Superchemistry: dynamics of coupled atomic and molecular Bose condensates

We analyze the dynamics of a dilute, trapped Bose-condensed atomic gas coupled to a diatomic molecular Bose gas by coherent Raman transitions. This system is shown to result in a new type of “superchemistry,” in which giant collective oscillations between the atomic and the molecular gas can occur. The phenomenon is caused by stimulated emission of bosonic atoms or molecules into their condensate phases.

Quantum dynamics of three coupled atomic Bose-Einstein condensates

The simplest model of three coupled Bose-Einstein condensates is investigated using a group theoretical method. The stationary solutions are determined using the SU(3) group under the mean-field approximation. This semiclassical analysis, using system symmetries, shows a transition in the dynamics of the system from self trapping to delocalization at a critical value for the coupling between the condensates. The global dynamics are investigated by examination of the stable points, and our analysis shows that the structure of the stable points depends on the ratio of the condensate coupling to the particle-particle interaction, and undergoes bifurcations as this ratio is varied. This semiclassical model is compared to a full quantum treatment, which also displays a dynamical transition. The quantum case has collapse and revival sequences superimposed on the semiclassical dynamics, reflecting the underlying discreteness of the spectrum. Nonzero circular current states are also demonstrated as one of the higher-dimensional effects displayed in this system.

Substructure and Dynamics of the Fornax Cluster

We present the first dynamical analysis of a galaxy cluster to include a large fraction of dwarf galaxies. Our sample of 108 Fornax Cluster members measured with the UK Schmidt Telescope FLAIR-II spectrograph contains 55 dwarf galaxies (15.5 > b(j) > 18.0 or -16 > M-B > -13.5). H alpha emission shows that of the dwarfs are star forming, twice the fraction implied by morphological classifications. The total sample has a mean velocity of 1493 +/- 36 kms s(-1) and a velocity dispersion of 374 +/- 26 km s(-1). The dwarf galaxies form a distinct population: their velocity dispersion (429 +/- 41 km s(-1)) is larger than that of the giants () at the 98% confidence level. This suggests that the dwarf population is dominated by infalling objects whereas the giants are virialized. The Fornax system has two components, the main Fornax Cluster centered on NGC 1399 with cz = 1478 km s(-1) and sigma (cz) = 370 km s(-1) and a subcluster centered 3 degrees to the southwest including NGC 1316 with cz = 1583 km s(-1) and sigma (cz) = 377 km s(-1). This partition is preferred over a single cluster at the 99% confidence level. The subcluster, a site of intense star formation, is bound to Fornax and probably infalling toward the cluster core for the first time. We discuss the implications of this substructure for distance estimates of the Fornax Cluster. We determine the cluster mass profile using the method of Diaferio, which does not assume a virialized sample. The mass within a projected radius of 1.4 Mpc is (7 +/- 2) x 10(13) M-., and the mass-to-light ratio is 300 +/- 100 M-./L-.. The mass is consistent with values derived from the projected mass virial estimator and X-ray measurements at smaller radii.

Evolution of Wolbachia pipientis transmission dynamics in insects

Wolbachia pipientis is an intracellular bacterial parasite of arthropods that enhances its transmission by manipulating host reproduction, most commonly by inducing cytoplasmic incompatibility. The discovery of isolates with modified cytoplasmic incompatibility phenotypes and others with novel virulence properties is an indication of the potential breadth of evolutionary strategies employed by Wolbachia.

Non-Equilibrium Reaction Rates in the Macroscopic Chemistry Method for DSMC Calculations

The Direct Simulation Monte Carlo (DSMC) method is used to simulate the flow of rarefied gases. In the Macroscopic Chemistry Method (MCM) for DSMC, chemical reaction rates calculated from local macroscopic flow properties are enforced in each cell. Unlike the standard total collision energy (TCE) chemistry model for DSMC, the new method is not restricted to an Arrhenius form of the reaction rate coefficient, nor is it restricted to a collision cross-section which yields a simple power-law viscosity. For reaction rates of interest in aerospace applications, chemically reacting collisions are generally infrequent events and, as such, local equilibrium conditions are established before a significant number of chemical reactions occur. Hence, the reaction rates which have been used in MCM have been calculated from the reaction rate data which are expected to be correct only for conditions of thermal equilibrium. Here we consider artificially high reaction rates so that the fraction of reacting collisions is not small and propose a simple method of estimating the rates of chemical reactions which can be used in the Macroscopic Chemistry Method in both equilibrium and non-equilibrium conditions. Two tests are presented: (1) The dissociation rates under conditions of thermal non-equilibrium are determined from a zero-dimensional Monte-Carlo sampling procedure which simulates ‘intra-modal’ non-equilibrium; that is, equilibrium distributions in each of the translational, rotational and vibrational modes but with different temperatures for each mode; (2) The 2-D hypersonic flow of molecular oxygen over a vertical plate at Mach 30 is calculated. In both cases the new method produces results in close agreement with those given by the standard TCE model in the same highly nonequilibrium conditions. We conclude that the general method of estimating the non-equilibrium reaction rate is a simple means by which information contained within non-equilibrium distribution functions predicted by the DSMC method can be included in the Macroscopic Chemistry Method.

