where the sums are over the partcle labels. In general H = p2 2m + V s(r ) V j = V nt (jr, r j j) (5) where V s s the sngle-partcle potental and V nt

Size: px
Start display at page:

Download "where the sums are over the partcle labels. In general H = p2 2m + V s(r ) V j = V nt (jr, r j j) (5) where V s s the sngle-partcle potental and V nt"

Transcription

1 Physcs 543 Quantum Mechancs II Fall 998 Hartree-Fock and the Self-consstent Feld Varatonal Methods In the dscusson of statonary perturbaton theory, I mentoned brey the dea of varatonal approxmaton schemes. The basc dea here s that the varatonal prncple: (h jhj,e h j ) = () j s equvalent to the Schrodnger equaton. In other words, the states j that satsfy ths equaton are egenstates of the Hamltonan. To be more explct, consder the quantty hh, E = d 3 r (r) H (r), E d 3 r (r) (r): (2) If we take the functonal dervatve of ths quantty wth respect to the functon (r) wehave hh, E = [H (r), E (r)] (r) d 3 r: (3) In order for ths to be statonary (the varaton s zero for all possble forms of (r)) we see that (r) must be a soluton to the Schrodnger equaton. We can also take the varaton wth respect to (r) and wll nd that hh,e s also statonary f (r) s a soluton to Schrodnger's equaton. The goal of these notes s to explore some of the mplcatons of usng a varatonal scheme to solve for the states of a many-partcle system. ote: I wll not go through the mathematcs of varatonal physcs and fucntonal dervatves. The methods used here are very smlar to those used n classcal physcs. Ths materal s avalable n most texts on Mathematcal Physcs. 2 Hartree Approxmaton Frst, consder a multpartcle system where the partcles are all dstngushable. The Hamltonan of the system has the form, H = H + 2 6=j V j (4)

2 where the sums are over the partcle labels. In general H = p2 2m + V s(r ) V j = V nt (jr, r j j) (5) where V s s the sngle-partcle potental and V nt s the nteracton potental. In the Hartee Approxmaton, one assumes that the egenstates of the total Hamltonan can be wrtten as a product of sngle partcle states. The varatonal approxmaton s then used to derve an equaton for these sngle partcle states. Explctly, assume the egenstates of H can be wrtten = (r) 2(r2) 3(r3) ::: (6) In ths case, the expectaton value of the total Hamltonan wll have the form d 3 r d 3 r2 d 3 r3::: (r) H + 2 6=j V j A (r) 2(r2):::: (7) Lookng at each of the terms, we see that the sums over H and V j reduce to one and two partcle expectaton values respectvely. For example, consder d 3 r d 3 r2 d 3 r3::: (r) 2(r2)::: H (r) 2(r2):::: (8) Because H only operates on the th partcle, the ntegrals over all the other partcle coordnates are equal to one. Therefore, ths term s smply, d 3 r (r ) H (r ): (9) Smlarly, the nteracton term reduces to d 3 r d 3 r j (r ) j (r j ) V j (r ) j (r j ): () Ths means that the quantty we wsh to be statonary s d 3 r (r ) H (r ) + 2 6=j, E Y d 3 r d 3 r j (r ) j (r j ) V j (r ) j (r j ) d 3 r (r ) (r ): () Takng the functonal dervatve of ths equaton wth respect to m(rm) for the m th partcle, we nd d 3 r Hm m (r m )+ 6=m d 3 r (r ) V m (r ) m (r m ), E m (r m ) A m (rm) =: (2) 2

