QUANTUM AND THERMAL MOTION IN MOLECULES FROM FIRST-PRINCIPLES
|
|
- Wilfrid Johns
- 6 years ago
- Views:
Transcription
1 QUANTUM AND THERMAL MOTION IN MOLECULES FROM FIRST-PRINCIPLES 1 Tapio T. Rantala, Department of Physics, Tampere University of Technology CONTENTS MOTIVATION PATH INTEGRAL APPROACH TO QUANTUM DYNAMICS AND STATISTICAL PHYSICS USING MONTE CARLO TECHNIQUE DEMONSTRATING THE FEATURES WITH H3 + BEYOND BORN OPPENHEIMER and ZERO-POINT MOTION ELECTRON NUCLEI COUPLING & ENERGETICS CHEMICAL REACTIONS and RELATED... STRONG CORRELATION
2 MOTIVATION 2 Electronic structure is the key quantity to materials properties and related phenomena: Mechanical, thermal, electrical, optical,.... Conventional ab initio / first-principles type methods suffer from laborious description of electron electron correlations (CI, MCHF, DFT-functionals) typically ignore nuclear quantum and thermal dynamics and coupling of electron nuclei dynamics (Born Oppenheimer approximation) give the zero-kelvin description, only.
3 PATH-INTEGRAL DESCRIPTION OF QUANTUM DYNAMICS 3 Time evolution of the wave function is laborious and challenging : path sampling interference of paths... Propagation in imaginary time is in better control. It can be used as formulation quantum statistical physics, and thus, finding the finite temperature electronic structure. Feynman Hibbs, Quantum Mechanics and Path Integrals, (McGraw-Hill, 1965) Feynman R.P., Rev. Mod. Phys. 20, (1948) Feynman R.P., Statistical Mechanics (Westview, Advance Book Classics, 1972)
4 CLASSICAL PATH 4 Let us consider particle dynamics from a to b. Lagrangian formulation of classical mechanics for finding the path/trajectory leads to equations of motion from minimization (extremum) of action Δx, Δt b = (x b,t b ) t b S= L( x,x,t)dt, where the Lagrangian L = T V. a = (x a,t a ) t a δs = 0 d L => => dt x L x = 0 V x =m x For example, the classical action of the free-particle is S = 1 2 m Δx 2 " % $ ' Δt = 12 # Δt & m (x x b a )2 (t b t a ).
5 5 QUANTUM PATH Usually, the most probable quantum path is the classical one, but other paths contribute, too, with a certain probability. Quantum probability of the particle propagation from a to b is P(b,a) = K(b,a) 2, the absolute square of the probability amplitude K. The probability amplitude is the sum over all oscillating phase factors φ of the paths xab as K(b,a) = φ[x ab ] all x ab where the phase is proportional to the action # φ[x(t)] = A exp% i $ S[x(t)] & ( ' PARTICLE Class. mech. Quantum mechanics λ PHOTON Geom. optics Physical optics
6 6 PATH-INTEGRAL Now, let us define the sum over all paths as a path-integral K(b,a) = e (i/ )S[b,a ] b a We call this kernel or propagator or Green s function. In terms of stationary eigenstates it can be written as K(b,a) = φ * n(a)φ n (b) e (i/ )E n (t b t a ) n Dx(t). The kernel satisfies the free-particle Schrödinger equation in space-time {xb,tb}. For example, the free-particle propagator takes now the form 1/2 # m & # K 0 (b,a) = % ( exp im(x x & b a % )2 (. $ 2πi (t b t a )' $ 2 (t b t a ) '
7 7 MIXED STATE DENSITY MATRIX Considering all states φ n (x) of the particle, for the probability p(x) of finding the particle/system in configuration space at x, we have P(x) = p n (x) = 1 φ * n(x)φ n (x) e βe n n Z. β = 1 n kt Now, define the mixed state density matrix (in position presentation) ρ(x',x) = φ * n(x')φ n (x) e βe n. n ( ρ(β) = e βh ). Thus, we find P(x) = 1 Z ρ(x,x) and normalization implies Z = ρ(x,x) dx = Tr(ρ). Expectation values evaluated from A = Tr(ρA) / Z
8 PATH-INTEGRAL EVALUATION OF DENSITY MATRIX 8 Now, compare ρ(x',x) = n φ * n(x')φ n (x) e βe n K(b,a) = φ * n(a)φ n (b) e (i/ )E n (t b t a ) n and in equilibrium (time independent hamiltonian) and tb > ta. Replacing (tb ta) = u by i β or β = i(tb ta)/ (imaginary time period) we obtain ρ(b,a), for which ρ(b,a) β for a time independent hamiltonian. Thus, we can evaluate the density matrix from a path-integral similarly where the imaginary time action is = H b ρ(b,a). Cf. ρ(x b,x a ;β) = S[x(u);β,0]= b 0 K(b,a) t b all x(u) = i H b K(b,a) e ( i/ )S[β,0] Dx(u),! m $ 2 x2 (u)+v(x(u)) " # % & du.
9 MONTE CARLO SAMPLING OF IMAGINARY TIME PATHS 9 For the density operator we can write ρ(β) = e βh = e β/2h e β/2h, if the kinetic and potential energies in the hamiltonian H = T + V commute. This becomes exact at the limit of imaginary time period goes to zero, the high temperature limit, because the potential energy approaches constant in position representation for each imaginary time step. Thus, we can write ρ(r 0,r M ;β) = ρ(r 0,r 1 ; τ)ρ(r 1,r 2 ; τ)... ρ(r M 1,r M ; τ) dr 1 dr 2...dr M 1, where τ = β / M, β = 1 kt and M is called the Trotter number. This allows numerical sampling of the imaginary time paths with a Monte Carlo method.
10 EVALUATION WITH MONTE CARLO 10 Monte Carlo allows straightforward numerical procedure for evaluation of multidimensional integrals. Metropolis Monte Carlo NVT (equilibrium) ensemble now yields mixed state density matrix with almost classical transparency Metropolis N. et al., J. Chem. Phys. 21, 1087, (1953). Ceperley D.M., Rev. Mod. Phys. 67, 279, (1995) and in Monte Carlo and... (Eds. K. Binder and G. Ciccotti, Editrice Compositori, Bologna, Italy 1996) Storer R.G., J. Math. Phys. 9, 964, (1968)
11 A-FEW-QUANTUM-PARTICLES SYSTEMS 11 A COUPLE OF ELECTRONS IN QUANTUM DOTS M. Leino & TTR, Physica Scripta T114, 44 (2004) M. Leino & TTR, Few Body Systems 40, 237 (2007) HYDROGEN ATOMS ON Ni SURFACE M. Leino, J. Nieminen & TTR, Surf. Sci. 600,1860 (2006) M. Leino, I. Kylänpää & TTR, Surf. Sci. 601, 1246 (2007) ELECTRONS AND NUCLEI QUANTUM DYNAMICS I. Kylänpää, M. Leino & TTR, PRA 76, (2007) I. Kylänpää, TTR, J.Chem.Phys. 133, (2010) I. Kylänpää, TTR, J.Chem.Phys. 135, (2011) THREE AND FOUR PARTICLE MOLECULES I. Kylänpää, TTR, PRA 80, , (2009) I. Kylänpää, TTR and DM. Ceperley, PRA 86, , (2012)
12 QUANTUM STATISTICAL PHYSICS PATH INTEGRAL MONTE CARLO APPROACH 12 An ab initio electronic structure approach with FULL ACCOUNT OF CORRELATION, the van der Waals interaction, for example! TEMPERATURE DEPENDENCE BEYOND BORN OPPENHEIMER APPROXIMATION INTERPRETATION OF EQUILIBRIUM DISSOCIATION REACTION Without the exchange interaction, now
13 SIZE AND TEMPERATURE SCALES 13 T / K PIMC PATH INTEGRAL MONTE CARLO APPROACH RT QC Ab initio- MOLDY MOLDY DFT size/atoms
14 QUANTUM / CLASSICAL APPROACHES TO DYNAMICS 14 time dependent MOLECULAR Wave packet approaches! TDDFT Car Parrinello and T > 0 equilibrium DYNAMICS! Metropolis Monte Carlo Rovibrational ab initio MOLDY T = 0 Molecular mechanics approaches DMC ab initio Quantum Chemistry / DFT / semiemp. electronic dyn.: nuclear dyn.: Classical Q Q Q Q Q Classical
15 15 H 3 + MOLECULE Quantum statistics of two electrons and three nuclei (five-particle system) as a function of temperature: Structure and energetics: quantum nature of nuclei pair correlation functions, contact densities,... dissociation temperature Comparison to the data from conventional quantum chemistry. p + e p + p + me, mp e
16 TOTAL ENERGY: FINITE NUCLEAR MASS AND ZERO-POINT ENERGY 16 Total energy of the H3 + ion up to the dissociation temperature. Born Oppenheimer approximation, classical nuclei and quantum nuclei.
17 MOLECULAR GEOMETRY AT LOW TEMPERATURE: ZERO-POINT MOTION 17 Internuclear distance. Quantum nuclei, classical nuclei, with FWHM Snapshot from simulation, projection to xy-plane. Trotter number 2 16.
18 PARTICLE PARTICLE CORRELATIONS: ELECTRON NUCLEI COUPLING 18 Pair correlation functions. Quantum p (solid), classical p (dashed) and Born Oppenheimer (dash-dotted). e e p p p e
19 ENERGETICS AT HIGH TEMPERATURES: DISSOCIATION RECOMBINATION 19 Total energy of the H3 + beyond dissociation temperature. Lowest density, mid density, highest density H+H + H 2 + +H H 2 +H H 3 + Barrier To Linearity T (K)
20 DISSOCIATION RECOMBINATION EQUILIBRIUM REACTION 20 The molecule and its fragments: Number of Monte Carlo Blocks H 3 + BTL H 2 +H + H 2 + +H 2H+H + We also evaluate molecular free energy, entropy and heat capacity! E (units of Hartree) The equilibrium composition of fragments at about 5000 K.
21 21 SOME THERMODYNAMICS Expected total energy or internal energy is E = U = 1 Z (Z) β and Z = e βe n = e βf or F = 1 ( ) β lnz β = 1 kt n All standard thermodynamic quantities and relations can be derived from the partition function Z or free energy F. F V = P F T = S F = U TS du = PdV+ dq dq = TdS E = U = 1 Z n E n e βe n = kt 2 (lnz) T = (βf) β. Department of Physics, TCOMP-EST Tapio Rantala: QUANTUM STATISTICAL PHYSICS APPROACH WITH PIMC NEST Aug 2012
22 PARTITION FUNCTION 22 Numerical integration of lnz(t) = lnz(t 1 )+ E dt k B T 2 gives the partition function. We use the boundary condition Z(0) = 1. T T 1
23 FREE ENERGY 23 F = 1 β lnz Free energy (units of Hartree) T (K) FIG. 4. Helmholtz free energy from Eq. (5)in the units of Hartree. Notations are the same as in Fig. 3.
24 ENTROPY 24 S = U F T, (11) 35 Entropy (units of k B ) T (K) FIG. 5. Entropy from Eq. (11) in the units of k B. Notations are the same as in Fig. 3.
25 MOLECULAR HEAT CAPACITY 25 C V = E T Heat capacity (units of k B ) 9/2 6/2 3/2 Vibration Rotation T (K) FIG. 6. Molecular heat capacity as a function of temperature calculated using the analytical model of this work. The values on the y-axis are given in units of the Boltzmann constant k B.
26 SOME FINAL NOTES AND THOUGHTS 26 Dynamics is always present in molecules and all the constituent particles participate: in zero Kelvin and in finite temperature In finite temperature there are no stable molecules classical dynamics emerge e + p + me, mp, m e + Essay title: Ehrenfest theorem and decoherence Instruction: Find definition/explanation of both from the literature and compare these two as a way from quantum mechanics to classical mechanics
Ayan Chattopadhyay Mainak Mustafi 3 rd yr Undergraduates Integrated MSc Chemistry IIT Kharagpur
Ayan Chattopadhyay Mainak Mustafi 3 rd yr Undergraduates Integrated MSc Chemistry IIT Kharagpur Under the supervision of: Dr. Marcel Nooijen Associate Professor Department of Chemistry University of Waterloo
More informationMD Thermodynamics. Lecture 12 3/26/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
MD Thermodynamics Lecture 1 3/6/18 1 Molecular dynamics The force depends on positions only (not velocities) Total energy is conserved (micro canonical evolution) Newton s equations of motion (second order
More information221A Lecture Notes Path Integral
A Lecture Notes Path Integral Feynman s Path Integral Formulation Feynman s formulation of quantum mechanics using the so-called path integral is arguably the most elegant. It can be stated in a single
More informationModels for Time-Dependent Phenomena
Models for Time-Dependent Phenomena I. Phenomena in laser-matter interaction: atoms II. Phenomena in laser-matter interaction: molecules III. Model systems and TDDFT Manfred Lein p.1 Outline Phenomena
More informationAb initio molecular dynamics and nuclear quantum effects
Ab initio molecular dynamics and nuclear quantum effects Luca M. Ghiringhelli Fritz Haber Institute Hands on workshop density functional theory and beyond: First principles simulations of molecules and
More informationFeynman Path Integrals in Quantum Mechanics
Feynman Path Integrals in Quantum Mechanics Christian Egli October, 2004 Abstract This text is written as a report to the seminar course in theoretical physics at KTH, Stockholm. The idea of this work
More informationNew Models for Aqueous Systems: Construction of Vibrational Wave Functions for use in Monte Carlo Simulations.
New Models for Aqueous Systems: Construction of Vibrational Wave Functions for use in Monte Carlo Simulations. Maria A. Gomez and Lawrence R. Pratt T-12 and CNLS Theoretical Division Los Alamos National
More information3.320: Lecture 19 (4/14/05) Free Energies and physical Coarse-graining. ,T) + < σ > dµ
3.320: Lecture 19 (4/14/05) F(µ,T) = F(µ ref,t) + < σ > dµ µ µ ref Free Energies and physical Coarse-graining T S(T) = S(T ref ) + T T ref C V T dt Non-Boltzmann sampling and Umbrella sampling Simple
More informationAb initio molecular dynamics. Simone Piccinin CNR-IOM DEMOCRITOS Trieste, Italy. Bangalore, 04 September 2014
Ab initio molecular dynamics Simone Piccinin CNR-IOM DEMOCRITOS Trieste, Italy Bangalore, 04 September 2014 What is MD? 1) Liquid 4) Dye/TiO2/electrolyte 2) Liquids 3) Solvated protein 5) Solid to liquid
More informationAb initio Molecular Dynamics Born Oppenheimer and beyond
Ab initio Molecular Dynamics Born Oppenheimer and beyond Reminder, reliability of MD MD trajectories are chaotic (exponential divergence with respect to initial conditions), BUT... With a good integrator
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 informationHigh Energy D 2 Bond from Feynman s Integral Wave Equation
Applying the Scientific Method to Understanding Anomalous Heat Effects: Opportunities and Challenges High Energy D Bond from Feynman s Integral Wave Equation By: Thomas Barnard Sponsored by: Coolesence
More information1 Imaginary Time Path Integral
1 Imaginary Time Path Integral For the so-called imaginary time path integral, the object of interest is exp( τh h There are two reasons for using imaginary time path integrals. One is that the application
More informationPhase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany
Phase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany Preliminaries Learning Goals Phase Equilibria Phase diagrams and classical thermodynamics
More informationCHEM3023: Spins, Atoms and Molecules
CHEM3023: Spins, Atoms and Molecules Lecture 3 The Born-Oppenheimer approximation C.-K. Skylaris Learning outcomes Separate molecular Hamiltonians to electronic and nuclear parts according to the Born-Oppenheimer
More informationHigh-order Chin actions in path integral Monte Carlo
High-order Chin actions in path integral Monte Carlo High-order actions and their applications, Barcelona 2009 Jordi Boronat Departament de Física i Enginyeria Nuclear Universitat Politècnica de Catalunya
More information3. RATE LAW AND STOICHIOMETRY
Page 1 of 39 3. RATE LAW AND STOICHIOMETRY Professional Reference Shelf R3.2 Abbreviated Lecture Notes Full Lecture Notes I. Overview II. Introduction A. The Transition State B. Procedure to Calculate
More informationTSTC Lectures: Theoretical & Computational Chemistry
TSTC Lectures: Theoretical & Computational Chemistry Rigoberto Hernandez Part I: ( background ) This Lecture Post-Modern Classical Mechanics Hamiltonian Mechanics Numerical Considerations Thermodynamics,
More informationExp. 4. Quantum Chemical calculation: The potential energy curves and the orbitals of H2 +
Exp. 4. Quantum Chemical calculation: The potential energy curves and the orbitals of H2 + 1. Objectives Quantum chemical solvers are used to obtain the energy and the orbitals of the simplest molecules
More informationG : Statistical Mechanics
G5.651: Statistical Mechanics Notes for Lecture 1 I. DERIVATION OF THE DISCRETIZED PATH INTEGRAL We begin our discussion of the Feynman path integral with the canonical ensemble. The epressions for the
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 informationLow-temperature isomers of the water hexamer
Low-temperature isomers of the water hexamer Volodymyr Babin and Francesco Paesani Department of Chemistry & Biochemistry, University of California, San Diego, La Jolla, CA 92093-0314 Hundreds of years
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 informationIntroduction to Path Integrals
Introduction to Path Integrals Consider ordinary quantum mechanics of a single particle in one space dimension. Let s work in the coordinate space and study the evolution kernel Ut B, x B ; T A, x A )
More information1 Interaction of Quantum Fields with Classical Sources
1 Interaction of Quantum Fields with Classical Sources A source is a given external function on spacetime t, x that can couple to a dynamical variable like a quantum field. Sources are fundamental in the
More informationPart II: Statistical Physics
Chapter 6: Boltzmann Statistics SDSMT, Physics Fall Semester: Oct. - Dec., 2013 1 Introduction: Very brief 2 Boltzmann Factor Isolated System and System of Interest Boltzmann Factor The Partition Function
More informationProblem #1 30 points Problem #2 30 points Problem #3 30 points Problem #4 30 points Problem #5 30 points
Name ME 5 Exam # November 5, 7 Prof. Lucht ME 55. POINT DISTRIBUTION Problem # 3 points Problem # 3 points Problem #3 3 points Problem #4 3 points Problem #5 3 points. EXAM INSTRUCTIONS You must do four
More information3.320 Lecture 18 (4/12/05)
3.320 Lecture 18 (4/12/05) Monte Carlo Simulation II and free energies Figure by MIT OCW. General Statistical Mechanics References D. Chandler, Introduction to Modern Statistical Mechanics D.A. McQuarrie,
More informationMolecular Dynamics. Park City June 2005 Tully
Molecular Dynamics John Lance Natasa Vinod Xiaosong Dufie Priya Sharani Hongzhi Group: August, 2004 Prelude: Classical Mechanics Newton s equations: F = ma = mq = p Force is the gradient of the potential:
More informationCircumventing the pathological behavior of path-integral Monte Carlo for systems with Coulomb potentials
Circumventing the pathological behavior of path-integral Monte Carlo for systems with Coulomb potentials M. H. Müser and B. J. Berne Department of Chemistry, Columbia University, New York, New York 10027
More informationPath Integrals. Andreas Wipf Theoretisch-Physikalisches-Institut Friedrich-Schiller-Universität, Max Wien Platz Jena
Path Integrals Andreas Wipf Theoretisch-Physikalisches-Institut Friedrich-Schiller-Universität, Max Wien Platz 1 07743 Jena 5. Auflage WS 2008/09 1. Auflage, SS 1991 (ETH-Zürich) I ask readers to report
More informationPhysics 127c: Statistical Mechanics. Application of Path Integrals to Superfluidity in He 4
Physics 17c: Statistical Mechanics Application of Path Integrals to Superfluidity in He 4 The path integral method, and its recent implementation using quantum Monte Carlo methods, provides both an intuitive
More informationPart II: Statistical Physics
Chapter 6: Boltzmann Statistics SDSMT, Physics Fall Semester: Oct. - Dec., 2014 1 Introduction: Very brief 2 Boltzmann Factor Isolated System and System of Interest Boltzmann Factor The Partition Function
More informationIntroduction Statistical Thermodynamics. Monday, January 6, 14
Introduction Statistical Thermodynamics 1 Molecular Simulations Molecular dynamics: solve equations of motion Monte Carlo: importance sampling r 1 r 2 r n MD MC r 1 r 2 2 r n 2 3 3 4 4 Questions How can
More information1 Equal-time and Time-ordered Green Functions
1 Equal-time and Time-ordered Green Functions Predictions for observables in quantum field theories are made by computing expectation values of products of field operators, which are called Green functions
More information(a) What are the probabilities associated with finding the different allowed values of the z-component of the spin after time T?
1. Quantum Mechanics (Fall 2002) A Stern-Gerlach apparatus is adjusted so that the z-component of the spin of an electron (spin-1/2) transmitted through it is /2. A uniform magnetic field in the x-direction
More informationThe Path Integral Formulation of Quantum Mechanics
Based on Quantum Mechanics and Path Integrals by Richard P. Feynman and Albert R. Hibbs, and Feynman s Thesis The Path Integral Formulation Vebjørn Gilberg University of Oslo July 14, 2017 Contents 1 Introduction
More informationEnergy and Forces in DFT
Energy and Forces in DFT Total Energy as a function of nuclear positions {R} E tot ({R}) = E DF T ({R}) + E II ({R}) (1) where E DF T ({R}) = DFT energy calculated for the ground-state density charge-density
More informationThermodynamics & Statistical Mechanics SCQF Level 9, U03272, PHY-3-ThermStat. Thursday 24th April, a.m p.m.
College of Science and Engineering School of Physics H T O F E E U D N I I N V E B R U S I R T Y H G Thermodynamics & Statistical Mechanics SCQF Level 9, U03272, PHY-3-ThermStat Thursday 24th April, 2008
More information2. Thermodynamics. Introduction. Understanding Molecular Simulation
2. Thermodynamics Introduction Molecular Simulations Molecular dynamics: solve equations of motion r 1 r 2 r n Monte Carlo: importance sampling r 1 r 2 r n How do we know our simulation is correct? Molecular
More informationIV. Classical Molecular Dynamics
IV. Classical Molecular Dynamics Basic Assumptions: 1. Born-Oppenheimer Approximation 2. Classical mechanical nuclear motion Unavoidable Additional Approximations: 1. Approximate potential energy surface
More informationmacroscopic view (phenomenological) microscopic view (atomistic) computing reaction rate rate of reactions experiments thermodynamics
Rate Theory (overview) macroscopic view (phenomenological) rate of reactions experiments thermodynamics Van t Hoff & Arrhenius equation microscopic view (atomistic) statistical mechanics transition state
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 informationp-adic Feynman s path integrals
p-adic Feynman s path integrals G.S. Djordjević, B. Dragovich and Lj. Nešić Abstract The Feynman path integral method plays even more important role in p-adic and adelic quantum mechanics than in ordinary
More informationLecture 1: The Equilibrium Green Function Method
Lecture 1: The Equilibrium Green Function Method Mark Jarrell April 27, 2011 Contents 1 Why Green functions? 2 2 Different types of Green functions 4 2.1 Retarded, advanced, time ordered and Matsubara
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 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 informationTransient Phenomena in Quantum Bound States Subjected to a Sudden Perturbation
Symmetry, Integrability and Geometry: Methods and Applications Vol. (5), Paper 3, 9 pages Transient Phenomena in Quantum Bound States Subjected to a Sudden Perturbation Marcos MOSHINSKY and Emerson SADURNÍ
More informationDepartment of Physics PRELIMINARY EXAMINATION 2014 Part I. Short Questions
Department of Physics PRELIMINARY EXAMINATION 2014 Part I. Short Questions Thursday May 15th, 2014, 14-17h Examiners: Prof. A. Clerk, Prof. M. Dobbs, Prof. G. Gervais (Chair), Prof. T. Webb, Prof. P. Wiseman
More informationLECTURES ON QUANTUM MECHANICS
LECTURES ON QUANTUM MECHANICS GORDON BAYM Unitsersity of Illinois A II I' Advanced Bock Progrant A Member of the Perseus Books Group CONTENTS Preface v Chapter 1 Photon Polarization 1 Transformation of
More informationQuantum Mechanics + Open Systems = Thermodynamics? Jochen Gemmer Tübingen,
Quantum Mechanics + Open Systems = Thermodynamics? Jochen Gemmer Tübingen, 09.02.2006 Motivating Questions more fundamental: understand the origin of thermodynamical behavior (2. law, Boltzmann distribution,
More informationS.K. Saikin May 22, Lecture 13
S.K. Saikin May, 007 13 Decoherence I Lecture 13 A physical qubit is never isolated from its environment completely. As a trivial example, as in the case of a solid state qubit implementation, the physical
More informationMonte Carlo simulations of harmonic and anharmonic oscillators in discrete Euclidean time
Monte Carlo simulations of harmonic and anharmonic oscillators in discrete Euclidean time DESY Summer Student Programme, 214 Ronnie Rodgers University of Oxford, United Kingdom Laura Raes University of
More informationFeynman s path integral approach to quantum physics and its relativistic generalization
Feynman s path integral approach to quantum physics and its relativistic generalization Jürgen Struckmeier j.struckmeier@gsi.de, www.gsi.de/ struck Vortrag im Rahmen des Winterseminars Aktuelle Probleme
More informationExpress the transition state equilibrium constant in terms of the partition functions of the transition state and the
Module 7 : Theories of Reaction Rates Lecture 33 : Transition State Theory Objectives After studying this Lecture you will be able to do the following. Distinguish between collision theory and transition
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 informationMCB100A/Chem130 MidTerm Exam 2 April 4, 2013
MCB1A/Chem13 MidTerm Exam 2 April 4, 213 Name Student ID True/False (2 points each). 1. The Boltzmann constant, k b T sets the energy scale for observing energy microstates 2. Atoms with favorable electronic
More informationMulti-Scale Modeling from First Principles
m mm Multi-Scale Modeling from First Principles μm nm m mm μm nm space space Predictive modeling and simulations must address all time and Continuum Equations, densityfunctional space scales Rate Equations
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 information2m + U( q i), (IV.26) i=1
I.D The Ideal Gas As discussed in chapter II, micro-states of a gas of N particles correspond to points { p i, q i }, in the 6N-dimensional phase space. Ignoring the potential energy of interactions, the
More informationQUANTUM MECHANICS. Franz Schwabl. Translated by Ronald Kates. ff Springer
Franz Schwabl QUANTUM MECHANICS Translated by Ronald Kates Second Revised Edition With 122Figures, 16Tables, Numerous Worked Examples, and 126 Problems ff Springer Contents 1. Historical and Experimental
More informationPhysics 408 Final Exam
Physics 408 Final Exam Name You are graded on your work (with partial credit where it is deserved) so please do not just write down answers with no explanation (or skip important steps)! Please give clear,
More informationIntroduction to Theories of Chemical Reactions. Graduate Course Seminar Beate Flemmig FHI
Introduction to Theories of Chemical Reactions Graduate Course Seminar Beate Flemmig FHI I. Overview What kind of reactions? gas phase / surface unimolecular / bimolecular thermal / photochemical What
More informationMolecular Aggregation
Molecular Aggregation Structure Analysis and Molecular Simulation of Crystals and Liquids ANGELO GAVEZZOTTI University of Milano OXFORD UNIVERSITY PRESS Contents PART I FUNDAMENTALS 1 The molecule: structure,
More informationInternational Physics Course Entrance Examination Questions
International Physics Course Entrance Examination Questions (May 2010) Please answer the four questions from Problem 1 to Problem 4. You can use as many answer sheets you need. Your name, question numbers
More information(i) T, p, N Gibbs free energy G (ii) T, p, µ no thermodynamic potential, since T, p, µ are not independent of each other (iii) S, p, N Enthalpy H
Solutions exam 2 roblem 1 a Which of those quantities defines a thermodynamic potential Why? 2 points i T, p, N Gibbs free energy G ii T, p, µ no thermodynamic potential, since T, p, µ are not independent
More informationPhysics PhD Qualifying Examination Part I Wednesday, August 26, 2015
Physics PhD Qualifying Examination Part I Wednesday, August 26, 2015 Name: (please print) Identification Number: STUDENT: Designate the problem numbers that you are handing in for grading in the appropriate
More informationNanoscale Energy Transport and Conversion A Parallel Treatment of Electrons, Molecules, Phonons, and Photons
Nanoscale Energy Transport and Conversion A Parallel Treatment of Electrons, Molecules, Phonons, and Photons Gang Chen Massachusetts Institute of Technology OXFORD UNIVERSITY PRESS 2005 Contents Foreword,
More informationAndré Schleife Department of Materials Science and Engineering
André Schleife Department of Materials Science and Engineering Yesterday you (should have) learned this: http://upload.wikimedia.org/wikipedia/commons/e/ea/ Simple_Harmonic_Motion_Orbit.gif 1. deterministic
More informationPlasma Spectroscopy Inferences from Line Emission
Plasma Spectroscopy Inferences from Line Emission Ø From line λ, can determine element, ionization state, and energy levels involved Ø From line shape, can determine bulk and thermal velocity and often
More informationQuantum Physics III (8.06) Spring 2007 FINAL EXAMINATION Monday May 21, 9:00 am You have 3 hours.
Quantum Physics III (8.06) Spring 2007 FINAL EXAMINATION Monday May 21, 9:00 am You have 3 hours. There are 10 problems, totalling 180 points. Do all problems. Answer all problems in the white books provided.
More informationSimulations of Dense Atomic Hydrogen in the Wigner Crystal Phase
Published in Journal of Physics and Chemistry of Solids, 67 (2006) 2136 Simulations of Dense Atomic Hydrogen in the Wigner Crystal Phase Burkhard Militzer Geophysical Laboratory, Carnegie Institution of
More informationComputational methods: Coupled Electron Ion Monte Carlo Path Integral Monte Carlo Examples: Electron gas Hydrogen
! Computational methods: Coupled Electron Ion Monte Carlo Path Integral Monte Carlo Examples: Electron gas Hydrogen WHO DID THE WORK? Miguel Morales, Livermore Carlo Pierleoni: L Aquila, Italy Jeff McMahon
More informationData Provided: A formula sheet and table of physical constants are attached to this paper.
Data Provided: A formula sheet and table of physical constants are attached to this paper. DEPARTMENT OF PHYSICS AND ASTRONOMY Spring Semester (2016-2017) From Thermodynamics to Atomic and Nuclear Physics
More information9.1 System in contact with a heat reservoir
Chapter 9 Canonical ensemble 9. System in contact with a heat reservoir We consider a small system A characterized by E, V and N in thermal interaction with a heat reservoir A 2 characterized by E 2, V
More informationICCP Project 2 - Advanced Monte Carlo Methods Choose one of the three options below
ICCP Project 2 - Advanced Monte Carlo Methods Choose one of the three options below Introduction In statistical physics Monte Carlo methods are considered to have started in the Manhattan project (1940
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 informationRole of closed paths in the path integral approach of statistical thermodynamics
Journal of Physics: Conference Series PAPER OPEN ACCESS Role of closed paths in the path integral approach of statistical thermodynamics To cite this article: Jean-Pierre Badiali 215 J. Phys.: Conf. Ser.
More information5.74 Introductory Quantum Mechanics II
MIT OpenCourseWare http://ocw.mit.edu 5.74 Introductory Quantum Mechanics II Spring 009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Andrei Tokmakoff,
More information140a Final Exam, Fall 2007., κ T 1 V P. (? = P or V ), γ C P C V H = U + PV, F = U TS G = U + PV TS. T v. v 2 v 1. exp( 2πkT.
4a Final Exam, Fall 27 Data: P 5 Pa, R = 8.34 3 J/kmol K = N A k, N A = 6.2 26 particles/kilomole, T C = T K 273.5. du = TdS PdV + i µ i dn i, U = TS PV + i µ i N i Defs: 2 β ( ) V V T ( ) /dq C? dt P?
More informationCharge and Energy Transfer Dynamits in Molecular Systems
Volkhard May, Oliver Kühn Charge and Energy Transfer Dynamits in Molecular Systems Second, Revised and Enlarged Edition WILEY- VCH WILEY-VCH Verlag GmbH & Co. KGaA Contents 1 Introduction 19 2 Electronic
More information1 Foundations of statistical physics
1 Foundations of statistical physics 1.1 Density operators In quantum mechanics we assume that the state of a system is described by some vector Ψ belonging to a Hilbert space H. If we know the initial
More informationChemistry 2000 Lecture 1: Introduction to the molecular orbital theory
Chemistry 2000 Lecture 1: Introduction to the molecular orbital theory Marc R. Roussel January 5, 2018 Marc R. Roussel Introduction to molecular orbitals January 5, 2018 1 / 24 Review: quantum mechanics
More informationIndependent electrons in an effective potential
ABC of DFT Adiabatic approximation Independent electrons in an effective potential Hartree Fock Density Functional Theory MBPT - GW Density Functional Theory in a nutshell Every observable quantity of
More informationNon equilibrium thermodynamics: foundations, scope, and extension to the meso scale. Miguel Rubi
Non equilibrium thermodynamics: foundations, scope, and extension to the meso scale Miguel Rubi References S.R. de Groot and P. Mazur, Non equilibrium Thermodynamics, Dover, New York, 1984 J.M. Vilar and
More informationSupporting Online Material (1)
Supporting Online Material The density functional theory (DFT) calculations were carried out using the dacapo code (http://www.fysik.dtu.dk/campos), and the RPBE (1) generalized gradient correction (GGA)
More informationTransition Theory Abbreviated Derivation [ A - B - C] # E o. Reaction Coordinate. [ ] # æ Æ
Transition Theory Abbreviated Derivation A + BC æ Æ AB + C [ A - B - C] # E A BC D E o AB, C Reaction Coordinate A + BC æ æ Æ æ A - B - C [ ] # æ Æ æ A - B + C The rate of reaction is the frequency of
More informationHandout 10. Applications to Solids
ME346A Introduction to Statistical Mechanics Wei Cai Stanford University Win 2011 Handout 10. Applications to Solids February 23, 2011 Contents 1 Average kinetic and potential energy 2 2 Virial theorem
More informationMCB100A/Chem130 MidTerm Exam 2 April 4, 2013
MCBA/Chem Miderm Exam 2 April 4, 2 Name Student ID rue/false (2 points each).. he Boltzmann constant, k b sets the energy scale for observing energy microstates 2. Atoms with favorable electronic configurations
More informationElectron in a Box. A wave packet in a square well (an electron in a box) changing with time.
Electron in a Box A wave packet in a square well (an electron in a box) changing with time. Last Time: Light Wave model: Interference pattern is in terms of wave intensity Photon model: Interference in
More informationPreliminary Examination - Day 2 August 16, 2013
UNL - Department of Physics and Astronomy Preliminary Examination - Day August 16, 13 This test covers the topics of Quantum Mechanics (Topic 1) and Thermodynamics and Statistical Mechanics (Topic ). Each
More informationA Brief Introduction to Statistical Mechanics
A Brief Introduction to Statistical Mechanics E. J. Maginn, J. K. Shah Department of Chemical and Biomolecular Engineering University of Notre Dame Notre Dame, IN 46556 USA Monte Carlo Workshop Universidade
More informationFIRST PUBLIC EXAMINATION SUBJECT 3: PHYSICAL CHEMISTRY
CCHE 4273 FIRST PUBLIC EXAMINATION Trinity Term 2005 Preliminary Examination in Chemistry SUBJECT 3: PHYSICAL CHEMISTRY Wednesday June 8 th 2005, 9:30am Time allowed: 2 ½ hours Candidates should answer
More information2. The Schrödinger equation for one-particle problems. 5. Atoms and the periodic table of chemical elements
1 Historical introduction The Schrödinger equation for one-particle problems 3 Mathematical tools for quantum chemistry 4 The postulates of quantum mechanics 5 Atoms and the periodic table of chemical
More informationEffective Actions Approach for Improving Numerical Calculation of Path Integrals
Effective Actions Approach for Improving Numerical Calculation of Path Integrals A. Balaž, A. Bogojević, I. Vidanović, A. Belić Scientific Computing Laboratory Institute of Physics Belgrade Pregrevica
More informationAb-initio modeling of opto-electronic properties of molecules in solvents and in proximity to a semiconductor nanoparticle
Ab-initio modeling of opto-electronic properties of molecules in solvents and in proximity to a semiconductor nanoparticle Alain Delgado (a,b), Stefano Corni (b), Carlo Andrea Rozzi (b) Stefano Pittalis
More informationThe role of specific cations and water entropy on the stability of branched DNA motif structures
The role of specific cations and water entropy on the stability of branched DNA motif structures Tod A Pascal, William A Goddard III, Prabal K Maiti and Nagarajan Vaidehi 1. Figures Supporting Information
More information4.1 LAWS OF MECHANICS - Review
4.1 LAWS OF MECHANICS - Review Ch4 9 SYSTEM System: Moving Fluid Definitions: System is defined as an arbitrary quantity of mass of fixed identity. Surrounding is everything external to this system. Boundary
More informationCHAPTER 6 Quantum Mechanics II
CHAPTER 6 Quantum Mechanics II 6.1 6.2 6.3 6.4 6.5 6.6 6.7 The Schrödinger Wave Equation Expectation Values Infinite Square-Well Potential Finite Square-Well Potential Three-Dimensional Infinite-Potential
More informationModels for Time-Dependent Phenomena
Models for Time-Dependent Phenomena I. Phenomena in laser-matter interaction: atoms II. Phenomena in laser-matter interaction: molecules III. Model systems and TDDFT Manfred Lein p. Outline Phenomena in
More information