Quantum dynamics with stochastic gauge simulations

The general idea of a stochastic gauge representation is introduced and compared with more traditional phase-space expansions, like the Wigner expansion. Stochastic gauges can be used to obtain an infinite class of positive-definite stochastic time-evolution equations, equivalent to master equations, for many systems including quantum time evolution. The method is illustrated with a variety of simple examples ranging from astrophysical molecular hydrogen production, through to the topical problem of Bose-Einstein condensation in an optical trap and the resulting quantum dynamics.

Gauge Poisson representations for birth/death master equations

Poisson representation techniques provide a powerful method for mapping master equations for birth/death processes -- found in many fields of physics, chemistry and biology -- into more tractable stochastic differential equations. However, the usual expansion is not exact in the presence of boundary terms, which commonly occur when the differential equations are nonlinear. In this paper, a gauge Poisson technique is introduced that eliminates boundary terms, to give an exact representation as a weighted rate equation with stochastic terms. These methods provide novel techniques for calculating and understanding the effects of number correlations in systems that have a master equation description. As examples, correlations induced by strong mutations in genetics, and the astrophysical problem of molecule formation on microscopic grain surfaces are analyzed. Exact analytic results are obtained that can be compared with numerical simulations, demonstrating that stochastic gauge techniques can give exact results where standard Poisson expansions are not able to.

Polymorphic CAG repeat length in the androgen receptor gene and association with neurodegeneration in a heterozygous female carrier of Kennedy's disease

Kennedy's disease (spinobulbar muscular atrophy) is an X-linked form of motor neuron disease affecting adult males carrying a CAG trinucleotide repeat expansion within the androgen receptor gene. While expression of Kennedy's disease is thought to be confined to males carrying the causative mutation, subclinical manifestations have been reported in a few female carriers of the disease. The reasons that females are protected from the disease are not clear, especially given that all other diseases caused by CAG expansions display dominant expression. In the current study, we report the identification of a heterozygote female carrying the Kennedy's disease mutation who was clinically diagnosed with motor neuron disease. We describe analysis of CAG repeat number in this individual as well as 33 relatives within the pedigree, including two male carriers of the Kennedy's mutation. The female heterozygote carried one expanded allele of the androgen receptor gene with CAG repeats numbering in the Kennedy's disease range (44 CAGs), with the normal allele numbering in the upper-normal range (28 CAGs). The subject has two sons, one of whom carries the mutant allele of the gene and has been clinically diagnosed with Kennedy's disease, whilst the other son carries the second allele of the gene with CAGs numbering in the upper normal range and displays a normal phenotype. This coexistence of motor neuron disease and the presence of one expanded allele and one allele at the upper limit of the normal range may be a coincidence. However, we hypothesize that the expression of the Kennedy's disease mutation combined with a second allele with a large but normal CAG repeat sequence may have contributed to the motor neuron degeneration displayed in the heterozygote female and discuss the possible reasons for phenotypic expression in particular individuals.

Green-Ampt approximations

The solution to the Green and Ampt infiltration equation is expressible in terms of the Lambert W-1 function. Approximations for Green and Ampt infiltration are thus derivable from approximations for the W-1 function and vice versa. An infinite family of asymptotic expansions to W-1 is presented. Although these expansions do not converge near the branch point of the W function (corresponds to Green-Ampt infiltration with immediate ponding), a method is presented for approximating W-1 that is exact at the branch point and asymptotically, with interpolation between these limits. Some existing and several new simple and compact yet robust approximations applicable to Green-Ampt infiltration and flux are presented, the most accurate of which has a maximum relative error of 5 x 10(-5)%. This error is orders of magnitude lower than any existing analytical approximations. (c) 2005 Elsevier Ltd. All rights reserved.

Steepness expansion for free surface flows in coastal aquifers

Free surface flow of groundwater in aquifers has been studied since the early 1960s. Previous investigations have been based on the Boussinesq equation, derived from the non-linear kinematic boundary condition. In fact, the Boussinesq equation is the zeroth-order equation in the shallow-water expansion. A key assumption in this expansion is that the mean thickness of the aquifer is small compared with a reference length, normally taken to be the linear decay length. In this study, we re-examine the expansion scheme for free surface groundwater flows, and propose a new expansion wherein the shallow-water assumption is replaced by a steepness assumption. A comparison with experimental data shows that the new model provides a better prediction of water table levels than the conventional shallow-water expansion. The applicable ranges of the two expansions are exhibited. (c) 2004 Elsevier B.V. All rights reserved.

Molecular analysis reveals tighter social regulation of immigration in patrilocal populations than in matrilocal populations

Human social organization can deeply affect levels of genetic diversity. This fact implies that genetic information can be used to study social structures, which is the basis of ethnogenetics. Recently, methods have been developed to extract this information from genetic data gathered from subdivided populations that have gone through recent spatial expansions, which is typical of most human populations. Here, we perform a Bayesian analysis of mitochondrial and Y chromosome diversity in three matrilocal and three patrilocal groups from northern Thailand to infer the number of males and females arriving in these populations each generation and to estimate the age of their range expansion. We find that the number of male immigrants is 8 times smaller in patrilocal populations than in matrilocal populations, whereas women move 2.5 times more in patrilocal populations than in matrilocal populations. In addition to providing genetic quantification of sex-specific dispersal rates in human populations, we show that although men and women are exchanged at a similar rate between matrilocal populations, there are far fewer men than women moving into patrilocal populations. This finding is compatible with the hypothesis that men are strictly controlling male immigration and promoting female immigration in patrilocal populations and that immigration is much less regulated in matrilocal populations.

Adverse feedback sequences in exploited marine systems: Are deliberate interruptive actions warranted?

