Chapter 2 Quantum chemistry using auxiliary field Monte Carlo
|
|
- Anis Brooks
- 5 years ago
- Views:
Transcription
1 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 of applications 1
2 Acknowledgment Shlomit Jacoby Prof. Daniel Neuhauser Prof. Anna Krylov Prof. Jan Martin Support: Israel Science Foundation Lise Meitner Center 2
3 Overview New quantum chemistry method: Works in terms of basis sets A continuous range of approximations from Hartree-Fock to full CI Based on Monte Carlo, but - No fixed nodes! Discuss the formalism Discuss application issues Show few results 3
4 Many-Fermion ground state (2 nd quantization) Hˆ = h T ˆ + ˆT V ˆ 1 2 ρ ρ ρ Electron density (matrix) 4
5 Goals Ground state energy Low lying excited states properties (correlation functions) How to do that rigorously??? 5
6 e β Ĥ φ Boltzmann is Hot! φ = a ψ + a ψ + a ψ = E0ψ 0 0 ae β E ψ 2 ae β β βe ψ 0 0 ae 0 E ψ ae β + βh 0 βh φ He φ E φ e φ 6
7 Groundstate is Cool Energy (Eh) STO-3G Full CI ˆ ˆ βh Φ He Φ E( ) ˆ β β = Hˆ = H E β β Φ e Φ d E ( β) = ( H ˆ E ( β) ) 2 < 0 dβ β β (au) 7 gs HF Determinant Variational!
8 How to do it??? If electrons were non-interacting easy: h T ˆ 1 N 1 Slater Det. ( ) ( ) e β ρ φ φ = φ β φ β N New Slater Det. 8
9 Proof (BCH) T T βh ˆ ρ βh ˆ ρ e φφ ˆˆ ˆ 1 2 φn e c1c2 cn 0 Slater Det. T T T T T T βh ˆ ρ βh ˆ ρ βh ˆ ρ βh ˆ ρ βh ˆ ρ βh ˆ ρ = e cˆ ˆ ˆ 1e e c2e e cne 0 dˆ dˆ dˆ 1 2 N ( ) ( ) ( ) φ β φ β φ β 1 2 N Slater Det. N T T βh ρˆ βh ρˆ βh βh ˆi = ( ) ˆj φi ( β) = ( ) φ ij ij j j j= 1 e c e e c e 9
10 2-body = very hard ˆ T V ˆ e βρ ρ φ φ N = huge # of det. 1 10
11 Slicing time Divide to N slices: Use Thus: βhˆ βhˆ βhˆ βhˆ e = e e e N β = ˆ T ˆ T β ρ V ρ ˆ βh βt ρˆ 2 2 e = e e + O ( β ) N β 1 T 1 T 1 T ˆ β ρˆ Vρˆ β ρˆ Vρˆ β ρˆ Vρˆ ˆ ˆ ˆ T T T βh βt ρ 2 βt ρ 2 βt ρ 2 e = e e e e e e + O( β) 1 2 Goal: Represent 2-body exponential as sum of 1-body exponentials N 11
12 Hubbard-Stratonovich Strategy Fourier transform of a Gaussian is a Gaussian: aρ 2 2 a x iρx e e e dx An operator! a = V β x = σv β V σ V σ β iσ V ρ β 1ρ ρ β e e e d σ 12
13 Problem in evaluation of integrals Consider a multidimensional integral f W 1 T σ W σ 2 (... ) e f d 1 T σ W σ N σ1 σn σ 2 N e d σ Can t do on a grid! Why? Most points are on grid boundaries 1D 2:1 2D 8:1 3D 26:1 1D 2:2 2D 12:4 3D 56:8 13
14 Solution: Monte Carlo sampling The integral f W 1 T σ W σ 2 (... ) e f d 1 T σ W σ N σ1 σn σ 2 N e d σ Is calculated by producing random numbers ( σ1... σ N ) that are Gaussian distributed and summing f W ( σ... σ ) 1 N f ( σ... σ ) 1 N 14
15 Monte Carlo scheme Φ e Hˆβ Φ = Φ U σ I ( β ) Φ W T T T N ( t iv ) ( t iv ) ( t iv ) 1 2 U σ ( β) Φ = e e e I β + σ ρˆ β + σ ρˆ β + σ ρˆ Φ W N 1 2 σ e n = 1 { } = T n β σ V σ n G. Sugiyama and S. E. Koonin, Ann. Phys. N.Y. 168, 1 (1986). 15
16 Exercise A 1D particle is in a potential well V(x) ˆ 1 ˆ 2 H = p + V ( x) 2 Based on the Hubbard-Stratonovich transformation derive a Monte Carlo method to calculate its GS energy Hint: pˆ β σ β 2 ˆ ( ) 2 iσp β = ( ) e ψ x A e e ψ x dσ 1 σ 2 β 2 = A e ψ( x σ β) dσ 16
17 Results for H Correlation energy (ev) β =4 au β =1 au β =2 au Sylvestrelli, Baroni & Car PRL 71, Iterations 17
18 What s causing the trouble? e Hˆ β = e ( h+ ivσ ) T ˆ ρ β W { σ } = 1 i σ T V ˆ ρ β O W ( β ) σ T Vσ β 1 Noise ( β ½ ) Signal ( β) Signal/Noise β ½ 0 Baer, Head-Gordon & Neuhauser, JCP 109, 6219 (1998) 18
19 19 Solution: contour-shift invariance ( ) W H U e β σ β = ( ) ( ) ( ) τ α τ σ τ σ i ( ) W ia d V i d V H U e e e T T β σ τ α σ τ α α β β τ τ β τ τ = Rom, Charutz & Neuhauser, CPL 270, 382 (1997)
20 τ τ τ Time slice stabilization Φ e βh Φ = Φ β τ e Φ τ τ ( ) H τ ( ) { } ( ) is τ = D W e σ τ σ τ S τ = i ln ( K iv( i )) ( ) T + σ α ρ τ β τ e ( τ τ) τ τ Φ Φ τ T σ Vα τ τ δ S τ = 0 δσ τ 20
21 The stabilizing shift α ( τ ) = Φ ( β τ ) ˆ ρ Φ( τ ) Φ ( β ) Φ( 0) when β, τ are very large: α( τ ) = Φ ˆ ρ Φ gs gs Need to know exact GS density! Use HF density as partial stabilizer 21 Baer, Head-Gordon & Neuhauser, JCP 109, 6219 (1998)
22 Successful Application to H 2 Correlation energy (ev) β= 1au β= 2au β= 4au β= 8au E E ( SC AFMC ) corr ( cc pvqz ) corr = 099. ( 2) ev = 110. ev 1.4 au Iterations 22
23 Things to note No uncontrolled approximation is made (no fixed nodes etc.) The auxiliary fields are sampled from a universal distribution W ( ) ( ) = 1 β T σ τ Vσ τ σ e 2 d 0 { } τ 23
24 Some application issues Two codes: Plane waves/pseudopotentials based code Gaussian based code (Shlomit Jacoby): Reads RHF/UHF GAMESS output Can get data from CASSCF/GVB or spin flip methods to perform multireference calculation 24
25 Inversion Barrier of Water H H O H H 25
26 Inversion Barrier of Water Barrier Correlation Energy (ev) MP5 L=12 au L=8 au β (au) MP2 Tarczay et al JCP 110, (1999) Baer, CPL 324, 101 (1999) 26
27 Electronic force on nuclei de ( + δ ) ( ) E R R E R F = lim dr δ R 0 R δ Force variance (statistical error) is infinite: σ ( F ) ( E) 2σ = lim = δ R 0 δ R 27
28 Solution: correlated sampling Auxiliary fields are sampled from a universal Coulomb-Gauss distribution: W 1 β T σ ( τ ) Vσ ( τ ) σ 2 d 0 { } = e τ Straightforward correlated sampling: Accurate potential surface Geometries Vibrational frequencies 28
29 Force variance is finite! St. Dev. Force (au) dr=1e-3 au dr=1e-4 au dr=1e-5 au β (au) 29
30 N 2 : Bond length R e (A) Experiment L=8 au L=12 au β (au) 30
31 N 2 : Harmonic frequency ω e (cm -1 ) L=12 au Experiment L=8 au β (au) 31
32 N 2 Heat of Formation D(kcal/mol) Experiment L=8 au L=12 au β (au) 32
33 The N 2 Potential Curve Potential Energy (ev) R (au) Hartree Fock β = 0.7 au β = 1.3 au β = 2.5 au Gdanitz CPL (1998) Baer, JCP 113, 473 (2000) 33
34 Multireference AFMC ψ gs Ĥ e β n C n Φ n J ( C) = C HC ε { C SC 1} S mm ' β Hˆ β Hˆ = Φm' e Φm H m' m = Φm' He Φm ˆ HC = SCE 34
35 Singlet-Triplet Splitting of CH 2 30 J (kcal/mole) singlet states 1-singlet state β (au) 35
36 Sources of Error in AFMC E(β) - E FCI (E h ) β = 1 au β = 0.5 au β = 0.2 au β = 0.1 au β = 0.01 au H 2 ccp-vdz I=10 5 β > 0 SE β < Also: Basis-set Frozen core Rid of small eigenvalues of V ijkl β (au) 36
37 Well-tempered AFMC E(β) - E FCI (E h ) K = 6 K = 1 H 2 ccp-vdz I = 5x10 5 R HH = 1.2 A Using several determinants considerably improves the performance β (au) 37
38 Simple example: excited states of H 2 Lines: FullCI #Dets: 6 o β = 0 au o β = 5 au 38
39 H 2 states: close up 39
40 Torsional Ethylene bond rapture Potential energy (E h ) β = 0 (RHF) β = 0.1 au β = 0.5 au β = 1.0 au β = 2.0 au β = 3.0 au β = 4.5 au C 2 H G Θ (deg) 40
41 Singlet-Triplet states of Ethylene Potential energy (ev) B 1u 6-31G β = 4.5 Eh A g Torsional angle Θ (deg) 41
42 Single bond-breaking 7 HF H+F Energy above min. (ev) FCI b = 0.1 au b = 1 au b = 2 au b = 3 au b = 4 au Basis set: 6-31** Full CI taken from Dutta and Sherrill, J. Chem. Phys. 118, 1610 (2003) R (Angstrom) 42
43 H 2 O PES for double bond breaking E-EFCI-min (ev) Basis set: DVZ 8 Detereminants β = 0 β = 0.5 au β = 1.5 au β = 3 au Full-CI AFMC results at various values of β, vs. Full-CI R eq =0.96Å R/Req 43
44 Summary New electronic structure method Formally exact, but has statistical error Can give PES s, break bonds, compute excited states No fixed nodes approximation Balance mutireference and MC: method to select a multi-reference space 44
Coupled-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 informationAuxiliary-field Monte Carlo methods in Fock space: sign problems and methods to circumvent them
Auxiliary-field Monte Carlo methods in Fock space: sign problems and methods to circumvent them Introduction Yoram Alhassid (Yale University) Finite-temperature auxiliary-field Monte Carlo methods in Fock
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 informationOVERVIEW OF QUANTUM CHEMISTRY METHODS
OVERVIEW OF QUANTUM CHEMISTRY METHODS Outline I Generalities Correlation, basis sets Spin II Wavefunction methods Hartree-Fock Configuration interaction Coupled cluster Perturbative methods III Density
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 informationElectronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory
Electronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory MARTIN HEAD-GORDON, Department of Chemistry, University of California, and Chemical Sciences Division, Lawrence Berkeley National
More informationElectron Correlation - Methods beyond Hartree-Fock
Electron Correlation - Methods beyond Hartree-Fock how to approach chemical accuracy Alexander A. Auer Max-Planck-Institute for Chemical Energy Conversion, Mülheim September 4, 2014 MMER Summerschool 2014
More informationStatistical properties of nuclei by the shell model Monte Carlo method
Statistical properties of nuclei by the shell model Monte Carlo method Introduction Yoram Alhassid (Yale University) Shell model Monte Carlo (SMMC) method Circumventing the odd particle-number sign problem
More informationAuxiliary-field quantum Monte Carlo methods in heavy nuclei
Mika Mustonen and Yoram Alhassid (Yale University) Introduction Auxiliary-field quantum Monte Carlo methods in heavy nuclei Auxiliary-field Monte Carlo (AFMC) methods at finite temperature Sign problem
More informationThe nature of superfluidity in the cold atomic unitary Fermi gas
The nature of superfluidity in the cold atomic unitary Fermi gas Introduction Yoram Alhassid (Yale University) Finite-temperature auxiliary-field Monte Carlo (AFMC) method The trapped unitary Fermi gas
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 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 informationFermions in the unitary regime at finite temperatures from path integral auxiliary field Monte Carlo simulations
Fermions in the unitary regime at finite temperatures from path integral auxiliary field Monte Carlo simulations Aurel Bulgac,, Joaquin E. Drut and Piotr Magierski University of Washington, Seattle, WA
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 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 information4πε. me 1,2,3,... 1 n. H atom 4. in a.u. atomic units. energy: 1 a.u. = ev distance 1 a.u. = Å
H atom 4 E a me =, n=,,3,... 8ε 0 0 π me e e 0 hn ε h = = 0.59Å E = me (4 πε ) 4 e 0 n n in a.u. atomic units E = r = Z n nao Z = e = me = 4πε = 0 energy: a.u. = 7. ev distance a.u. = 0.59 Å General results
More informationShifted-contour auxiliary-field Monte Carlo for molecular electronic structure
JOURNAL OF CHEMICAL PHYSICS VOLUME 109, NUMBER 19 15 NOVEMBER 1998 Shifted-contour auxiliary-field Monte Carlo for molecular electronic structure Naomi Rom, Eyal Fattal, Ashish K. Gupta, Emily A. Carter,
More informationQuantum Theory of Many-Particle Systems, Phys. 540
Quantum Theory of Many-Particle Systems, Phys. 540 Questions about organization Second quantization Questions about last class? Comments? Similar strategy N-particles Consider Two-body operators in Fock
More informationChemistry 3502/4502. Final Exam Part I. May 14, 2005
Advocacy chit Chemistry 350/450 Final Exam Part I May 4, 005. For which of the below systems is = where H is the Hamiltonian operator and T is the kinetic-energy operator? (a) The free particle
More informationA microscopic approach to nuclear dynamics. Cédric Simenel CEA/Saclay, France
A microscopic approach to nuclear dynamics Cédric Simenel CEA/Saclay, France Introduction Quantum dynamics of complex systems (nuclei, molecules, BEC, atomic clusters...) Collectivity: from vibrations
More informationAuxiliary-field quantum Monte Carlo methods for nuclei and cold atoms
Introduction Auxiliary-field quantum Monte Carlo methods for nuclei and cold atoms Yoram Alhassid (Yale University) Auxiliary-field Monte Carlo (AFMC) methods at finite temperature Sign problem and good-sign
More informationTDDFT as a tool in biophysics
TDDFT as a tool in biophysics The primary event in vision Robert Send Universität Karlsruhe 09.09.08 Robert Send TDDFT as a tool in biophysics 09.09.08 1 / 28 Outline 1 Human vision 2 The methods 3 The
More informationFaddeev Random Phase Approximation (FRPA) Application to Molecules
Faddeev Random Phase Approximation (FRPA) Application to Molecules Matthias Degroote Center for Molecular Modeling (CMM) Ghent University INT 2011 Spring Program Fermions from Cold Atoms to Neutron Stars:
More information( R)Ψ el ( r;r) = E el ( R)Ψ el ( r;r)
Born Oppenheimer Approximation: Ĥ el ( R)Ψ el ( r;r) = E el ( R)Ψ el ( r;r) For a molecule with N electrons and M nuclei: Ĥ el What is E el (R)? s* potential surface Reaction Barrier Unstable intermediate
More informationarxiv: v1 [nucl-th] 18 Jan 2018
Nuclear deformation in the configuration-interaction shell model arxiv:181.6175v1 [nucl-th] 18 Jan 218 Y. Alhassid, 1 G.F. Bertsch 2,3 C.N. Gilbreth, 2 and M.T. Mustonen 1 1 Center for Theoretical Physics,
More informationAN INTRODUCTION TO QUANTUM CHEMISTRY. Mark S. Gordon Iowa State University
AN INTRODUCTION TO QUANTUM CHEMISTRY Mark S. Gordon Iowa State University 1 OUTLINE Theoretical Background in Quantum Chemistry Overview of GAMESS Program Applications 2 QUANTUM CHEMISTRY In principle,
More informationMolecular Simulation I
Molecular Simulation I Quantum Chemistry Classical Mechanics E = Ψ H Ψ ΨΨ U = E bond +E angle +E torsion +E non-bond Jeffry D. Madura Department of Chemistry & Biochemistry Center for Computational Sciences
More informationQuantum Monte Carlo. QMC methods in the continuum
Quantum Monte Carlo Premise: need to use simulation techniques to solve manybody quantum problems just as you need them classically. Both the wavefunction and expectation values are determined by the simulations.
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 informationLecture 5: More about one- Final words about the Hartree-Fock theory. First step above it by the Møller-Plesset perturbation theory.
Lecture 5: More about one- determinant wave functions Final words about the Hartree-Fock theory. First step above it by the Møller-Plesset perturbation theory. Items from Lecture 4 Could the Koopmans theorem
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 informationIntroduction to Density Functional Theory
Introduction to Density Functional Theory S. Sharma Institut für Physik Karl-Franzens-Universität Graz, Austria 19th October 2005 Synopsis Motivation 1 Motivation : where can one use DFT 2 : 1 Elementary
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 informationDiffusion Monte Carlo
Diffusion Monte Carlo Notes for Boulder Summer School 2010 Bryan Clark July 22, 2010 Diffusion Monte Carlo The big idea: VMC is a useful technique, but often we want to sample observables of the true ground
More informationElectron-Proton Correlation, Theory, and Tunneling Splittings. Sharon Hammes-Schiffer Penn State University
Nuclear-Electronic Orbital Approach: Electron-Proton Correlation, Multicomponent Density Functional Theory, and Tunneling Splittings Sharon Hammes-Schiffer Penn State University Nuclear Quantum Effects
More informationAn Introduction to Quantum Chemistry and Potential Energy Surfaces. Benjamin G. Levine
An Introduction to Quantum Chemistry and Potential Energy Surfaces Benjamin G. Levine This Week s Lecture Potential energy surfaces What are they? What are they good for? How do we use them to solve chemical
More informationBuilding a wavefunction within the Complete-Active. Cluster with Singles and Doubles formalism: straightforward description of quasidegeneracy
Building a wavefunction within the Complete-Active Active-Space Coupled-Cluster Cluster with Singles and Doubles formalism: straightforward description of quasidegeneracy Dmitry I. Lyakh (Karazin Kharkiv
More informationMulti-reference Density Functional Theory. COLUMBUS Workshop Argonne National Laboratory 15 August 2005
Multi-reference Density Functional Theory COLUMBUS Workshop Argonne National Laboratory 15 August 2005 Capt Eric V. Beck Air Force Institute of Technology Department of Engineering Physics 2950 Hobson
More informationShell model Monte Carlo level density calculations in the rare-earth region
Shell model Monte Carlo level density calculations in the rare-earth region Kadir Has University Istanbul, Turkey Workshop on Gamma Strength and Level Density in Nuclear Physics and Nuclear Technology
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 informationMethods for Treating Electron Correlation CHEM 430
Methods for Treating Electron Correlation CHEM 430 Electron Correlation Energy in the Hartree-Fock approximation, each electron sees the average density of all of the other electrons two electrons cannot
More informationIntroduction to Hartree-Fock Molecular Orbital Theory
Introduction to Hartree-Fock Molecular Orbital Theory C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology Origins of Mathematical Modeling in Chemistry Plato (ca. 428-347
More informationWave function methods for the electronic Schrödinger equation
Wave function methods for the electronic Schrödinger equation Zürich 2008 DFG Reseach Center Matheon: Mathematics in Key Technologies A7: Numerical Discretization Methods in Quantum Chemistry DFG Priority
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 informationDensity Func,onal Theory (Chapter 6, Jensen)
Chem 580: DFT Density Func,onal Theory (Chapter 6, Jensen) Hohenberg- Kohn Theorem (Phys. Rev., 136,B864 (1964)): For molecules with a non degenerate ground state, the ground state molecular energy and
More informationQuantum Chemistry Methods
1 Quantum Chemistry Methods T. Helgaker, Department of Chemistry, University of Oslo, Norway The electronic Schrödinger equation Hartree Fock theory self-consistent field theory basis functions and basis
More informationIs a system of fermions in the crossover BCS-BEC. BEC regime a new type of superfluid?
Is a system of fermions in the crossover BCS-BEC BEC regime a new type of superfluid? Finite temperature properties of a Fermi gas in the unitary regime Aurel Bulgac,, Joaquin E. Drut, Piotr Magierski
More informationSecond quantization. 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. Institut
More informationLecture 5. Hartree-Fock Theory. WS2010/11: Introduction to Nuclear and Particle Physics
Lecture 5 Hartree-Fock Theory WS2010/11: Introduction to Nuclear and Particle Physics Particle-number representation: General formalism The simplest starting point for a many-body state is a system of
More informationTime-dependent DMRG:
The time-dependent DMRG and its applications Adrian Feiguin Time-dependent DMRG: ^ ^ ih Ψ( t) = 0 t t [ H ( t) E ] Ψ( )... In a truncated basis: t=3 τ t=4 τ t=5τ t=2 τ t= τ t=0 Hilbert space S.R.White
More informationQuantum Chemical and Dynamical Tools for Solving Photochemical Problems
2.165430 3.413060 3.889592 9 H 3.413060 2.165430 1.099610 2.165430 3.413060 10 H 3.889592 3.413060 2.165430 1.099610 2.165430 11 H 3.413060 3.889592 3.413060 2.165430 1.099610 12 H 2.165430 3.413060 3.889592
More informationPseudopotentials for hybrid density functionals and SCAN
Pseudopotentials for hybrid density functionals and SCAN Jing Yang, Liang Z. Tan, Julian Gebhardt, and Andrew M. Rappe Department of Chemistry University of Pennsylvania Why do we need pseudopotentials?
More informationSpeed of light c = m/s. x n e a x d x = 1. 2 n+1 a n π a. He Li Ne Na Ar K Ni 58.
Physical Chemistry II Test Name: KEY CHEM 464 Spring 18 Chapters 7-11 Average = 1. / 16 6 questions worth a total of 16 points Planck's constant h = 6.63 1-34 J s Speed of light c = 3. 1 8 m/s ħ = h π
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 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 informationChemistry 3502/4502. Final Exam Part I. May 14, 2005
Chemistry 3502/4502 Final Exam Part I May 14, 2005 1. For which of the below systems is = where H is the Hamiltonian operator and T is the kinetic-energy operator? (a) The free particle (e) The
More informationCE 530 Molecular Simulation
1 CE 530 Molecular Simulation Lecture 14 Molecular Models David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Review Monte Carlo ensemble averaging, no dynamics easy
More informationH 2 in the minimal basis
H 2 in the minimal basis Alston J. Misquitta Centre for Condensed Matter and Materials Physics Queen Mary, University of London January 27, 2016 Overview H 2 : The 1-electron basis. The two-electron basis
More informationRecent advances in quantum Monte Carlo for quantum chemistry: optimization of wave functions and calculation of observables
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
More information6. Auxiliary field continuous time quantum Monte Carlo
6. Auxiliary field continuous time quantum Monte Carlo The purpose of the auxiliary field continuous time quantum Monte Carlo method 1 is to calculate the full Greens function of the Anderson impurity
More information2.4. Quantum Mechanical description of hydrogen atom
2.4. Quantum Mechanical description of hydrogen atom Atomic units Quantity Atomic unit SI Conversion Ang. mom. h [J s] h = 1, 05459 10 34 Js Mass m e [kg] m e = 9, 1094 10 31 kg Charge e [C] e = 1, 6022
More informationThe Hartree-Fock Method
1/41 18. December, 2015 Problem Statement First Quantization Second Quantization References Overview 1. Problem Statement 2. First Quantization 3. Second Quantization 4. 5. References 2/41 Problem Statement
More informationElectron Correlation Methods
Electron Correlation Methods HF method: electron-electron interaction is replaced by an average interaction E HF c = E 0 E HF E 0 exact ground state energy E HF HF energy for a given basis set HF E c
More informationElectron Correlation
Electron Correlation Levels of QM Theory HΨ=EΨ Born-Oppenheimer approximation Nuclear equation: H n Ψ n =E n Ψ n Electronic equation: H e Ψ e =E e Ψ e Single determinant SCF Semi-empirical methods Correlation
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 informationResonating Valence Bond point of view in Graphene
Resonating Valence Bond point of view in Graphene S. A. Jafari Isfahan Univ. of Technology, Isfahan 8456, Iran Nov. 29, Kolkata S. A. Jafari, Isfahan Univ of Tech. RVB view point in graphene /2 OUTLINE
More informationChemistry 3502 Physical Chemistry II (Quantum Mechanics) 3 Credits Fall Semester 2003 Christopher J. Cramer. Lecture 25, November 5, 2003
Chemistry 3502 Physical Chemistry II (Quantum Mechanics) 3 Credits Fall Semester 2003 Christopher J. Cramer Lecture 25, November 5, 2003 (Some material in this lecture has been adapted from Cramer, C.
More informationIntroduction to Vibrational Spectroscopy
Introduction to Vibrational Spectroscopy Harmonic oscillators The classical harmonic oscillator The uantum mechanical harmonic oscillator Harmonic approximations in molecular vibrations Vibrational spectroscopy
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 information2. The interaction between an electron gas and the potential field is given by
68A. Solutions to Exercises II. April 5. Carry out a Hubbard Stratonovich transformation for the following interaction Hamiltonians, factorizing the terms in brackets, and write one or two sentences interpreting
More informationDensity Functional Theory for Electrons in Materials
Density Functional Theory for Electrons in Materials Richard M. Martin Department of Physics and Materials Research Laboratory University of Illinois at Urbana-Champaign 1 Density Functional Theory 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 informationPredictive Computing for Solids and Liquids
Predictive Computing for Solids and Liquids So Hirata Department of Chemistry May 214 Blue Waters Symposium 1 Schrödinger equation for a water molecule 1-particle, 3-dimensional partial differential equation
More informationAb initio calculations for potential energy surfaces. D. Talbi GRAAL- Montpellier
Ab initio calculations for potential energy surfaces D. Talbi GRAAL- Montpellier A theoretical study of a reaction is a two step process I-Electronic calculations : techniques of quantum chemistry potential
More informationQuantum Monte Carlo Methods for Strongly Correlated Electron Systems
This is page 1 Printer: Opaque this Quantum Monte Carlo Methods for Strongly Correlated Electron Systems Shiwei Zhang ABSTRACT We review some of the recent development in quantum Monte Carlo (QMC) methods
More informationConvergence properties of the coupled-cluster method: the accurate calculation of molecular properties for light systems
1 Convergence properties of the coupled-cluster method: the accurate calculation of molecular properties for light systems T. Helgaker Centre for Theoretical and Computational Chemistry, Department of
More informationMO Calculation for a Diatomic Molecule. /4 0 ) i=1 j>i (1/r ij )
MO Calculation for a Diatomic Molecule Introduction The properties of any molecular system can in principle be found by looking at the solutions to the corresponding time independent Schrodinger equation
More informationHartree, Hartree-Fock and post-hf methods
Hartree, Hartree-Fock and post-hf methods MSE697 fall 2015 Nicolas Onofrio School of Materials Engineering DLR 428 Purdue University nonofrio@purdue.edu 1 The curse of dimensionality Let s consider a multi
More informationQuantum chemistry and vibrational spectra
Chapter 3 Quantum chemistry and vibrational spectra This chapter presents the quantum chemical results for the systems studied in this work, FHF (Section 3.) and OHF (Section 3.3). These triatomic anions
More informationINTRODUCTION TO QUANTUM ELECTRODYNAMICS by Lawrence R. Mead, Prof. Physics, USM
INTRODUCTION TO QUANTUM ELECTRODYNAMICS by Lawrence R. Mead, Prof. Physics, USM I. The interaction of electromagnetic fields with matter. The Lagrangian for the charge q in electromagnetic potentials V
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 informationNew Frontiers in Nuclear Structure Theory
New Frontiers in Nuclear Structure Theory From Realistic Interactions to the Nuclear Chart Robert Roth Institut für Kernphysik Technical University Darmstadt Overview Motivation Nucleon-Nucleon Interactions
More informationThe successful wavefunction can be written as a determinant: # 1 (2) # 2 (2) Electrons. This can be generalized to our 2N-electron wavefunction:
T2. CNDO to AM1: The Semiempirical Molecular Orbital Models The discussion in sections T2.1 T2.3 applies also to ab initio molecular orbital calculations. T2.1 Slater Determinants Consider the general
More informationThe shape distribution of nuclear level densities in the shell model Monte Carlo method
The shape distribution of nuclear level densities in the shell model Monte Carlo method Introduction Yoram Alhassid (Yale University) Shell model Monte Carlo (SMMC) method and level densities Nuclear deformation
More informationd 3 r d 3 vf( r, v) = N (2) = CV C = n where n N/V is the total number of molecules per unit volume. Hence e βmv2 /2 d 3 rd 3 v (5)
LECTURE 12 Maxwell Velocity Distribution Suppose we have a dilute gas of molecules, each with mass m. If the gas is dilute enough, we can ignore the interactions between the molecules and the energy will
More informationIntroduction to computational chemistry Exercise I: Structure and electronic energy of a small molecule. Vesa Hänninen
Introduction to computational chemistry Exercise I: Structure and electronic energy of a small molecule Vesa Hänninen 1 Introduction In this exercise the equilibrium structure and the electronic energy
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 informationHomework Problem Set 5 Solutions. E e + H corr (a.u.) a.) Determine the bond dissociation enthalpy of ethane in kcal/mol (1 a.u. = kcal/mol).
Chemistry 380.37 Dr. Jean M. Standard Homework Problem Set 5 Solutions 1. Given below are the sum of electronic and thermal enthalpies, E e + H corr, from Hartree-Fock calculations using a 6-31G(d) basis
More informationBridging Scales Through Wavefunction Analysis
Bridging Scales Through Wavefunction Analysis Felix Plasser Institute for Theoretical Chemistry, University of Vienna Excited States Bridging Scales Marseille, November 7 10, 2016 F. Plasser Wavefunction
More informationPH 451/551 Quantum Mechanics Capstone Winter 201x
These are the questions from the W7 exam presented as practice problems. The equation sheet is PH 45/55 Quantum Mechanics Capstone Winter x TOTAL POINTS: xx Weniger 6, time There are xx questions, for
More informationLecture 5: Harmonic oscillator, Morse Oscillator, 1D Rigid Rotor
Lecture 5: Harmonic oscillator, Morse Oscillator, 1D Rigid Rotor It turns out that the boundary condition of the wavefunction going to zero at infinity is sufficient to quantize the value of energy that
More informationVibronic quantum dynamics of exciton relaxation/trapping in molecular aggregates
Symposium, Bordeaux Vibronic quantum dynamics of exciton relaxation/trapping in molecular aggregates Alexander Schubert Institute of Physical and Theoretical Chemistry, University of Würzburg November
More informationThe Rigorous Calculation of Molecular Properties to Chemical Accuracy. T. Helgaker, Department of Chemistry, University of Oslo, Norway
1 The Rigorous Calculation of Molecular Properties to Chemical Accuracy T. Helgaker, Department of Chemistry, University of Oslo, Norway A. C. Hennum and T. Ruden, University of Oslo, Norway S. Coriani,
More informationOptimization 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 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 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 informationStudy of Ozone in Tribhuvan University, Kathmandu, Nepal. Prof. S. Gurung Central Department of Physics, Tribhuvan University, Kathmandu, Nepal
Study of Ozone in Tribhuvan University, Kathmandu, Nepal Prof. S. Gurung Central Department of Physics, Tribhuvan University, Kathmandu, Nepal 1 Country of the Mt Everest 2 View of the Mt Everest 3 4 5
More informationMean-Field Theory: HF and BCS
Mean-Field Theory: HF and BCS Erik Koch Institute for Advanced Simulation, Jülich N (x) = p N! ' (x ) ' 2 (x ) ' N (x ) ' (x 2 ) ' 2 (x 2 ) ' N (x 2 )........ ' (x N ) ' 2 (x N ) ' N (x N ) Slater determinant
More informationMean-field concept. (Ref: Isotope Science Facility at Michigan State University, MSUCL-1345, p. 41, Nov. 2006) 1/5/16 Volker Oberacker, Vanderbilt 1
Mean-field concept (Ref: Isotope Science Facility at Michigan State University, MSUCL-1345, p. 41, Nov. 2006) 1/5/16 Volker Oberacker, Vanderbilt 1 Static Hartree-Fock (HF) theory Fundamental puzzle: The
More informationΨ g = N g (1s a + 1s b ) Themes. Take Home Message. Energy Decomposition. Multiple model wave fns possible. Begin describing brain-friendly MO
Themes Multiple model wave fns possible Different explanations of bonding Some you have seen before Errors vary (not consistently better or worse) Convenience varies Brain v. algebra v. computer Where
More information