Simulating icosahedral gold clusters
Revision as of 05:29, 1 April 2022 by Miguel Caro (talk | contribs) (Created page with "In this tutorial we are going to use a GAP to optimize the structure of an icosahedral gold cluster and estimate its formation energy. In this tutorial we will make use of sta...")
In this tutorial we are going to use a GAP to optimize the structure of an icosahedral gold cluster and estimate its formation energy. In this tutorial we will make use of static (single-point) calculations (turbogap predict) and geometry optimization.
Prerequisites for this tutorial
- A TurboGAP installation
- An ASE installation
- Numpy
- A utility that can perform simple least-squares regression (I will use gnuplot, which lets me plot the fit at the same time)
Computing the equilibrium lattice parameter of the reference structure
The thermodynamically stable crystal structure of gold is fcc. When we compute the formation energy of a gold nanocluster we will be referring its energy per atom to that of the fcc structure. Therefore, the first step in this tutorial is to find out the lattice parameter of fcc gold by computing the energy vs lattice parameter curve, and finding the minimum.
