Round Robin Study: Molecular Simulation of Thermodynamic Properties from Models with Internal Degrees of Freedom.

Thermodynamic properties are often modeled by classical force fields which describe the interactions on the atomistic scale. Molecular simulations are used for retrieving thermodynamic data from such models, and many simulation techniques and computer codes are available for that purpose. In the present round robin study, the following fundamental question is addressed: Will different user groups working with different simulation codes obtain coinciding results within the statistical uncertainty of their data? A set of 24 simple simulation tasks is defined and solved by five user groups working with eight molecular simulation codes: DL_POLY, GROMACS, IMC, LAMMPS, ms2, NAMD, Tinker, and TOWHEE.

ls1 mardyn: The Massively Parallel Molecular Dynamics Code for Large Systems.

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 that were previously out of scope for molecular dynamics simulation.

Stability of xenon oxides at high pressures.

Xenon, which is quite inert under ambient conditions, may become reactive under pressure. The possibility of the formation of stable xenon oxides and silicates in the interior of the Earth could explain the atmospheric missing xenon paradox. Using an ab initio evolutionary algorithm, we predict the existence of thermodynamically stable Xe-O compounds at high pressures (XeO, XeO(2) and XeO(3) become stable at pressures above 83, 102 and 114 GPa, respectively).

Constrained evolutionary algorithm for structure prediction of molecular crystals: methodology and applications.

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.

Ionic high-pressure form of elemental boron.

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% of boron, and it was not until 1909 that 99% 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.

Evolutionary crystal structure prediction as a tool in materials design.

Ab initio methods allow a more or less straightforward prediction of numerous physical properties of solids, but require the knowledge of their crystal structure. The evolutionary algorithm USPEX, developed by us in 2004-2006, enables reliable prediction of the stable crystal structure without relying on any experimental data. Numerous tests (mostly for systems with up to 28 atoms in the unit cell, and a few tests with up to 128 atoms/cell) showed a success rate of nearly 100%.

Crystal structure prediction using ab initio evolutionary techniques: principles and applications.

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 (nearly 100%) 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.