https://chemwiki.ch.ic.ac.uk/api.php?action=feedcontributions&feedformat=atom&user=Js6313ChemWiki - User contributions [en]2026-05-03T08:01:40ZUser contributionsMediaWiki 1.43.0https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553573Rep:Mod:js6313 mgo2016-03-11T11:48:43Z<p>Js6313: </p>
<hr />
<div>=== Introduction ===<br />
[[File:Js6313Lattice primative.png|thumb|'''''Fig 1. '''''The primative unit cell of MgO]]<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. <br />
<br />
When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges. The unit cell is a parallelepiped with <br />
<br />
<math>\textstyle \alpha = \beta = \gamma = 60^{\circ} </math> and <math>\textstyle a = b = c = 2.9783 \AA </math> which gives a conventional lattice constant of <math>\textstyle 2.10597 \AA </math>.<br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== Overview of The First Brillouin Zone, Dispersion Curves and the Density of States (DOS) ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|'''''Fig 2. '''''Dispersion curve taken with a grid size of 50]]<br />
In a crystal with lattice vector ''a'', the wavevector <math>\textstyle k = \frac{2 \pi}{\lambda} </math>. Only wavevectors between <math>\textstyle -\frac{\pi}{a} \text{ and } \frac{\pi}{a} </math> are unique; these are in what is known as the first Brillouin zone. Every k-point in this domain is unique, and has it's own energy level which under ideality is given by:<br />
<br />
<math> E(k) = \frac{\hbar^2 k^2}{2m} + E_{pot} </math><br />
<br />
Note that this only holds true within the stated boundary conditions for the first Brillouin zone.<br />
<br />
Dispersion curves are a plot of frequency against k, within the first Brillouin zone. Each line is a phonon - at any point, a maximum of six phonons is observed. This is because there are 2 atoms in each unit cell, each with 3 vibrational degrees of freedom. At some points (where lines converge), phonons are degenerate.<br />
[[File:Js6313dos.png|thumb|'''''Fig 3. '''''Density of phonon states along k-point L (left) and along 50 k-points (right)]]<br />
A density of states (DOS) graph shows the frequency of phonons along a given number of k-points; setting this to 1 gives a view along the L k-point shown below. Note that the convergence of the lower energy curves corresponds to a higher DOS. As the number of sampled k-points increases, a more representative view of the density of states is gained.<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation vs Molecular Dynamics ===<br />
For the quasi-harmonic approximation, Every MgO bond can be thought of as a mass on a spring, the motion of which is governed by Hooke's law:<br />
<br />
<math>F = -k(x - \bar{x})</math><br />
<br />
This is a good approximation at very low temperatures, where bond lengths do not change significantly, however it breaks down once vibrations start to stretch bonds closer to their breaking points. The QMA is a variation on this model whereby k is dependent on the volume, and that the and that the harmonic approximation works for all volumes.<br />
<br />
Molecular dynamics takes a set of starting conditions, and takes a series of snapshots at regular intervals (5 ps in this instance), until the properties reach equilibrium. A random selection of velocities is selected such that the overall temperature of the 32 unit cells is that specified beforehand; from here, the forces are computed, and then:<br />
<br />
<math>a = \frac{F}{m}</math><br />
<br />
<math>v_{\text{new}} = v_{\text{old}} + a \delta t</math><br />
<br />
<math>r_{\text{new}} = r_{\text{old}} + v_{\text{new}} \delta t </math><br />
<br />
where r is the position, v is the velocity and a is the acceleration.<br />
<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly set the temperature of 300 K and the pressure to 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brillouin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation. A comparison of free energies can be seen in '''''Fig 5. '''''Accuracies to within 1 meV, 0.5 meV and 0.1 meV require grid sizes of 2, 2 and 3 respectively.<br />
<br />
=== Thermal Expansion of MgO Using the Quasi-Harmonic Approximation and Molecular Dynamics ===<br />
[[File:Js6313Lattice cube.png|left|thumb|'''''Fig 4. '''''The conventional unit cell for MgO]]<br />
[[File:Temperturevsfreeenergy.png|thumb|'''''Fig 5. '''''Graph showing temperature vs free energy for the QMA (green) and MD (red)]]<br />
The calculations were run at temperatures ranging from 0 K to 6000 K; past 2000 K when using the QMA, DLV gave error messages, possibly due to the dissociation of the lattice. The results are summarised in '''''Fig 5. '''''The two models give similar results at low temperatures, but diverge at T > 450 K. This is because anharmonicity increaces with increacing temperature, so the QM model breaks down (it predicts that free energy will fall with tempreature, whereas the MD model uses the formula <math>\textstyle E_k = \frac{k_B T}{2} </math>. This means that kinetic energy of the individual ions and thus the free energy increases with temperature, as shown in the graph and proven by experiment.<ref><nowiki>http://www.asminternational.org/documents/10192/3449965/ACFAAD6.pdf/2d574bfc-e104-48c5-8d8e-c33ffe91c3a0</nowiki></ref><br />
<br />
==== The Thermal Expansion Coefficient ====<br />
[[File:Volume vs temperature2.png|thumb|'''''Fig 6. '''''Primative cell volume as a function of temperature for the QMA (green) and using MD (red)]]<br />
The thermal expansion coefficient, <math>\textstyle \alpha_V </math>, was calculated using by comparison of the cell volumes over temperatures ranging from 100 K to 1000 K; this can be seen in '''''Fig 6. '''''From this, it can be seen that while MD shows linearity all the way down to 100 K, the QMA deviates at temperatures below 200 K. Because <math>\textstyle \alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math>, the units of volume cancel out and the coefficient is the gradient of the graph - for the QMA this is found to be <math>\textstyle 4.47 \times 10^{-4} K^{-1}</math>, and for MD it is <math>\textstyle 5.81 \times 10^{-4} K^{-1}</math>. This is in good agreement with literature.<ref><nowiki>NIELSEN, T. H., and M. H. LEIPOLD. "Thermal Expansion In Air Of Ceramic Oxides To 2200Oc". Journal of the American Ceramic Society 46.8 (1963): 381-387.</nowiki></ref><br />
<br />
The MD model is, in this case, more accurate. The QMA model uses only a single unit cell, and assumes periodicity of vibrations which means that the amplitude of the phonons must not cause interraction with adjacent cells. The MD model makes no such assumptions, and does not in fact consider the unit cell at all; rather, it computes the motion and position of every ion individually with respect to their immediate surroundings, to give a more physically realistic picture.<br />
<br />
=== Conclusions ===<br />
Both of the models used were accurate for determining the thermal expansion coefficient, whereas the QMA deviated enormously from experiment (and theory) at temperatures above 500 K. This can be explained by the fact that at higher temperatures, anharmonicity increaces greatly which causes innacuracies in the QMA. In order to improve the models, the ions could take into account orbitals rather than being approximated to round charges.<br />
<br />
=== References ===<br />
<references /></div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553571Rep:Mod:js6313 mgo2016-03-11T11:47:34Z<p>Js6313: </p>
<hr />
<div>=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''[[File:Js6313Lattice primative.png|thumb|'''''Fig 1. '''''The primative unit cell of MgO]]<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. <br />
<br />
When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges. The unit cell is a parallelepiped with <br />
<br />
<math>\textstyle \alpha = \beta = \gamma = 60^{\circ} </math> and <math>\textstyle a = b = c = 2.9783 \AA </math> which gives a conventional lattice constant of <math>\textstyle 2.10597 \AA </math>.<br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== Overview of The First Brillouin Zone, Dispersion Curves and the Density of States (DOS) ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|'''''Fig 2. '''''Dispersion curve taken with a grid size of 50]]<br />
In a crystal with lattice vector ''a'', the wavevector <math>\textstyle k = \frac{2 \pi}{\lambda} </math>. Only wavevectors between <math>\textstyle -\frac{\pi}{a} \text{ and } \frac{\pi}{a} </math> are unique; these are in what is known as the first Brillouin zone. Every k-point in this domain is unique, and has it's own energy level which under ideality is given by:<br />
<br />
<math> E(k) = \frac{\hbar^2 k^2}{2m} + E_{pot} </math><br />
<br />
Note that this only holds true within the stated boundary conditions for the first Brillouin zone.<br />
<br />
Dispersion curves are a plot of frequency against k, within the first Brillouin zone. Each line is a phonon - at any point, a maximum of six phonons is observed. This is because there are 2 atoms in each unit cell, each with 3 vibrational degrees of freedom. At some points (where lines converge), phonons are degenerate.<br />
[[File:Js6313dos.png|thumb|'''''Fig 3. '''''Density of phonon states along k-point L (left) and along 50 k-points (right)]]<br />
A density of states (DOS) graph shows the frequency of phonons along a given number of k-points; setting this to 1 gives a view along the L k-point shown below. Note that the convergence of the lower energy curves corresponds to a higher DOS. As the number of sampled k-points increases, a more representative view of the density of states is gained.<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation vs Molecular Dynamics ===<br />
For the quasi-harmonic approximation, Every MgO bond can be thought of as a mass on a spring, the motion of which is governed by Hooke's law:<br />
<br />
<math>F = -k(x - \bar{x})</math><br />
<br />
This is a good approximation at very low temperatures, where bond lengths do not change significantly, however it breaks down once vibrations start to stretch bonds closer to their breaking points. The QMA is a variation on this model whereby k is dependent on the volume, and that the and that the harmonic approximation works for all volumes.<br />
<br />
Molecular dynamics takes a set of starting conditions, and takes a series of snapshots at regular intervals (5 ps in this instance), until the properties reach equilibrium. A random selection of velocities is selected such that the overall temperature of the 32 unit cells is that specified beforehand; from here, the forces are computed, and then:<br />
<br />
<math>a = \frac{F}{m}</math><br />
<br />
<math>v_{\text{new}} = v_{\text{old}} + a \delta t</math><br />
<br />
<math>r_{\text{new}} = r_{\text{old}} + v_{\text{new}} \delta t </math><br />
<br />
where r is the position, v is the velocity and a is the acceleration.<br />
<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly set the temperature of 300 K and the pressure to 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brillouin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation. A comparison of free energies can be seen in '''''Fig 5. '''''Accuracies to within 1 meV, 0.5 meV and 0.1 meV require grid sizes of 2, 2 and 3 respectively.<br />
<br />
=== Thermal Expansion of MgO Using the Quasi-Harmonic Approximation and Molecular Dynamics ===<br />
[[File:Js6313Lattice cube.png|left|thumb|'''''Fig 4. '''''The conventional unit cell for MgO]]<br />
[[File:Temperturevsfreeenergy.png|thumb|'''''Fig 5. '''''Graph showing temperature vs free energy for the QMA (green) and MD (red)]]<br />
The calculations were run at temperatures ranging from 0 K to 6000 K; past 2000 K when using the QMA, DLV gave error messages, possibly due to the dissociation of the lattice. The results are summarised in '''''Fig 5. '''''The two models give similar results at low temperatures, but diverge at T > 450 K. This is because anharmonicity increaces with increacing temperature, so the QM model breaks down (it predicts that free energy will fall with tempreature, whereas the MD model uses the formula <math>\textstyle E_k = \frac{k_B T}{2} </math>. This means that kinetic energy of the individual ions and thus the free energy increases with temperature, as shown in the graph and proven by experiment.<ref><nowiki>http://www.asminternational.org/documents/10192/3449965/ACFAAD6.pdf/2d574bfc-e104-48c5-8d8e-c33ffe91c3a0</nowiki></ref><br />
<br />
==== The Thermal Expansion Coefficient ====<br />
[[File:Volume vs temperature2.png|thumb|'''''Fig 6. '''''Primative cell volume as a function of temperature for the QMA (green) and using MD (red)]]<br />
The thermal expansion coefficient, <math>\textstyle \alpha_V </math>, was calculated using by comparison of the cell volumes over temperatures ranging from 100 K to 1000 K; this can be seen in '''''Fig 6. '''''From this, it can be seen that while MD shows linearity all the way down to 100 K, the QMA deviates at temperatures below 200 K. Because <math>\textstyle \alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math>, the units of volume cancel out and the coefficient is the gradient of the graph - for the QMA this is found to be <math>\textstyle 4.47 \times 10^{-4} K^{-1}</math>, and for MD it is <math>\textstyle 5.81 \times 10^{-4} K^{-1}</math>. This is in good agreement with literature.<ref><nowiki>NIELSEN, T. H., and M. H. LEIPOLD. "Thermal Expansion In Air Of Ceramic Oxides To 2200Oc". Journal of the American Ceramic Society 46.8 (1963): 381-387.</nowiki></ref><br />
<br />
The MD model is, in this case, more accurate. The QMA model uses only a single unit cell, and assumes periodicity of vibrations which means that the amplitude of the phonons must not cause interraction with adjacent cells. The MD model makes no such assumptions, and does not in fact consider the unit cell at all; rather, it computes the motion and position of every ion individually with respect to their immediate surroundings, to give a more physically realistic picture.<br />
<br />
=== Conclusions ===<br />
Both of the models used were accurate for determining the thermal expansion coefficient, whereas the QMA deviated enormously from experiment (and theory) at temperatures above 500 K. This can be explained by the fact that at higher temperatures, anharmonicity increaces greatly which causes innacuracies in the QMA. In order to improve the models, the ions could take into account orbitals rather than being approximated to round charges.<br />
<br />
=== References ===<br />
<references /></div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553570Rep:Mod:js6313 mgo2016-03-11T11:46:51Z<p>Js6313: </p>
<hr />
<div>=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''[[File:Js6313Lattice primative.png|thumb|'''''Fig 1. '''''The primative unit cell of MgO]]<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. <br />
<br />
When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges. The unit cell is a parallelepiped with <br />
<br />
<math>\textstyle \alpha = \beta = \gamma = 60^{\circ} </math> and <math>\textstyle a = b = c = 2.9783 \AA </math> which gives a conventional lattice constant of <math>\textstyle 2.10597 \AA </math>.<br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== Overview of The First Brillouin Zone, Dispersion Curves and the Density of States (DOS) ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|'''''Fig 2. '''''Dispersion curve taken with a grid size of 50]]<br />
In a crystal with lattice vector ''a'', the wavevector <math>\textstyle k = \frac{2 \pi}{\lambda} </math>. Only wavevectors between <math>\textstyle -\frac{\pi}{a} \text{ and } \frac{\pi}{a} </math> are unique; these are in what is known as the first Brillouin zone. Every k-point in this domain is unique, and has it's own energy level which under ideality is given by:<br />
<br />
<math> E(k) = \frac{\hbar^2 k^2}{2m} + E_{pot} </math><br />
<br />
Note that this only holds true within the stated boundary conditions for the first Brillouin zone.<br />
<br />
Dispersion curves are a plot of frequency against k, within the first Brillouin zone. Each line is a phonon - at any point, a maximum of six phonons is observed. This is because there are 2 atoms in each unit cell, each with 3 vibrational degrees of freedom. At some points (where lines converge), phonons are degenerate.<br />
[[File:Js6313dos.png|thumb|'''''Fig 3. '''''Density of phonon states along k-point L (left) and along 50 k-points (right)]]<br />
A density of states (DOS) graph shows the frequency of phonons along a given number of k-points; setting this to 1 gives a view along the L k-point shown below. Note that the convergence of the lower energy curves corresponds to a higher DOS. As the number of sampled k-points increases, a more representative view of the density of states is gained.<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation vs Molecular Dynamics ===<br />
For the quasi-harmonic approximation, Every MgO bond can be thought of as a mass on a spring, the motion of which is governed by Hooke's law:<br />
<br />
<math>F = -k(x - \bar{x})</math><br />
<br />
This is a good approximation at very low temperatures, where bond lengths do not change significantly, however it breaks down once vibrations start to stretch bonds closer to their breaking points. The QMA is a variation on this model whereby k is dependent on the volume, and that the and that the harmonic approximation works for all volumes.<br />
<br />
Molecular dynamics takes a set of starting conditions, and takes a series of snapshots at regular intervals (5 ps in this instance), until the properties reach equilibrium. A random selection of velocities is selected such that the overall temperature of the 32 unit cells is that specified beforehand; from here, the forces are computed, and then:<br />
<br />
<math>a = \frac{F}{m}</math><br />
<br />
<math>v_{\text{new}} = v_{\text{old}} + a \delta t</math><br />
<br />
<math>r_{\text{new}} = r_{\text{old}} + v_{\text{new}} \delta t </math><br />
<br />
where r is the position, v is the velocity and a is the acceleration.<br />
<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly set the temperature of 300 K and the pressure to 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brillouin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation. A comparison of free energies can be seen in '''''Fig 5. '''''Accuracies to within 1 meV, 0.5 meV and 0.1 meV require grid sizes of 2, 2 and 3 respectively.<br />
<br />
=== Thermal Expansion of MgO Using the Quasi-Harmonic Approximation and Molecular Dynamics ===<br />
[[File:Js6313Lattice cube.png|left|thumb|'''''Fig 4. '''''The conventional unit cell for MgO]]<br />
[[File:Temperturevsfreeenergy.png|thumb|'''''Fig 5. '''''Graph showing temperature vs free energy for the QMA (green) and MD (red)]]<br />
The calculations were run at temperatures ranging from 0 K to 6000 K; past 2000 K when using the QMA, DLV gave error messages, possibly due to the dissociation of the lattice. The results are summarised in '''''Fig 5. '''''The two models give similar results at low temperatures, but diverge at T > 450 K. This is because anharmonicity increaces with increacing temperature, so the QM model breaks down (it predicts that free energy will fall with tempreature, whereas the MD model uses the formula <math>\textstyle E_k = \frac{k_B T}{2} </math>. This means that kinetic energy of the individual ions and thus the free energy increases with temperature, as shown in the graph and proven by experiment.<ref><nowiki>http://www.asminternational.org/documents/10192/3449965/ACFAAD6.pdf/2d574bfc-e104-48c5-8d8e-c33ffe91c3a0</nowiki></ref><br />
<br />
==== The Thermal Expansion Coefficient ====<br />
[[File:Volume vs temperature2.png|thumb|'''''Fig 6. '''''Primative cell volume as a function of temperature for the QMA (green) and using MD (red)]]<br />
The thermal expansion coefficient, <math>\textstyle \alpha_V </math>, was calculated using by comparison of the cell volumes over temperatures ranging from 100 K to 1000 K; this can be seen in '''''Fig 6. '''''From this, it can be seen that while MD shows linearity all the way down to 100 K, the QMA deviates at temperatures below 200 K. Because <math>\textstyle \alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math>, the units of volume cancel out and the coefficient is the gradient of the graph - for the QMA this is found to be <math>\textstyle 4.47 \times 10^{-4} K^{-1}</math>, and for MD it is <math>\textstyle 5.81 \times 10^{-4} K^{-1}</math>. This is in good agreement with literature.<ref><nowiki>NIELSEN, T. H., and M. H. LEIPOLD. "Thermal Expansion In Air Of Ceramic Oxides To 2200Oc". Journal of the American Ceramic Society 46.8 (1963): 381-387.</nowiki></ref><br />
<br />
The MD model is, in this case, more accurate. The QMA model uses only a single unit cell, and assumes periodicity of vibrations which means that the amplitude of the phonons must not cause interraction with adjacent cells. The MD model makes no such assumptions, and does not in fact consider the unit cell at all; rather, it computes the motion and position of every ion individually with respect to their immediate surroundings, to give a more physically realistic picture.<br />
<br />
=== Conclusions ===<br />
Both of the models used were accurate for determining the thermal expansion coefficient, whereas the QMA deviated enormously from experiment (and theory) at temperatures above 500 K. This can be explained by the fact that at higher temperatures, anharmonicity increaces greatly which causes innacuracies in the QMA. In order to improve the models, the ions could take into account orbitals rather than being approximated to round charges.</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553567Rep:Mod:js6313 mgo2016-03-11T11:45:01Z<p>Js6313: </p>
<hr />
<div>=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''[[File:Js6313Lattice primative.png|thumb|'''''Fig 1. '''''The primative unit cell of MgO]]<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. <br />
<br />
When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges. The unit cell is a parallelepiped with <br />
<br />
<math>\textstyle \alpha = \beta = \gamma = 60^{\circ} </math> and <math>\textstyle a = b = c = 2.9783 \AA </math> which gives a conventional lattice constant of <math>\textstyle 2.10597 \AA </math>.<br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== Overview of The First Brillouin Zone, Dispersion Curves and the Density of States (DOS) ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|'''''Fig 2. '''''Dispersion curve taken with a grid size of 50]]<br />
In a crystal with lattice vector ''a'', the wavevector <math>\textstyle k = \frac{2 \pi}{\lambda} </math>. Only wavevectors between <math>\textstyle -\frac{\pi}{a} \text{ and } \frac{\pi}{a} </math> are unique; these are in what is known as the first Brillouin zone. Every k-point in this domain is unique, and has it's own energy level which under ideality is given by:<br />
<br />
<math> E(k) = \frac{\hbar^2 k^2}{2m} + E_{pot} </math><br />
<br />
Note that this only holds true within the stated boundary conditions for the first Brillouin zone.<br />
<br />
Dispersion curves are a plot of frequency against k, within the first Brillouin zone. Each line is a phonon - at any point, a maximum of six phonons is observed. This is because there are 2 atoms in each unit cell, each with 3 vibrational degrees of freedom. At some points (where lines converge), phonons are degenerate.<br />
[[File:Js6313dos.png|thumb|'''''Fig 3. '''''Density of phonon states along k-point L (left) and along 50 k-points (right)]]<br />
A density of states (DOS) graph shows the frequency of phonons along a given number of k-points; setting this to 1 gives a view along the L k-point shown below. Note that the convergence of the lower energy curves corresponds to a higher DOS. As the number of sampled k-points increases, a more representative view of the density of states is gained.<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation vs Molecular Dynamics ===<br />
For the quasi-harmonic approximation, Every MgO bond can be thought of as a mass on a spring, the motion of which is governed by Hooke's law:<br />
<br />
<math>F = -k(x - \bar{x})</math><br />
<br />
This is a good approximation at very low temperatures, where bond lengths do not change significantly, however it breaks down once vibrations start to stretch bonds closer to their breaking points. The QMA is a variation on this model whereby k is dependent on the volume, and that the and that the harmonic approximation works for all volumes.<br />
<br />
Molecular dynamics takes a set of starting conditions, and takes a series of snapshots at regular intervals (5 ps in this instance), until the properties reach equilibrium. A random selection of velocities is selected such that the overall temperature of the 32 unit cells is that specified beforehand; from here, the forces are computed, and then:<br />
<br />
<math>a = \frac{F}{m}</math><br />
<br />
<math>v_{\text{new}} = v_{\text{old}} + a \delta t</math><br />
<br />
<math>r_{\text{new}} = r_{\text{old}} + v_{\text{new}} \delta t </math><br />
<br />
where r is the position, v is the velocity and a is the acceleration.<br />
<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly set the temperature of 300 K and the pressure to 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brillouin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation. A comparison of free energies can be seen in '''''Fig 5. '''''Accuracies to within 1 meV, 0.5 meV and 0.1 meV require grid sizes of 2, 2 and 3 respectively.<br />
<br />
=== Thermal Expansion of MgO Using the Quasi-Harmonic Approximation and Molecular Dynamics ===<br />
[[File:Js6313Lattice cube.png|left|thumb|'''''Fig 4. '''''The conventional unit cell for MgO]]<br />
[[File:Temperturevsfreeenergy.png|thumb|'''''Fig 5. '''''Graph showing temperature vs free energy for the QMA (green) and MD (red)]]<br />
The calculations were run at temperatures ranging from 0 K to 6000 K; past 2000 K when using the QMA, DLV gave error messages, possibly due to the dissociation of the lattice. The results are summarised in '''''Fig 5. '''''The two models give similar results at low temperatures, but diverge at T > 450 K. This is because anharmonicity increaces with increacing temperature, so the QM model breaks down (it predicts that free energy will fall with tempreature, whereas the MD model uses the formula <math>\textstyle E_k = \frac{k_B T}{2} </math>. This means that kinetic energy of the individual ions and thus the free energy increases with temperature, as shown in the graph and proven by experiment.<ref><nowiki>http://www.asminternational.org/documents/10192/3449965/ACFAAD6.pdf/2d574bfc-e104-48c5-8d8e-c33ffe91c3a0</nowiki></ref><br />
<br />
==== The Thermal Expansion Coefficient ====<br />
[[File:Volume vs temperature2.png|thumb|'''''Fig 6. '''''Primative cell volume as a function of temperature for the QMA (green) and using MD (red)]]<br />
The thermal expansion coefficient, <math>\textstyle \alpha_V </math>, was calculated using by comparison of the cell volumes over temperatures ranging from 100 K to 1000 K; this can be seen in '''''Fig 6. '''''From this, it can be seen that while MD shows linearity all the way down to 100 K, the QMA deviates at temperatures below 200 K. Because <math>\textstyle \alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math>, the units of volume cancel out and the coefficient is the gradient of the graph - for the QMA this is found to be <math>\textstyle 4.47 \times 10^{-4} K^{-1}</math>, and for MD it is <math>\textstyle 5.81 \times 10^{-4} K^{-1}</math>. This is in good agreement with literature.<ref><nowiki>NIELSEN, T. H., and M. H. LEIPOLD. "Thermal Expansion In Air Of Ceramic Oxides To 2200Oc". Journal of the American Ceramic Society 46.8 (1963): 381-387.</nowiki></ref><br />
<br />
The MD model is, in this case, more accurate. The QMA model uses only a single unit cell, and assumes periodicity of vibrations which means that the amplitude of the phonons must not cause interraction with adjacent cells. The MD model makes no such assumptions, and does not in fact consider the unit cell at all; rather, it computes the motion and position of every ion individually with respect to their immediate surroundings, to give a more physically realistic picture.<br />
<br />
=== Conclusion ===<br />
Both of the models used were accurate for determining the thermal expansion coefficient, whereas the QMA deviated enormously from experiment (and theory) at temperatures above 500 K. This can be explained by the fact that at higher temperatures, anharmonicity increaces greatly which causes innacuracies in the QMA. In order to improve the models, the ions could take into account orbitals rather than being approximated to round charges.</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553559Rep:Mod:js6313 mgo2016-03-11T11:40:15Z<p>Js6313: </p>
<hr />
<div>=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''[[File:Js6313Lattice primative.png|thumb|'''''Fig 1. '''''The primative unit cell of MgO]]<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. <br />
<br />
When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges. The unit cell is a parallelepiped with <br />
<br />
<math>\textstyle \alpha = \beta = \gamma = 60^{\circ} </math> and <math>\textstyle a = b = c = 2.9783 \AA </math> which gives a conventional lattice constant of <math>\textstyle 2.10597 \AA </math>.<br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== Overview of The First Brillouin Zone, Dispersion Curves and the Density of States (DOS) ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|'''''Fig 2. '''''Dispersion curve taken with a grid size of 50]]<br />
In a crystal with lattice vector ''a'', the wavevector <math>\textstyle k = \frac{2 \pi}{\lambda} </math>. Only wavevectors between <math>\textstyle -\frac{\pi}{a} \text{ and } \frac{\pi}{a} </math> are unique; these are in what is known as the first Brillouin zone. Every k-point in this domain is unique, and has it's own energy level which under ideality is given by:<br />
<br />
<math> E(k) = \frac{\hbar^2 k^2}{2m} + E_{pot} </math><br />
<br />
Note that this only holds true within the stated boundary conditions for the first Brillouin zone.<br />
<br />
Dispersion curves are a plot of frequency against k, within the first Brillouin zone. Each line is a phonon - at any point, a maximum of six phonons is observed. This is because there are 2 atoms in each unit cell, each with 3 vibrational degrees of freedom. At some points (where lines converge), phonons are degenerate.<br />
[[File:Js6313dos.png|thumb|'''''Fig 3. '''''Density of phonon states along k-point L (left) and along 50 k-points (right)]]<br />
A density of states (DOS) graph shows the frequency of phonons along a given number of k-points; setting this to 1 gives a view along the L k-point shown below. Note that the convergence of the lower energy curves corresponds to a higher DOS. As the number of sampled k-points increases, a more representative view of the density of states is gained.<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation vs Molecular Dynamics ===<br />
For the quasi-harmonic approximation, Every MgO bond can be thought of as a mass on a spring, the motion of which is governed by Hooke's law:<br />
<br />
<math>F = -k(x - \bar{x})</math><br />
<br />
This is a good approximation at very low temperatures, where bond lengths do not change significantly, however it breaks down once vibrations start to stretch bonds closer to their breaking points. The QMA is a variation on this model whereby k is dependent on the volume, and that the and that the harmonic approximation works for all volumes.<br />
<br />
Molecular dynamics takes a set of starting conditions, and takes a series of snapshots at regular intervals (5 ps in this instance), until the properties reach equilibrium. A random selection of velocities is selected such that the overall temperature of the 32 unit cells is that specified beforehand; from here, the forces are computed, and then:<br />
<br />
<math>a = \frac{F}{m}</math><br />
<br />
<math>v_{\text{new}} = v_{\text{old}} + a \delta t</math><br />
<br />
<math>r_{\text{new}} = r_{\text{old}} + v_{\text{new}} \delta t </math><br />
<br />
where r is the position, v is the velocity and a is the acceleration.<br />
<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly set the temperature of 300 K and the pressure to 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brillouin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation. A comparison of free energies can be seen in '''''Fig 5. '''''Accuracies to within 1 meV, 0.5 meV and 0.1 meV require grid sizes of 2, 2 and 3 respectively.<br />
<br />
=== Thermal Expansion of MgO Using the Quasi-Harmonic Approximation and Molecular Dynamics ===<br />
[[File:Js6313Lattice cube.png|left|thumb|'''''Fig 4. '''''The conventional unit cell for MgO]]<br />
[[File:Temperturevsfreeenergy.png|thumb|'''''Fig 5. '''''Graph showing temperature vs free energy for the QMA (green) and MD (red)]]<br />
The calculations were run at temperatures ranging from 0 K to 6000 K; past 2000 K when using the QMA, DLV gave error messages, possibly due to the dissociation of the lattice. The results are summarised in '''''Fig 5. '''''The two models give similar results at low temperatures, but diverge at T > 450 K. This is because anharmonicity increaces with increacing temperature, so the QM model breaks down (it predicts that free energy will fall with tempreature, whereas the MD model uses the formula <math>\textstyle E_k = \frac{k_B T}{2} </math>. This means that kinetic energy of the individual ions and thus the free energy increases with temperature, as shown in the graph and proven by experiment.<ref><nowiki>http://www.asminternational.org/documents/10192/3449965/ACFAAD6.pdf/2d574bfc-e104-48c5-8d8e-c33ffe91c3a0</nowiki></ref><br />
<br />
==== The Thermal Expansion Coefficient ====<br />
[[File:Volume vs temperature2.png|thumb|'''''Fig 6. '''''Primative cell volume as a function of temperature for the QMA (green) and using MD (red)]]<br />
The thermal expansion coefficient, <math>\textstyle \alpha_V </math>, was calculated using by comparison of the cell volumes over temperatures ranging from 100 K to 1000 K; this can be seen in '''''Fig 6. '''''From this, it can be seen that while MD shows linearity all the way down to 100 K, the QMA deviates at temperatures below 200 K. Because <math>\textstyle \alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math>, the units of volume cancel out and the coefficient is the gradient of the graph - for the QMA this is found to be <math>\textstyle 4.47 \times 10^{-4} K^{-1}</math>, and for MD it is <math>\textstyle 5.81 \times 10^{-4} K^{-1}</math>. This is in good agreement with literature.<ref><nowiki>NIELSEN, T. H., and M. H. LEIPOLD. "Thermal Expansion In Air Of Ceramic Oxides To 2200Oc". Journal of the American Ceramic Society 46.8 (1963): 381-387.</nowiki></ref><br />
<br />
The MD model is, in this case, more accurate. The QMA model uses only a single unit cell, and assumes periodicity of vibrations which means that the amplitude of the phonons must not cause interraction with adjacent cells. The MD model makes no such assumptions, and does not in fact consider the unit cell at all; rather, it computes the motion and position of every ion individually with respect to their immediate surroundings, to give a more physically realistic picture.<br />
<br />
=== Conclusion ===<br />
Both of the models used were accurate for determining the thermal expansion coefficient, whereas the QMA deviated enormously from experiment (and theory) at temperatures above 500 K. This can be explained by the fact that at higher temperatures, anharmonicity increaces greatly which causes innacuracies in the QMA, whereas</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553533Rep:Mod:js6313 mgo2016-03-11T11:24:38Z<p>Js6313: </p>
<hr />
<div>=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''[[File:Js6313Lattice primative.png|thumb|'''''Fig 1. '''''The primative unit cell of MgO]]<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. <br />
<br />
When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges. The unit cell is a parallelepiped with <br />
<br />
<math>\textstyle \alpha = \beta = \gamma = 60^{\circ} </math> and <math>\textstyle a = b = c = 2.9783 \AA </math> which gives a conventional lattice constant of <math>\textstyle 2.10597 \AA </math>.<br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== Overview of The First Brillouin Zone, Dispersion Curves and the Density of States (DOS) ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|'''''Fig 2. '''''Dispersion curve taken with a grid size of 50]]<br />
In a crystal with lattice vector ''a'', the wavevector <math>\textstyle k = \frac{2 \pi}{\lambda} </math>. Only wavevectors between <math>\textstyle -\frac{\pi}{a} \text{ and } \frac{\pi}{a} </math> are unique; these are in what is known as the first Brillouin zone. Every k-point in this domain is unique, and has it's own energy level which under ideality is given by:<br />
<br />
<math> E(k) = \frac{\hbar^2 k^2}{2m} + E_{pot} </math><br />
<br />
Note that this only holds true within the stated boundary conditions for the first Brillouin zone.<br />
<br />
Dispersion curves are a plot of frequency against k, within the first Brillouin zone. Each line is a phonon - at any point, a maximum of six phonons is observed. This is because there are 2 atoms in each unit cell, each with 3 vibrational degrees of freedom. At some points (where lines converge), phonons are degenerate.<br />
[[File:Js6313dos.png|thumb|'''''Fig 3. '''''Density of phonon states along k-point L (left) and along 50 k-points (right)]]<br />
A density of states (DOS) graph shows the frequency of phonons along a given number of k-points; setting this to 1 gives a view along the L k-point shown below. Note that the convergence of the lower energy curves corresponds to a higher DOS. As the number of sampled k-points increases, a more representative view of the density of states is gained.<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation vs Molecular Dynamics ===<br />
For the quasi-harmonic approximation, Every MgO bond can be thought of as a mass on a spring, the motion of which is governed by Hooke's law:<br />
<br />
<math>F = -k(x - \bar{x})</math><br />
<br />
This is a good approximation at very low temperatures, where bond lengths do not change significantly, however it breaks down once vibrations start to stretch bonds closer to their breaking points. The QMA is a variation on this model whereby k is dependent on the volume, and that the and that the harmonic approximation works for all volumes.<br />
<br />
Molecular dynamics takes a set of starting conditions, and takes a series of snapshots at regular intervals (5 ps in this instance), until the properties reach equilibrium. A random selection of velocities is selected such that the overall temperature of the 32 unit cells is that specified beforehand; from here, the forces are computed, and then:<br />
<br />
<math>a = \frac{F}{m}</math><br />
<br />
<math>v_{\text{new}} = v_{\text{old}} + a \delta t</math><br />
<br />
<math>r_{\text{new}} = r_{\text{old}} + v_{\text{new}} \delta t </math><br />
<br />
where r is the position, v is the velocity and a is the acceleration.<br />
<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly set the temperature of 300 K and the pressure to 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brillouin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation. A comparison of free energies can be seen in '''''Fig 5. '''''Accuracies to within 1 meV, 0.5 meV and 0.1 meV require grid sizes of 2, 2 and 3 respectively.<br />
<br />
=== Thermal Expansion of MgO Using the Quasi-Harmonic Approximation and Molecular Dynamics ===<br />
[[File:Js6313Lattice cube.png|left|thumb|'''''Fig 4. '''''The conventional unit cell for MgO]]<br />
[[File:Temperturevsfreeenergy.png|thumb|'''''Fig 5. '''''Graph showing temperature vs free energy for the QMA (green) and MD (red)]]<br />
The calculations were run at temperatures ranging from 0 K to 6000 K; past 2000 K when using the QMA, DLV gave error messages, possibly due to the dissociation of the lattice. The results are summarised in '''''Fig 5. '''''The two models give similar results at low temperatures, but diverge at T > 450 K. This is because anharmonicity increaces with increacing temperature, so the QM model breaks down (it predicts that free energy will fall with tempreature, whereas the MD model uses the formula <math>\textstyle E_k = \frac{k_B T}{2} </math>. This means that kinetic energy of the individual ions and thus the free energy increases with temperature, as shown in the graph and proven by experiment.<ref><nowiki>http://www.asminternational.org/documents/10192/3449965/ACFAAD6.pdf/2d574bfc-e104-48c5-8d8e-c33ffe91c3a0</nowiki></ref><br />
<br />
==== The Thermal Expansion Coefficient ====<br />
[[File:Volume vs temperature2.png|thumb|'''''Fig 6. '''''Primative cell volume as a function of temperature for the QMA (green) and using MD (red)]]<br />
The thermal expansion coefficient, <math>\textstyle \alpha_V </math>, was calculated using by comparison of the cell volumes over temperatures ranging from 100 K to 1000 K; this can be seen in '''''Fig 6. '''''From this, it can be seen that while MD shows linearity all the way down to 100 K, the QMA deviates at temperatures below 200 K. Because <math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math>, the units of volume cancel out and the coefficient is the gradient of the graph - for the QMA this is found to be</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553518Rep:Mod:js6313 mgo2016-03-11T11:19:35Z<p>Js6313: </p>
<hr />
<div>=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''[[File:Js6313Lattice primative.png|thumb|'''''Fig 1. '''''The primative unit cell of MgO]]<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. <br />
<br />
When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges. The unit cell is a parallelepiped with <br />
<br />
<math>\textstyle \alpha = \beta = \gamma = 60^{\circ} </math> and <math>\textstyle a = b = c = 2.9783 \AA </math> which gives a conventional lattice constant of <math>\textstyle 2.10597 \AA </math>.<br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== Overview of The First Brillouin Zone, Dispersion Curves and the Density of States (DOS) ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|'''''Fig 2. '''''Dispersion curve taken with a grid size of 50]]<br />
In a crystal with lattice vector ''a'', the wavevector <math>\textstyle k = \frac{2 \pi}{\lambda} </math>. Only wavevectors between <math>\textstyle -\frac{\pi}{a} \text{ and } \frac{\pi}{a} </math> are unique; these are in what is known as the first Brillouin zone. Every k-point in this domain is unique, and has it's own energy level which under ideality is given by:<br />
<br />
<math> E(k) = \frac{\hbar^2 k^2}{2m} + E_{pot} </math><br />
<br />
Note that this only holds true within the stated boundary conditions for the first Brillouin zone.<br />
<br />
Dispersion curves are a plot of frequency against k, within the first Brillouin zone. Each line is a phonon - at any point, a maximum of six phonons is observed. This is because there are 2 atoms in each unit cell, each with 3 vibrational degrees of freedom. At some points (where lines converge), phonons are degenerate.<br />
[[File:Js6313dos.png|thumb|'''''Fig 3. '''''Density of phonon states along k-point L (left) and along 50 k-points (right)]]<br />
A density of states (DOS) graph shows the frequency of phonons along a given number of k-points; setting this to 1 gives a view along the L k-point shown below. Note that the convergence of the lower energy curves corresponds to a higher DOS. As the number of sampled k-points increases, a more representative view of the density of states is gained.<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation vs Molecular Dynamics ===<br />
For the quasi-harmonic approximation, Every MgO bond can be thought of as a mass on a spring, the motion of which is governed by Hooke's law:<br />
<br />
<math>F = -k(x - \bar{x})</math><br />
<br />
This is a good approximation at very low temperatures, where bond lengths do not change significantly, however it breaks down once vibrations start to stretch bonds closer to their breaking points. The QMA is a variation on this model whereby k is dependent on the volume, and that the and that the harmonic approximation works for all volumes.<br />
<br />
Molecular dynamics takes a set of starting conditions, and takes a series of snapshots at regular intervals (5 ps in this instance), until the properties reach equilibrium. A random selection of velocities is selected such that the overall temperature of the 32 unit cells is that specified beforehand; from here, the forces are computed, and then:<br />
<br />
<math>a = \frac{F}{m}</math><br />
<br />
<math>v_{\text{new}} = v_{\text{old}} + a \delta t</math><br />
<br />
<math>r_{\text{new}} = r_{\text{old}} + v_{\text{new}} \delta t </math><br />
<br />
where r is the position, v is the velocity and a is the acceleration.<br />
<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly set the temperature of 300 K and the pressure to 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brillouin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation. A comparison of free energies can be seen in '''''Fig 5. '''''Accuracies to within 1 meV, 0.5 meV and 0.1 meV require grid sizes of 2, 2 and 3 respectively.<br />
<br />
=== Thermal Expansion of MgO Using the Quasi-Harmonic Approximation and Molecular Dynamics ===<br />
[[File:Js6313Lattice cube.png|left|thumb|'''''Fig 4. '''''The conventional unit cell for MgO]]<br />
[[File:Temperturevsfreeenergy.png|thumb|'''''Fig 5. '''''Graph showing temperature vs free energy for the QMA (green) and MD (red)]]<br />
The calculations were run at temperatures ranging from 0 K to 6000 K; past 2000 K when using the QMA, DLV gave error messages, possibly due to the dissociation of the lattice. The results are summarised in '''''Fig 5. '''''The two models give similar results at low temperatures, but diverge at T > 450 K. This is because anharmonicity increaces with increacing temperature, so the QM model breaks down (it predicts that free energy will fall with tempreature, whereas the MD model uses the formula <math>\textstyle E_k = \frac{k_B T}{2} </math>. This means that kinetic energy of the individual ions and thus the free energy increases with temperature, as shown in the graph and proven by experiment.<ref><nowiki>http://www.asminternational.org/documents/10192/3449965/ACFAAD6.pdf/2d574bfc-e104-48c5-8d8e-c33ffe91c3a0</nowiki></ref><br />
<br />
==== The Thermal Expansion Coefficient ====<br />
[[File:Volume vs temperature2.png|thumb|'''''Fig 6. '''''Primative cell volume as a function of temperature for the QMA (green) and using MD (red)]]<br />
The thermal expansion coefficient, <math>\textstyle \alpha_V </math>, was calculated using by comparison of the cell volumes over temperatures ranging from 100 K to 1000 K; this can be seen in '''''Fig 6. '''''From this, it can be seen that while MD shows linearity all the way down to 100 K, the QMA deviates at temperatures below</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553516Rep:Mod:js6313 mgo2016-03-11T11:18:12Z<p>Js6313: </p>
<hr />
<div>=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''[[File:Js6313Lattice primative.png|thumb|'''''Fig 1. '''''The primative unit cell of MgO]]<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. <br />
<br />
When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges. The unit cell is a parallelepiped with <br />
<br />
<math>\textstyle \alpha = \beta = \gamma = 60^{\circ} </math> and <math>\textstyle a = b = c = 2.9783 \AA </math> which gives a conventional lattice constant of <math>\textstyle 2.10597 \AA </math>.<br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== Overview of The First Brillouin Zone, Dispersion Curves and the Density of States (DOS) ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|'''''Fig 2. '''''Dispersion curve taken with a grid size of 50]]<br />
In a crystal with lattice vector ''a'', the wavevector <math>\textstyle k = \frac{2 \pi}{\lambda} </math>. Only wavevectors between <math>\textstyle -\frac{\pi}{a} \text{ and } \frac{\pi}{a} </math> are unique; these are in what is known as the first Brillouin zone. Every k-point in this domain is unique, and has it's own energy level which under ideality is given by:<br />
<br />
<math> E(k) = \frac{\hbar^2 k^2}{2m} + E_{pot} </math><br />
<br />
Note that this only holds true within the stated boundary conditions for the first Brillouin zone.<br />
<br />
Dispersion curves are a plot of frequency against k, within the first Brillouin zone. Each line is a phonon - at any point, a maximum of six phonons is observed. This is because there are 2 atoms in each unit cell, each with 3 vibrational degrees of freedom. At some points (where lines converge), phonons are degenerate.<br />
[[File:Js6313dos.png|thumb|'''''Fig 3. '''''Density of phonon states along k-point L (left) and along 50 k-points (right)]]<br />
A density of states (DOS) graph shows the frequency of phonons along a given number of k-points; setting this to 1 gives a view along the L k-point shown below. Note that the convergence of the lower energy curves corresponds to a higher DOS. As the number of sampled k-points increases, a more representative view of the density of states is gained.<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation vs Molecular Dynamics ===<br />
For the quasi-harmonic approximation, Every MgO bond can be thought of as a mass on a spring, the motion of which is governed by Hooke's law:<br />
<br />
<math>F = -k(x - \bar{x})</math><br />
<br />
This is a good approximation at very low temperatures, where bond lengths do not change significantly, however it breaks down once vibrations start to stretch bonds closer to their breaking points. The QMA is a variation on this model whereby k is dependent on the volume, and that the and that the harmonic approximation works for all volumes.<br />
<br />
Molecular dynamics takes a set of starting conditions, and takes a series of snapshots at regular intervals (5 ps in this instance), until the properties reach equilibrium. A random selection of velocities is selected such that the overall temperature of the 32 unit cells is that specified beforehand; from here, the forces are computed, and then:<br />
<br />
<math>a = F/m</math><br />
<br />
<math>v_{\text{new}} = v_{\text{old}} + a \delta t</math><br />
<br />
<math>r_{\text{new}} = r_{\text{old}} + v_{\text{new}} \delta t </math><br />
<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly set the temperature of 300 K and the pressure to 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brillouin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation. A comparison of free energies can be seen in '''''Fig 5. '''''Accuracies to within 1 meV, 0.5 meV and 0.1 meV require grid sizes of 2, 2 and 3 respectively.<br />
<br />
=== Thermal Expansion of MgO Using the Quasi-Harmonic Approximation and Molecular Dynamics ===<br />
[[File:Js6313Lattice cube.png|left|thumb|'''''Fig 4. '''''The conventional unit cell for MgO]]<br />
[[File:Temperturevsfreeenergy.png|thumb|'''''Fig 5. '''''Graph showing temperature vs free energy for the QMA (green) and MD (red)]]<br />
The calculations were run at temperatures ranging from 0 K to 6000 K; past 2000 K when using the QMA, DLV gave error messages, possibly due to the dissociation of the lattice. The results are summarised in '''''Fig 5. '''''The two models give similar results at low temperatures, but diverge at T > 450 K. This is because anharmonicity increaces with increacing temperature, so the QM model breaks down (it predicts that free energy will fall with tempreature, whereas the MD model uses the formula <math>\textstyle E_k = \frac{k_B T}{2} </math>. This means that kinetic energy of the individual ions and thus the free energy increases with temperature, as shown in the graph and proven by experiment.<ref><nowiki>http://www.asminternational.org/documents/10192/3449965/ACFAAD6.pdf/2d574bfc-e104-48c5-8d8e-c33ffe91c3a0</nowiki></ref><br />
<br />
==== The Thermal Expansion Coefficient ====<br />
[[File:Volume vs temperature2.png|thumb|'''''Fig 6. '''''Primative cell volume as a function of temperature for the QMA (green) and using MD (red)]]<br />
The thermal expansion coefficient, <math>\textstyle \alpha_V </math>, was calculated using by comparison of the cell volumes over temperatures ranging from 100 K to 1000 K; this can be seen in '''''Fig 6. '''''From this, it can be seen that while MD shows linearity all the way down to 100 K, the QMA deviates at temperatures below</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553515Rep:Mod:js6313 mgo2016-03-11T11:17:53Z<p>Js6313: </p>
<hr />
<div>=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''[[File:Js6313Lattice primative.png|thumb|'''''Fig 1. '''''The primative unit cell of MgO]]<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. <br />
<br />
When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges. The unit cell is a parallelepiped with <br />
<br />
<math>\textstyle \alpha = \beta = \gamma = 60^{\circ} </math> and <math>\textstyle a = b = c = 2.9783 \AA </math> which gives a conventional lattice constant of <math>\textstyle 2.10597 \AA </math>.<br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== Overview of The First Brillouin Zone, Dispersion Curves and the Density of States (DOS) ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|'''''Fig 2. '''''Dispersion curve taken with a grid size of 50]]<br />
In a crystal with lattice vector ''a'', the wavevector <math>\textstyle k = \frac{2 \pi}{\lambda} </math>. Only wavevectors between <math>\textstyle -\frac{\pi}{a} \text{ and } \frac{\pi}{a} </math> are unique; these are in what is known as the first Brillouin zone. Every k-point in this domain is unique, and has it's own energy level which under ideality is given by:<br />
<br />
<math> E(k) = \frac{\hbar^2 k^2}{2m} + E_{pot} </math><br />
<br />
Note that this only holds true within the stated boundary conditions for the first Brillouin zone.<br />
<br />
Dispersion curves are a plot of frequency against k, within the first Brillouin zone. Each line is a phonon - at any point, a maximum of six phonons is observed. This is because there are 2 atoms in each unit cell, each with 3 vibrational degrees of freedom. At some points (where lines converge), phonons are degenerate.<br />
[[File:Js6313dos.png|thumb|'''''Fig 3. '''''Density of phonon states along k-point L (left) and along 50 k-points (right)]]<br />
A density of states (DOS) graph shows the frequency of phonons along a given number of k-points; setting this to 1 gives a view along the L k-point shown below. Note that the convergence of the lower energy curves corresponds to a higher DOS. As the number of sampled k-points increases, a more representative view of the density of states is gained.<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation vs Molecular Dynamics ===<br />
For the quasi-harmonic approximation, Every MgO bond can be thought of as a mass on a spring, the motion of which is governed by Hooke's law:<br />
<br />
<math>F = -k(x - \bar{x})</math><br />
<br />
This is a good approximation at very low temperatures, where bond lengths do not change significantly, however it breaks down once vibrations start to stretch bonds closer to their breaking points. The QMA is a variation on this model whereby k is dependent on the volume, and that the and that the harmonic approximation works for all volumes.<br />
<br />
Molecular dynamics takes a set of starting conditions, and takes a series of snapshots at regular intervals (5 ps in this instance), until the properties reach equilibrium. A random selection of velocities is selected such that the overall temperature of the 32 unit cells is that specified beforehand; from here, the forces are computed, and then:<br />
<br />
<math>a = F/m</math><br />
<br />
<math>v_{\text{new}} = v_{\text{old}} + a \delta t</math><br />
<br />
<math>r_{\text{new}} = r_{\text{old}} + v_{\text{new}} \delta t<br />
<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly set the temperature of 300 K and the pressure to 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brillouin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation. A comparison of free energies can be seen in '''''Fig 5. '''''Accuracies to within 1 meV, 0.5 meV and 0.1 meV require grid sizes of 2, 2 and 3 respectively.<br />
<br />
=== Thermal Expansion of MgO Using the Quasi-Harmonic Approximation and Molecular Dynamics ===<br />
[[File:Js6313Lattice cube.png|left|thumb|'''''Fig 4. '''''The conventional unit cell for MgO]]<br />
[[File:Temperturevsfreeenergy.png|thumb|'''''Fig 5. '''''Graph showing temperature vs free energy for the QMA (green) and MD (red)]]<br />
The calculations were run at temperatures ranging from 0 K to 6000 K; past 2000 K when using the QMA, DLV gave error messages, possibly due to the dissociation of the lattice. The results are summarised in '''''Fig 5. '''''The two models give similar results at low temperatures, but diverge at T > 450 K. This is because anharmonicity increaces with increacing temperature, so the QM model breaks down (it predicts that free energy will fall with tempreature, whereas the MD model uses the formula <math>\textstyle E_k = \frac{k_B T}{2} </math>. This means that kinetic energy of the individual ions and thus the free energy increases with temperature, as shown in the graph and proven by experiment.<ref><nowiki>http://www.asminternational.org/documents/10192/3449965/ACFAAD6.pdf/2d574bfc-e104-48c5-8d8e-c33ffe91c3a0</nowiki></ref><br />
<br />
==== The Thermal Expansion Coefficient ====<br />
[[File:Volume vs temperature2.png|thumb|'''''Fig 6. '''''Primative cell volume as a function of temperature for the QMA (green) and using MD (red)]]<br />
The thermal expansion coefficient, <math>\textstyle \alpha_V </math>, was calculated using by comparison of the cell volumes over temperatures ranging from 100 K to 1000 K; this can be seen in '''''Fig 6. '''''From this, it can be seen that while MD shows linearity all the way down to 100 K, the QMA deviates at temperatures below</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553510Rep:Mod:js6313 mgo2016-03-11T11:15:23Z<p>Js6313: </p>
<hr />
<div>=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''[[File:Js6313Lattice primative.png|thumb|'''''Fig 1. '''''The primative unit cell of MgO]]<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. <br />
<br />
When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges. The unit cell is a parallelepiped with <br />
<br />
<math>\textstyle \alpha = \beta = \gamma = 60^{\circ} </math> and <math>\textstyle a = b = c = 2.9783 \AA </math> which gives a conventional lattice constant of <math>\textstyle 2.10597 \AA </math>.<br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== Overview of The First Brillouin Zone, Dispersion Curves and the Density of States (DOS) ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|'''''Fig 2. '''''Dispersion curve taken with a grid size of 50]]<br />
In a crystal with lattice vector ''a'', the wavevector <math>\textstyle k = \frac{2 \pi}{\lambda} </math>. Only wavevectors between <math>\textstyle -\frac{\pi}{a} \text{ and } \frac{\pi}{a} </math> are unique; these are in what is known as the first Brillouin zone. Every k-point in this domain is unique, and has it's own energy level which under ideality is given by:<br />
<br />
<math> E(k) = \frac{\hbar^2 k^2}{2m} + E_{pot} </math><br />
<br />
Note that this only holds true within the stated boundary conditions for the first Brillouin zone.<br />
<br />
Dispersion curves are a plot of frequency against k, within the first Brillouin zone. Each line is a phonon - at any point, a maximum of six phonons is observed. This is because there are 2 atoms in each unit cell, each with 3 vibrational degrees of freedom. At some points (where lines converge), phonons are degenerate.<br />
[[File:Js6313dos.png|thumb|'''''Fig 3. '''''Density of phonon states along k-point L (left) and along 50 k-points (right)]]<br />
A density of states (DOS) graph shows the frequency of phonons along a given number of k-points; setting this to 1 gives a view along the L k-point shown below. Note that the convergence of the lower energy curves corresponds to a higher DOS. As the number of sampled k-points increases, a more representative view of the density of states is gained.<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation vs Molecular Dynamics ===<br />
For the quasi-harmonic approximation, Every MgO bond can be thought of as a mass on a spring, the motion of which is governed by Hooke's law:<br />
<br />
<math>F = -k(x - \bar{x})</math><br />
<br />
This is a good approximation at very low temperatures, where bond lengths do not change significantly, however it breaks down once vibrations start to stretch bonds closer to their breaking points. The QMA is a variation on this model whereby k is dependent on the volume, and that the and that the harmonic approximation works for all volumes.<br />
<br />
Molecular dynamics takes a set of starting conditions, and takes a series of snapshots at regular intervals (5 ps in this instance), until the properties reach equilibrium. A random selection of velocities is selected such that the overall temperature of the 32 unit cells is that specified beforehand; from here, the forces are computed, and then:<br />
<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly set the temperature of 300 K and the pressure to 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brillouin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation. A comparison of free energies can be seen in '''''Fig 5. '''''Accuracies to within 1 meV, 0.5 meV and 0.1 meV require grid sizes of 2, 2 and 3 respectively.<br />
<br />
=== Thermal Expansion of MgO Using the Quasi-Harmonic Approximation and Molecular Dynamics ===<br />
[[File:Js6313Lattice cube.png|left|thumb|'''''Fig 4. '''''The conventional unit cell for MgO]]<br />
[[File:Temperturevsfreeenergy.png|thumb|'''''Fig 5. '''''Graph showing temperature vs free energy for the QMA (green) and MD (red)]]<br />
The calculations were run at temperatures ranging from 0 K to 6000 K; past 2000 K when using the QMA, DLV gave error messages, possibly due to the dissociation of the lattice. The results are summarised in '''''Fig 5. '''''The two models give similar results at low temperatures, but diverge at T > 450 K. This is because anharmonicity increaces with increacing temperature, so the QM model breaks down (it predicts that free energy will fall with tempreature, whereas the MD model uses the formula <math>\textstyle E_k = \frac{k_B T}{2} </math>. This means that kinetic energy of the individual ions and thus the free energy increases with temperature, as shown in the graph and proven by experiment.<ref><nowiki>http://www.asminternational.org/documents/10192/3449965/ACFAAD6.pdf/2d574bfc-e104-48c5-8d8e-c33ffe91c3a0</nowiki></ref><br />
<br />
==== The Thermal Expansion Coefficient ====<br />
[[File:Volume vs temperature2.png|thumb|'''''Fig 6. '''''Primative cell volume as a function of temperature for the QMA (green) and using MD (red)]]<br />
The thermal expansion coefficient, <math>\textstyle \alpha_V </math>, was calculated using by comparison of the cell volumes over temperatures ranging from 100 K to 1000 K; this can be seen in '''''Fig 6. '''''From this, it can be seen that while MD shows linearity all the way down to 100 K, the QMA deviates at temperatures below</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553496Rep:Mod:js6313 mgo2016-03-11T11:05:18Z<p>Js6313: /* The Thermal Expansion Coefficient */</p>
<hr />
<div>= Using DLVisualize and GULP to Explore the Free Energy and Thermal Expansion of MgO =<br />
<br />
=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''[[File:Js6313Lattice primative.png|thumb|'''''Fig 1. '''''The primative unit cell of MgO]]<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. <br />
<br />
When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges. The unit cell is a parallelepiped with <br />
<br />
<math>\textstyle \alpha = \beta = \gamma = 60^{\circ} </math> and <math>\textstyle a = b = c = 2.9783 \AA </math> which gives a conventional lattice constant of <math>\textstyle 2.10597 \AA </math>.<br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== Overview of The First Brillouin Zone, Dispersion Curves and the Density of States (DOS) ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|'''''Fig 2. '''''Dispersion curve taken with a grid size of 50]]<br />
In a crystal with lattice vector ''a'', the wavevector <math>\textstyle k = \frac{2 \pi}{\lambda} </math>. Only wavevectors between <math>\textstyle -\frac{\pi}{a} \text{ and } \frac{\pi}{a} </math> are unique; these are in what is known as the first Brillouin zone. Every k-point in this domain is unique, and has it's own energy level which under ideality is given by:<br />
<br />
<math> E(k) = \frac{\hbar^2 k^2}{2m} + E_{pot} </math><br />
<br />
Note that this only holds true within the stated boundary conditions for the first Brillouin zone.<br />
<br />
Dispersion curves are a plot of frequency against k, within the first Brillouin zone. Each line is a phonon - at any point, a maximum of six phonons is observed. This is because there are 2 atoms in each unit cell, each with 3 vibrational degrees of freedom. At some points (where lines converge), phonons are degenerate.<br />
[[File:Js6313dos.png|thumb|'''''Fig 3. '''''Density of phonon states along k-point L (left) and along 50 k-points (right)]]<br />
A density of states (DOS) graph shows the frequency of phonons along a given number of k-points; setting this to 1 gives a view along the L k-point shown below. Note that the convergence of the lower energy curves corresponds to a higher DOS. As the number of sampled k-points increases, a more representative view of the density of states is gained.<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation ===<br />
Every MgO bond can be thought of as a mass on a spring, the motion of which is governed by Hooke's law:<br />
<br />
<math>F = -k(x - \bar{x})</math><br />
<br />
This is a good approximation at very low temperatures, where bond lengths do not change significantly, however it breaks down once vibrations start to stretch bonds closer to their breaking points. The QMA is a variation on this model whereby k is dependent on the volume, and that the and that the harmonic approximation works for all volumes.<br />
<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly set the temperature of 300 K and the pressure to 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brillouin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation. A comparison of free energies can be seen in '''''Fig 5. '''''Accuracies to within 1 meV, 0.5 meV and 0.1 meV require grid sizes of 2, 2 and 3 respectively.<br />
<br />
=== Thermal Expansion of MgO Using the Quasi-Harmonic Approximation and Molecular Dynamics ===<br />
[[File:Js6313Lattice cube.png|left|thumb|'''''Fig 4. '''''The conventional unit cell for MgO]]<br />
[[File:Temperturevsfreeenergy.png|thumb|'''''Fig 5. '''''Graph showing temperature vs free energy for the QMA (green) and MD (red)]]<br />
The calculations were run at temperatures ranging from 0 K to 6000 K; past 2000 K when using the QMA, DLV gave error messages, possibly due to the dissociation of the lattice. The results are summarised in '''''Fig 5. '''''The two models give similar results at low temperatures, but diverge at T > 450 K. This is because anharmonicity increaces with increacing temperature, so the QM model breaks down (it predicts that free energy will fall with tempreature, whereas the MD model uses the formula <math>\textstyle E_k = \frac{k_B T}{2} </math>. This means that kinetic energy of the individual ions and thus the free energy increases with temperature, as shown in the graph and proven by experiment.<ref><nowiki>http://www.asminternational.org/documents/10192/3449965/ACFAAD6.pdf/2d574bfc-e104-48c5-8d8e-c33ffe91c3a0</nowiki></ref><br />
<br />
==== The Thermal Expansion Coefficient ====<br />
[[File:Volume vs temperature2.png|thumb|'''''Fig 6. '''''Primative cell volume as a function of temperature for the QMA (green) and using MD (red)]]<br />
The thermal expansion coefficient, <math>\textstyle \alpha_V </math>, was calculated using by comparison of the cell volumes over temperatures ranging from 100 K to 1000 K; this can be seen in '''''Fig 6. '''''From this, it can be seen that</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=File:Volume_vs_temperature2.png&diff=553492File:Volume vs temperature2.png2016-03-11T11:04:40Z<p>Js6313: </p>
<hr />
<div></div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=File:Volume_vs_temperature.png&diff=553487File:Volume vs temperature.png2016-03-11T11:03:03Z<p>Js6313: </p>
<hr />
<div></div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553472Rep:Mod:js6313 mgo2016-03-11T10:50:10Z<p>Js6313: </p>
<hr />
<div>= Using DLVisualize and GULP to Explore the Free Energy and Thermal Expansion of MgO =<br />
<br />
=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''[[File:Js6313Lattice primative.png|thumb|'''''Fig 1. '''''The primative unit cell of MgO]]<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. <br />
<br />
When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges. The unit cell is a parallelepiped with <br />
<br />
<math>\textstyle \alpha = \beta = \gamma = 60^{\circ} </math> and <math>\textstyle a = b = c = 2.9783 \AA </math> which gives a conventional lattice constant of <math>\textstyle 2.10597 \AA </math>.<br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== Overview of The First Brillouin Zone, Dispersion Curves and the Density of States (DOS) ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|'''''Fig 2. '''''Dispersion curve taken with a grid size of 50]]<br />
In a crystal with lattice vector ''a'', the wavevector <math>\textstyle k = \frac{2 \pi}{\lambda} </math>. Only wavevectors between <math>\textstyle -\frac{\pi}{a} \text{ and } \frac{\pi}{a} </math> are unique; these are in what is known as the first Brillouin zone. Every k-point in this domain is unique, and has it's own energy level which under ideality is given by:<br />
<br />
<math> E(k) = \frac{\hbar^2 k^2}{2m} + E_{pot} </math><br />
<br />
Note that this only holds true within the stated boundary conditions for the first Brillouin zone.<br />
<br />
Dispersion curves are a plot of frequency against k, within the first Brillouin zone. Each line is a phonon - at any point, a maximum of six phonons is observed. This is because there are 2 atoms in each unit cell, each with 3 vibrational degrees of freedom. At some points (where lines converge), phonons are degenerate.<br />
[[File:Js6313dos.png|thumb|'''''Fig 3. '''''Density of phonon states along k-point L (left) and along 50 k-points (right)]]<br />
A density of states (DOS) graph shows the frequency of phonons along a given number of k-points; setting this to 1 gives a view along the L k-point shown below. Note that the convergence of the lower energy curves corresponds to a higher DOS. As the number of sampled k-points increases, a more representative view of the density of states is gained.<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation ===<br />
Every MgO bond can be thought of as a mass on a spring, the motion of which is governed by Hooke's law:<br />
<br />
<math>F = -k(x - \bar{x})</math><br />
<br />
This is a good approximation at very low temperatures, where bond lengths do not change significantly, however it breaks down once vibrations start to stretch bonds closer to their breaking points. The QMA is a variation on this model whereby k is dependent on the volume, and that the and that the harmonic approximation works for all volumes.<br />
<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly set the temperature of 300 K and the pressure to 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brillouin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation. A comparison of free energies can be seen in '''''Fig 5. '''''Accuracies to within 1 meV, 0.5 meV and 0.1 meV require grid sizes of 2, 2 and 3 respectively.<br />
<br />
=== Thermal Expansion of MgO Using the Quasi-Harmonic Approximation and Molecular Dynamics ===<br />
[[File:Js6313Lattice cube.png|left|thumb|'''''Fig 4. '''''The conventional unit cell for MgO]]<br />
[[File:Temperturevsfreeenergy.png|thumb|'''''Fig 5. '''''Graph showing temperature vs free energy for the QMA (green) and MD (red)]]<br />
The calculations were run at temperatures ranging from 0 K to 6000 K; past 2000 K when using the QMA, DLV gave error messages, possibly due to the dissociation of the lattice. The results are summarised in '''''Fig 5. '''''The two models give similar results at low temperatures, but diverge at T > 450 K. This is because anharmonicity increaces with increacing temperature, so the QM model breaks down (it predicts that free energy will fall with tempreature, whereas the MD model uses the formula <math>\textstyle E_k = \frac{k_B T}{2} </math>. This means that kinetic energy of the individual ions and thus the free energy increases with temperature, as shown in the graph and proven by experiment.<ref><nowiki>http://www.asminternational.org/documents/10192/3449965/ACFAAD6.pdf/2d574bfc-e104-48c5-8d8e-c33ffe91c3a0</nowiki></ref><br />
<br />
==== The Thermal Expansion Coefficient ====<br />
The thermal expansion coefficient, <math>\textstyle \alpha_V </math>, was calculated using by comparison of the cell volumes over temperatures ranging from 100 K to 1000 K; this can be seen in '''''Fig 6.'''''.</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553471Rep:Mod:js6313 mgo2016-03-11T10:48:23Z<p>Js6313: </p>
<hr />
<div>= Using DLVisualize and GULP to Explore the Free Energy and Thermal Expansion of MgO =<br />
<br />
=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''[[File:Js6313Lattice primative.png|thumb|'''''Fig 1. '''''The primative unit cell of MgO]]<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. <br />
<br />
When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges. The unit cell is a parallelepiped with <br />
<br />
<math>\textstyle \alpha = \beta = \gamma = 60^{\circ} </math> and <math>\textstyle a = b = c = 2.9783 \AA </math> which gives a conventional lattice constant of <math>\textstyle 2.10597 \AA </math>.<br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== Overview of The First Brillouin Zone, Dispersion Curves and the Density of States (DOS) ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|'''''Fig 2. '''''Dispersion curve taken with a grid size of 50]]<br />
In a crystal with lattice vector ''a'', the wavevector <math>\textstyle k = \frac{2 \pi}{\lambda} </math>. Only wavevectors between <math>\textstyle -\frac{\pi}{a} \text{ and } \frac{\pi}{a} </math> are unique; these are in what is known as the first Brillouin zone. Every k-point in this domain is unique, and has it's own energy level which under ideality is given by:<br />
<br />
<math> E(k) = \frac{\hbar^2 k^2}{2m} + E_{pot} </math><br />
<br />
Note that this only holds true within the stated boundary conditions for the first Brillouin zone.<br />
<br />
Dispersion curves are a plot of frequency against k, within the first Brillouin zone. Each line is a phonon - at any point, a maximum of six phonons is observed. This is because there are 2 atoms in each unit cell, each with 3 vibrational degrees of freedom. At some points (where lines converge), phonons are degenerate.<br />
[[File:Js6313dos.png|thumb|'''''Fig 3. '''''Density of phonon states along k-point L (left) and along 50 k-points (right)]]<br />
A density of states (DOS) graph shows the frequency of phonons along a given number of k-points; setting this to 1 gives a view along the L k-point shown below. Note that the convergence of the lower energy curves corresponds to a higher DOS. As the number of sampled k-points increases, a more representative view of the density of states is gained.<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation ===<br />
Every MgO bond can be thought of as a mass on a spring, the motion of which is governed by Hooke's law:<br />
<br />
<math>F = -k(x - \bar{x})</math><br />
<br />
This is a good approximation at very low temperatures, where bond lengths do not change significantly, however it breaks down once vibrations start to stretch bonds closer to their breaking points. The QMA is a variation on this model whereby k is dependent on the volume, and that the and that the harmonic approximation works for all volumes.<br />
<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly set the temperature of 300 K and the pressure to 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brillouin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation. A comparison of free energies can be seen in '''''Fig 5. '''''Accuracies to within 1 meV, 0.5 meV and 0.1 meV require grid sizes of 2, 2 and 3 respectively.<br />
<br />
=== Thermal Expansion of MgO Using the Quasi-Harmonic Approximation and Molecular Dynamics ===<br />
[[File:Js6313Lattice cube.png|left|thumb|'''''Fig 4. '''''The conventional unit cell for MgO]]<br />
[[File:Temperturevsfreeenergy.png|thumb|'''''Fig 5. '''''Graph showing temperature vs free energy for the QMA (green) and MD (red)]]<br />
The calculations were run at temperatures ranging from 0 K to 6000 K; past 2000 K when using the QMA, DLV gave error messages, possibly due to the dissociation of the lattice. The results are summarised in '''''Fig 5. '''''The two models give similar results at low temperatures, but diverge at T > 450 K. This is because anharmonicity increaces with increacing temperature, so the QM model breaks down (it predicts that free energy will fall with tempreature, whereas the MD model uses the formula <math>\textstyle E_k = \frac{k_B T}{2} </math>. This means that kinetic energy of the individual ions and thus the free energy increases with temperature, as shown in the graph and proven by experiment.<ref><nowiki>http://www.asminternational.org/documents/10192/3449965/ACFAAD6.pdf/2d574bfc-e104-48c5-8d8e-c33ffe91c3a0</nowiki></ref><br />
<br />
==== The Thermal Expansion Coefficient ====<br />
The thermal expansion coefficient,</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553470Rep:Mod:js6313 mgo2016-03-11T10:46:37Z<p>Js6313: </p>
<hr />
<div>= Using DLVisualize and GULP to Explore the Free Energy and Thermal Expansion of MgO =<br />
<br />
=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''[[File:Js6313Lattice primative.png|thumb|'''''Fig 1. '''''The primative unit cell of MgO]]<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. <br />
<br />
When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges. The unit cell is a parallelepiped with <br />
<br />
<math>\textstyle \alpha = \beta = \gamma = 60^{\circ} </math> and <math>\textstyle a = b = c = 2.9783 \AA </math> which gives a conventional lattice constant of <math>\textstyle 2.10597 \AA </math>.<br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== Overview of The First Brillouin Zone, Dispersion Curves and the Density of States (DOS) ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|'''''Fig 2. '''''Dispersion curve taken with a grid size of 50]]<br />
In a crystal with lattice vector ''a'', the wavevector <math>\textstyle k = \frac{2 \pi}{\lambda} </math>. Only wavevectors between <math>\textstyle -\frac{\pi}{a} \text{ and } \frac{\pi}{a} </math> are unique; these are in what is known as the first Brillouin zone. Every k-point in this domain is unique, and has it's own energy level which under ideality is given by:<br />
<br />
<math> E(k) = \frac{\hbar^2 k^2}{2m} + E_{pot} </math><br />
<br />
Note that this only holds true within the stated boundary conditions for the first Brillouin zone.<br />
<br />
Dispersion curves are a plot of frequency against k, within the first Brillouin zone. Each line is a phonon - at any point, a maximum of six phonons is observed. This is because there are 2 atoms in each unit cell, each with 3 vibrational degrees of freedom. At some points (where lines converge), phonons are degenerate.<br />
[[File:Js6313dos.png|thumb|'''''Fig 3. '''''Density of phonon states along k-point L (left) and along 50 k-points (right)]]<br />
A density of states (DOS) graph shows the frequency of phonons along a given number of k-points; setting this to 1 gives a view along the L k-point shown below. Note that the convergence of the lower energy curves corresponds to a higher DOS. As the number of sampled k-points increases, a more representative view of the density of states is gained.<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation ===<br />
Every MgO bond can be thought of as a mass on a spring, the motion of which is governed by Hooke's law:<br />
<br />
<math>F = -k(x - \bar{x})</math><br />
<br />
This is a good approximation at very low temperatures, where bond lengths do not change significantly, however it breaks down once vibrations start to stretch bonds closer to their breaking points. The QMA is a variation on this model whereby k is dependent on the volume, and that the and that the harmonic approximation works for all volumes.<br />
<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly set the temperature of 300 K and the pressure to 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brillouin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation. A comparison of free energies can be seen in '''''Fig 5. '''''Accuracies to within 1 meV, 0.5 meV and 0.1 meV require grid sizes of 2, 2 and 3 respectively.<br />
<br />
=== Thermal Expansion of MgO Using the Quasi-Harmonic Approximation and Molecular Dynamics ===<br />
[[File:Js6313Lattice cube.png|left|thumb|'''''Fig 4. '''''The conventional unit cell for MgO]]<br />
[[File:Temperturevsfreeenergy.png|thumb|'''''Fig 5. '''''Graph showing temperature vs free energy for the QMA (green) and MD (red)]]<br />
The calculations were run at temperatures ranging from 0 K to 6000 K; past 2000 K when using the QMA, DLV gave error messages, possibly due to the dissociation of the lattice. The results are summarised in '''''Fig 5. '''''The two models give similar results at low temperatures, but diverge at T > 450 K. This is because anharmonicity increaces with increacing temperature, so the QM model breaks down (it predicts that free energy will fall with tempreature, whereas the MD model uses the formula <math>\textstyle E_k = \frac{k_B T}{2} </math>. This means that kinetic energy of the individual ions and thus the free energy increases with temperature, as shown in the graph and proven by experiment.</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553468Rep:Mod:js6313 mgo2016-03-11T10:42:32Z<p>Js6313: </p>
<hr />
<div>= Using DLVisualize and GULP to Explore the Free Energy and Thermal Expansion of MgO =<br />
<br />
=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''[[File:Js6313Lattice primative.png|thumb|'''''Fig 1. '''''The primative unit cell of MgO]]<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. <br />
<br />
When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges. The unit cell is a parallelepiped with <br />
<br />
<math>\textstyle \alpha = \beta = \gamma = 60^{\circ} </math> and <math>\textstyle a = b = c = 2.9783 \AA </math> which gives a conventional lattice constant of <math>\textstyle 2.10597 \AA </math>.<br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== Overview of The First Brillouin Zone, Dispersion Curves and the Density of States (DOS) ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|'''''Fig 2. '''''Dispersion curve taken with a grid size of 50]]<br />
In a crystal with lattice vector ''a'', the wavevector <math>\textstyle k = \frac{2 \pi}{\lambda} </math>. Only wavevectors between <math>\textstyle -\frac{\pi}{a} \text{ and } \frac{\pi}{a} </math> are unique; these are in what is known as the first Brillouin zone. Every k-point in this domain is unique, and has it's own energy level which under ideality is given by:<br />
<br />
<math> E(k) = \frac{\hbar^2 k^2}{2m} + E_{pot} </math><br />
<br />
Note that this only holds true within the stated boundary conditions for the first Brillouin zone.<br />
<br />
Dispersion curves are a plot of frequency against k, within the first Brillouin zone. Each line is a phonon - at any point, a maximum of six phonons is observed. This is because there are 2 atoms in each unit cell, each with 3 vibrational degrees of freedom. At some points (where lines converge), phonons are degenerate.<br />
[[File:Js6313dos.png|thumb|'''''Fig 3. '''''Density of phonon states along k-point L (left) and along 50 k-points (right)]]<br />
A density of states (DOS) graph shows the frequency of phonons along a given number of k-points; setting this to 1 gives a view along the L k-point shown below. Note that the convergence of the lower energy curves corresponds to a higher DOS. As the number of sampled k-points increases, a more representative view of the density of states is gained.<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation ===<br />
Every MgO bond can be thought of as a mass on a spring, the motion of which is governed by Hooke's law:<br />
<br />
<math>F = -k(x - \bar{x})</math><br />
<br />
This is a good approximation at very low temperatures, where bond lengths do not change significantly, however it breaks down once vibrations start to stretch bonds closer to their breaking points. The QMA is a variation on this model whereby k is dependent on the volume, and that the and that the harmonic approximation works for all volumes.<br />
<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly set the temperature of 300 K and the pressure to 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brillouin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation. A comparison of free energies can be seen in '''''Fig 5. '''''Accuracies to within 1 meV, 0.5 meV and 0.1 meV require grid sizes of 2, 2 and 3 respectively.<br />
<br />
=== Thermal Expansion of MgO Using the Quasi-Harmonic Approximation and Molecular Dynamics ===<br />
[[File:Js6313Lattice cube.png|left|thumb|'''''Fig 4. '''''The conventional unit cell for MgO]]<br />
[[File:Temperturevsfreeenergy.png|thumb|'''''Fig 5. '''''Graph showing temperature vs free energy for the QMA (green) and MD (red)]]<br />
The calculations were run at temperatures ranging from 0 K to 6000 K; past 2000 K when using the QMA, DLV gave error messages, possibly due to the dissociation of the lattice. The results are summarised in '''''Fig 5. '''''The two models give similar results at low temperatures, but diverge at T > 450 K. This is because anharmonicity increaces with increacing temperature, so the QM model breaks down (it predicts that free energy will fall with tempreature, whereas the MD model uses the formula</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=File:Temperturevsfreeenergy.png&diff=553462File:Temperturevsfreeenergy.png2016-03-11T10:35:13Z<p>Js6313: </p>
<hr />
<div></div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553444Rep:Mod:js6313 mgo2016-03-11T09:58:04Z<p>Js6313: </p>
<hr />
<div>= Using DLVisualize and GULP to Explore the Free Energy and Thermal Expansion of MgO =<br />
<br />
=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''[[File:Js6313Lattice primative.png|thumb|The primative unit cell of MgO]]<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. <br />
<br />
When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges. The unit cell is a parallelepiped with <br />
<br />
<math>\textstyle \alpha = \beta = \gamma = 60^{\circ} </math> and <math>\textstyle a = b = c = 2.9783 \AA </math> which gives a conventional lattice constant of <math>\textstyle 2.10597 \AA </math>.<br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== Overview of The First Brillouin Zone, Dispersion Curves and the Density of States (DOS) ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|Dispersion curve taken with a grid size of 50]]<br />
In a crystal with lattice vector ''a'', the wavevector <math>\textstyle k = \frac{2 \pi}{\lambda} </math>. Only wavevectors between <math>\textstyle -\frac{\pi}{a} \text{ and } \frac{\pi}{a} </math> are unique; these are in what is known as the first Brillouin zone. Every k-point in this domain is unique, and has it's own energy level which under ideality is given by:<br />
<br />
<math> E(k) = \frac{\hbar^2 k^2}{2m} + E_{pot} </math><br />
<br />
Note that this only holds true within the stated boundary conditions for the first Brillouin zone.<br />
<br />
Dispersion curves are a plot of frequency against k, within the first Brillouin zone. Each line is a phonon - at any point, a maximum of six phonons is observed. This is because there are 2 atoms in each unit cell, each with 3 vibrational degrees of freedom. At some points (where lines converge), phonons are degenerate.<br />
[[File:Js6313dos.png|thumb|Density of phonon states along k-point L (left) and along 50 k-points (right)]]<br />
A density of states (DOS) graph shows the frequency of phonons along a given number of k-points; setting this to 1 gives a view along the L k-point shown below. Note that the convergence of the lower energy curves corresponds to a higher DOS. As the number of sampled k-points increases, a more representative view of the density of states is gained.<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation ===<br />
Every MgO bond can be thought of as a mass on a spring, the motion of which is governed by Hooke's law:<br />
<br />
<math>F = -k(x - \bar{x})</math><br />
<br />
This is a good approximation at very low temperatures, where bond lengths do not change significantly, however it breaks down once vibrations start to stretch bonds closer to their breaking points. The QMA is a variation on this model whereby k is dependent on the volume, and that the and that the harmonic approximation works for all volumes.<br />
<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly set the temperature of 300 K and the pressure to 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brillouin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation.<br />
<br />
=== Thermal Expansion of MgO Using the Quasi-Harmonic Approximation and Molecular Dynamics ===<br />
[[File:Js6313Lattice cube.png|left|thumb|The conventional unit cell for MgO]]The calculations were run at temperatures ranging from 0 K to 6000 K; past 2000 K when using the QMA, DLV gave error messages, possibly due to the dissociation of the lattice. The results are summarised in the graph below:</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=File:Js6313dos.png&diff=553442File:Js6313dos.png2016-03-11T09:53:27Z<p>Js6313: </p>
<hr />
<div></div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=File:Js6313002_phonon_DOS_shrink1.png&diff=553436File:Js6313002 phonon DOS shrink1.png2016-03-11T09:45:08Z<p>Js6313: </p>
<hr />
<div></div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553419Rep:Mod:js6313 mgo2016-03-11T09:30:38Z<p>Js6313: </p>
<hr />
<div>= Using DLVisualize and GULP to Explore the Free Energy and Thermal Expansion of MgO =<br />
<br />
=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''[[File:Js6313Lattice primative.png|thumb|The primative unit cell of MgO]]<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. <br />
<br />
When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges. The unit cell is a parallelepiped with <math>\textstyle \alpha = \beta = \gamma = 60^{\circ} </math> and <math>\textstyle a = b = c = 2.9783 \AA </math> which gives a conventional lattice constant of <math>\textstyle 2.10597 \AA </math>.<br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== The First Brillouin Zone ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|Dispersion curve taken with a grid size of 50]]<br />
In a crystal with lattice vector ''a'', the wavevector <math>\textstyle k = \frac{2 \pi}{\lambda} </math>. Only wavevectors between <math>\textstyle -\frac{\pi}{a} \text{ and } \frac{\pi}{a} </math> are unique; these are in what is known as the first Brillouin zone. Every k-point in this domain is unique, and has it's own energy level which under ideality is given by:<br />
<br />
<math> E(k) = \frac{\hbar^2 k^2}{2m} + E_{pot} </math><br />
<br />
Note that this only holds true within the stated boundary conditions for the first Brillouin zone.<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation ===<br />
[[File:Js6313Lattice cube.png|left|thumb|The conventional unit cell for MgO]]Every MgO bond can be thought of as a mass on a spring, the motion of which is governed by Hooke's law:<br />
<br />
<math>F = -k(x - \bar{x})</math><br />
<br />
This is a good approximation at very low temperatures, where bond lengths do not change significantly, however it breaks down once vibrations start to stretch bonds closer to their breaking points. The QMA is a variation on this model whereby k is dependent on the volume, and that the and that the harmonic approximation works for all volumes.<br />
<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly set the temperature of 300 K and the pressure to 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brillouin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation.</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553418Rep:Mod:js6313 mgo2016-03-11T09:29:58Z<p>Js6313: </p>
<hr />
<div>= Using DLVisualize and GULP to Explore the Free Energy and Thermal Expansion of MgO =<br />
<br />
=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''[[File:Js6313Lattice primative.png|thumb|The primative unit cell of MgO]]<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. <br />
<br />
When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges. The unit cell is a parallelepiped with <math>\textstyle \alpha = \beta = \gamma = 60^{\circ} </math> and <math>\textstyle a = b = c = 2.9783 \AA </math> which gives a conventional lattice constant of <math>\textstyle 2.10597 \AA </math>.<br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== The First Brillouin Zone ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|Dispersion curve taken with a grid size of 50]]<br />
In a crystal with lattice vector ''a'', the wavevector <math>\textstyle k = \frac{2 \pi}{\lambda} </math>. Only wavevectors between <math>\textstyle -\frac{\pi}{a} \text{ and } \frac{\pi}{a} </math> are unique; these are in what is known as the first Brillouin zone. Every k-point in this domain is unique, and has it's own energy level which under ideality is given by:<br />
<br />
<math> E(k) = \frac{\hbar^2 k^2}{2m} + E_{pot} </math><br />
<br />
Note that this only holds true within the stated boundary conditions for the first Brillouin zone.<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation ===<br />
[[File:Js6313Lattice cube.png|left|thumb|The conventional unit cell for MgO]]Every MgO bond can be thought of as a mass on a spring, the motion of which is governed by Hooke's law:<br />
<br />
<math>F = -k(x - \xbar)</math><br />
<br />
This is a good approximation at very low temperatures, where bond lengths do not change significantly, however it breaks down once vibrations start to stretch bonds closer to their breaking points. The QMA is a variation on this model whereby k is dependent on the volume, and that the and that the harmonic approximation works for all volumes.<br />
<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly set the temperature of 300 K and the pressure to 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brillouin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation.</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553411Rep:Mod:js6313 mgo2016-03-11T09:23:10Z<p>Js6313: </p>
<hr />
<div>= Using DLVisualize and GULP to Explore the Free Energy and Thermal Expansion of MgO =<br />
<br />
=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''[[File:Js6313Lattice primative.png|thumb|The primative unit cell of MgO]]<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. <br />
<br />
When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges. The unit cell is a parallelepiped with <math>\textstyle \alpha = \beta = \gamma = 60^{\circ} </math> and <math>\textstyle a = b = c = 2.9783 \AA </math> which gives a conventional lattice constant of <math>\textstyle 2.10597 \AA </math>.<br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== The First Brillouin Zone ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|Dispersion curve taken with a grid size of 50]]<br />
In a crystal with lattice vector ''a'', the wavevector <math>\textstyle k = \frac{2 \pi}{\lambda} </math>. Only wavevectors between <math>\textstyle -\frac{\pi}{a} \text{ and } \frac{\pi}{a} </math> are unique; these are in what is known as the first Brillouin zone. Every k-point in this domain is unique, and has it's own energy level which under ideality is given by:<br />
<br />
<math> E(k) = \frac{\hbar^2 k^2}{2m} + E_{pot} </math><br />
<br />
Note that this only holds true within the stated boundary conditions for the first Brillouin zone.<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation ===<br />
[[File:Js6313Lattice cube.png|left|thumb|The conventional unit cell for MgO]]Every MgO bond can be thought of as a mass on a spring, the motion of which is governed by Hooke's law:<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly set the temperature of 300 K and the pressure to 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brillouin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation.</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553402Rep:Mod:js6313 mgo2016-03-11T09:17:19Z<p>Js6313: /* The First Brouillin Zone */</p>
<hr />
<div>= Using DLVisualize and GULP to Explore the Free Energy and Thermal Expansion of MgO =<br />
<br />
=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''[[File:Js6313Lattice primative.png|thumb|The primative unit cell of MgO]]<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. [[File:Js6313Lattice cube.png|left|thumb|The conventional unit cell for MgO]]When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges. The unit cell is a parallelepiped with <math>\textstyle \alpha = \beta = \gamma = 60^{\circ} </math> and <math>\textstyle a = b = c = 2.9783 \AA </math> which gives a conventional lattice constant of <math>\textstyle 2.10597 \AA </math>.<br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== The First Brouillin Zone ===<br />
In a crystal with lattice vector ''a'', the wavevector <math>\textstyle k = \frac{2 \pi}{\lambda} </math>. Only wavevectors between <math>\textstyle -\frac{\pi}{a} \text{ and } \frac{\pi}{a} </math> are unique; these are in what is known as the first Brouillin zone. Every k-point in this domain is unique, and has it's own energy level which is given by:<br />
<br />
<math> E(k) = \frac{\hbar^2 k^2}{2m} + E_{pot} </math><br />
<br />
Note that this only holds true within the stated boundary conditions for the first Brouillin zone.<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|Dispersion curve taken with a grid size of 50]]<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly set the temperature of 300 K and the pressure to 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brouillin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation.</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553398Rep:Mod:js6313 mgo2016-03-11T09:15:24Z<p>Js6313: </p>
<hr />
<div>= Using DLVisualize and GULP to Explore the Free Energy and Thermal Expansion of MgO =<br />
<br />
=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''[[File:Js6313Lattice primative.png|thumb|The primative unit cell of MgO]]<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. [[File:Js6313Lattice cube.png|left|thumb|The conventional unit cell for MgO]]When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges. The unit cell is a parallelepiped with <math>\textstyle \alpha = \beta = \gamma = 60^{\circ} </math> and <math>\textstyle a = b = c = 2.9783 \AA </math> which gives a conventional lattice constant of <math>\textstyle 2.10597 \AA </math>.<br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== The First Brouillin Zone ===<br />
In a crystal with lattice vector ''a'', the wavevector <math>\textstyle k = \frac{2 \pi}{\lambda} </math>. Only wavevectors between <math>\textstyle -\frac{\pi}{a} \text{ and } \frac{\pi}{a} </math> are unique; these are in what is known as the first Brouillin zone. Every k-point in this domain is unique, and has it's own energy level which is given by:<br />
<br />
<math> E(k) = \frac{\hbar^2 k^2}{2m} + E_{pot} </math><br />
<br />
<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|Dispersion curve taken with a grid size of 50]]<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly set the temperature of 300 K and the pressure to 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brouillin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation.</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553396Rep:Mod:js6313 mgo2016-03-11T09:15:02Z<p>Js6313: </p>
<hr />
<div>= Using DLVisualize and GULP to Explore the Free Energy and Thermal Expansion of MgO =<br />
<br />
=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''[[File:Js6313Lattice primative.png|thumb|The primative unit cell of MgO]]<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. [[File:Js6313Lattice cube.png|left|thumb|The conventional unit cell for MgO]]When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges. The unit cell is a parallelepiped with <math>\textstyle \alpha = \beta = \gamma = 60^{\circ} </math> and <math>\textstyle a = b = c = 2.9783 \AA </math> which gives a conventional lattice constant of <math>\textstyle 2.10597 \AA </math>.<br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== The First Brouillin Zone ===<br />
In a crystal with lattice vector ''a'', the wavevector <math>\textstyle k = \frac{2 \pi}{\lambda} </math>. Only wavevectors between <math>\textstyle -\frac{pi}{a} \text{and} \frac{pi}{a} </math> are unique; these are in what is known as the first Brouillin zone. Every k-point in this domain is unique, and has it's own energy level which is given by:<br />
<br />
<math> E(k) = \frac{\hbar^2 k^2}{2m} + E_{pot} </math><br />
<br />
<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|Dispersion curve taken with a grid size of 50]]<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly set the temperature of 300 K and the pressure to 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brouillin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation.</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553395Rep:Mod:js6313 mgo2016-03-11T09:14:52Z<p>Js6313: </p>
<hr />
<div>= Using DLVisualize and GULP to Explore the Free Energy and Thermal Expansion of MgO =<br />
<br />
=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''[[File:Js6313Lattice primative.png|thumb|The primative unit cell of MgO]]<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. [[File:Js6313Lattice cube.png|left|thumb|The conventional unit cell for MgO]]When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges. The unit cell is a parallelepiped with <math>\textstyle \alpha = \beta = \gamma = 60^{\circ} </math> and <math>\textstyle a = b = c = 2.9783 \AA </math> which gives a conventional lattice constant of <math>\textstyle 2.10597 \AA </math>.<br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== The First Brouillin Zone ===<br />
In a crystal with lattice vector ''a'', the wavevector <math>\textstyle k = \frac{2 \pi}{\lambda} </math>. Only wavevectors between <math>\textstyle -\frac{pi}{a} \text{and} \frac{pi}{a} </math> are unique; these are in what is known as the first Brouillin zone. Every k-point in this domain is unique, and has it's own energy level which is given by:<br />
<br />
<math> E(k) = \frac{\hbar^2 k^2}{2m} + E_{pot}<br />
<br />
<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|Dispersion curve taken with a grid size of 50]]<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly set the temperature of 300 K and the pressure to 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brouillin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation.</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553387Rep:Mod:js6313 mgo2016-03-11T09:03:16Z<p>Js6313: </p>
<hr />
<div>= Using DLVisualize and GULP to Explore the Free Energy and Thermal Expansion of MgO =<br />
<br />
=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''[[File:Js6313Lattice primative.png|thumb|The primative unit cell of MgO]]<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. [[File:Js6313Lattice cube.png|left|thumb|The conventional unit cell for MgO]]When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges. The unit cell is a parallelepiped with <math>\textstyle \alpha = \beta = \gamma = 60^{\circ} </math> and <math>\textstyle a = b = c = 2.9783 \AA </math> which gives a conventional lattice constant of <math> 2.10597 \AA </math>.<br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== The First Brouillin Zone ===<br />
In a crystal with lattice vector ''a'', the wavevector <br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|Dispersion curve taken with a grid size of 50]]<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly set the temperature of 300 K and the pressure to 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brouillin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation.</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553385Rep:Mod:js6313 mgo2016-03-11T09:00:18Z<p>Js6313: </p>
<hr />
<div>= Using DLVisualize and GULP to Explore the Free Energy and Thermal Expansion of MgO =<br />
<br />
=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''[[File:Js6313Lattice primative.png|thumb|The primative unit cell of MgO]]<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. [[File:Js6313Lattice cube.png|left|thumb|The conventional unit cell for MgO]]When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges. The unit cell is a parallelepiped with <math>\textstyle \alpha = \beta = \gamma = 60^{\circ} </math> and <math>\textstyle a = b = c = 2.9783 \AA </math> which gives a conventional lattice constant of <math> 2.10597 \AA </math>.<br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|Dispersion curve taken with a grid size of 50]]<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly set the temperature of 300 K and the pressure to 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brouillin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation.</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553368Rep:Mod:js6313 mgo2016-03-11T08:31:20Z<p>Js6313: </p>
<hr />
<div>= Using DLVisualize and GULP to Explore the Free Energy and Thermal Expansion of MgO =<br />
<br />
=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''[[File:Js6313Lattice primative.png|thumb|The primative unit cell of MgO]]<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. [[File:Js6313Lattice cube.png|left|thumb|The conventional unit cell for MgO]]When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges. The unit cell is a parallelepiped with <math>\textstyle \alpha = \beta = \gamma = 60^{\circ} </math> and <math>\textstyle a = b = c = 2.9783 \AA which gives a conventional lattice constant of 2.10597 \AA </math>.<br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|Dispersion curve taken with a grid size of 50]]<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly set the temperature of 300 K and the pressure to 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brouillin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation.</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553363Rep:Mod:js6313 mgo2016-03-11T08:24:24Z<p>Js6313: </p>
<hr />
<div>= Using DLVisualize and GULP to Explore the Free Energy and Thermal Expansion of MgO =<br />
<br />
=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''[[File:Js6313Lattice primative.png|thumb|The primative unit cell of MgO]]<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. [[File:Js6313Lattice cube.png|left|thumb]]When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges. The unit cell is a parallelepiped with <math>\textstyle \alpha = \beta = \gamma = 2.10597 \AA </math> and <math>\textstyle a = b = c = 60^{\circ} </math><br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|Dispersion curve taken with a grid size of 50]]<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly use a temperature of 300 K and a pressure of 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brouillin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation.</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553360Rep:Mod:js6313 mgo2016-03-11T08:23:24Z<p>Js6313: </p>
<hr />
<div>= Using DLVisualize and GULP to Explore the Free Energy and Thermal Expansion of MgO =<br />
<br />
=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''[[File:Js6313Lattice primative.png|thumb|The primative unit cell of MgO]]<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. <br />
<br />
When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
[[File:Js6313Lattice cube.png|left|thumb]]<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges. The unit cell is a parallelepiped with <math>\textstyle \alpha = \beta = \gamma = 2.10597 \AA </math> and <math>\textstyle a = b = c = 60^{\circ} </math><br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|Dispersion curve taken with a grid size of 50]]<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly use a temperature of 300 K and a pressure of 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brouillin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation.</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553359Rep:Mod:js6313 mgo2016-03-11T08:22:53Z<p>Js6313: </p>
<hr />
<div>[[File:Js6313Lattice primative.png|thumb|The primative unit cell of MgO]]<br />
<br />
= Using DLVisualize and GULP to Explore the Free Energy and Thermal Expansion of MgO =<br />
<br />
=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. <br />
<br />
When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
[[File:Js6313Lattice cube.png|left|thumb]]<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges. The unit cell is a parallelepiped with <math>\textstyle \alpha = \beta = \gamma = 2.10597 \AA </math> and <math>\textstyle a = b = c = 60^{\circ} </math><br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|Dispersion curve taken with a grid size of 50]]<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly use a temperature of 300 K and a pressure of 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brouillin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation.</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553358Rep:Mod:js6313 mgo2016-03-11T08:22:26Z<p>Js6313: /* Using DLVisualize and GULP to Explore the Free Energy and Thermal Expansion of MgO */</p>
<hr />
<div>[[File:Js6313Lattice primative.png|thumb|The primative unit cell of MgO]]<br />
<br />
= Using DLVisualize and GULP to Explore the Free Energy and Thermal Expansion of MgO =<br />
<br />
=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. <br />
<br />
When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
[[File:Js6313Lattice cube.png|left|thumb]]<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges. The unit cell is a parallelepiped with <math>\textstyle \alpha = \beta = \gamma = 2.10597 \AA </math> and <math>\textstyle a = b = c = 60^{\circ} </math><br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|Dispersion curve taken with a grid size of 50][File:Js6313Lattice cube.png|left|thumb]]<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly use a temperature of 300 K and a pressure of 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brouillin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation.</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=File:Js6313Lattice_primative.png&diff=553353File:Js6313Lattice primative.png2016-03-11T08:20:00Z<p>Js6313: </p>
<hr />
<div></div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=File:Js6313Lattice_cube.png&diff=553351File:Js6313Lattice cube.png2016-03-11T08:18:31Z<p>Js6313: </p>
<hr />
<div></div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553345Rep:Mod:js6313 mgo2016-03-11T08:11:31Z<p>Js6313: </p>
<hr />
<div>= Using DLVisualize and GULP to Explore the Free Energy and Thermal Expansion of MgO =<br />
<br />
=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. <br />
<br />
When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges. The unit cell is a parallelepiped with <math>\textstyle \alpha = \beta = \gamma = 2.10597 \AA </math> and <math>\textstyle a = b = c = 60^{\circ} </math><br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|Dispersion curve taken with a grid size of 50]]<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly use a temperature of 300 K and a pressure of 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brouillin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation.</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553344Rep:Mod:js6313 mgo2016-03-11T08:08:29Z<p>Js6313: </p>
<hr />
<div>= Using DLVisualize and GULP to Explore the Free Energy and Thermal Expansion of MgO =<br />
<br />
=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. <br />
<br />
When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges. The unit cell is a parallelepiped with <math>\textstyle \alpha = \beta = \gamma = 2.1060 \AA </math> and <math>\textstyle a = b = c = 60^{\circ} </math><br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|Dispersion curve taken with a grid size of 50]]<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly use a temperature of 300 K and a pressure of 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brouillin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation.</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553342Rep:Mod:js6313 mgo2016-03-11T08:07:02Z<p>Js6313: /* Introduction */</p>
<hr />
<div>= Using DLVisualize and GULP to Explore the Free Energy and Thermal Expansion of MgO =<br />
<br />
=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. <br />
<br />
When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges. The unit cell is a parallelepiped with <math>\textstyle \alpha = \beta = \gamma = 2.1060 \AA </math> and <math>\textstyle a = b = c = 60 \degree </math><br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|Dispersion curve taken with a grid size of 50]]<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly use a temperature of 300 K and a pressure of 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brouillin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation.</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553341Rep:Mod:js6313 mgo2016-03-11T08:06:46Z<p>Js6313: </p>
<hr />
<div>= Using DLVisualize and GULP to Explore the Free Energy and Thermal Expansion of MgO =<br />
<br />
=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. <br />
<br />
When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges. The unit cell is a parallelepiped with <math>\textstyle \alpha = \beta = \gamma = 2.1060 \AA </math> and <math>\textstyle a = b = c = 60 \degrees </math><br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|Dispersion curve taken with a grid size of 50]]<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly use a temperature of 300 K and a pressure of 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brouillin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation.</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553339Rep:Mod:js6313 mgo2016-03-11T08:06:12Z<p>Js6313: </p>
<hr />
<div>= Using DLVisualize and GULP to Explore the Free Energy and Thermal Expansion of MgO =<br />
<br />
=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. <br />
<br />
When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges. The unit cell is a parallelepiped with <math>\textstyle \alpha = \beta = \gamma = 2.1060 \angstrom</math> and <math>\textstyle a = b = c = 60 \degree</math><br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|Dispersion curve taken with a grid size of 50]]<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly use a temperature of 300 K and a pressure of 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brouillin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation.</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553333Rep:Mod:js6313 mgo2016-03-11T08:01:52Z<p>Js6313: /* Using DLVisualize and GULP to Explore the Free Energy and Thermal Expansion of MgO */</p>
<hr />
<div>= Using DLVisualize and GULP to Explore the Free Energy and Thermal Expansion of MgO =<br />
<br />
=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. <br />
<br />
When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>, which relates how much the crystal expands when subjected to a change in temperature <math>\textstyle \delta V</math>:<br />
<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C. Its high melting point can be attributed to the high lattice enthalpy that results from small ions with high charges.<br />
<br />
In general, both the QMA and MD models are in good agreement with experiment, although the thermal expansion of MgO at very low temperatures deviates from literature when using the former.<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|Dispersion curve taken with a grid size of 50]]<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly use a temperature of 300 K and a pressure of 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brouillin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation.</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553330Rep:Mod:js6313 mgo2016-03-11T07:55:42Z<p>Js6313: </p>
<hr />
<div>= Using DLVisualize and GULP to Explore the Free Energy and Thermal Expansion of MgO =<br />
<br />
=== Introduction ===<br />
'''NEED RESULTS HERE<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<'''<br />
<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. <br />
<br />
When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>:<br />
<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
Magnesium oxide is a basic cubic lattice at temperatures below 2800 °C<br />
<br />
=== Computing the Free Energy - The Quasi-Harmonic Approximation ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|Dispersion curve taken with a grid size of 50]]<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly use a temperature of 300 K and a pressure of 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brouillin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation.</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553315Rep:Mod:js6313 mgo2016-03-11T07:47:06Z<p>Js6313: </p>
<hr />
<div>= Using DLVisualize and GULP to Explore the Free Energy and Thermal Expansion of MgO =<br />
<br />
=== Introduction ===<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. <br />
<br />
When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>:<br />
<br />
<math>\alpha_V = \frac{1}{V} \left (\frac{\delta V}{\delta T} \right )_p </math><br />
<br />
=== Computing the Free Energy - The Harmonic Approximation ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|Dispersion curve taken with a grid size of 50]]<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly use a temperature of 300 K and a pressure of 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brouillin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation.</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553314Rep:Mod:js6313 mgo2016-03-11T07:46:41Z<p>Js6313: </p>
<hr />
<div>= Using DLVisualize and GULP to Explore the Free Energy and Thermal Expansion of MgO =<br />
<br />
=== Introduction ===<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. <br />
<br />
When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>:<br />
<br />
<math>\alpha_V = \frac{1}{V} \left (frac{\delta V}{\delta T} \right )_p </math><br />
<br />
=== Computing the Free Energy - The Harmonic Approximation ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|Dispersion curve taken with a grid size of 50]]<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly use a temperature of 300 K and a pressure of 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brouillin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation.</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553312Rep:Mod:js6313 mgo2016-03-11T07:45:57Z<p>Js6313: /* Using DLVisualize and GULP to Explore the Free Energy and Thermal Expansion of MgO */</p>
<hr />
<div>= Using DLVisualize and GULP to Explore the Free Energy and Thermal Expansion of MgO =<br />
<br />
=== Introduction ===<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. <br />
<br />
When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>:<br />
<br />
<math>\alpha_V = \frac{1}{V} \left (frac{\delta V}{\delta T} )\right_p </math><br />
<br />
=== Computing the Free Energy - The Harmonic Approximation ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|Dispersion curve taken with a grid size of 50]]<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly use a temperature of 300 K and a pressure of 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brouillin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation.</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=553311Rep:Mod:js6313 mgo2016-03-11T07:45:37Z<p>Js6313: </p>
<hr />
<div>= Using DLVisualize and GULP to Explore the Free Energy and Thermal Expansion of MgO =<br />
<br />
=== Introduction ===<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. <br />
<br />
When any crystal is heated, it expands. The extent to which this happens is encapsulated by the thermal expansion coefficient, <math>\textstyle \alpha_V</math>:<br />
<br />
<math>\alpha_V = \frac{1}{V} \left (frac{\delta V}{\delta T} )\right_p<br />
<br />
=== Computing the Free Energy - The Harmonic Approximation ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|Dispersion curve taken with a grid size of 50]]<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly use a temperature of 300 K and a pressure of 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brouillin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation.</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=Rep:Mod:js6313_mgo&diff=552849Rep:Mod:js6313 mgo2016-03-10T19:20:39Z<p>Js6313: </p>
<hr />
<div>= Using DLVisualize and GULP to Explore the Free Energy and Thermal Expansion of MgO =<br />
<br />
=== Introduction ===<br />
At any temperature above 0K, atoms within a crystal vibrate; the modes of these vibrations are called phonons. These can be modelled in a number of ways, including the quasi-harmonic approximation (QMA), where each bond is equivalent to a harmonic oscillator, and equilibrium bond lengths are temperature dependent. A more accurate model hostilities molecular dynamics (MD), which is based upon allowing the dynamic development of forces in agreement with Newton's second law. In order to ascertain a completely representative view of all physical properties of a lattice, the sum over all possible values of k would have to be taken - this would take infinitely long, so a compromise is made between accuracy and computation cost. For the QMA, this was found to be a grid size of 20. <br />
<br />
=== Computing the Free Energy - The Harmonic Approxomation ===<br />
[[File:Js6313 001 phonons 50npoints.png|thumb|Dispersion curve taken with a grid size of 50]]<br />
In order to ascertain the minimum shrinking values required to get an accurate (converged) value for the free energy, jobs were using values of 1, 2, 3, 4, 5, 10 and 20. All calculations here explicitly use a temperature of 300 K and a pressure of 0 GPa. Note how the computation time increases exponentially with shrinking value:<br />
{| class="wikitable"<br />
!Grid size (n)<br />
!Number of Brouillin zone points sampled (<math>\textstyle \frac{n^3}{2}</math>for n > 5)<br />
!Helmholtz free energy (eV)<br />
!Difference between FE at grid size = n and FE at grid size = 50 (meV)<br />
!Calculation time (ms)<br />
|-<br />
|1<br />
|1<br />
|<nowiki>-40.930301</nowiki><br />
|3.818<br />
|2.0<br />
|-<br />
|2<br />
|4<br />
|<nowiki>-40.926609</nowiki><br />
|0.126<br />
|2.8<br />
|-<br />
|3<br />
|18<br />
|<nowiki>-40.926432</nowiki><br />
|0.051<br />
|6.3<br />
|-<br />
|4<br />
|32<br />
|<nowiki>-40.926450</nowiki><br />
|0.033<br />
|9.1<br />
|-<br />
|5<br />
|75<br />
|<nowiki>-40.926463</nowiki><br />
|0.020<br />
|20.9<br />
|-<br />
|10<br />
|500<br />
|<nowiki>-40.926480</nowiki><br />
|0.003<br />
|136.7<br />
|-<br />
|20<br />
|4000<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|1479.8<br />
|-<br />
|50<br />
|62500<br />
|<nowiki>-40.926483</nowiki><br />
|0<br />
|107751.6<br />
|}<br />
<br />
From this, it was decided that a grid size of 20 was a good compromise between calculation time and accuracy; this was used for subsequent calculations using the quasi-harmonic approximation.</div>Js6313https://chemwiki.ch.ic.ac.uk/index.php?title=File:Js6313_001_phonons_50npoints.png&diff=552847File:Js6313 001 phonons 50npoints.png2016-03-10T19:18:25Z<p>Js6313: </p>
<hr />
<div></div>Js6313