Rights statement: © 2016 American Physical Society
Accepted author manuscript, 293 KB, PDF document
Research output: Contribution to Journal/Magazine › Journal article › peer-review
Research output: Contribution to Journal/Magazine › Journal article › peer-review
}
TY - JOUR
T1 - Trail-Needs pseudopotentials in quantum Monte Carlo calculations with plane-wave/blip basis sets
AU - Drummond, Neil David
AU - Trail, J. R.
AU - Needs, Richard
N1 - © 2016 American Physical Society
PY - 2016/10/15
Y1 - 2016/10/15
N2 - We report a systematic analysis of the performance of a widely used set of Dirac-Fock pseudopotentials for quantum Monte Carlo (QMC) calculations. We study each atom in the periodic table from hydrogen (Z = 1) to mercury (Z = 80), with the exception of the 4f elements (57 ≤ Z ≤ 70).We demonstrate that ghost states are a potentially serious problem when plane-wave basis sets are used in density functional theory (DFT) orbital-generation calculations, but that this problem can be almost entirely eliminated by choosing the s channel to be local in the DFT calculation; thed channel can then be chosen to be local in subsequent QMC calculations, which generally leadsto more accurate results. We investigate the achievable energy variance per electron with differentlevels of trial wave function and we determine appropriate plane-wave cutoff energies for DFT calculations for each pseudopotential. We demonstrate that the so-called “T-move” scheme in diffusion Monte Carlo is essential for many elements. We investigate the optimal choice of spherical integration rule for pseudopotential projectors in QMC calculations. The information reported here will prove crucial in the planning and execution of QMC projects involving beyond-first-row elements.
AB - We report a systematic analysis of the performance of a widely used set of Dirac-Fock pseudopotentials for quantum Monte Carlo (QMC) calculations. We study each atom in the periodic table from hydrogen (Z = 1) to mercury (Z = 80), with the exception of the 4f elements (57 ≤ Z ≤ 70).We demonstrate that ghost states are a potentially serious problem when plane-wave basis sets are used in density functional theory (DFT) orbital-generation calculations, but that this problem can be almost entirely eliminated by choosing the s channel to be local in the DFT calculation; thed channel can then be chosen to be local in subsequent QMC calculations, which generally leadsto more accurate results. We investigate the achievable energy variance per electron with differentlevels of trial wave function and we determine appropriate plane-wave cutoff energies for DFT calculations for each pseudopotential. We demonstrate that the so-called “T-move” scheme in diffusion Monte Carlo is essential for many elements. We investigate the optimal choice of spherical integration rule for pseudopotential projectors in QMC calculations. The information reported here will prove crucial in the planning and execution of QMC projects involving beyond-first-row elements.
M3 - Journal article
VL - 94
JO - Physical review B
JF - Physical review B
SN - 1098-0121
IS - 16
M1 - 165170
ER -