On the use of molecular dynamics simulation to calculate X-ray thermal diffuse scattering from molecular crystals

Molecular dynamics simulation of organic crystals: introducing the CLP-dyncry environment

The CLP-dyncry molecular dynamics (MD) program suite and force field environment is introduced and validated with its ad hoc features for the treatment of organic crystalline matter. The package, stemming from a preliminary implementation on organic liquids (Gavezzotti & Lo Presti, 2019), includes modules for the preliminary generation of molecular force field files from ab initio derived force constants, and for the preparation of crystalline simulation boxes from general crystallographic information, including Cambridge Structural Database CIFs. The intermolecular potential is the atom–atom Coulomb–London–Pauli force field, well tested as calibrated on sublimation enthalpies of organic crystals. These products are then submitted to a main MD module that drives the time integration and produces dynamic information in the form of coordinate and energy trajectories, which are in turn processed by several kinds of crystal-oriented analytic modules. The whole setup is tested on a variety of bulk crystals of rigid, non-rigid and hydrogen-bonded compounds for the reproduction of radial distribution functions and of crystal-specific collective orientational variables against X-ray data. In a series of parallel tests, some advantages of a dedicated program as opposed to software more oriented to biomolecular simulation (Gromacs) are highlighted. The different and improved view of crystal packing that results from joining static structural information from X-ray analysis with dynamic upgrades is also pointed out. The package is available for free distribution with I/O examples and Fortran source codes.

Evaluation of General and Tailor Made Force Fields via X-ray Thermal Diffuse Scattering Using Molecular Dynamics and Monte Carlo Simulations of Crystalline Aspirin

We have performed a comparison of the experimental thermal diffuse scattering (TDS) from crystalline Aspirin (form I) to that calculated from molecular dynamics (MD) simulations based on a variety of general force fields and a tailor-made force field (TMFF). A comparison is also made with Monte Carlo (MC) simulations which use a "harmonic network" approach to describe the intermolecular interactions. These comparisons were based on the hypothesis that TDS could be a useful experimental data in validation of such simulation parameter sets, especially when calculations of dynamical properties (e.g., thermodynamic free energies) from molecular crystals are concerned. Currently such a validation of force field parameters against experimental data is often limited to calculation of specific physical properties, e.g., absolute lattice energies usually at 0 K or heat capacity measurements. TDS harvested from in-house or synchrotron experiments comprises highly detailed structural information representative of the dynamical motions of the crystal lattice. Thus, TDS is a well-suited experimental data-driven means of cross validating theoretical approaches targeted at understanding dynamical properties of crystals. We found from the results of our investigation that the TMFF and COMPASS (from the commercial software "Materials Studio") parameter sets gave the best agreement with experiment. From our homologous MC simulation analysis we are able to show that force constants associated with the molecular torsion angles are likely to be a strong contributing factor for the apparent reason why these aforementioned force fields performed better.

A method to perform modulated structure studies using the programZMC

The Monte Carlo simulation toolkitZMCprovides a general method for Monte Carlo simulation of disordered molecular displacements, orientations and site occupations. It is part of a suite of programs that allows convenient calculation of three-dimensional reciprocal space diffraction intensities. Such intensities are useful for extensive structure interpretations from diffraction experiments. Some physical properties can also be calculated from the simulation. This article reports the implementation of an integratedZMCmodule that facilitates the construction of modulated molecular crystals and allows the corresponding calculation of modulation satellites. Use of the module is explained and examples illustrated. The method is ideally suited for novice treatment of modulated structure features and also as a basic platform to explore complex structure/properties relationships using atomistic simulation techniques.