Using the level set method to model endogenous lava dome growth

Modeling volcanic phenomena is complicated by free-surfaces often supporting large rheological gradients. Analytical solutions and analogue models provide explanations for fundamental characteristics of lava flows. But more sophisticated models are needed, incorporating improved physics and rheology to capture realistic events. To advance our understanding of the flow dynamics of highly viscous lava in Peléean lava dome formation, axi-symmetrical Finite Element Method (FEM) models of generic endogenous dome growth have been developed. We use a novel technique, the level-set method, which tracks a moving interface, leaving the mesh unaltered. The model equations are formulated in an Eulerian framework. In this paper we test the quality of this technique in our numerical scheme by considering existing analytical and experimental models of lava dome growth which assume a constant Newtonian viscosity. We then compare our model against analytical solutions for real lava domes extruded on Soufrière, St. Vincent, W.I. in 1979 and Mount St. Helens, USA in October 1980 using an effective viscosity. The level-set method is found to be computationally light and robust enough to model the free-surface of a growing lava dome. Also, by modeling the extruded lava with a constant pressure head this naturally results in a drop in extrusion rate with increasing dome height, which can explain lava dome growth observables more appropriately than when using a fixed extrusion rate. From the modeling point of view, the level-set method will ultimately provide an opportunity to capture more of the physics while benefiting from the numerical robustness of regular grids.

Dependence of transient dynamics in a class-C laser upon variation of inversion with time

The transient statistics of a gain-switched coherently pumped class-C laser displays a linear correlation between the first passage time and subsequent peak intensity. Measurements are reported showing a positive or negative sign of this linear correlation, controlled through the switching time and the laser detuning. Further measurements of the small-signal laser gain combined with calculations involving a three-level laser model indicate that this sign fundamentally depends upon the way the laser inversion varies during the gain switching, despite the added dynamics of the laser polarization in the class-C laser. [S1050-2947(97)07112-6].

Regulation of photosynthetic pigments in micro-algae by multiple environmental factors: a dynamic balance hypothesis

Environmental effects on the concentration of photosynthetic pigments in micro-algae can be explained by dynamics of photosystem synthesis and deactivation. A model that couples photosystem losses to the relative cellular rates of energy harvesting (light absorption) and assimilation predicts optimal concentrations of light-harvesting pigments and balanced energy flow under environmental conditions that affect light availability and metabolic rates. Effects of light intensity, nutrient supply and temperature on growth rate and pigment levels were similar to general patterns observed across diverse micro-algal taxa. Results imply that dynamic behaviour associated with photophysical stress, and independent of gene regulation, might constitute one mechanism for photo-acclimation of photosynthesis.

Spin-glass dynamics

Spin glasses are magnetic systems with conflicting and random interactions between the individual spins. The dynamics of spin glasses, as of structural glasses, reflect their complexity. Both in experimental and numerical work the relaxation below the freezing temperature depends strongly on the annealing conditions (aging) and, above the freezing point, relaxation in equilibrium is slow and non-exponential, In this Forum, observed characteristics of the dynamics were summarized and the physical models proposed to explain them were outlined. (C) 1998 Elsevier Science B.V. All rights reserved.

A study of the stability of subcycling algorithms in structural dynamics

Algorithms for explicit integration of structural dynamics problems with multiple time steps (subcycling) are investigated. Only one such algorithm, due to Smolinski and Sleith has proved to be stable in a classical sense. A simplified version of this algorithm that retains its stability is presented. However, as with the original version, it can be shown to sacrifice accuracy to achieve stability. Another algorithm in use is shown to be only statistically stable, in that a probability of stability can be assigned if appropriate time step limits are observed. This probability improves rapidly with the number of degrees of freedom in a finite element model. The stability problems are shown to be a property of the central difference method itself, which is modified to give the subcycling algorithm. A related problem is shown to arise when a constraint equation in time is introduced into a time-continuous space-time finite element model. (C) 1998 Elsevier Science S.A.

Ecosystem dynamics at disturbed and undisturbed sites in North Queensland wet tropical rain forest. III. Nutrient returns to the forest floor through litterfall

