First-principles quantum dynamics in interacting Bose gases II: stochastic gauges

First principles simulations of the quantum dynamics of interacting Bose gases using the stochastic gauge representation are analysed. In a companion paper, we showed how the positive-P representation can be applied to these problems using stochastic differential equations. That method, however, is limited by increased sampling error as time evolves. Here, we show how the sampling error can be greatly reduced and the simulation time significantly extended using stochastic gauges. In particular, local stochastic gauges (a subset) are investigated. Improvements are confirmed in numerical calculations of single-, double- and multi-mode systems in the weak-mode coupling regime. Convergence issues are investigated, including the recognition of two modes by which stochastic equations produced by phase-space methods in general can diverge: movable singularities and a noise-weight relationship. The example calculated here displays wave-like behaviour in spatial correlation functions propagating in a uniform 1D gas after a sudden change in the coupling constant. This could in principle be tested experimentally using Feshbach resonance methods.

Experimental tests of quantum nonlinear dynamics in atom optics

Cold atoms in optical potentials provide an ideal test bed to explore quantum nonlinear dynamics. Atoms are prepared in a magneto-optic trap or as a dilute Bose-Einstein condensate and subjected to a far detuned optical standing wave that is modulated. They exhibit a wide range of dynamics, some of which can be explained by classical theory while other aspects show the underlying quantum nature of the system. The atoms have a mixed phase space containing regions of regular motion which appear as distinct peaks in the atomic momentum distribution embedded in a sea of chaos. The action of the atoms is of the order of Planck's constant, making quantum effects significant. This tutorial presents a detailed description of experiments measuring the evolution of atoms in time-dependent optical potentials. Experimental methods are developed providing means for the observation and selective loading of regions of regular motion. The dependence of the atomic dynamics on the system parameters is explored and distinct changes in the atomic momentum distribution are observed which are explained by the applicable quantum and classical theory. The observation of a bifurcation sequence is reported and explained using classical perturbation theory. Experimental methods for the accurate control of the momentum of an ensemble of atoms are developed. They use phase space resonances and chaotic transients providing novel ensemble atomic beamsplitters. The divergence between quantum and classical nonlinear dynamics is manifest in the experimental observation of dynamical tunnelling. It involves no potential barrier. However a constant of motion other than energy still forbids classically this quantum allowed motion. Atoms coherently tunnel back and forth between their initial state of oscillatory motion and the state 180 out of phase with the initial state.

Fock-state dynamics in Raman photoassociation of Bose-Einstein condensates

By stochastic modeling of the process of Raman photoassociation of Bose-Einstein condensates, we show that, the farther the initial quantum state is from a coherent state, the farther the one-dimensional predictions are from those of the commonly used zero-dimensional approach. We compare the dynamics of condensates, initially in different quantum states, finding that, even when the quantum prediction for an initial coherent state is relatively close to the Gross-Pitaevskii prediction, an initial Fock state gives qualitatively different predictions. We also show that this difference is not present in a single-mode type of model, but that the quantum statistics assume a more important role as the dimensionality of the model is increased. This contrasting behavior in different dimensions, well known with critical phenomena in statistical mechanics, makes itself plainly visible here in a mesoscopic system and is a strong demonstration of the need to consider physically realistic models of interacting condensates.

Green turtle somatic growth dynamics in a spatially disjunct Great Barrier Reef metapopulation

The growth dynamics of green sea turtles resident in four separate foraging grounds of the southern Great Barrier Reef genetic stock were assessed using a nonparametric regression modeling approach. Juveniles recruit to these grounds at the same size, but grow at foraging-ground-dependent rates that result in significant differences in expected size- or age-at-maturity. Mean age-at-maturity was estimated to vary from 25-50 years depending on the ground. This stock comprises mainly the same mtDNA haplotype, so geographic variability might be due to local environmental conditions rather than genetic factors, although the variability was not a function of latitudinal variation in environmental conditions or whether the food stock was seagrass or algae. Temporal variability in growth rates was evident in response to local environmental stochasticity, so geographic variability might be due to local food stock dynamics. Despite such variability, the expected size-specific growth rate function at all grounds displayed a similar nonmonotonic growth pattern with a juvenile growth spurt at 60-70 cm curved carapace length, (CCL) or 15-20 years of age. Sex-specific growth differences were also evident with females tending to grow faster than similar-sized males after the Juvenile growth spurt. It is clear that slow sex-specific growth displaying both spatial and temporal variability and a juvenile growth spurt are distinct growth behaviors of green turtles from this stock.