3 The factor of two n front of the nteracton term drops out because there are two equal terms from the dervatve, m = and m = j. Ths means that each of the sngle partcle states n our guess at the total state must satsfy, H m m (r m 6=m d 3 r (r ) V m (r ) A m (r m )=E m (r m ) (3) Ths s usually called the \Hartree equaton". 3 Hartree Potental and Self-Consstency What we have shown above s that the Hartree approxmaton reduced our -partcle problem to a set of sngle partcle equatons that we know (n theory) how to solve. The nteracton between the partcles s reduced to a sngle potental term of the form V H (r m )= d 3 r j (r )j 2 V nt (jr, r m j) (4) 6=m Ths s the \Hartree potental". One usually thnks about ths potental as sayng that the m th partcle nteracts wth the probablty densty of all the other partcles, 6= m, (m) (r) = 6=m j (r)j 2 (5) through the potental V nt (jr, r2j). If we consder, for example, a coulomb nteracton between each of the partcles, the Hartree potental looks lke the nteracton wth the charge densty due to all the other partcles. V H (r m )=q d 3 r q(m) (r) jr m, rj : (6) The dculty wth solvng ths problem s that the Hartree potental connects the solutons of each of the sngle-partcle states. In other words, ths potental for partcle m depends on the states of all the other - partcles. For ths reason, t s sometmes referred to as the \Self-Consstent" potental. The dculty n solvng for the -partcle state s smply that of ndng a set of sngle partcle states (r) that together are solutons to the Hartree equaton. The approach to solvng ths problem numercally s essentally:. Guess a set of sngle partcle states,. 2. Compute the Hartree potental, V H, for each of the partcles usng the states guessed n (). 3. Solve the sngle-partcle (Hartee) equatons for the states usng V H. 4. If the states from (3) are the same as the states from (), you are done. 3

4 5. If the states from (3) are derent than the states from (), use these as a new guess at the states. 6. Repeat untl convergence (solutons = guesses). There are some subtletes n the process. The most mportant s that after solvng for the sngle partcle states n (3), the new guess at the states s a mxture of the old guess from () and the new solutons. If one uses the new states as the guesses, wld oscllatons n the numercs often occur. Of course, the soluton to the Hartree equaton s not exact. The approxmaton used s that the soluton to the problem can be wrtten as a product of sngle partcle states. Ths s the usual varatonal smplcatons where the regon of Hlbert space consdered s lmted by some guess as to the form of the states. \Exact" solutons to the -partcle problem nvolve lnear combnatons of these Hartee-lke bass states. Ths sort of soluton method has been able to produce exact dagonalzaton of Hamltonans of up to about dozen partcles, dependng on the physcal system. Of course, ths requres a lot of computer. 4 Hartee-Fock and Exchange The Hartree-Fock approxmaton follows drectly from what was done above for the Hatree approxmaton, only takng nto account the need for ant-symmetrc states for Fermons. The Hamltonan for the system s the same as for the Hartree approxmaton. The exchange requrement on the system of partcles wll add an extra nteracton term n the nal result. We start wth the Slater determnant state = p! (n ;n 2 ;:::) n n 2 n 3 ::: n (r) n2 (r2) n3 (r3) ::: (7) where the sum s over all possble permutatons of the ndces (n;n2;n3; :::). (ote that here I'm permutng the ndces labelng the sngle-partcle state functons rather than the ndces of the partcle coordnates as was done prevously. You mght want to convnce yourself that one gets the same states.) Frst, consder the matrx element of one of the sngle-partcle operators, H. Ths wll have the form, n n (n ;n 2 ;:::) (n ;:::) 2 n 3 ::: n n 2 n 3 ::: d 3 r d 3 r2 d 3 r3::: n (r ) n (r 2)::: H s n (r) n2 (r2):::: 2 ;n 2 (8) Because the operator only operates on the th state, ths can be rewrtten as (n ;n 2 ;:::) j6= (n ;n 2 ;:::) n n 2 n 3 ;::: n n 2 n 3 ;::: Y d 3 r j n j (r j ) nj (r j ) d 3 r s n (r ) H s n (r ): (9) 4

5 The sngle partcle states are orthogonal. Ths means that for there to be a non-zero matrx element, the sum on permutatons must have n j = n j for j 6= etc. But ths also means that n = n because that's the only state left. Ths gves only the dagonal terms (! of them) n the double sum over permutatons, and the product n n 2 n 3 ;::: n n 2 n 3 ;::: =. Therefore, the Hartree-Fock result for the sngle-partcle terms s the same as the Hartree result, d 3 r (r ) H (r ): (2) The nteracton potental term s a lttle more problematc. In ths case, we're lookng at a matrx element of the form, n n (n ;n 2 ;:::) (n ;:::) 2 n 3 ;::: n n 2 n 3 ;::: d 3 r d 3 r2 d 3 r3::: n (r ) n (r 2)::: V ;j n (r) n2 (r2):::: 2 ;n 2 (2) Agan, the ntegrals over the r k for k 6= ; j wll lmt the sum over the permutatons to those where n k = n k. For the and j states, there are two possbltes, ether n = n and n j = n j or n = n j and n j = n. Ths gves the matrx element, (n ;n 2 ;:::) d 3 r d 3 r2 V j [ n :::n :::n j ::: n :::n :::n j ::: n (r j ) n (r ) nj (r j ) + n :::n :::n j ::: n :::n j :::n ::: n (r j ) nj (r ) n (r j )]: (22) The product of the 's n the rst term wll equal because the order of the n k 's s the same n each. For the second term the product of the 's wll equal - because two ndces (n and n j )have been swtched. The sum over permutatons just gves a factor of!, so our nteracton potental term reduces to: d 3 r d 3 r j V j j n (r )j 2 j nj (r j )j 2, d 3 r d 3 r j V j n (r j ) nj (r ) n (r j ): (23) The rst term s the same as the Hartree analyss, the second term (sometmes called the Exchange or Fock term) arses from the Paul excluson prncple. The quantty that we need to mnmze n ths case s d 3 r (r ) H (r ) + 2, 2, E 6=j Y 6=j d 3 r d 3 r j n (r j ) V j n (r ) nj (r j ) d 3 r d 3 r j n (r j ) V j nj (r ) n (r j ) d 3 r (r ) (r ): (24) Performng the same mnmzaton process as above, we'll take the functonal dervatve of the terms above wth respect to some state m(r m ). The results, smlar to the Hartree case, are the sngle partcle equatons H m m (r m ) + 6=m d 3 r (r ) V m (r ) m (r m ) 5

6 , 6=m d 3 r (r ) V m (r m ) m (r )=E m (r m ) (25) The nal term s sometmes wrtten, d 3 6=m (r ) V m (r m ) m (r ) A (26) Ths s a \non-local nteracton" term for m. It depends on the value of m everywhere. Usually, nteracton terms are local, of the form V (r) (r). A non-local potental term has the form V (r) (r)! d 3 r U(r; r ) (r ) (27) Because of the ant-symmetrzaton requrement on Fermons, that the states \avod" each other, ths exchange energy depends on the value of the state everywhere. 6

A how to guide to second quantization method.

A how to guide to second quantization method. Phys. 67 (Graduate Quantum Mechancs Sprng 2009 Prof. Pu K. Lam. Verson 3 (4/3/2009 A how to gude to second quantzaton method. -> Second quantzaton s a mathematcal notaton desgned to handle dentcal partcle

More information

The Feynman path integral

The Feynman path integral The Feynman path ntegral Aprl 3, 205 Hesenberg and Schrödnger pctures The Schrödnger wave functon places the tme dependence of a physcal system n the state, ψ, t, where the state s a vector n Hlbert space

More information

Poisson brackets and canonical transformations

Poisson brackets and canonical transformations rof O B Wrght Mechancs Notes osson brackets and canoncal transformatons osson Brackets Consder an arbtrary functon f f ( qp t) df f f f q p q p t But q p p where ( qp ) pq q df f f f p q q p t In order

More information

Some Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS)

Some Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS) Some Comments on Acceleratng Convergence of Iteratve Sequences Usng Drect Inverson of the Iteratve Subspace (DIIS) C. Davd Sherrll School of Chemstry and Bochemstry Georga Insttute of Technology May 1998

More information

12. The Hamilton-Jacobi Equation Michael Fowler

12. The Hamilton-Jacobi Equation Michael Fowler 1. The Hamlton-Jacob Equaton Mchael Fowler Back to Confguraton Space We ve establshed that the acton, regarded as a functon of ts coordnate endponts and tme, satsfes ( ) ( ) S q, t / t+ H qpt,, = 0, and

More information

PHYS 705: Classical Mechanics. Calculus of Variations II

PHYS 705: Classical Mechanics. Calculus of Variations II 1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary

More information

Density matrix. c α (t)φ α (q)

Density matrix. c α (t)φ α (q) Densty matrx Note: ths s supplementary materal. I strongly recommend that you read t for your own nterest. I beleve t wll help wth understandng the quantum ensembles, but t s not necessary to know t n

More information

THEOREMS OF QUANTUM MECHANICS

THEOREMS OF QUANTUM MECHANICS THEOREMS OF QUANTUM MECHANICS In order to develop methods to treat many-electron systems (atoms & molecules), many of the theorems of quantum mechancs are useful. Useful Notaton The matrx element A mn

More information

Lecture 20: Noether s Theorem

Lecture 20: Noether s Theorem Lecture 20: Noether s Theorem In our revew of Newtonan Mechancs, we were remnded that some quanttes (energy, lnear momentum, and angular momentum) are conserved That s, they are constant f no external

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng

More information

CME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 13

CME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 13 CME 30: NUMERICAL LINEAR ALGEBRA FALL 005/06 LECTURE 13 GENE H GOLUB 1 Iteratve Methods Very large problems (naturally sparse, from applcatons): teratve methods Structured matrces (even sometmes dense,

More information

PHYS 215C: Quantum Mechanics (Spring 2017) Problem Set 3 Solutions

PHYS 215C: Quantum Mechanics (Spring 2017) Problem Set 3 Solutions PHYS 5C: Quantum Mechancs Sprng 07 Problem Set 3 Solutons Prof. Matthew Fsher Solutons prepared by: Chatanya Murthy and James Sully June 4, 07 Please let me know f you encounter any typos n the solutons.

More information

PHYS 705: Classical Mechanics. Hamilton-Jacobi Equation

PHYS 705: Classical Mechanics. Hamilton-Jacobi Equation 1 PHYS 705: Classcal Mechancs Hamlton-Jacob Equaton Hamlton-Jacob Equaton There s also a very elegant relaton between the Hamltonan Formulaton of Mechancs and Quantum Mechancs. To do that, we need to derve

More information

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve

More information

Advanced Quantum Mechanics

Advanced Quantum Mechanics Advanced Quantum Mechancs Rajdeep Sensarma! sensarma@theory.tfr.res.n ecture #9 QM of Relatvstc Partcles Recap of ast Class Scalar Felds and orentz nvarant actons Complex Scalar Feld and Charge conjugaton

More information

Lagrangian Field Theory

Lagrangian Field Theory Lagrangan Feld Theory Adam Lott PHY 391 Aprl 6, 017 1 Introducton Ths paper s a summary of Chapter of Mandl and Shaw s Quantum Feld Theory [1]. The frst thng to do s to fx the notaton. For the most part,

More information

The non-negativity of probabilities and the collapse of state

The non-negativity of probabilities and the collapse of state The non-negatvty of probabltes and the collapse of state Slobodan Prvanovć Insttute of Physcs, P.O. Box 57, 11080 Belgrade, Serba Abstract The dynamcal equaton, beng the combnaton of Schrödnger and Louvlle

More information

Workshop: Approximating energies and wave functions Quantum aspects of physical chemistry

Workshop: Approximating energies and wave functions Quantum aspects of physical chemistry Workshop: Approxmatng energes and wave functons Quantum aspects of physcal chemstry http://quantum.bu.edu/pltl/6/6.pdf Last updated Thursday, November 7, 25 7:9:5-5: Copyrght 25 Dan Dll (dan@bu.edu) Department

More information

ψ ij has the eigenvalue

ψ ij has the eigenvalue Moller Plesset Perturbaton Theory In Moller-Plesset (MP) perturbaton theory one taes the unperturbed Hamltonan for an atom or molecule as the sum of the one partcle Foc operators H F() where the egenfunctons

More information

PHYS 705: Classical Mechanics. Newtonian Mechanics

PHYS 705: Classical Mechanics. Newtonian Mechanics 1 PHYS 705: Classcal Mechancs Newtonan Mechancs Quck Revew of Newtonan Mechancs Basc Descrpton: -An dealzed pont partcle or a system of pont partcles n an nertal reference frame [Rgd bodes (ch. 5 later)]

More information

14 The Postulates of Quantum mechanics

14 The Postulates of Quantum mechanics 14 The Postulates of Quantum mechancs Postulate 1: The state of a system s descrbed completely n terms of a state vector Ψ(r, t), whch s quadratcally ntegrable. Postulate 2: To every physcally observable

More information

Robert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations

Robert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations Quantum Physcs 量 理 Robert Esberg Second edton CH 09 Multelectron atoms ground states and x-ray exctatons 9-01 By gong through the procedure ndcated n the text, develop the tme-ndependent Schroednger equaton

More information

Snce h( q^; q) = hq ~ and h( p^ ; p) = hp, one can wrte ~ h hq hp = hq ~hp ~ (7) the uncertanty relaton for an arbtrary state. The states that mnmze t

Snce h( q^; q) = hq ~ and h( p^ ; p) = hp, one can wrte ~ h hq hp = hq ~hp ~ (7) the uncertanty relaton for an arbtrary state. The states that mnmze t 8.5: Many-body phenomena n condensed matter and atomc physcs Last moded: September, 003 Lecture. Squeezed States In ths lecture we shall contnue the dscusson of coherent states, focusng on ther propertes

More information

Lecture 12: Discrete Laplacian

Lecture 12: Discrete Laplacian Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly

More information

Lecture 14: Forces and Stresses

Lecture 14: Forces and Stresses The Nuts and Bolts of Frst-Prncples Smulaton Lecture 14: Forces and Stresses Durham, 6th-13th December 2001 CASTEP Developers Group wth support from the ESF ψ k Network Overvew of Lecture Why bother? Theoretcal

More information

EPR Paradox and the Physical Meaning of an Experiment in Quantum Mechanics. Vesselin C. Noninski

EPR Paradox and the Physical Meaning of an Experiment in Quantum Mechanics. Vesselin C. Noninski EPR Paradox and the Physcal Meanng of an Experment n Quantum Mechancs Vesseln C Nonnsk vesselnnonnsk@verzonnet Abstract It s shown that there s one purely determnstc outcome when measurement s made on

More information

XII. The Born-Oppenheimer Approximation

XII. The Born-Oppenheimer Approximation X. The Born-Oppenhemer Approxmaton The Born- Oppenhemer (BO) approxmaton s probably the most fundamental approxmaton n chemstry. From a practcal pont of vew t wll allow us to treat the ectronc structure

More information

A particle in a state of uniform motion remain in that state of motion unless acted upon by external force.

A particle in a state of uniform motion remain in that state of motion unless acted upon by external force. The fundamental prncples of classcal mechancs were lad down by Galleo and Newton n the 16th and 17th centures. In 1686, Newton wrote the Prncpa where he gave us three laws of moton, one law of gravty,

More information

Non-interacting Spin-1/2 Particles in Non-commuting External Magnetic Fields

Non-interacting Spin-1/2 Particles in Non-commuting External Magnetic Fields EJTP 6, No. 0 009) 43 56 Electronc Journal of Theoretcal Physcs Non-nteractng Spn-1/ Partcles n Non-commutng External Magnetc Felds Kunle Adegoke Physcs Department, Obafem Awolowo Unversty, Ile-Ife, Ngera

More information

International Journal of Pure and Applied Sciences and Technology

International Journal of Pure and Applied Sciences and Technology Int. J. Pure Appl. Sc. Technol., 4() (03), pp. 5-30 Internatonal Journal of Pure and Appled Scences and Technology ISSN 9-607 Avalable onlne at www.jopaasat.n Research Paper Schrödnger State Space Matrx

More information

Lagrange Multipliers. A Somewhat Silly Example. Monday, 25 September 2013

Lagrange Multipliers. A Somewhat Silly Example. Monday, 25 September 2013 Lagrange Multplers Monday, 5 September 013 Sometmes t s convenent to use redundant coordnates, and to effect the varaton of the acton consstent wth the constrants va the method of Lagrange undetermned

More information

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1 P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the

More information

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law: CE304, Sprng 2004 Lecture 4 Introducton to Vapor/Lqud Equlbrum, part 2 Raoult s Law: The smplest model that allows us do VLE calculatons s obtaned when we assume that the vapor phase s an deal gas, and

More information

Quantum Statistical Mechanics

Quantum Statistical Mechanics Chapter 8 Quantum Statstcal Mechancs 8.1 Mcrostates For a collecton of N partcles the classcal mcrostate of the system s unquely specfed by a vector (q, p) (q 1...q N, p 1...p 3N )(q 1...q 3N,p 1...p 3N

More information

Errors in Nobel Prize for Physics (7) Improper Schrodinger Equation and Dirac Equation

Errors in Nobel Prize for Physics (7) Improper Schrodinger Equation and Dirac Equation Errors n Nobel Prze for Physcs (7) Improper Schrodnger Equaton and Drac Equaton u Yuhua (CNOOC Research Insttute, E-mal:fuyh945@sna.com) Abstract: One of the reasons for 933 Nobel Prze for physcs s for

More information

Probabilistic method to determine electron correlation energy

Probabilistic method to determine electron correlation energy Probablstc method to determne electron elaton energy T.R.S. Prasanna Department of Metallurgcal Engneerng and Materals Scence Indan Insttute of Technology, Bombay Mumba 400076 Inda A new method to determne

More information

Advanced Circuits Topics - Part 1 by Dr. Colton (Fall 2017)

Advanced Circuits Topics - Part 1 by Dr. Colton (Fall 2017) Advanced rcuts Topcs - Part by Dr. olton (Fall 07) Part : Some thngs you should already know from Physcs 0 and 45 These are all thngs that you should have learned n Physcs 0 and/or 45. Ths secton s organzed

More information

Physics 181. Particle Systems

Physics 181. Particle Systems Physcs 181 Partcle Systems Overvew In these notes we dscuss the varables approprate to the descrpton of systems of partcles, ther defntons, ther relatons, and ther conservatons laws. We consder a system

More information

Yukawa Potential and the Propagator Term

Yukawa Potential and the Propagator Term PHY304 Partcle Physcs 4 Dr C N Booth Yukawa Potental an the Propagator Term Conser the electrostatc potental about a charge pont partcle Ths s gven by φ = 0, e whch has the soluton φ = Ths escrbes the

More information

5.03, Inorganic Chemistry Prof. Daniel G. Nocera Lecture 2 May 11: Ligand Field Theory

5.03, Inorganic Chemistry Prof. Daniel G. Nocera Lecture 2 May 11: Ligand Field Theory 5.03, Inorganc Chemstry Prof. Danel G. Nocera Lecture May : Lgand Feld Theory The lgand feld problem s defned by the followng Hamltonan, h p Η = wth E n = KE = where = m m x y z h m Ze r hydrogen atom

More information

Integrals and Invariants of Euler-Lagrange Equations

Integrals and Invariants of Euler-Lagrange Equations Lecture 16 Integrals and Invarants of Euler-Lagrange Equatons ME 256 at the Indan Insttute of Scence, Bengaluru Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng,

More information

PHYS 705: Classical Mechanics. Canonical Transformation II

PHYS 705: Classical Mechanics. Canonical Transformation II 1 PHYS 705: Classcal Mechancs Canoncal Transformaton II Example: Harmonc Oscllator f ( x) x m 0 x U( x) x mx x LT U m Defne or L p p mx x x m mx x H px L px p m p x m m H p 1 x m p m 1 m H x p m x m m

More information

The Second Anti-Mathima on Game Theory

The Second Anti-Mathima on Game Theory The Second Ant-Mathma on Game Theory Ath. Kehagas December 1 2006 1 Introducton In ths note we wll examne the noton of game equlbrum for three types of games 1. 2-player 2-acton zero-sum games 2. 2-player

More information

1 Derivation of Point-to-Plane Minimization

1 Derivation of Point-to-Plane Minimization 1 Dervaton of Pont-to-Plane Mnmzaton Consder the Chen-Medon (pont-to-plane) framework for ICP. Assume we have a collecton of ponts (p, q ) wth normals n. We want to determne the optmal rotaton and translaton

More information

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and

More information

Introduction to Density Functional Theory. Jeremie Zaffran 2 nd year-msc. (Nanochemistry)

Introduction to Density Functional Theory. Jeremie Zaffran 2 nd year-msc. (Nanochemistry) Introducton to Densty Functonal Theory Jereme Zaffran nd year-msc. (anochemstry) A- Hartree appromatons Born- Oppenhemer appromaton H H H e The goal of computatonal chemstry H e??? Let s remnd H e T e

More information

Difference Equations

Difference Equations Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1

More information

LAGRANGIAN MECHANICS

LAGRANGIAN MECHANICS LAGRANGIAN MECHANICS Generalzed Coordnates State of system of N partcles (Newtonan vew): PE, KE, Momentum, L calculated from m, r, ṙ Subscrpt covers: 1) partcles N 2) dmensons 2, 3, etc. PE U r = U x 1,

