CE 530 Molecular Simulation
|
|
- Eric Wiggins
- 5 years ago
- Views:
Transcription
1 CE 530 Molecular Simulatio Lecture 5 Log-rage forces ad Ewald sum David A. Kofke Departmet of Chemical Egieerig SUNY Buffalo kofke@eg.buffalo.edu
2 Review Itermolecular forces arise from quatum mechaics too complex to iclude i legthy simulatios of bulk phases Empirical forms give simple formulas to approximate behavior itramolecular forms: bed, stretch, torsio itermolecular: va der Waals, electrostatics, polarizatio Ulike-atom iteractios weak lik i quatitative work
3 3 Trucatig the Potetial Bulk system modeled via periodic boudary coditio ot feasible to iclude iteractios with all images must trucate potetial at half the box legth (at most) to have all separatios treated cosistetly Cotributios from distat separatios may be importat These two are same distace from cetral atom, yet: Black atom iteracts Gree atom does ot These two are earest images for cetral atom Oly iteractios cosidered
4 Trucatig the Potetial 4 Potetial trucatio itroduces discotiuity Correspods to a ifiite force Problematic for MD simulatios ruis eergy coservatio Shifted potetials Removes ifiite force Still discotiuity i force Shifted-force potetials Routiely used i MD ur () ur ( c) r r () r = 0 r > r For quatitative work eed to re-itroduce log-rage iteractios u s du ur () ur ( ) ( r r) r r usf () r = dr 0 r > r c c c c c c
5 Leard-Joes example r c =.5σ Trucatig the Potetial 5 u s () r ur () ur ( c) r r = 0 r > r c c du ur () ur ( ) ( r r) r r usf () r = dr 0 r > r c c c c u(r) f(r) u(r) u-shift u(r) f-shift f(r) f-shift..6.0 Separatio, r/σ Separatio, r/σ.4.5
6 6 Radial Distributio Fuctio Radial distributio fuctio, g(r) key quatity i statistical mechaics quatifies correlatio betwee atom pairs Defiitio gr () = ρ() rdr ρ id () rdr Number of atoms at r i actual system dr 4 3 Hard-sphere g(r) Low desity High desity Number of atoms at r for ideal gas id N ρ () rdr = dr V Here s a applet that computes g(r)
7 8 Radial Distributio Fuctio. Java Code public class MeterRDF exteds MeterFuctio /** * Computes RDF for the curret cofiguratio */ public double[] curretvalue() { iterator.reset(); //prepare iterator of atom pairs for(it i=0; i<poits; i++) {y[i] = 0.0;} //zero histogram while(iterator.hasnext()) { //loop over all pairs i phase double r = Math.sqrt(iterator.ext().r()); //get pair separatio if(r < xmax) { it idex = (it)(r/delr); //determie histogram idex y[idex]+=; //add oce for each atom } } it = phase.atomcout(); //compute ormalizatio: divide by double orm = */phase.volume(); //, ad desity*(volume of shell) for(it i=0; i<poits; i++) {y[i] /= (orm*vshell[i]);} retur y; }
8 9 Simple Log-Rage Correctio Approximate distat iteractios by assumig uiform distributio beyod cutoff: g(r) = r > r cut Correctios to thermodyamic properties Iteral eergy Virial N Ulrc = ρ u()4 r πr dr cut Chemical potetial r du Plrc = ρ r 4πr dr 6 dr r cut U µ lrc = ρ ur ()4πrdr= N r cut lrc U Expressio for Leard-Joes model P LJ lrc LJ lrc σ σ πnρσ ε 3 r c r c 8 = 9 3 πρ σ ε σ 3 σ = 9 r c r c For r c /σ =.5, these are about 5-0% of the total values
9 Coulombic Log-Rage Correctio Coulombic iteractios must be treated specially very log rage /r form does ot die off as quickly as volume grows r c fiite oly because + ad cotributios cacel Methods Full lattice sum Here is a applet demostratig direct approach Ewald sum 4 r π rdr = Treat surroudigs as dielectric cotiuum Leard-Joes Coulomb 4 0
10 Aside: Fourier Series Cosider periodic fuctio o - L/, +L/ A Fourier series provides a equivalet represetatio of the fuctio f (x) = a 0 + ( a cosx + b si x) = The coefficiets are + L / / + L / / a = f ( x)cos( π x / L) dx b = f ( x)si( π x / L) dx f(x) Oe period L/ x +L/
11 Fourier Series Example f(x) is a square wave f(x) + L / / a = f ( x)cos( π x / L) dx 0 + L / πx L dx / 0 = cos( / ) cos( πx / L) dx = 0 + L / / b = f ( x)si( π x / L) dx = 4 π odd 0 eve x
12 Fourier Series Example 3 f(x) is a square wave = f(x) + L / / a = f ( x)cos( π x / L) dx 0 + L / πx L dx / 0 = cos( / ) cos( πx / L) dx = 0 + L / / b = f ( x)si( π x / L) dx = 4 π odd 0 eve x
13 Fourier Series Example 4 f(x) is a square wave =, 3 f(x) + L / / a = f ( x)cos( π x / L) dx 0 + L / πx L dx / 0 = cos( / ) cos( πx / L) dx = 0 + L / / b = f ( x)si( π x / L) dx = 4 π odd 0 eve x
14 Fourier Series Example 5 f(x) is a square wave =, 3, 5 f(x) + L / / a = f ( x)cos( π x / L) dx 0 + L / πx L dx / 0 = cos( / ) cos( πx / L) dx = 0 + L / / b = f ( x)si( π x / L) dx = 4 π odd 0 eve x
15 Fourier Series Example 6 f(x) is a square wave =, 3, 5, 7 f(x) + L / / a = f ( x)cos( π x / L) dx 0 + L / πx L dx / 0 = cos( / ) cos( πx / L) dx = 0 + L / / b = f ( x)si( π x / L) dx = 4 π odd 0 eve x
16 7 Fourier Represetatio The set of Fourier-space coefficiets b cotai complete iformatio about the fuctio f (x) = a 0 + Although f(x) is periodic to ifiity, b is o-egligible over oly a fiite rage Sometimes the Fourier represetatio is more coveiet to use. = ( ) a cosx + b si x b
17 8 Covergece of Fourier Sum If f(x) = si(πkx/l), trasform is simple b = for = k Coverges very quickly! b = 0 otherwise f(x) b x
18 9 Observatios o Fourier Sum Smooth fuctios f(x) require few coefficiets b Sharp fuctios (square wave) require more coefficiets Large- coefficiets describe high-frequecy behavior of f(x) large = short wavelegth Small- coefficiets describe low-frequecy behavior small = log wavelegth e.g., = 0 coefficiet is simple average of f(x)
19 Fourier Trasform As L icreases, f(x) becomes less periodic Fourier trasform arrives i limit of L à Compact form obtaied with expoetial form of cos/si f (x) = a 0 + a = π b = π + L/ L/ + L/ L/ e = iθ ( + b ) a cos πx L f (x)cos(π x / L) dx f (x)si(π x / L) dx Useful relatios derivative = cosθ + isiθ covolutio si πx L iverse forward + ˆ π ikx f( x) = f( k) e dk + ˆ( ) ( ) π ikx f k = f x e dx a = real part of trasform b = imagiary part f ( x) ( k) = ( πik) f( k) f() t g( x t) dt ( ) ˆ k = f( k) gˆ ( k) ( m) m ˆ 0
20 Fourier Trasform Example Gaussia α x f( x) = exp α π Trasform is also a Gaussia! f(x) / ˆ( α k f k ) = exp π α x Width of trasform is reciprocal of width of fuctio k-space is reciprocal space sharp f(x) requires more values of F(k) for good represetatio δ(x-x o ) trasforms ito a sie/cosie wave of frequecy x o : ˆ( ) o δ k = e πikx F(k) x k
21 Fourier Trasform Relevace May features of statistical-mechaical systems are described i k-space structure trasport behavior electrostatics This descriptio focuses o the correlatios show over a particular legth scale (depedig o k) Macroscopic observables are recovered i the k à 0 limit Correspodig treatmet is applied i the time/frequecy domais
22 Review of Basic Electrostatics 3 Force betwee charges I terms of electric field Static electric field satisfies Er () = 4 πρ() r Er () = 0 Charge desity ρ(r) for poit charge q : Electrostatic potetial zero curl implies E ca be writte potetial eergy of charge q at r, relative to positio at ifiity u() r = qφ() r Poisso s equatio φ = 4πρ qq F= r r ˆ Fr () = q Er () ρ() r = q δ() r Er () = φ() r Mass aalogy: u(z) = m gz = mφ(z)
23 4 Ewald Sum We wat to sum the iteractio eergy of each charge i the cetral volume with all images of the other charges express i terms of electostatic potetial U q = qφ i ri charge i i cetral volume ( ) the charge desity creatig the potetial is ρ(r) = q δ (r r ),image vectors i = q δ r (r + L) ( ) this is a periodic fuctio (of period L), but it is very sharp Fourier represetatio would ever coverge
24 5 Ewald Sum: Fourier. Compute field istead by smearig all the charges 3/ ρ() r = q( α/ π) exp α r ( r + L) Electrostatic potetial via Poisso equatio direct space form reciprocal space iclude = 0 φ() r = 4 πρ() r k φ( k) = 4 πρ( k) Discrete Fourier trasform the charge desity V V ρ( k) = de r ρ( r) = V qe ikr ikr k e /4α Large α takes ρ back to δ fuctio a = L b = L + L/ / + L/ / f (x)cos(π x / L) dx f (x)si(π x / L) dx
25 6 Ewald Sum. Fourier. Use Poisso s equatio for electrostatic potetial 4π φ( k) = ρ( k) k Ivert trasform to recover real-space potetial φ() r = φ( k) e k 0 = V k 0 ikr 4π q e k ik ( r r ) k /4α e f (x) = a 0 + ( a cosx+b si x) i priciple requires sum over ifiite umber of wave vectors k but reciprocal Gaussia goes to zero quickly if α is small (broad Gaussia, large smearig of charge) =
26 Ewald Sum. Fourier 3. 7 The electrostatic eergy ca ow be obtaied for poit charges i potetial of smeared charges U q = q iφ ri i = = Two correctios are eeded self iteractio k /4α i k 0 k i, V k 0 k ( ) 4πV e 4πV e k correct for smearig /4α qq ρ( k) e ik ( r r ) i product of idetical sums i ρ( k) = qe kr V φ(r) = V k 0 4πq k e ik (r r ) e k /4α
27 Ewald Sum. Self Iteractio. 8 I Ewald sum, each poit charge is replaced by smeared Gaussia cetered o that charge this is doe to estimate the electrostatic potetial field x x All poit charges iteract with the resultig field to yield the potetial eergy This meas that the poit charge iteracts with its smeared represetatio We eed to subtract this x
28 Ewald Sum. Self Iteractio. We work i real space to deal with the self term Poisso s equatio for the electrostatic potetial due to a sigle smeared charge 9 φ() r = 4 πρ() r The solutio is q ( α ) φ() r = erf r r ρ() r = ( α/ π) exp α r r 3/ q = r r r r I particular, at r = 0 φ(0) = ( α/ π) q / The self-correctio subtracts this for each charge U self = = α ( ) π q φ(0) q idepedet of cofiguratio
29 30 Ewald Sum. Smearig Correctio. We add the correct field ad subtract the approximate oe to correct for the smearig p G Δ φ () r = φ () r φ () r q q = erf r r r r q = erfc r r ( α r r ) ( α r r ) This field is short raged for large α (arrow Gaussias) ca view as poit charges surrouded by shieldig coutercharge distributio x
30 3 Ewald Sum. Smearig Correctio. Sum iteractio of all charges with field correctio coveiet to stay i real space Usually α is chose so that sum coverges withi cetral image Total Coulomb eergy Δ U = qδφ ( r ) = i i i i qq i erfc r i ( α r i ) U = U ( α) U ( α) +Δ U( α) c q self each term depeds o α, but the sum is idepedet of it if eough lattice vectors are used i the reciprocal- ad real-space sums Here is a applet that demostrate the Ewald method
31 3 Ewald Method. Commets Basic form requires a O(N ) calculatio efficiecy ca be itroduced to reduce to O(N 3/ ) good value of α is 5L, but should check for give applicatio ca be exteded to sum poit dipoles Other methods are i commo use reactio field particle-particle/particle mesh fast multipole
32 33 Summary Cotributios from distat iteractios caot be eglected potetial trucated at o more tha half box legth treat log-rage assumig uiform radial distributio fuctio Coulombic iteractios require explicit summig of images too costly to perform direct sum Ewald method is more efficiet smear charges to approximate electrostatic field simple correctio for self iteractio real-space correctio for smearig
62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +
62. Power series Defiitio 16. (Power series) Give a sequece {c }, the series c x = c 0 + c 1 x + c 2 x 2 + c 3 x 3 + is called a power series i the variable x. The umbers c are called the coefficiets of
More informationMicroscopic Theory of Transport (Fall 2003) Lecture 6 (9/19/03) Static and Short Time Properties of Time Correlation Functions
.03 Microscopic Theory of Trasport (Fall 003) Lecture 6 (9/9/03) Static ad Short Time Properties of Time Correlatio Fuctios Refereces -- Boo ad Yip, Chap There are a umber of properties of time correlatio
More information17 Phonons and conduction electrons in solids (Hiroshi Matsuoka)
7 Phoos ad coductio electros i solids Hiroshi Matsuoa I this chapter we will discuss a miimal microscopic model for phoos i a solid ad a miimal microscopic model for coductio electros i a simple metal.
More informationEwald Summation for Coulomb Interactions in a Periodic Supercell
Ewald Summatio for Coulomb Iteractios i a Periodic Supercell Har Lee ad Wei Cai Departmet of Mechaical Egieerig, Staford Uiversity, CA 9435-44 Jauary, 29 Cotets Problem Statemet 2 Charge Distributio Fuctio
More informationChapter 9: Numerical Differentiation
178 Chapter 9: Numerical Differetiatio Numerical Differetiatio Formulatio of equatios for physical problems ofte ivolve derivatives (rate-of-chage quatities, such as velocity ad acceleratio). Numerical
More informationLecture 4 Conformal Mapping and Green s Theorem. 1. Let s try to solve the following problem by separation of variables
Lecture 4 Coformal Mappig ad Gree s Theorem Today s topics. Solvig electrostatic problems cotiued. Why separatio of variables does t always work 3. Coformal mappig 4. Gree s theorem The failure of separatio
More information1.3 Convergence Theorems of Fourier Series. k k k k. N N k 1. With this in mind, we state (without proof) the convergence of Fourier series.
.3 Covergece Theorems of Fourier Series I this sectio, we preset the covergece of Fourier series. A ifiite sum is, by defiitio, a limit of partial sums, that is, a cos( kx) b si( kx) lim a cos( kx) b si(
More information3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,
3. Z Trasform Referece: Etire Chapter 3 of text. Recall that the Fourier trasform (FT) of a DT sigal x [ ] is ω ( ) [ ] X e = j jω k = xe I order for the FT to exist i the fiite magitude sese, S = x [
More informationL = n i, i=1. dp p n 1
Exchageable sequeces ad probabilities for probabilities 1996; modified 98 5 21 to add material o mutual iformatio; modified 98 7 21 to add Heath-Sudderth proof of de Fietti represetatio; modified 99 11
More informationFilter banks. Separately, the lowpass and highpass filters are not invertible. removes the highest frequency 1/ 2and
Filter bas Separately, the lowpass ad highpass filters are ot ivertible T removes the highest frequecy / ad removes the lowest frequecy Together these filters separate the sigal ito low-frequecy ad high-frequecy
More informationStatistical Pattern Recognition
Statistical Patter Recogitio Classificatio: No-Parametric Modelig Hamid R. Rabiee Jafar Muhammadi Sprig 2014 http://ce.sharif.edu/courses/92-93/2/ce725-2/ Ageda Parametric Modelig No-Parametric Modelig
More informationSignal Processing in Mechatronics. Lecture 3, Convolution, Fourier Series and Fourier Transform
Sigal Processig i Mechatroics Summer semester, 1 Lecture 3, Covolutio, Fourier Series ad Fourier rasform Dr. Zhu K.P. AIS, UM 1 1. Covolutio Covolutio Descriptio of LI Systems he mai premise is that the
More informationPhysics 324, Fall Dirac Notation. These notes were produced by David Kaplan for Phys. 324 in Autumn 2001.
Physics 324, Fall 2002 Dirac Notatio These otes were produced by David Kapla for Phys. 324 i Autum 2001. 1 Vectors 1.1 Ier product Recall from liear algebra: we ca represet a vector V as a colum vector;
More informationMatsubara-Green s Functions
Matsubara-Gree s Fuctios Time Orderig : Cosider the followig operator If H = H the we ca trivially factorise this as, E(s = e s(h+ E(s = e sh e s I geeral this is ot true. However for practical applicatio
More information4. Partial Sums and the Central Limit Theorem
1 of 10 7/16/2009 6:05 AM Virtual Laboratories > 6. Radom Samples > 1 2 3 4 5 6 7 4. Partial Sums ad the Cetral Limit Theorem The cetral limit theorem ad the law of large umbers are the two fudametal theorems
More informationLesson 10: Limits and Continuity
www.scimsacademy.com Lesso 10: Limits ad Cotiuity SCIMS Academy 1 Limit of a fuctio The cocept of limit of a fuctio is cetral to all other cocepts i calculus (like cotiuity, derivative, defiite itegrals
More informationThe z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j
The -Trasform 7. Itroductio Geeralie the complex siusoidal represetatio offered by DTFT to a represetatio of complex expoetial sigals. Obtai more geeral characteristics for discrete-time LTI systems. 7.
More informationChapter 4. Fourier Series
Chapter 4. Fourier Series At this poit we are ready to ow cosider the caoical equatios. Cosider, for eample the heat equatio u t = u, < (4.) subject to u(, ) = si, u(, t) = u(, t) =. (4.) Here,
More informationChapter 7: The z-transform. Chih-Wei Liu
Chapter 7: The -Trasform Chih-Wei Liu Outlie Itroductio The -Trasform Properties of the Regio of Covergece Properties of the -Trasform Iversio of the -Trasform The Trasfer Fuctio Causality ad Stability
More informationDiscrete Mathematics for CS Spring 2008 David Wagner Note 22
CS 70 Discrete Mathematics for CS Sprig 2008 David Wager Note 22 I.I.D. Radom Variables Estimatig the bias of a coi Questio: We wat to estimate the proportio p of Democrats i the US populatio, by takig
More informationCHAPTER 10 INFINITE SEQUENCES AND SERIES
CHAPTER 10 INFINITE SEQUENCES AND SERIES 10.1 Sequeces 10.2 Ifiite Series 10.3 The Itegral Tests 10.4 Compariso Tests 10.5 The Ratio ad Root Tests 10.6 Alteratig Series: Absolute ad Coditioal Covergece
More information7.1 Convergence of sequences of random variables
Chapter 7 Limit Theorems Throughout this sectio we will assume a probability space (, F, P), i which is defied a ifiite sequece of radom variables (X ) ad a radom variable X. The fact that for every ifiite
More informationLecture 5-2: Polytropes. Literature: MWW chapter 19
Lecture 5-2: Polytropes Literature: MWW chapter 9!" Preamble The 4 equatios of stellar structure divide ito two groups: Mass ad mometum describig the mechaical structure ad thermal equilibrium ad eergy
More informationApproximations and more PMFs and PDFs
Approximatios ad more PMFs ad PDFs Saad Meimeh 1 Approximatio of biomial with Poisso Cosider the biomial distributio ( b(k,,p = p k (1 p k, k λ: k Assume that is large, ad p is small, but p λ at the limit.
More informationInfinite Sequences and Series
Chapter 6 Ifiite Sequeces ad Series 6.1 Ifiite Sequeces 6.1.1 Elemetary Cocepts Simply speakig, a sequece is a ordered list of umbers writte: {a 1, a 2, a 3,...a, a +1,...} where the elemets a i represet
More informationQuestion 1: The magnetic case
September 6, 018 Corell Uiversity, Departmet of Physics PHYS 337, Advace E&M, HW # 4, due: 9/19/018, 11:15 AM Questio 1: The magetic case I class, we skipped over some details, so here you are asked to
More informationPH 411/511 ECE B(k) Sin k (x) dk (1)
Fall-26 PH 4/5 ECE 598 A. La Rosa Homework-2 Due -3-26 The Homework is iteded to gai a uderstadig o the Heiseberg priciple, based o a compariso betwee the width of a pulse ad the width of its spectral
More informationEE / EEE SAMPLE STUDY MATERIAL. GATE, IES & PSUs Signal System. Electrical Engineering. Postal Correspondence Course
Sigal-EE Postal Correspodece Course 1 SAMPLE STUDY MATERIAL Electrical Egieerig EE / EEE Postal Correspodece Course GATE, IES & PSUs Sigal System Sigal-EE Postal Correspodece Course CONTENTS 1. SIGNAL
More informationSECTION 2 Electrostatics
SECTION Electrostatics This sectio, based o Chapter of Griffiths, covers effects of electric fields ad forces i static (timeidepedet) situatios. The topics are: Electric field Gauss s Law Electric potetial
More informationx a x a Lecture 2 Series (See Chapter 1 in Boas)
Lecture Series (See Chapter i Boas) A basic ad very powerful (if pedestria, recall we are lazy AD smart) way to solve ay differetial (or itegral) equatio is via a series expasio of the correspodig solutio
More information(A sequence also can be thought of as the list of function values attained for a function f :ℵ X, where f (n) = x n for n 1.) x 1 x N +k x N +4 x 3
MATH 337 Sequeces Dr. Neal, WKU Let X be a metric space with distace fuctio d. We shall defie the geeral cocept of sequece ad limit i a metric space, the apply the results i particular to some special
More informationPHYC - 505: Statistical Mechanics Homework Assignment 4 Solutions
PHYC - 55: Statistical Mechaics Homewor Assigmet 4 Solutios Due February 5, 14 1. Cosider a ifiite classical chai of idetical masses coupled by earest eighbor sprigs with idetical sprig costats. a Write
More informationIntroduction to Extreme Value Theory Laurens de Haan, ISM Japan, Erasmus University Rotterdam, NL University of Lisbon, PT
Itroductio to Extreme Value Theory Laures de Haa, ISM Japa, 202 Itroductio to Extreme Value Theory Laures de Haa Erasmus Uiversity Rotterdam, NL Uiversity of Lisbo, PT Itroductio to Extreme Value Theory
More informationSequences and Series of Functions
Chapter 6 Sequeces ad Series of Fuctios 6.1. Covergece of a Sequece of Fuctios Poitwise Covergece. Defiitio 6.1. Let, for each N, fuctio f : A R be defied. If, for each x A, the sequece (f (x)) coverges
More informationTopic 9: Sampling Distributions of Estimators
Topic 9: Samplig Distributios of Estimators Course 003, 2016 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be
More informationTIME-CORRELATION FUNCTIONS
p. 8 TIME-CORRELATION FUNCTIONS Time-correlatio fuctios are a effective way of represetig the dyamics of a system. They provide a statistical descriptio of the time-evolutio of a variable for a esemble
More informationPhys 6303 Final Exam Solutions December 19, 2012
Phys 633 Fial Exam s December 19, 212 You may NOT use ay book or otes other tha supplied with this test. You will have 3 hours to fiish. DO YOUR OWN WORK. Express your aswers clearly ad cocisely so that
More informationLecture #5: Begin Quantum Mechanics: Free Particle and Particle in a 1D Box
561 Fall 013 Lecture #5 page 1 Last time: Lecture #5: Begi Quatum Mechaics: Free Particle ad Particle i a 1D Box u 1 u 1-D Wave equatio = x v t * u(x,t): displacemets as fuctio of x,t * d -order: solutio
More information8. IRREVERSIBLE AND RANDOM PROCESSES Concepts and Definitions
8. IRREVERSIBLE ND RNDOM PROCESSES 8.1. Cocepts ad Defiitios I codesed phases, itermolecular iteractios ad collective motios act to modify the state of a molecule i a time-depedet fashio. Liquids, polymers,
More informationAssignment 2 Solutions SOLUTION. ϕ 1 Â = 3 ϕ 1 4i ϕ 2. The other case can be dealt with in a similar way. { ϕ 2 Â} χ = { 4i ϕ 1 3 ϕ 2 } χ.
PHYSICS 34 QUANTUM PHYSICS II (25) Assigmet 2 Solutios 1. With respect to a pair of orthoormal vectors ϕ 1 ad ϕ 2 that spa the Hilbert space H of a certai system, the operator  is defied by its actio
More informationChapter 7 z-transform
Chapter 7 -Trasform Itroductio Trasform Uilateral Trasform Properties Uilateral Trasform Iversio of Uilateral Trasform Determiig the Frequecy Respose from Poles ad Zeros Itroductio Role i Discrete-Time
More informationRandom Variables, Sampling and Estimation
Chapter 1 Radom Variables, Samplig ad Estimatio 1.1 Itroductio This chapter will cover the most importat basic statistical theory you eed i order to uderstad the ecoometric material that will be comig
More informationLecture 18: Sampling distributions
Lecture 18: Samplig distributios I may applicatios, the populatio is oe or several ormal distributios (or approximately). We ow study properties of some importat statistics based o a radom sample from
More informationPH 411/511 ECE B(k) Sin k (x) dk (1)
Fall-27 PH 4/5 ECE 598 A. La Rosa Homework-3 Due -7-27 The Homework is iteded to gai a uderstadig o the Heiseberg priciple, based o a compariso betwee the width of a pulse ad the width of its spectral
More informationExpectation and Variance of a random variable
Chapter 11 Expectatio ad Variace of a radom variable The aim of this lecture is to defie ad itroduce mathematical Expectatio ad variace of a fuctio of discrete & cotiuous radom variables ad the distributio
More informationJacob Hays Amit Pillay James DeFelice 4.1, 4.2, 4.3
No-Parametric Techiques Jacob Hays Amit Pillay James DeFelice 4.1, 4.2, 4.3 Parametric vs. No-Parametric Parametric Based o Fuctios (e.g Normal Distributio) Uimodal Oly oe peak Ulikely real data cofies
More informationSalmon: Lectures on partial differential equations. 3. First-order linear equations as the limiting case of second-order equations
3. First-order liear equatios as the limitig case of secod-order equatios We cosider the advectio-diffusio equatio (1) v = 2 o a bouded domai, with boudary coditios of prescribed. The coefficiets ( ) (2)
More informationAnalysis of composites with multiple rigid-line reinforcements by the BEM
Aalysis of composites with multiple rigid-lie reiforcemets by the BEM Piotr Fedeliski* Departmet of Stregth of Materials ad Computatioal Mechaics, Silesia Uiversity of Techology ul. Koarskiego 18A, 44-100
More informationLecture 8: Solving the Heat, Laplace and Wave equations using finite difference methods
Itroductory lecture otes o Partial Differetial Equatios - c Athoy Peirce. Not to be copied, used, or revised without explicit writte permissio from the copyright ower. 1 Lecture 8: Solvig the Heat, Laplace
More informationThe time evolution of the state of a quantum system is described by the time-dependent Schrödinger equation (TDSE): ( ) ( ) 2m "2 + V ( r,t) (1.
Adrei Tokmakoff, MIT Departmet of Chemistry, 2/13/2007 1-1 574 TIME-DEPENDENT QUANTUM MECHANICS 1 INTRODUCTION 11 Time-evolutio for time-idepedet Hamiltoias The time evolutio of the state of a quatum system
More informationTopic 9: Sampling Distributions of Estimators
Topic 9: Samplig Distributios of Estimators Course 003, 2018 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be
More information1 of 7 7/16/2009 6:06 AM Virtual Laboratories > 6. Radom Samples > 1 2 3 4 5 6 7 6. Order Statistics Defiitios Suppose agai that we have a basic radom experimet, ad that X is a real-valued radom variable
More informationNonequilibrium Excess Carriers in Semiconductors
Lecture 8 Semicoductor Physics VI Noequilibrium Excess Carriers i Semicoductors Noequilibrium coditios. Excess electros i the coductio bad ad excess holes i the valece bad Ambiolar trasort : Excess electros
More informationRay Optics Theory and Mode Theory. Dr. Mohammad Faisal Dept. of EEE, BUET
Ray Optics Theory ad Mode Theory Dr. Mohammad Faisal Dept. of, BUT Optical Fiber WG For light to be trasmitted through fiber core, i.e., for total iteral reflectio i medium, > Ray Theory Trasmissio Ray
More informationTopic 9: Sampling Distributions of Estimators
Topic 9: Samplig Distributios of Estimators Course 003, 2018 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be
More information7-1. Chapter 4. Part I. Sampling Distributions and Confidence Intervals
7-1 Chapter 4 Part I. Samplig Distributios ad Cofidece Itervals 1 7- Sectio 1. Samplig Distributio 7-3 Usig Statistics Statistical Iferece: Predict ad forecast values of populatio parameters... Test hypotheses
More informationFFTs in Graphics and Vision. The Fast Fourier Transform
FFTs i Graphics ad Visio The Fast Fourier Trasform 1 Outlie The FFT Algorithm Applicatios i 1D Multi-Dimesioal FFTs More Applicatios Real FFTs 2 Computatioal Complexity To compute the movig dot-product
More informationPHYS-3301 Lecture 10. Wave Packet Envelope Wave Properties of Matter and Quantum Mechanics I CHAPTER 5. Announcement. Sep.
Aoucemet Course webpage http://www.phys.ttu.edu/~slee/3301/ PHYS-3301 Lecture 10 HW3 (due 10/4) Chapter 5 4, 8, 11, 15, 22, 27, 36, 40, 42 Sep. 27, 2018 Exam 1 (10/4) Chapters 3, 4, & 5 CHAPTER 5 Wave
More information5.61 Fall 2013 Problem Set #3
5.61 Fall 013 Problem Set #3 1. A. McQuarrie, page 10, #3-3. B. McQuarrie, page 10, #3-4. C. McQuarrie, page 18, #4-11.. McQuarrie, pages 11-1, #3-11. 3. A. McQuarrie, page 13, #3-17. B. McQuarrie, page
More informationLecture Chapter 6: Convergence of Random Sequences
ECE5: Aalysis of Radom Sigals Fall 6 Lecture Chapter 6: Covergece of Radom Sequeces Dr Salim El Rouayheb Scribe: Abhay Ashutosh Doel, Qibo Zhag, Peiwe Tia, Pegzhe Wag, Lu Liu Radom sequece Defiitio A ifiite
More informationAppendix: The Laplace Transform
Appedix: The Laplace Trasform The Laplace trasform is a powerful method that ca be used to solve differetial equatio, ad other mathematical problems. Its stregth lies i the fact that it allows the trasformatio
More informationRandom Models. Tusheng Zhang. February 14, 2013
Radom Models Tusheg Zhag February 14, 013 1 Radom Walks Let me describe the model. Radom walks are used to describe the motio of a movig particle (object). Suppose that a particle (object) moves alog the
More informationFIR Filter Design: Part II
EEL335: Discrete-Time Sigals ad Systems. Itroductio I this set of otes, we cosider how we might go about desigig FIR filters with arbitrary frequecy resposes, through compositio of multiple sigle-peak
More informationMath 155 (Lecture 3)
Math 55 (Lecture 3) September 8, I this lecture, we ll cosider the aswer to oe of the most basic coutig problems i combiatorics Questio How may ways are there to choose a -elemet subset of the set {,,,
More informationUnbiased Estimation. February 7-12, 2008
Ubiased Estimatio February 7-2, 2008 We begi with a sample X = (X,..., X ) of radom variables chose accordig to oe of a family of probabilities P θ where θ is elemet from the parameter space Θ. For radom
More informationSignal Processing in Mechatronics
Sigal Processig i Mechatroics Zhu K.P. AIS, UM. Lecture, Brief itroductio to Sigals ad Systems, Review of Liear Algebra ad Sigal Processig Related Mathematics . Brief Itroductio to Sigals What is sigal
More informationMa 530 Infinite Series I
Ma 50 Ifiite Series I Please ote that i additio to the material below this lecture icorporated material from the Visual Calculus web site. The material o sequeces is at Visual Sequeces. (To use this li
More informationCS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 5
CS434a/54a: Patter Recogitio Prof. Olga Veksler Lecture 5 Today Itroductio to parameter estimatio Two methods for parameter estimatio Maimum Likelihood Estimatio Bayesia Estimatio Itroducto Bayesia Decisio
More informationEcon 325/327 Notes on Sample Mean, Sample Proportion, Central Limit Theorem, Chi-square Distribution, Student s t distribution 1.
Eco 325/327 Notes o Sample Mea, Sample Proportio, Cetral Limit Theorem, Chi-square Distributio, Studet s t distributio 1 Sample Mea By Hiro Kasahara We cosider a radom sample from a populatio. Defiitio
More informationTrue Nature of Potential Energy of a Hydrogen Atom
True Nature of Potetial Eergy of a Hydroge Atom Koshu Suto Key words: Bohr Radius, Potetial Eergy, Rest Mass Eergy, Classical Electro Radius PACS codes: 365Sq, 365-w, 33+p Abstract I cosiderig the potetial
More informationradians A function f ( x ) is called periodic if it is defined for all real x and if there is some positive number P such that:
Fourier Series. Graph of y Asix ad y Acos x Amplitude A ; period 36 radias. Harmoics y y six is the first harmoic y y six is the th harmoics 3. Periodic fuctio A fuctio f ( x ) is called periodic if it
More informationMathematical Methods for Physics and Engineering
Mathematical Methods for Physics ad Egieerig Lecture otes Sergei V. Shabaov Departmet of Mathematics, Uiversity of Florida, Gaiesville, FL 326 USA CHAPTER The theory of covergece. Numerical sequeces..
More informationLecture 9: Diffusion, Electrostatics review, and Capacitors. Context
EECS 5 Sprig 4, Lecture 9 Lecture 9: Diffusio, Electrostatics review, ad Capacitors EECS 5 Sprig 4, Lecture 9 Cotext I the last lecture, we looked at the carriers i a eutral semicoductor, ad drift currets
More informationFrequency Domain Filtering
Frequecy Domai Filterig Raga Rodrigo October 19, 2010 Outlie Cotets 1 Itroductio 1 2 Fourier Represetatio of Fiite-Duratio Sequeces: The Discrete Fourier Trasform 1 3 The 2-D Discrete Fourier Trasform
More informationExercises and Problems
HW Chapter 4: Oe-Dimesioal Quatum Mechaics Coceptual Questios 4.. Five. 4.4.. is idepedet of. a b c mu ( E). a b m( ev 5 ev) c m(6 ev ev) Exercises ad Problems 4.. Model: Model the electro as a particle
More informationSolutions to Final Exam Review Problems
. Let f(x) 4+x. Solutios to Fial Exam Review Problems Math 5C, Witer 2007 (a) Fid the Maclauri series for f(x), ad compute its radius of covergece. Solutio. f(x) 4( ( x/4)) ( x/4) ( ) 4 4 + x. Sice the
More informationUniversity of Washington Department of Chemistry Chemistry 453 Winter Quarter 2015
Uiversity of Wasigto Departmet of Cemistry Cemistry 453 Witer Quarter 15 Lecture 14. /11/15 Recommeded Text Readig: Atkis DePaula: 9.1, 9., 9.3 A. Te Equipartitio Priciple & Eergy Quatizatio Te Equipartio
More information7.1 Convergence of sequences of random variables
Chapter 7 Limit theorems Throughout this sectio we will assume a probability space (Ω, F, P), i which is defied a ifiite sequece of radom variables (X ) ad a radom variable X. The fact that for every ifiite
More informationPRELIM PROBLEM SOLUTIONS
PRELIM PROBLEM SOLUTIONS THE GRAD STUDENTS + KEN Cotets. Complex Aalysis Practice Problems 2. 2. Real Aalysis Practice Problems 2. 4 3. Algebra Practice Problems 2. 8. Complex Aalysis Practice Problems
More informationChapter 6 Infinite Series
Chapter 6 Ifiite Series I the previous chapter we cosidered itegrals which were improper i the sese that the iterval of itegratio was ubouded. I this chapter we are goig to discuss a topic which is somewhat
More informationBasics of Probability Theory (for Theory of Computation courses)
Basics of Probability Theory (for Theory of Computatio courses) Oded Goldreich Departmet of Computer Sciece Weizma Istitute of Sciece Rehovot, Israel. oded.goldreich@weizma.ac.il November 24, 2008 Preface.
More informationA Lattice Green Function Introduction. Abstract
August 5, 25 A Lattice Gree Fuctio Itroductio Stefa Hollos Exstrom Laboratories LLC, 662 Nelso Park Dr, Logmot, Colorado 853, USA Abstract We preset a itroductio to lattice Gree fuctios. Electroic address:
More informationDeterministic Model of Multipath Fading for Circular and Parabolic Reflector Patterns
To appear i the Proceedigs of the 5 IEEE outheastco, (Ft. Lauderdale, FL), April 5 Determiistic Model of Multipath Fadig for Circular ad Parabolic Reflector Patters Dwight K. Hutcheso dhutche@clemso.edu
More information2. Fourier Series, Fourier Integrals and Fourier Transforms
Mathematics IV -. Fourier Series, Fourier Itegrals ad Fourier Trasforms The Fourier series are used for the aalysis of the periodic pheomea, which ofte appear i physics ad egieerig. The Fourier itegrals
More informationLecture 2: Poisson Sta*s*cs Probability Density Func*ons Expecta*on and Variance Es*mators
Lecture 2: Poisso Sta*s*cs Probability Desity Fuc*os Expecta*o ad Variace Es*mators Biomial Distribu*o: P (k successes i attempts) =! k!( k)! p k s( p s ) k prob of each success Poisso Distributio Note
More informationHE ATOM & APPROXIMATION METHODS MORE GENERAL VARIATIONAL TREATMENT. Examples:
5.6 4 Lecture #3-4 page HE ATOM & APPROXIMATION METHODS MORE GENERAL VARIATIONAL TREATMENT Do t restrict the wavefuctio to a sigle term! Could be a liear combiatio of several wavefuctios e.g. two terms:
More informationOffice: JILA A709; Phone ;
Office: JILA A709; Phoe 303-49-7841; email: weberjm@jila.colorado.edu Problem Set 5 To be retured before the ed of class o Wedesday, September 3, 015 (give to me i perso or slide uder office door). 1.
More informationWe are mainly going to be concerned with power series in x, such as. (x)} converges - that is, lims N n
Review of Power Series, Power Series Solutios A power series i x - a is a ifiite series of the form c (x a) =c +c (x a)+(x a) +... We also call this a power series cetered at a. Ex. (x+) is cetered at
More informationECE Spring Prof. David R. Jackson ECE Dept. Notes 20
ECE 6341 Sprig 016 Prof. David R. Jackso ECE Dept. Notes 0 1 Spherical Wave Fuctios Cosider solvig ψ + k ψ = 0 i spherical coordiates z φ θ r y x Spherical Wave Fuctios (cot.) I spherical coordiates we
More informationMathematical Statistics - MS
Paper Specific Istructios. The examiatio is of hours duratio. There are a total of 60 questios carryig 00 marks. The etire paper is divided ito three sectios, A, B ad C. All sectios are compulsory. Questios
More information1 Introduction: within and beyond the normal approximation
Tel Aviv Uiversity, 205 Large ad moderate deviatios Itroductio: withi ad beyod the ormal approximatio a Mathematical prelude................ b Physical prelude................... 3 a Mathematical prelude
More informationPrinciple Of Superposition
ecture 5: PREIMINRY CONCEP O RUCUR NYI Priciple Of uperpositio Mathematically, the priciple of superpositio is stated as ( a ) G( a ) G( ) G a a or for a liear structural system, the respose at a give
More informationChapter 10 Partial Differential Equations and Fourier Series
Math-33 Chapter Partial Differetial Equatios November 6, 7 Chapter Partial Differetial Equatios ad Fourier Series Math-33 Chapter Partial Differetial Equatios November 6, 7. Boudary Value Problems for
More informationMATH4822E FOURIER ANALYSIS AND ITS APPLICATIONS
MATH48E FOURIER ANALYSIS AND ITS APPLICATIONS 7.. Cesàro summability. 7. Summability methods Arithmetic meas. The followig idea is due to the Italia geometer Eresto Cesàro (859-96). He shows that eve if
More informationSOLUTIONS: ECE 606 Homework Week 7 Mark Lundstrom Purdue University (revised 3/27/13) e E i E T
SOUIONS: ECE 606 Homework Week 7 Mark udstrom Purdue Uiversity (revised 3/27/13) 1) Cosider a - type semicoductor for which the oly states i the badgap are door levels (i.e. ( E = E D ). Begi with the
More information1 Introduction to reducing variance in Monte Carlo simulations
Copyright c 010 by Karl Sigma 1 Itroductio to reducig variace i Mote Carlo simulatios 11 Review of cofidece itervals for estimatig a mea I statistics, we estimate a ukow mea µ = E(X) of a distributio by
More informationMechanics Physics 151
Mechaics Physics 151 Lecture 4 Cotiuous Systems ad Fields (Chapter 13) What We Did Last Time Built Lagragia formalism for cotiuous system Lagragia L = L dxdydz d L L Lagrage s equatio = dx η, η Derived
More informationStatisticians use the word population to refer the total number of (potential) observations under consideration
6 Samplig Distributios Statisticias use the word populatio to refer the total umber of (potetial) observatios uder cosideratio The populatio is just the set of all possible outcomes i our sample space
More informationhttp://www.xelca.l/articles/ufo_ladigsbaa_houte.aspx imulatio Output aalysis 3/4/06 This lecture Output: A simulatio determies the value of some performace measures, e.g. productio per hour, average queue
More informationThe Riemann Zeta Function
Physics 6A Witer 6 The Riema Zeta Fuctio I this ote, I will sketch some of the mai properties of the Riema zeta fuctio, ζ(x). For x >, we defie ζ(x) =, x >. () x = For x, this sum diverges. However, we
More information