Undular tidal bore dynamics in the Daly Estuary, Northern Australia

Measurements in the macro-tidal Daly Estuary show that the presence of an undular tidal bore contributed negligibly to the dissipation of tidal energy. No recirculation bubble was observed between a trough and the following wave crest in the lee waves following the undular bore. This differs to stationary undular bores in laboratory experiments at larger Froude numbers where a recirculation bubble exists. Secondary motions and the turbulence generated by the undular bore had no measurable influence on the sediment transport. This situation contrasts with the intense sediment resuspension observed in breaking tidal bores. The tidally averaged sediment budget in the Daly Estuary was controlled by the asymmetry of tidal currents. The undular bore may widen the river by breaking along the banks that it undercuts, leading to bank slippage. A patch of river-wide macro-turbulence of 3-min duration occurred about 20 min after the passage of the bore during accelerating tidal currents. (C) 2004 Elsevier Ltd. All rights reserved.

Modelling, nitrogen dynamics in sugarcane systems: Recent advances and applications

Substantial amounts of nitrogen (N) fertiliser are necessary for commercial sugarcane production because of the large biomass produced by sugarcane crops. Since this fertiliser is a substantial input cost and has implications if N is lost to the environment, there are pressing needs to optimise the supply of N to the crops' requirements. The complexity of the N cycle and the strong influence of climate, through its moderation of N transformation processes in the soil and its impact on N uptake by crops, make simulation-based approaches to this N management problem attractive. In this paper we describe the processes to be captured in modelling soil and plant N dynamics in sugarcane systems, and review the capability for modelling these processes. We then illustrate insights gained into improved management of N through simulation-based studies for the issues of crop residue management, irrigation management and greenhouse gas emissions. We conclude by identifying processes not currently represented in the models used for simulating N cycling in sugarcane production systems, and illustrate ways in which these can be partially overcome in the short term. (c) 2005 Elsevier B.V. All rights reserved.

Patch dynamics in grazed subtropical native pastures in south-east Queensland

Patch formation is common in grazed grasslands but the mechanisms involved in the formation and maintenance of patches are not clear. To increase our knowledge on this subject we examined possible reasons for patch formation and the influence of management on changes between patch states in three experiments in native pasture communities in the Crows Nest district, south-east Queensland. In these communities, small-scale patches (tall grassland (dominated by large and medium tussock grasses), short swards (dominated by short tussock grasses and sedges), and lawns (dominated by stoloniferous and/or rhizomatous grasses)) are readily apparent. We hypothesized that the formation of short sward and lawn patches in areas of tall grassland was due to combinations of grazing and soil fertility effects. This was tested in Experiment 1 by applying a factorial combination of defoliation, nutrient application and transplants of short tussock and stoloniferous species to a uniform area of tall grassland. Total species density declined during the experiment, was lower with high nutrient applications, but was not affected by defoliation. There were significant changes in abundance of species that provided support for our hypotheses. With light defoliation and low nutrients, the tall grassland remained dominated by large tussock grasses and contained considerable amounts of forbs. With heavy defoliation, the pastures were dominated by medium tussock grasses and there were significant decreases in forbs and increases in sedges (mainly with low nutrients) and stoloniferous grasses (mainly with high nutrients). Total germinable seed densities and those of most species groups were significantly lower in the heavy defoliation than the light defoliation plots. Total soil seed numbers were not affected by nutrient application but there were fewer seeds of the erect forbs and more sedge seeds in plots with high nutrients. The use of resting from grazing and fire to manage transitions between patches was tested. In Experiment 2, changes in species density and abundance were measured for 5 years in the three patch types with and without grazing. Experiment 3 examined the effects of fire, grazing and resting on short sward patches over 4 years. In Experiment 2, total species density was lower in lawn than short sward or tall grassland patches, and there were more species of erect forbs than other plant groups in all patch types. The lawn patches were originally dominated by Cynodon spp. This dominance continued with grazing but in ungrazed patches the abundance of Cynodon spp. declined and that of forbs increased. In the short sward patches, dominance of short tussock grasses continued with grazing but in ungrazed plots their abundance declined while that of large tussock grasses increased. The tall grassland patches remained dominated by large and medium tussock species. In Experiment 3, fire had no effect on species abundance. On the grazed plots the short tussock grasses remained dominant but where the plots were rested from grazing the small tussock grasses declined and the large tussock grasses increased in abundance. The slow and relatively small changes in these experiments over 4 or 5 years showed how stable the composition of these pastures is, and that rapid changes between patch types are unlikely.

