Graphitization simulation with van der Waals corrections
WARNING: TUTORIAL UNDER CONSTRUCTION!!!!
In this tutorial we are going to carry out simulations of carbon graphitization similar to those in the simple molecular dynamics tutorial, but we will also add a barostat so that the simulation box size changes to mimick the effect of coupling to an external pressure. We will also look at how to use van der Waals (vdW) corrections on top of a GAP force field. Note that vdW corrections are generally not available for all potentials. Those potentials without "proper" vdW corrections (but again, not all) may include fixed pair-wise long-range interactions that allow you to incorporate dispersion effects in a limited way (but we will not cover those in this tutorial).
It is strongly recommended that you do the simple molecular dynamics tutorial before this one.
Prerequisites for this tutorial
- A TurboGAP installation
Optional
- gnuplot (for plotting)
- An ASE installation
- VMD (for visualization, ASE can handle visualization but it's slow)
Initial configuration
We will start our simulation from a melted carbon sample with 216 atoms, generated as in the simple molecular dynamics tutorial. For convenience, you can download the file from here.
