A complete analytic theory for structure and dynamics of populations and communities spanning wide ranges in body size

The prediction and management of ecosystem responses to global environmental change would profit from a clearer understanding of the mechanisms determining the structure and dynamics of ecological communities. The analytic theory presented here develops a causally closed picture for the mechanisms controlling community and population size structure, in particular community size spectra, and their dynamic responses to perturbations, with emphasis on marine ecosystems. Important implications are summarised in non-technical form. These include the identification of three different responses of community size spectra to size-specific pressures (of which one is the classical trophic cascade), an explanation for the observed slow recovery of fish communities from exploitation, and clarification of the mechanism controlling predation mortality rates. The theory builds on a community model that describes trophic interactions among size-structured populations and explicitly represents the full life cycles of species. An approximate time-dependent analytic solution of the model is obtained by coarse graining over maturation body sizes to obtain a simple description of the model steady state, linearising near the steady state, and then eliminating intraspecific size structure by means of the quasi-neutral approximation. The result is a convolution equation for trophic interactions among species of different maturation body sizes, which is solved analytically using a novel technique based on a multiscale expansion.

Relativistic theory for localized electrostatic excitations in degenerate electron-ion plasmas

A self-consistent relativistic two-fluid model is proposed for electron-ion plasma dynamics. A one-dimensional geometry is adopted. Electrons are treated as a relativistically degenerate fluid, governed by an appropriate equation of state. The ion fluid is also allowed to be relativistic, but is cold, nondegenerate, and subject only to an electrostatic potential. Exact stationary-profile solutions are sought, at the ionic scale, via the Sagdeev pseudopotential method. The analysis provides the pulse existence region, in terms of characteristic relativistic parameters, associated with the (ultrahigh) particle density.

Testing Wavefunction Collapse Models using Parametric Heating of a Trapped Nanosphere

We propose a mechanism for testing the theory of collapse models such as continuous spontaneous localization (CSL) by examining the parametric heating rate of a trapped nanosphere. The random localizations of the center-of-mass for a given particle predicted by the CSL model can be understood as a stochastic force embodying a source of heating for the nanosphere. We show that by utilising a Paul trap to levitate the particle and optical cooling, it is possible to reduce environmental decoher- ence to such a level that CSL dominates the dynamics and contributes the main source of heating. We show that this approach allows measurements to be made on the timescale of seconds, and that the free parameter λcsl which characterises the model ought to be testable to values as low as 10^{−12} Hz.

Spin and spatial dynamics in electron-impact scattering off S-wave He using R-matrix with Time-Dependence theory

R-matrix with time-dependence theory is applied to electron-impact ionisation processes for He in the S-wave model. Cross sections for electron-impact excitation, ionisation and ionisation with excitation for impact energies between 25 and 225 eV are in excellent agreement with benchmark cross sections. Ultra-fast dynamics induced by a scattering event is observed through time-dependent signatures associated with autoionisation from doubly excited states. Further insight into dynamics can be obtained through examination of the spin components of the time-dependent wavefunction.

Particle Image Velocimetry Investigation of a finite amplitude pressure wave

Particle image velocimetry is used to study the motion of gas within a duct subject to the passage of a finite amplitude pressure wave. The wave is representative of the pressure waves found in the exhaust systems of internal combustion engines. Gas particles are accelerated from stationary to 150 m/s and then back to stationary in 8 ms. It is demonstrated that gas particles at the head of the wave travel at the same velocity across the duct cross section at a given point in time. Towards the tail of the wave viscous effects are plainly evident causing the flow profile to tend towards parabolic. However, the instantaneous mean particle velocity across the section is shown to match well with the velocity calculated from a corresponding measured pressure history using 1D gas dynamic theory. The measured pressure history at a point in the duct was acquired using a high speed pressure transducer of the type typically used for engine research in intake and exhaust systems. It is demonstrated that these are unable to follow the rapid changes in pressure accurately and that they are prone to resonate under certain circumstances.

Free energy and molecular dynamics calculations for the cubic-tetragonal phase transition in zirconia

The high-temperature cubic-tetragonal phase transition of pure stoichiometric zirconia is studied by molecular dynamics (MD) simulations and within the framework of the Landau theory of phase transformations. The interatomic forces are calculated using an empirical, self-consistent, orthogonal tight-binding model, which includes atomic polarizabilities up to the quadrupolar level. A first set of standard MD calculations shows that, on increasing temperature, one particular vibrational frequency softens. The temperature evolution of the free-energy surfaces around the phase transition is then studied with a second set of calculations. These combine the thermodynamic integration technique with constrained MD simulations. The results seem to support the thesis of a second-order phase transition but with unusual, very anharmonic behavior above the transition temperature.

Unexpected higher-order effects in charged particle impact ionization at high energies

