1 / 30

Order(N) methods in QMC Dario Alfè [d.alfe@ucl.ac.uk]

2007 Summer School on Computational Materials Science Quantum Monte Carlo: From Minerals and Materials to Molecules July 9 –19, 2007 • University of Illinois at Urbana–Champaign http://www.mcc.uiuc.edu/summerschool/2007/qmc/. Order(N) methods in QMC Dario Alfè [d.alfe@ucl.ac.uk].

jewell
Download Presentation

Order(N) methods in QMC Dario Alfè [d.alfe@ucl.ac.uk]

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. 2007 Summer School on Computational Materials ScienceQuantum Monte Carlo: From Minerals and Materials to Molecules July 9 –19, 2007 • University of Illinois at Urbana–Champaign http://www.mcc.uiuc.edu/summerschool/2007/qmc/ Order(N) methods in QMC Dario Alfè [d.alfe@ucl.ac.uk]

  2. Monte Carlo methods • Evaluation of multidimensional integrals: • Quadrature methods, error ~ M-4/d • Monte Carlo methods, error ~ M-0.5

  3. Variational Monte Carlo

  4. Diffusion Monte Carlo Extracting the ground state: substitute  = it Fixed nodes:

  5. Cost of evaluating T • N electrons • N single particle orbitals i • ~ N basis functions for each i if Cost proportional to N3

  6. Cost of evaluating Energy E Variance of E is proportional to N  Cost of total energy proportional to N4 HOWEVER Variance of E/atom is proportional to 1/N  Cost of energy/atom proportional to N2 (relevant to free energies in phase transitions, surface energies, ….)

  7. Cost of evaluating T • N electrons • N single particle orbitals i • ~ N basis functions for each i if Cost proportional to N3

  8. Cost of evaluating T • N electrons • N single particle orbitals i • ~ o(1) basis functions for each i if if fl(r) is localised Cost proportional to N2

  9. Cost of evaluating Energy E Variance of E is proportional to N  Cost of total energy proportional to N3 HOWEVER Variance of E/atom is proportional to 1/N  Cost of energy/atom proportional to N (relevant to free energies in phase transitions, surface energies, ….)

  10. B-splines: Localised functions sitting at the points of a uniform grid E. Hernàndez, M. J. Gillan and C. M. Goringe, Phys. Rev. B, 20 (1997) D. Alfè  and M. J. Gillan, Phys. Rev. B, Rapid Comm., 70, 161101, (2004)

  11. Cost of evaluating T • N electrons • N single particle orbitals i • ~ o(1) basis functions for each i if Cost proportional to N2

  12. Cost of evaluating T • N electrons • ~ o(1) single particle localised orbitals i,if i is localised • ~ o(1) basis functions for each i if Cost proportional to N(Linear scaling)

  13. Cost of evaluating Energy E Variance of E is proportional to N  Cost of total energy proportional to N2 HOWEVER Variance of E/atom is proportional to 1/N  Cost of energy/atom independent on N! (relevant to free energies in phase transitions, surface energies, ….)

  14. B-splines (blips) • Defined on a uniform grid • Localised: f(x)=0 for |x|>2 • Continuous with first and second derivative continuous Grid spacing 1: Grid spacing a:

  15. Blips in three dimensions: Single-particle orbital representation: For each position r there are only 64 blip functions that are non zero, for any size and any grid spacing a. By contrast, in a plane wave representation the number of plane waves is proportional to the size of the system, and also depends on the plane wave cutoff: For example, the number of plane waves in MgO is ~ 3000 per atom! (HF pseudopotentials, 100 Ha PW cutoff)

  16. Approximate equivalence between blips and plane waves One dimension:

  17. Now consider: If nand nwere proportional then n= 1, and blip waves would be identical to plane waves. Let’s see what value of n has for a few values on n. Take N=12, for example Therefore if we replace nwithn we make a small error at low k but maybe a significant error for large values ofk. Let’s do it anyway, then we will come back to the size of the error and how to reduce it.

  18. Representing single particle orbitals in 3 d

  19. Improving the quality of the B-spline representation Reduce a and achieve better quality May not be good enough Blips are systematically improvable

  20. Example: Silicon in the ß-Sn structure, 16 atoms, 1st orbital

  21. Energies in QMC (eV/atom)(Silicon inß-Sn structure, 16 atoms, 15 Ry PW cutoff)

  22. Energies in QMC (eV/atom)(MgO inNaCl structure, 8 atoms, 200 Ry PW cutoff) Blips and plane waves give identical results, blips are 2 orders of magnitude faster on 8 atoms, 3 and 4 order of magnitudes faster on 80 ad 800 atoms respectively (on 800 atoms: Blips 1 day, PW 38 years)

  23. Achieving linear scaling Maximally localised Wannier functions(Marzari Vanderbilt) [Williamson et al., PRL, 87, 246406 (2001)] Unitary transformation: Linear scaling obtained by truncating orbitals to zero outside their localisation region. And:

  24. The VMC and DMC total energies are exactly invariant with respect to arbitrary (non singular) linear combinations of single-electron orbitals! Transformation does not have to be unitary. Arbitrary linear combination: And: Provided

  25. We have set of orbitals extending over region . Want to make orbital maximally localised in sub-region . Maximise “localisation weight” P: Linear scaling obtained by truncating localised orbitals to zero outside their localisation region . Then for given r number of non-zero orbitals is O(1). Express P as: where: Weight P is maximal when cneigenvector of associated with maximum eigenvalue , and

  26. Localisation weight in MgO Convergence of localisation weight Q=1-P in MgO. Squares: present method (D. Alfè  and M. J. Gillan, J. Phys. Cond. Matter 16, L305-L311 (2004)) Diamonds: MLWF (Williamson et al. PRL, 87, 246406 (2001)

  27. Localisation weight in MgO Convergence of localisation weight Q=1-P in MgO. Squares: present method (D. Alfè  and M. J. Gillan, J. Phys. Cond. Matter 16, L305-L311 (2004)) Diamonds: MLWF (Williamson et al. PRL, 87, 246406 (2001)

  28. VMC and DMC energies in MgO Convergence of VMC energy. Squares: Non-orthogonal orbitals. Diamons: MLWF Convergence of DMC energy, non-orthogonal orbitals

  29. VMC and DMC energies in MgO Convergence of VMC energy. Squares: Non-orthogonal orbitals. Diamons: MLWF Convergence of DMC energy, non-orthogonal orbitals Spherical cutoff at 6 a.u on a 7.8 a.u. cubic box (64 atoms): - WF evaluation 4 times faster - Memory occupancy 4 times smaller

More Related