MD Algorithms and the Properties of Water¶
In this exercise the student will explore the behaviour of several molecular dynamics algorithms by examining their effects on the properties of a simple model for water: the SPC, or Simple Point Charge, model.
The SPC model for water, despite its simplicity, finds use in many simulations, particularly of biological systems. It is a rigid model, with three charged sites and Lennard-Jones interactions. Its properties (and limitations) are well known. The purpose of this exercise is to determine what effects, if any, the choice of molecular dynamics algorithm has on the statistical properties of this system.
DL_POLY offers a wide selection of MD algorithms including:
Fixed volume, fixed energy (NVE);
Fixed volume, constant temperature (NVT-Berendsen);
Fixed volume, constant temperature (NVT-Andersen);
Fixed volume, constant temperature (NVT-Langevin);
Fixed volume, constant temperature (NVT-Evans);
Fixed volume, constant temperature (NVT-Nosé-Hoover);
Fixed volume, constant temperature (NVT-Gentle-Stochastic);
Constant pressure, constant temperature (NPT-Berendsen);
Constant pressure, constant temperature (NPT-Nosé-Hoover);
Constant pressure, constant temperature (NPT-Martyna-Tuckerman-Klein);
Constant pressure, constant temperature (NPT-Langevin);
The DL_POLY_4 manual contains the complete list.
Where Notice that we have omitted the constant stress algorithms from the list - can you think why? DL_POLY_4 also offers two possible ways of handling the rigid body dynamics of the water molecules:
Rigid-body quaternion dynamics;
Constraint bond dynamics.
All of these algorithms are described in the DL_POLY_4 manual and in the references cited therein. You can select any of these in the CONTROL file panel of the GUI, which creates the run CONTROL file.
Our objective in this exercise is to perform a single state point simulation of water, using as many different algorithms as possible (up to 12 in all) and see what differences are manifest in the results.
Download Exercise4.tar.xz to $DL_POLY/data directory and unpack:
cd $DL_POLY/data wget https://ccp5.gitlab.io/dlpoly/Exercise4.tar.xz
Now go to the $DL_POLY/execute directory and start up the GUI:
Copy the contents of the subdirectory Exercise4 into the execute subdirectory. To do this using the GUI, select from the main menu Execute > Store/Fetch Data and in the Data Archive window type Exercise4 in the Fetch box and click Fetch. You will obtain the files CONTROL, FIELD and CONFIG. Proceed as follows.
Using the FIELD file without modification, simulate the water system using the NVE ensemble. Run, see , for at least 1000 timesteps. Examine the final results in terms of the mean thermodynamic properties (system energy and so on), structure (as revealed by the RDF) and mean-square displacement. This will be the reference system.
Repeat the study under fixed volume, constant temperature conditions, using one of the thermostated algorithms mentioned above. Set the system temperature to equal the mean temperature obtained from step 1.
Repeat the study under constant pressure, constant temperature conditions using one (or more) of the above mentioned algorithms. Set the system temperature and pressure to correspond to the mean values obtained in step 1.
Compare closely the results you obtain. Note any differences you feel are significant and try to work out what may have caused them. Broadly speaking, one should get similar results, but we are attempting to gain some feel for any systematic differences between the methods.
If you have the stamina, you may edit the FIELD file and switch the quaternion and constraint algorithms and repeat the study!
In performing this exercise, you will need to find adequate values for the ‘relaxation parameters’ of the thermostat and barostat. Clearly the choice of these represents additional degrees of freedom in the model, which may or may not also have an impact on the results. What do you think?
If all this is too boring for you, you may consider picking one particular algorithm and making a study of Hydrogen bonding in this system. First of all, decide for yourself how H-bonding may be identified and how you may quantify the degree of H-bonding as a function of temperature (say).