An introduction to classical molecular dynamics simulation for experimental scattering users

TL;DR

This paper introduces an open educational resource to help scattering experiment users understand classical molecular dynamics simulations, emphasizing accessibility and practical application in data analysis.

Contribution

It provides a freely accessible, interactive web-based educational tool and open source software to teach classical molecular dynamics to scattering users.

Findings

Educational resource enhances understanding of molecular dynamics

Open source tools facilitate learning and customization

Improves integration of simulations in scattering data analysis

Abstract

Classical molecular dynamics simulations are a common component of multi-modal analyses from scattering measurements, such as small-angle scattering and diffraction. Users of these experimental techniques often have no formal training in the theory and practice of molecular dynamics simulation, leading to the possibility of these simulations being treated as a "black box" analysis technique. In this article, we describe an open educational resource (OER) designed to introduce classical molecular dynamics to users of scattering methods. This resource is available as a series of interactive web pages, which can be easily accessed by students, and as an open source software repository, which can be freely copied, modified, and redistributed by educators. The topic covered in this OER includes classical atomistic modelling, parameterising interatomic potentials, molecular dynamics simulations, typical sources of error, and some of the approaches to using simulations in the analysis of scattering data.

Figures

Full text

Usage

Electronic Supplementary Information (ESI) available: All analysis/plotting scripts and figure files, allowing for a fully reproducible, and automated, analysis workflow for the work presented is available at https://github.com/arm61/sim_and_scat_paper (DOI: 10.5281/zenodo.2556826) under a CC BY-SA 4.0 license.

I Introduction

The use of molecular dynamics simulations to help analyse and interpret experimental data from small-angle scattering and diffraction studies has grown significantly over the past ten years Pan et al. (2012); Boldon et al. (2015); Hub (2018); Ivanović et al. (2018); East et al. (2016); Wall et al. (2014); Wall (2018); Satoh et al. (2015), with the percentage of small-angle scattering publications that mention molecular dynamics reaching more than $20\text{\,}\mathrm{\char 37\relax}$ in 2018 (Fig. 1). Users of scattering and diffraction techniques often have backgrounds in experimental science and may have received little formal training in the theory or practice of computational modelling. This can lead to molecular dynamics simulations being used as a “black box”, without understanding the underlying methodologies, or considering possible sources of error. The use of molecular simulations without any technical understanding can lead to unintended but severe errors. To help support researchers use of molecular dynamics simulation in their analysis of scattering data while reducing this risk of modelling errors, a number of software tools, such as WAXSiS and SASSIE Chen and Hub (2014); Knight and Hub (2015); Perkins et al. (2016), have been developed that present easy-to-use, graphical, web-based user interfaces.

A complementary approach is to organise educational activities, such as lectures or workshops, tailored to introduce molecular simulation techniques to audiences of scattering and diffraction users. One example is the annual ISIS Neutron Training Course, which includes a module titled “An Introduction to Molecular Dynamics for Neutron Scattering”. This module covers the fundamentals of classical molecular dynamics simulation, presents applications of these methods in neutron science, and gives students practical hands-on experience with the SASSIE software package Perkins et al. (2016).

While lectures and workshops are an effective tool for education and training, participation can be limited due to difficulties attending in person (due to location or cost) or physical limits on student numbers. An alternative educational strategy gaining popularity within scientific and engineering communities is the publication of technology-enhanced open educational resources (OERs). These are courses, lectures, or learning resources published online that are freely available for use by anyone. In addition to their broad accessibility, these resources have permissive “open” copyright licences that allow their use in the “5R activities”: retain, reuse, revise, remix, and redistribute Wiley . Publishing a resource as an OER increases its reach and impact, because others may use it in their own teaching not only in its original form, but are free to modify, and redistribute, the material to better suit their aims. The Jupyter Notebook framework Kluyver et al. (2016) has become a popular platform for OERs that teach computational skills, because it allows authors to include instructional text, images, and other media, alongside executable, editable code, in an example of “literate programming” Knuth (1984). This format encourages students to directly interact with code examples by running, editing, and rerunning these within the source document Barba et al. (2017), supporting exploratory experiential learning Papert (1993).

Here, we present an online, open-source, interactive learning resource written to introduce members of the scattering and diffraction community to molecular dynamics simulations. This OER comprises five lessons that introduce classical molecular dynamics methods and show how these can be used to assist in the analysis of experimental scattering data by the calculation of a simulated scattering profile from the molecular dynamics simulation. We use the open-source Python library pylj McCluskey et al. (2018, 2019a) to provide simple, but computationally authentic, examples of simulations, that demonstrate visually and programmatically the conceptual relationships between simulation and scattering techniques. In this article, we discuss the structure of the resource and describe how a student may get the most from the resource 111In this work the student refers to anyone working through the OER, regardless of career position..

II Assumed prior knowledge

The OER, entitled “The Interaction Between Simulation and Scattering” makes use of the Python programming language to include interactive examples of mathematical and algorithmic content. To be able to fully use this resource, therefore, requires some knowledge of, or willingness to learn, Python. We have, however, attempted to design the resource in such a fashion that an in-depth knowledge of Python is not a prerequisite. It is anticipated that the users of this resource have some experience of at least undergraduate level chemistry or physics, and a commensurate level of mathematical understanding. This background knowledge will be particularly important for the sections dealing with the functional form and chemical basis of classical interaction potentials.

III Resource construction

