DL_POLY_3 is a general-purpose massively parallel molecular dynamics simulation package embedding a highly efficient set of methods and algorithms such as: Domain Decomposition (DD), Linked Cells (LC), Daresbury Advanced Fourier Transform (DAFT), Trotter derived Velocity Verlet (VV) integration and RATTLE. Written to support academic research, it has a wide range of applications and can run on a wide range of computers; from single processor workstations to multi-processor computers. The code development has placed particular emphasis on the efficient utilization of multi-processor power by optimised memory workload and distribution, which makes it possible to simulate systems of the order of tens of millions of particles and beyond. In this paper we discuss the new DL_POLY_3 design, and report on the performance, capability and scalability. We also discuss new features implemented to simulate highly non-equilibrium processes of radiation damage and analyse the structural damage during such processes. have introduced the Domain Decomposition (DD) 10 version, DL_POLY_3, to permit simulation of systems of the order of tens of millions of atoms and beyond. As we shall see in the Performance and discussion section, DL_POLY_3's inherent parallelism allows close to perfect parallelisation up to impressively high processor counts.
Scientific research relies on computer software, yet software is not always developed following practices that ensure its quality and sustainability. This manuscript does not aim to propose new software development best practices, but rather to provide simple recommendations that encourage the adoption of existing best practices. Software development best practices promote better quality software, and better quality software improves the reproducibility and reusability of research. These recommendations are designed around Open Source values, and provide practical suggestions that contribute to making research software and its source code more discoverable, reusable and transparent. This manuscript is aimed at developers, but also at organisations, projects, journals and funders that can increase the quality and sustainability of research software by encouraging the adoption of these recommendations.
Advanced potential energy surfaces are defined as theoretical models that explicitly include many-body effects that transcend the standard fixed-charge, pairwise-additive paradigm typically used in molecular simulation. However, several factors relating to their software implementation have precluded their widespread use in condensed-phase simulations: the computational cost of the theoretical models, a paucity of approximate models and algorithmic improvements that can ameliorate their cost, underdeveloped interfaces and limited dissemination in computational code bases that are widely used in the computational chemistry community, and software implementations that have not kept pace with modern high-performance computing (HPC) architectures, such as multicore CPUs and modern graphics processing units (GPUs). In this Feature article we review recent progress made in these areas, including well-defined polarization approximations and new multipole electrostatic formulations, novel methods for solving the mutual polarization equations and increasing the MD time step, combining linear scaling electronic structure methods with new QM/MM methods that account for mutual polarization between the two regions, and the greatly improved software deployment of these models and methods onto GPU and CPU hardware platforms. We have now approached an era where multipole-based polarizable force fields can be routinely used to obtain computational results comparable to state-of-the-art density functional theory while reaching sampling statistics that are acceptable when compared to that obtained from simpler fixed partial charge force fields.
The DL_POLY package provides a set of classical molecular dynamics programs that have application over a wide range of atomic and molecular systems. Written for parallel computers they offer capabilities stretching from small systems consisting of a few hundred atoms running on a single processor, up to systems of several million atoms running on massively parallel computers with thousands of processors. In this article we describe the structure of the programs and some applications.
DL_POLY_3 is a general-purpose molecular-dynamics simulation package embedding a highly efficient domain decomposition (DD) parallelization strategy. It was developed at Daresbury Laboratory under the auspices of the Engineering and Physical Sciences Research Council. Written to support academic research, it has a wide range of applications and will run on a wide range of computers; from single-processor workstations to multi-processor computers, with accent on the efficient use of multi-processor power. A new DD adaptation of the smoothed particle mesh Ewald method for calculating long-range forces in molecular simulations, incorporating a novel three-dimensional fast Fourier transform (the Daresbury Advanced Fourier Transform), makes it possible to simulate systems of the order of one million particles and beyond. DL_POLY_3 structure, functionality, performance and availability are described in this feature paper.
Understanding and predicting a material's performance in response to high-energy radiation damage, as well as designing future materials to be used in intense radiation environments, requires the knowledge of the structure, morphology and amount of radiation-induced structural changes 1-5 . We report the results of molecular dynamics simulations of high-energy radiation damage in iron in the range 0.2-0.5 MeV. We analyze and quantify the nature of collision cascades both at the global and local scale. We find that the structure of high-energy collision cascades becomes increasingly continuous as opposed to showing sub-cascade branching reported previously. At the local length scale, we find large defect clusters and novel small vacancy and interstitial clusters. These features form the basis for physical models aimed at understanding the effects of high energy radiation damage in structural materials.
We use molecular-dynamics simulations to study processes related to resistance to amorphization by radiation damage. We simulate high-energy radiation events in SiO 2 , GeO 2 , TiO 2 , Al 2 O 3 , and MgO, and find that simulation results match the experiments. We discuss the difference between the damage that the structures along this series can support. We find that for the same material, activation barriers for damage recovery can strongly depend on the degree of structural damage. We observe that the effect of resistance to amorphization is primarily governed by the relaxation processes at the time scales of several picoseconds. On this time scale, we observe two distinct relaxation processes, reversible elastic deformation around the radiation cascade and recovery of the in-cascade damage of high topological disorder. Finally, we discuss how resistance to amorphization is related to interatomic interactions and to the nature of the chemical bond.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.