Recent advances in quantum Monte Carlo for quantum chemistry: optimization of wave functions and calculation of observables
|
|
- James Cameron
- 5 years ago
- Views:
Transcription
1 Recent advances in quantum Monte Carlo for quantum chemistry: optimization of wave functions and calculation of observables Julien Toulouse 1, Cyrus J. Umrigar 2, Roland Assaraf 1 1 Laboratoire de Chimie Théorique, Université Pierre et Marie Curie - CNRS, Paris, France. 2 Laboratory of Atomic and Solid State Physics, Cornell University, Ithaca, New York, USA. julien.toulouse@upmc.fr Web page: March 2009
2 1 Optimization of wave functions 2 Calculation of observables
3 1 Optimization of wave functions 2 Calculation of observables
4 Trial wave function Jastrow-Slater wave function Ψ(p) = Ĵ(α) N CSF i=1 c i C i Ĵ(α) = Jastrow factor (with e-e, e-n, e-e-n terms) C i = Configuration state function (CSF) = linear combination of Slater determinants of given symmetry.
5 Trial wave function Jastrow-Slater wave function Ψ(p) = Ĵ(α) N CSF i=1 c i C i Ĵ(α) = Jastrow factor (with e-e, e-n, e-e-n terms) C i = Configuration state function (CSF) = linear combination of Slater determinants of given symmetry. The Slater determinants are made of orbitals expanded on a Slater basis: φ k (r) = N basis µ=1 λ kµ χ µ (r) χ(r) = N(ζ)r n 1 e ζr S l,m (θ, φ)
6 Trial wave function Jastrow-Slater wave function Ψ(p) = Ĵ(α) N CSF i=1 c i C i Ĵ(α) = Jastrow factor (with e-e, e-n, e-e-n terms) C i = Configuration state function (CSF) = linear combination of Slater determinants of given symmetry. The Slater determinants are made of orbitals expanded on a Slater basis: φ k (r) = N basis µ=1 λ kµ χ µ (r) χ(r) = N(ζ)r n 1 e ζr S l,m (θ, φ) Parameters to optimize p = {α, c, λ, ζ}: Jastrow parameters α, CSF coefficients c, orbital coefficients λ and basis exponents ζ
7 Wave function optimization: why and how? Important for both VMC and DMC in order to reduce the systematic error reduce the statistical uncertainty
8 Wave function optimization: why and how? Important for both VMC and DMC in order to reduce the systematic error reduce the statistical uncertainty How to optimize? Until recently: minimization of the variance of the energy OK for the few Jastrow parameters but does not work well for the many CSF and orbital parameters Since recently: minimization of the energy (+ possibly a small fraction of variance) works well for all the parameters the energy is a better criterion
9 Optimization method: principle Expansion of the wave function around p 0 to linear order in p = p p 0 : Ψ [1] (p) = Ψ 0 + j p j Ψ j where Ψ 0 = Ψ(p 0 ) and Ψ j = Ψ(p0 )) p j.
10 Optimization method: principle Expansion of the wave function around p 0 to linear order in p = p p 0 : Ψ [1] (p) = Ψ 0 + j p j Ψ j where Ψ 0 = Ψ(p 0 ) and Ψ j = Ψ(p0 )) p j. Normalization of wave function chosen so that the derivatives Ψ j are orthogonal to Ψ 0.
11 Optimization method: principle Expansion of the wave function around p 0 to linear order in p = p p 0 : Ψ [1] (p) = Ψ 0 + j p j Ψ j where Ψ 0 = Ψ(p 0 ) and Ψ j = Ψ(p0 )) p j. Normalization of wave function chosen so that the derivatives Ψ j are orthogonal to Ψ 0. Minimization of the energy = generalized eigenvalue equation: ( E0 g T )( ) ( )( ) / T 1 = E g/2 H p lin 0 S p where E 0 = Ψ 0 Ĥ Ψ 0, g i = E(p0 ) p i, H ij = Ψ i Ĥ Ψ j, S ij = Ψ i Ψ j.
12 Optimization method: principle Expansion of the wave function around p 0 to linear order in p = p p 0 : Ψ [1] (p) = Ψ 0 + j p j Ψ j where Ψ 0 = Ψ(p 0 ) and Ψ j = Ψ(p0 )) p j. Normalization of wave function chosen so that the derivatives Ψ j are orthogonal to Ψ 0. Minimization of the energy = generalized eigenvalue equation: ( E0 g T )( ) ( )( ) / T 1 = E g/2 H p lin 0 S p where E 0 = Ψ 0 Ĥ Ψ 0, g i = E(p0 ) p i, H ij = Ψ i Ĥ Ψ j, S ij = Ψ i Ψ j. Update of the parameters: p 0 p 0 + p.
13 Optimization method: robustness The linear method is equivalent to a stabilized Newton method: ( E0 g T )( ) ( )( ) / T 1 = E g/2 H p lin 0 S p { (h + 2 E S) p = g 2 E = g T p where h = 2(H E 0 S) is an approximate Hessian, and E = E 0 E lin > 0 is the energy stabilization. = more robust than Newton method
14 Optimization method: robustness The linear method is equivalent to a stabilized Newton method: ( E0 g T )( ) ( )( ) / T 1 = E g/2 H p lin 0 S p { (h + 2 E S) p = g 2 E = g T p where h = 2(H E 0 S) is an approximate Hessian, and E = E 0 E lin > 0 is the energy stabilization. = more robust than Newton method In quantum chemistry, it is known as super-ci method or augmented Hessian method.
15 Optimization method: robustness The linear method is equivalent to a stabilized Newton method: ( E0 g T )( ) ( )( ) / T 1 = E g/2 H p lin 0 S p { (h + 2 E S) p = g 2 E = g T p where h = 2(H E 0 S) is an approximate Hessian, and E = E 0 E lin > 0 is the energy stabilization. = more robust than Newton method In quantum chemistry, it is known as super-ci method or augmented Hessian method. Additional stabilization: H ij H ij + a δ ij where a 0.
16 Optimization method: on a finite VMC sample The generalized eigenvalue equation is estimated as ( E0 gr T/2 )( ) ( ) ( T 1 = E g L /2 H p lin 0 S p with Ψi (R) g L,i /2 = Ψ 0 (R) Ψi (R) H ij = Ψ 0 (R) H(R)Ψ 0 (R) Ψ 0 (R) H(R)Ψ j (R) Ψ 0 (R) non-symmetric! Ψ 2 0 Ψ 2 0 Ψ0 (R) H(R)Ψ j (R) and g R,j /2 = Ψ 0 (R) Ψ 0 (R) Ψi (R) Ψ j (R) and S ij = Ψ 0 (R) Ψ 0 (R) ) Ψ 2 0 Ψ 2 0
17 Optimization method: on a finite VMC sample The generalized eigenvalue equation is estimated as ( E0 gr T/2 )( ) ( ) ( T 1 = E g L /2 H p lin 0 S p with Ψi (R) g L,i /2 = Ψ 0 (R) Ψi (R) H ij = Ψ 0 (R) H(R)Ψ 0 (R) Ψ 0 (R) H(R)Ψ j (R) Ψ 0 (R) non-symmetric! Ψ 2 0 Ψ 2 0 Ψ0 (R) H(R)Ψ j (R) and g R,j /2 = Ψ 0 (R) Ψ 0 (R) Ψi (R) Ψ j (R) and S ij = Ψ 0 (R) Ψ 0 (R) = Zero-variance principle of Nightingale et al. (PRL 2001): If there is some p so that Ψ 0 (R) + j p j Ψ j (R) = Ψ exact (R) then p is found with zero variance. In practice, these non-symmetric estimators reduce the fluctuations on p by 1 or 2 orders of magnitude. ) Ψ 2 0 Ψ 2 0
18 Optimization method: mixing a fraction of variance How to minimize the energy variance with the linear method? { V = min V 0 + gv T p + 1 } p 2 pt h V p
19 Optimization method: mixing a fraction of variance How to minimize the energy variance with the linear method? { V = min V 0 + gv T p + 1 } p 2 pt h V p V = min p ( 1 p T ) ( V 0 gv T/2 )( 1 g V /2 h V /2 + V 0 S p ( 1 p T ) ( 1 0 T )( ) 1 0 S p )
20 Optimization method: mixing a fraction of variance How to minimize the energy variance with the linear method? { V = min V 0 + gv T p + 1 } p 2 pt h V p V = min p ( 1 p T ) ( V 0 gv T/2 )( 1 g V /2 h V /2 + V 0 S p ( 1 p T ) ( 1 0 T )( ) 1 0 S p ) ( V0 g T V /2 g V /2 h V /2 + V 0 S ) ( 1 p ) ( 1 0 T = V 0 S )( 1 p )
21 Optimization method: mixing a fraction of variance How to minimize the energy variance with the linear method? { V = min V 0 + gv T p + 1 } p 2 pt h V p V = min p ( 1 p T ) ( V 0 gv T/2 )( 1 g V /2 h V /2 + V 0 S p ( 1 p T ) ( 1 0 T )( ) 1 0 S p ) ( V0 g T V /2 g V /2 h V /2 + V 0 S ) ( 1 p matrix to add to the energy matrix ) ( 1 0 T = V 0 S )( 1 p )
22 Simultaneous optimization of all parameters Optimization of 149 parameters = 24 (Jastrow) + 49 (CSF) + 64 (orbitals) + 12 (exponents) for C 2 molecule : Energy (Hartree) Energy (Hartree) Iterations Iterations = Energy converges up to 1 mhartree in a few iterations
23 Systematic improvement in QMC For C 2 molecule: total energies for a series of fully optimized Jastrow-Slater wave functions: Energy (Hartree) VMC CCSD(T)/cc-pVQZ Exact J*SD J*CAS(8,5) J*CAS(8,7) J*CAS(8,8) J*RAS(8,26) Wave function = Systematic improvement in VMC
24 Systematic improvement in QMC For C 2 molecule: total energies for a series of fully optimized Jastrow-Slater wave functions: Energy (Hartree) VMC DMC CCSD(T)/cc-pVQZ Exact J*SD J*CAS(8,5) J*CAS(8,7) J*CAS(8,8) J*RAS(8,26) Wave function = Systematic improvement in VMC and DMC!
25 Potential energy curve of C 2 molecule ( 1 Σ + g ) Jastrow single determinant wave function VMC J SD : Energy (Hartree) Morse potential size-consistency error Interatomic distance (Bohr)
26 Potential energy curve of C 2 molecule ( 1 Σ + g ) Jastrow single determinant wave function VMC J SD : Energy (Hartree) DMC J SD Morse potential size-consistency error Interatomic distance (Bohr) = Single-determinant DMC is size-consistent but with broken spin symmetry at dissociation, Ψ DMC Ŝ 2 Ψ DMC = 2
27 Potential energy curve of C 2 molecule ( 1 Σ + g ) Jastrow multideterminant wave function: Energy (Hartree) VMC J CAS(8,8) DMC J CAS(8,8) Morse potential Interatomic distance (Bohr)
28 Potential energy curve of C 2 molecule ( 1 Σ + g ) Jastrow multideterminant wave function: Energy (Hartree) VMC J CAS(8,8) DMC J CAS(8,8) Morse potential Interatomic distance (Bohr) = DMC gives dissociation energy with chemical accuracy (1 kcal/mol 0.04 ev): D DMC = 6.482(3) vs D exact = 6.44(2) ev
29 Dissociation energies of diatomic molecules Jastrow multideterminant (full valence CAS) wave functions: Errors on dissociation energy (ev) DMC J CAS VMC J CAS MCSCF CAS Li 2 Be 2 B 2 C 2 N 2 O 2 F 2 Ne 2 Molecules = Near chemical accuracy in DMC
30 Example of application Binding energy of 2 NO 2 to a fragment of carbon nanotube: Estimates for the full nanotube (9,0): B3LYP calculations: no binding QMC calculations: weak binding ( 10 kcal/mol) Lawson, Bauschlicher, Toulouse, Filippi, Umrigar, Chem. Phys. Lett., 466, 170 (2008)
31 1 Optimization of wave functions 2 Calculation of observables
32 Calculation of an observable in VMC Energy Estimator: E L (R) = H(R)Ψ(R) Ψ(R) Systematic error: δe = O(δΨ 2 ) Variance: σ 2 (E L ) = O(δΨ 2 ) } Quadratic Zero-Variance Zero-Bias property
33 Calculation of an observable in VMC Energy Estimator: E L (R) = H(R)Ψ(R) Ψ(R) Systematic error: δe = O(δΨ 2 ) Variance: σ 2 (E L ) = O(δΨ 2 ) } Quadratic Zero-Variance Zero-Bias property Arbitrary observable Ô (which does not commute with Ĥ) Estimator: O L (R) = O(R)Ψ(R) Ψ(R) } Systematic error: δo = O(δΨ) Quadratic Zero-Variance Variance: σ 2 (O L ) = O(1) Zero-Bias property
34 Zero-Variance Zero-Bias estimators (Assaraf & Caffarel) Based on the Hellmann-Feynman theorem ( ) de λ Ô = dλ λ=0 where E λ = Ψ λ Ĥ + λô Ψ λ
35 Zero-Variance Zero-Bias estimators (Assaraf & Caffarel) Based on the Hellmann-Feynman theorem ( ) de λ Ô = dλ λ=0 one can define an improved estimator: O improved (R) = O(R)Ψ(R) Ψ(R) where E λ = Ψ λ Ĥ + λô Ψ λ + O ZV (R) + O ZB (R) [ H(R)Ψ ] (R) Ψ (R) with the ZV term: O ZV (R) = Ψ E L (R) (R) Ψ(R) and the ZB term: O ZB (R) = 2 [E L (R) E] Ψ (R) Ψ(R)
36 Zero-Variance Zero-Bias estimators (Assaraf & Caffarel) Based on the Hellmann-Feynman theorem ( ) de λ Ô = dλ λ=0 one can define an improved estimator: O improved (R) = O(R)Ψ(R) Ψ(R) where E λ = Ψ λ Ĥ + λô Ψ λ + O ZV (R) + O ZB (R) [ H(R)Ψ ] (R) Ψ (R) with the ZV term: O ZV (R) = Ψ E L (R) (R) Ψ(R) and the ZB term: O ZB (R) = 2 [E L (R) E] Ψ (R) Ψ(R) Quadratic Zero-Variance Zero-Bias property Systematic error: δo improved = O(δΨ 2 + δψ δψ ) Variance: σ 2 (O improved ) = O(δΨ 2 + δψ 2 + δψ δψ )
37 Example of improved QMC estimators Dipole moment of CH molecule ( 2 Π) in VMC: Dipole moment (Debye) 1.9 usual estimator Interatomic distance (Bohr)
38 Example of improved QMC estimators Dipole moment of CH molecule ( 2 Π) in VMC: Dipole moment (Debye) usual estimator improved estimator Interatomic distance (Bohr) = Reduction of statistical uncertainty!
39 Example of improved QMC estimators Correlation hole of C 2 molecule in VMC: 0.5 usual histogram estimator Correlation hole (Bohr -1 ) Interelectronic distance (Bohr)
40 Example of improved QMC estimators Correlation hole of C 2 molecule in VMC: Correlation hole (Bohr -1 ) usual histogram estimator improved estimator Interelectronic distance (Bohr) = Reduction of statistical uncertainty!
41 Summary and perspectives Summary efficient wave function optimization method by minimization of VMC energy near chemical accuracy with compact wave functions improved estimators for observables in QMC Toulouse, Umrigar, JCP 126, (2007) Umrigar, Toulouse, Filippi, Sorella, Hennig, PRL 98, (2007) Toulouse, Assaraf, Umrigar, JCP 26, (2007) Perspectives Toulouse, Umrigar, JCP 128, (2008) optimization by minimization of DMC energy optimization of molecular geometry excited states
Optimization of quantum Monte Carlo wave functions by energy minimization
Optimization of quantum Monte Carlo wave functions by energy minimization Julien Toulouse, Roland Assaraf, Cyrus J. Umrigar Laboratoire de Chimie Théorique, Université Pierre et Marie Curie and CNRS, Paris,
More informationOptimization of quantum Monte Carlo (QMC) wave functions by energy minimization
Optimization of quantum Monte Carlo (QMC) wave functions by energy minimization Julien Toulouse, Cyrus Umrigar, Roland Assaraf Cornell Theory Center, Cornell University, Ithaca, New York, USA. Laboratoire
More informationQuantum Monte Carlo wave functions and their optimization for quantum chemistry
Quantum Monte Carlo wave functions and their optimization for quantum chemistry Julien Toulouse Université Pierre & Marie Curie and CNRS, Paris, France CEA Saclay, SPhN Orme des Merisiers April 2015 Outline
More informationTime-dependent linear-response variational Monte Carlo.
Time-dependent linear-response variational Monte Carlo. Bastien Mussard bastien.mussard@colorado.edu https://mussard.github.io/ Julien Toulouse julien.toulouse@upmc.fr Sorbonne University, Paris (web)
More informationWave Function Optimization and VMC
Wave Function Optimization and VMC Jeremy McMinis The University of Illinois July 26, 2012 Outline Motivation History of Wave Function Optimization Optimization in QMCPACK Multideterminant Example Motivation:
More informationSize-extensive wave functions for QMC A linear-scaling GVB approach
Size-extensive wave functions for QMC A linear-scaling GVB approach Claudia Filippi, University of Twente, The Netherlands Francesco Fracchia, University of Pisa, Italy Claudio Amovilli, University of
More informationTime-dependent linear-response variational Monte Carlo
Time-dependent linear-response variational Monte Carlo Bastien Mussard 1,2, Emanuele Coccia 1,3, Roland Assaraf 1, Matt Otten 4, C. J. Umrigar 4, and Julien Toulouse 1 1 Laboratoire de Chimie Théorique,
More informationQuantum Monte Carlo backflow calculations of benzene dimers
Quantum Monte Carlo backflow calculations of benzene dimers Kathleen Schwarz*, Cyrus Umrigar**, Richard Hennig*** *Cornell University Department of Chemistry, **Cornell University Department of Physics,
More informationDensity-functional theory at noninteger electron numbers
Density-functional theory at noninteger electron numbers Tim Gould 1, Julien Toulouse 2 1 Griffith University, Brisbane, Australia 2 Université Pierre & Marie Curie and CRS, Paris, France July 215 Introduction
More informationAb-initio molecular dynamics for High pressure Hydrogen
Ab-initio molecular dynamics for High pressure Hydrogen Claudio Attaccalite Institut d'electronique, Microélectronique et Nanotechnologie (IEMN), Lille Outline A brief introduction to Quantum Monte Carlo
More informationFixed-Node quantum Monte Carlo for Chemistry
Fixed-Node quantum Monte Carlo for Chemistry Michel Caffarel Lab. Physique et Chimie Quantiques, CNRS-IRSAMC, Université de Toulouse e-mail : caffarel@irsamc.ups-tlse.fr. p.1/29 The N-body problem of Chemistry
More informationarxiv:physics/ v2 [physics.chem-ph] 9 Apr 2005
Electronic Quantum Monte Carlo Calculations of Atomic Forces, Vibrations, and Anharmonicities arxiv:physics/0411209v2 [physics.chem-ph] 9 Apr 2005 Myung Won Lee a), Massimo Mella b), and Andrew M. Rappe
More informationZero-variance zero-bias quantum Monte Carlo estimators of the spherically and system-averaged pair density
THE JOURNAL OF CHEMICAL PHYSICS 126, 244112 2007 Zero-variance zero-bias quantum Monte Carlo estimators of the spherically and system-averaged pair density Julien Toulouse a Cornell Theory Center, Cornell
More informationPseudo-Hermitian eigenvalue equations in linear-response electronic-structure theory
1/11 Pseudo-Hermitian eigenvalue equations in linear-response electronic-structure theory Julien Toulouse Université Pierre & Marie Curie and CNRS, 4 place Jussieu, Paris, France Web page: www.lct.jussieu.fr/pagesperso/toulouse/
More informationResonating Valence Bond wave function with molecular orbitals: application to diatomic molecules
Resonating Valence Bond wave function with molecular orbitals: application to diatomic molecules M. Marchi 1,2, S. Azadi 2, M. Casula 3, S. Sorella 1,2 1 DEMOCRITOS, National Simulation Center, 34014,
More informationRange-separated density-functional theory with long-range random phase approximation
Range-separated density-functional theory with long-range random phase approximation Julien Toulouse 1 Wuming Zhu 1, Andreas Savin 1, János Ángyán2 1 Laboratoire de Chimie Théorique, UPMC Univ Paris 6
More informationChapter 19 Quantum Monte Carlo Calculations of Electronic Excitation Energies: The Case of the Singlet n π (CO) Transition in Acrolein
Chapter 19 Quantum Monte Carlo Calculations of Electronic Excitation Energies: The Case of the Singlet n π (CO) Transition in Acrolein Julien Toulouse, Michel Caffarel, Peter Reinhardt, Philip E. Hoggan,
More informationAdiabatic-connection fluctuation-dissipation density-functional theory based on range separation
Adiabatic-connection fluctuation-dissipation density-functional theory based on range separation Julien Toulouse 1 I. Gerber 2, G. Jansen 3, A. Savin 1, W. Zhu 1, J. Ángyán 4 1 Laboratoire de Chimie Théorique,
More informationVariational Monte Carlo Optimization and Excited States
Variational Monte Carlo Optimization and Excited States Eric Neuscamman August 9, 2018 motivation charge transfer core spectroscopy double excitations the menu aperitif: number counting Jastrows main course:
More informationarxiv: v2 [physics.chem-ph] 11 Mar 2011
Quantum Monte Carlo with Jastrow-Valence-Bond wave functions Benoît Braïda 1, Julien Toulouse 1, Michel Caffarel 2, and C. J. Umrigar 3 1 Laboratoire de Chimie Théorique, Université Pierre et Marie Curie
More informationThéorie de la fonctionnnelle de la densité avec séparation de portée pour les forces de van der Waals
Théorie de la fonctionnnelle de la densité avec séparation de portée pour les forces de van der Waals Julien Toulouse 1 Iann Gerber 2, Georg Jansen 3, Andreas Savin 1, János Ángyán 4 1 Laboratoire de Chimie
More informationMulticonfiguration wave functions for quantum Monte Carlo calculations of first-row diatomic molecules
Multiconfiguration wave functions for quantum Monte Carlo calculations of first-row diatomic molecules Claudia Filippi Laboratory of Atomic and Solid State Physics and Theory Center, Cornell University,
More informationIntroduction to multiconfigurational quantum chemistry. Emmanuel Fromager
Institut de Chimie, Strasbourg, France Page 1 Emmanuel Fromager Institut de Chimie de Strasbourg - Laboratoire de Chimie Quantique - Université de Strasbourg /CNRS M2 lecture, Strasbourg, France. Notations
More informationQMC dissociation energy of the water dimer: Time step errors and backflow calculations
QMC dissociation energy of the water dimer: Time step errors and backflow calculations Idoia G. de Gurtubay and Richard J. Needs TCM group. Cavendish Laboratory University of Cambridge Idoia G. de Gurtubay.
More informationAssessment of range-separated time-dependent density-functional theory for calculating C 6 dispersion coefficients
1/10 Assessment of range-separated time-dependent density-functional theory for calculating C 6 dispersion coefficients Julien Toulouse 1,2, Elisa Rebolini 1, Tim Gould 3, John F. Dobson 3, Prasenjit Seal
More informationDFT basée sur le théorème de fluctuation-dissipation avec séparation de portée pour les interactions de van der Waals
DFT basée sur le théorème de fluctuation-dissipation avec séparation de portée pour les interactions de van der Waals Julien Toulouse 1 Iann Gerber 2, Georg Jansen 3, Andreas Savin 1, János Ángyán 4 1
More informationExtending Kohn-Sham density-functional theory
Extending Kohn-Sham density-functional theory Julien Toulouse Université Pierre & Marie Curie and CNRS, Paris, France Email: julien.toulouse@upmc.fr Web page: www.lct.jussieu.fr/pagesperso/toulouse/ January
More informationOverview of variational and projector Monte Carlo methods
Overview of variational and projector Monte Carlo methods Cyrus Umrigar Physics Department, Cornell University, Ithaca. Email: CyrusUmrigar@cornell.edu GDR Correlation, Marseille, 8-10 July, 2015 Outline
More informationQuantum Monte Carlo methods
Quantum Monte Carlo methods Lubos Mitas North Carolina State University Urbana, August 2006 Lubos_Mitas@ncsu.edu H= 1 2 i i 2 i, I Z I r ii i j 1 r ij E ion ion H r 1, r 2,... =E r 1, r 2,... - ground
More informationIntroduction to Computational Chemistry
Introduction to Computational Chemistry Vesa Hänninen Laboratory of Physical Chemistry Chemicum 4th floor vesa.hanninen@helsinki.fi September 10, 2013 Lecture 3. Electron correlation methods September
More informationQuantum Monte Carlo for excited state calculations
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
More informationSemistochastic Quantum Monte Carlo A Hybrid of Exact Diagonalization and QMC Methods and Optimization of FN-PMC energies and FN-PMC forces
Semistochastic Quantum Monte Carlo A Hybrid of Exact Diagonalization and QMC Methods and Optimization of FN-PMC energies and FN-PMC forces Cyrus Umrigar Physics Department, Cornell University, Ithaca.
More informationQuantum Monte Carlo study of the ground state of the two-dimensional Fermi fluid
Quantum Monte Carlo study of the ground state of the two-dimensional Fermi fluid N. D. Drummond and R. J. Needs TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge
More informationUncertainty in Molecular Photoionization!
Uncertainty in Molecular Photoionization! Robert R. Lucchese! Department of Chemistry! Texas A&M University Collaborators:! At Texas A&M: R. Carey, J. Lopez, J. Jose! At ISMO, Orsay, France: D. Dowek and
More informationElectronic structure quantum Monte Carlo methods and variable spins: beyond fixedphase/node
Electronic structure quantum Monte Carlo methods and variable spins: beyond fixedphase/node approximations Cody Melton, M. Chandler Bennett, L. Mitas, with A. Ambrosetti, F. Pederiva North Carolina State
More informationQMC dissociation energies of three-electron hemibonded radical cation dimers... and water clusters
QMC dissociation energies of three-electron hemibonded radical cation dimers... and water clusters Idoia G. de Gurtubay TCM group. Cavendish Laboratory University of Cambridge Idoia G. de Gurtubay. Quantum
More informationMethods for calculating forces within quantum Monte Carlo
Methods for calculating forces within quantum Monte Carlo A Badinski 1,2, P D Haynes 1,3, J R Trail 1,4, and R J Needs 1 1 Theory of Condensed Matter Group, Cavendish Laboratory, Cambridge CB3 0HE, UK
More informationLecture 9. Hartree Fock Method and Koopman s Theorem
Lecture 9 Hartree Fock Method and Koopman s Theorem Ψ(N) is approximated as a single slater determinant Φ of N orthogonal One electron spin-orbitals. One electron orbital φ i = φ i (r) χ i (σ) χ i (σ)
More informationOrbital-dependent backflow transformations in quantum Monte Carlo
transformations in quantum Monte Carlo P. Seth, P. López Ríos, and R. J. Needs TCM group, Cavendish Laboratory, University of Cambridge 5 December 2012 VMC and DMC Optimization Wave functions Variational
More informationQuantum Monte Carlo Benchmarks Density Functionals: Si Defects
Quantum Monte Carlo Benchmarks Density Functionals: Si Defects K P Driver, W D Parker, R G Hennig, J W Wilkins (OSU) C J Umrigar (Cornell), R Martin, E Batista, B Uberuaga (LANL), J Heyd, G Scuseria (Rice)
More informationIntroduction to Path Integral Monte Carlo. Part I.
Introduction to Path Integral Monte Carlo. Part I. Alexey Filinov, Jens Böning, Michael Bonitz Institut für Theoretische Physik und Astrophysik, Christian-Albrechts-Universität zu Kiel, D-24098 Kiel, Germany
More informationExercise 1: Structure and dipole moment of a small molecule
Introduction to computational chemistry Exercise 1: Structure and dipole moment of a small molecule Vesa Hänninen 1 Introduction In this exercise the equilibrium structure and the dipole moment of a small
More informationDepartment of Physics and NCSA University of Illinois, Urbana-Champaign, IL, 61801, USA
UNDERSTANDING ELECTRONIC WAVE FUNCTIONS D. M. Ceperley Department of Physics and NCSA University of Illinois, Urbana-Champaign, IL, 61801, USA INTRODUCTION In this article I discuss some aspects of what
More informationAdvanced Electronic Structure Theory Density functional theory. Dr Fred Manby
Advanced Electronic Structure Theory Density functional theory Dr Fred Manby fred.manby@bris.ac.uk http://www.chm.bris.ac.uk/pt/manby/ 6 Strengths of DFT DFT is one of many theories used by (computational)
More informationCoupled-Cluster Perturbative Triples for Bond Breaking
Coupled-Cluster Perturbative Triples for Bond Breaking Andrew G. Taube and Rodney J. Bartlett Quantum Theory Project University of Florida INT CC Meeting Seattle July 8, 2008 Why does chemistry need triples?
More informationTowards accurate all-electron quantum Monte Carlo calculations of transition-metal systems: Spectroscopy of the copper atom
THE JOURNAL OF CHEMICAL PHYSICS 123, 094102 2005 Towards accurate all-electron quantum Monte Carlo calculations of transition-metal systems: Spectroscopy of the copper atom Michel Caffarel, a Jean-Pierre
More informationLinear response time-dependent density functional theory
Linear response time-dependent density functional theory Emmanuel Fromager Laboratoire de Chimie Quantique, Université de Strasbourg, France fromagere@unistra.fr Emmanuel Fromager (UdS) Cours RFCT, Strasbourg,
More informationNodal surfaces and interdimensional degeneracies
THE JOURNAL OF CHEMICAL PHYSICS 142, 214112 (2015) Nodal surfaces and interdimensional degeneracies Pierre-François Loos 1,a) and Dario Bressanini 2,b) 1 Research School of Chemistry, Australian National
More informationNoncollinear spins in QMC: spiral Spin Density Waves in the HEG
Noncollinear spins in QMC: spiral Spin Density Waves in the HEG Zoltán Radnai and Richard J. Needs Workshop at The Towler Institute July 2006 Overview What are noncollinear spin systems and why are they
More informationVan der Waals Interactions Between Thin Metallic Wires and Layers
Van der Waals Interactions Between Thin Metallic Wires and Layers N. D. Drummond and R. J. Needs TCM Group, Cavendish Laboratory, University of Cambridge, Cambridge CB3 0HE, United Kingdom Quantum Monte
More informationQuantum Monte Carlo calculations of two neutrons in finite volume
Quantum Monte Carlo calculations of two neutrons in finite volume Philipp Klos with J. E. Lynn, I. Tews, S. Gandolfi, A. Gezerlis, H.-W. Hammer, M. Hoferichter, and A. Schwenk Nuclear Physics from Lattice
More informationMETALLIZATION AND DISSOCIATION IN HIGH PRESSURE LIQUID HYDROGEN BY AN EFFICIENT MOLECULAR DYNAMICS WITH QUANTUM MONTE CARLO
METALLIZATION AND DISSOCIATION IN HIGH PRESSURE LIQUID HYDROGEN BY AN EFFICIENT MOLECULAR DYNAMICS WITH QUANTUM MONTE CARLO Author: Guglielmo Mazzola Supervisor: Prof. Sandro Sorella A thesis submitted
More informationMulticonfigurational Quantum Chemistry. Björn O. Roos as told by RL Department of Theoretical Chemistry Chemical Center Lund University Sweden
Multiconfigurational Quantum Chemistry Björn O. Roos as told by RL Department of Theoretical Chemistry Chemical Center Lund University Sweden April 20, 2009 1 The Slater determinant Using the spin-orbitals,
More informationAb-initio simulation of liquid water by quantum Monte Carlo
Ab-initio simulation of liquid water by quantum Monte Carlo Sandro Sorella G. Mazzola & Y. Luo SISSA, IOM DEMOCRITOS, Trieste A. Zen, L. Guidoni U. of L Aquila, L Aquila 28 July 2014, Mike Towler Institute,
More informationAll-electron quantum Monte Carlo calculations for the noble gas atoms He to Xe
All-electron quantum Monte Carlo calculations for the noble gas atoms He to Xe A. Ma, N. D. Drummond, M. D. Towler, and R. J. Needs Theory of Condensed Matter Group, Cavendish Laboratory, University of
More informationQuantum Mechanical Simulations
Quantum Mechanical Simulations Prof. Yan Wang Woodruff School of Mechanical Engineering Georgia Institute of Technology Atlanta, GA 30332, U.S.A. yan.wang@me.gatech.edu Topics Quantum Monte Carlo Hartree-Fock
More informationAuxiliary-field quantum Monte Carlo calculations of excited states and strongly correlated systems
Auxiliary-field quantum Monte Carlo calculations of excited states and strongly correlated systems Formally simple -- a framework for going beyond DFT? Random walks in non-orthogonal Slater determinant
More informationCalculating potential energy curves with fixed-node diffusion Monte Carlo: CO and N 2
1 Calculating potential energy curves with fixed-node diffusion Monte Carlo: CO and N 2 Andrew D. Powell and Richard Dawes* Department of Chemistry, Missouri University of Science and Technology, Rolla,
More informationDensity Functional Theory
Chemistry 380.37 Fall 2015 Dr. Jean M. Standard October 28, 2015 Density Functional Theory What is a Functional? A functional is a general mathematical quantity that represents a rule to convert a function
More informationCombining density-functional theory and many-body methods
Combining density-functional theory and many-body methods Julien Toulouse Université Pierre & Marie Curie and CNRS, Paris, France Vrije Universiteit Amsterdam, Netherlands November 2017 Outline 2/23 1
More informationChemistry 483 Lecture Topics Fall 2009
Chemistry 483 Lecture Topics Fall 2009 Text PHYSICAL CHEMISTRY A Molecular Approach McQuarrie and Simon A. Background (M&S,Chapter 1) Blackbody Radiation Photoelectric effect DeBroglie Wavelength Atomic
More informationJack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah
1. Born-Oppenheimer approx.- energy surfaces 2. Mean-field (Hartree-Fock) theory- orbitals 3. Pros and cons of HF- RHF, UHF 4. Beyond HF- why? 5. First, one usually does HF-how? 6. Basis sets and notations
More informationVariational Methods for Electronic Structure
Variational Methods for Electronic Structure The hydrogen atom is a two-body system consisting of a proton and an electron. If spin and relativistic effects are ignored, then the Schrödinger equation for
More informationGround-state properties, excitations, and response of the 2D Fermi gas
Ground-state properties, excitations, and response of the 2D Fermi gas Introduction: 2D FG and a condensed matter perspective Auxiliary-field quantum Monte Carlo calculations - exact* here Results on spin-balanced
More informationThe Nature of the Interlayer Interaction in Bulk. and Few-Layer Phosphorus
Supporting Information for: The Nature of the Interlayer Interaction in Bulk and Few-Layer Phosphorus L. Shulenburger, A.D. Baczewski, Z. Zhu, J. Guan, and D. Tománek, Sandia National Laboratories, Albuquerque,
More informationThe Overhauser Instability
The Overhauser Instability Zoltán Radnai and Richard Needs TCM Group ESDG Talk 14th February 2007 Typeset by FoilTEX Introduction Hartree-Fock theory and Homogeneous Electron Gas Noncollinear spins and
More informationarxiv: v1 [cond-mat.str-el] 17 Aug 2016
arxiv:1608.04901v1 [cond-mat.str-el] 17 Aug 2016 Evolutionary algorithm based configuration interaction approach Rahul Chakraborty 1 1, a) and Debashree Ghosh Physical and Materials Chemistry Division,
More informationContinuum variational and diffusion quantum Monte Carlo calculations
Continuum variational and diffusion quantum Monte Carlo calculations R J Needs, M D Towler, N D Drummond and P López Ríos Theory of Condensed Matter Group, Cavendish Laboratory, Cambridge CB3 0HE, UK Abstract.
More informationBasis sets for electron correlation
Basis sets for electron correlation Trygve Helgaker Centre for Theoretical and Computational Chemistry Department of Chemistry, University of Oslo, Norway The 12th Sostrup Summer School Quantum Chemistry
More informationRandom Estimates in QMC
Page 1 Random Estimates in QMC J.R. Trail and R.J. Needs Theory of Condensed Matter Group, Cavendish Laboratory, University of Cambridge, UK February 2007 Page 2 Numerical Integration Unkown function,
More informationSpin-crossover molecules: puzzling systems for electronic structure methods. Andrea Droghetti School of Physics and CRANN, Trinity College Dublin
Spin-crossover molecules: puzzling systems for electronic structure methods Andrea Droghetti School of Physics and CRANN, Trinity College Dublin Introduction Introduction Can we read a molecule spin? Can
More informationSet the initial conditions r i. Update neighborlist. Get new forces F i
v Set the initial conditions r i ( t 0 ), v i ( t 0 ) Update neighborlist Quantum mechanical models Get new forces F i ( r i ) Solve the equations of motion numerically over time step Δt : r i ( t n )
More informationChemistry 881 Lecture Topics Fall 2001
Chemistry 881 Lecture Topics Fall 2001 Texts PHYSICAL CHEMISTRY A Molecular Approach McQuarrie and Simon MATHEMATICS for PHYSICAL CHEMISTRY, Mortimer i. Mathematics Review (M, Chapters 1,2,3 & 4; M&S,
More information(1/2) M α 2 α, ˆTe = i. 1 r i r j, ˆV NN = α>β
Chemistry 26 Spectroscopy Week # The Born-Oppenheimer Approximation, H + 2. Born-Oppenheimer approximation As for atoms, all information about a molecule is contained in the wave function Ψ, which is the
More informationAccurate barrier heights using diffusion Monte Carlo
Accurate barrier heights using diffusion Monte Carlo Kittithat Krongchon, Brian Busemeyer, and Lucas K. Wagner Department of Physics, University of Illinois at Urbana-Champaign (Dated: March 9, 217) Fixed
More informationAn Overview of Quantum Monte Carlo Methods. David M. Ceperley
An Overview of Quantum Monte Carlo Methods David M. Ceperley Department of Physics and National Center for Supercomputing Applications University of Illinois Urbana-Champaign Urbana, Illinois 61801 In
More informationElectronic structure theory: Fundamentals to frontiers. VI. Analysis and more.
Electronic structure theory: Fundamentals to frontiers. VI. Analysis and more. MARTIN HEAD-GORDON Department of Chemistry, University of California, Berkeley, and, Chemical Sciences Division, Lawrence
More information3: Many electrons. Orbital symmetries. l =2 1. m l
3: Many electrons Orbital symmetries Atomic orbitals are labelled according to the principal quantum number, n, and the orbital angular momentum quantum number, l. Electrons in a diatomic molecule experience
More informationCorrelation in correlated materials (mostly transition metal oxides) Lucas K. Wagner University of Illinois at Urbana-Champaign
Correlation in correlated materials (mostly transition metal oxides) Lucas K. Wagner University of Illinois at Urbana-Champaign Understanding of correlated materials is mostly phenomenological FN- DMC
More informationJack Simons Henry Eyring Scientist and Professor Chemistry Department University of Utah
1. Born-Oppenheimer approx.- energy surfaces 2. Mean-field (Hartree-Fock) theory- orbitals 3. Pros and cons of HF- RHF, UHF 4. Beyond HF- why? 5. First, one usually does HF-how? 6. Basis sets and notations
More informationPairing in Cold Atoms and other Applications for Quantum Monte Carlo methods
Physics Procedia 3 (2010) 1397 1410 www.elsevier.com/locate/procedia Pairing in Cold Atoms and other Applications for Quantum Monte Carlo methods M. Bajdich a, J. Kolorenč a, L. Mitas a and P. J. Reynolds
More informationWave-Function Optimization by Least-Squares Fitting of the Exact Wave Function Sampled by Quantum Monte Carlo
Wave-Function Optimization by Least-Squares Fitting of the Exact Wave Function Sampled by Quantum Monte Carlo R. BIANCHI, D. BRESSANINI, P. CREMASCHI, M. MELLA, AND G. MOROSI Dipartimento di Chimica Fisica
More informationThe wavefunction that describes a bonding pair of electrons:
4.2. Molecular Properties from VB Theory a) Bonding and Bond distances The wavefunction that describes a bonding pair of electrons: Ψ b = a(h 1 ) + b(h 2 ) where h 1 and h 2 are HAOs on adjacent atoms
More informationRenner-Teller Effect in Tetra-Atomic Molecules
Groupe de Chimie Théorique du MSME Renner-Teller Effect in Tetra-Atomic Molecules Laurent Jutier, G. Dhont, H. Khalil and C. Léonard jutier@univ-mlv.fr (non linear) Outline General Presentation Structure
More informationthe abdus salam international centre for theoretical physics
united nations educational, scientific and cultural organization the abdus salam international centre for theoretical physics international atomic energy agency SMR 1595-16 Joint DEMOCRITOS - ICTP School
More informationThe search for the nodes of the fermionic ground state
The search for the nodes of the fermionic ground state ψ ( R) = 0 + - + - - + - + Some believe that unicorns are just shy and difficult to find. However, most people believe blue ones do not exist. We
More informationarxiv: v1 [cond-mat.mes-hall] 10 Feb 2010
Quantum Monte Carlo study of the ground state of the two-dimensional Fermi fluid N. D. Drummond and R. J. Needs TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge
More informationSession 1. Introduction to Computational Chemistry. Computational (chemistry education) and/or (Computational chemistry) education
Session 1 Introduction to Computational Chemistry 1 Introduction to Computational Chemistry Computational (chemistry education) and/or (Computational chemistry) education First one: Use computational tools
More informationFractional-charge and Fractional-spin errors in Range-Separated Density-Functional Theory with Random Phase Approximation.
Fractional-charge and Fractional-spin errors in Range-Separated Density-Functional Theory with Random Phase Approximation. Bastien Mussard, Julien Toulouse Laboratoire de Chimie Théorique Institut du Calcul
More informationAnion-π and π-π cooperative interactions
Chapter 5 Anion-π and π-π cooperative interactions 5.1 Introduction The design of selective receptors of anionic species is a very active area of research within supramolecular chemistry due to the potential
More informationQUACCS 2014: Quantum Monte Carlo for Electrons
QUACCS 2014: Quantum Monte Carlo for Electrons Pierre-François Loos Research School of Chemistry, Australian National University, Canberra, ACT 0200, Australia pf.loos@anu.edu.au Abstract The purpose of
More informationThe QMC Petascale Project
The QMC Petascale Project Richard G. Hennig What will a petascale computer look like? What are the limitations of current QMC algorithms for petascale computers? How can Quantum Monte Carlo algorithms
More informationChapter 2 Quantum chemistry using auxiliary field Monte Carlo
Chapter 2 Quantum chemistry using auxiliary field Monte Carlo 1. The Hubbard-Stratonovich Transformation 2. Neuhauser s shifted contour 3. Calculation of forces and PESs 4. Multireference AFMC 5. Examples
More informationCombining density-functional theory and many-body methods
Combining density-functional theory and many-body methods Julien Toulouse Université Pierre & Marie Curie and CNRS, Paris, France Albuquerque New Mexico, USA June 2016 Outline 2/23 1 A brief overview of
More informationDevelopments in Electronic Structure Theory.
Developments in Electronic Structure Theory. Bastien Mussard Laboratoire de Chimie Théorique Institut du Calcul et de la Simulation Sorbonne Universités, Université Pierre et Marie Curie bastien.mussard@upmc.fr
More informationSection 3 Electronic Configurations, Term Symbols, and States
Section 3 Electronic Configurations, Term Symbols, and States Introductory Remarks- The Orbital, Configuration, and State Pictures of Electronic Structure One of the goals of quantum chemistry is to allow
More informationYingwei Wang Computational Quantum Chemistry 1 Hartree energy 2. 2 Many-body system 2. 3 Born-Oppenheimer approximation 2
Purdue University CHM 67300 Computational Quantum Chemistry REVIEW Yingwei Wang October 10, 2013 Review: Prof Slipchenko s class, Fall 2013 Contents 1 Hartree energy 2 2 Many-body system 2 3 Born-Oppenheimer
More information1 Rayleigh-Schrödinger Perturbation Theory
1 Rayleigh-Schrödinger Perturbation Theory All perturbative techniques depend upon a few simple assumptions. The first of these is that we have a mathematical expression for a physical quantity for which
More informationElectronic structure quantum Monte Carlo methods with variable spins and fixed-phase/node approximations
Electronic structure quantum Monte Carlo methods with variable spins and fixed-phase/node approximations C. Melton, M.C. Bennett, L. Mitas with A. Ambrosetti, F. Pederiva North Carolina State University
More informationKevin Driver 1 Shuai Zhang 1 Burkhard Militzer 1 R. E. Cohen 2.
Quantum Monte Carlo Simulations of a Single Iron Impurity in MgO Kevin Driver 1 Shuai Zhang 1 Burkhard Militzer 1 R. E. Cohen 2 1 Department of Earth & Planetary Science University of California, Berkeley
More information