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

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

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

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

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

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

The mean free path for electron conduction in metallic fullerenes

We calculate the electrical resistivity due to electron-phonon scattering for a model of A3C60 (A= K, Rb), using an essentially exact quantum Monte-Carlo calculation. In agreement with experiment, we obtain exceptionally large metallic resistivities at large temperatures T. This illustrates that the apparent mean free path can be much shorter than the separation of the molecules. An interpretation of this result is given. The calculation also explains the linear behavior in T at small T.Comment: 4 pages, RevTeX, 3 eps figure, additional material available at http://www.mpi-stuttgart.mpg.de/docs/ANDERSEN/fullerene

Ewald method for polytropic potentials in arbitrary dimensionality

The Ewald summation technique is generalised to power-law 1/|r|^k potentials in three-, two- and one-dimensional geometries with explicit formulae for all the components of the sums. The cases of short-range, long-range and "marginal" interactions are treated separately. The jellium model, as a particular case of a charge-neutral system, is discussed and the explicit forms of the Ewald sums for such system are presented. A generalised form of the Ewald sums for a noncubic (nonsquare) simulation cell for three- (two-) dimensional geometry is obtained and its possible field of application is discussed. A procedure for the optimisation of the involved parameters in actual simulations is developed and an example of its application is presented.Comment: 41 pages, 3 figure