MP2 ground state calculation - Tutorial: Difference between revisions

Revision as of 10:24, 21 November 2018

Overview > MP2 > LTMP2 > stochastic LTMP2 > High energy contributions using stochastic LTMP2 > List of tutorials

(UNDER CONSTRUCTION)

This tutorial introduces how to calculate the ground state energy using second order Møller-Plesset perturbation theory (MP2) with VASP. Currently there are three implementations available:

MP2^[1]: this is the first implementation in VASP and is very efficient for small systems, i.e. systems with less than ~32 valence electrons per unit cell or unit cells smaller than ~150 Å³. The system size scaling of this algorithm is N⁵.
LTMP2^[2]: the Laplace transformed MP2 algorithm has a lower scaling (N⁴) than the previous MP2 algorithm and is therefore efficient for large systems with unit cells larger than ~150 Å³ and more than ~32 valence electrons.
stochastic LTMP2^[3]: faster calculations at the price of statistical noise can be achieved with the stochastic MP2 algorithm. It is an optimal choice for very large systems where only relative errors per valence electron (say 1 meV per valence electron) are relevant. Keeping the absolute error fixed, the algorithm exhibits a cubic scaling with the system size, N³, whereas for a fixed relative error, a linear scaling, N¹, can be achieved. Note that there is no k-point sampling and no spin polarization implemented for this algorithm.

Both LTMP2 as well as stochastic LTMP2 are high performance algorithms that can parallelize the MP2 calculation over thousands of CPUs.

At first, one should select the best algorithm according to the considered system size. In the following, a step by step instruction for each algorithm is presented.

Preparation: the Hartree-Fock ground state

In order to calculate the Hartree-Fock ground state, use the following INCAR file

ISMEAR = 0 ; SIGMA = 0.05 
ALGO = A
LHFCALC = .TRUE. ; AEXX = 1.0
EDIFF = 1E-6
ENCUT = # 10-20% larger than ENMAX in the POTCAR file
LORBITALREAL = .TRUE. # only necessary for LTMP2 and stochastic LTMP2

Keep the OUTCAR file to read-out the Hartree-Fock ground state energy later.

Calculating the unoccupied Hartree-Fock orbitals

We also need the unoccupied/virtual Hartree-Fock orbitals to perform MP2 calculations. The number of necessary orbitals should be equal to the number of plane-waves, that can be found via

nplw=`awk '/number of plane-waves:/ {print $5} ' < OUTCAR_OF_HARTREE_FOCK

For the Gamma-only version of VASP, twice the number of plane-waves have to be used.

Set the INCAR file to

ISMEAR = 0 ; SIGMA = 0.05
ALGO = Exact
LHFCALC = .TRUE. ; AEXX = 1.0
NELM = 1
NBANDS = # number of plane-waves (favorably a multiple of the used mpi-ranks)
ENCUT = # same value as in the Hartree-Fock step
LORBITALREAL = .TRUE. # only necessary for LTMP2 and stochastic LTMP2

Make sure that VASP reads the WAVECAR file from the previous Hartree-Fock step.

Actual MP2 calculations

Depending on your choice, please switch to the corresponding page.

References

[marsman-1] M. Marsman, A. Grüneis, J. Paier, and G. Kresse, J. Chem. Phys. 130, 184103 (2009).

[schaefer2017-2] T. Schäfer, B. Ramberger, and G. Kresse, J. Chem. Phys. 146, 104101 (2017).

[schaefer2018-3] T. Schäfer, B. Ramberger, and G. Kresse, J. Chem. Phys. 148, 064103 (2018).

[1]

[2]

[3]

@@ Line 5: / Line 5: @@
 This tutorial introduces how to calculate the ground state energy using second order Møller-Plesset perturbation theory (MP2) with VASP. Currently there are three implementations available:
-* '''MP2'''<ref name="marsman"/>: this was the first implementation in VASP and is very efficient for small systems, i.e. systems with less than ~32 valence electrons per unit cell or unit cells smaller than ~150 Å&#179;. The system size scaling of this algorithm is N&#8309;.
+* '''MP2'''<ref name="marsman"/>: this is the first implementation in VASP and is very efficient for small systems, i.e. systems with less than ~32 valence electrons per unit cell or unit cells smaller than ~150 Å&#179;. The system size scaling of this algorithm is N&#8309;.
 * '''LTMP2'''<ref name="schaefer2017"/>: the Laplace transformed MP2 algorithm has a lower scaling (N&#8308;) than the previous MP2 algorithm and is therefore efficient for large systems with unit cells larger than ~150 Å&#179; and more than ~32 valence electrons.
 * '''stochastic LTMP2'''<ref name="schaefer2018"/>: faster calculations at the price of statistical noise can be achieved with the stochastic MP2 algorithm. It is an optimal choice for very large systems where only relative errors per valence electron (say 1 meV per valence electron) are relevant. Keeping the absolute error fixed, the algorithm exhibits a cubic scaling with the system size, N&#179;, whereas for a fixed relative error, a linear scaling, N&#185;, can be achieved. Note that there is no k-point sampling and no spin polarization implemented for this algorithm.