Inseting this into the left hand side of the equation of motion above gives the most commonly used algoithm in classical molecula dynamics simulations
|
|
- Irma Nash
- 5 years ago
- Views:
Transcription
1 Chem465 in 2000 Univesity of Washington Lectue notes Hannes Jonsson Classical dynamics When we ae dealing with heavy atoms and high enough enegy o tempeatue, it is often suciently accuate to neglect quantum eects and calculate the time evolution of a system (such as an atom o a collection of atoms) using the classical equation of motion, Newton's second law, F ~ = d ~P dt. Hee F ~ is the foce acting on the classical paticle and P ~ is the momentum given by P ~ = m~v. Since the velocity is the st deivative of position with espect to time, the equation of motion is a second ode dieential equation fo the position of the paticle as a function of time (the second deivative being the acceleation) d 2 x(t) dt 2 = F=m : In most cases of inteest in chemisty, the foce is consevative, i.e. it can be witten as the dieential of a scala function, the potential enegy, V (x). In one dimensional systems F =,dv=dx. The foce is then only a function of the coodinate of the paticle, F (x), and not, fo example, dependent on the velocity. The poblem of calculating the time evolution of such a system, theefoe, educes to the poblem of nding a good epesentation of the potential enegy as a function of the coodinates of the atoms and then solving the second ode dieential equation fo the classical tajectoies of the atoms. We will st focus on the second aspect of this poblem. We will use a nite dieence appoximation to calculate numeically the tajectoy fo a given initial condition and a given potential enegy function. The nite dieence algoithm is called the Velet algoithm. Then, we will examine a few potential functions which ae commonly used to epesent appoximately the inteaction between atoms. Simple deivation of the Velet algoithm Hee, time will be discetized, i.e. time evolves in discete steps t = h, but the spacial coodinates will not be discetized. To begin with, we will focus on a one-dimensional poblem fo simplicity, and then genealize the algoithm to thee-dimensional space. Let x k be the position afte k timesteps, i.e. x k = x(t = kh). The cental nite dieence appoximation fo the second deivative is x k = x k+1 + x k,1, 2x k h 2 : 16
2 Inseting this into the left hand side of the equation of motion above gives the most commonly used algoithm in classical molecula dynamics simulations, the Velet algoithm x k+1 = 2x k, x k,1 + h 2 F (x k )=m : Since the dieential equation is second ode, we need to know two points to get the ecusion stated. Deivation using Taylo's expansions The Velet algoithm can also be deived using Taylo expansions. Then it becomes cleae what the accuacy of the algoithm is, in paticula how the eo scales with the size of the timestep. Taylo expanding x(t + h) about x(t) gives x(t + h) = x(t) + h _x(t) + h2 2 x(t) + h3 3 x(3) + O(h 4 ) : Expanding x(t, h) gives a simila expession, except the odd tems have a minus sign x(t, h) = x(t), h _x(t) + h2 2 When the two expessions ae added, the odd tems cancel x(t), h3 3 x(3) + O(h 4 ) : x(t + h) + x(t, h) = 2x(t) + h 2 x(t) + O(h 4 ) : This expession is geneal, valid fo any function x(t) that is dieentiable enough times. The equation of motion fo classical dynamics can now be used to eliminate the second deivative x(t + h) + x(t, h) = 2x(t) + h 2 F (x(t))=m + O(h 4 ) : The Velet algoithm is obtained by neglecting the fouth and highe ode tems. It tuns out that this algoithm is `stable', i.e. the eo made at any given step tends to decay athe than magnify late on. Thee is no guaantee fo such behaviou in nite dieence algoithms. Some ae stable and some ae not. The easiest way to nd out whethe an algoithm is stable is to ty it and see if it woks. Deivation of the Velocity Velet algoithm In classical mechanics it is possible to know simultaneously both the position and velocity of paticles (paticles being atoms, fo example). A classical system is often epesented as a point in `phase space', a space of all the coodinates and momenta of paticles in the system. Knowing the position and velocity (o momentum) at one point 17
3 in time completely species the state of the system and, though the equation of motion, the futue (and past) motion of the paticles, as long as the foce acting on the paticles ae known. The velocity, howeve, does not appea explicitly in the Velet algoithm. If it is needed, fo example, to calculate the kinetic enegy o to implement some kind of tempeatue contol (modications of the kinetic enegy), a cental nite dieence estimate can be used, _x k =(x k+1, x k,1 )=2h. It is a bit awkwad that _x k cannot be found until x k+1 has been calculated. In stating up a calculation, we typically know whee the paticle is and we have some idea of what velocity it has (fo example by knowing the kinetic enegy), but this is not enough to stat the Velet algoithm because the positions at two adjacent timesteps ae needed. Given the coodinate x(0) and velocity _x(0), the coodinate at a pevious step x(,h) needs to be constucted. A dieent algoithm, the Velocity Velet algoithm, explicitly includes the velocity at each step and is `self-stating' fom the position and velocity at the initial time. It is mathematically identical to the oiginal Velet algoithm in the sense that it geneates the same tajectoy in the absence of oundo eos in the compute. Any n,th ode dieential equation can be educed to a set of n st ode dieential equations. In paticula, the classical equation of motion x(t) = F (x)=m is second ode and can be educed to two st ode equations. Let v(t) = _x(t) = dx=dt and the two equations ae _x(t) = v(t) (1) and _v(t) = F (x(t))=m (2): A nite dieence algoithm can again be deived using Taylo expansions. Stating with x(t + h) x(t + h) = x(t) + h _x(t) + h2 2 x(t) + O(h3 ) : Using v(t) to eliminate _x and F=m to eliminate x gives x(t + h) = x(t) + hv(t) + h2 2 Then, expanding the second function, v(t + h) F (x(t)) m + O(h 3 ) (1b): v(t + h) = v(t) + h _v(t) + h2 2 v(t) + O(h3 ) : We can use F=m to eliminate _v, but we need to develop an expession fo v in tems of known quantities. This can be done by expanding _v(t + h) _v(t + h) = _v(t) + hv(t) + O(h 2 ) : 18
4 It is enough to go up to ode h 2 hee because we only need an appoximation that is good to ode h 3 to the quantity h2 v(t). Multiplying by h=2 and eaanging, gives 2 so, the expession fo v(t + h) becomes h 2 2 v(t) = h 2 (_v(t + h), _v(t)) + O(h3 ) v(t + h) = v(t) + h _v(t) + h 2 (_v(t + h), _v(t)) + O(h3 ) : Using the equation of motion, this can nally be ewitten as v(t + h) = v(t) + h 2m (F (x(t + h)) + F (x(t))) + O(h3 ) (2b): Schematically, the velocity Velet algoithm is as follows Given x k and v k and an expession fo F (x) Step 1: calculate x k+1 = x k + hv k + h 2 F (x k )=2m Step 2: evaluate F (x k+1 ) Step 3: calculate v k+1 = v k + h 2m (F (x k)+f (x k+1 )) Now all quantities fo the new step, k +1, have been found, go back to step 1. Genealization to many atoms in thee dimensions: The Velet algoithm can easily be genealized to highe dimensions and many atoms. Fo N atoms in thee dimensional space, the potential function depends on 3N coodinates, V (x1;y1;z1;x2;y2;z2;x3;:::;xn;yn;zn) = V (~ 1 ;~ 2 ;:::;~ N ). The foce on atom numbe i is ~F i =, i V (~ 1 ;~ 2 ;:::;~ i ;:::;~ N ) whee the subscipt i on the gadient opeato means that the deivatives should be taken with espect to coodinates of atom i. Witing each catesian component sepaately Fxi =, d dxi V (x1;y1;z1;x2;y2;:::;xi;:::) 19
5 Fyi =, d dyi V (x1;y1;z1;x2;y2;:::;xi;:::) Fzi =, d dzi V (x1;y1;z1;x2;y2;:::;xi;:::) : A ecusion elation given by the Velet algoithm o the velocity Velet algoithm can be used fo each [coodinate,velocity] pai, but they all get coupled togethe though the foce, which geneally depends on coodinates of all the atoms. Let x1 k be the x coodinate of atom 1 at timestep k and y1 k be the y coodinate of atom 1 at timestep k, etc. Futhemoe, let v x denote the x component of the velocity and v y denote the y component and v z the z component. The foce components ae Fxi k Fxi(x1 k ;y1 k ;:::;yn k ): Then a schematic algoithm fo N atoms in thee dimensions can be witten as follows: Given x1 k ;y1 k ;:::;yn k and u1 k ;v1 k ;:::;vn k and an expession fo Fx1 k ;Fy1 k ;:::;FyN k Step 1: calculate x1 k+1 = x1 k + hv x 1 k + h 2 Fx1 k =2m y1 k+1 = y1 k + hv y 1 k + h 2 Fy1 k =2m z1 k+1 = z1 k + hv z 1 k + h 2 Fz1 k =2m ::: zn k+1 = zn k + hv z N k + h 2 FzN k =2m Step 2: evaluate Fx1 k+1 ;Fy1 k+1 :::;FzN k+1 Step 3: calculate v x 1 k+1 = v x 1 k + h 2m (Fx1 k + Fx1 k+1 ) v y 1 k+1 = v y 1 k + h 2m (Fy1 k + Fy1 k+1 ) ::: v z N k+1 = v z N k + h 2m (FzN k + FzN k+1 ) Now all quantities fo the new step, k +1, have been found, go back to step 1. Potential enegy functions In calculating the dynamics of atoms and molecules, it is vey impotant to choose a potential function that accuately mimics the system of inteest. The potential function descibes how the potential enegy of a system of atoms depends on the coodinates of the 20
6 atoms. It is assumed hee that the electons adjust to new atomic positions much faste than the motion of the atomic nuclei (this is called the Bon-Oppenheime appoximation) and the potential function that is needed fo the nuclea motion is stictly the enegy obtained afte calculating the electonic wavefunction keeping the atomic coodinates xed. The calculation of the electonic wavefunction fo a system of many atoms is vey dif- cult. Most often, a simple functional fom is assumed fo the potential function and the paametes adjusted to epoduce some expeimental o theoetical data. Those ae called empiical o semi-empiical potential functions. Although we use the oppotunity hee to biey discuss potential functions in the context of classical simulations, the same consideations and functions apply when the motion of the atoms and molecules is teated quantum mechanically. In pinciple, the inteaction potential of N atoms can be expanded in a many-body expansion V ( 1 ; 2 ;:::; N ) = X i v i (~ i )+ X i X j>i v 2 (~ i ;~ j )+ X i X j>i X k>j>i v 3 (~ i ;~ j ;~ k )+ ::: whee the st sum is ove one-body tems, the second sum ove paiwise inteactions, the thid one ove thee-body contibutions, etc. The one-body tems aise if an extenal eld is applied (fo example, the inteaction with an extenal electical eld, o the potential descibing the wall of a containe). The paiwise inteactions ae usually most impotant and ae at shot ange epulsive due to epulsion of the two electon clouds and at long ange attactive due to the induced-dipole/induced-dipole inteactions. The thee-body tems aise because the inteaction of a pai of atoms is modied by the pesence of a thid atom. Fo ae gases, the paiwise potentials alone descibe quite well the potential enegy function, but even thee the thee-body coections ae signicant (amounting to ca. 10% in the binding enegy of the cystals of heavie ae gases). Fo systems with stong, bonding inteactions, the many-body expansion is a vey poo appoach, because the expansion conveges vey slowly. The most commonly used pai potential is the Lennad-Jones potential v() = , whee the st tem descibes the epulsion at shot ange as the electon clouds ovelap moe than is optimal and the second tem epesents the long ange attaction. The fom of the second tem,6 has the coect behaviou fo the induced-dipole/induced-dipole inteaction. The fom of the epulsive tem is, howeve, not suppoted by theoetical calculations and is simply chosen fo convenience. A moe accuate pai potential fom, fo example, fo descibing the inteaction of ae gas atoms is v() = Ae,, f () C 6 : 21
7 The ovelap of the closed shell electon clouds is usually well descibed with an exponential fom. The function f () is a switching function, which is 1 at long ange but goes to zeo at shot ange and pevents the attactive enegy tem fom diveging. Bonded inteactions, whee the attactive inteaction comes fom the fomation of a chemical bond, is bette descibed with a Mose potential. v() = D e,2(,0), 2e,(, 0) : The attaction decays exponentially as the ovelap of the electon clouds deceases. In addition to this exponential attaction, a longe ange van de Waals attaction should be added, but in compaison to chemical bond enegies, this can often be neglected. When bonding constaints need to be taken into account, fo example, when the numbe of bonded neighbos is small and the angle between bonds has a pefeed value, the many-body tems ae vey impotant and the potential function becomes quite complicated. One good example is cabon, which only bonds to at most 4 othe atoms if sp 3 hybidized (as opposed to ae gas atoms that can suound themselves with up to 12 nea neighbos), but in dieent hybidization states will pefe 2 o 3 neighbos. A geat deal of eot in theoetical chemisty goes into the development of accuate and yet convenient functional foms that take into account many-body eects. One functional fom that has been used extensively in modeling chemical eactions involving fomation and beaking of chemical bonds is the LEPS fom. As an example, the `extended LEPS' potential fo a thee atom system A; B and C, capable of descibing a substitution eaction A + BC! AB + C is of the fom E( ab ; bc ; ac ) = Q ab 1+a + Q bc 1+b + Q ac 1+c, J 2 ab (1 + a) 2 + J 2 bc (1 + b) 2 + J 2 ac (1 + c) 2, J ab J bc (1 + a)(1 + b), J bc J ac (1 + b)(1 + c), J ab J ac (1 + a)(1 + c) The paametes can be chosen in such away that the enegy is lage when all thee atoms ae close to each othe, i.e. when only one bond is allowed at a time. The paametes a; b and c can be adjusted to change the height and location of the potential baie fo the eaction. The functional fom is inspied by appoximate calculations of the inteaction of thee H atoms. The Qs stand fo Coulomb inteactions between the electon clouds and the nuclei and the Js stand fo exchange inteactions (tems that aise because of the indistinguishability of the quantum mechanical electons). A typical fom fo the distance dependence of the Q and J tems is simila to the Mose potential Q() = D e,2(, 0), e,(, 0)
8 and J() = D 4 e,2(,0), 6e,(, 0) : The bonded inteaction of two atoms is Q+J (in the singlet state, whee unpaied electons on each atom have opposite spin and can pai up as the atoms appoach each othe) and the non-bonded inteaction is Q, J (puely epulsive aswould be the case if the unpaied electons have the same spin, a tiplet). Deivation of the foce fo L-J potentials: In calculating the classical dynamics of atoms it is necessay to be able to evaluate the foce acting on the atoms. Given a potential enegy suface, the task is to calculate the gadient ~F i =, ~ i V (~ 1 ;~ 2 ;:::;~ i ;:::) with espect to the coodinates of each atom. This equation says that the foce acting on atom i is obtained by taking the potential function, which depends on the coodinates of all the atoms in the system, and dieentiate it with espect to the coodinates of atom i. This can be a bit involved when the potential function is complicated. Fotunately, Mathematica and simila pogams can be of geat help. As a simple example, a deivation of the foce acting on atoms inteacting via a paiwise additive potential is given below. The total potential enegy of the system is a sum ove a potential, v(), fo each distinct pai of atoms V (~ 1 ;~ 2 ;:::;~ N ) = NX i NX j<i v( ij ) whee ij is the distance between atoms i and j. 23
9 The distance can be obtained fom the coodinates of the atoms (which efe to some abitay oigin of the coodinate system) q q ij = (~ i, ~ j ) (~ i, ~ j ) = x 2 ij + y2 ij + z2 ij whee x ij = x i, x j. The foce on atom k is a vecto pointing in the diection of the steepest descent of the potential enegy. The expession is ~F k =, ~ k V (~ 1 ;~ 2 ;:::;~ k ;:::) =, X j6=k =, X k ~ k v( kj ) k k v( kj kj Using the chain ule fo dieentiation, we can see that k = 1. Theefoe, the deivative of the pai potential above can be witten in tems of v 0, i.e. the deivative with espect to the agument of the function. Fo the x component of the foce Fo all thee kj kj = v 0 ( kj kj ~ k v( kj ) = v0 ( kj ) kj = x kj kj v 0 ( kj ) : (x kj ^x + y kj ^y + z kj ^z) = (~ k, ~ j ) kj v 0 ( kj ) = ^ kj v 0 ( kj ) whee ^ kj is a unit vecto in the diection fom atom j to k. we get In the paticulaly simple case whee the pai poential is a Lennad-Jones potential The foce on atom k becomes v() = 4 v 0 () = 4 ~F k =, X j6=k 12 6,, (~ k, ~ j ) kj kj : : 24
Nuclear and Particle Physics - Lecture 20 The shell model
1 Intoduction Nuclea and Paticle Physics - Lectue 0 The shell model It is appaent that the semi-empiical mass fomula does a good job of descibing tends but not the non-smooth behaviou of the binding enegy.
More informationProblem Set 10 Solutions
Chemisty 6 D. Jean M. Standad Poblem Set 0 Solutions. Give the explicit fom of the Hamiltonian opeato (in atomic units) fo the lithium atom. You expession should not include any summations (expand them
More information7.2. Coulomb s Law. The Electric Force
Coulomb s aw Recall that chaged objects attact some objects and epel othes at a distance, without making any contact with those objects Electic foce,, o the foce acting between two chaged objects, is somewhat
More information(read nabla or del) is defined by, k. (9.7.1*)
9.7 Gadient of a scala field. Diectional deivative Some of the vecto fields in applications can be obtained fom scala fields. This is vey advantageous because scala fields can be handled moe easily. The
More informationQuestion 1: The dipole
Septembe, 08 Conell Univesity, Depatment of Physics PHYS 337, Advance E&M, HW #, due: 9/5/08, :5 AM Question : The dipole Conside a system as discussed in class and shown in Fig.. in Heald & Maion.. Wite
More informationASTR415: Problem Set #6
ASTR45: Poblem Set #6 Cuan D. Muhlbege Univesity of Mayland (Dated: May 7, 27) Using existing implementations of the leapfog and Runge-Kutta methods fo solving coupled odinay diffeential equations, seveal
More informationInteratomic Forces. Overview
Inteatomic Foces Oveview an de Walls (shot ange ~1/ 6, weak ~0.010.1 e) Ionic (long ange, ~1/, stong ~510 e) Metallic (no simple dependence, ~0.1e) Covalent (no simple dependence, diectional,~3 e) Hydogen
More information5.61 Physical Chemistry Lecture #23 page 1 MANY ELECTRON ATOMS
5.6 Physical Chemisty Lectue #3 page MAY ELECTRO ATOMS At this point, we see that quantum mechanics allows us to undestand the helium atom, at least qualitatively. What about atoms with moe than two electons,
More informationQualifying Examination Electricity and Magnetism Solutions January 12, 2006
1 Qualifying Examination Electicity and Magnetism Solutions Januay 12, 2006 PROBLEM EA. a. Fist, we conside a unit length of cylinde to find the elationship between the total chage pe unit length λ and
More informationPROBLEM SET #1 SOLUTIONS by Robert A. DiStasio Jr.
POBLM S # SOLUIONS by obet A. DiStasio J. Q. he Bon-Oppenheime appoximation is the standad way of appoximating the gound state of a molecula system. Wite down the conditions that detemine the tonic and
More informationEM Boundary Value Problems
EM Bounday Value Poblems 10/ 9 11/ By Ilekta chistidi & Lee, Seung-Hyun A. Geneal Desciption : Maxwell Equations & Loentz Foce We want to find the equations of motion of chaged paticles. The way to do
More informationanubhavclasses.wordpress.com CBSE Solved Test Papers PHYSICS Class XII Chapter : Electrostatics
CBS Solved Test Papes PHYSICS Class XII Chapte : lectostatics CBS TST PAPR-01 CLASS - XII PHYSICS (Unit lectostatics) 1. Show does the foce between two point chages change if the dielectic constant of
More informationPhysics 235 Chapter 5. Chapter 5 Gravitation
Chapte 5 Gavitation In this Chapte we will eview the popeties of the gavitational foce. The gavitational foce has been discussed in geat detail in you intoductoy physics couses, and we will pimaily focus
More informationLecture 8 - Gauss s Law
Lectue 8 - Gauss s Law A Puzzle... Example Calculate the potential enegy, pe ion, fo an infinite 1D ionic cystal with sepaation a; that is, a ow of equally spaced chages of magnitude e and altenating sign.
More informationLet us next consider how one can calculate one- electron elements of the long- range nuclear- electron potential:
Hatee- Fock Long Range In the following assume the patition of 1/ into shot ange and long ance components 1 efc( / a) ef ( / a) + Hee a is a paamete detemining the patitioning between shot ange and long
More informationAST 121S: The origin and evolution of the Universe. Introduction to Mathematical Handout 1
Please ead this fist... AST S: The oigin and evolution of the Univese Intoduction to Mathematical Handout This is an unusually long hand-out and one which uses in places mathematics that you may not be
More informationElectrostatics. 1. Show does the force between two point charges change if the dielectric constant of the medium in which they are kept increase?
Electostatics 1. Show does the foce between two point chages change if the dielectic constant of the medium in which they ae kept incease? 2. A chaged od P attacts od R whee as P epels anothe chaged od
More informationarxiv: v1 [physics.gen-ph] 18 Aug 2018
Path integal and Sommefeld quantization axiv:1809.04416v1 [physics.gen-ph] 18 Aug 018 Mikoto Matsuda 1, and Takehisa Fujita, 1 Japan Health and Medical technological college, Tokyo, Japan College of Science
More informationON THE TWO-BODY PROBLEM IN QUANTUM MECHANICS
ON THE TWO-BODY PROBLEM IN QUANTUM MECHANICS L. MICU Hoia Hulubei National Institute fo Physics and Nuclea Engineeing, P.O. Box MG-6, RO-0775 Buchaest-Maguele, Romania, E-mail: lmicu@theoy.nipne.o (Received
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department. Problem Set 10 Solutions. r s
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Depatment Physics 8.033 Decembe 5, 003 Poblem Set 10 Solutions Poblem 1 M s y x test paticle The figue above depicts the geomety of the poblem. The position
More informationQuantum Chemistry Notes:
Quantum Chemisty otes: Hatee-Fock equations The Hatee-Fock method is the undelying appoximation to nealy all methods of computational chemisty, both ab initio and semi-empiical. Theefoe, a clea undestanding
More informationMany Electron Atoms. Electrons can be put into approximate orbitals and the properties of the many electron systems can be catalogued
Many Electon Atoms The many body poblem cannot be solved analytically. We content ouselves with developing appoximate methods that can yield quite accuate esults (but usually equie a compute). The electons
More informationPhysics 107 TUTORIAL ASSIGNMENT #8
Physics 07 TUTORIAL ASSIGNMENT #8 Cutnell & Johnson, 7 th edition Chapte 8: Poblems 5,, 3, 39, 76 Chapte 9: Poblems 9, 0, 4, 5, 6 Chapte 8 5 Inteactive Solution 8.5 povides a model fo solving this type
More informationPhysics 505 Homework No. 9 Solutions S9-1
Physics 505 Homewok No 9 s S9-1 1 As pomised, hee is the tick fo summing the matix elements fo the Stak effect fo the gound state of the hydogen atom Recall, we need to calculate the coection to the gound
More information( ) [ ] [ ] [ ] δf φ = F φ+δφ F. xdx.
9. LAGRANGIAN OF THE ELECTROMAGNETIC FIELD In the pevious section the Lagangian and Hamiltonian of an ensemble of point paticles was developed. This appoach is based on a qt. This discete fomulation can
More informationPhysics 221 Lecture 41 Nonlinear Absorption and Refraction
Physics 221 Lectue 41 Nonlinea Absoption and Refaction Refeences Meye-Aendt, pp. 97-98. Boyd, Nonlinea Optics, 1.4 Yaiv, Optical Waves in Cystals, p. 22 (Table of cystal symmeties) 1. Intoductoy Remaks.
More informationCharges, Coulomb s Law, and Electric Fields
Q&E -1 Chages, Coulomb s Law, and Electic ields Some expeimental facts: Expeimental fact 1: Electic chage comes in two types, which we call (+) and (). An atom consists of a heavy (+) chaged nucleus suounded
More informationAs is natural, our Aerospace Structures will be described in a Euclidean three-dimensional space R 3.
Appendix A Vecto Algeba As is natual, ou Aeospace Stuctues will be descibed in a Euclidean thee-dimensional space R 3. A.1 Vectos A vecto is used to epesent quantities that have both magnitude and diection.
More informationNewton s Laws, Kepler s Laws, and Planetary Orbits
Newton s Laws, Keple s Laws, and Planetay Obits PROBLEM SET 4 DUE TUESDAY AT START OF LECTURE 28 Septembe 2017 ASTRONOMY 111 FALL 2017 1 Newton s & Keple s laws and planetay obits Unifom cicula motion
More informationRigid Body Dynamics 2. CSE169: Computer Animation Instructor: Steve Rotenberg UCSD, Winter 2018
Rigid Body Dynamics 2 CSE169: Compute Animation nstucto: Steve Rotenbeg UCSD, Winte 2018 Coss Poduct & Hat Opeato Deivative of a Rotating Vecto Let s say that vecto is otating aound the oigin, maintaining
More informationCHAPTER 25 ELECTRIC POTENTIAL
CHPTE 5 ELECTIC POTENTIL Potential Diffeence and Electic Potential Conside a chaged paticle of chage in a egion of an electic field E. This filed exets an electic foce on the paticle given by F=E. When
More informationLecture 23: Central Force Motion
Lectue 3: Cental Foce Motion Many of the foces we encounte in natue act between two paticles along the line connecting the Gavity, electicity, and the stong nuclea foce ae exaples These types of foces
More informationChem 453/544 Fall /08/03. Exam #1 Solutions
Chem 453/544 Fall 3 /8/3 Exam # Solutions. ( points) Use the genealized compessibility diagam povided on the last page to estimate ove what ange of pessues A at oom tempeatue confoms to the ideal gas law
More informationHomework # 3 Solution Key
PHYSICS 631: Geneal Relativity Homewok # 3 Solution Key 1. You e on you hono not to do this one by hand. I ealize you can use a compute o simply look it up. Please don t. In a flat space, the metic in
More informationChapter 13 Gravitation
Chapte 13 Gavitation In this chapte we will exploe the following topics: -Newton s law of gavitation, which descibes the attactive foce between two point masses and its application to extended objects
More informationVectors, Vector Calculus, and Coordinate Systems
Apil 5, 997 A Quick Intoduction to Vectos, Vecto Calculus, and Coodinate Systems David A. Randall Depatment of Atmospheic Science Coloado State Univesity Fot Collins, Coloado 80523. Scalas and vectos Any
More informationNuclear size corrections to the energy levels of single-electron atoms
Nuclea size coections to the enegy levels of single-electon atoms Babak Nadii Nii a eseach Institute fo Astonomy and Astophysics of Maagha (IAAM IAN P. O. Box: 554-44. Abstact A study is made of nuclea
More informationA Relativistic Electron in a Coulomb Potential
A Relativistic Electon in a Coulomb Potential Alfed Whitehead Physics 518, Fall 009 The Poblem Solve the Diac Equation fo an electon in a Coulomb potential. Identify the conseved quantum numbes. Specify
More informationComputational Methods of Solid Mechanics. Project report
Computational Methods of Solid Mechanics Poject epot Due on Dec. 6, 25 Pof. Allan F. Bowe Weilin Deng Simulation of adhesive contact with molecula potential Poject desciption In the poject, we will investigate
More information2. Electrostatics. Dr. Rakhesh Singh Kshetrimayum 8/11/ Electromagnetic Field Theory by R. S. Kshetrimayum
2. Electostatics D. Rakhesh Singh Kshetimayum 1 2.1 Intoduction In this chapte, we will study how to find the electostatic fields fo vaious cases? fo symmetic known chage distibution fo un-symmetic known
More informationEN40: Dynamics and Vibrations. Midterm Examination Thursday March
EN40: Dynamics and Vibations Midtem Examination Thusday Mach 9 2017 School of Engineeing Bown Univesity NAME: Geneal Instuctions No collaboation of any kind is pemitted on this examination. You may bing
More informationAbsorption Rate into a Small Sphere for a Diffusing Particle Confined in a Large Sphere
Applied Mathematics, 06, 7, 709-70 Published Online Apil 06 in SciRes. http://www.scip.og/jounal/am http://dx.doi.og/0.46/am.06.77065 Absoption Rate into a Small Sphee fo a Diffusing Paticle Confined in
More informationLecture 3.7 ELECTRICITY. Electric charge Coulomb s law Electric field
Lectue 3.7 ELECTRICITY Electic chage Coulomb s law Electic field ELECTRICITY Inteaction between electically chages objects Many impotant uses Light Heat Rail tavel Computes Cental nevous system Human body
More informationBut for simplicity, we ll define significant as the time it takes a star to lose all memory of its original trajectory, i.e.,
Stella elaxation Time [Chandasekha 1960, Pinciples of Stella Dynamics, Chap II] [Ostike & Davidson 1968, Ap.J., 151, 679] Do stas eve collide? Ae inteactions between stas (as opposed to the geneal system
More informationChapter 22: Electric Fields. 22-1: What is physics? General physics II (22102) Dr. Iyad SAADEDDIN. 22-2: The Electric Field (E)
Geneal physics II (10) D. Iyad D. Iyad Chapte : lectic Fields In this chapte we will cove The lectic Field lectic Field Lines -: The lectic Field () lectic field exists in a egion of space suounding a
More informationb Ψ Ψ Principles of Organic Chemistry lecture 22, page 1
Pinciples of Oganic Chemisty lectue, page. Basis fo LCAO and Hückel MO Theoy.. Souces... Hypephysics online. http://hypephysics.phy-ast.gsu.edu/hbase/quantum/qm.html#c... Zimmeman, H. E., Quantum Mechanics
More informationNumerical solution of diffusion mass transfer model in adsorption systems. Prof. Nina Paula Gonçalves Salau, D.Sc.
Numeical solution of diffusion mass tansfe model in adsoption systems Pof., D.Sc. Agenda Mass Tansfe Mechanisms Diffusion Mass Tansfe Models Solving Diffusion Mass Tansfe Models Paamete Estimation 2 Mass
More informationNuclear reactions of heavy ions
Autho: Facultat de Física, Univesitat de Bacelona, Diagonal 645, 08028 Bacelona, Spain. Adviso: Xavie Vinyes Abstact: In this wok nuclea eactions of heavy ions ae studied, focusing on elastic scatteing.
More informationTo Feel a Force Chapter 7 Static equilibrium - torque and friction
To eel a oce Chapte 7 Chapte 7: Static fiction, toque and static equilibium A. Review of foce vectos Between the eath and a small mass, gavitational foces of equal magnitude and opposite diection act on
More informationKEPLER S LAWS OF PLANETARY MOTION
EPER S AWS OF PANETARY MOTION 1. Intoduction We ae now in a position to apply what we have leaned about the coss poduct and vecto valued functions to deive eple s aws of planetay motion. These laws wee
More informationF(r) = r f (r) 4.8. Central forces The most interesting problems in classical mechanics are about central forces.
4.8. Cental foces The most inteesting poblems in classical mechanics ae about cental foces. Definition of a cental foce: (i) the diection of the foce F() is paallel o antipaallel to ; in othe wods, fo
More informationCOMP Parallel Computing SMM (3) OpenMP Case Study: The Barnes-Hut N-body Algorithm
COMP 633 - Paallel Computing Lectue 8 Septembe 14, 2017 SMM (3) OpenMP Case Study: The Banes-Hut N-body Algoithm Topics Case study: the Banes-Hut algoithm Study an impotant algoithm in scientific computing»
More informationVoltage ( = Electric Potential )
V-1 of 10 Voltage ( = lectic Potential ) An electic chage altes the space aound it. Thoughout the space aound evey chage is a vecto thing called the electic field. Also filling the space aound evey chage
More informationBasic Principles in Molecular Modeling
Basic Pinciples in Molecula Modeling Univesity of Heidelbeg Molecula Biotechnology (WS 003/004) Mathematical Methods in Bioinfomatics D. M. Diehl by Bama Mostofian (85360) and Filipp Fan (03674) Contents.
More informationd 2 x 0a d d =0. Relative to an arbitrary (accelerating frame) specified by x a = x a (x 0b ), the latter becomes: d 2 x a d 2 + a dx b dx c
Chapte 6 Geneal Relativity 6.1 Towads the Einstein equations Thee ae seveal ways of motivating the Einstein equations. The most natual is pehaps though consideations involving the Equivalence Pinciple.
More informationElectromagnetism Physics 15b
lectomagnetism Physics 15b Lectue #20 Dielectics lectic Dipoles Pucell 10.1 10.6 What We Did Last Time Plane wave solutions of Maxwell s equations = 0 sin(k ωt) B = B 0 sin(k ωt) ω = kc, 0 = B, 0 ˆk =
More informationis the instantaneous position vector of any grid point or fluid
Absolute inetial, elative inetial and non-inetial coodinates fo a moving but non-defoming contol volume Tao Xing, Pablo Caica, and Fed Sten bjective Deive and coelate the govening equations of motion in
More informationAlgebra-based Physics II
lgebabased Physics II Chapte 19 Electic potential enegy & The Electic potential Why enegy is stoed in an electic field? How to descibe an field fom enegetic point of view? Class Website: Natual way of
More informationDonnishJournals
DonnishJounals 041-1189 Donnish Jounal of Educational Reseach and Reviews. Vol 1(1) pp. 01-017 Novembe, 014. http:///dje Copyight 014 Donnish Jounals Oiginal Reseach Pape Vecto Analysis Using MAXIMA Savaş
More informationFZX: Personal Lecture Notes from Daniel W. Koon St. Lawrence University Physics Department CHAPTER 7
FZX: Pesonal Lectue Notes fom Daniel W. Koon St. Lawence Univesity Physics Depatment CHAPTER 7 Please epot any glitches, bugs o eos to the autho: dkoon at stlawu.edu. 7. Momentum and Impulse Impulse page
More information1D2G - Numerical solution of the neutron diffusion equation
DG - Numeical solution of the neuton diffusion equation Y. Danon Daft: /6/09 Oveview A simple numeical solution of the neuton diffusion equation in one dimension and two enegy goups was implemented. Both
More informationQuantum Mechanics II
Quantum Mechanics II Pof. Bois Altshule Apil 25, 2 Lectue 25 We have been dicussing the analytic popeties of the S-matix element. Remembe the adial wave function was u kl () = R kl () e ik iπl/2 S l (k)e
More informationRydberg-Rydberg Interactions
Rydbeg-Rydbeg Inteactions F. Robicheaux Aubun Univesity Rydbeg gas goes to plasma Dipole blockade Coheent pocesses in fozen Rydbeg gases (expts) Theoetical investigation of an excitation hopping though
More informationConservative Averaging Method and its Application for One Heat Conduction Problem
Poceedings of the 4th WSEAS Int. Conf. on HEAT TRANSFER THERMAL ENGINEERING and ENVIRONMENT Elounda Geece August - 6 (pp6-) Consevative Aveaging Method and its Application fo One Heat Conduction Poblem
More informationAn Exact Solution of Navier Stokes Equation
An Exact Solution of Navie Stokes Equation A. Salih Depatment of Aeospace Engineeing Indian Institute of Space Science and Technology, Thiuvananthapuam, Keala, India. July 20 The pincipal difficulty in
More informationPhysics 121 Hour Exam #5 Solution
Physics 2 Hou xam # Solution This exam consists of a five poblems on five pages. Point values ae given with each poblem. They add up to 99 points; you will get fee point to make a total of. In any given
More informationSection 1: Main results of Electrostatics and Magnetostatics. Electrostatics
Chage density ection 1: ain esults of Electostatics and agnetostatics Electostatics The most fundamental quantity of electostatics is electic chage. Chage comes in two vaieties, which ae called positive
More information1 Spherical multipole moments
Jackson notes 9 Spheical multipole moments Suppose we have a chage distibution ρ (x) wheeallofthechageiscontained within a spheical egion of adius R, as shown in the diagam. Then thee is no chage in the
More informationDOING PHYSICS WITH MATLAB COMPUTATIONAL OPTICS
DOING PHYIC WITH MTLB COMPUTTIONL OPTIC FOUNDTION OF CLR DIFFRCTION THEORY Ian Coope chool of Physics, Univesity of ydney ian.coope@sydney.edu.au DOWNLOD DIRECTORY FOR MTLB CRIPT View document: Numeical
More informationPhysics 161 Fall 2011 Extra Credit 2 Investigating Black Holes - Solutions The Following is Worth 50 Points!!!
Physics 161 Fall 011 Exta Cedit Investigating Black Holes - olutions The Following is Woth 50 Points!!! This exta cedit assignment will investigate vaious popeties of black holes that we didn t have time
More informationPhysics 506 Winter 2006 Homework Assignment #9 Solutions
Physics 506 Winte 2006 Homewok Assignment #9 Solutions Textbook poblems: Ch. 12: 12.2, 12.9, 12.13, 12.14 12.2 a) Show fom Hamilton s pinciple that Lagangians that diffe only by a total time deivative
More information17.1 Electric Potential Energy. Equipotential Lines. PE = energy associated with an arrangement of objects that exert forces on each other
Electic Potential Enegy, PE Units: Joules Electic Potential, Units: olts 17.1 Electic Potential Enegy Electic foce is a consevative foce and so we can assign an electic potential enegy (PE) to the system
More informationAppendix B The Relativistic Transformation of Forces
Appendix B The Relativistic Tansfomation of oces B. The ou-foce We intoduced the idea of foces in Chapte 3 whee we saw that the change in the fou-momentum pe unit time is given by the expession d d w x
More informationMath Notes on Kepler s first law 1. r(t) kp(t)
Math 7 - Notes on Keple s fist law Planetay motion and Keple s Laws We conside the motion of a single planet about the sun; fo simplicity, we assign coodinates in R 3 so that the position of the sun is
More informationINTRODUCTION. 2. Vectors in Physics 1
INTRODUCTION Vectos ae used in physics to extend the study of motion fom one dimension to two dimensions Vectos ae indispensable when a physical quantity has a diection associated with it As an example,
More informationTHE LAPLACE EQUATION. The Laplace (or potential) equation is the equation. u = 0. = 2 x 2. x y 2 in R 2
THE LAPLACE EQUATION The Laplace (o potential) equation is the equation whee is the Laplace opeato = 2 x 2 u = 0. in R = 2 x 2 + 2 y 2 in R 2 = 2 x 2 + 2 y 2 + 2 z 2 in R 3 The solutions u of the Laplace
More informationBetween any two masses, there exists a mutual attractive force.
YEAR 12 PHYSICS: GRAVITATION PAST EXAM QUESTIONS Name: QUESTION 1 (1995 EXAM) (a) State Newton s Univesal Law of Gavitation in wods Between any two masses, thee exists a mutual attactive foce. This foce
More informationChapter 5 Force and Motion
Chapte 5 Foce and Motion In Chaptes 2 and 4 we have studied kinematics, i.e., we descibed the motion of objects using paametes such as the position vecto, velocity, and acceleation without any insights
More informationElectric Field. y s +q. Point charge: Uniformly charged sphere: Dipole: for r>>s :! ! E = 1. q 1 r 2 ˆr. E sphere. at <0,r,0> at <0,0,r>
Electic Field Point chage: E " ˆ Unifomly chaged sphee: E sphee E sphee " Q ˆ fo >R (outside) fo >s : E " s 3,, at z y s + x Dipole moment: p s E E s "#,, 3 s "#,, 3 at
More informationNuclear models: Shell model
Lectue 3 Nuclea models: Shell model WS0/3: Intoduction to Nuclea and Paticle Physics,, Pat I Nuclea models Nuclea models Models with stong inteaction between the nucleons Liquid dop model α-paticle model
More informationSingle Particle State AB AB
LECTURE 3 Maxwell Boltzmann, Femi, and Bose Statistics Suppose we have a gas of N identical point paticles in a box of volume V. When we say gas, we mean that the paticles ae not inteacting with one anothe.
More information11) A thin, uniform rod of mass M is supported by two vertical strings, as shown below.
Fall 2007 Qualifie Pat II 12 minute questions 11) A thin, unifom od of mass M is suppoted by two vetical stings, as shown below. Find the tension in the emaining sting immediately afte one of the stings
More informationChapter 5 Force and Motion
Chapte 5 Foce and Motion In chaptes 2 and 4 we have studied kinematics i.e. descibed the motion of objects using paametes such as the position vecto, velocity and acceleation without any insights as to
More informationI. CONSTRUCTION OF THE GREEN S FUNCTION
I. CONSTRUCTION OF THE GREEN S FUNCTION The Helmohltz equation in 4 dimensions is 4 + k G 4 x, x = δ 4 x x. In this equation, G is the Geen s function and 4 efes to the dimensionality. In the vey end,
More informationThree dimensional flow analysis in Axial Flow Compressors
1 Thee dimensional flow analysis in Axial Flow Compessos 2 The ealie assumption on blade flow theoies that the flow inside the axial flow compesso annulus is two dimensional means that adial movement of
More informationTeoría del Funcional de la Densidad (Density Functional Theory)
Teoía del Funcional de la Densidad (Density Functional Theoy) Motivation: limitations of the standad appoach based on the wave function. The electonic density n() as the key vaiable: Functionals & Thomas-Femi
More informationAb Initio Calculations of Intermolecular Interactions. calculating dispersion energies is hard; (BSSE)
V()/k B / K Ab Initio Calculations of Intemolecula Inteactions 100 80 0 40 20 0 Calculated Ne 2 Potentials basis=aug-cc-vqz HF MP2 QCISD(T) B3LYP V() / k B /K 100 80 0 40 20 0 Ne 2 Potentials - Calc &
More informationScattering in Three Dimensions
Scatteing in Thee Dimensions Scatteing expeiments ae an impotant souce of infomation about quantum systems, anging in enegy fom vey low enegy chemical eactions to the highest possible enegies at the LHC.
More information3. Electromagnetic Waves II
Lectue 3 - Electomagnetic Waves II 9 3. Electomagnetic Waves II Last time, we discussed the following. 1. The popagation of an EM wave though a macoscopic media: We discussed how the wave inteacts with
More informationyou of a spring. The potential energy for a spring is given by the parabola U( x)
Small oscillations The theoy of small oscillations is an extemely impotant topic in mechanics. Conside a system that has a potential enegy diagam as below: U B C A x Thee ae thee points of stable equilibium,
More informationFall 2016 Semester METR 3113 Atmospheric Dynamics I: Introduction to Atmospheric Kinematics and Dynamics
Fall 06 Semeste METR 33 Atmospheic Dynamics I: Intoduction to Atmospheic Kinematics Dynamics Lectue 7 Octobe 3 06 Topics: Scale analysis of the equations of hoizontal motion Geostophic appoximation eostophic
More informationPhysics 202, Lecture 2
Physics 202, Lectue 2 Todays Topics Electic Foce and Electic Fields Electic Chages and Electic Foces Coulomb's Law Physical Field The Electic Field Electic Field Lines Motion of Chaged Paticle in Electic
More informationAP Physics - Coulomb's Law
AP Physics - oulomb's Law We ve leaned that electons have a minus one chage and potons have a positive one chage. This plus and minus one business doesn t wok vey well when we go in and ty to do the old
More informationDoublet structure of Alkali spectra:
Doublet stuctue of : Caeful examination of the specta of alkali metals shows that each membe of some of the seies ae closed doublets. Fo example, sodium yellow line, coesponding to 3p 3s tansition, is
More informationPhysics 111. Ch 12: Gravity. Newton s Universal Gravity. R - hat. the equation. = Gm 1 m 2. F g 2 1. ˆr 2 1. Gravity G =
ics Announcements day, embe 9, 004 Ch 1: Gavity Univesal Law Potential Enegy Keple s Laws Ch 15: Fluids density hydostatic equilibium Pascal s Pinciple This week s lab will be anothe physics wokshop -
More information2 Lecture 2: The Bohr atom (1913) and the Schrödinger equation (1925)
1 Lectue 1: The beginnings of quantum physics 1. The Sten-Gelach expeiment. Atomic clocks 3. Planck 1900, blackbody adiation, and E ω 4. Photoelectic effect 5. Electon diffaction though cystals, de Boglie
More informationReview: Electrostatics and Magnetostatics
Review: Electostatics and Magnetostatics In the static egime, electomagnetic quantities do not vay as a function of time. We have two main cases: ELECTROSTATICS The electic chages do not change postion
More informationLecture 4 Povh Krane Enge Williams
Lectue 4 Povh Kane Enge Williams the Deuteon 6. Ch. 4 Ch. Ch 3 d-wave admixtue 4..6 3.5 tenso foce 4..6 3.5 missing S state 4.4.5 3.4 isospin.3 6.7 3.4 Poblems on Lectue 4 What is the minimum photon enegy
More informationModeling Fermi Level Effects in Atomistic Simulations
Mat. Res. Soc. Symp. Poc. Vol. 717 Mateials Reseach Society Modeling Femi Level Effects in Atomistic Simulations Zudian Qin and Scott T. Dunham Depatment of Electical Engineeing, Univesity of Washington,
More informationMechanics Physics 151
Mechanics Physics 151 Lectue 5 Cental Foce Poblem (Chapte 3) What We Did Last Time Intoduced Hamilton s Pinciple Action integal is stationay fo the actual path Deived Lagange s Equations Used calculus
More information