Quantum Monte Carlo for excited state calculations
Claudia Filippi
MESA+ Institute for Nanotechnology, Universiteit Twente, The Netherlands
Winter School in Theoretical Chemistry, Helsinki, Finland, Dec 2013 Outline
I Lecture 1: Continuum quantum Monte Carlo methods
− Variational Monte Carlo − Jastrow-Slater wave functions − Diffusion Monte Carlo
I Lecture 2: Quantum Monte Carlo and excited states
− Benchmark and validation: A headache − FCI with random walks: A remarkable development Simple examples: Ethene and butadiene − Assessment of conventional QMC From the gas phase to the protein environment − Prospects and conclusions e.g. Geophysics: Bulk Fe with 96 atoms/cell (Alf´e2009)
Electronic structure methods Scaling − Density functional theory methods N3
− Traditional post Hartree-Fock methods > N6
− Quantum Monte Carlo techniques N4
Stochastic solution of Schr¨odingerequation Accurate calculations for medium-large systems → Molecules of typically 10-50 1st/2nd-row atoms → Relatively little experience with transition metals → Solids (Si, C ... Fe, MnO, FeO) Electronic structure methods Scaling − Density functional theory methods N3
− Traditional post Hartree-Fock methods > N6
− Quantum Monte Carlo techniques N4
Stochastic solution of Schr¨odingerequation Accurate calculations for medium-large systems → Molecules of typically 10-50 1st/2nd-row atoms → Relatively little experience with transition metals → Solids (Si, C ... Fe, MnO, FeO) e.g. Geophysics: Bulk Fe with 96 atoms/cell (Alf´e2009) Quantum Monte Carlo as a wave function based method
Ψ(r1,..., rN ) → CASSCF, MRCI, CASPT2, CC ...
Can we construct an accurate but more compact Ψ?
Explicit dependence on the inter-electronic distances rij Use of non-orthogonal orbitals ...
→ Quantum Monte Carlo as alternative route A different way of writing the expectation values
Consider the expectation value of the Hamiltonian on Ψ
hΨ|H|Ψi R dR Ψ∗(R)HΨ(R) E = = ≥ E V hΨ|Ψi R dR Ψ∗(R)Ψ(R) 0
Z HΨ(R) |Ψ(R)|2 = dR Ψ(R) R dR|Ψ(R)|2
Z ←− = dR EL(R) ρ(R)= hEL(R)iρ
HΨ(R) ρ(R) is a probability density and E (R) = the local energy L Ψ(R) Variational Monte Carlo: a random walk of the electrons
Use Monte Carlo integration to compute expectation values
. Sample R from ρ(R) using Metropolis algorithm HΨ(R) . Average local energy E (R) = to obtain E as L Ψ(R) V
M 1 X E = hE (R)i ≈ E (R ) V L ρ M L i i=1 R
Random walk in 3N dimensions, R = (r1,..., rN )
Just a trick to evaluate integrals in many dimensions Is it really “just” a trick?
Si21H22
Number of electrons 4 × 21 + 22 = 106 Number of dimensions 3 × 106 = 318
Integral on a grid with 10 points/dimension → 10318 points!
MC is a powerful trick ⇒ Freedom in form of the wave function Ψ Are there any conditions on many-body Ψ to be used in VMC?
Within VMC, we can use any “computable” wave function if
. Continuous, normalizable, proper symmetry
. Finite variance
hΨ|(H − E )2|Ψi σ2 = V = h(E (R) − E )2i hΨ|Ψi L V ρ
σ since the Monte Carlo error goes as err(EV ) ∼ √ M
Zero variance principle: if Ψ → Ψ0, EL(R) does not fluctuate Expectation values in variational Monte Carlo: The energy
hΨ|H|Ψi Z HΨ(R) |Ψ(R)|2 Z E = = dR = dR E (R) ρ(R) V hΨ|Ψi Ψ(R) R dR|Ψ(R)|2 L Z EV = dR EL(R) ρ(R) Z 2 2 σ = dR(EL(R) − EV ) ρ(R)
Estimate EV and σ from M independent samples as
M 1 X E¯ = E (R ) V M L i i=1 M 1 X σ¯2 = (E (R ) − E¯ )2 M L i V i=1 Typical VMC run
Example: Local energy and average energy of acetone (C3H6O) -34
-35
-36 σ VMC
-37 Energy (Hartree) -38
-39 0 500 1000 1500 2000 MC step
EVMC = hEL(R)iρ = −36.542 ± 0.001 Hartree (40×20000 steps)
2 σVMC = h(EL(R) − EVMC) iρ = 0.90 Hartree Variational Monte Carlo and the generalized Metropolis algorithm
|Ψ(R)|2 How do we sample distribution function ρ(R) = ? R dR|Ψ(R)|2
Aim → Obtain a set of {R1, R2,..., RM } distributed as ρ(R)
Generate a Markov chain
. Start from arbitrary initial state Ri
. Use stochastic transition matrix P(Rf |Ri) X P(Rf |Ri) ≥ 0 P(Rf |Ri) = 1.
Rf R as probability of making transition Ri → Rf
. Evolve the system by repeated application of P How do we sample ρ(R)?
To sample ρ, use P which satisfies stationarity condition :
X P(Rf |Ri) ρ(Ri) = ρ(Rf ) ∀ Rf i
⇒ If we start with ρ, we continue to sample ρ
Stationarity condition + stochastic property of P + ergodicity
⇒ Any initial distribution will evolve to ρ A more stringent condition
In practice, we impose detailed balance condition