Monte Carlo and molecular dynamics computer simulations are now able to provide detailed information concerning the structure, dynamics, and thermodynamics of pore fluids in 2:1 clays. This article will discuss interparticle interaction potentials currently available for atomistic simulations of clay-water systems, and will describe how computational techniques can be applied to modelling of cl...

We use numerical simulation to examine the possibility of a reversible liquid-liquid transition in supercooled water and related systems. In particular, for two atomistic models of water, we have computed free energies as functions of multiple order parameters, where one is density and another distinguishes crystal from liquid. For a range of temperatures and pressures, separate free energy bas...

Results showed that the way the fluid-wall interaction is modeled in molecular dynamics simulations has a strong effect on the resultant simulation of liquid injection into a gas. It was found that modeling the wall as individual atoms (atomistic model) interacting with the fluid resulted in the fluid remaining in an injection tube absent any pressure or body force to force it into a gaseous re...

Our atomistic molecular dynamics simulations reveal the existence of bound and free water molecules in the hydration layer of an aqueous micelle. The bound water molecules can be either singly or doubly hydrogen bonded to the polar head group on the surface of the micelle. The ratio of bound to free water is found to be approximately equal to 9:1 at 300 K. email:[email protected] email:bbagchi@...

Very few works exist to date on development of a consistent energy-based coupling of atomistic and continuum models of materials in more than one dimension. The difficulty in constructing such a coupling consists in defining a coupled energy whose minimizers are free from uncontrollable errors on the atomistic/continuum interface. In this paper a consistent coupling in three dimensions is propo...

This paper extends our earlier studies of free energy functions of density and crystalline order parameters for models of supercooled water, which allows us to examine the possibility of two distinct metastable liquid phases [D. T. Limmer and D. Chandler, J. Chem. Phys. 135, 134503 (2011) and preprint arXiv:1107.0337 (2011)]. Low-temperature reversible free energy surfaces of several different ...

Gay-Berne anisotropic potential has been widely used to evaluate the non-bonded interactions between coarse-grained particles being described as elliptical rigid bodies. In this paper, we are presenting a coarse-grained model for twenty kinds of amino acids and proteins, based on the anisotropic Gay-Berne and point electric multipole (EMP) potentials. We demonstrate that the anisotropic coarse-...

The formation and motion of lattice defects such as cracks, dislocations, or grain boundaries, occurs when the lattice configuration loses stability, that is, when an eigenvalue of the Hessian of the lattice energy functional becomes negative. When the atomistic energy is approximated by a hybrid energy that couples atomistic and continuum models, the accuracy of the approximation can only be g...