Most of the experimental and theoretical studies of electron-impact ionization of atoms, referred to as (e, 2e), have concentrated on the scattering plane. The assumption has been that all the important physical effects will be observable in the scattering plane. However, very recently it has been shown that, for C6+-helium ionization, experiment and theory are in nice agreement in the scattering plane and in very bad agreement out of the scattering plane. This lack of agreement between experiment and theory has been explained in terms of higher-order scattering effects between the projectile and target ion. We have examined electron-impact ionization of magnesium and have observed similar higher-order effects. The results of the electron-impact ionization of magnesium indicate the possible deficiencies in the calculation of fully differential cross sections in previous heavy particle ionization work.

Ab initio molecular dynamics simulation of a room temperature ionic liquid

Ab initio molecular dynamics simulations have been performed for the first time on the room-temperature organic ionic liquid dimethyl imidazolium chloride [DMIM][Cl] using density functional theory. The aim is to compare the local liquid structure with both that obtained from two different classical force fields and from neutron scattering experiments. The local structure around the cation shows significant differences compared to both the classical calculations and the neutron results. In particular, and unlike in the gas-phase ion pair, chloride ions tend to be located near a ring C-H proton in a position suggesting hydrogen bonding. The results are used to suggest ways in which the classical potentials may be improved.

Restricted rotation dynamics and far-infrared spectra of liquid water governed by elastic interactions

A semi-phenomenological model describing wideband dielectric and far-infrared spectra of liquid water was proposed recently by the same authors [J. Mol. Struct. 606 (2002) 9], where a small dipole-moment component changing harmonically with time determines a weak absorption band (termed here the R-band) centred at the wavenumber v similar to 200 cm(-1). In the present work, a rough molecular theory of the R-band based on the concept of elastic interactions is given. Stretching and bending of hydrogen bonds cause restricted rotation (RR) of a polar water molecule in terms of a dimer comprising the H- bonded molecules. Analytical expression for the RR frequency nu(str) is derived as a function of the RR amplitude, geometrical parameters and force constants. The density g(nu(str)) of frequency distribution is shown to be centred in the R-band. The spectrum of the dipolar auto-correlation function calculated for this structural-dynamical model is found. A composite model comprising two intermolecular potentials is proposed, which yields for water a good description of the experimental wideband (from 0 to 1000 cm(- 1)) spectra of complex permittivity and of absorption coefficient. The presented interpretation of these spectra is based on a concept that water presents a two-component solution, with components differing by the types of molecular rotation. (C) 2003 Elsevier B.V. All rights reserved.

Semiclassical treatment of Wannier's theory when the exponent diverges

We consider a non-standard application of the Wannier model. A physical example is the single ionization of a hydrogenic beryllium ion with a fully stripped beryllium ion, where the ratio of the charge of the third particle to the charges of the escaping particles is 1/4; we investigate the single ionization by an electron of an atom comprising an electron and a nucleus of charge 1/4. An infinite exponent is obtained suggesting that this process is not tractable within the Wannier model. A modified version of Crothers' uniform semiclassical wavefunction for the outgoing particles has been adopted, since the Wannier exponents and are infinite for an effective charge of Z = 1/4. We use Bessel functions to describe the Peterkop functions u and u and derive a new turning point ?. Since u is well behaved at infinity, there exists only the singularity in u at infinity, thus we employ a one- (rather than two-) dimensional change of dependent variable, ensuring that a uniform solution is obtained that avoids semiclassical breakdown on the Wannier ridge. The regularized final-state asymptotic wavefunction is employed, along with a continuum-distorted-wave approximation for the initial-state wavefunction to obtain total cross sections on an absolute scale. © 2006 IOP Publishing Ltd.

Kinetics of phase transformations in a model with metastable fluid-fluid separation: A molecular dynamics study

By molecular dynamics (MD) simulations we study the crystallization process in a model system whose particles interact by a spherical pair potential with a narrow and deep attractive well adjacent to a hard repulsive core. The phase diagram of the model displays a solid-fluid equilibrium, with a metastable fluid-fluid separation. Our computations are restricted to fairly small systems (from 2592 to 10368 particles) and cover long simulation times, with constant energy trajectories extending up to 76x10(6) MD steps. By progressively reducing the system temperature below the solid-fluid line, we first observe the metastable fluid-fluid separation, occurring readily and almost reversibly upon crossing the corresponding line in the phase diagram. The nucleation of the crystal phase takes place when the system is in the two-fluid metastable region. Analysis of the temperature dependence of the nucleation time allows us to estimate directly the nucleation free energy barrier. The results are compared with the predictions of classical nucleation theory. The critical nucleus is identified, and its structure is found to be predominantly fcc. Following nucleation, the solid phase grows steadily across the system, incorporating a large number of localized and extended defects. We discuss the relaxation processes taking place both during and after the crystallization stage. The relevance of our simulation for the kinetics of protein crystallization under normal experimental conditions is discussed. (C) 2002 American Institute of Physics.