More information

HW #6, due Oct Toy Dirac Model, Wick s theorem, LSZ reduction formula. Consider the following quantum mechanics Lagrangian,

HW #6, due Oct Toy Dirac Model, Wick s theorem, LSZ reduction formula. Consider the following quantum mechanics Lagrangian, HW #6, due Oct 5. Toy Drac Model, Wck s theorem, LSZ reducton formula. Consder the followng quantum mechancs Lagrangan, L ψ(σ 3 t m)ψ, () where σ 3 s a Paul matrx, and ψ s defned by ψ ψ σ 3. ψ s a twocomponent

More information

Physics 5153 Classical Mechanics. Principle of Virtual Work-1

Physics 5153 Classical Mechanics. Principle of Virtual Work-1 P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal

More information

Linear Approximation with Regularization and Moving Least Squares

Linear Approximation with Regularization and Moving Least Squares Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...

More information

Solutions to Problems Fundamentals of Condensed Matter Physics

Solutions to Problems Fundamentals of Condensed Matter Physics Solutons to Problems Fundamentals of Condensed Matter Physcs Marvn L. Cohen Unversty of Calforna, Berkeley Steven G. Loue Unversty of Calforna, Berkeley c Cambrdge Unversty Press 016 1 Acknowledgement

More information

Mechanics Physics 151

