Evaluation of Coulomb potential in a triclinic cell with periodic boundary conditions

Lekner and Sperb's work on the evaluation of Coulomb energy and forces under periodic boundary conditions is generalized that makes it possible to use a triclinic unit cell in simulations in 3D rather than just an orthorhombic cell. The expressions obtained are in a similar form as previously obtained by Lekner and Sperb for the especial case of orthorhombic cell

Effect of the Berendsen thermostat on dynamical properties of water

The effect of the Berendsen thermostat on the dynamical properties of bulk SPC/E water is tested by generating power spectra associated with fluctuations in various observables. The Berendsen thermostat is found to be very effective in preserving temporal correlations in fluctuations of tagged particle quantities over a very wide range of frequencies. Even correlations in fluctuations of global properties, such as the total potential energy, are well-preserved for time periods shorter than the thermostat time constant. Deviations in dynamical behaviour from the microcanonical limit do not, however, always decrease smoothly with increasing values of the thermostat time constant but may be somewhat larger for some intermediate values of $\tau_B$, specially in the supercooled regime, which are similar to time scales for slow relaxation processes in bulk water.Comment: 21 pages, 5 figures, To be published in Mol. Phy

Expressions for forces and torques in molecular simulations using rigid bodies

Expressions for intermolecular forces and torques, derived from pair potentials between rigid non-spherical units, are presented. The aim is to give compact and clear expressions, which are easily generalised, and which minimise the risk of error in writing molecular dynamics simulation programs. It is anticipated that these expressions will be useful in the simulation of liquid crystalline systems, and in coarse-grained modelling of macromolecules

Selective-pivot sampling of radial distribution functions in asymmetric liquid mixtures

We present a Monte Carlo algorithm for selectively sampling radial distribution functions and effective interaction potentials in asymmetric liquid mixtures. We demonstrate its efficiency for hard-sphere mixtures, and for model systems with more general interactions, and compare our simulations with several analytical approximations. For interaction potentials containing a hard-sphere contribution, the algorithm yields the contact value of the radial distribution function.Comment: 5 pages, 5 figure

Atomistic studies of transformation pathways and energetics in plutonium

One of the most challenging problems in understanding the structural phase transformations in Pu is to determine the energetically favored, continuous atomic pathways from one crystal symmetry to another. This problem involves enumerating candidate pathways and studying their energetics to garner insight into instabilities and energy barriers. The purpose of this work is to investigate the energetics of two transformation pathways for the delta to alpha' transformation in Pu that were recently proposed [Lookman et al., Phys. Rev. Lett. 100:145504, 2008] on the basis of symmetry. These pathways require the presence of either an intermediate hexagonal closed-packed (hcp) structure or a simple hexagonal (sh) structure. A subgroup of the parent fcc and the intermediate hexagonal structure, which has trigonal symmetry, facilitates the transformation to the intermediate hcp or sh structure. Phonons then break the translational symmetry from the intermediate hcp or sh structure to the final monoclinic symmetry of the alpha' structure. We perform simulations using the modified embedded atom method (MEAM) for Pu to investigate these candidate pathways. Our main conclusion is that the path via hcp is energetically favored and the volume change for both pathways essentially occurs in the second step of the transformation, i.e. from the intermediate sh or hcp to the monoclinic structure. Our work also highlights the deficiency of the current state-of-the-art MEAM potential in capturing the anisotropy associated with the lower symmetry monoclinic structure.Comment: 12 pages, 5 figures, accepted for publication in Philos. Ma

Frame dragging with optical vortices

General Relativistic calculations in the linear regime have been made for electromagnetic beams of radiation known as optical vortices. These exotic beams of light carry a physical quantity known as optical orbital angular momentum (OAM). It is found that when a massive spinning neutral particle is placed along the optical axis, a phenomenon known as inertial frame dragging occurs. Our results are compared with those found previously for a ring laser and an order of magnitude estimate of the laser intensity needed for a precession frequency of 1 Hz is given for these "steady" beams of light.Comment: 13 pages, 2 figure

On the interplay between sedimentation and phase separation phenomena in two-dimensional colloidal fluids

Colloidal particles that are confined to an interface effectively form a two-dimensional fluid. We examine the dynamics of such colloids when they are subject to a constant external force, which drives them in a particular direction over the surface. Such a situation occurs, for example, for colloidal particles that have settled to the bottom of their container, when the container is tilted at an angle, so that they `sediment' to the lower edge of the surface. We focus in particular on the case when there are attractive forces between the colloids which causes them to phase separate into regions of high density and low density and we study the influence of this phase separation on the sedimentation process. We model the colloids as Brownian particles and use both Brownian dynamics computer simulations and dynamical density functional theory (DDFT) to obtain the time evolution of the ensemble average one-body density profiles of the colloids. We consider situations where the external potential varies only in one direction so that the ensemble average density profiles vary only in this direction. We solve the DDFT in one-dimension, by assuming that the density profile only varies in one direction. However, we also solve the DDFT in two-dimensions, allowing the fluid density profile to vary in both the $x$- and $y$-directions. We find that in certain situations the two-dimensional DDFT is clearly superior to its one-dimensional counterpart when compared with the simulations and we discuss this issue.Comment: 17 pages, 10 figures, submitted to Molecular Physic

Path integral Monte Carlo simulations for rigid rotors and their application to water

In this work the path integral formulation for rigid rotors, proposed by M\"user and Berne [Phys. Rev. Lett. {\bf 77}, 2638 (1996)], is described in detail. It is shown how this formulation can be used to perform Monte Carlo simulations of water. Our numerical results show that whereas some properties of water can be accurately reproduced using classical simulations with an empirical potential which, implicitly, includes quantum effects, other properties can only be described quantitatively when quantum effects are explicitly incorporated. In particular, quantum effects are extremely relevant when it comes to describing the equation of state of the ice phases at low temperatures, the structure of the ices at low temperatures, and the heat capacity of both liquid water and the ice phases. They also play a minor role in the relative stability of the ice phases.Comment: to appear in Molecular Physics (2011

Void Growth in BCC Metals Simulated with Molecular Dynamics using the Finnis-Sinclair Potential

The process of fracture in ductile metals involves the nucleation, growth, and linking of voids. This process takes place both at the low rates involved in typical engineering applications and at the high rates associated with dynamic fracture processes such as spallation. Here we study the growth of a void in a single crystal at high rates using molecular dynamics (MD) based on Finnis-Sinclair interatomic potentials for the body-centred cubic (bcc) metals V, Nb, Mo, Ta, and W. The use of the Finnis-Sinclair potential enables the study of plasticity associated with void growth at the atomic level at room temperature and strain rates from 10^9/s down to 10^6/s and systems as large as 128 million atoms. The atomistic systems are observed to undergo a transition from twinning at the higher end of this range to dislocation flow at the lower end. We analyze the simulations for the specific mechanisms of plasticity associated with void growth as dislocation loops are punched out to accommodate the growing void. We also analyse the process of nucleation and growth of voids in simulations of nanocrystalline Ta expanding at different strain rates. We comment on differences in the plasticity associated with void growth in the bcc metals compared to earlier studies in face-centred cubic (fcc) metals.Comment: 24 pages, 12 figure