CE 530 Molecular Simulation

Size: px
Start display at page:

Download "CE 530 Molecular Simulation"

Transcription

1 CE 53 Molcular Simulation Lctur 8 Fr-nrgy calculations David A. Kofk Dpartmnt of Chmical Enginring SUNY Buffalo kofk@ng.buffalo.du

2 2 Fr-Enrgy Calculations Uss of fr nrgy Phas quilibria Raction quilibria Solvation Stability Kintics Calculation mthods Fr-nrgy prturbation Thrmodynamic intgration Paramtr-hopping Histogram intrpolation

3 3 Ensmbl Avrags Simpl nsmbl avrags ar of th form M dπ ( M( To valuat: sampl points in phas spac with probability π( at ach point, valuat M( simpl avrag of all valus givs <M> Prvious xampl man squar distanc from origin in rgion R 2 s( 2 ds( r d r ( insid R s outsid R sampl only points in R, avrag r 2 Principl applis to both MD and MC

4 Ensmbl Volums Entropy and fr nrgy rlat to th siz of th nsmbl.g., S k lnω(e,v,n No ffctiv way to masur th siz of th nsmbl no phas-spac function that givs siz of R whil sampling only R imagin bing plac rpatdly at random points on an island what could you masur at ach point to dtrmin th siz of th island? Volum of nsmbl is numrically unwildy.g. for hard sphrs r., Ω 5 33 Ω numbr of stats of givn E,V,N 4 r.5, Ω 3 7 r.9, Ω 5-42 Shap of important rgion is vry complx cannot apply mthods that xploit som simpl gomtric pictur

5 Rfrnc Systms All fr-nrgy mthods ar basd on calculation of frnrgy diffrncs Exampl volum of R can b masurd as a fraction of th total volum 5 Ω R s( Ω sampl th rfrnc systm kp an avrag of th fraction of tim occupying targt systm what w gt is th diffrnc SR S kln ( ΩR / Ω Usfulnss of fr-nrgy diffrnc it may b th quantity of intrst anyway if rfrnc is simpl, its absolut fr nrgy can b valuatd analytically.g., idal gas, harmonic crystal

6 Hard Sphr Chmical Potntial Chmical potntial is an ntropy diffrnc S/ k βµ + N UV, For hard sphrs, th nrgy is zro or infinity any chang in N that dos not caus ovrlap will b chang at constant U To gt ntropy diffrnc [ SUV (,, N SUV (,, N ] simulat a systm of N+ sphrs, on non-intracting ghost occasionally s if th ghost sphr ovrlaps anothr rcord th fraction of th tim it dos not ovrlap S/ k V Δ βµ Ω N + V f 3 V 3 NΛ Ω 3 N NΛ NΛ non-ovrlap Hr is an applt dmonstrating this calculation N+ N 6

7 Fr-Enrgy Prturbation Widom mthod is an xampl of a fr-nrgy prturbation (FEP tchniqu FEP givs fr-nrgy diffrnc btwn two systms labld, Working quation ( A A Q Q d d π ( d d ( U U U d ( U U ( U U U U U Fr-nrgy diffrnc is a ratio of partition functions 7

8 Fr-Enrgy Prturbation Widom mthod is an xampl of a fr-nrgy prturbation (FEP tchniqu FEP givs fr-nrgy diffrnc btwn two systms labld, Working quation ( A A Q Q d d d ( U U d π ( d ( U U ( U U U U U U Add and subtract rfrnc-systm nrgy 8

9 Fr-Enrgy Prturbation Widom mthod is an xampl of a fr-nrgy prturbation (FEP tchniqu FEP givs fr-nrgy diffrnc btwn two systms labld, Working quation ( A A Q Q d ( U U U ( U U d π ( ( U U d d d U U U 9 Idntify rfrncsystm probability distribution

10 Fr-Enrgy Prturbation Widom mthod is an xampl of a fr-nrgy prturbation (FEP tchniqu FEP givs fr-nrgy diffrnc btwn two systms labld, Working quation ( A A Q Q Sampl th rgion important to systm, masur proprtis of systm d d d ( U U U d ( U U d π ( ( U U U U U Writ as rfrncsystm nsmbl avrag

11 Chmical potntial For chmical potntial, U - U is th nrgy of turning on th ghost particl call this u t, th tst-particl nrgy ( A A µ V NΛ 3 u t tst-particl position may b slctd at random in simulation volum for hard sphrs, -βu t is for ovrlap, othrwis thn (as bfor avrag is th fraction of configurations with no ovrlap This is known as Widom s insrtion mthod

12 Dltion Mthod Th FEP formula may b usd also with th rols of th rfrnc and targt systm rvrsd Original: à Modifid: à ( A A ( U U sampl th systm, valuat proprtis of systm Considr application to hard sphrs + β( A A + β( U U 3 + βµ N + β V βu t is infinity for ovrlap, othrwis but ovrlaps ar nvr sampld tru avrag is product of tchnically, formula is corrct u t in practic simulation avrag is always zro mthod is flawd in application many tims th flaw with dltion is not as obvious as this

13 4 Othr Typs of Prturbation Many typs of fr-nrgy diffrncs can b computd Thrmodynamic stat tmpratur, dnsity, mixtur composition Hamiltonian for a singl molcul or for ntir systm.g., valuat fr nrgy diffrnc for hard sphrs with and without lctrostatic dipol momnt Configuration distanc/orintation btwn two soluts.g, protin and ligand Ordr paramtr idntifying phass ordr paramtr is a quantity that can b usd to idntify th thrmodynamic phas a systm is in.g, crystal structur, orintational ordr, magntization

14 Gnral Numrical Problms 5 Sampling problms limit rang of FEP calculations Targt systm configurations must b ncountrd whn sampling rfrnc systm Two typs of problm aris targt-systm spac vry small targt systm outsid of rfrnc first situation is mor common although dltion FEP provids an avoidabl xampl of th lattr

15 Staging Mthods 6 Multistag FEP can b usd to rmdy th sampling problm dfin a potntial U w intrmdiat btwn and systms valuat total fr-nrgy diffrnc as Each stag may b sampld in ithr dirction yilding four staging schms choos to avoid dltion calculation A A ( A A + ( A A w w W Umbrlla sampling W Bnntt's mthod W Stagd dltion W Stagd insrtion Us stagd insrtion Us umbrlla sampling Us Bnntt s mthod W W W

16 Exampl of Staging Mthod Hard-sphr chmical potntial Us small-diamtr sphr as intrmdiat β Δ( A A β u( σ w N t In first stag, masur fraction of tim random insrtion of small sphr finds no ovrlap N ( A N Aw ( u( u( t Δ + β σ σ In scond stag, small sphr movs around with othrs. Masur fraction of tim no ovrlap is found whn it is grown to full-siz sphr w 7 ( A N A N u ( t ( u( u( t β Δ + β σ σ σ N w

17 8 Multipl Stags W 3 W 2 W 2 W 3 W W 2 W W 3 W Multistag insrtion Multistag umbrlla sampling Multistag Bnntt s mthod W W2 W3 W2 W W3 W2 W W3

18 9 Non-Boltzmann Sampling Th FEP mthods ar an instanc of a mor gnral tchniqu that aims to improv sampling Unlik biasing mthods, improvmnt ntails a chang in th limiting distribution Apply a formula to rcovr th corrct avrag Q M dm( Q ( U U U d M( Q Q W W W M W ( U U W ( U U W U W W W

19 2 Thrmodynamic Intgration. Thrmodynamics givs formulas for variation of fr nrgy with stat d( βa Udβ βpdv + βµ dn βa βa U P β V VN, TN, Ths can b intgratd to obtain a fr-nrgy diffrnc drivativs can b masurd as normal nsmbl avrags βa( V βa( V P( V dv 2 V V 2 this is usually how fr nrgis ar masurd xprimntally

20 Thrmodynamic Intgration 2. TI can b xtndd to follow uncommon (or unphysical intgration paths much lik FEP, can b applid for any typ of fr-nrgy chang Formalism Lt λ b a paramtr dscribing th path th potntial nrgy is a function of λ nsmbl formula for th drivativ thn β A ln Q N N U( r ; λ dr 3N Λ N! λ λ Q λ Q N N U( r ; λ N + dr βu ( r ; λ 3N Λ βu λ N! 2 βu βa( λ2 βa( λ + dλ λ λ λ λ 2

21 Thrmodynamic Intgration Exampl 22 Th soft-sphr pair potntial is givn by n ur ( σ ε r Exhibits simplifying bhavior bcaus εσ n is th only potntial paramtr Softnss and rang varis with n larg n limit lads to hard sphrs small n lads to Coulombic bhavior Thrmodynamic intgration can b usd to masur fr nrgy as a function of softnss s /n Intgrand is βa s βε ur (ln( σ / r 2 s U/ε Sparation, r/σ n 2 n 2 n 6 n 3 2.5

22 Paramtr Hopping. Thory 23 Viw fr-nrgy paramtr λ as anothr dimnsion in phas spac N E E( p, r, λ Partition function N Q dp dr λ N N E( p, r, λ Mont Carlo trials includ changs in λ Probability that systm has λ λ or λ λ N N N N E N N E( p, r, λ N N β ( p, r, λ dp dr + dp dr Q + Q π( λ π( λ N N N N N E( p, r, λ Q d d Q+ Q p r Q+ Q N N N N E( p, r, λ Q d d Q+ Q p r Q+ Q N λ λ

23 Paramtr Hopping. Implmntation Mont Carlo simulation in which l-chang trials ar attmptd Accpt trials as usual, with probability min[, -βδu ] 24 Rcord fractions f, f of configurations spnt in λ λ and λ λ Fr nrgy is givn by ratio ( A A Q Q ( Q + Q f In practic, systm may spnd almost no tim in on of th valus Can apply wighting function w(λ to ncourag it to sampl both Accpt trials with probability min[,(w n /w o -βδu ] Fr nrgy is Q Q ( Q + Q f Good choic for w has f f Multivalu xtnsion is particularly ffctiv l taks on a continuum of valus ( A A w f wf

24 25 Summary Fr nrgy calculations ar ndd to modl th most intrsting physical bhaviors All usful mthods ar basd on computing fr-nrgy diffrnc Four gnral approachs Fr-nrgy prturbation Thrmodynamic intgration Paramtr hopping Distribution-function mthods FEP is asymmtric Dltion mthod is awful Four approachs to basic multistaging Umbrlla sampling, Bnntt s mthod, stagd insrtion/dltion Non-Boltzmann mthods improv sampling

Lecture 19: Free Energies in Modern Computational Statistical Thermodynamics: WHAM and Related Methods

Lecture 19: Free Energies in Modern Computational Statistical Thermodynamics: WHAM and Related Methods Statistical Thrmodynamics Lctur 19: Fr Enrgis in Modrn Computational Statistical Thrmodynamics: WHAM and Rlatd Mthods Dr. Ronald M. Lvy ronlvy@tmpl.du Dfinitions Canonical nsmbl: A N, V,T = k B T ln Q

More information

Chemical Physics II. More Stat. Thermo Kinetics Protein Folding...

Chemical Physics II. More Stat. Thermo Kinetics Protein Folding... Chmical Physics II Mor Stat. Thrmo Kintics Protin Folding... http://www.nmc.ctc.com/imags/projct/proj15thumb.jpg http://nuclarwaponarchiv.org/usa/tsts/ukgrabl2.jpg http://www.photolib.noaa.gov/corps/imags/big/corp1417.jpg

More information

Elements of Statistical Thermodynamics

Elements of Statistical Thermodynamics 24 Elmnts of Statistical Thrmodynamics Statistical thrmodynamics is a branch of knowldg that has its own postulats and tchniqus. W do not attmpt to giv hr vn an introduction to th fild. In this chaptr,

More information

Lecture 13: Conformational Sampling: MC and MD

Lecture 13: Conformational Sampling: MC and MD Statistical Thrmodynamics Lctur 13: Conformational Sampling: MC and MD Dr. Ronald M. Lvy ronlvy@tmpl.du Contributions from Mik Andrc and Danil Winstock Importanc Sampling and Mont Carlo Mthods Enrgy functions

More information

Ch. 24 Molecular Reaction Dynamics 1. Collision Theory

Ch. 24 Molecular Reaction Dynamics 1. Collision Theory Ch. 4 Molcular Raction Dynamics 1. Collision Thory Lctur 16. Diffusion-Controlld Raction 3. Th Matrial Balanc Equation 4. Transition Stat Thory: Th Eyring Equation 5. Transition Stat Thory: Thrmodynamic

More information

Hydrogen Atom and One Electron Ions

Hydrogen Atom and One Electron Ions Hydrogn Atom and On Elctron Ions Th Schrödingr quation for this two-body problm starts out th sam as th gnral two-body Schrödingr quation. First w sparat out th motion of th cntr of mass. Th intrnal potntial

More information

Brief Introduction to Statistical Mechanics

Brief Introduction to Statistical Mechanics Brif Introduction to Statistical Mchanics. Purpos: Ths nots ar intndd to provid a vry quick introduction to Statistical Mchanics. Th fild is of cours far mor vast than could b containd in ths fw pags.

More information

The van der Waals interaction 1 D. E. Soper 2 University of Oregon 20 April 2012

The van der Waals interaction 1 D. E. Soper 2 University of Oregon 20 April 2012 Th van dr Waals intraction D. E. Sopr 2 Univrsity of Orgon 20 pril 202 Th van dr Waals intraction is discussd in Chaptr 5 of J. J. Sakurai, Modrn Quantum Mchanics. Hr I tak a look at it in a littl mor

More information

Coupled Pendulums. Two normal modes.

Coupled Pendulums. Two normal modes. Tim Dpndnt Two Stat Problm Coupld Pndulums Wak spring Two normal mods. No friction. No air rsistanc. Prfct Spring Start Swinging Som tim latr - swings with full amplitud. stationary M +n L M +m Elctron

More information

Lecture 37 (Schrödinger Equation) Physics Spring 2018 Douglas Fields

Lecture 37 (Schrödinger Equation) Physics Spring 2018 Douglas Fields Lctur 37 (Schrödingr Equation) Physics 6-01 Spring 018 Douglas Filds Rducd Mass OK, so th Bohr modl of th atom givs nrgy lvls: E n 1 k m n 4 But, this has on problm it was dvlopd assuming th acclration

More information

SECTION where P (cos θ, sin θ) and Q(cos θ, sin θ) are polynomials in cos θ and sin θ, provided Q is never equal to zero.

SECTION where P (cos θ, sin θ) and Q(cos θ, sin θ) are polynomials in cos θ and sin θ, provided Q is never equal to zero. SETION 6. 57 6. Evaluation of Dfinit Intgrals Exampl 6.6 W hav usd dfinit intgrals to valuat contour intgrals. It may com as a surpris to larn that contour intgrals and rsidus can b usd to valuat crtain

More information

Data Assimilation 1. Alan O Neill National Centre for Earth Observation UK

Data Assimilation 1. Alan O Neill National Centre for Earth Observation UK Data Assimilation 1 Alan O Nill National Cntr for Earth Obsrvation UK Plan Motivation & basic idas Univariat (scalar) data assimilation Multivariat (vctor) data assimilation 3d-Variational Mthod (& optimal

More information

Middle East Technical University Department of Mechanical Engineering ME 413 Introduction to Finite Element Analysis

Middle East Technical University Department of Mechanical Engineering ME 413 Introduction to Finite Element Analysis Middl East Tchnical Univrsity Dpartmnt of Mchanical Enginring ME 43 Introduction to Finit Elmnt Analysis Chaptr 3 Computr Implmntation of D FEM Ths nots ar prpard by Dr. Cünyt Srt http://www.m.mtu.du.tr/popl/cunyt

More information

Fourier Transforms and the Wave Equation. Key Mathematics: More Fourier transform theory, especially as applied to solving the wave equation.

Fourier Transforms and the Wave Equation. Key Mathematics: More Fourier transform theory, especially as applied to solving the wave equation. Lur 7 Fourir Transforms and th Wav Euation Ovrviw and Motivation: W first discuss a fw faturs of th Fourir transform (FT), and thn w solv th initial-valu problm for th wav uation using th Fourir transform

More information

On the Hamiltonian of a Multi-Electron Atom

On the Hamiltonian of a Multi-Electron Atom On th Hamiltonian of a Multi-Elctron Atom Austn Gronr Drxl Univrsity Philadlphia, PA Octobr 29, 2010 1 Introduction In this papr, w will xhibit th procss of achiving th Hamiltonian for an lctron gas. Making

More information

MCB137: Physical Biology of the Cell Spring 2017 Homework 6: Ligand binding and the MWC model of allostery (Due 3/23/17)

MCB137: Physical Biology of the Cell Spring 2017 Homework 6: Ligand binding and the MWC model of allostery (Due 3/23/17) MCB37: Physical Biology of th Cll Spring 207 Homwork 6: Ligand binding and th MWC modl of allostry (Du 3/23/7) Hrnan G. Garcia March 2, 207 Simpl rprssion In class, w drivd a mathmatical modl of how simpl

More information

3 Finite Element Parametric Geometry

3 Finite Element Parametric Geometry 3 Finit Elmnt Paramtric Gomtry 3. Introduction Th intgral of a matrix is th matrix containing th intgral of ach and vry on of its original componnts. Practical finit lmnt analysis rquirs intgrating matrics,

More information

MCE503: Modeling and Simulation of Mechatronic Systems Discussion on Bond Graph Sign Conventions for Electrical Systems

MCE503: Modeling and Simulation of Mechatronic Systems Discussion on Bond Graph Sign Conventions for Electrical Systems MCE503: Modling and Simulation o Mchatronic Systms Discussion on Bond Graph Sign Convntions or Elctrical Systms Hanz ichtr, PhD Clvland Stat Univrsity, Dpt o Mchanical Enginring 1 Basic Assumption In a

More information

Mor Tutorial at www.dumblittldoctor.com Work th problms without a calculator, but us a calculator to chck rsults. And try diffrntiating your answrs in part III as a usful chck. I. Applications of Intgration

More information

Addition of angular momentum

Addition of angular momentum Addition of angular momntum April, 0 Oftn w nd to combin diffrnt sourcs of angular momntum to charactriz th total angular momntum of a systm, or to divid th total angular momntum into parts to valuat th

More information

Quasi-Classical States of the Simple Harmonic Oscillator

Quasi-Classical States of the Simple Harmonic Oscillator Quasi-Classical Stats of th Simpl Harmonic Oscillator (Draft Vrsion) Introduction: Why Look for Eignstats of th Annihilation Oprator? Excpt for th ground stat, th corrspondnc btwn th quantum nrgy ignstats

More information

u r du = ur+1 r + 1 du = ln u + C u sin u du = cos u + C cos u du = sin u + C sec u tan u du = sec u + C e u du = e u + C

u r du = ur+1 r + 1 du = ln u + C u sin u du = cos u + C cos u du = sin u + C sec u tan u du = sec u + C e u du = e u + C Tchniqus of Intgration c Donald Kridr and Dwight Lahr In this sction w ar going to introduc th first approachs to valuating an indfinit intgral whos intgrand dos not hav an immdiat antidrivativ. W bgin

More information

u x v x dx u x v x v x u x dx d u x v x u x v x dx u x v x dx Integration by Parts Formula

u x v x dx u x v x v x u x dx d u x v x u x v x dx u x v x dx Integration by Parts Formula 7. Intgration by Parts Each drivativ formula givs ris to a corrsponding intgral formula, as w v sn many tims. Th drivativ product rul yilds a vry usful intgration tchniqu calld intgration by parts. Starting

More information

AS 5850 Finite Element Analysis

AS 5850 Finite Element Analysis AS 5850 Finit Elmnt Analysis Two-Dimnsional Linar Elasticity Instructor Prof. IIT Madras Equations of Plan Elasticity - 1 displacmnt fild strain- displacmnt rlations (infinitsimal strain) in matrix form

More information

Title: Vibrational structure of electronic transition

Title: Vibrational structure of electronic transition Titl: Vibrational structur of lctronic transition Pag- Th band spctrum sn in th Ultra-Violt (UV) and visibl (VIS) rgions of th lctromagntic spctrum can not intrprtd as vibrational and rotational spctrum

More information

Bifurcation Theory. , a stationary point, depends on the value of α. At certain values

Bifurcation Theory. , a stationary point, depends on the value of α. At certain values Dnamic Macroconomic Thor Prof. Thomas Lux Bifurcation Thor Bifurcation: qualitativ chang in th natur of th solution occurs if a paramtr passs through a critical point bifurcation or branch valu. Local

More information

Dealing with quantitative data and problem solving life is a story problem! Attacking Quantitative Problems

Dealing with quantitative data and problem solving life is a story problem! Attacking Quantitative Problems Daling with quantitati data and problm soling lif is a story problm! A larg portion of scinc inols quantitati data that has both alu and units. Units can sa your butt! Nd handl on mtric prfixs Dimnsional

More information

Derivation of Electron-Electron Interaction Terms in the Multi-Electron Hamiltonian

Derivation of Electron-Electron Interaction Terms in the Multi-Electron Hamiltonian Drivation of Elctron-Elctron Intraction Trms in th Multi-Elctron Hamiltonian Erica Smith Octobr 1, 010 1 Introduction Th Hamiltonian for a multi-lctron atom with n lctrons is drivd by Itoh (1965) by accounting

More information

orbiting electron turns out to be wrong even though it Unfortunately, the classical visualization of the

orbiting electron turns out to be wrong even though it Unfortunately, the classical visualization of the Lctur 22-1 Byond Bohr Modl Unfortunatly, th classical visualization of th orbiting lctron turns out to b wrong vn though it still givs us a simpl way to think of th atom. Quantum Mchanics is ndd to truly

More information

Molecular Orbitals in Inorganic Chemistry

Molecular Orbitals in Inorganic Chemistry Outlin olcular Orbitals in Inorganic Chmistry Dr. P. Hunt p.hunt@imprial.ac.uk Rm 167 (Chmistry) http://www.ch.ic.ac.uk/hunt/ octahdral complxs forming th O diagram for Oh colour, slction ruls Δoct, spctrochmical

More information

6.1 Integration by Parts and Present Value. Copyright Cengage Learning. All rights reserved.

6.1 Integration by Parts and Present Value. Copyright Cengage Learning. All rights reserved. 6.1 Intgration by Parts and Prsnt Valu Copyright Cngag Larning. All rights rsrvd. Warm-Up: Find f () 1. F() = ln(+1). F() = 3 3. F() =. F() = ln ( 1) 5. F() = 6. F() = - Objctivs, Day #1 Studnts will b

More information

Addition of angular momentum

Addition of angular momentum Addition of angular momntum April, 07 Oftn w nd to combin diffrnt sourcs of angular momntum to charactriz th total angular momntum of a systm, or to divid th total angular momntum into parts to valuat

More information

Outline. Thanks to Ian Blockland and Randy Sobie for these slides Lifetimes of Decaying Particles Scattering Cross Sections Fermi s Golden Rule

Outline. Thanks to Ian Blockland and Randy Sobie for these slides Lifetimes of Decaying Particles Scattering Cross Sections Fermi s Golden Rule Outlin Thanks to Ian Blockland and andy obi for ths slids Liftims of Dcaying Particls cattring Cross ctions Frmi s Goldn ul Physics 424 Lctur 12 Pag 1 Obsrvabls want to rlat xprimntal masurmnts to thortical

More information

5.80 Small-Molecule Spectroscopy and Dynamics

5.80 Small-Molecule Spectroscopy and Dynamics MIT OpnCoursWar http://ocw.mit.du 5.80 Small-Molcul Spctroscopy and Dynamics Fall 008 For information about citing ths matrials or our Trms of Us, visit: http://ocw.mit.du/trms. Lctur # 3 Supplmnt Contnts

More information

MATH 1080 Test 2-SOLUTIONS Spring

MATH 1080 Test 2-SOLUTIONS Spring MATH Tst -SOLUTIONS Spring 5. Considr th curv dfind by x = ln( 3y + 7) on th intrval y. a. (5 points) St up but do not simplify or valuat an intgral rprsnting th lngth of th curv on th givn intrval. =

More information

Where k is either given or determined from the data and c is an arbitrary constant.

Where k is either given or determined from the data and c is an arbitrary constant. Exponntial growth and dcay applications W wish to solv an quation that has a drivativ. dy ky k > dx This quation says that th rat of chang of th function is proportional to th function. Th solution is

More information

COMPUTER GENERATED HOLOGRAMS Optical Sciences 627 W.J. Dallas (Monday, April 04, 2005, 8:35 AM) PART I: CHAPTER TWO COMB MATH.

COMPUTER GENERATED HOLOGRAMS Optical Sciences 627 W.J. Dallas (Monday, April 04, 2005, 8:35 AM) PART I: CHAPTER TWO COMB MATH. C:\Dallas\0_Courss\03A_OpSci_67\0 Cgh_Book\0_athmaticalPrliminaris\0_0 Combath.doc of 8 COPUTER GENERATED HOLOGRAS Optical Scincs 67 W.J. Dallas (onday, April 04, 005, 8:35 A) PART I: CHAPTER TWO COB ATH

More information

Exam 1. It is important that you clearly show your work and mark the final answer clearly, closed book, closed notes, no calculator.

Exam 1. It is important that you clearly show your work and mark the final answer clearly, closed book, closed notes, no calculator. Exam N a m : _ S O L U T I O N P U I D : I n s t r u c t i o n s : It is important that you clarly show your work and mark th final answr clarly, closd book, closd nots, no calculator. T i m : h o u r

More information

Search sequence databases 3 10/25/2016

Search sequence databases 3 10/25/2016 Sarch squnc databass 3 10/25/2016 Etrm valu distribution Ø Suppos X is a random variabl with probability dnsity function p(, w sampl a larg numbr S of indpndnt valus of X from this distribution for an

More information

Finite element discretization of Laplace and Poisson equations

Finite element discretization of Laplace and Poisson equations Finit lmnt discrtization of Laplac and Poisson quations Yashwanth Tummala Tutor: Prof S.Mittal 1 Outlin Finit Elmnt Mthod for 1D Introduction to Poisson s and Laplac s Equations Finit Elmnt Mthod for 2D-Discrtization

More information

Principles of Humidity Dalton s law

Principles of Humidity Dalton s law Principls of Humidity Dalton s law Air is a mixtur of diffrnt gass. Th main gas componnts ar: Gas componnt volum [%] wight [%] Nitrogn N 2 78,03 75,47 Oxygn O 2 20,99 23,20 Argon Ar 0,93 1,28 Carbon dioxid

More information

Integration by Parts

Integration by Parts Intgration by Parts Intgration by parts is a tchniqu primarily for valuating intgrals whos intgrand is th product of two functions whr substitution dosn t work. For ampl, sin d or d. Th rul is: u ( ) v'(

More information

Byeong-Joo Lee

Byeong-Joo Lee OSECH - MSE calphad@postch.ac.kr Equipartition horm h avrag nrgy o a particl pr indpndnt componnt o motion is ε ε ' ε '' ε ''' U ln Z Z ε < ε > U ln Z β ( ε ' ε '' ε ''' / Z' Z translational kintic nrgy

More information

Probability and Stochastic Processes: A Friendly Introduction for Electrical and Computer Engineers Roy D. Yates and David J.

Probability and Stochastic Processes: A Friendly Introduction for Electrical and Computer Engineers Roy D. Yates and David J. Probability and Stochastic Procsss: A Frindly Introduction for Elctrical and Computr Enginrs Roy D. Yats and David J. Goodman Problm Solutions : Yats and Goodman,4.3. 4.3.4 4.3. 4.4. 4.4.4 4.4.6 4.. 4..7

More information

Statistical Thermodynamics: Sublimation of Solid Iodine

Statistical Thermodynamics: Sublimation of Solid Iodine c:374-7-ivap-statmch.docx mar7 Statistical Thrmodynamics: Sublimation of Solid Iodin Chm 374 For March 3, 7 Prof. Patrik Callis Purpos:. To rviw basic fundamntals idas of Statistical Mchanics as applid

More information

2008 AP Calculus BC Multiple Choice Exam

2008 AP Calculus BC Multiple Choice Exam 008 AP Multipl Choic Eam Nam 008 AP Calculus BC Multipl Choic Eam Sction No Calculator Activ AP Calculus 008 BC Multipl Choic. At tim t 0, a particl moving in th -plan is th acclration vctor of th particl

More information

Forces. Quantum ElectroDynamics. α = = We have now:

Forces. Quantum ElectroDynamics. α = = We have now: W hav now: Forcs Considrd th gnral proprtis of forcs mdiatd by xchang (Yukawa potntial); Examind consrvation laws which ar obyd by (som) forcs. W will nxt look at thr forcs in mor dtail: Elctromagntic

More information

Cramér-Rao Inequality: Let f(x; θ) be a probability density function with continuous parameter

Cramér-Rao Inequality: Let f(x; θ) be a probability density function with continuous parameter WHEN THE CRAMÉR-RAO INEQUALITY PROVIDES NO INFORMATION STEVEN J. MILLER Abstract. W invstigat a on-paramtr family of probability dnsitis (rlatd to th Parto distribution, which dscribs many natural phnomna)

More information

Pipe flow friction, small vs. big pipes

Pipe flow friction, small vs. big pipes Friction actor (t/0 t o pip) Friction small vs larg pips J. Chaurtt May 016 It is an intrsting act that riction is highr in small pips than largr pips or th sam vlocity o low and th sam lngth. Friction

More information

ECE507 - Plasma Physics and Applications

ECE507 - Plasma Physics and Applications ECE507 - Plasma Physics and Applications Lctur 7 Prof. Jorg Rocca and Dr. Frnando Tomasl Dpartmnt of Elctrical and Computr Enginring Collisional and radiativ procsss All particls in a plasma intract with

More information

The Matrix Exponential

The Matrix Exponential Th Matrix Exponntial (with xrciss) by D. Klain Vrsion 207.0.05 Corrctions and commnts ar wlcom. Th Matrix Exponntial For ach n n complx matrix A, dfin th xponntial of A to b th matrix A A k I + A + k!

More information

EE243 Advanced Electromagnetic Theory Lec # 23 Scattering and Diffraction. Reading: Jackson Chapter , lite

EE243 Advanced Electromagnetic Theory Lec # 23 Scattering and Diffraction. Reading: Jackson Chapter , lite Applid M Fall 6, Nuruthr Lctur #3 Vr /5/6 43 Advancd lctromagntic Thory Lc # 3 cattring and Diffraction calar Diffraction Thory Vctor Diffraction Thory Babint and Othr Principls Optical Thorm ading: Jackson

More information

Math 34A. Final Review

Math 34A. Final Review Math A Final Rviw 1) Us th graph of y10 to find approimat valus: a) 50 0. b) y (0.65) solution for part a) first writ an quation: 50 0. now tak th logarithm of both sids: log() log(50 0. ) pand th right

More information

Constants and Conversions:

Constants and Conversions: EXAM INFORMATION Radial Distribution Function: P 2 ( r) RDF( r) Br R( r ) 2, B is th normalization constant. Ordr of Orbital Enrgis: Homonuclar Diatomic Molculs * * * * g1s u1s g 2s u 2s u 2 p g 2 p g

More information

Introduction to Condensed Matter Physics

Introduction to Condensed Matter Physics Introduction to Condnsd Mattr Physics pcific hat M.P. Vaughan Ovrviw Ovrviw of spcific hat Hat capacity Dulong-Ptit Law Einstin modl Dby modl h Hat Capacity Hat capacity h hat capacity of a systm hld at

More information

Chapter 8: Electron Configurations and Periodicity

Chapter 8: Electron Configurations and Periodicity Elctron Spin & th Pauli Exclusion Principl Chaptr 8: Elctron Configurations and Priodicity 3 quantum numbrs (n, l, ml) dfin th nrgy, siz, shap, and spatial orintation of ach atomic orbital. To xplain how

More information

Rational Approximation for the one-dimensional Bratu Equation

Rational Approximation for the one-dimensional Bratu Equation Intrnational Journal of Enginring & Tchnology IJET-IJES Vol:3 o:05 5 Rational Approximation for th on-dimnsional Bratu Equation Moustafa Aly Soliman Chmical Enginring Dpartmnt, Th British Univrsity in

More information

22/ Breakdown of the Born-Oppenheimer approximation. Selection rules for rotational-vibrational transitions. P, R branches.

22/ Breakdown of the Born-Oppenheimer approximation. Selection rules for rotational-vibrational transitions. P, R branches. Subjct Chmistry Papr No and Titl Modul No and Titl Modul Tag 8/ Physical Spctroscopy / Brakdown of th Born-Oppnhimr approximation. Slction ruls for rotational-vibrational transitions. P, R branchs. CHE_P8_M

More information

ME469A Numerical Methods for Fluid Mechanics

ME469A Numerical Methods for Fluid Mechanics ME469A Numrical Mthods for Fluid Mchanics Handout #5 Gianluca Iaccarino Finit Volum Mthods Last tim w introducd th FV mthod as a discrtization tchniqu applid to th intgral form of th govrning quations

More information

cycle that does not cross any edges (including its own), then it has at least

cycle that does not cross any edges (including its own), then it has at least W prov th following thorm: Thorm If a K n is drawn in th plan in such a way that it has a hamiltonian cycl that dos not cross any dgs (including its own, thn it has at last n ( 4 48 π + O(n crossings Th

More information

High Energy Physics. Lecture 5 The Passage of Particles through Matter

High Energy Physics. Lecture 5 The Passage of Particles through Matter High Enrgy Physics Lctur 5 Th Passag of Particls through Mattr 1 Introduction In prvious lcturs w hav sn xampls of tracks lft by chargd particls in passing through mattr. Such tracks provid som of th most

More information

Chapter 13 GMM for Linear Factor Models in Discount Factor form. GMM on the pricing errors gives a crosssectional

Chapter 13 GMM for Linear Factor Models in Discount Factor form. GMM on the pricing errors gives a crosssectional Chaptr 13 GMM for Linar Factor Modls in Discount Factor form GMM on th pricing rrors givs a crosssctional rgrssion h cas of xcss rturns Hors rac sting for charactristic sting for pricd factors: lambdas

More information

Difference -Analytical Method of The One-Dimensional Convection-Diffusion Equation

Difference -Analytical Method of The One-Dimensional Convection-Diffusion Equation Diffrnc -Analytical Mthod of Th On-Dimnsional Convction-Diffusion Equation Dalabav Umurdin Dpartmnt mathmatic modlling, Univrsity of orld Economy and Diplomacy, Uzbistan Abstract. An analytical diffrncing

More information

4037 ADDITIONAL MATHEMATICS

4037 ADDITIONAL MATHEMATICS CAMBRIDGE INTERNATIONAL EXAMINATIONS GCE Ordinary Lvl MARK SCHEME for th Octobr/Novmbr 0 sris 40 ADDITIONAL MATHEMATICS 40/ Papr, maimum raw mark 80 This mark schm is publishd as an aid to tachrs and candidats,

More information

4 x 4, and. where x is Town Square

4 x 4, and. where x is Town Square Accumulation and Population Dnsity E. A city locatd along a straight highway has a population whos dnsity can b approimatd by th function p 5 4 th distanc from th town squar, masurd in mils, whr 4 4, and

More information

5.62 Physical Chemistry II Spring 2008

5.62 Physical Chemistry II Spring 2008 MIT OpnCoursWar http://ocw.mit.du 5.62 Physical Chmistry II Spring 2008 For information about citing ths matrials or our Trms of Us, visit: http://ocw.mit.du/trms. 5.62 Lctur #7: Translational Part of

More information

Lecture 28 Title: Diatomic Molecule : Vibrational and Rotational spectra

Lecture 28 Title: Diatomic Molecule : Vibrational and Rotational spectra Lctur 8 Titl: Diatomic Molcul : Vibrational and otational spctra Pag- In this lctur w will undrstand th molcular vibrational and rotational spctra of diatomic molcul W will start with th Hamiltonian for

More information

Give the letter that represents an atom (6) (b) Atoms of A and D combine to form a compound containing covalent bonds.

Give the letter that represents an atom (6) (b) Atoms of A and D combine to form a compound containing covalent bonds. 1 Th diagram shows th lctronic configurations of six diffrnt atoms. A B C D E F (a) You may us th Priodic Tabl on pag 2 to hlp you answr this qustion. Answr ach part by writing on of th lttrs A, B, C,

More information

1 Isoparametric Concept

1 Isoparametric Concept UNIVERSITY OF CALIFORNIA BERKELEY Dpartmnt of Civil Enginring Spring 06 Structural Enginring, Mchanics and Matrials Profssor: S. Govindj Nots on D isoparamtric lmnts Isoparamtric Concpt Th isoparamtric

More information

Full Waveform Inversion Using an Energy-Based Objective Function with Efficient Calculation of the Gradient

Full Waveform Inversion Using an Energy-Based Objective Function with Efficient Calculation of the Gradient Full Wavform Invrsion Using an Enrgy-Basd Objctiv Function with Efficint Calculation of th Gradint Itm yp Confrnc Papr Authors Choi, Yun Sok; Alkhalifah, ariq Ali Citation Choi Y, Alkhalifah (217) Full

More information

Differential Equations

Differential Equations Prfac Hr ar m onlin nots for m diffrntial quations cours that I tach hr at Lamar Univrsit. Dspit th fact that ths ar m class nots, th should b accssibl to anon wanting to larn how to solv diffrntial quations

More information

Collisions between electrons and ions

Collisions between electrons and ions DRAFT 1 Collisions btwn lctrons and ions Flix I. Parra Rudolf Pirls Cntr for Thortical Physics, Unirsity of Oxford, Oxford OX1 NP, UK This rsion is of 8 May 217 1. Introduction Th Fokkr-Planck collision

More information

MAE4700/5700 Finite Element Analysis for Mechanical and Aerospace Design

MAE4700/5700 Finite Element Analysis for Mechanical and Aerospace Design MAE4700/5700 Finit Elmnt Analysis for Mchanical and Arospac Dsign Cornll Univrsity, Fall 2009 Nicholas Zabaras Matrials Procss Dsign and Control Laboratory Sibly School of Mchanical and Arospac Enginring

More information

Electrochemical Energy Systems Spring 2014 MIT, M. Z. Bazant. Midterm Exam

Electrochemical Energy Systems Spring 2014 MIT, M. Z. Bazant. Midterm Exam 10.66 Elctrochmical Enrgy Systms Spring 014 MIT, M. Z. Bazant Midtrm Exam Instructions. This is a tak-hom, opn-book xam du in Lctur. Lat xams will not b accptd. You may consult any books, handouts, or

More information

hep-lat/ Dec 93

hep-lat/ Dec 93 GLUON VERSUS MESON EXCHANGE IN HADRON-HADRON SYSTEMS ON THE LATTICE 1 H. MARKUM, K. RABITSCH, W. SAKULER Institut fur Krnphysik, Tchnisch Univrsitat Win A-100 Vinna, Austria hp-lat/931059 15 Dc 93 Th intraction

More information

MATH 319, WEEK 15: The Fundamental Matrix, Non-Homogeneous Systems of Differential Equations

MATH 319, WEEK 15: The Fundamental Matrix, Non-Homogeneous Systems of Differential Equations MATH 39, WEEK 5: Th Fundamntal Matrix, Non-Homognous Systms of Diffrntial Equations Fundamntal Matrics Considr th problm of dtrmining th particular solution for an nsmbl of initial conditions For instanc,

More information

Solution: APPM 1360 Final (150 pts) Spring (60 pts total) The following parts are not related, justify your answers:

Solution: APPM 1360 Final (150 pts) Spring (60 pts total) The following parts are not related, justify your answers: APPM 6 Final 5 pts) Spring 4. 6 pts total) Th following parts ar not rlatd, justify your answrs: a) Considr th curv rprsntd by th paramtric quations, t and y t + for t. i) 6 pts) Writ down th corrsponding

More information

4.2 Design of Sections for Flexure

4.2 Design of Sections for Flexure 4. Dsign of Sctions for Flxur This sction covrs th following topics Prliminary Dsign Final Dsign for Typ 1 Mmbrs Spcial Cas Calculation of Momnt Dmand For simply supportd prstrssd bams, th maximum momnt

More information

Function Spaces. a x 3. (Letting x = 1 =)) a(0) + b + c (1) = 0. Row reducing the matrix. b 1. e 4 3. e 9. >: (x = 1 =)) a(0) + b + c (1) = 0

Function Spaces. a x 3. (Letting x = 1 =)) a(0) + b + c (1) = 0. Row reducing the matrix. b 1. e 4 3. e 9. >: (x = 1 =)) a(0) + b + c (1) = 0 unction Spacs Prrquisit: Sction 4.7, Coordinatization n this sction, w apply th tchniqus of Chaptr 4 to vctor spacs whos lmnts ar functions. Th vctor spacs P n and P ar familiar xampls of such spacs. Othr

More information

Abstract Interpretation. Lecture 5. Profs. Aiken, Barrett & Dill CS 357 Lecture 5 1

Abstract Interpretation. Lecture 5. Profs. Aiken, Barrett & Dill CS 357 Lecture 5 1 Abstract Intrprtation 1 History On brakthrough papr Cousot & Cousot 77 (?) Inspird by Dataflow analysis Dnotational smantics Enthusiastically mbracd by th community At last th functional community... At

More information

Toward the understanding of QCD phase structure at finite temperature and density

Toward the understanding of QCD phase structure at finite temperature and density Toward th undrstanding o QCD phas structur at init tmpratur and dnsity Shinji Ejiri iigata Univrsity HOT-QCD collaboration S. Ejiri1 S. Aoki T. Hatsuda3 K. Kanaya Y. akagawa1 H. Ohno H. Saito and T. Umda5

More information

ME 321 Kinematics and Dynamics of Machines S. Lambert Winter 2002

ME 321 Kinematics and Dynamics of Machines S. Lambert Winter 2002 3.4 Forc Analysis of Linkas An undrstandin of forc analysis of linkas is rquird to: Dtrmin th raction forcs on pins, tc. as a consqunc of a spcifid motion (don t undrstimat th sinificanc of dynamic or

More information

Basic Polyhedral theory

Basic Polyhedral theory Basic Polyhdral thory Th st P = { A b} is calld a polyhdron. Lmma 1. Eithr th systm A = b, b 0, 0 has a solution or thr is a vctorπ such that π A 0, πb < 0 Thr cass, if solution in top row dos not ist

More information

Middle East Technical University Department of Mechanical Engineering ME 413 Introduction to Finite Element Analysis

Middle East Technical University Department of Mechanical Engineering ME 413 Introduction to Finite Element Analysis Middl East Tchnical Univrsity Dpartmnt of Mchanical Enginring ME Introduction to Finit Elmnt Analysis Chaptr 5 Two-Dimnsional Formulation Ths nots ar prpard by Dr. Cünyt Srt http://www.m.mtu.du.tr/popl/cunyt

More information

Electrochemistry L E O

Electrochemistry L E O Rmmbr from CHM151 A rdox raction in on in which lctrons ar transfrrd lctrochmistry L O Rduction os lctrons xidation G R ain lctrons duction W can dtrmin which lmnt is oxidizd or rducd by assigning oxidation

More information

Lecture 2: Discrete-Time Signals & Systems. Reza Mohammadkhani, Digital Signal Processing, 2015 University of Kurdistan eng.uok.ac.

Lecture 2: Discrete-Time Signals & Systems. Reza Mohammadkhani, Digital Signal Processing, 2015 University of Kurdistan eng.uok.ac. Lctur 2: Discrt-Tim Signals & Systms Rza Mohammadkhani, Digital Signal Procssing, 2015 Univrsity of Kurdistan ng.uok.ac.ir/mohammadkhani 1 Signal Dfinition and Exampls 2 Signal: any physical quantity that

More information

Electromagnetic scattering. Graduate Course Electrical Engineering (Communications) 1 st Semester, Sharif University of Technology

Electromagnetic scattering. Graduate Course Electrical Engineering (Communications) 1 st Semester, Sharif University of Technology Elctromagntic scattring Graduat Cours Elctrical Enginring (Communications) 1 st Smstr, 1388-1389 Sharif Univrsity of Tchnology Contnts of lctur 8 Contnts of lctur 8: Scattring from small dilctric objcts

More information

There is an arbitrary overall complex phase that could be added to A, but since this makes no difference we set it to zero and choose A real.

There is an arbitrary overall complex phase that could be added to A, but since this makes no difference we set it to zero and choose A real. Midtrm #, Physics 37A, Spring 07. Writ your rsponss blow or on xtra pags. Show your work, and tak car to xplain what you ar doing; partial crdit will b givn for incomplt answrs that dmonstrat som concptual

More information

2. Background Material

2. Background Material S. Blair Sptmbr 3, 003 4. Background Matrial Th rst of this cours dals with th gnration, modulation, propagation, and ction of optical radiation. As such, bic background in lctromagntics and optics nds

More information

Utilizing exact and Monte Carlo methods to investigate properties of the Blume Capel Model applied to a nine site lattice.

Utilizing exact and Monte Carlo methods to investigate properties of the Blume Capel Model applied to a nine site lattice. Utilizing xat and Mont Carlo mthods to invstigat proprtis of th Blum Capl Modl applid to a nin sit latti Nik Franios Writing various xat and Mont Carlo omputr algorithms in C languag, I usd th Blum Capl

More information

Classical Magnetic Dipole

Classical Magnetic Dipole Lctur 18 1 Classical Magntic Dipol In gnral, a particl of mass m and charg q (not ncssarily a point charg), w hav q g L m whr g is calld th gyromagntic ratio, which accounts for th ffcts of non-point charg

More information

Estimation of apparent fraction defective: A mathematical approach

Estimation of apparent fraction defective: A mathematical approach Availabl onlin at www.plagiarsarchlibrary.com Plagia Rsarch Library Advancs in Applid Scinc Rsarch, 011, (): 84-89 ISSN: 0976-8610 CODEN (USA): AASRFC Estimation of apparnt fraction dfctiv: A mathmatical

More information

The Matrix Exponential

The Matrix Exponential Th Matrix Exponntial (with xrciss) by Dan Klain Vrsion 28928 Corrctions and commnts ar wlcom Th Matrix Exponntial For ach n n complx matrix A, dfin th xponntial of A to b th matrix () A A k I + A + k!

More information

EXST Regression Techniques Page 1

EXST Regression Techniques Page 1 EXST704 - Rgrssion Tchniqus Pag 1 Masurmnt rrors in X W hav assumd that all variation is in Y. Masurmnt rror in this variabl will not ffct th rsults, as long as thy ar uncorrlatd and unbiasd, sinc thy

More information

Why is a E&M nature of light not sufficient to explain experiments?

Why is a E&M nature of light not sufficient to explain experiments? 1 Th wird world of photons Why is a E&M natur of light not sufficint to xplain xprimnts? Do photons xist? Som quantum proprtis of photons 2 Black body radiation Stfan s law: Enrgy/ ara/ tim = Win s displacmnt

More information

Higher order derivatives

Higher order derivatives Robrto s Nots on Diffrntial Calculus Chaptr 4: Basic diffrntiation ruls Sction 7 Highr ordr drivativs What you nd to know alrady: Basic diffrntiation ruls. What you can larn hr: How to rpat th procss of

More information

CHAPTER 16 HW: CONJUGATED SYSTEMS

CHAPTER 16 HW: CONJUGATED SYSTEMS APTER 6 W: JUGATED SYSTEMS NAMING PLYENES. Giv th IUPA nam for ach compound, including cis/trans or E/Z dsignations whr ndd. ompound no E/Z trans or E 2 3 4 3 Nam trans-2-mthyl-2,4-hxadin 2-mthoxy-,3-cyclohptadin

More information

ELECTRON-MUON SCATTERING

ELECTRON-MUON SCATTERING ELECTRON-MUON SCATTERING ABSTRACT Th lctron charg is considrd to b distributd or xtndd in spac. Th diffrntial of th lctron charg is st qual to a function of lctron charg coordinats multiplid by a four-dimnsional

More information

A central nucleus. Protons have a positive charge Electrons have a negative charge

A central nucleus. Protons have a positive charge Electrons have a negative charge Atomic Structur Lss than ninty yars ago scintists blivd that atoms wr tiny solid sphrs lik minut snookr balls. Sinc thn it has bn discovrd that atoms ar not compltly solid but hav innr and outr parts.

More information

Human vision is determined based on information theory:

Human vision is determined based on information theory: Human vision is dtrmind basd on information thory: Supplmntary Information Alfonso Dlgado-Bonal,2 and F. Javir Martn Torrs,3 [] Instituto Andaluz d Cincias d la Tirra CSIC-UGR, Avda. d Las Palmras n 4,

More information