The MOLDY short-range molecular dynamics package

Abstract

We describe a parallelised version of the MOLDY molecular dynamics program. This Fortran code is aimed at systems which may be described by shortrange potentials and specifically those which may be addressed with the embedded atom method. This includes a wide range of transition metals and alloys. MOLDY provides a range of options in terms of the molecular dynamics ensemble used and the boundary conditions which may be applied. A number of standard potentials are provided, and the modular structure of the code allows new potentials to be added easily. The code is parallelised using OpenMP and can therefore be run on shared memory systems, including modern multicore processors. Particular attention is paid to the updates required in the main force loop, where synchronisation is often required in OpenMP implementations of molecular dynamics. We examine the performance of the parallel code in detail and give some examples of applications to realistic problems, including the dynamic compression of copper and carbon migration in an iron-carbon alloy. Moldy addresses the problem of many atoms (of order 10 6 ) interacting via a classical interatomic potential on a timescale of microseconds. It is designed for problems where statistics must be gathered over a number of equivalent runs, such as measuring thermodynamic properities, diffusion, radiation damage, fracture, twinning deformation, nucleation and growth of phase transitions, sputtering etc. In the vast majority of materials, the interactions are non-pairwise, and the code must be able to deal with many-body forces. Solution method: Molecular dynamics involves integrating Newton's equations of motion. MOLDY uses verlet (for good energy conservation) or predictor-corrector (for accurate trajectories) algorithms. It is parallelised using open MP. It also includes a static minimisation routine to find the lowest energy structure. Boundary conditions for surfaces, clusters, grain boundaries, thermostat (Nose), barostat (ParrinelloRahman), and externally applied strain are provided. The initial configuration can be either a repeated unit cell or have all atoms given explictly. Initial velocities are generated internally, but it is also possible to specify the velocity of a particular atom. A wide range of interatomic force models are implemented, including embedded atom, morse or Lennard Jones. Thus the program is especially well suited to calculations of metals. KeywordsRestrictions: The code is designed for short-ranged potentials, and there is no Ewald sum. Thus for long range interactions where all particles interact with all others, the order-N scaling will fail. Different interatomic potential forms require recompilation of the code. Additional comments: There is a set of associated open-source analysis software for postprocessing and visualisation. This includes local crystal structure recognition and identification of topological defects. Running time: A set of test modules for running time are provided. The code scales as order N. The parallelisation ...

Parallelization of Molecular-Dynamics Simulations Using Tasks

“…A similar idea is described in Ref. [3]. While the implementation in the latter article uses a locking/unlocking mechanism to restrict access to particles in the same cell, the cell method uses a conditional task schedule to achieve the same effect without explicit locking.…”

Study of Liquid Lead Corrosion of Fast Nuclear Reactorand Its Mitigation by Using Molecular Dynamics Method

“…2 (11) We can see the graph is linearly enough with correlation coefficient 2 cor R = 0.977. By using (10) and (11) …”

“…After complete MD simulation then for every temperature we calculated all properties: MSD, D and D(T) by using (8)- (10). Fig.…”

“…In this present work we study the corrosion by employing the molecular dynamics (MD) simulation method that is a powerful technique for material computation [9]- [10]. The application of MD simulation method to study the corrosion of liquid metal is relatively new and an interesting research topic.…”

Thermal conductivity of argon: Modeling and comparison with experimental data

“…The heat transfer is because of the anharmonicity of the potential in the crystal, which is known as phonon conductivity, studied by the Monte Carlo method [2]. The MD simulation is for now a well-established method to calculate several transport and thermal properties in fluids [3]. The Lennard-Jones potential has been considered a pivotal model [4], which takes the molecules as hard spheres containing two main characteristics: repulsive at short distance and attractive at intermediate and long distances [5,6].…”