Mechanics Physics 151 Mechancs Physcs 5 Lecture 3 Contnuous Systems an Fels (Chapter 3) Where Are We Now? We ve fnshe all the essentals Fnal wll cover Lectures through Last two lectures: Classcal Fel Theory Start wth wave equatons

More information

CHARACTERISTICS OF COMPLEX SEPARATION SCHEMES AND AN ERROR OF SEPARATION PRODUCTS OUTPUT DETERMINATION

CHARACTERISTICS OF COMPLEX SEPARATION SCHEMES AND AN ERROR OF SEPARATION PRODUCTS OUTPUT DETERMINATION Górnctwo Geonżynera Rok 0 Zeszyt / 006 Igor Konstantnovch Mladetskj * Petr Ivanovch Plov * Ekaterna Nkolaevna Kobets * Tasya Igorevna Markova * CHARACTERISTICS OF COMPLEX SEPARATION SCHEMES AND AN ERROR

More information

Multi-electron atoms (11) 2010 update Extend the H-atom picture to more than 1 electron: H-atom sol'n use for N-elect., assume product wavefct.

Multi-electron atoms (11) 2010 update Extend the H-atom picture to more than 1 electron: H-atom sol'n use for N-elect., assume product wavefct. Mult-electron atoms (11) 2010 update Extend the H-atom pcture to more than 1 electron: VII 33 H-atom sol'n use for -elect., assume product wavefct. n ψ = φn l m where: ψ mult electron w/fct φ n l m one

