-
ls1 mardyn: The massively parallel molecular dynamics code for large systems
Authors:
Christoph Niethammer,
Stefan Becker,
Martin Bernreuther,
Martin Buchholz,
Wolfgang Eckhardt,
Alexander Heinecke,
Stephan Werth,
Hans-Joachim Bungartz,
Colin W. Glass,
Hans Hasse,
Jadran Vrabec,
Martin Horsch
Abstract:
The molecular dynamics simulation code ls1 mardyn is presented. It is a highly scalable code, optimized for massively parallel execution on supercomputing architectures, and currently holds the world record for the largest molecular simulation with over four trillion particles. It enables the application of pair potentials to length and time scales which were previously out of scope for molecular…
▽ More
The molecular dynamics simulation code ls1 mardyn is presented. It is a highly scalable code, optimized for massively parallel execution on supercomputing architectures, and currently holds the world record for the largest molecular simulation with over four trillion particles. It enables the application of pair potentials to length and time scales which were previously out of scope for molecular dynamics simulation. With an efficient dynamic load balancing scheme, it delivers high scalability even for challenging heterogeneous configurations. Presently, multi-center rigid potential models based on Lennard-Jones sites, point charges and higher-order polarities are supported. Due to its modular design, ls1 mardyn can be extended to new physical models, methods, and algorithms, allowing future users to tailor it to suit their respective needs. Possible applications include scenarios with complex geometries, e.g. for fluids at interfaces, as well as non-equilibrium molecular dynamics simulation of heat and mass transfer.
△ Less
Submitted 20 August, 2014;
originally announced August 2014.
-
Stability of xenon oxides at high pressures
Authors:
Qiang Zhu,
Daniel Y. Jung,
Artem R. Oganov,
Colin W. Glass,
Carlo Gatti,
Andriy O. Lyakhov
Abstract:
Xenon, which is quite inert under ambient conditions, may become reactive under pressure. The possibility of formation of stable xenon oxides and silicates in the interior of the Earth could explain the atmospheric missing xenon paradox. Using the ab initio evolutionary algorithm, we predict the thermodynamical stabilization of Xe-O compounds at high pressures (XeO, XeO2 and XeO3 at pressures abov…
▽ More
Xenon, which is quite inert under ambient conditions, may become reactive under pressure. The possibility of formation of stable xenon oxides and silicates in the interior of the Earth could explain the atmospheric missing xenon paradox. Using the ab initio evolutionary algorithm, we predict the thermodynamical stabilization of Xe-O compounds at high pressures (XeO, XeO2 and XeO3 at pressures above 83, 102 and 114 GPa, respectively). Our calculations indicate large charge transfer in these oxides, suggesting that large electronegativity difference and pressure are the key factors favoring the formation of xenon compounds. Xenon compounds in the Earth's mantle, however, cannot directly explain the missing xenon paradox: xenon oxides are unstable in equilibrium with metallic iron in the Earth's lower mantle, while xenon silicates are predicted to spontaneously decompose at all mantle pressures (<136 GPa). This does not preclude Xe atoms from being retained in defects of mantle silicates and oxides.
△ Less
Submitted 29 November, 2012; v1 submitted 28 November, 2012;
originally announced November 2012.
-
Constrained evolutionary algorithm for structure prediction of molecular crystals: methodology and applications
Authors:
Qiang Zhu,
Artem R. Oganov,
Colin W. Glass,
Harold T. Stokes
Abstract:
Evolutionary crystal structure prediction proved to be a powerful approach for studying a wide range of materials. Here, we present a specifically designed algorithm for the prediction of the structure of complex crystals consisting of well-defined molecular units. The main feature of this new approach is that each unit is treated as a whole body, which drastically reduces the search space and imp…
▽ More
Evolutionary crystal structure prediction proved to be a powerful approach for studying a wide range of materials. Here, we present a specifically designed algorithm for the prediction of the structure of complex crystals consisting of well-defined molecular units. The main feature of this new approach is that each unit is treated as a whole body, which drastically reduces the search space and improves the efficiency, but necessitates the introduction of new variation operators described here. To increase diversity of the population of structures, the initial population and part($\scriptsize{\sim}$20%) of the new generations are generated using space group symmetry combined with random cell parameters and random positions and orientations of molecular units. We illustrate the efficiency and reliability of this approach by number of tests (ice, ammonia, carbon dioxide, methane, benzene, glycine and butane-1,4-diammonium dibromide). This approach easily predicts the crystal structure of methane \emph{A} containing 21 methane molecules (105 atoms) per unit cell. We demonstrate that this new approach has also a high potential for the study of complex inorganic crystals on the examples of a complex hydrogen storage material Mg(BH$_4$)$_2$ and elemental boron.
△ Less
Submitted 17 May, 2012; v1 submitted 20 April, 2012;
originally announced April 2012.
-
Static and dynamic properties of curved vapour-liquid interfaces by massively parallel molecular dynamics simulation
Authors:
Martin T. Horsch,
Svetlana K. Miroshnichenko,
Jadran Vrabec,
Colin W. Glass,
Christoph Niethammer,
Martin F. Bernreuther,
Erich A. Müller,
George Jackson
Abstract:
Curved fluid interfaces are investigated on the nanometre length scale by molecular dynamics simulation. Thereby, droplets surrounded by a metastable vapour phase are stabilized in the canonical ensemble. Analogous simulations are conducted for cylindrical menisci separating vapour and liquid phases under confinement in planar nanopores. Regarding the emergence of nanodroplets during nucleation, a…
▽ More
Curved fluid interfaces are investigated on the nanometre length scale by molecular dynamics simulation. Thereby, droplets surrounded by a metastable vapour phase are stabilized in the canonical ensemble. Analogous simulations are conducted for cylindrical menisci separating vapour and liquid phases under confinement in planar nanopores. Regarding the emergence of nanodroplets during nucleation, a non-equilibrium phenomenon, both the non-steady dynamics of condensation processes and stationary quantities related to supersaturated vapours are considered. Results for the truncated and shifted Lennard-Jones fluid and for mixtures of quadrupolar fluids confirm the applicability of the capillarity approximation and the classical nucleation theory.
△ Less
Submitted 20 October, 2011;
originally announced October 2011.
-
Ionic high-pressure form of elemental boron
Authors:
Artem R. Oganov,
Jiuhua Chen,
Carlo Gatti,
Yanzhang Ma,
Yanming Ma,
Colin W. Glass,
Zhenxian Liu,
Tony Yu,
Oleksandr O. Kurakevych,
Vladimir L. Solozhenko
Abstract:
Boron is an element of fascinating chemical complexity. Controversies have shrouded this element since its discovery was announced in 1808: the new 'element' turned out to be a compound containing less than 60-70 percent of boron, and it was not until 1909 that 99-percent pure boron was obtained. And although we now know of at least 16 polymorphs, the stable phase of boron is not yet experimenta…
▽ More
Boron is an element of fascinating chemical complexity. Controversies have shrouded this element since its discovery was announced in 1808: the new 'element' turned out to be a compound containing less than 60-70 percent of boron, and it was not until 1909 that 99-percent pure boron was obtained. And although we now know of at least 16 polymorphs, the stable phase of boron is not yet experimentally established even at ambient conditions. Boron's complexities arise from frustration: situated between metals and insulators in the periodic table, boron has only three valence electrons, which would favour metallicity, but they are sufficiently localized that insulating states emerge. However, this subtle balance between metallic and insulating states is easily shifted by pressure, temperature and impurities. Here we report the results of high-pressure experiments and ab initio evolutionary crystal structure predictions that explore the structural stability of boron under pressure and, strikingly, reveal a partially ionic high-pressure boron phase. This new phase is stable between 19 and 89 GPa, can be quenched to ambient conditions, and has a hitherto unknown structure (space group Pnnm, 28 atoms in the unit cell) consisting of icosahedral B12 clusters and B2 pairs in a NaCl-type arrangement. We find that the ionicity of the phase affects its electronic bandgap, infrared adsorption and dielectric constants, and that it arises from the different electronic properties of the B2 pairs and B12 clusters and the resultant charge transfer between them.
△ Less
Submitted 18 November, 2009; v1 submitted 16 November, 2009;
originally announced November 2009.
-
Crystal structure prediction using ab initio evolutionary techniques: principles and applications
Authors:
A. R. Oganov,
C. W. Glass
Abstract:
We have developed an efficient and reliable methodology for crystal structure prediction, merging ab initio total-energy calculations and a specifically devised evolutionary algorithm. This method allows one to predict the most stable crystal structure and a number of low-energy metastable structures for a given compound at any P-T conditions without requiring any experimental input. Extremely h…
▽ More
We have developed an efficient and reliable methodology for crystal structure prediction, merging ab initio total-energy calculations and a specifically devised evolutionary algorithm. This method allows one to predict the most stable crystal structure and a number of low-energy metastable structures for a given compound at any P-T conditions without requiring any experimental input. Extremely high success rate has been observed in a few tens of tests done so far, including ionic, covalent, metallic, and molecular structures with up to 40 atoms in the unit cell. We have been able to resolve some important problems in high-pressure crystallography and report a number of new high-pressure crystal structures. Physical reasons for the success of this methodology are discussed.
△ Less
Submitted 16 November, 2009;
originally announced November 2009.
-
New high-pressure form of boron is significantly ionic
Authors:
Artem R. Oganov,
Jiuhua Chen,
Carlo Gatti,
Yanzhang Ma,
Yanming Ma,
Colin W. Glass,
Zhenxian Liu,
Tony Yu,
Oleksandr O. Kurakevych,
Vladimir L. Solozhenko
Abstract:
The comment of Dubrovinskaia et al. is scientifically flawed. The high-pressure form of boron, discovered by Oganov et al., is indeed new and its bonding has a significant ionic character, as demonstrated in Ref. 1.
The comment of Dubrovinskaia et al. is scientifically flawed. The high-pressure form of boron, discovered by Oganov et al., is indeed new and its bonding has a significant ionic character, as demonstrated in Ref. 1.
△ Less
Submitted 27 December, 2010; v1 submitted 4 August, 2009;
originally announced August 2009.