The goal of this lab is to become more comfortable with running simulations in GROMACS and to learn how to create your own input files for basic simulations. The system of study is a single acetone molecule (CH3)2CO solvated in SPC/E water.
The lab procedure assumes independent work with GROMACS manual. Refer to Lab 2b for general information about running GROMACS.
We will use the acetone model from the paper “Kirkwood-Buff derived force field for mixtures of acetone and water” by S. Weerasinghe and P. E. Smith (J. Chem. Phys. 118, p 10663, 2003). This model corresponds to the united atom type: the methyl groups are represented by a single interaction site.
Write a topology file (topol.top) for a system containing only one acetone molecule and no water. Table I in the paper is useful for specifying geometry of acetone. Check the structure of the topology file from previous examples. Look through Chapter 5: Topologies and read carefully Section 5.7: File Formats of GROMACS 5.0 manual in PDF format here . Find the combination rules that were used in the paper and invoke those rules in the
[default] section, see Table II in particular.
Create a coordinate file (conf.gro) for a single acetone molecule contained in a simulation box of approximately 1.8 x 1.8 x 1.8 nm3. You don't need to have a perfect structure.
IMPORTANT: You can get the geometry of acetone from the Computational Chemistry Comparison and Benchmark Database.
WARNING: The coordinate file
conf.gro has a fixed format. The columns are in fixed positions. Refer to section 5.7.6.
Check that you can actually run a short simulation of this isolated acetone molecule, 10 steps are enough. Remember that to run the simulation you must load the GROMACS module first, then execute grompp, and finally mdrun. Read the output of these two programs, and pay attention to error and warning messages.
Check the trajectory from the simulation with VMD or PyMOL
You may continue with the following only after steps 2 through 5 were completed successfully.
(A) Include [moleculetype] SOL for the SPC/E water model in the topology file.
GROMACS provides an auxiliary script genbox to fill the simulation box with water molecules. Check
gmx genbox -h and carefully read the program help and manual. Use
gmx genbox to solvate your initial simulation box that has the solitary acetone molecule.
(B) Equilibrate a system to T=300 K and P=1 atm. For that, you need to equilibration stages: a 50 ps of NVT, and 100 ps of NPT (with stochastic dynamics, integrator=sd).
(C) Include the harmonic bonds, angle, and improper dihedral in topology file. You may use information in Table 1. Equilibrate the system and perform calculations in NPT ensemble.