More information

Mathematical Preparations

Mathematical Preparations 1 Introducton Mathematcal Preparatons The theory of relatvty was developed to explan experments whch studed the propagaton of electromagnetc radaton n movng coordnate systems. Wthn expermental error the

More information

Inexact Newton Methods for Inverse Eigenvalue Problems

Inexact Newton Methods for Inverse Eigenvalue Problems Inexact Newton Methods for Inverse Egenvalue Problems Zheng-jan Ba Abstract In ths paper, we survey some of the latest development n usng nexact Newton-lke methods for solvng nverse egenvalue problems.

More information

2. Hartree-Fock formalism

2. Hartree-Fock formalism 2. Hartree-Fock formalsm 2.1 Slater determnants One speaks of the Hartree-Fock (HF) approxmaton and of the HF method. The former s also called the self-consstent feld approxmaton, or mean feld approxmaton.

More information

Limited Dependent Variables

Limited Dependent Variables Lmted Dependent Varables. What f the left-hand sde varable s not a contnuous thng spread from mnus nfnty to plus nfnty? That s, gven a model = f (, β, ε, where a. s bounded below at zero, such as wages

More information

CHAPTER 14 GENERAL PERTURBATION THEORY

CHAPTER 14 GENERAL PERTURBATION THEORY CHAPTER 4 GENERAL PERTURBATION THEORY 4 Introducton A partcle n orbt around a pont mass or a sphercally symmetrc mass dstrbuton s movng n a gravtatonal potental of the form GM / r In ths potental t moves

More information

Rate of Absorption and Stimulated Emission

Rate of Absorption and Stimulated Emission MIT Department of Chemstry 5.74, Sprng 005: Introductory Quantum Mechancs II Instructor: Professor Andre Tokmakoff p. 81 Rate of Absorpton and Stmulated Emsson The rate of absorpton nduced by the feld

More information

8.323 Relativistic Quantum Field Theory I

8.323 Relativistic Quantum Field Theory I MI OpenCourseWare http://ocw.mt.edu 8.323 Relatvstc Quantum Feld heory I Sprng 2008 For nformaton about ctng these materals or our erms of Use, vst: http://ocw.mt.edu/terms. MASSACHUSES INSIUE OF ECHNOLOGY

More information

Lecture 10. Reading: Notes and Brennan Chapter 5

Lecture 10. Reading: Notes and Brennan Chapter 5 Lecture tatstcal Mechancs and Densty of tates Concepts Readng: otes and Brennan Chapter 5 Georga Tech C 645 - Dr. Alan Doolttle C 645 - Dr. Alan Doolttle Georga Tech How do electrons and holes populate

More information

Classical Field Theory

Classical Field Theory Classcal Feld Theory Before we embark on quantzng an nteractng theory, we wll take a dverson nto classcal feld theory and classcal perturbaton theory and see how far we can get. The reader s expected to

More information

P A = (P P + P )A = P (I P T (P P ))A = P (A P T (P P )A) Hence if we let E = P T (P P A), We have that

P A = (P P + P )A = P (I P T (P P ))A = P (A P T (P P )A) Hence if we let E = P T (P P A), We have that Backward Error Analyss for House holder Reectors We want to show that multplcaton by householder reectors s backward stable. In partcular we wsh to show fl(p A) = P (A) = P (A + E where P = I 2vv T s the

More information

Moments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is.

Moments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is. Moments of Inerta Suppose a body s movng on a crcular path wth constant speed Let s consder two quanttes: the body s angular momentum L about the center of the crcle, and ts knetc energy T How are these

More information

University of Washington Department of Chemistry Chemistry 453 Winter Quarter 2015

University of Washington Department of Chemistry Chemistry 453 Winter Quarter 2015 Lecture 2. 1/07/15-1/09/15 Unversty of Washngton Department of Chemstry Chemstry 453 Wnter Quarter 2015 We are not talkng about truth. We are talkng about somethng that seems lke truth. The truth we want

More information

1. Mean-Field Theory. 2. Bjerrum length

1. Mean-Field Theory. 2. Bjerrum length 1. Mean-Feld Theory Contnuum models lke the Posson-Nernst-Planck equatons are mean-feld approxmatons whch descrbe how dscrete ons are affected by the mean concentratons c and potental φ. Each on mgrates

More information

Explicit constructions of all separable two-qubits density matrices and related problems for three-qubits systems

Explicit constructions of all separable two-qubits density matrices and related problems for three-qubits systems Explct constructons of all separable two-qubts densty matrces and related problems for three-qubts systems Y. en-ryeh and. Mann Physcs Department, Technon-Israel Insttute of Technology, Hafa 2000, Israel

More information

Dynamics of a Superconducting Qubit Coupled to an LC Resonator

Dynamics of a Superconducting Qubit Coupled to an LC Resonator Dynamcs of a Superconductng Qubt Coupled to an LC Resonator Y Yang Abstract: We nvestgate the dynamcs of a current-based Josephson juncton quantum bt or qubt coupled to an LC resonator. The Hamltonan of

More information

Notes on Analytical Dynamics

Notes on Analytical Dynamics Notes on Analytcal Dynamcs Jan Peters & Mchael Mstry October 7, 004 Newtonan Mechancs Basc Asssumptons and Newtons Laws Lonely pontmasses wth postve mass Newtons st: Constant velocty v n an nertal frame

More information

Effects of Ignoring Correlations When Computing Sample Chi-Square. John W. Fowler February 26, 2012

Effects of Ignoring Correlations When Computing Sample Chi-Square. John W. Fowler February 26, 2012 Effects of Ignorng Correlatons When Computng Sample Ch-Square John W. Fowler February 6, 0 It can happen that ch-square must be computed for a sample whose elements are correlated to an unknown extent.

More information

OPTIMISATION. Introduction Single Variable Unconstrained Optimisation Multivariable Unconstrained Optimisation Linear Programming

OPTIMISATION. Introduction Single Variable Unconstrained Optimisation Multivariable Unconstrained Optimisation Linear Programming OPTIMIATION Introducton ngle Varable Unconstraned Optmsaton Multvarable Unconstraned Optmsaton Lnear Programmng Chapter Optmsaton /. Introducton In an engneerng analss, sometmes etremtes, ether mnmum or

More information

763622S ADVANCED QUANTUM MECHANICS Solution Set 1 Spring c n a n. c n 2 = 1.

763622S ADVANCED QUANTUM MECHANICS Solution Set 1 Spring c n a n. c n 2 = 1. 7636S ADVANCED QUANTUM MECHANICS Soluton Set 1 Sprng 013 1 Warm-up Show that the egenvalues of a Hermtan operator  are real and that the egenkets correspondng to dfferent egenvalues are orthogonal (b)

More information

STAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 16

STAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 16 STAT 39: MATHEMATICAL COMPUTATIONS I FALL 218 LECTURE 16 1 why teratve methods f we have a lnear system Ax = b where A s very, very large but s ether sparse or structured (eg, banded, Toepltz, banded plus

More information

DETAILED DERIVATION OF THE GENERAL MASTER EQUATION IN QUANTUM OPTICS

DETAILED DERIVATION OF THE GENERAL MASTER EQUATION IN QUANTUM OPTICS DETAILED DERIVATION OF THE GENERAL MASTER EQUATION IN QUANTUM OPTICS TIMOTHY JONES Abstract. An unorgnal but more detaled dervaton that typcal s presented of the master equaton. 1. Assumptons We begn by

More information

10. Canonical Transformations Michael Fowler

10. Canonical Transformations Michael Fowler 10. Canoncal Transformatons Mchael Fowler Pont Transformatons It s clear that Lagrange s equatons are correct for any reasonable choce of parameters labelng the system confguraton. Let s call our frst

More information

5.04, Principles of Inorganic Chemistry II MIT Department of Chemistry Lecture 32: Vibrational Spectroscopy and the IR

5.04, Principles of Inorganic Chemistry II MIT Department of Chemistry Lecture 32: Vibrational Spectroscopy and the IR 5.0, Prncples of Inorganc Chemstry II MIT Department of Chemstry Lecture 3: Vbratonal Spectroscopy and the IR Vbratonal spectroscopy s confned to the 00-5000 cm - spectral regon. The absorpton of a photon

More information

PHY688, Statistical Mechanics

PHY688, Statistical Mechanics Department of Physcs & Astronomy 449 ESS Bldg. Stony Brook Unversty January 31, 2017 Nuclear Astrophyscs James.Lattmer@Stonybrook.edu Thermodynamcs Internal Energy Densty and Frst Law: ε = E V = Ts P +

More information

LECTURE 9 CANONICAL CORRELATION ANALYSIS

LECTURE 9 CANONICAL CORRELATION ANALYSIS LECURE 9 CANONICAL CORRELAION ANALYSIS Introducton he concept of canoncal correlaton arses when we want to quantfy the assocatons between two sets of varables. For example, suppose that the frst set of

More information

Mechanics Physics 151

Mechanics Physics 151 Mechancs Physcs 151 Lecture 3 Lagrange s Equatons (Goldsten Chapter 1) Hamlton s Prncple (Chapter 2) What We Dd Last Tme! Dscussed mult-partcle systems! Internal and external forces! Laws of acton and

More information

Thermodynamics and statistical mechanics in materials modelling II

Thermodynamics and statistical mechanics in materials modelling II Course MP3 Lecture 8/11/006 (JAE) Course MP3 Lecture 8/11/006 Thermodynamcs and statstcal mechancs n materals modellng II A bref résumé of the physcal concepts used n materals modellng Dr James Ellott.1

More information

Solutions to Problem Set 6

Solutions to Problem Set 6 Solutons to Problem Set 6 Problem 6. (Resdue theory) a) Problem 4.7.7 Boas. n ths problem we wll solve ths ntegral: x sn x x + 4x + 5 dx: To solve ths usng the resdue theorem, we study ths complex ntegral:

More information

From Biot-Savart Law to Divergence of B (1)

From Biot-Savart Law to Divergence of B (1) From Bot-Savart Law to Dvergence of B (1) Let s prove that Bot-Savart gves us B (r ) = 0 for an arbtrary current densty. Frst take the dvergence of both sdes of Bot-Savart. The dervatve s wth respect to

More information

Integrals and Invariants of

Integrals and Invariants of Lecture 16 Integrals and Invarants of Euler Lagrange Equatons NPTEL Course Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng, Indan Insttute of Scence, Banagalore

More information

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2 Salmon: Lectures on partal dfferental equatons 5. Classfcaton of second-order equatons There are general methods for classfyng hgher-order partal dfferental equatons. One s very general (applyng even to

More information

Relaxation Methods for Iterative Solution to Linear Systems of Equations

Relaxation Methods for Iterative Solution to Linear Systems of Equations Relaxaton Methods for Iteratve Soluton to Lnear Systems of Equatons Gerald Recktenwald Portland State Unversty Mechancal Engneerng Department gerry@pdx.edu Overvew Techncal topcs Basc Concepts Statonary

More information

Lecture 6/7 (February 10/12, 2014) DIRAC EQUATION. The non-relativistic Schrödinger equation was obtained by noting that the Hamiltonian 2

Lecture 6/7 (February 10/12, 2014) DIRAC EQUATION. The non-relativistic Schrödinger equation was obtained by noting that the Hamiltonian 2 P470 Lecture 6/7 (February 10/1, 014) DIRAC EQUATION The non-relatvstc Schrödnger equaton was obtaned by notng that the Hamltonan H = P (1) m can be transformed nto an operator form wth the substtutons

More information

Lecture 16 Statistical Analysis in Biomaterials Research (Part II)

Lecture 16 Statistical Analysis in Biomaterials Research (Part II) 3.051J/0.340J 1 Lecture 16 Statstcal Analyss n Bomaterals Research (Part II) C. F Dstrbuton Allows comparson of varablty of behavor between populatons usng test of hypothess: σ x = σ x amed for Brtsh statstcan

More information

Mechanics Physics 151

Mechanics Physics 151 Mechancs Physcs 5 Lecture 7 Specal Relatvty (Chapter 7) What We Dd Last Tme Worked on relatvstc knematcs Essental tool for epermental physcs Basc technques are easy: Defne all 4 vectors Calculate c-o-m

More information

Physics 607 Exam 1. ( ) = 1, Γ( z +1) = zγ( z) x n e x2 dx = 1. e x2

Physics 607 Exam 1. ( ) = 1, Γ( z +1) = zγ( z) x n e x2 dx = 1. e x2 Physcs 607 Exam 1 Please be well-organzed, and show all sgnfcant steps clearly n all problems. You are graded on your wor, so please do not just wrte down answers wth no explanaton! Do all your wor on

More information

Tracking with Kalman Filter

Tracking with Kalman Filter Trackng wth Kalman Flter Scott T. Acton Vrgna Image and Vdeo Analyss (VIVA), Charles L. Brown Department of Electrcal and Computer Engneerng Department of Bomedcal Engneerng Unversty of Vrgna, Charlottesvlle,

More information

Lecture 5.8 Flux Vector Splitting

Lecture 5.8 Flux Vector Splitting Lecture 5.8 Flux Vector Splttng 1 Flux Vector Splttng The vector E n (5.7.) can be rewrtten as E = AU (5.8.1) (wth A as gven n (5.7.4) or (5.7.6) ) whenever, the equaton of state s of the separable form

More information

1 (1 + ( )) = 1 8 ( ) = (c) Carrying out the Taylor expansion, in this case, the series truncates at second order:

1 (1 + ( )) = 1 8 ( ) = (c) Carrying out the Taylor expansion, in this case, the series truncates at second order: 68A Solutons to Exercses March 05 (a) Usng a Taylor expanson, and notng that n 0 for all n >, ( + ) ( + ( ) + ) We can t nvert / because there s no Taylor expanson around 0 Lets try to calculate the nverse

More information

THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions

THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructons by George Hardgrove Chemstry Department St. Olaf College Northfeld, MN 55057 hardgrov@lars.acc.stolaf.edu Copyrght George

More information

ON MECHANICS WITH VARIABLE NONCOMMUTATIVITY

ON MECHANICS WITH VARIABLE NONCOMMUTATIVITY ON MECHANICS WITH VARIABLE NONCOMMUTATIVITY CIPRIAN ACATRINEI Natonal Insttute of Nuclear Physcs and Engneerng P.O. Box MG-6, 07725-Bucharest, Romana E-mal: acatrne@theory.npne.ro. Receved March 6, 2008

More information

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,

More information

This chapter illustrates the idea that all properties of the homogeneous electron gas (HEG) can be calculated from electron density.

This chapter illustrates the idea that all properties of the homogeneous electron gas (HEG) can be calculated from electron density. 1 Unform Electron Gas Ths chapter llustrates the dea that all propertes of the homogeneous electron gas (HEG) can be calculated from electron densty. Intutve Representaton of Densty Electron densty n s

More information

Lagrangian Theory. Several-body Systems

Lagrangian Theory. Several-body Systems Lagrangan Theory of Several-body Systems Ncholas Wheeler, Reed College Physcs Department November 995 Introducton. Let the N-tuple of 3-vectors {x (t) : =, 2,..., N} descrbe, relatve to an nertal frame,

More information

Supplementary information. Coulomb oscillations in a gate-controlled few-layer graphene quantum dot

Supplementary information. Coulomb oscillations in a gate-controlled few-layer graphene quantum dot Supplementary nformaton Coulomb oscllatons n a gate-controlled few-layer graphene quantum dot Ypu Song * Haonan Xong Wentao Jang Hongy Zhang Xao Xue Cheng Ma Yuln Ma Luyan Sun Hayan Wang and Lumng Duan

More information