Molecular transport in nanopores

Simulation of the transport of methane in cylindrical silica mesopores have been performed using equilibrium and nonequilibrium molecular dynamics (NEMD) as well as dual control volume grand canonical molecular dynamics methods. It is demonstrated that all three techniques yield the same transport coefficient even in the presence of viscous flow. A modified locally averaged density model for viscous flow, combined with consideration of wall slip through a frictional condition, gives a convincing interpretation of the variation of the transport coefficient over a wide range of densities, and for various pore sizes and temperatures. Wall friction coefficients extracted from NEMD simulations are found to be consistent with momentum transfer arguments, and the approach is shown to be more meaningful than the classical slip length concept. (C) 2003 American Institute of Physics.

Hydrodynamic origin of diffusion in nanopores

We study the transport of a subcritical Lennard-Jones fluid in a cylindrical nanopore, using a combination of equilibrium and nonequilibrium as well as dual control volume grand canonical molecular dynamics methods. We show that all three techniques yield the same value of the transport coefficient for diffusely reflecting pore walls, even in the presence of viscous transport. We also demonstrate that the classical Knudsen mechanism is not manifested, and that a combination of viscous flow and momentum exchange at the pore wall governs the transport over a wide range of densities.

Bayesian inferences on the recent island colonization history by the bird Zosterops lateralis lateralis

The founding of new populations by small numbers of colonists has been considered a potentially important mechanism promoting evolutionary change in island populations. Colonizing species, such as members of the avian species complex Zosterops lateralis, have been used to support this idea. A large amount of background information on recent colonization history is available for one Zosterops subspecies, Z. lateralis lateralis, providing the opportunity to reconstruct the population dynamics of its colonization sequence. We used a Bayesian approach to combine historical and demographic information available on Z. l. lateralis with genotypic data from six microsatellite loci, and a rejection algorithm to make simultaneous inferences on the demographic parameters describing the recent colonization history of this subspecies in four southwest Pacific islands. Demographic models assuming mutation–drift equilibrium or a large number of founders were better supported than models assuming founder events for three of four recently colonized island populations. Posterior distributions of demographic parameters supported (i) a large stable effective population size of several thousands individuals with point estimates around 4000–5000; (ii) a founder event of very low intensity with a large effective number of founders around 150–200 individuals for each island in three of four islands, suggesting the colonization of those islands by one flock of large size or several flocks of average size; and (iii) a founder event of higher intensity on Norfolk Island with an effective number of founders around 20 individuals, suggesting colonization by a single flock of moderate size. Our inferences on demographic parameters, especially those on the number of founders, were relatively insensitive to the precise choice of prior distributions for microsatellite mutation processes and demographic parameters, suggesting that our analysis provides a robust description of the recent colonization history of the subspecies.

Effector dynamics of rhythmic wrist activity and its implications for (modeling) bimanual coordination

To examine the role of the effector dynamics of the wrist in the production of rhythmic motor activity, we estimated the phase shifts between the EMG and the task-related output for a rhythmic isometric torque production task and an oscillatory movement, and found a substantial difference (45-52degrees) between the two. For both tasks, the relation between EMG and task-related output (torque or displacement) was adequately reproduced with a physiologically motivated musculoskeletal model. The model simulations demonstrated the importance of the contribution of passive structures to the overall dynamics and provided an account for the observed phase shifts in the dynamic task. Additional simulations of the musculoskeletal model with added load suggested that particular changes in the phase relation between EMG and movement may follow largely from the intrinsic muscle dynamics, rather than being the result of adaptations in the neural control of joint stiffness. The implications of these results are discussed in relation to (models of) interlimb coordination in rhythmic tasks. (C) 2004 Elsevier B.V. All rights reserved.

Tractable molecular theory of transport of Lennard-Jones fluids in nanopores

