Prospects for release-node quantum Monte Carlo.

We perform release-node quantum Monte Carlo simulations on the first row diatomic molecules in order to assess how accurately their ground-state energies can be obtained. An analysis of the fermion-boson energy difference is shown to be strongly dependent on the nuclear charge, Z, which in turn determines the growth of variance of the release-node energy. It is possible to use maximum entropy analysis to extrapolate to ground-state energies only for the low Z elements. For the higher Z dimers beyond boron, the error growth is too large to allow accurate data for long enough imaginary times. Within the limit of our statistics we were able to estimate, in atomic units, the ground-state energy of Li(2) (-14.9947(1)), Be(2) (-29.3367(7)), and B(2)(-49.410(2)).

A semiclassical study of the thermal conductivity of low temperature liquids.

The conventional classical energy current auto-correlation function has been extended into a quantum mechanical version and then approximated by the linearized semiclassical initial value representation approach. Comparison of the thermal conductivity to simulation results shows that about 15% quantum correction to the classical molecular dynamics results for liquid neon are quantitatively predicted. For liquid para-hydrogen the quantum effects are sufficiently large that the linearized semiclassical approach is only 20% accurate, while for both liquid He(4) and He(3) the thermal conductivity disagrees by a factor of 2, although exchange effects appear to play a minor role.

Atomistic methods in fluid simulation.

Atomistic methods, such as molecular dynamics and direct simulation Monte Carlo, constitute a powerful and growing set of techniques for fluid-dynamics simulation. The more fundamental nature of such methods, which exhibit nonlinear transport effects and small-scale fluctuations, extends their modelling accuracy to a significantly wider range of scales and regimes than the more traditional Navier-Stokes-based continuum fluid-simulation techniques. In this paper, we describe the current state of the art in atomistic fluid simulation, from both a theoretical and a computational standpoint, and outline the advantages and limitations of such methods. In addition, we present an overview of some recent atomistic-simulation results on fluid instabilities and on the physical scaling of atomistic techniques. Finally, we suggest possible avenues of future research in the field.

Tethered DNA dynamics in shear flow.

We study the cyclic dynamics of a single polymer tethered to a hard wall in shear flow using Brownian dynamics, the lattice Boltzmann method, and a recent stochastic event-driven molecular dynamics algorithm. We focus on the dynamics of the free end (last bead) of the tethered chain and we examine the cross-correlation function and power spectral density of the chain extensions in the flow and gradient directions as a function of chain length N and dimensionless shear rate Wi. Extensive simulation results suggest a classical fluctuation-dissipation stochastic process and question the existence of periodicity of the cyclic dynamics, as previously claimed. We support our numerical findings with a simple analytical calculation for a harmonic dimer in shear flow.

Scaling of atomistic fluid dynamics simulations.

We have performed a series of large-scale atomistic simulations of the Rayleigh-Taylor instability including up to 5.7 x 10(9) particles and spanning time and length scales of up to 170 ns and 45 microm , respectively. The results suggest that atomistic fluid dynamics simulations exhibit the same scaling as solutions of the continuum Navier-Stokes equations. Furthermore, a comparison with macroscopic Rayleigh-Taylor experiments suggests that the results of such atomistic simulations can, in fact, be scaled up to macroscopic dimensions, even for complex, nonstationary flows.

Stochastic hard-sphere dynamics for hydrodynamics of nonideal fluids.

A novel stochastic fluid model is proposed with a nonideal structure factor consistent with compressibility, and adjustable transport coefficients. This stochastic hard-sphere dynamics (SHSD) algorithm is a modification of the direct simulation Monte Carlo algorithm and has several computational advantages over event-driven hard-sphere molecular dynamics. Surprisingly, SHSD results in an equation of state and a pair correlation function identical to that of a deterministic Hamiltonian system of penetrable spheres interacting with linear core pair potentials. The fluctuating hydrodynamic behavior of the SHSD fluid is verified for the Brownian motion of a nanoparticle suspended in a compressible solvent.

The importance of fluctuations in fluid mixing.

A ubiquitous example of fluid mixing is the Rayleigh-Taylor instability, in which a heavy fluid initially sits atop a light fluid in a gravitational field. The subsequent development of the unstable interface between the two fluids is marked by several stages. At first, each interface mode grows exponentially with time before transitioning to a nonlinear regime characterized by more complex hydrodynamic mixing. Unfortunately, traditional continuum modeling of this process has generally been in poor agreement with experiment. Here, we indicate that the natural, random fluctuations of the flow field present in any fluid, which are neglected in continuum models, can lead to qualitatively and quantitatively better agreement with experiment. We performed billion-particle atomistic simulations and magnetic levitation experiments with unprecedented control of initial interface conditions. A comparison between our simulations and experiments reveals good agreement in terms of the growth rate of the mixing front as well as the new observation of droplet breakup at later times. These results improve our understanding of many fluid processes, including interface phenomena that occur, for example, in supernovae, the detachment of droplets from a faucet, and ink jet printing. Such instabilities are also relevant to the possible energy source of inertial confinement fusion, in which a millimeter-sized capsule is imploded to initiate nuclear fusion reactions between deuterium and tritium. Our results suggest that the applicability of continuum models would be greatly enhanced by explicitly including the effects of random fluctuations.

Nanohydrodynamics simulations: an atomistic view of the Rayleigh-Taylor instability.

Nanohydrodynamics simulations, hydrodynamics on the nanometer and nanosecond scale by molecular dynamics simulations for up to 100 million particles, are performed on the latest generation of supercomputers. Such simulations exhibit Rayleigh-Taylor instability, the mixing of a heavy fluid on top of a light in the presence of a gravitational field, initiated by thermal fluctuations at the interface, leading to the chaotic regime in the long-time evolution of the mixing process. The early-time behavior is in general agreement with linear analysis of continuum theory (Navier-Stokes), and the late-time behavior agrees quantitatively with experimental observations. Nanohydrodynamics provides insights into the turbulent mixing process that are inaccessible to either continuum calculations or to experiment.