Analysis of articulatory dynamics in dysarthria following brain injury in childhood using electromagnetic articulography and electropalatography

Two physiological assessments, electromagnetic articulography (EMA) and electropalatography (EPG), were used simultaneously to investigate the articulatory dynamics in an 18-year-old male with dysarthria 9 years following traumatic brain injury (TBI). Eight words consisting of /t, s, integral, k/ in word initial and word final positions were produced up to 10 times. A nonneurologically impaired male served as a control subject. Six parameters were analyzed using EMA: velocity, acceleration, deceleration, distance, duration, and motion path of tongue movements. Using EPG, the pattern and amount of tongue-to-palate contact and the duration of the closure/constriction phase of each consonant produced were assessed. Timing disturbances in the TBI speaker's speech were highlighted in perceptual assessments in the form of prolonged phonemes and a reduced speech rate. EMA analysis revealed that the approach and release phase durations of the consonant productions were within normal limits. Kinematic strategies such as decreased velocity and decreased distances traveled by the tongue, however, may have counterbalanced each other to produce these appropriate results. EPG examination revealed significantly longer closure/constriction phase periods, which may have contributed to the prolonged phonemes and reduced speech rate observed. The implications of these findings for the development of treatment programs for dysarthria subsequent to TBI will be highlighted.

Universality of quantum critical dynamics in a planar optical parametric oscillator

We analyze the critical quantum fluctuations in a coherently driven planar optical parametric oscillator. We show that the presence of transverse modes combined with quantum fluctuations changes the behavior of the quantum image critical point. This zero-temperature nonequilibrium quantum system has the same universality class as a finite-temperature magnetic Lifshitz transition.

First-principles quantum dynamics in interacting Bose gases: I. The positive P representation

The performance of the positive P phase-space representation for exact many- body quantum dynamics is investigated. Gases of interacting bosons are considered, where the full quantum equations to simulate are of a Gross-Pitaevskii form with added Gaussian noise. This method gives tractable simulations of many-body systems because the number of variables scales linearly with the spatial lattice size. An expression for the useful simulation time is obtained, and checked in numerical simulations. The dynamics of first-, second- and third-order spatial correlations are calculated for a uniform interacting 1D Bose gas subjected to a change in scattering length. Propagation of correlations is seen. A comparison is made with other recent methods. The positive P method is particularly well suited to open systems as no conservation laws are hard-wired into the calculation. It also differs from most other recent approaches in that there is no truncation of any kind.

Network topology and the evolution of dynamics in an artificial genetic regulatory network model created by whole genome duplication and divergence

Topological measures of large-scale complex networks are applied to a specific artificial regulatory network model created through a whole genome duplication and divergence mechanism. This class of networks share topological features with natural transcriptional regulatory networks. Specifically, these networks display scale-free and small-world topology and possess subgraph distributions similar to those of natural networks. Thus, the topologies inherent in natural networks may be in part due to their method of creation rather than being exclusively shaped by subsequent evolution under selection. The evolvability of the dynamics of these networks is also examined by evolving networks in simulation to obtain three simple types of output dynamics. The networks obtained from this process show a wide variety of topologies and numbers of genes indicating that it is relatively easy to evolve these classes of dynamics in this model. (c) 2006 Elsevier Ireland Ltd. All rights reserved.

