Example: Graphene sheet with dispersion correction

Download Graphene_Dispersion.run

A normal GGA would give only negligible interaction between two graphene sheets.

Use the dispersion option in the XC key block.

In the first run we use BP86-D, in the second BLYP-D3 and in the third run BLYP-D3(BJ).

# First run: dispersion default

"$ADFBIN/band" <<eor
TITLE Dispersion energy with two parallel graphene sheets

XC
   gga scf bp86
   dispersion default
End

NumericalQuality Basic

GeometryFile $ADFHOME/examples/band/Graphene_Dispersion/Graphene_double_layer.xyz

BasisDefaults
   BasisType TZP
   Core Large
End

end input
eor

rm RUNKF

# Second run: dispersion Grimme3

"$ADFBIN/band" <<eor
TITLE Grimme3 dispersion energy with two parallel graphene sheets

XC
   gga scf blyp
   dispersion Grimme3
end

Gradients
End

NumericalQuality Basic

GeometryFile $ADFHOME/examples/band/Graphene_Dispersion/Graphene_double_layer.xyz

BasisDefaults
   BasisType TZP
   Core Large
End

end input
eor

rm RUNKF

# Third run: dispersion Grimme3 bjdamp

"$ADFBIN/band" <<eor

XC
   gga scf blyp
   dispersion Grimme3 bjdamp
end

Gradients
End

NumericalQuality Basic

GeometryFile $ADFHOME/examples/band/Graphene_Dispersion/Graphene_double_layer.xyz

BasisDefaults
   BasisType TZP
   Core Large
End

end input
eor