The resource is available in two main formats. First, as a series of web pages, hosted at http://pythoninchemistry.org/sim_and_scat. Secondly as the source-code repository used to build these webpages, hosted at https://github.com/pythoninchemistry/sim_and_scat McCluskey et al. (2019b). The source content consists of a set of Jupyter Notebooks and Markdown files, which are automatically compiled using the jupyter-book tool Lau and Holdgraf (2019) to generate the web version. This system allows the resulting webpages to include text, equations, and figures, which we use to describe key concepts and to explain details of algorithms, as well as Python code blocks, which we use to provide specific examples. The web pages have Thebelab and BinderHub integrations Ragan-Kelley et al. (2019a, b); Jupyter and others (2018), which allow students to launch interactive versions of these web pages that allow execution and editing of the included Python code. We believe that the ability to read the resource as an “interactive document” will help students to develop a deeper understanding of the materials than they would from a traditional static document. The resource is provided under a CC-BY license Commons (2019), while the jupyter-book software is shared under an MIT license Source (2019), both of which are open and highly permissive. This allows readers to reuse or remix the material to enhance their own educational platform, and for secondary authors to contribute to improving the source material.

IV Resource outline

The resource follows a simple outline that introduces key aspects of molecular dynamics simulations.

IV.1 Home

The welcome page introduces the resource and gives the student information about how the resource may be used, including details of the Thebelab and BinderHub integrations. This page also contains details about the permitted use and sharing of the content of the resource, and about the resource license. This page also includes a list of authors and contributors.

IV.2 Classical methods

After a brief introduction, this section introduces concepts related to classical simulation methods. This includes the use of interatomic potential functions, alongside some examples, such as the Lennard-Jones and Buckingham potential models Lennard-Jones (1924); Buckingham (1938). We then introduce the problem of parameterising a potential model, including the use of higher accuracy quantum mechanical calculations to do so. The presence of off-the-shelf, general potential models are discussed; with the caveat that they may still require system specific optimisation. Finally, we mention mixing rules; again discussing the possible problems that a user may encounter if applying these blindly to specific systems.

IV.3 Molecular dynamics

Having introduced the concept of classical interatomic potentials, we then discuss how these are used in molecular dynamics simulations. We work through how a one-dimensional NVE (constant number, volume, and energy) molecular dynamics simulation may be built, using the Velocity-Verlet algorithm and the Lennard-Jones potential model Swope et al. (1982); Lennard-Jones (1924). The Velocity-Verlet algorithm is introduced in terms of Newton’s laws of motion and the generalised equations of motion. Finally, we discuss a range of key factors that can affect molecular dynamics simulations, including choice of simulation ensembles, setting the range cut-off for an interatomic potential, and the use of periodic boundary conditions.

IV.4 pylj and the interaction with scattering

The final aspect of the resource covers using molecular dynamics simulations to understand scattering profiles. This is presented as a practical example, using the open-source pylj package McCluskey et al. (2018, 2019a). We demonstrate a two-dimensional molecular dynamics simulation of argon particles interacting through a Lennard-Jones potential. The student is first shown a working pylj simulation and invited to interact with the simulation and the custom plotting functionality of pylj. The concept of a radial distribution function (RDF) is then shown, and the students are given the opportunity to run some pylj simulations with the RDF being output alongside the simulation window. Next we present the Debye equation Debye (1915) and show how it may be used to calculate scattering data from a simulation. The student is invited to observe the effect of simulation temperature on the resulting scattering profile. We finish by discussing alternative, faster, algorithms for calculating scattering profiles, such as the Fibonacci Sequence or Golden Vectors methods Svergun (1994); Watson and Curtis (2013).

IV.5 “Real” simulation and scattering

Having shown the development of a scattering profile from an idealised system, we then direct the student to a popular resource for the GROMACS Berendsen et al. (1995) molecular dynamics software. This resource gives a quick introduction to using GROMACS to simulate a lyzosyme molecule in buffer Lemkul , the student may then use their own simulated trajectory or one that can be downloaded from the OER. We show how the system may be visualised, introduce the MDAnalysis Python package for the analysis of molecular dynamics trajectories Michaud-Agrawal et al. (2011); Gowers et al. (2016), and show the scattering profile developed from the lysozyme simulation compared with experimental data Franke et al. (2015). The module finishes by pointing the student to resources to more easily resolve scattering data from molecular simulation, such as SASSIE and CRYSOL Perkins et al. (2016); Svergun et al. (1995). The focus of this module is to introduce simulation methodologies to users of scattering to aid their understanding, not to derive the exact mechanics of the calculation of scattering from simulation. Resources for this purpose already exist and have well developed tutorials, so it is not necessary to recreate such software here Perkins et al. (2016); Svergun et al. (1995).

V Future outlook

The aim of this resource is benefit students and researchers interested in combining molecular dynamics analysis with scattering and diffraction experiments. Our intention in publishing this material as an OER is to promote interest from other parties towards reuse, and remixing the material for their own educational purposes. We plan to implement this material within training at the ISIS Neutron and Muon Source as well as Diamond Light Source. Finally, it is hoped, via student and community feedback, that the implementation, content, and pedagogical approach of this resource can be increased over time.

VI Author contributions

The open education resource was developed, and the manuscript written, by A. R. M. with input from all authors.

Acknowledgements.

A. R. M. is grateful to the University of Bath and Diamond Light Source for co-funding a studentship (Studentship No. STU0149). B. J. M. acknowledges support from the Royal Society (Grant No. UF130329).