Several mechanisms for self-enhancing feedback instabilities in marine ecosystems are identified and briefly elaborated. It appears that adverse phases of operation may be abruptly triggered by explosive breakouts in abundance of one or more previously suppressed populations. Moreover, an evident capacity of marine organisms to accomplish extensive geographic habitat expansions may expand and perpetuate a breakout event. This set of conceptual elements provides a framework for interpretation of a sequence of events that has occurred in the Northern Benguela Current Large Marine Ecosystem (off south-western Africa). This history can illustrate how multiple feedback loops might interact with one another in unanticipated and quite malignant ways, leading not only to collapse of customary resource stocks but also to degradation of the ecosystem to such an extent that disruption of customary goods and services may go beyond fisheries alone to adversely affect other major global ecosystem concerns (e.g. proliferations of jellyfish and other slimy, stingy, toxic and/or noxious organisms, perhaps even climate change itself, etc.). The wisdom of management interventions designed to interrupt an adverse mode of feedback operation is pondered. Research pathways are proposed that may lead to improved insights needed: (i) to avoid potential 'triggers' that might set adverse phases of feedback loop operation into motion; and (ii) to diagnose and properly evaluate plausible actions to reverse adverse phases of feedback operation that might already have been set in motion. These pathways include the drawing of inferences from available 'quasi-experiments' produced either by short-term climatic variation or inadvertently in the course of biased exploitation practices, and inter-regional applications of the comparative method of science.

Anomalous Excitation Spectra of Frustrated Quantum Antiferromagnets

We use series expansions to study the excitation spectra of spin-1/2 antiferromagnets on anisotropic triangular lattices. For the isotropic triangular lattice model (TLM), the high-energy spectra show several anomalous features that differ strongly from linear spin-wave theory (LSWT). Even in the Neel phase, the deviations from LSWT increase sharply with frustration, leading to rotonlike minima at special wave vectors. We argue that these results can be interpreted naturally in a spinon language and provide an explanation for the previously observed anomalous finite-temperature properties of the TLM. In the coupled-chains limit, quantum renormalizations strongly enhance the one-dimensionality of the spectra, in agreement with experiments on Cs2CuCl4.

Critical behavior of one-particle spectral weights in the transverse Ising model

We investigate the critical behavior of the spectral weight of a single quasiparticle, one of the key observables in experiment, for the particular case of the transverse Ising model. Series expansions are calculated for the linear chain and the square and simple cubic lattices. For the chain model, a conjectured exact result is discovered. For the square and simple cubic lattices, series analyses are used to estimate the critical exponents. The results agree with the general predictions of Sachdev [Quantum Phase Transitions (Cambridge University Press, Cambridge, England, 1999)].

On functionally-fitted Runge-Kutta methods

Functionally-fitted methods are generalizations of collocation techniques to integrate an equation exactly if its solution is a linear combination of a chosen set of basis functions. When these basis functions are chosen as the power functions, we recover classical algebraic collocation methods. This paper shows that functionally-fitted methods can be derived with less restrictive conditions than previously stated in the literature, and that other related results can be derived in a much more elegant way. The novelty in our approach is to fully retain the collocation framework without reverting back into derivations based on cumbersome Taylor series expansions.

Excitation spectra of the spin-1/2 triangular-lattice Heisenberg antiferromagnet

We use series expansion methods to calculate the dispersion relation of the one-magnon excitations for the spin-(1)/(2) triangular-lattice nearest-neighbor Heisenberg antiferromagnet above a three-sublattice ordered ground state. Several striking features are observed compared to the classical (large-S) spin-wave spectra. Whereas, at low energies the dispersion is only weakly renormalized by quantum fluctuations, significant anomalies are observed at high energies. In particular, we find rotonlike minima at special wave vectors and strong downward renormalization in large parts of the Brillouin zone, leading to very flat or dispersionless modes. We present detailed comparison of our calculated excitation energies in the Brillouin zone with the spin-wave dispersion to order 1/S calculated recently by Starykh, Chubukov, and Abanov [Phys. Rev. B74, 180403(R) (2006)]. We find many common features but also some quantitative and qualitative differences. We show that at temperatures as low as 0.1J the thermally excited rotons make a significant contribution to the entropy. Consequently, unlike for the square lattice model, a nonlinear sigma model description of the finite-temperature properties is only applicable at temperatures < 0.1J. Finally, we review recent NMR measurements on the organic compound kappa-(BEDT-TTF)(2)Cu-2(CN)(3). We argue that these are inconsistent with long-range order and a description of the low-energy excitations in terms of interacting magnons, and that therefore a Heisenberg model with only nearest-neighbor exchange does not offer an adequate description of this material.

Accelerated leap methods for simulating discrete stochastic chemical kinetics

Biologists are increasingly conscious of the critical role that noise plays in cellular functions such as genetic regulation, often in connection with fluctuations in small numbers of key regulatory molecules. This has inspired the development of models that capture this fundamentally discrete and stochastic nature of cellular biology - most notably the Gillespie stochastic simulation algorithm (SSA). The SSA simulates a temporally homogeneous, discrete-state, continuous-time Markov process, and of course the corresponding probabilities and numbers of each molecular species must all remain positive. While accurately serving this purpose, the SSA can be computationally inefficient due to very small time stepping so faster approximations such as the Poisson and Binomial τ-leap methods have been suggested. This work places these leap methods in the context of numerical methods for the solution of stochastic differential equations (SDEs) driven by Poisson noise. This allows analogues of Euler-Maruyuma, Milstein and even higher order methods to be developed through the Itô-Taylor expansions as well as similar derivative-free Runge-Kutta approaches. Numerical results demonstrate that these novel methods compare favourably with existing techniques for simulating biochemical reactions by more accurately capturing crucial properties such as the mean and variance than existing methods.