Lecture 14: Forces and Stresses
|
|
- Marlene Lane
- 5 years ago
- Views:
Transcription
1 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
2 Overvew of Lecture Why bother? Theoretcal background CASTEP detals Symmetry and User Constrants Concluson Nuts and Bolts 2001 Lecture 14: Forces and Stresses 2
3 Why bother? (I) Structure optmsaton Mnmum energy corresponds to zero force Much more effcent than just usng energy alone Equlbrum bond lengths, angles, etc. Mnmum enthalpy corresponds to zero force and stress Can therefore mnmse enthalpy w.r.t. supercell shape due to nternal stress and external pressure Pressure-drven phase transtons Nuts and Bolts 2001 Lecture 14: Forces and Stresses 3
4 Why bother? (II) Molecular dynamcs Can do classcal dynamcs of ons usng forces derved from ab nto electronc structure Copes wth unusual geometry, bond-breakng, chemcal reactons, catalyss, dffuson, etc Incorporates effects of fnte temperature of ons Can generate thermodynamc nformaton from ensemble averagng Tme dependent phenomena Temperature drven phase transtons Nuts and Bolts 2001 Lecture 14: Forces and Stresses 4
5 Theoretcal Background Hellman-Feynman Theorem basc Quantum Mechancs Densty Functonal Theory how t apples n DFT Nuts and Bolts 2001 Lecture 14: Forces and Stresses 5
6 Hellman-Feynman Theorem (I) Classcally we have the force F at poston R s determned from the potental energy as F = U R ( R) Quantum mechancally we therefore expect F = R E where E = H Nuts and Bolts 2001 Lecture 14: Forces and Stresses 6
7 Hellman-Feynman Theorem (II) If we wrte the three unt cell vectors a, b, c as the columns of a matrx h then the effect of an appled stran s to change the shape of the unt cell: h = ( I + å)h We then have the stress tensor s related to the stran tensor e by: where Ù = σ αβ = 1 Ω ε αβ s the volume of the unt cell. a b c ( ) E Nuts and Bolts 2001 Lecture 14: Forces and Stresses 7
8 Stress and stran n acton c c b α b α γ β a c s xx γ β a a+δa s xy b γ γ+δγ α β a NB Much messer f non-orthogonal cell Nuts and Bolts 2001 Lecture 14: Forces and Stresses 8
9 Nuts and Bolts 2001 Lecture 14: Forces and Stresses 9 Hellman-Feynman Theorem (III) The Hellman-Feynman Theorem states that for any perturbaton λ we have whch obvously ncludes the case we are nterested n. We have assumed that the wavefuncton s properly normalsed and s an exact egenstate of H. = + = + + = + + = λ λ λ λ λ λ λ λ λ λ H H E H E H H H E
10 Hellman-Feynman Theorem (IV) To evaluate <E> for an unknown wavefuncton we frst expand t n terms of a complete set of fxed bass functons ϕ = c ϕ and then use the Varatonal Prncple to fnd the set of complex coeffcents c that mnmse the energy. If the bass set s ncomplete then we arrve at an upperbound for the energy. Nuts and Bolts 2001 Lecture 14: Forces and Stresses 10
11 Hellman-Feynman Theorem (V) If we have an approxmate egenstate, for example from usng an ncomplete bass set, then we must keep all 3 terms n the general expresson. If our bass set depends upon the onc postons, such as atomc centred Gaussans, then the other dervatves n the general expresson wll contrbute so-called Pulay forces (stresses). Note that Pulay forces (stresses) wll vansh n the lmt of a complete bass set, but that ths s never realzed n practce, or f poston ndependent bass functons, such as plane-waves, are used. Nuts and Bolts 2001 Lecture 14: Forces and Stresses 11
12 Nuts and Bolts 2001 Lecture 14: Forces and Stresses 12 Hellman-Feynman Theorem (VI) If we choose plane-waves as our bass functons, then because these functons are ndependent of the onc coordnates, t can easly seen that the general expresson for the forces becomes: ( ) ( ) ( ) ( ) = = j j j j j j r R H r c c r c H r c R R E, * * * * ϕ ϕ ϕ ϕ
13 Hellman-Feynman Theorem (VII) That s, we can calculate the forces usng the same expanson coeffcents as we used to varatonally mnmse the energy, usng matrx elements of the onc dervatve of the Hamltonan. Ths makes calculaton of the forces relatvely cheap once the varatonal energy mnmsaton has been completed f we are usng a plane-wave bass set. Smlar expressons can also be derved for the stresses. Nuts and Bolts 2001 Lecture 14: Forces and Stresses 13
14 Densty Functonal Theory (I) In DFT we have the Kohn-Sham Hamltonan: (, R ) = + V ( r) + V ( r, R ) + V ( r ) V ( R ) ˆ 2 r r e e on e XC on on H Therefore we only get contrbutons to the forces from the electron-on (pseudo)potental and the on-on Coulomb nteracton (the Ewald sum). Also contrbuton from exchange-correlaton potental f usng non-lnear core correctons. However, for the stresses, we also get a contrbuton from the knetc energy and Hartree terms. Nuts and Bolts 2001 Lecture 14: Forces and Stresses 14
15 Densty Functonal Theory (II) As we do not have a complete bass, the wavefuncton wll not be exact even wthn DFT. If we use a varatonal method to mnmze the total energy, then we know that the energy and hence the wavefuncton wll be correct to second order errors. However, the forces wll only be correct to frst order need a larger bass set for accurate forces than for energes. Nuts and Bolts 2001 Lecture 14: Forces and Stresses 15
16 Nuts and Bolts 2001 Lecture 14: Forces and Stresses 16
17 Densty Functonal Theory (III) However, f we use a non-varatonal mnmsaton technque, such as densty mxng, then such statements cannot be made. We can no longer guarantee that the energy found s an upper-bound on the true ground state energy. Ths complcates the applcaton of the Hellman- Feynman theorem. Consequently, non-varatonal forces and stresses are less relable. Nuts and Bolts 2001 Lecture 14: Forces and Stresses 17
18 CASTEP detals (I) Forces and stresses are almost the hghest level functonalty Sngle subroutne call to frst dervatves module s all that s requred to return the forces for the current model f ground state s already known. Dtto stress. frstd puts together the dfferent contrbutons to the force (stress) usng other functonal modules so the physcs s obvous. These n turn call down to operatons on charge denstes and potentals. Even at ths level the physcs s obvous and the detals of USPs etc are hdden. Nuts and Bolts 2001 Lecture 14: Forces and Stresses 18
19 CASTEP detals (II) The use of Ultra-Soft Pseudopotentals further complcates thngs, as there are now addtonal contrbutons to both the forces and the stresses from the charge augmentaton. However, the modular desgn of new CASTEP completely hdes ths from the hgher level programmer. Nuts and Bolts 2001 Lecture 14: Forces and Stresses 19
20 CASTEP detals (III) There s a problem wth the applcaton of the Hellman-Feynman theorem wth non-varatonal mnmsers Consequently the CASTEP code contans a frst-order correcton to the forces derved from densty mxng. However, the correspondng correcton to the stresses s not known. Ths has mplcatons for structure optmsaton and molecular dynamcs wth densty mxng Therefore the more recent Ensemble DFT approach (whch s fully varatonal) s to be preferred. Nuts and Bolts 2001 Lecture 14: Forces and Stresses 20
21 Nuts and Bolts 2001 Lecture 14: Forces and Stresses 21
22 CASTEP detals (IV) If the unt cell changes shape then the number of plane-waves and the FFT grd wll, at some pont, change dscontnuously. Consequently, t becomes dffcult to compare results at dfferent cell szes at the same nomnal bass set sze (cut-off energy) as the effectve qualty of the bass set s not the same, unless the bass set s fully converged (mpossble). Ths can be countered by usng the Fnte Bass Set Correcton, whch calculates the change n total energy upon changng the bass set sze at a fxed cell sze, and then uses ths to correct the total energy and stress at nearby cell szes. Nuts and Bolts 2001 Lecture 14: Forces and Stresses 22
23 Nuts and Bolts 2001 Lecture 14: Forces and Stresses 23
24 Symmetry If the symmetry of the system has been calculated (controlled by keyword n nput fle) then ths can be used to symmetrse the forces and stresses. Ths ensures that the forces (stresses) have the same symmetry as the model. Consequently, the symmetry of the system wll be preserved n any structural relaxaton or dynamcs. Nuts and Bolts 2001 Lecture 14: Forces and Stresses 24
25 User Constrants (I) Sometmes t s desred to mpose addtonal constrants on any structural relaxaton or dynamcs. Currently, CASTEP can apply any arbtrary number of lnear constrants on the atomc coordnates, up to the number of degrees of freedom. E.g. fxng an atom, constranng an atom to move n a lne or plane, fxng the relatve postons of pars of atoms, fxng the centre of mass of the system, etc. Nuts and Bolts 2001 Lecture 14: Forces and Stresses 25
26 User Constrants (II) Both force and stress constrants work n the same way by Lagrange multplers n the extended Lagrangan of the system: d dt L cons ( q, q& ) = L ( q q& ) λs ( q ) L q& cons = 0, L q 0 λ S m q&& = F λ q So for a gven set of constrants S(q) we need to know the dervatves of the constrants and then f we can determne λ we have the constrant force. S q Nuts and Bolts 2001 Lecture 14: Forces and Stresses 26
27 Nuts and Bolts 2001 Lecture 14: Forces and Stresses 27 User Constrants (III) Lnear constrants for the onc moton makes the matrx of constrant dervatves trval and therefore determnng λ: ( ) = = = a m a F m a m a m F a q a q S λ λ 0
28 User Constrants (IV) If we have multple constrants, S(q), R(q), etc. then we may satsfy each constrant smultaneously f the constrant matrces are all mutually orthogonal and so we use Gram-Schmdt on the coeffcents. Wth non-lnear constrants, both the constrants and ther dervatves need to be specfed, and t s not possble n general to satsfy them all smultaneously. Consequently, teratve procedures such as SHAKE or RATTLE are requred. Nuts and Bolts 2001 Lecture 14: Forces and Stresses 28
29 User Constrants (V) Constrants can also be appled to the unt cell lengths and angles (gvng 6 degrees of freedom). Any length (angle) can be held constant, or ted to one or both of the others. Ths s most useful f only a subset of the symmetres of the orgnal unt cell s to be enforced. Nuts and Bolts 2001 Lecture 14: Forces and Stresses 29
30 User Constrants (VI) Cell constrants are mplemented n a slghtly dfferent way. The stress tensor s transformed from ts normal symmetrc representaton n Cartesan coordnates nto the space of cell lengths and angles. There s then a 1:1 correspondence between the stress components and the cell degrees of freedom and so the constrants may be trvally appled to the stress and hence to the evoluton of the system. Nuts and Bolts 2001 Lecture 14: Forces and Stresses 30
31 Concluson Hellman-Feynman Theorem gves us a smple recpe for calculatng ab nto forces and stresses plane-wave bass has bg advantage DFT mplementaton Can be combned wth symmetry and/or constrants Major use n structural relaxaton and molecular dynamcs Nuts and Bolts 2001 Lecture 14: Forces and Stresses 31
Army Ants Tunneling for Classical Simulations
Electronc Supplementary Materal (ESI) for Chemcal Scence. Ths journal s The Royal Socety of Chemstry 2014 electronc supplementary nformaton (ESI) for Chemcal Scence Army Ants Tunnelng for Classcal Smulatons
More informationRobert 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 informationPhysics 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 informationCanonical transformations
Canoncal transformatons November 23, 2014 Recall that we have defned a symplectc transformaton to be any lnear transformaton M A B leavng the symplectc form nvarant, Ω AB M A CM B DΩ CD Coordnate transformatons,
More informationPHYS 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 informationPhysics 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 informationThe 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 informationNUMERICAL DIFFERENTIATION
NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the
More informationTHEOREMS 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 informationA 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 informationIndeterminate pin-jointed frames (trusses)
Indetermnate pn-jonted frames (trusses) Calculaton of member forces usng force method I. Statcal determnacy. The degree of freedom of any truss can be derved as: w= k d a =, where k s the number of all
More informationwhere 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
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.
More information10. 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 informationAdvanced 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 informationPhysics 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 informationLecture 10 Support Vector Machines II
Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed
More informationCSci 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 informationNon-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 information5. Response properties in ab initio schemes
5. Response propertes n ab nto schemes A number of mportant physcal observables s expressed va dervatves of total energy (or free energy) E. Examples are: E R 2 E R a R b forces on the nucle; crtcal ponts
More informationMechanics 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 informationPolynomial Regression Models
LINEAR REGRESSION ANALYSIS MODULE XII Lecture - 6 Polynomal Regresson Models Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur Test of sgnfcance To test the sgnfcance
More information1 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations
Physcs 171/271 -Davd Klenfeld - Fall 2005 (revsed Wnter 2011) 1 Dervaton of Rate Equatons from Sngle-Cell Conductance (Hodgkn-Huxley-lke) Equatons We consder a network of many neurons, each of whch obeys
More informationCHAPTER 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 informationLagrange 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 informationDensity 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 information5.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 informationProbabilistic 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 informationPhysics 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 informationThermodynamics 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 informationSome 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 informationCausal Diamonds. M. Aghili, L. Bombelli, B. Pilgrim
Causal Damonds M. Aghl, L. Bombell, B. Plgrm Introducton The correcton to volume of a causal nterval due to curvature of spacetme has been done by Myrhem [] and recently by Gbbons & Solodukhn [] and later
More information763622S 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 informationImplicit Integration Henyey Method
Implct Integraton Henyey Method In realstc stellar evoluton codes nstead of a drect ntegraton usng for example the Runge-Kutta method one employs an teratve mplct technque. Ths s because the structure
More informationMechanics Physics 151
Mechancs Physcs 5 Lecture 0 Canoncal Transformatons (Chapter 9) What We Dd Last Tme Hamlton s Prncple n the Hamltonan formalsm Dervaton was smple δi δ p H(, p, t) = 0 Adonal end-pont constrants δ t ( )
More informationELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM
ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look
More informationψ = i c i u i c i a i b i u i = i b 0 0 b 0 0
Quantum Mechancs, Advanced Course FMFN/FYSN7 Solutons Sheet Soluton. Lets denote the two operators by  and ˆB, the set of egenstates by { u }, and the egenvalues as  u = a u and ˆB u = b u. Snce the
More informationcoordinates. Then, the position vectors are described by
Revewng, what we have dscussed so far: Generalzed coordnates Any number of varables (say, n) suffcent to specfy the confguraton of the system at each nstant to tme (need not be the mnmum number). In general,
More informationEcon107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)
I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes
More informationIntroduction 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 informationEEE 241: Linear Systems
EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they
More informationMechanics 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 informationSIO 224. m(r) =(ρ(r),k s (r),µ(r))
SIO 224 1. A bref look at resoluton analyss Here s some background for the Masters and Gubbns resoluton paper. Global Earth models are usually found teratvely by assumng a startng model and fndng small
More informationLINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity
LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 30 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 2 Remedes for multcollnearty Varous technques have
More informationStatistical mechanics handout 4
Statstcal mechancs handout 4 Explan dfference between phase space and an. Ensembles As dscussed n handout three atoms n any physcal system can adopt any one of a large number of mcorstates. For a quantum
More informationSupplementary Notes for Chapter 9 Mixture Thermodynamics
Supplementary Notes for Chapter 9 Mxture Thermodynamcs Key ponts Nne major topcs of Chapter 9 are revewed below: 1. Notaton and operatonal equatons for mxtures 2. PVTN EOSs for mxtures 3. General effects
More informationTHE 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 informationGrover s Algorithm + Quantum Zeno Effect + Vaidman
Grover s Algorthm + Quantum Zeno Effect + Vadman CS 294-2 Bomb 10/12/04 Fall 2004 Lecture 11 Grover s algorthm Recall that Grover s algorthm for searchng over a space of sze wors as follows: consder the
More informationLagrangian 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 informationOpen Systems: Chemical Potential and Partial Molar Quantities Chemical Potential
Open Systems: Chemcal Potental and Partal Molar Quanttes Chemcal Potental For closed systems, we have derved the followng relatonshps: du = TdS pdv dh = TdS + Vdp da = SdT pdv dg = VdP SdT For open systems,
More informationPhysics 106a, Caltech 11 October, Lecture 4: Constraints, Virtual Work, etc. Constraints
Physcs 106a, Caltech 11 October, 2018 Lecture 4: Constrants, Vrtual Work, etc. Many, f not all, dynamcal problems we want to solve are constraned: not all of the possble 3 coordnates for M partcles (or
More informationNegative Binomial Regression
STATGRAPHICS Rev. 9/16/2013 Negatve Bnomal Regresson Summary... 1 Data Input... 3 Statstcal Model... 3 Analyss Summary... 4 Analyss Optons... 7 Plot of Ftted Model... 8 Observed Versus Predcted... 10 Predctons...
More informationQuantum Field Theory III
Quantum Feld Theory III Prof. Erck Wenberg February 16, 011 1 Lecture 9 Last tme we showed that f we just look at weak nteractons and currents, strong nteracton has very good SU() SU() chral symmetry,
More informationCHAPTER 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 informationSnce 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 information5.60 Thermodynamics & Kinetics Spring 2008
MIT OpenCourseWare http://ocw.mt.edu 5.60 Thermodynamcs & Knetcs Sprng 2008 For nformaton about ctng these materals or our Terms of Use, vst: http://ocw.mt.edu/terms. 5.60 Sprng 2008 Lecture #29 page 1
More information7. Products and matrix elements
7. Products and matrx elements 1 7. Products and matrx elements Based on the propertes of group representatons, a number of useful results can be derved. Consder a vector space V wth an nner product ψ
More informationKernel Methods and SVMs Extension
Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general
More information12. 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 informationNumerical 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 informationWorkshop: 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
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 informationLesson 5: Kinematics and Dynamics of Particles
Lesson 5: Knematcs and Dynamcs of Partcles hs set of notes descrbes the basc methodology for formulatng the knematc and knetc equatons for multbody dynamcs. In order to concentrate on the methodology and
More informationLecture 4. Macrostates and Microstates (Ch. 2 )
Lecture 4. Macrostates and Mcrostates (Ch. ) The past three lectures: we have learned about thermal energy, how t s stored at the mcroscopc level, and how t can be transferred from one system to another.
More informationReport on Image warping
Report on Image warpng Xuan Ne, Dec. 20, 2004 Ths document summarzed the algorthms of our mage warpng soluton for further study, and there s a detaled descrpton about the mplementaton of these algorthms.
More informationCHAPTER 6. LAGRANGE S EQUATIONS (Analytical Mechanics)
CHAPTER 6 LAGRANGE S EQUATIONS (Analytcal Mechancs) 1 Ex. 1: Consder a partcle movng on a fxed horzontal surface. r P Let, be the poston and F be the total force on the partcle. The FBD s: -mgk F 1 x O
More informationIntegrals 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 informationProf. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model
EXACT OE-DIMESIOAL ISIG MODEL The one-dmensonal Isng model conssts of a chan of spns, each spn nteractng only wth ts two nearest neghbors. The smple Isng problem n one dmenson can be solved drectly n several
More informationLinear Feature Engineering 11
Lnear Feature Engneerng 11 2 Least-Squares 2.1 Smple least-squares Consder the followng dataset. We have a bunch of nputs x and correspondng outputs y. The partcular values n ths dataset are x y 0.23 0.19
More informationIn this section is given an overview of the common elasticity models.
Secton 4.1 4.1 Elastc Solds In ths secton s gven an overvew of the common elastcty models. 4.1.1 The Lnear Elastc Sold The classcal Lnear Elastc model, or Hooean model, has the followng lnear relatonshp
More informationElectronic Structure for Excited States (multiconfigurational methods) Spiridoula Matsika
Electronc Structure for Excted States (multconfguratonal methods) Sprdoula Matska Excted Electronc States Theoretcal treatment of excted states s needed for: UV/Vs electronc spectroscopy Photochemstry
More informationLAGRANGIAN 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 informationMoments 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 informationProgramming Project 1: Molecular Geometry and Rotational Constants
Programmng Project 1: Molecular Geometry and Rotatonal Constants Center for Computatonal Chemstry Unversty of Georga Athens, Georga 30602 Summer 2012 1 Introducton Ths programmng project s desgned to provde
More informationSupporting Information
Supportng Informaton The neural network f n Eq. 1 s gven by: f x l = ReLU W atom x l + b atom, 2 where ReLU s the element-wse rectfed lnear unt, 21.e., ReLUx = max0, x, W atom R d d s the weght matrx to
More informationTracking 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 informationn α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0
MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector
More informationLecture 2: Prelude to the big shrink
Lecture 2: Prelude to the bg shrnk Last tme A slght detour wth vsualzaton tools (hey, t was the frst day... why not start out wth somethng pretty to look at?) Then, we consdered a smple 120a-style regresson
More informationMathematical 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 informationDUE: WEDS FEB 21ST 2018
HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant
More informationPHYS 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 informationOn balancing multiple video streams with distributed QoS control in mobile communications
On balancng multple vdeo streams wth dstrbuted QoS control n moble communcatons Arjen van der Schaaf, José Angel Lso Arellano, and R. (Inald) L. Lagendjk TU Delft, Mekelweg 4, 68 CD Delft, The Netherlands
More informationNP-Completeness : Proofs
NP-Completeness : Proofs Proof Methods A method to show a decson problem Π NP-complete s as follows. (1) Show Π NP. (2) Choose an NP-complete problem Π. (3) Show Π Π. A method to show an optmzaton problem
More informationIntroduction 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 informationLecture 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 informationSupplemental document
Electronc Supplementary Materal (ESI) for Physcal Chemstry Chemcal Physcs. Ths journal s the Owner Socetes 01 Supplemental document Behnam Nkoobakht School of Chemstry, The Unversty of Sydney, Sydney,
More informationQuantum Mechanics for Scientists and Engineers. David Miller
Quantum Mechancs for Scentsts and Engneers Davd Mller Types of lnear operators Types of lnear operators Blnear expanson of operators Blnear expanson of lnear operators We know that we can expand functons
More informationUncertainty in measurements of power and energy on power networks
Uncertanty n measurements of power and energy on power networks E. Manov, N. Kolev Department of Measurement and Instrumentaton, Techncal Unversty Sofa, bul. Klment Ohrdsk No8, bl., 000 Sofa, Bulgara Tel./fax:
More informationand Statistical Mechanics Material Properties
Statstcal Mechancs and Materal Propertes By Kuno TAKAHASHI Tokyo Insttute of Technology, Tokyo 15-855, JAPA Phone/Fax +81-3-5734-3915 takahak@de.ttech.ac.jp http://www.de.ttech.ac.jp/~kt-lab/ Only for
More informationRate 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 informationIntroduction to Statistical Methods
Introducton to Statstcal Methods Physcs 4362, Lecture #3 hermodynamcs Classcal Statstcal Knetc heory Classcal hermodynamcs Macroscopc approach General propertes of the system Macroscopc varables 1 hermodynamc
More informationPhysics 53. Rotational Motion 3. Sir, I have found you an argument, but I am not obliged to find you an understanding.
Physcs 53 Rotatonal Moton 3 Sr, I have found you an argument, but I am not oblged to fnd you an understandng. Samuel Johnson Angular momentum Wth respect to rotatonal moton of a body, moment of nerta plays
More informationModule 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur
Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:
More informationarxiv:cond-mat/ v2 [cond-mat.mes-hall] 3 Jan 2006
arxv:cond-mat/0210519v2 [cond-mat.mes-hall] 3 Jan 2006 Non Equlbrum Green s Functons for Dummes: Introducton to the One Partcle NEGF equatons Magnus Paulsson Dept. of mcro- and nano-technology, NanoDTU,
More informationFeb 14: Spatial analysis of data fields
Feb 4: Spatal analyss of data felds Mappng rregularly sampled data onto a regular grd Many analyss technques for geophyscal data requre the data be located at regular ntervals n space and/or tme. hs s
More informationEEL 6266 Power System Operation and Control. Chapter 3 Economic Dispatch Using Dynamic Programming
EEL 6266 Power System Operaton and Control Chapter 3 Economc Dspatch Usng Dynamc Programmng Pecewse Lnear Cost Functons Common practce many utltes prefer to represent ther generator cost functons as sngle-
More informationTransfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system
Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng
More informationKinematics of Fluids. Lecture 16. (Refer the text book CONTINUUM MECHANICS by GEORGE E. MASE, Schaum s Outlines) 17/02/2017
17/0/017 Lecture 16 (Refer the text boo CONTINUUM MECHANICS by GEORGE E. MASE, Schaum s Outlnes) Knematcs of Fluds Last class, we started dscussng about the nematcs of fluds. Recall the Lagrangan and Euleran
More informationClassical 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 informationWeek 9 Chapter 10 Section 1-5
Week 9 Chapter 10 Secton 1-5 Rotaton Rgd Object A rgd object s one that s nondeformable The relatve locatons of all partcles makng up the object reman constant All real objects are deformable to some extent,
More informationLecture 3. Ax x i a i. i i
18.409 The Behavor of Algorthms n Practce 2/14/2 Lecturer: Dan Spelman Lecture 3 Scrbe: Arvnd Sankar 1 Largest sngular value In order to bound the condton number, we need an upper bound on the largest
More information