The nutrient contents and accessions in litterfall over a period of 3 y are reported for undisturbed areas and at two sites disturbed by selective harvesting in tropical rain forest in North Queensland, Australia. Mean concentrations (mg g(-1) dry weight) of nutrients in litterfall ranged from 10 to 12 for nitrogen; 0.33 to 0.43 for phosphorus; 3.6 to 4.3 for potassium; 6.0 to 10.5 for calcium and 1.7 to 2.6 for magnesium. These concentrations are in the middle to lower part of the spectrum of values recorded for tropical forest. Accessions of nutrients in litterfall (kg ha(-1) y(-1)) ranged from 59 to 64 N; 1.9 to 2.4 P; 20 to 24K; 34 to 63 Ca; and 9 to 16 Mg. These rates, particularly for IN and P, are among the lowest recorded for tropical forests. There were no consistent between-site differences in total nutrient accessions in small litterfall. In terms of the contribution of litterfall to the accessions of nutrients to the forest floor, this suggests that the logged sites have recovered from the effects of selective harvesting within 25 y. Nutrient accessions at each site were distinctly seasonal, with maximum accessions occurring in the late dry season to late in the wet season. Leaf-fall accounted for the largest proportion of nutrient accessions over the study period, although at certain times accessions in both reproductive material and wood were significant. A cyclone which crossed the coast near the study sites resulted in large nutrient accessions over a short period but had little effect on the total annual accession. A comparison with previous studies of litterfall in Australian tropical rainforests indicates that nutrient return in litterfall is directly related to soil fertility.

Factors influencing biodiversity and distributional gradients in mangroves

Numerous factors affect the distribution of mangrove plants. Most mangrove species are typically dispersed by water-buoyant propagules, allowing them to lake advantage of estuarine, coastal and ocean currents both to replenish existing stands and to establish new ones. The direction they travel depends on sea currents and land barriers, but the dispersal distance depends on the time that propagules remain buoyant and viable. This is expected to differ for each species. Similarly, each species will also differ in establishment success and growth development rate, and each has tolerance limits and growth responses which are apparently unique. Such attributes are presumably responsible for the characteristic distributional ranges of each species, as each responds to the environmental, physical and biotic settings they might occupy. In practice, species are often ordered by the interplay of different factors along environmental gradients, and these may conveniently be considered at four geographic scales-global, regional, estuarine and intertidal. We believe these influencing factors act similarly around the world, and to demonstrate this point, we present examples of distributional gradients from the two global biogeographic regions, the Atlantic East Pacific and the Indo-West Pacific.

The influence of joint position on the dynamics of perception-action coupling

Six right-handed subjects performed rhythmic flexion and extension movements of the index finger in time with an auditory metronome. On each block of trials, the wrist of the response hand was placed in a extended, neutral or flexed position. In the flex-on-the-beat condition, subjects were instructed to coordinate maximum excursion in the direction of finger flexion with each beat of the metronome. In the extend-on-the-beat condition, subjects were instructed to coordinate maximum excursion in the direction of finger extension with each beat of the metronome. The frequency of the metronome was increased from 2.00 Hz to 3.75 Hz in 8 steps (8 s epochs) of 0.25 Hz. During trials prepared in the extend-on-the-beat pattern, all subjects exhibited transitions to either a flex-on-the-beat pattern or to phase wandering as the frequency of pacing was increased. The time at which these transitions occurred was reliably influenced by the position of the wrist. Four subjects exhibited qualitative departures from the flex-on-the-beat pattern at pacing frequencies that were greater than those at which the extend-on-the-beat pattern could be maintained. The lime at which these departures occurred was not influenced by the position of the wrist. These results are discussed with reference to the constraints imposed on the coordination dynamics by the intrinsic properties of the neuromuscular-skeletal system.

Relationships between blood pressure and heart rate in the saltwater crocodile Crocodylus porosus

The cardiac limb of the baroreflex loop was studied in the saltwater crocodile Crocodylus porosus, The classical pharmacological methodology using phenylephrine and sodium nitroprusside was used to trigger blood pressure changes, and the resulting alterations in heart rate were analysed quantitatively using a logistic function. Interindividual differences in resting heart rates and blood pressures were observed, but all seven animals displayed clear baroreflex responses. Atropine and sotalol greatly attenuated the response. A maximal baroreflex gain of 7.2 beats min(-1) kPa(-1) was found at a mean aortic pressure of 6.1 kPa, indicating the active role of the baroreflex in a wide pressure range encompassing hypotensive and hypertensive states. At the lowest mean aortic pressures (5.0 kPa), the synergistic role of the pulmonary-to-systemic shunt in buffering the blood pressure drop also contributes to blood pressure regulation, Pulse pressure showed a better correlation,vith heart rate and also a higher gain than mean aortic, systolic or diastolic pressures, and this is taken as an indicator of the existence of a differential control element working simultaneously with a linear proportional element.