Basics of thermal conductivity calculations
|
|
- Cori Poole
- 6 years ago
- Views:
Transcription
1 Bascs of thermal conductvty calculatons L Ju sometme n Equlbrum Correlaton The bass of calculaton of thermal conductvty and other transport coeffcents usng MD s lnear response theory (Green 1954, Kubo 1957, Mor 1958, Kadanoff and Martn The thermal conductvty ˆκ can be expressed n terms of equlbrum tme autocorrelaton functon of the heat current J q : ˆκ = 1 k B T 2 Ω dτ < J q (J q (τ > (1 Here Ω s the total volume. The nstantaneous heat current s gven by J q = d (E hr (2 dt where E s the mcroscopc energy for each partcle and h s the enthalpy per partcle of the system. If the nteracton s of par potental type, then E = p2 2m V (r j (3 j 1
2 from whch (2 can be wrtten as J q = d (E hr dt = ( (E hv + r [ p F + 1 V j (v v j ] m 2 r j = [(E h 1 (r r j V j ] v 2 r j = (E hv 1 (r r j [ V j (v + v j ] (4 2 r pars On the other hand, when the nteracton s of many-body nature, how to devde the potental energy among the nteractng partcles s not so obvous, and n fact can be qute arbtary. Yet when t s done J q can be derved from the general formalsm: J q = d (E hr dt = ( (E hv + r (v F + E v j r j j = ( (E hv + r ( v E j E + r v j r j r j j = ( (E hv + r j ( E v j r j j (5 (5 provdes an obvous method to calculate the thermal conductvty n a model system, snce all the dynamcal varables n (5 are explctly known n a molecular dynamcs smulaton. Although (1 was derved for a canoncal ensemble, mcrocanoncal ensemble s usually used n MD for t s much easer to mplement, the dfference s only trval snce the phase ponts generated n MD would be n the peak center of densty dstrbuton of the correspondng canoncal ensemble. We wll call ths method equlbrum correlaton method. 2
3 2 Nonequlbrum MD Non-equlbrum molecular dynamcs (NEMD s a formalsm whch bypasses (1, and gve the correct results. It s also based on the lnear response theory, but wth consderable orgnalty and theoretcal nterest n tself, manly, the concept of equvalence of equatons of motons. We know that when a physcal system has Hamltonan H(q, p, the equaton of moton can be readly wrtten down: The Louvlle operator s defned as LA(p, q = so that from Louvlle theorem = q = H/ p (6 p = H/ q ( A q + A p q p ( A H A H q p p q = {A, H} (7 (ṗ, q Tr( (p, q = = Df Dt f Lf (8 t there s f(t = exp(tlf(, where f(t s the densty dstrbuton functonal of phase ponts n Γ space. It can be shown that L s a Hermtan operator. Because the measurement of A(p, q s Ā(t = = Γ Γ dγf(ta(p, q dγ[exp(tlf(]a(p, q (9 3
4 so dā(t dt = = = Γ Γ Γ dγ[l exp(tlf(]a(p, q dγ exp(tlf([ LA(p, q] dγ[exp(tlf(]ȧ(p, q (1 where we have defned the tme dervatve of an arbtary dynamcal varable to be A(p, q = LA(p, q ( A = q + A p q p (11 where the second row conforms both to (7 and the physcal nterpretaton of tme dervatve. From here we can derve the lnear response formulae: suppose the system was orgnally governed by H and s already n thermal equlbrum: f = exp( βh f / t = β exp( βh (L H (12 At t = the perturbaton H was turned on, so H = H + H (13 and L = L + L, f(t = f + f (t (14 The flow equaton (8 remans vald as long as the system s descrbed by a Hamltonan, f t = Lf 4
5 and to frst order n perturbaton strength t can be lnearzed, Multplyng both sde by exp( tl, we have f t = L f + L f (15 [exp( tl f ] t = exp( tl L f whch can be ntergrated to gve f (t = = β = β = β dτ exp((t τl [L f ] dτ exp((t τl f [L H ] dτ exp((t τl f [L H ] dτ exp((t τl f Ḣ (τ (16 For an arbtary dynamcal varable A, the measurement of response to such perturbaton at t > s Ã(t Ā = dγf (ta(p, q Γ = β dτ dγf Ḣ (τ[exp( (t τl A(p, q] = β Ths s the central result of lnear response theory. Γ dτ < A(t τḣ (τ > (17 3 NEMD scheme by Cccott et al As to how (17 can be used to get ˆκ wthout resortng to (1 explctly, let s compare the two equatons: suppose we let A = J q x (p, q n (17, then the rght hand sde of (17 and (1 would look almost smlar, except that one has Ḣ (τ, the other has J q (τ. However, we can 5
6 choose a perturbaton Hamltonan to satsfy Ḣ = λj q x (18 where λ s a very small constant, such that the RHS of (17 and (1 are dentcal except for a constant. Under such a perturbaton, the thermal conductvty ˆκ s smply ( q J ˆκ = lm lm x (t t λ λωt The above formulae assume ˆκ to be an sotropc tensor,.e., ˆκ = κi. When ˆκ sn t sotropc, we can let (19 Ḣ = λj q (t (2 and ( q J j κ j = lm lm (t t λ λωt (21 Whether ˆκ s sotropc or not depend on the local symmetry of the crystal. When the local symmetry s cubc (sold argon or tetrahedral (SC, we can be sure that ˆκ s sotropc from group theory. The problem mmedately follows that how can we choose a H to satsfy (18. However, the answer s obvous once we note that J q s defned by the equaton n (2, so that J q = d (E hr dt H = b (E hr (22 would satsfy (18 automatcally, where b s a small constant. Up to now an deal NEMD scheme had been completely set up. To mplement t we need to put H nto the equaton of moton: r = (p /m (1 + b r p = F (1 + b r ( (E h 1 2 j (r r j V j r b (23 6
7 For sotropc systems, ( b κ = lm Jq (t t b 2 ΩT (24 Ths scheme, however, s only deal and cannot be mplemented on the computer because of the fnte perodcty of the smulaton cell. Imagne the stuaton descrbed by (22 and (23: n order to generate a steady heat current n the system, a constant feld had to be mposed on the whole space,.e., the potental (22 descrbng ths drvng force should be a lnear functon of r (n strct mathematcal sense, ts value rangng from to +. The problem s that on a computer, nstead of a straght lne from to +, we can only generate potental functons of jgjaw shape, wth perodcty of the smulaton cell: dscontnuty would happen on the boundary wth fnte jumps. In other word, there can be no absolute measures of r s that appear n (23. The measures n a peorodc cell can only be of exp(k r, wth k commesurate to the cell, whch means we can only have the fnte k values of the thermal conductvty f we stck to the method. Cccott et al (1978 show that we can study the fnte k thermal conductvty of the system by takng H = λ wth the fnte k heat flux defned by J q k = dentfy A n (17 by J q k k, we have (E hφ, φ = exp(k r (25 ( exp( k r (E hi 1 2 j V j r j p (26 r m J q λβ t k (t k = dτ k < J q k Ω (Jq k (τ > k (27 from where we can get the fnte k autocorrelaton functon of heat flux. And by extrapolatng the fnte k results onto k = we wll arrve at what we want: the macroscopc thermal conductvty. The equaton of moton descrbed by (25 s r = (p /m (1 + λφ (28 ( p = F (1 + λφ λ (29 (E h φ r 1 2 j (φ φ j V j r 7
8 Cccott s scheme s a practcal method to calculate kk. The shortcomng les n the nontrval extrapolaton procedure from fnte k to k =. Another strange feature about (28,29 s that the veloctes r are not smply gven by the momenta p devded by mass. 4 Gllan and Dxon: Equvalence of Equaton of Motons Look at (27 and (28, 29, and let s thnk carefully of where the dffculty les: obvously when we take the lmt k, there must be J q k k n order to satsfy (27. The second term on RHS of (29 ( λ (E h φ 1 r 2 j (φ φ j V j r satsfes ths requrement explctly, snce φ / r, (φ φ j k n k lmt: and we know that the response J q k s proportonal to perturbaton strength appearng n the equaton of moton (n lnear regme. The effect of other two terms, λφ p /m and λφ F, on the other hand, are not so obvous, because there s no drect proportonalty wth k nsde. To study them, let s seperate them out. Consder the followng equaton of moton r = (p /m (1 + λφ p = F (1 + λφ (3 The tme evoluton of an arbtary dynamcal varable A under ths equaton of moton would be where the operators Λ and Λ 1 are Λ Λ 1 DA Dt = (Λ + λλ 1 A (31 = ( p + F m r p = ( p φ + F m r p Note that (3 s not beng descrbed by an Hamltonan, thus although Λ and Λ 1 looks (32 8
9 smlar to the Louvlle operator, they lack the Hermtan property of L. Nevertheless, we can formulate a lnear response theory lke (17 because Hermcty s not requred n the dervaton. Snce flow equaton (8 s no longer vald here, we should use Hesenberg representaton to do the ensemble average for the response to the perturbaton. Take A(t = A (t + A (t, where A (t satsfy DA /Dt = Λ A, the tme evoluton of A n the correspondng unperturbed system. Lnearze (31 nto multply both sde by exp( tλ, there s DA Dt = Λ A + λλ 1 A (t (33 D[exp( tλ A ] Dt = exp( tλ [λλ 1 A (t] whch can be ntergrated to gve So, to lnear order n λ, the measurement Ã(t A (t = λ dτ exp((t τλ [Λ 1 A (τ] Ã(t = Ā(t + λ dτ dγf( exp((t τλ [Λ 1 A (τ] (34 Γ If at t = the system s already n equlbrum, exp((t τλ can be replaced by unty, whch reads Ã(t = Ā(t + λ dτ dγf [Λ 1 A (τ] = Ā(t + λ Γ dτ < Λ 1 A (τ > (35 Plug n the expresson for Λ 1 Ã(t = Ā(t + λ dτ dγf [ Γ = Ā(t λ dτ dγf ( Γ φ ( p m p m φ r A (τ] + F r p A (τ (36 9
10 snce f φ s a constant we should get zero by turnng around both / r and / p, and φ s only a functon of r. Because f p m = k B T f p (36 can be rewrtten as Ã(t = Ā(t + λk B T = Ā(t λk B T dτ dγ Γ dτ dγf Γ f p φ r A (τ φ r A (τ p = Ā(t + λ dτ < Λ 2 A (τ > (37 wth Λ 2 = k B T φ (38 r p Compare (37 wth (35, we see that although perturbaton operators Λ 1 and Λ 2 are algebracally dfferent, the responses they would nduce for an arbtary dynamcal varable are totally dentcal after ensemble averagng. Ths s a truly amazng concluson, snce t proves the exstence of a famly of dynamcal systems wth dfferent equaton of motons, but wth the same macroscopc propertes (n lnear regme. For nstance, Λ 2 can be generated by the followng equaton of moton r p = p /m = F λk B T φ r (39 Although (3 and (39 are dfferent equaton of motons, they are actually equvalent. Now let s go back to thermal conductvty and see how old problems can easly be solved by transferrng the equaton of moton. Snce perturbatons are addtve n lnear regme, (28, 29 are equvalent to r p = p /m ( = F λ (E + k B T h φ 1 r 2 j (φ φ j V j r (4 1
11 As perturbatons terms are all explctly proportonal to k now, we can take the k lmt drectly, and get r = p /m p = F + b ( (E + k B T hi 1 2 wth the k = heat current J q (t gven by (4, and ( b < κ = lm Jq (t > t b 2 ΩT j V j r j r (41 (42 Another way to look at t s to take φ = b r and do equvalence transformaton for (23. All the formalsms above were frst gven by Gllan and Dxon (1983. We know that < J q > = for an equlbrated system. However, Jq (τ and J q (τ are hghly correlated random varables, so nstead of calculatng < J q (t > drectly, we calculate the ensemble average of ther dfference < J q (t J q >, whch must possess much hgher statstcal accuracy as ensemble fluctuatons largely cancel each other. 5 Implementaton on Argon As a model system for the calculaton of thermal conductvty, we choose LJ6-12 argon. There are several purposes n dong ths: To verfy the NEMD algorthm theoretcally derved by Gllan and Dxon. To compare the accuracy between equlbrum correlaton method and NEMD. To see the dfference between lqud state and sold state calculatons. To check the senstvty of results wth respect to cell sze,.e., how many k ponts are enough. In the possble future, put defects nto the system and get emprcal functonal dependence between defect concentraton and thermal conductvty. 11
12 Followng the work by Gllan and Dxon (1983, we mplemented both the equlbrum correlaton and NEMD algorthms on LJ6-12 argon systems. In order to get reasonable results, we must ensure the phase trajectores we generated by MD are hghly smooth, so that the dffernce between the perturbed trajectory and the unperturbed trajectory, whch s very small, won t be nundated by numercal nose. Ths requres us to use a smooth truncated LJ6-12 potental, as any dscontnuty n the potental functon would cause jumps n the phase trajectory: V (r = 4ɛ[(σ/r 12 (σ/r 6 ] r < r γ(r 1 r 2 r < r < r 1 r 1 < r (43 where, for a gven r 1, the constants r and γ are chosen to ensure the contnuty of V and ts frst dervatve at r = r. We ntroduced the cutoff at r 1 = 2σ, same as Gllan and Dxon (1983, and get r = σ, γ = ɛ σ 2 (44 All varables are of double precson n the program to ensure numercal accuracy. Fg.(1 shows a run at T R =.69, DR =.84, NP = 256. The thermal conductvty s n reduced unt, the X axs showng the tme steps, wth one step t = 1ps. A total number of 2745 NEMD runs had been taken wthn the range of 34 runnng steps. For the equlbrum correlaton method, dfferent orgns had been taken (each runnng step after equlbraton as a new orgn, wth correlaton length = 12 steps: the number of effectve orgns s much smaller than the number of as these orgns are hghly correlated, but snce t s a long run we are qute sure that the phase space was fully sampled. The calculated κ for lqud argon at ths condton s.99 ±.3, close to the publshd results by Gllan and Dxon. As we can see from Fg.(1, the comparson between equlbrum correlaton method and the NEMD method by Gllan and Dxon should be consdered as qute satsfactory: not only are the asymptotc value for thermal conductvty close, but also for the whole tme range of the autocorrelaton ntegral. Ths drectly proved the theory developed n prevous sectons of the equvalence transformaton of dynamcal systems, based on whch ths NEMD algorthm 12
13 was put forward. Fg.(2 shows one run at T R =.2983, DR =.9, NP = 256 for sold argon. A total number of 1887 NEMD runs had been taken wthn the range of one mllon steps, the correlaton length beng 45, whch s almost 4 tmes the length needed for the lqud to converge. Snce thermal conductvty n a crystal manly results from the aharmoncty of the lattce,.e., nelastc phonon-phonon collsons, the sze of the smulaton cell s an mportant ssue here, as t drectly controls how many k ponts we have. A run had been taken wth exactly the same condton as Fg.(2, except that NP = 18. The result s shown n Fg.(3: We can see that the equlbrum correlaton curve changed very lttle, but the NEMD curve became shaky after about 15 steps. For overall consderatons we decded that t does not make much dfference between 18 and 256 partcles f we use equlbrum correlaton method. There are several subtle ponts concernng the accuracy of the NEMD method: 1. We must ensure that at t = the perturbaton force be nstantaneously added onto the copy of the orgnal system, not a step later or somethng (for mplementaton n a 5th order Gear method see ball.f. 2. As numercal error (n the fnte-step ntergraton process wll set n for both the perturbed and unperturbed systems, ther trajectores would eventually drft away what they should be. Ths makes the delcate busness of takng the dfferences between two trajectores more and more precarous as tme goes on. We therefore do not recommend NEMD for a very long correlaton run f numercal accuracy can not be guaranteed. The best way to see ths s to test the reversblty of the system after a certan tme, by checkng how much the heat current has changed and compare to the dfference n heat current generated by the perturbaton b. 3. To reduce the error dscussed above, we can ether shorten the ntergraton step, whch s costly, or to ncrease perturbaton strength b to suppress the nose (however, ths wll undermne the valdty of the entre algorthm snce t only works n lnear regme. There would be an optmum value dependng on how many NEMD runs one s gong to take (the more you re gong to take, the smaller the optmum b should be. Usually b s around /σ for argon systems. 13
14 A The Argon Code: ball.f ball.f (65 lnes s a smple molecular dynamcs code for LJ6-12 system wth nteracton potental of Eq.(43 type, able to perform both equlbrum correlaton method and NEMD to get the thermal conductvty of the model system. It uses the temperature couplng method proposed by Berendsen et al (1984, wth couplng tme constant of about 3 ps. Be aware of the fact that temperature couplng must be stopped before NEMD starts, so that the smoothness of the phase trajectory could be guaranteed. ball.f usually generates 5 output fles: temp.out: Ths s the overall record of nstantaneous thermodynamcal varables durng the run, such as the potental energy, pressure and MSD. gr.out: Radal dstrbuton functon. Sometmes you ll fnd B n the fle whch means paragraph breakng, after whch the counter s reset and G(r s re-measured for a new perod of tme. confg: Confguraton fle whch saves the postons and veloctes of partcles, so that later programs can contnue from prevous runs. corr.out: Autocorrelaton functon of heat current n the system, startng from the second row: frst column s the steps, second column s the < J q ( J q (t > functonal, thrd column s the tme ntegral of t. The frst row s the system parameters: temperature, densty and t. nemd.out: Heat current of the perturbed system, startng from the thrd row: ( J x q, J y q, J z q. The frst row s system parameters: temperature, densty and t, the second row s perturbaton strength b : (b x, b y, b z. The nput parameters are put n a fle (suppose called con, whch usually looks lke ths: E 8.. The frst row means 256 partcles wth T R =.2983, DR =.9, the means not readng from the confg fles (f t s 1, then the confg fle of prevous run should be appended behnd 14
15 Thermal conductvty of LJ6 12 system at TR =.69 DR = dashed lne: NEMD algorthm; sold lne: equlbrum correlaton Fgure 1: Flud near trple pont, NP =
16 1.2 Thermal conductvty of LJ6 12 system at TR =.2983 DR = dashed lne: NEMD algorthm; sold lne: equlbrum correlaton Fgure 2: Sold argon, NP =
17 1.2 Thermal conductvty of LJ6 12 system at TR =.2983 DR = dashed lne: NEMD algorthm; sold lne: equlbrum correlaton Fgure 3: Sold argon, NP = 18 17
18 the 4th row. The second row means the total run s 1 steps wth the frst 15 steps as temperature couplng (the couplng constant s 3. ps, and wrtng on temp.out thermodynamcal varables every 5 steps. The thrd row: the frst 1 means equlbrum correlaton turned on, second 1 means NEMD turned on, 45 s correlaton length of both methods (they are the same. The fourth row s perturbaton strength b : (b x, b y, b z. To run, just type: % ball < con & where con s the above condton fle. References [1] H. J. C. Berendsen et al, J. Chem. Phys. 81, 3684 (1984. [2] G. Cccott, G. Jacucc and I. R. McDonald, J. Phys. C 11, (1978. [3] G. Cccott et al, J. Statst. Phys. 21, 1-22 (1979. [4] M. J. Gllan and M. Dxon, J. Phys. C 16, , (1983. [5] M. S. Green, J. Chem. Phys. 22, 398, (1954. [6] L. P. Kadanoff and P. C. Martn, Ann. Phys. (NY 24, (1963. [7] R. Kubo, J. Phys. Soc. Japan 12, 57 (1957. [8] R. Kubo, Rep. Prog. Phys. 29, (1966. [9] H. Mor, Phys. Rev. 112, (
12. The Hamilton-Jacobi Equation Michael Fowler
1. The Hamlton-Jacob Equaton Mchael Fowler Back to Confguraton Space We ve establshed that the acton, regarded as a functon of ts coordnate endponts and tme, satsfes ( ) ( ) S q, t / t+ H qpt,, = 0, and
More informationAdvanced Quantum Mechanics
Advanced Quantum Mechancs Rajdeep Sensarma! sensarma@theory.tfr.res.n ecture #9 QM of Relatvstc Partcles Recap of ast Class Scalar Felds and orentz nvarant actons Complex Scalar Feld and Charge conjugaton
More informationLecture 12: Discrete Laplacian
Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly
More informationLecture 14: Forces and Stresses
The Nuts and Bolts of Frst-Prncples Smulaton Lecture 14: Forces and Stresses Durham, 6th-13th December 2001 CASTEP Developers Group wth support from the ESF ψ k Network Overvew of Lecture Why bother? Theoretcal
More informationCHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng
More information10. Canonical Transformations Michael Fowler
10. Canoncal Transformatons Mchael Fowler Pont Transformatons It s clear that Lagrange s equatons are correct for any reasonable choce of parameters labelng the system confguraton. Let s call our frst
More informationCSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography
CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve
More informationNUMERICAL DIFFERENTIATION
NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the
More informationOpen Systems: Chemical Potential and Partial Molar Quantities Chemical Potential
Open Systems: Chemcal Potental and Partal Molar Quanttes Chemcal Potental For closed systems, we have derved the followng relatonshps: du = TdS pdv dh = TdS + Vdp da = SdT pdv dg = VdP SdT For open systems,
More informationPHYS 705: Classical Mechanics. Canonical Transformation II
1 PHYS 705: Classcal Mechancs Canoncal Transformaton II Example: Harmonc Oscllator f ( x) x m 0 x U( x) x mx x LT U m Defne or L p p mx x x m mx x H px L px p m p x m m H p 1 x m p m 1 m H x p m x m m
More informationDifference Equations
Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1
More informationPhysics 5153 Classical Mechanics. Principle of Virtual Work-1
P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal
More information1 Matrix representations of canonical matrices
1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:
More informationThe Feynman path integral
The Feynman path ntegral Aprl 3, 205 Hesenberg and Schrödnger pctures The Schrödnger wave functon places the tme dependence of a physcal system n the state, ψ, t, where the state s a vector n Hlbert space
More informationCHAPTER 14 GENERAL PERTURBATION THEORY
CHAPTER 4 GENERAL PERTURBATION THEORY 4 Introducton A partcle n orbt around a pont mass or a sphercally symmetrc mass dstrbuton s movng n a gravtatonal potental of the form GM / r In ths potental t moves
More informationELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM
ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look
More information3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X
Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number
More informationCanonical transformations
Canoncal transformatons November 23, 2014 Recall that we have defned a symplectc transformaton to be any lnear transformaton M A B leavng the symplectc form nvarant, Ω AB M A CM B DΩ CD Coordnate transformatons,
More informationProf. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model
EXACT OE-DIMESIOAL ISIG MODEL The one-dmensonal Isng model conssts of a chan of spns, each spn nteractng only wth ts two nearest neghbors. The smple Isng problem n one dmenson can be solved drectly n several
More informationPhysics 181. Particle Systems
Physcs 181 Partcle Systems Overvew In these notes we dscuss the varables approprate to the descrpton of systems of partcles, ther defntons, ther relatons, and ther conservatons laws. We consder a system
More informationTransfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system
Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng
More informationPhysics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1
P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the
More informationOne-sided finite-difference approximations suitable for use with Richardson extrapolation
Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,
More informationThe non-negativity of probabilities and the collapse of state
The non-negatvty of probabltes and the collapse of state Slobodan Prvanovć Insttute of Physcs, P.O. Box 57, 11080 Belgrade, Serba Abstract The dynamcal equaton, beng the combnaton of Schrödnger and Louvlle
More informationSnce h( q^; q) = hq ~ and h( p^ ; p) = hp, one can wrte ~ h hq hp = hq ~hp ~ (7) the uncertanty relaton for an arbtrary state. The states that mnmze t
8.5: Many-body phenomena n condensed matter and atomc physcs Last moded: September, 003 Lecture. Squeezed States In ths lecture we shall contnue the dscusson of coherent states, focusng on ther propertes
More informationSolution Thermodynamics
Soluton hermodynamcs usng Wagner Notaton by Stanley. Howard Department of aterals and etallurgcal Engneerng South Dakota School of nes and echnology Rapd Cty, SD 57701 January 7, 001 Soluton hermodynamcs
More informationLecture Note 3. Eshelby s Inclusion II
ME340B Elastcty of Mcroscopc Structures Stanford Unversty Wnter 004 Lecture Note 3. Eshelby s Incluson II Chrs Wenberger and We Ca c All rghts reserved January 6, 004 Contents 1 Incluson energy n an nfnte
More informationNumerical Heat and Mass Transfer
Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and
More informationEcon107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)
I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes
More informationDensity matrix. c α (t)φ α (q)
Densty matrx Note: ths s supplementary materal. I strongly recommend that you read t for your own nterest. I beleve t wll help wth understandng the quantum ensembles, but t s not necessary to know t n
More informationWeek3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity
Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle
More informationThe equation of motion of a dynamical system is given by a set of differential equations. That is (1)
Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence
More informationFoundations of Arithmetic
Foundatons of Arthmetc Notaton We shall denote the sum and product of numbers n the usual notaton as a 2 + a 2 + a 3 + + a = a, a 1 a 2 a 3 a = a The notaton a b means a dvdes b,.e. ac = b where c s an
More informationLectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix
Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could
More information= z 20 z n. (k 20) + 4 z k = 4
Problem Set #7 solutons 7.2.. (a Fnd the coeffcent of z k n (z + z 5 + z 6 + z 7 + 5, k 20. We use the known seres expanson ( n+l ( z l l z n below: (z + z 5 + z 6 + z 7 + 5 (z 5 ( + z + z 2 + z + 5 5
More informationModule 1 : The equation of continuity. Lecture 1: Equation of Continuity
1 Module 1 : The equaton of contnuty Lecture 1: Equaton of Contnuty 2 Advanced Heat and Mass Transfer: Modules 1. THE EQUATION OF CONTINUITY : Lectures 1-6 () () () (v) (v) Overall Mass Balance Momentum
More informationCausal Diamonds. M. Aghili, L. Bombelli, B. Pilgrim
Causal Damonds M. Aghl, L. Bombell, B. Plgrm Introducton The correcton to volume of a causal nterval due to curvature of spacetme has been done by Myrhem [] and recently by Gbbons & Solodukhn [] and later
More informationThermodynamics and statistical mechanics in materials modelling II
Course MP3 Lecture 8/11/006 (JAE) Course MP3 Lecture 8/11/006 Thermodynamcs and statstcal mechancs n materals modellng II A bref résumé of the physcal concepts used n materals modellng Dr James Ellott.1
More informationNumerical Transient Heat Conduction Experiment
Numercal ransent Heat Conducton Experment OBJECIVE 1. o demonstrate the basc prncples of conducton heat transfer.. o show how the thermal conductvty of a sold can be measured. 3. o demonstrate the use
More informationEPR Paradox and the Physical Meaning of an Experiment in Quantum Mechanics. Vesselin C. Noninski
EPR Paradox and the Physcal Meanng of an Experment n Quantum Mechancs Vesseln C Nonnsk vesselnnonnsk@verzonnet Abstract It s shown that there s one purely determnstc outcome when measurement s made on
More informationRepresentation theory and quantum mechanics tutorial Representation theory and quantum conservation laws
Representaton theory and quantum mechancs tutoral Representaton theory and quantum conservaton laws Justn Campbell August 1, 2017 1 Generaltes on representaton theory 1.1 Let G GL m (R) be a real algebrac
More informationSupplementary Notes for Chapter 9 Mixture Thermodynamics
Supplementary Notes for Chapter 9 Mxture Thermodynamcs Key ponts Nne major topcs of Chapter 9 are revewed below: 1. Notaton and operatonal equatons for mxtures 2. PVTN EOSs for mxtures 3. General effects
More informationChapter 12. Ordinary Differential Equation Boundary Value (BV) Problems
Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(
More information1 GSW Iterative Techniques for y = Ax
1 for y = A I m gong to cheat here. here are a lot of teratve technques that can be used to solve the general case of a set of smultaneous equatons (wrtten n the matr form as y = A), but ths chapter sn
More informationWorkshop: Approximating energies and wave functions Quantum aspects of physical chemistry
Workshop: Approxmatng energes and wave functons Quantum aspects of physcal chemstry http://quantum.bu.edu/pltl/6/6.pdf Last updated Thursday, November 7, 25 7:9:5-5: Copyrght 25 Dan Dll (dan@bu.edu) Department
More informationC/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1
C/CS/Phy9 Problem Set 3 Solutons Out: Oct, 8 Suppose you have two qubts n some arbtrary entangled state ψ You apply the teleportaton protocol to each of the qubts separately What s the resultng state obtaned
More informationThe Order Relation and Trace Inequalities for. Hermitian Operators
Internatonal Mathematcal Forum, Vol 3, 08, no, 507-57 HIKARI Ltd, wwwm-hkarcom https://doorg/0988/mf088055 The Order Relaton and Trace Inequaltes for Hermtan Operators Y Huang School of Informaton Scence
More informationTHE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions
THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructons by George Hardgrove Chemstry Department St. Olaf College Northfeld, MN 55057 hardgrov@lars.acc.stolaf.edu Copyrght George
More informationA particle in a state of uniform motion remain in that state of motion unless acted upon by external force.
The fundamental prncples of classcal mechancs were lad down by Galleo and Newton n the 16th and 17th centures. In 1686, Newton wrote the Prncpa where he gave us three laws of moton, one law of gravty,
More informationModule 9. Lecture 6. Duality in Assignment Problems
Module 9 1 Lecture 6 Dualty n Assgnment Problems In ths lecture we attempt to answer few other mportant questons posed n earler lecture for (AP) and see how some of them can be explaned through the concept
More informationPHYS 705: Classical Mechanics. Newtonian Mechanics
1 PHYS 705: Classcal Mechancs Newtonan Mechancs Quck Revew of Newtonan Mechancs Basc Descrpton: -An dealzed pont partcle or a system of pont partcles n an nertal reference frame [Rgd bodes (ch. 5 later)]
More information2 Finite difference basics
Numersche Methoden 1, WS 11/12 B.J.P. Kaus 2 Fnte dfference bascs Consder the one- The bascs of the fnte dfference method are best understood wth an example. dmensonal transent heat conducton equaton T
More informationPHYS 705: Classical Mechanics. Calculus of Variations II
1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary
More informationImplicit Integration Henyey Method
Implct Integraton Henyey Method In realstc stellar evoluton codes nstead of a drect ntegraton usng for example the Runge-Kutta method one employs an teratve mplct technque. Ths s because the structure
More informationRigid body simulation
Rgd bod smulaton Rgd bod smulaton Once we consder an object wth spacal etent, partcle sstem smulaton s no longer suffcent Problems Problems Unconstraned sstem rotatonal moton torques and angular momentum
More informationMOLECULAR DYNAMICS ,..., What is it? 2 = i i
MOLECULAR DYNAMICS What s t? d d x t 2 m 2 = F ( x 1,..., x N ) =1,,N r ( x1 ( t),..., x ( t)) = v = ( x& 1 ( t ),..., x& ( t )) N N What are some uses of molecular smulatons and modelng? Conformatonal
More informationThis column is a continuation of our previous column
Comparson of Goodness of Ft Statstcs for Lnear Regresson, Part II The authors contnue ther dscusson of the correlaton coeffcent n developng a calbraton for quanttatve analyss. Jerome Workman Jr. and Howard
More informationUniversity of Washington Department of Chemistry Chemistry 453 Winter Quarter 2015
Lecture 2. 1/07/15-1/09/15 Unversty of Washngton Department of Chemstry Chemstry 453 Wnter Quarter 2015 We are not talkng about truth. We are talkng about somethng that seems lke truth. The truth we want
More informationSTATISTICAL MECHANICAL ENSEMBLES 1 MICROSCOPIC AND MACROSCOPIC VARIABLES PHASE SPACE ENSEMBLES. CHE 524 A. Panagiotopoulos 1
CHE 54 A. Panagotopoulos STATSTCAL MECHACAL ESEMBLES MCROSCOPC AD MACROSCOPC ARABLES The central queston n Statstcal Mechancs can be phrased as follows: f partcles (atoms, molecules, electrons, nucle,
More informationSection 8.3 Polar Form of Complex Numbers
80 Chapter 8 Secton 8 Polar Form of Complex Numbers From prevous classes, you may have encountered magnary numbers the square roots of negatve numbers and, more generally, complex numbers whch are the
More informationPHYS 705: Classical Mechanics. Hamilton-Jacobi Equation
1 PHYS 705: Classcal Mechancs Hamlton-Jacob Equaton Hamlton-Jacob Equaton There s also a very elegant relaton between the Hamltonan Formulaton of Mechancs and Quantum Mechancs. To do that, we need to derve
More information(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate
Internatonal Journal of Mathematcs and Systems Scence (018) Volume 1 do:10.494/jmss.v1.815 (Onlne Frst)A Lattce Boltzmann Scheme for Dffuson Equaton n Sphercal Coordnate Debabrata Datta 1 *, T K Pal 1
More informationThis model contains two bonds per unit cell (one along the x-direction and the other along y). So we can rewrite the Hamiltonian as:
1 Problem set #1 1.1. A one-band model on a square lattce Fg. 1 Consder a square lattce wth only nearest-neghbor hoppngs (as shown n the fgure above): H t, j a a j (1.1) where,j stands for nearest neghbors
More informationHomework Assignment 3 Due in class, Thursday October 15
Homework Assgnment 3 Due n class, Thursday October 15 SDS 383C Statstcal Modelng I 1 Rdge regresson and Lasso 1. Get the Prostrate cancer data from http://statweb.stanford.edu/~tbs/elemstatlearn/ datasets/prostate.data.
More informationBézier curves. Michael S. Floater. September 10, These notes provide an introduction to Bézier curves. i=0
Bézer curves Mchael S. Floater September 1, 215 These notes provde an ntroducton to Bézer curves. 1 Bernsten polynomals Recall that a real polynomal of a real varable x R, wth degree n, s a functon of
More informationRobert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations
Quantum Physcs 量 理 Robert Esberg Second edton CH 09 Multelectron atoms ground states and x-ray exctatons 9-01 By gong through the procedure ndcated n the text, develop the tme-ndependent Schroednger equaton
More informationEEE 241: Linear Systems
EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they
More informationA how to guide to second quantization method.
Phys. 67 (Graduate Quantum Mechancs Sprng 2009 Prof. Pu K. Lam. Verson 3 (4/3/2009 A how to gude to second quantzaton method. -> Second quantzaton s a mathematcal notaton desgned to handle dentcal partcle
More informationMMA and GCMMA two methods for nonlinear optimization
MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons
More informationAppendix B. The Finite Difference Scheme
140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton
More informationSupplemental document
Electronc Supplementary Materal (ESI) for Physcal Chemstry Chemcal Physcs. Ths journal s the Owner Socetes 01 Supplemental document Behnam Nkoobakht School of Chemstry, The Unversty of Sydney, Sydney,
More informationprinceton univ. F 17 cos 521: Advanced Algorithm Design Lecture 7: LP Duality Lecturer: Matt Weinberg
prnceton unv. F 17 cos 521: Advanced Algorthm Desgn Lecture 7: LP Dualty Lecturer: Matt Wenberg Scrbe: LP Dualty s an extremely useful tool for analyzng structural propertes of lnear programs. Whle there
More informationBezier curves. Michael S. Floater. August 25, These notes provide an introduction to Bezier curves. i=0
Bezer curves Mchael S. Floater August 25, 211 These notes provde an ntroducton to Bezer curves. 1 Bernsten polynomals Recall that a real polynomal of a real varable x R, wth degree n, s a functon of the
More information1 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations
Physcs 171/271 -Davd Klenfeld - Fall 2005 (revsed Wnter 2011) 1 Dervaton of Rate Equatons from Sngle-Cell Conductance (Hodgkn-Huxley-lke) Equatons We consder a network of many neurons, each of whch obeys
More informationRate of Absorption and Stimulated Emission
MIT Department of Chemstry 5.74, Sprng 005: Introductory Quantum Mechancs II Instructor: Professor Andre Tokmakoff p. 81 Rate of Absorpton and Stmulated Emsson The rate of absorpton nduced by the feld
More informationMore metrics on cartesian products
More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of
More informationPES 1120 Spring 2014, Spendier Lecture 6/Page 1
PES 110 Sprng 014, Spender Lecture 6/Page 1 Lecture today: Chapter 1) Electrc feld due to charge dstrbutons -> charged rod -> charged rng We ntroduced the electrc feld, E. I defned t as an nvsble aura
More informationAnalytical Chemistry Calibration Curve Handout
I. Quck-and Drty Excel Tutoral Analytcal Chemstry Calbraton Curve Handout For those of you wth lttle experence wth Excel, I ve provded some key technques that should help you use the program both for problem
More informationLecture Notes on Linear Regression
Lecture Notes on Lnear Regresson Feng L fl@sdueducn Shandong Unversty, Chna Lnear Regresson Problem In regresson problem, we am at predct a contnuous target value gven an nput feature vector We assume
More informationAffine transformations and convexity
Affne transformatons and convexty The purpose of ths document s to prove some basc propertes of affne transformatons nvolvng convex sets. Here are a few onlne references for background nformaton: http://math.ucr.edu/
More informationChap 5. Dynamics in Condensed Phases
PG/KA/Chap 5-1 Chap 5. Dynamcs n Condensed Phases Generalzed Langevn Equaton Phenomenologcal ntroducton m v = m Γ(t τ)v(τ)dτ + R(t) Γ(t) : frcton kernel frcton depends on the past ( = memory effect : delayed
More informationLinear Feature Engineering 11
Lnear Feature Engneerng 11 2 Least-Squares 2.1 Smple least-squares Consder the followng dataset. We have a bunch of nputs x and correspondng outputs y. The partcular values n ths dataset are x y 0.23 0.19
More informationwhere the sums are over the partcle labels. In general H = p2 2m + V s(r ) V j = V nt (jr, r j j) (5) where V s s the sngle-partcle potental and V nt
Physcs 543 Quantum Mechancs II Fall 998 Hartree-Fock and the Self-consstent Feld Varatonal Methods In the dscusson of statonary perturbaton theory, I mentoned brey the dea of varatonal approxmaton schemes.
More informationInductance Calculation for Conductors of Arbitrary Shape
CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors
More informationMonte Carlo method II
Course MP3 Lecture 5 14/11/2006 Monte Carlo method II How to put some real physcs nto the Monte Carlo method Dr James Ellott 5.1 Monte Carlo method revsted In lecture 4, we ntroduced the Monte Carlo (MC)
More informationStatistical mechanics handout 4
Statstcal mechancs handout 4 Explan dfference between phase space and an. Ensembles As dscussed n handout three atoms n any physcal system can adopt any one of a large number of mcorstates. For a quantum
More information9 Characteristic classes
THEODORE VORONOV DIFFERENTIAL GEOMETRY. Sprng 2009 [under constructon] 9 Characterstc classes 9.1 The frst Chern class of a lne bundle Consder a complex vector bundle E B of rank p. We shall construct
More informationInner Product. Euclidean Space. Orthonormal Basis. Orthogonal
Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,
More informationEntropy generation in a chemical reaction
Entropy generaton n a chemcal reacton E Mranda Área de Cencas Exactas COICET CCT Mendoza 5500 Mendoza, rgentna and Departamento de Físca Unversdad aconal de San Lus 5700 San Lus, rgentna bstract: Entropy
More informationLinear Regression Analysis: Terminology and Notation
ECON 35* -- Secton : Basc Concepts of Regresson Analyss (Page ) Lnear Regresson Analyss: Termnology and Notaton Consder the generc verson of the smple (two-varable) lnear regresson model. It s represented
More informationUncertainty and auto-correlation in. Measurement
Uncertanty and auto-correlaton n arxv:1707.03276v2 [physcs.data-an] 30 Dec 2017 Measurement Markus Schebl Federal Offce of Metrology and Surveyng (BEV), 1160 Venna, Austra E-mal: markus.schebl@bev.gv.at
More informationIntroduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:
CE304, Sprng 2004 Lecture 4 Introducton to Vapor/Lqud Equlbrum, part 2 Raoult s Law: The smplest model that allows us do VLE calculatons s obtaned when we assume that the vapor phase s an deal gas, and
More informationAdvanced Circuits Topics - Part 1 by Dr. Colton (Fall 2017)
Advanced rcuts Topcs - Part by Dr. olton (Fall 07) Part : Some thngs you should already know from Physcs 0 and 45 These are all thngs that you should have learned n Physcs 0 and/or 45. Ths secton s organzed
More informationLinear Approximation with Regularization and Moving Least Squares
Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...
More informationDETAILED DERIVATION OF THE GENERAL MASTER EQUATION IN QUANTUM OPTICS
DETAILED DERIVATION OF THE GENERAL MASTER EQUATION IN QUANTUM OPTICS TIMOTHY JONES Abstract. An unorgnal but more detaled dervaton that typcal s presented of the master equaton. 1. Assumptons We begn by
More informationON MECHANICS WITH VARIABLE NONCOMMUTATIVITY
ON MECHANICS WITH VARIABLE NONCOMMUTATIVITY CIPRIAN ACATRINEI Natonal Insttute of Nuclear Physcs and Engneerng P.O. Box MG-6, 07725-Bucharest, Romana E-mal: acatrne@theory.npne.ro. Receved March 6, 2008
More informationLecture 10 Support Vector Machines II
Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed
More informationGeorgia Tech PHYS 6124 Mathematical Methods of Physics I
Georga Tech PHYS 624 Mathematcal Methods of Physcs I Instructor: Predrag Cvtanovć Fall semester 202 Homework Set #7 due October 30 202 == show all your work for maxmum credt == put labels ttle legends
More informationIntegrals and Invariants of Euler-Lagrange Equations
Lecture 16 Integrals and Invarants of Euler-Lagrange Equatons ME 256 at the Indan Insttute of Scence, Bengaluru Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng,
More informationTHE CURRENT BALANCE Physics 258/259
DSH 1988, 005 THE CURRENT BALANCE Physcs 58/59 The tme average force between two parallel conductors carryng an alternatng current s measured by balancng ths force aganst the gravtatonal force on a set
More information