We present here a tractable theory of transport of simple fluids in cylindrical nanopores, which is applicable over a wide range of densities and pore sizes. In the Henry law low-density region the theory considers the trajectories of molecules oscillating between diffuse wall collisions, while at higher densities beyond this region the contribution from viscous flow becomes significant and is included through our recent approach utilizing a local average density model. The model is validated by means of equilibrium as well nonequilibrium molecular dynamics simulations of supercritical methane transport in cylindrical silica pores over a wide range of temperature, density, and pore size. The model for the Henry law region is exact and found to yield an excellent match with simulations at all conditions, including the single-file region of very small pore size where it is shown to provide the density-independent collective transport coefficient. It is also shown that in the absence of dispersive interactions the model reduces to the classical Knudsen result, but in the presence of such interactions the latter model drastically overpredicts the transport coefficient. For larger micropores beyond the single-file region the transport coefficient is reduced at high density because of intermolecular interactions and hindrance to particle crossings leading to a large decrease in surface slip that is not well represented by the model. However, for mesopores the transport coefficient increases monotonically with density, over the range studied, and is very well predicted by the theory, though at very high density the contribution from surface slip is slightly overpredicted. It is also seen that the concept of activated diffusion, commonly associated with diffusion in small pores, is fundamentally invalid for smooth pores, and the apparent activation energy is not simply related to the minimum pore potential or the adsorption energy as generally assumed. (C) 2004 American Institute of Physics.

Modeling molecular transport in slit pores

We examine the transport of methane in microporous carbon by performing equilibrium and nonequilibrium molecular dynamics simulations over a range of pore sizes, densities, and temperatures. We interpret these simulation results using two models of the transport process. At low densities, we consider a molecular flow model, in which intermolecular interactions are neglected, and find excellent agreement between transport diffusion coefficients determined from simulation, and those predicted by the model. Simulation results indicate that the model can be applied up to fluid densities of the order to 0.1-1 nm(-3). Above these densities, we consider a slip flow model, combining hydrodynamic theory with a slip condition at the solid-fluid interface. As the diffusion coefficient at low densities can be accurately determined by the molecular flow model, we also consider a model where the slip condition is supplied by the molecular flow model. We find that both density-dependent models provide a useful means of estimating the transport coefficient that compares well with simulation. (C) 2004 American Institute of Physics.

Critical fluctuations and entanglement in the nondegenerate parametric oscillator

We present a fully quantum mechanical treatment of the nondegenerate optical parametric oscillator both below and near threshold. This is a nonequilibrium quantum system with a critical point phase transition, that is also known to exhibit strong yet easily observed squeezing and quantum entanglement. Our treatment makes use of the positive P representation and goes beyond the usual linearized theory. We compare our analytical results with numerical simulations and find excellent agreement. We also carry out a detailed comparison of our results with those obtained from stochastic electrodynamics, a theory obtained by truncating the equation of motion for the Wigner function, with a view to locating regions of agreement and disagreement between the two. We calculate commonly used measures of quantum behavior including entanglement, squeezing, and Einstein-Podolsky-Rosen (EPR) correlations as well as higher order tripartite correlations, and show how these are modified as the critical point is approached. These results are compared with those obtained using two degenerate parametric oscillators, and we find that in the near-critical region the nondegenerate oscillator has stronger EPR correlations. In general, the critical fluctuations represent an ultimate limit to the possible entanglement that can be achieved in a nondegenerate parametric oscillator.

Comparisons of diffusive and viscous contributions to transport coefficients of light gases in single-walled carbon nanotubes

We examine here the relative importance of different contributions to transport of light gases in single walled carbon nanotubes, using methane and hydrogen as examples. Transport coefficients at 298 K are determined using molecular dynamics simulation with atomistic models of the nanotube wall, from which the diffusive and viscous contributions are resolved using a recent approach that provides an explicit expression for the latter. We also exploit an exact theory for the transport of Lennard-Jones fluids at low density considering diffuse reflection at the tube wall, thereby permitting the estimation of Maxwell coefficients for the wall reflection. It is found that reflection from the carbon nanotube wall is nearly specular, as a result of which slip flow dominates, and the viscous contribution is small in comparison, even for a tube as large as 8.1 nm in diameter. The reflection coefficient for hydrogen is 3-6 times as large as that for methane in tubes of 1.36 nm diameter, indicating less specular reflection for hydrogen and greater sensitivity to atomic detail of the surface. This reconciles results showing that transport coefficients for hydrogen and methane, obtained in simulation, are comparable in tubes of this size. With increase in adsorbate density, the reflection coefficient increases, suggesting that adsorbate interactions near the wall serve to roughen the local potential energy landscape perceived by fluid molecules.

