
Full configuration interaction quantum Monte Carlo and coupled cluster Monte Carlo: a framework for stochastic quantum chemistry James Spencer1,2 1Thomas Young Centre, Dept. of Materials, Imperial College London 2Condensed Maer Theory, Dept. of Physics, Imperial College London August 4, 2013 Slater determinant space Expand wavefunction in a basis: Ψ(τ) = X f(i) Di (1) | i | i i Slater determinants are a good choice for Fermionic systems: Di = Di1,i2,...,iN | i | i φi1 (1) φi1 (2) φi1 (N) ··· 1 φi2 (1) φi2 (2) φi2 (N) = . ··· . (2) √ . .. N! . φiN (1) φiN (2) φiN (N) ··· (Sadly the Hilbert space grows factorially...) Hamiltonian projection For a given wavefunction ansatz, f(i) satisfy: (Hˆ E) Ψ = 0 (3) − | i Solve coupled equations: Di (Hˆ E) Ψ = 0 (4) h | − | i or Di 1 δτ(Hˆ E) Ψ = Di Ψ (5) h | − − | i h | i Stochastic algorithm energy contribution: for each occupied φT Hˆ Ψ(τ) site i h | | i φT Ψ(τ) h | i for each psip in Sample Di Ψ(τ) wavefunction h | i annihilate parent and child psips select a random site, j spawn new psip on j: Dj Hˆ Di h | | i Sample p(j i) Hamiltonian | sign( Di Ψ(τ) Dj Hˆ Di ) − h | ih | | i FCIQMC Use linear expansion of determinants (G.H. Booth, A.J.W. Thom and A. Alavi, JCP 131 054106 (2009)): Ψ(τ) = X ci(τ) Di (6) | i | i i Particularly eicient for (some) quantum systems. ◮ Applications to atoms, molecules and solids. ◮ Allowed access to FCI energies for Hilbert spaces orders of magnitude greater than conventional diagonalisation. ◮ Recent developments: ◮ Semi-stochastic: Petruzielo et al, PRL 109, 230201 (2012). ◮ Approximations: Roggero et al, arXiv:1304.1549. ◮ Properties via Hellmann–Feynman: JSS, WMC Foulkes (in preparation). ◮ Density Matrix antum Monte Carlo: N. Blunt, T. Rogers, JSS, WMC Foulkes, arXiv:1303.5007. Example: 2D Hubbard model 18 sites, half filling, U/t = 4. 8 10 2.5 (a) / 2.0 1.5 1.0 0.5 No. of psips −14 (b) S(τ) −15 E(τ) − /t 16 EFCI −17 −18 Energy −19 −20 0 5 10 15 20 Uτ Negative sign problem Without annihilation, both positive and negative psips can exist on the same determinant In-phase combination evolves according to H . Out-of-phase −| | combination evolves according to H. 0.00 −0.05 −0.10 −0.15 E(τ) − 0.20 λ(τ) −0.25 E0 λ0 0 50 100 150 200 250 τ/a.u. JSS, N Blunt, WMC Foulkes, JCP 136, 054110 (2012) Population dynamics: 1 component analogue Analytic solution. Captures− key features of FCIQMC simulations: Vmax/κn0 = 25.0; Vmax/Tmax = 30.0 Vmax/κn0 = 12.5; Vmax/Tmax = 15.0 Vmax/κn0 = 12.5; Vmax/Tmax = 30.0 100 80 0 60 /n ) τ ( p 40 20 0 0 20 40 60 80 100 Vmaxτ JSS, N Blunt, WMC Foulkes, JCP 136, 054110 (2012) But... Orders of magnitude improvement equates to only a handful more electrons. Coupled Cluster Monte Carlo Sample coupled cluster wavefunction instead of CI wavefunction: Tˆ Ψ(τ) = e DHF (7) | i | i ˆ a a ab ab T = X ti aˆi + X tij aˆij + (8) ··· ia ijab Highly parallel. Keeps advantanges of traditional CC: fast convergence with truncation and size consistency. Alex Thom, PRL 105, 263004, (2010). Stochastic quantum chemistry Many other quantum chemistry methods exist. Specific methods can be highly accurate for specific systems. Can we exploit the infrastructure required for FCIQMC and CCMC to sample these theories as well? Stochastic algorithm energy contribution: for each occupied φT Hˆ Ψ(τ) site i h | | i φT Ψ(τ) h | i for each psip in Sample Di Ψ(τ) wavefunction h | i annihilate parent and child psips select a random site, j spawn new psip on j: Dj Hˆ Di h | | i Sample p(j i) Hamiltonian | sign( Di Ψ(τ) Dj Hˆ Di ) − h | ih | | i Restricted active space configuration interaction Restrict space in which excitations are allowed to regions of (chemical) interest. Reduces the size of the Hilbert space without (hopefully) removing the important physics of the system of interest. Zimmerman et al, JCP 137 164110 (2012) RAS CI H O; 6-31G** basis set; CISDTQ; RAS=(4,3,8); RAS III 2. 2 ≤ −76.00 S(τ) −76.05 Ep (τ) RAS −76.10 −76.15 −76.20 −76.25 Energy / hartrees Energy −76.30 −76.35 −76.40 0 20 40 60 80 100 τ Spin-flip configuration interaction Truncate the configuration interaction expansion based upon a high spin reference determinant: Useful for studying radicals, high-spin complexes... “Size consistent” if unpaired electrons in reference are localized. Anna Krylov, CPL 338 375 (2001); CPL 350 522 (2001) SF CI Be–Ne; 100A; 6-31G basis set; Ms = 1 spin-flip reference. −1.43e2 0.00 S(τ) Ep (τ) SFCI −0.05 CISD −0.10 Energy / hartrees Energy −0.15 −0.20 0 200 400 600 800 1000 τ Work in progress: random phase approximation Direct RPA is exactly equivalent to direct ring coupled cluster doubles1. Tˆ2 Ansatz: Ψ = e DHF . | i | i Include only ring diagrams in action of Hamiltonian—must modify the generation of random excitations to exclude ladder and mosiac diagrams. 1Scuseria, Henderson, Sorenson, JCP 129 231101 (2008) Conclusions FCIQMC has allowed us to evaluate exact energies for previously impossible systems. Conclusions FCIQMC has allowed us to evaluate exact energies for previously impossible systems. But much, much more: Conclusions FCIQMC has allowed us to evaluate exact energies for previously impossible systems. But much, much more: ◮ General framework for stochastic quantum chemistry methods. Conclusions FCIQMC has allowed us to evaluate exact energies for previously impossible systems. But much, much more: ◮ General framework for stochastic quantum chemistry methods. ◮ Simple to develop and test new theories. Conclusions FCIQMC has allowed us to evaluate exact energies for previously impossible systems. But much, much more: ◮ General framework for stochastic quantum chemistry methods. ◮ Simple to develop and test new theories. ◮ “Baeries included”: ◮ parallelization ◮ lower memory demands ◮ faster ◮ initiator approximation ◮ ... Acknowledgments Alex Thom and Mahew Foulkes. Alavi group (Ali Alavi, George Booth, James Shepherd, Simon Smart, Catherine Overy, Nick Blunt). HANDE contributors: Alex Thom, Will Vigor, Nick Blunt, Tom Rogers, Will Handley. Imperial College High Performance Computing. PSI4 (hp://www.psicode.org)..
Details
-
File Typepdf
-
Upload Time-
-
Content LanguagesEnglish
-
Upload UserAnonymous/Not logged-in
-
File Pages23 Page
-
File Size-