Download - Full configuration interaction quantum Monte Carlo and coupled cluster … · 2013-08-14 · Full configuration interaction quantum Monte Carlo and coupled cluster Monte Carlo: a

Transcript
Page 1: Full configuration interaction quantum Monte Carlo and coupled cluster … · 2013-08-14 · Full configuration interaction quantum Monte Carlo and coupled cluster Monte Carlo: a

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 Ma�er Theory, Dept. of Physics, Imperial College London

August 4, 2013

Page 2: Full configuration interaction quantum Monte Carlo and coupled cluster … · 2013-08-14 · Full configuration interaction quantum Monte Carlo and coupled cluster Monte Carlo: a

Slater determinant space

Expand wavefunction in a basis:

|Ψ(τ)〉 =∑

i

f(i) |Di〉 (1)

Slater determinants are a good choice for Fermionic systems:

|Di〉 = |Di1,i2,...,iN 〉

=1√N !

φi1(1) φi1(2) · · · φi1(N)φi2(1) φi2(2) · · · φi2(N)

......

. . ....

φiN (1) φiN (2) · · · φiN (N)

(2)

(Sadly the Hilbert space grows factorially. . . )

Page 3: Full configuration interaction quantum Monte Carlo and coupled cluster … · 2013-08-14 · Full configuration interaction quantum Monte Carlo and coupled cluster Monte Carlo: a

Hamiltonian projection

For a given wavefunction ansatz, f(i) satisfy:

(H − E) |Ψ〉 = 0 (3)

Solve coupled equations:

〈Di|(H − E)|Ψ〉 = 0 (4)

or

〈Di|1− δτ(H − E)|Ψ〉 = 〈Di|Ψ〉 (5)

Page 4: Full configuration interaction quantum Monte Carlo and coupled cluster … · 2013-08-14 · Full configuration interaction quantum Monte Carlo and coupled cluster Monte Carlo: a

Stochastic algorithmfor each occupied

site i

energy contribution:

〈φT |H|Ψ(τ)〉〈φT |Ψ(τ)〉

for each psip in

〈Di|Ψ(τ)〉

select a random

site, j

annihilate parent

and child psips

spawn new psip on j:

〈Dj |H|Di〉p(j|i)

−sign(〈Di|Ψ(τ)〉 〈Dj |H|Di〉)

Sample

wavefunction

Sample

Hamiltonian

Page 5: Full configuration interaction quantum Monte Carlo and coupled cluster … · 2013-08-14 · Full configuration interaction quantum Monte Carlo and coupled cluster Monte Carlo: a

FCIQMC

Use linear expansion of determinants (G.H. Booth, A.J.W. Thom and

A. Alavi, JCP 131 054106 (2009)):

|Ψ(τ)〉 =∑

i

ci(τ) |Di〉 (6)

Particularly e�icient 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.

Page 6: Full configuration interaction quantum Monte Carlo and coupled cluster … · 2013-08-14 · Full configuration interaction quantum Monte Carlo and coupled cluster Monte Carlo: a

Example: 2D Hubbard model

18 sites, half filling, U/t = 4.

0.51.01.52.02.5

No.

of

psi

ps/108

(a)

0 5 10 15 20

−20

−19

−18

−17

−16

−15

−14

En

erg

y/t

(b) S(τ)

E(τ)EFCI

Page 7: Full configuration interaction quantum Monte Carlo and coupled cluster … · 2013-08-14 · Full configuration interaction quantum Monte Carlo and coupled cluster Monte Carlo: a

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-phasecombination evolves according toH .

0 50 100 150 200 250

τ/a.u.

−0.25

−0.20

−0.15

−0.10

−0.05

0.00

E(τ)

λ(τ)E0

λ0

JSS, N Blunt, WMC Foulkes, JCP 136, 054110 (2012)

Page 8: Full configuration interaction quantum Monte Carlo and coupled cluster … · 2013-08-14 · Full configuration interaction quantum Monte Carlo and coupled cluster Monte Carlo: a

Population dynamics: 1−component analogueAnalytic solution. Captures key features of FCIQMC simulations:

0 20 40 60 80 100

Vmaxτ

0

20

40

60

80

100

p(τ)/n0

Vmax/κn0 = 25.0; Vmax/Tmax = 30.0

Vmax/κn0 = 12.5; Vmax/Tmax = 15.0

Vmax/κn0 = 12.5; Vmax/Tmax = 30.0

JSS, N Blunt, WMC Foulkes, JCP 136, 054110 (2012)

Page 9: Full configuration interaction quantum Monte Carlo and coupled cluster … · 2013-08-14 · Full configuration interaction quantum Monte Carlo and coupled cluster Monte Carlo: a

But. . .

Orders of magnitude improvement equates to only a handful more

electrons.

Page 10: Full configuration interaction quantum Monte Carlo and coupled cluster … · 2013-08-14 · Full configuration interaction quantum Monte Carlo and coupled cluster Monte Carlo: a