Scaling and non-scaling of muscle activity, kinematics, and dynamics in sit-ups with different degrees of difficulty

The purpose of this study was to investigate how the CNS adjusts motor patterns for variants of a complex axial movement-the situp. Adjustments were induced by changing the support surface contact and mass distribution of the body. Healthy adults performed straight-legged sit-ups, 3 s in duration, with support added to or removed from the lumbar trunk, or with mass added to the head or to the legs. Each of these interventions either increased or decreased the difficulty of the task. The study addressed the extent to which changes in sit-up difficulty are compensated by scaling of muscle activity, kinematics, and dynamics versus the extent to which they are compensated by changing discretely the motor pattern. The analysis of muscle activity, kinematics, and dynamics focused on the first 30-40% of the sit-up-the trunk flexion phase-since this is the most critical part of the movement. Our results demonstrate that, in some respects, sit-up kinematics and dynamics scaled with difficulty, but in other respects, they did not. Muscle activity also scaled, in many respects, but in more difficult sit-ups, abdominal flexor activity decreased instead of increased. Non-scaling changes in these parameters suggest that complex movements, such as the sit-up, may require discrete changes in motor pattern in order to deal with large loads, which challenge the available leverage. (C) 2005 Elsevier Ltd. All rights reserved.

Spatial patterns of aerobic and anaerobic mineralization rates and oxygen penetration dynamics in coral reef sediments

Oxygen consumption rates (OCR), aerobic mineralization and sulfate reduction rates (SRR) were studied in the permeable carbonate reef sediments of Heron Reef, Australia. We selected 4 stations with different hydrodynamic regimes for this study. In situ oxygen penetration into the sediments was measured with an autonomous microsensor profiler. Areal OCR were quantified from the measured oxygen penetration depth and volumetric OCR. Oxygen penetration and dynamics (median penetration depths at the 4 stations ranged between 0.3 and 2.2 cm), OCR (median 57 to 196 mmol C m(-2) d(-1)), aerobic mineralization (median 24 to 176 mmol C m(-2) d(-1)) and SRR (median 9 to 42 mmol C m(-2) d(-1)) were highly variable between sites. The supply of oxygen by pore water advection was a major cause for high mineralization rates by stimulating aerobic mineralization at all sites. However, estimated bottom water filtration rates could not explain the differences in volumetric OCR and SRR between the 4 stations. This suggests that local mineralization rates are additionally controlled by factors other than current driven pore water advection, e.g. by the distribution of the benthic fauna or by local differences in labile organic carbon supply from sources such as benthic photosynthesis. Carbon mineralization rates were among the highest reported for coral reef sediments, stressing the role of these sediments in the functioning of the reef ecosystem.

Transcriptional network dynamics in macrophage activation

Transcriptional regulatory networks govern cell differentiation and the cellular response to external stimuli. However, mammalian model systems have not yet been accessible for network analysis. Here, we present a genome-wide network analysis of the transcriptional regulation underlying the mouse macrophage response to bacterial lipopolysaccharide (LPS). Key to uncovering the network structure is our combination of time-series cap analysis of gene expression with in silico prediction of transcription factor binding sites. By integrating microarray and qPCR time-series expression data with a promoter analysis, we find dynamic subnetworks that describe how signaling pathways change dynamically during the progress of the macrophage LPS response, thus defining regulatory modules characteristic of the inflammatory response. In particular, our integrative analysis enabled us to suggest novel roles for the transcription factors ATF-3 and NRF-2 during the inflammatory response. We believe that our system approach presented here is applicable to understanding cellular differentiation in higher eukaryotes. (c) 2006 Elsevier Inc. All rights reserved.