Purifying and Reversible Physical Processes

Starting from the observation that reversible processes cannot increase the purity of any input state, we study deterministic physical processes, which map a set of states to a set of pure states. Such a process must map any state to the same pure output, if purity is demanded for the input set of all states. But otherwise, when the input set is restricted, it is possible to find non-trivial purifying processes. For the most restricted case of only two input states, we completely characterize the output of any such map. We furthermore consider maps, which combine the property of purity and reversibility on a set of states, and we derive necessary and sufficient conditions on sets, which permit such processes.Comment: 5 pages, no figures, v2: only minimal change

On the Convergence of Ritz Pairs and Refined Ritz Vectors for Quadratic Eigenvalue Problems

For a given subspace, the Rayleigh-Ritz method projects the large quadratic eigenvalue problem (QEP) onto it and produces a small sized dense QEP. Similar to the Rayleigh-Ritz method for the linear eigenvalue problem, the Rayleigh-Ritz method defines the Ritz values and the Ritz vectors of the QEP with respect to the projection subspace. We analyze the convergence of the method when the angle between the subspace and the desired eigenvector converges to zero. We prove that there is a Ritz value that converges to the desired eigenvalue unconditionally but the Ritz vector converges conditionally and may fail to converge. To remedy the drawback of possible non-convergence of the Ritz vector, we propose a refined Ritz vector that is mathematically different from the Ritz vector and is proved to converge unconditionally. We construct examples to illustrate our theory.Comment: 20 page

Further results on non-diagonal Bianchi type III vacuum metrics

We present the derivation, for these vacuum metrics, of the Painlev\'e VI equation first obtained by Christodoulakis and Terzis, from the field equations for both minkowskian and euclidean signatures. This allows a complete discussion and the precise connection with some old results due to Kinnersley. The hyperk\"ahler metrics are shown to belong to the Multi-Centre class and for the cases exhibiting an integrable geodesic flow the relevant Killing tensors are given. We conclude by the proof that for the Bianchi B family, excluding type III, there are no hyperk\"ahler metrics.Comment: 21 pages, no figure

Fast linear algebra is stable

In an earlier paper, we showed that a large class of fast recursive matrix multiplication algorithms is stable in a normwise sense, and that in fact if multiplication of $n$-by-$n$ matrices can be done by any algorithm in $O(n^{\omega + \eta})$ operations for any $\eta > 0$, then it can be done stably in $O(n^{\omega + \eta})$ operations for any $\eta > 0$. Here we extend this result to show that essentially all standard linear algebra operations, including LU decomposition, QR decomposition, linear equation solving, matrix inversion, solving least squares problems, (generalized) eigenvalue problems and the singular value decomposition can also be done stably (in a normwise sense) in $O(n^{\omega + \eta})$ operations.Comment: 26 pages; final version; to appear in Numerische Mathemati

A Matrix Hyperbolic Cosine Algorithm and Applications

In this paper, we generalize Spencer's hyperbolic cosine algorithm to the matrix-valued setting. We apply the proposed algorithm to several problems by analyzing its computational efficiency under two special cases of matrices; one in which the matrices have a group structure and an other in which they have rank-one. As an application of the former case, we present a deterministic algorithm that, given the multiplication table of a finite group of size $n$, it constructs an expanding Cayley graph of logarithmic degree in near-optimal O(n^2 log^3 n) time. For the latter case, we present a fast deterministic algorithm for spectral sparsification of positive semi-definite matrices, which implies an improved deterministic algorithm for spectral graph sparsification of dense graphs. In addition, we give an elementary connection between spectral sparsification of positive semi-definite matrices and element-wise matrix sparsification. As a consequence, we obtain improved element-wise sparsification algorithms for diagonally dominant-like matrices.Comment: 16 pages, simplified proof and corrected acknowledging of prior work in (current) Section

Brane-Antibrane Inflation in Orbifold and Orientifold Models

We analyse the cosmological implications of brane-antibrane systems in string-theoretic orbifold and orientifold models. In a class of realistic models, consistency conditions require branes and antibranes to be stuck at different fixed points, and so their mutual attraction generates a potential for one of the radii of the underlying torus or the 4D string dilaton. Assuming that all other moduli have been fixed by string effects, we find that this potential leads naturally to a period of cosmic inflation with the radion or dilaton field as the inflaton. The slow-roll conditions are satisfied more generically than if the branes were free to move within the space. The appearance of tachyon fields at certain points in moduli space indicates the onset of phase transitions to different non-BPS brane systems, providing ways of ending inflation and reheating the corresponding observable brane universe. In each case we find relations between the inflationary parameters and the string scale to get the correct spectrum of density perturbations. In some examples the small numbers required as inputs are no smaller than 0.01, and are the same small quantities which are required to explain the gauge hierarchy.Comment: 30 pages, 2 figures. Substantial changes on version 1. New cosmological scenarios proposed including the dilaton as the inflaton. Main conclusions unchange