Projected Gross-Pitaevskii equation for harmonically confined Bose gases at finite temperature

We extend the projected Gross-Pitaevskii equation formalism of Davis [Phys. Rev. Lett. 87, 160402 (2001)] to the experimentally relevant case of thermal Bose gases in harmonic potentials and outline a robust and accurate numerical scheme that can efficiently simulate this system. We apply this method to investigate the equilibrium properties of the harmonically trapped three-dimensional projected Gross-Pitaevskii equation at finite temperature and consider the dependence of condensate fraction, position, and momentum distributions and density fluctuations on temperature. We apply the scheme to simulate an evaporative cooling process in which the preferential removal of high-energy particles leads to the growth of a Bose-Einstein condensate. We show that a condensate fraction can be inferred during the dynamics even in this nonequilibrium situation.

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.

Influence of adsorbate interaction on transport in confined spaces

This article provides a review of the recent theory of transport in nanopores developed in the author's laboratory. In particular the influence of fluid-solid interactions on the transport coefficient is examined, showing that such interactions reduce the value of the coefficient by almost an order of magnitude in comparison to the Knudsen theory for non-interacting systems. The activation energy and potential energy barriers for diffusion in smooth pores with a one-dimensional potential energy profile are also discussed, indicating the inadequacy of the commonly used assumption of proportionality between the activation energy and heat of adsorption or the minimum pore potential energy. A further feature affected by fluid-solid interactions is the nature of the reflection of fluid molecules colliding with a pore wall surface, varying from being nearly specular - such as in carbon nanotubes - to nearly diffuse for amorphous solids. Diffuse reflection leads to momentum loss and reduced transport coefficients. However, fluid-solid interactions do not affect the transport coefficient in the single-file diffusion regime when the surface reflection is diffuse, and the transport coefficient in this case is largely independent of the adsorbed density.

Modern developments in hypersonic wind tunnels

The development of new methods of producing hypersonic wind-tunnel flows at increasing velocities during the last few decades is reviewed with attention to airbreathing propulsion, hypervelocity aerodynamics and superorbital aerodynamics. The role of chemical reactions in these flows leads to use of a binary scaling simulation parameter, which can be related to the Reynolds number, and which demands that smaller wind tunnels require higher reservoir pressure levels for simulation of flight phenomena. The use of combustion heated vitiated wind tunnels for propulsive research is discussed, as well as the use of reflected shock tunnels for the same purpose. A flight experiment validating shock-tunnel results is described, and relevant developments in shock tunnel instrumentation are outlined. The use of shock tunnels for hypervelocity testing is reviewed, noting the role of driver gas contamination in determining test time, and presenting examples of air dissociation effects on model flows. Extending the hypervelocity testing range into the superorbital regime with useful test times is seen to be possible by use of expansion tube/tunnels with a free piston driver.

Kinetic energy conserving integrators for Gaussian thermostatted SLLOD

A new integration scheme is developed for nonequilibrium molecular dynamics simulations where the temperature is constrained by a Gaussian thermostat. The utility of the scheme is demonstrated by its application to the SLLOD algorithm which is the standard nonequilibrium molecular dynamics algorithm for studying shear flow. Unlike conventional integrators, the new integrators are constructed using operator-splitting techniques to ensure stability and that little or no drift in the kinetic energy occurs. Moreover, they require minimum computer memory and are straightforward to program. Numerical experiments show that the efficiency and stability of the new integrators compare favorably with conventional integrators such as the Runge-Kutta and Gear predictor-corrector methods. (C) 1999 American Institute of Physics. [S0021-9606(99)50125-6].