Coupled Cluster Monte Carlo

Sample coupled cluster wavefunction instead of CI wavefunction:

|Ψ(τ)〉 = eT |DHF〉 (7)

T =∑

ia

tai aai +

ijab

tabij aabij + · · · (8)

Highly parallel.

Keeps advantanges of traditional CC: fast convergence with

truncation and size consistency.

Alex Thom, PRL 105, 263004, (2010).

Page 11: Full configuration interaction quantum Monte Carlo and coupled cluster … · 2013-08-14 · Full configuration interaction quantum Monte Carlo and coupled cluster Monte Carlo: a

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?

Page 12: Full configuration interaction quantum Monte Carlo and coupled cluster … · 2013-08-14 · Full configuration interaction quantum Monte Carlo and coupled cluster Monte Carlo: a

Stochastic algorithmfor each occupied

site i

energy contribution:

〈φT |H|Ψ(τ)〉〈φT |Ψ(τ)〉

for each psip in

〈Di|Ψ(τ)〉

select a random

site, j

annihilate parent

and child psips

spawn new psip on j:

〈Dj |H|Di〉p(j|i)

−sign(〈Di|Ψ(τ)〉 〈Dj |H|Di〉)

Sample

wavefunction

Sample

Hamiltonian

Page 13: Full configuration interaction quantum Monte Carlo and coupled cluster … · 2013-08-14 · Full configuration interaction quantum Monte Carlo and coupled cluster Monte Carlo: a

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)

Page 14: Full configuration interaction quantum Monte Carlo and coupled cluster … · 2013-08-14 · Full configuration interaction quantum Monte Carlo and coupled cluster Monte Carlo: a

RAS CIH2O; 6-31G** basis set; CISDTQ; RAS=(4,3,8); RAS III ≤ 2.

0 20 40 60 80 100τ

−76.40

−76.35

−76.30

−76.25

−76.20

−76.15

−76.10

−76.05

−76.00

Energ

y /

hart

rees

S(τ)

Ep (τ)

RAS

Page 15: Full configuration interaction quantum Monte Carlo and coupled cluster … · 2013-08-14 · Full configuration interaction quantum Monte Carlo and coupled cluster Monte Carlo: a

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)

Page 16: Full configuration interaction quantum Monte Carlo and coupled cluster … · 2013-08-14 · Full configuration interaction quantum Monte Carlo and coupled cluster Monte Carlo: a

SF CIBe–Ne; 100A; 6-31G basis set; Ms = 1 spin-flip reference.

0 200 400 600 800 1000τ

−0.20

−0.15

−0.10

−0.05

0.00

Energ

y /

hart

rees

−1.43e2

S(τ)

Ep (τ)

SFCICISD

Page 17: Full configuration interaction quantum Monte Carlo and coupled cluster … · 2013-08-14 · Full configuration interaction quantum Monte Carlo and coupled cluster Monte Carlo: a

Work in progress: random phase approximation

Direct RPA is exactly equivalent to direct ring coupled cluster

doubles1.

Ansatz: |Ψ〉 = eT2 |DHF〉.

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)

Page 18: Full configuration interaction quantum Monte Carlo and coupled cluster … · 2013-08-14 · Full configuration interaction quantum Monte Carlo and coupled cluster Monte Carlo: a

Conclusions

FCIQMC has allowed us to evaluate exact energies for previously

impossible systems.

Page 19: Full configuration interaction quantum Monte Carlo and coupled cluster … · 2013-08-14 · Full configuration interaction quantum Monte Carlo and coupled cluster Monte Carlo: a

Conclusions

FCIQMC has allowed us to evaluate exact energies for previously

impossible systems.

But much, much more:

Page 20: Full configuration interaction quantum Monte Carlo and coupled cluster … · 2013-08-14 · Full configuration interaction quantum Monte Carlo and coupled cluster Monte Carlo: a

Conclusions

FCIQMC has allowed us to evaluate exact energies for previously

impossible systems.

But much, much more:

◮ General framework for stochastic quantum chemistry methods.

Page 21: Full configuration interaction quantum Monte Carlo and coupled cluster … · 2013-08-14 · Full configuration interaction quantum Monte Carlo and coupled cluster Monte Carlo: a

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.

Page 22: Full configuration interaction quantum Monte Carlo and coupled cluster … · 2013-08-14 · Full configuration interaction quantum Monte Carlo and coupled cluster Monte Carlo: a

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.

◮ “Ba�eries included”:◮ parallelization◮ lower memory demands◮ faster◮ initiator approximation◮ . . .

Page 23: Full configuration interaction quantum Monte Carlo and coupled cluster … · 2013-08-14 · Full configuration interaction quantum Monte Carlo and coupled cluster Monte Carlo: a

Acknowledgments

Alex Thom and Ma�hew 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 (h�p://www.psicode.org).