Lessons learned from infertility investigations in the public sector

Objectives. To determine the main factors causing infertility in an urban, tertiary hospital population. To establish if any such major causal factor could be used to rationalise and improve the service for infertile couples in the public sector.Design. A retrospective analysis of the hospital records of 206 women who had a tubal patency test (hysterosalpingogram) performed and the results of the investigations performed in the couples with infertility.Results. Of the 206 women 38 (18.5%) had normal fallopian tubes on hysterosalpingogram; 33 (16%) had unilateral obstruction and 135 (65.5%) had bilateral tubal obstruction. Of the latter group 81 (60%) had significant hydrosalpinges. Semen analysis results in 148 partners (71.8%) demonstrated a normal count in 85 (62%), normal motility in 70 (51%) and normal morphology in only 25 (18%). Testing for ovulation (mid-luteal phase progesterone) was positive in 91 of 124 women tested (73%). Compliance, technical and logistical problems were encountered with both semen analysis and mid-luteal phase progesterone tests.Conclusions. Infertility is a major problem in South Africa, with limited resources for investigation and treatment in the public sector. Tubal factor infertility was the most common cause of infertility demonstrated in this study. In the presence of bilateral tubal obstruction with hydrosalpinges the prognosis is so poor that unless assisted reproductive techniques are available and affordable, further infertility investigations do not seem justified. Recommendations on an approach to the infertile couple in the public sector is outlined

A weakly stable algorithm for general Toeplitz systems

We show that a fast algorithm for the QR factorization of a Toeplitz or Hankel matrix A is weakly stable in the sense that R^T.R is close to A^T.A. Thus, when the algorithm is used to solve the semi-normal equations R^T.Rx = A^Tb, we obtain a weakly stable method for the solution of a nonsingular Toeplitz or Hankel linear system Ax = b. The algorithm also applies to the solution of the full-rank Toeplitz or Hankel least squares problem.Comment: 17 pages. An old Technical Report with postscript added. For further details, see http://wwwmaths.anu.edu.au/~brent/pub/pub143.htm

Debris disk size distributions: steady state collisional evolution with P-R drag and other loss processes

We present a new scheme for determining the shape of the size distribution, and its evolution, for collisional cascades of planetesimals undergoing destructive collisions and loss processes like Poynting-Robertson drag. The scheme treats the steady state portion of the cascade by equating mass loss and gain in each size bin; the smallest particles are expected to reach steady state on their collision timescale, while larger particles retain their primordial distribution. For collision-dominated disks, steady state means that mass loss rates in logarithmic size bins are independent of size. This prescription reproduces the expected two phase size distribution, with ripples above the blow-out size, and above the transition to gravity-dominated planetesimal strength. The scheme also reproduces the expected evolution of disk mass, and of dust mass, but is computationally much faster than evolving distributions forward in time. For low-mass disks, P-R drag causes a turnover at small sizes to a size distribution that is set by the redistribution function (the mass distribution of fragments produced in collisions). Thus information about the redistribution function may be recovered by measuring the size distribution of particles undergoing loss by P-R drag, such as that traced by particles accreted onto Earth. Although cross-sectional area drops with 1/age^2 in the PR-dominated regime, dust mass falls as 1/age^2.8, underlining the importance of understanding which particle sizes contribute to an observation when considering how disk detectability evolves. Other loss processes are readily incorporated; we also discuss generalised power law loss rates, dynamical depletion, realistic radiation forces and stellar wind drag.Comment: Accepted for publication by Celestial Mechanics and Dynamical Astronomy (special issue on EXOPLANETS

Amides do not always work: observation of guest binding in an amide-functionalised porous host

An amide-functionalised metal organic frame-work (MOF) material, MFM-136, shows a high CO2 uptake of 12.6 mmol g-1 at 20 bar and 298 K. MFM-136 is the first example of acylamide pyrimidyl isophthalate MOF without open metal sites, and thus provides a unique platform to study guest bind-ing, particularly the role of free amides. Neutron diffraction reveals that, surprisingly, there is no direct binding between the adsorbed CO2/CH4 molecules and the pendant amide group in the pore. This observation has been confirmed un-ambiguously by inelastic neutron spectroscopy. This suggests that introduction of functional groups solely may not neces-sarily induce specific guest-host binding in porous materials, but it is a combination of pore size, geometry, and functional group that leads to enhanced gas adsorption properties