• Volume/Page
  • Keyword
  • DOI
  • Citation
  • Advanced
   
 
 
 

You Tube Flickr Twitter UniPHY Group iResearch App Facebook

J. Chem. Phys. 132, 154113 (2010); http://dx.doi.org/10.1063/1.3380831 (9 pages)

Size-consistent variational approaches to nonlocal pseudopotentials: Standard and lattice regularized diffusion Monte Carlo methods revisited

Michele Casula1, Saverio Moroni2,3, Sandro Sorella2,3, and Claudia Filippi4

1Centre de Physique Théorique, CNRS, École Polytechnique, 91128 Palaiseau Cedex, France
2International School for Advanced Studies (SISSA), Via Beirut 2-4, 34014 Trieste, Italy
3CNR-IOM Democritos National Simulation Center, Via Beirut 2-4, 34014 Trieste, Italy
4Faculty of Science and Technology and MESA+Research Institute, University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands

View MapView Map

(Received 1 February 2010; accepted 13 March 2010; published online 21 April 2010)

We propose improved versions of the standard diffusion Monte Carlo (DMC) and the lattice regularized diffusion Monte Carlo (LRDMC) algorithms. For the DMC method, we refine a scheme recently devised to treat nonlocal pseudopotential in a variational way. We show that such scheme—when applied to large enough systems—maintains its effectiveness only at correspondingly small enough time-steps, and we present two simple upgrades of the method which guarantee the variational property in a size-consistent manner. For the LRDMC method, which is size-consistent and variational by construction, we enhance the computational efficiency by introducing: (i) an improved definition of the effective lattice Hamiltonian which remains size-consistent and entails a small lattice-space error with a known leading term and (ii) a new randomization method for the positions of the lattice knots which requires a single lattice-space.

© 2010 American Institute of Physics

Article Outline

  1. INTRODUCTION
  2. DMC AND NONLOCAL PSEUDOPOTENTIALS
    1. Beyond the LA
  3. SIZE-CONSISTENCY
    1. Size-consistent variational formulations: SVDMC Version 1
    2. Size-consistent variational formulations: SVDMC Version 2
  4. LRDMC AND NONLOCAL PSEUDOPOTENTIALS
    1. Small a2 correction for good trial function
    2. Well defined lattice regularization
  5. PERFORMANCE OF THE PROPOSED METHODS
    1. Time-step error
    2. Relative efficiency
  6. CONCLUSIONS

RELATED DATABASES

To view database links for this article, you need to log in.

KEYWORDS and PACS

PACS

  • 71.15.Dx

    Computational methodology (Brillouin zone sampling, iterative diagonalization, pseudopotential construction)

  • 31.15.xt

    Variational techniques

ARTICLE DATA

PUBLICATION DATA

ISSN

0021-9606 (print)  
1089-7690 (online)

For access to fully linked references, you need to log in.
    W. M. C. Foulkes, L. Mitas, R. J. Needs, and G. Rajagopal, Rev. Mod. Phys. 73, 33 (2001).

    B. L. Hammond, P. J. Reynolds, and W. A. Lester, Jr., J. Chem. Phys. 87, 1130 (1987)JCPSA6000087000002001130000001.

    S. Fahy, X. W. Wang, and S. G. Louie, Phys. Rev. B 42, 3503 (1990).

    L. Mitáš, E. L. Shirley, and D. M. Ceperley, J. Chem. Phys. 95, 3467 (1991)JCPSA6000095000005003467000001.

    M. Casula, C. Filippi, and S. Sorella, Phys. Rev. Lett. 95, 100201 (2005).

    M. Casula, Phys. Rev. B 74, 161102(R) (2006).

    D. F. B. ten Haaf, H. J. M. van Bemmel, J. M. J. van Leeuwen, W. van Saarloos, and D. M. Ceperley, Phys. Rev. B 51, 13039 (1995).

    M. Burkatzki, C. Filippi, and M. Dolg, J. Chem. Phys. 126, 234105 (2007)JCPSA6000126000023234105000001.

    C. Filippi and C. J. Umrigar, J. Chem. Phys. 105, 213 (1996)JCPSA6000105000001000213000001. As Jastrow correlation factor, we use the exponential of the sum of two fifth-order polynomials of the electron-nuclear (e-n) and the electron-electron (e-e) distances, respectively. The Jastrow factor is adapted to deal with pseudoatoms, and the scaling factor kappa is set to 0.8 a.u. for the oxygen systems and 0.5 a.u. for oxirane.

    C. J. Umrigar, J. Toulouse, C. Filippi, S. Sorella, and R. G. Henning, Phys. Rev. Lett. 98, 110201 (2007).

    The short time propagator of DMC version 2 lends itself to an implementation of the Reptation quantum Monte Carlo method [S. Baroni and S. Moroni, Phys. Rev. Lett. 82, 4745 (1999)] with single-particle moves.

    S. Sorella and L. Capriotti, Phys. Rev. B 61, 2599 (2000).

    I. G. Gurtubay and R. J. Needs, J. Chem. Phys. 127, 124306 (2007)JCPSA6000127000012124306000001.


For access to citing articles, you need to log in.


Figures (4)

Access to article objects (figures, tables, multimedia) requires a subscription; log in to view available files.
(Access to supplementary files, where available, is free for this journal.)


Close
Google Calendar
ADVERTISEMENT

close