Search

Data and metadata documentation requirements for explainable-AI-ready (XAIR) models and data in physics-based simulation technology are discussed by analysing different perspectives from the literature on two core aspects: First, the scope of the simulation; this category is taken to include subject matter, the objective with which the simulation is conducted, and the object of reference, i.e., the simulated physical system or process. Second, the artefacts that need to be documented in order to make data and models XAIR, and modelling and simulation workflows explainable; two CEN workshop agreements, MODA and ModGra, are compared for this purpose. As a result, minimum requirements for an ontologization of the scope of simulation artefacts are formulated, and the object-objective abstractness diagram is proposed as a tool for visualizing the landscape of use cases for physics-based simulation.

It is shown that the residual entropy (entropy minus that of the ideal gas at the same temperature and density) is mostly synonymous with the independent variable of density scaling, identifying a direct link between these two approaches. The two-body residual entropy is demonstrated to not be a suitable surrogate for the total residual entropy in the gas phase. The residual entropy and the effective hardness of interaction (itself a derivative at constant residual entropy) are studied for the Lennard-Jones monomer and dimer as well as a range of rigid molecular models for carbon dioxide. It is observed that the density scaling exponent is related to the two-body interactions.

The Riemann problem is one of the basic building blocks for numerical methods in computational fluid mechanics. Nonetheless, there are still open questions and gaps in theory and modelling for situations with complex thermodynamic behavior. In this series, we compare numerical solutions of the macroscopic flow equations with molecular dynamics simulation data. To enable molecular dynamics for sufficiently large scales in time and space, we selected the truncated and shifted Lennard-Jones potential for which also highly accurate equations of state are available. A comparison of a two-phase Riemann problem is shown, which involves a liquid and a vapor phase, with an undergoing phase transition. The loss of hyperbolicity allows for the occurrence of anomalous wave structures. We successfully compare the molecular dynamics data with two macroscopic numerical solutions obtained by either assuming local phase equilibrium or by imposing a kinetic relation and allowing for metastable states.

The MolMod database is presented, which is openly accessible at http://molmod.boltzmann-zuse.de/ and contains presently intermolecular force fields for over 150 pure fluids. It was developed and is maintained by the Boltzmann-Zuse Society for Computational Molecular Engineering (BZS). The set of molecular models in the MolMod database provides a coherent framework for molecular simulations of fluids. The molecular models in the MolMod database consist of Lennard-Jones interaction sites, point charges, and point dipoles and quadrupoles, which can be equivalently represented by multiple point charges. The force fields can be exported as input files for the simulation programs ms2 and ls1 mardyn, Gromacs, and LAMMPS. To characterise the semantics associated with the numerical database content, a force-field nomenclature is introduced that can also be used in other contexts in materials modelling at the atomistic and mesoscopic levels. The models of the pure substances that are included in the data base were generally optimised such as to yield good representations of experimental data of the vapour-liquid equilibrium with a focus on the vapour pressure and the saturated liquid density. In many cases, the models also yield good predictions of caloric, transport, and interfacial properties of the pure fluids. For all models, references to the original works in which they were developed are provided. The models can be used straightforwardly for predictions of properties of fluid mixtures using established combination rules. Input errors are a major source of errors in simulations. The MolMod database contributes to reducing such errors.

In manuscript arXiv:1703.08719 [cond-mat.soft], it was claimed that the well-known deduction of Tolman's law is not rigorous, since Tolman's argument implies that two different definitions of the surface tension, called $\gamma$ and $\bar\gamma$ in the manuscript, coincide. This claim is retracted as it can be shown by free-energy minimization that $\gamma = \bar\gamma$ indeed holds for the Laplace radius. Joachim Gro\ss, Philipp Rehner, Carlos Vega, \O{}ivind Wilhelmsen, and the anonymous reviewers of The Journal of Chemical Physics contributed to finding the mistake in the manuscript., Comment: Manuscript arXiv:1703.08719 [cond-mat.soft] contains claims that were shown to be incorrect; the manuscript and the claims are retracted

A slab based long range correction approach for multi-site Lennard-Jones models is presented for systems with a planar film geometry that is based on the work by Janecek, J. Phys. Chem. B 110: 6264 (2006). It is efficient because it relies on a center-of-mass cutoff scheme and scales in terms of numerics almost perfectly with the molecule number. For validation, a series of simulations with the two-center Lennard-Jones model fluid, carbon dioxide and cyclohexane is carried out. The results of the present approach, a site-based long range correction and simulations without any long range correction are compared with respect to the saturated liquid density and the surface tension. The present simulation results exhibit only a weak dependence on the cutoff radius, indicating a high accuracy of the implemented long range correction.

A new version release (2.0) of the molecular simulation tool ms2 [S. Deublein et al., Comput. Phys. Commun. 182 (2011) 2350] is presented. Version 2.0 of ms2 features a hybrid parallelization based on MPI and OpenMP for molecular dynamics simulation to achieve higher scalability. Furthermore, the formalism by Lustig [R. Lustig, Mol. Phys. 110 (2012) 3041] is implemented, allowing for a systematic sampling of Massieu potential derivatives in a single simulation run. Moreover, the Green-Kubo formalism is extended for the sampling of the electric conductivity and the residence time. To remove the restriction of the preceding version to electro-neutral molecules, Ewald summation is implemented to consider ionic long range interactions. Finally, the sampling of the radial distribution function is added.

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.

The present level of development of molecular force field methods is assessed from the point of view of simulation-based engineering, outlining the immediate perspective for further development and highlighting the newly emerging discipline of Computational Molecular Engineering (CME) which makes basic research in soft matter physics fruitful for industrial applications. Within the coming decade, major breakthroughs can be reached if a research focus is placed on processes at interfaces, combining aspects where an increase in the accessible length and time scales due to massively parallel high-performance computing will lead to particularly significant improvements.

Massively-parallel molecular dynamics simulation is applied to systems containing electrolytes, vapour-liquid interfaces, and biomolecules in contact with water-oil interfaces. Novel molecular models of alkali halide salts are presented and employed for the simulation of electrolytes in aqueous solution. The enzymatically catalysed hydroxylation of oleic acid is investigated by molecular dynamics simulation taking the internal degrees of freedom of the macromolecules into account. Thereby, Ewald summation methods are used to compute the long range electrostatic interactions. In systems with a phase boundary, the dispersive interaction, which is modelled by the Lennard-Jones potential here, has a more significant long range contribution than in homogeneous systems. This effect is accounted for by implementing the Janecek cutoff correction scheme. On this basis, the HPC infrastructure at the Steinbuch Centre for Computing was accessed and efficiently used, yielding new insights on the molecular systems under consideration.