You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
Copy file name to clipboardExpand all lines: _build/classical_methods/intro.md
+2-2Lines changed: 2 additions & 2 deletions
Display the source diff
Display the rich diff
Original file line number
Diff line number
Diff line change
@@ -13,7 +13,7 @@ comment: "***PROGRAMMATICALLY GENERATED, DO NOT EDIT. SEE ORIGINAL FILES IN /con
13
13
# Classical methods
14
14
15
15
**Classical methods** is a phrase used to describe techniques that make use of a potential model (sometimes called a force-field) to simulated chemical systems.
16
-
These can be molecular dynamics (which we will cover in this module), Monte Carlo, Langevin dynamics, etc.
16
+
These can be molecular dynamics (which we will cover in this resource), Monte Carlo, Langevin dynamics, etc.
17
17
18
18
In order to simulate a **real** chemical system, it is necessary to model the electrons and their interactions.
19
19
This is achieved by using quantum mechanical calculations, where the energy of a chemical system is calculated by finding an approximate solution to the Schrödinger equation.
@@ -24,4 +24,4 @@ With this assumption envoked, it is possible to find the potential energy using
24
24
The first is some ground-state quantum mechanical method (e.g. density functional theory), however as mentioned above these are limited in system size.
25
25
The alternative involves modelling the electron distributions with some mathematical function to determine the potential energy.
26
26
These **potential models** are usually faster to calculate than the quantum mechanical method and therefore may be used on larger systems.
27
-
However, this **simplification** does have a drawback in that the correct potential models must be determined for each system to accurately determine the energy.
27
+
However, this **simplification** does have a drawback in that the correct potential models must be determined for each system to accurately determine the energy.
Copy file name to clipboardExpand all lines: _build/important_considerations/ensembles.md
+2-2Lines changed: 2 additions & 2 deletions
Display the source diff
Display the rich diff
Original file line number
Diff line number
Diff line change
@@ -12,7 +12,7 @@ comment: "***PROGRAMMATICALLY GENERATED, DO NOT EDIT. SEE ORIGINAL FILES IN /con
12
12
---
13
13
## Ensembles
14
14
15
-
The molecular dynamics algorithm outlined in the previous lesson makes use of the NVE ensemble (also known as the microcanonical ensemble), where the number of particles (N), volume of the system (V), and energy of the system (E) are all **kept constant**.
15
+
The molecular dynamics algorithm outlined in the previously makes use of the NVE ensemble (also known as the microcanonical ensemble), where the number of particles (N), volume of the system (V), and energy of the system (E) are all **kept constant**.
16
16
This is not the only, or the most accurate, ensemble that exists, there is also other such as:
17
17
- NVT (canonical): number of particles (N), volume of system (V), temperature of the simulation (T)
18
18
- NPT (isothermal-isobaric): number of particles (N), pressure of system (P), temperature of the simulation (T)
@@ -31,7 +31,7 @@ This means that the velocities of the particles may be rescaled by the following
where $T_{\text{target}}$ is the target temperature for the themostat, and $\bar{T}$ is the average simulation temperature.
34
-
pylj [[1,2](#references)], the software that you shall use in the next lesson uses this method for producing an NVT simulation, using the `heat_bath` function.
34
+
pylj [[1,2](#references)], the software that you shall use later uses this method for producing an NVT simulation, using the `heat_bath` function.
35
35
Various **other methods** for thermostatting exist, such as the Anderson, Nosé-Hoover, or the Berendsen [[3-6](#references)].
36
36
37
37
In order to achieve the NPT ensemble, it is necessary to use a **barostat** in addition to a thermostat.
In the above equation, $r_x$ is the distance between the two particles in the $x$-dimension and $\mathbf{r}$ is the overall distance vector.
86
86
The above equation must be determined to find the force in each dimension.
87
-
However, this lesson will **only** consider particles interacting in a one-dimensional space.
87
+
However, currently we will **only** consider particles interacting in a one-dimensional space.
88
88
89
89
The Python code below shows how to determine the acceleration on each atom of argon due to each other atom of argon.
90
90
Due to Newton's third law, we are able to **increase the efficiency** of this algorithm as the force on atom $i$ will be equal and opposite to the force on atom $j$.
Copy file name to clipboardExpand all lines: _build/parameterisation/intro.md
+3-1Lines changed: 3 additions & 1 deletion
Display the source diff
Display the rich diff
Original file line number
Diff line number
Diff line change
@@ -30,12 +30,14 @@ This means that the parameters should really be obtained by optimising them with
30
30
Commonly this involves either experimental measurements, e.g. X-ray crystallography, or quantum mechanical calculations; we will be focusing on the latter.
31
31
32
32
More can be found out about quantum mechanical calculations in the textbooks mentioned in the introduction (in particular Jeremy Harvey's Computational Chemistry Primer [[1](#references)]).
33
-
However, for our current purposes we only need to remember the previous lesson where it was introduced that quantum calculations are more accurate than classical simulations.
33
+
However, for our current purposes we only need to remember that quantum calculations are more accurate than classical simulations.
34
34
35
35
### Quantum mechanical calculations
36
36
37
37
These are more accurate then classical simulations. However, they are severely limited in the system size, with a maximum simulation size in the order of hundreds atoms.
38
38
39
+
### Parameterising a Lennard-Jones interaction
40
+
39
41
We will stick with the example of a Lennard-Jones interaction, however the arguments and methods discussed are **extensible to all different interaction types**.
40
42
To generate the potential energy model between two particles of argon, we could conduct quantum mechanical calculations at a range of inter-atom separations, from 3 to 8 Å, finding the energy between the two particles at each separation.
Copy file name to clipboardExpand all lines: content/classical_methods/intro.md
+2-2Lines changed: 2 additions & 2 deletions
Display the source diff
Display the rich diff
Original file line number
Diff line number
Diff line change
@@ -1,7 +1,7 @@
1
1
# Classical methods
2
2
3
3
**Classical methods** is a phrase used to describe techniques that make use of a potential model (sometimes called a force-field) to simulated chemical systems.
4
-
These can be molecular dynamics (which we will cover in this module), Monte Carlo, Langevin dynamics, etc.
4
+
These can be molecular dynamics (which we will cover in this resource), Monte Carlo, Langevin dynamics, etc.
5
5
6
6
In order to simulate a **real** chemical system, it is necessary to model the electrons and their interactions.
7
7
This is achieved by using quantum mechanical calculations, where the energy of a chemical system is calculated by finding an approximate solution to the Schrödinger equation.
@@ -12,4 +12,4 @@ With this assumption envoked, it is possible to find the potential energy using
12
12
The first is some ground-state quantum mechanical method (e.g. density functional theory), however as mentioned above these are limited in system size.
13
13
The alternative involves modelling the electron distributions with some mathematical function to determine the potential energy.
14
14
These **potential models** are usually faster to calculate than the quantum mechanical method and therefore may be used on larger systems.
15
-
However, this **simplification** does have a drawback in that the correct potential models must be determined for each system to accurately determine the energy.
15
+
However, this **simplification** does have a drawback in that the correct potential models must be determined for each system to accurately determine the energy.
Copy file name to clipboardExpand all lines: content/important_considerations/ensembles.md
+2-2Lines changed: 2 additions & 2 deletions
Display the source diff
Display the rich diff
Original file line number
Diff line number
Diff line change
@@ -1,6 +1,6 @@
1
1
## Ensembles
2
2
3
-
The molecular dynamics algorithm outlined in the previous lesson makes use of the NVE ensemble (also known as the microcanonical ensemble), where the number of particles (N), volume of the system (V), and energy of the system (E) are all **kept constant**.
3
+
The molecular dynamics algorithm outlined in the previously makes use of the NVE ensemble (also known as the microcanonical ensemble), where the number of particles (N), volume of the system (V), and energy of the system (E) are all **kept constant**.
4
4
This is not the only, or the most accurate, ensemble that exists, there is also other such as:
5
5
- NVT (canonical): number of particles (N), volume of system (V), temperature of the simulation (T)
6
6
- NPT (isothermal-isobaric): number of particles (N), pressure of system (P), temperature of the simulation (T)
@@ -19,7 +19,7 @@ This means that the velocities of the particles may be rescaled by the following
where $T_{\text{target}}$ is the target temperature for the themostat, and $\bar{T}$ is the average simulation temperature.
22
-
pylj [[1,2](#references)], the software that you shall use in the next lesson uses this method for producing an NVT simulation, using the `heat_bath` function.
22
+
pylj [[1,2](#references)], the software that you shall use later uses this method for producing an NVT simulation, using the `heat_bath` function.
23
23
Various **other methods** for thermostatting exist, such as the Anderson, Nosé-Hoover, or the Berendsen [[3-6](#references)].
24
24
25
25
In order to achieve the NPT ensemble, it is necessary to use a **barostat** in addition to a thermostat.
Copy file name to clipboardExpand all lines: content/molecular_dynamics/intro.ipynb
+1-1Lines changed: 1 addition & 1 deletion
Original file line number
Diff line number
Diff line change
@@ -102,7 +102,7 @@
102
102
"\n",
103
103
"In the above equation, $r_x$ is the distance between the two particles in the $x$-dimension and $\\mathbf{r}$ is the overall distance vector. \n",
104
104
"The above equation must be determined to find the force in each dimension. \n",
105
-
"However, this lesson will **only** consider particles interacting in a one-dimensional space. \n",
105
+
"However, currently we will **only** consider particles interacting in a one-dimensional space. \n",
106
106
"\n",
107
107
"The Python code below shows how to determine the acceleration on each atom of argon due to each other atom of argon. \n",
108
108
"Due to Newton's third law, we are able to **increase the efficiency** of this algorithm as the force on atom $i$ will be equal and opposite to the force on atom $j$. "
Copy file name to clipboardExpand all lines: content/parameterisation/intro.ipynb
+3-1Lines changed: 3 additions & 1 deletion
Original file line number
Diff line number
Diff line change
@@ -37,7 +37,7 @@
37
37
"Commonly this involves either experimental measurements, e.g. X-ray crystallography, or quantum mechanical calculations; we will be focusing on the latter. \n",
38
38
"\n",
39
39
"More can be found out about quantum mechanical calculations in the textbooks mentioned in the introduction (in particular Jeremy Harvey's Computational Chemistry Primer [[1](#references)]).\n",
40
-
"However, for our current purposes we only need to remember the previous lesson where it was introduced that quantum calculations are more accurate than classical simulations."
40
+
"However, for our current purposes we only need to remember that quantum calculations are more accurate than classical simulations."
41
41
]
42
42
},
43
43
{
@@ -53,6 +53,8 @@
53
53
"cell_type": "markdown",
54
54
"metadata": {},
55
55
"source": [
56
+
"### Parameterising a Lennard-Jones interaction\n",
57
+
"\n",
56
58
"We will stick with the example of a Lennard-Jones interaction, however the arguments and methods discussed are **extensible to all different interaction types**. \n",
57
59
"To generate the potential energy model between two particles of argon, we could conduct quantum mechanical calculations at a range of inter-atom separations, from 3 to 8 Å, finding the energy between the two particles at each separation.\n",
0 commit comments