Step 9a - Centre-of-Mass RDF
Examining the centre-of-mass RDF between molecules is often useful since it gives a general picture of the arrangement between molecule. A related quantity is the coordination number, i.e. the number of molecules within solvation shells that typically form around molecules in the liquid. Here we’ll set up calculation of the centre-of-mass RDF and the associated coordination number for the first solvation shell. Dissolve provides a predefined layer to calculate both at once, however first we need to define a suitable site - the centre of mass of the water molecule.
Define the Centre-of-Mass Site
If you don’t already have a site COM
representing the centre-of-mass of the molecule, perform the following steps:
Open the Water and go to the Analysis Sites section
Click
to create a new site
In the species viewer to the right, click-drag with the left mouse button and draw a box around all three atoms of the water molecule to select them.
Click the
button in the toolbar of the viewer to set the currently-selected atoms as defining the origin
Select the Weight by mass option in the Origin Atom(s) section
Double-click on the new site in the list and change its name to
COM
You’ll see the site represented in the species viewer as a 3D crosshair indicating the coordinates of the origin, surrounded by a small cube.
Create an Analysis Layer
We’ll now create one of Dissolve’s predefined layers to give us the processing modules we need:
Layer ⇨ Create ⇨ Analysis ⇨ RDF & Coordination Number
The new layer contains the following modules:
Module | Purpose |
---|---|
CalculateRDF | Calculates the radial distribution function between two specified sites |
CalculateCN | Calculates a coordination number by summing histogram data calculated by a CalculateRDF module |
We’ll need to set up both of these modules to calculate exactly what we need. First, CalculateRDF :
Select the CalculateRDF module to display its options
Open the Control settings group
Press the button for SiteA and choose the
COM
site
Select the
COM
site for SiteB as well
Enable the ExcludeSameMolecule option
Since we are calculating the RDF of a site around itself, the ExcludeSameMolecule option prevents consideration of the same site with itself (which would give a distance of zero and result in a large undesirable spike at the beginning of the RDF). The distance range over which to calculate the RDF can be set in the Control settings group, but the defaults (0 - 10 Å with a step size of 0.05 Å) are fine for what we need.
Now the CalculateCN module:
Select the CalculateCN module to display its options
Open the Ranges settings group
Set the maximum for RangeA to 3.4 Å
The CalculateCN module calculates up to three separate coordination numbers over different distance regions of the target RDF, but we will focus on the first coordination shell in this example. The target RDF is taken from a CalculateRDF module, specified by the SourceRDF keyword (this will have been set automatically to the CalculateRDF module in the current layer).