A Simple Aromatic
Structure of liquid benzene from isotopic neutron data
80 minutes
Summary
You might always be asking the question “How much does quantity X differ between the equilibrated simulation at the reference forcefield parameters, and the refined simulation at the empirical parameters?”. This is a very important question, and this example will show how to checkpoint simulation data for comparison later on. It will also introduce the utility of size factor scaling when dealing with molecular systems containing rings.
You will set up and equilibrate a small liquid benzene simulation before calculating a few properties of interest, and then refine the potential against experimental data and recalculate the properties of interest to see how much they change. The data are three neutron scattering datasets measured on the NIMROD instrument at the ISIS Pulsed Neutron and Muon Source in 2014. We’ll be referencing the paper by Headen et al. along the way - “Structure of π-π Interactions in Aromatic Liquids”, T. F. Headen, C. A. Howard, N. T. Skipper, M. A. Wilkinson, D. T. Bowron and A. K. Soper, J. Am. Chem. Soc. 132, 5735 (2010).
Preparation
Download the example data as a zip or tar.gz and unpack it if you haven’t done so already. Alternatively, download the data files listed below and save them to a location of your choice.
- C6H6: C6H6.mint01
- C6H6:C6D6 (50:50): 5050.mint01
- C6D6: C6D6.mint01
File ⇨ New
File ⇨ Save As…
Let's beginSave your input file under a sensible name in same directory as the example data