arxiv: v1 [cond-mat.stat-mech] 22 Dec 2017
|
|
- Shanon Lawrence
- 5 years ago
- Views:
Transcription
1 Notes on the Hybrid Monte Carlo Method arxiv: v1 [cond-at.stat-ech] 22 Dec 2017 Jerey C. Paler, 1, Air Haji-Akbari, 2, Rakesh S. Singh, 2 Fausto Martelli, 3 Roberto Car, 3 Athanassios Z. Panagiotopoulos, 2 and Pablo G. Debenedetti 2, 1 Departent of Cheical and Bioolecular Engineering, University of Houston, Houston, Texas 77204, USA 2 Departent of Cheical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA 3 Departent of Cheistry, Princeton University, Princeton, New Jersey 08544, USA Dated: Deceber 25, 2017 I. INTRODUCTION The hybrid Monte Carlo HMC ethod is a widely used olecular siulation technique for coputing the equilibriu properties of condensed atter systes [1 3]. The basic HMC algorith for sapling fro the canonical NVT enseble consists of three steps [3]: Step 1. Draw a coplete set of initial oenta p {p i } N i=1 fro the distribution P p. Step 2. Propagate a icrocanonical olecular dynaics MD trajectory using a tie-reversible and voluepreserving integrator e.g., the velocity-verlet algorith to take the syste fro state x to x, where x {x i } N i=1 is a coplete set of particle coordinates for the syste. The length of the trajectory t = n steps is specified by the nuber of integration steps n steps and tie step. Step 3. Accept or reject the new configuration x according to the Metropolis-Hastings criterion, P acc {x, p} {x, p } = in 1, e βux P p, 1 e βux P p where Ux is the potential energy of the syste. The ost coon choice for P p is the Maxwell-Boltzann MB distribution [1 3], PMBp std Ce βkp, 2 where β = k B T 1 and C is a teperature specific noralization constant. The kinetic energy Kp is defined using the standard expression fro classical echanics. For rigid bodies, for exaple, p {p co i, Ω i } N i=1, where p co i i v co i, v co i is the linear center of ass co velocity for olecule i, and ω i and Ω i denote olecule i s angular velocity and angular oentu, respectively. The j th coponent of olecule i s angular oentu j = 1, 2, 3 is Ω j = I j,k ω i,k k = 1, 2, 3, where I j,k denotes the generic olecule s i in this case inertia tensor. The kinetic energy for rigid bodies is thus given by Kp 1 N 2 i=1 [pco i p co i / i + ω i Ω i ]. Insertion of P MB p into Eq. 1 yields: Pacc std {x, p} {x, p } = in 1, e β[hx,p Hx,p], 3 where Hx, p Hx, p Ux Ux+Kp Kp is the difference between the classical Hailtonians Hx, p and Hx, p. The superscript std in Eqs. 2 and 3 denotes that these are the standard choices for P p and P acc, which are used in the overwheling ajority of HMC siulation studies of condensed atter systes. When used together, they ensure that the HMC algorith will satisfy detailed balance and thus asyptotically saple configurations fro the Boltzann distribution P eq x e βux. Sapling can also be perfored using other initial oentu distributions that are even functions of p i.e., P p = P p [2, 3]. As we deonstrate below, however, it is generally necessary to odify the acceptance criterion as prescribed by Eq. 1 to preserve detailed balance. Hence if P p P MB p, then using the standard acceptance criterion in Eq. 3 will generally result in detailed balance violations and concoitant sapling errors. Corresponding author:jcpaler@uh.edu Current address: Departent of Cheical and Environental Engineering, Yale University, New Haven, CT 06520, USA Corresponding author:pdebene@princeton.edu
2 2 II. DETAILED BALANCE The detailed balance condition is given by: P eq xγx x = P eq x Γx x, 4 where P eq x is the equilibriu distribution and Γx x is the Markov transition probability fro state x to x. This condition is sufficient to ensure that P eq x will be a stationary state for the Markov processes generated by an MC sapling algorith. Here, we exaine the constraints Eq. 4 iposes on the choice of P p in the HMC algorith under the following conditions: i The MD trajectories are propagated using a tie-reversible and volue-preserving integration schee. The latter properties ensures that dxdp = dx dp ; non-volue-preserving integrators ay be used in HMC, but a Jacobian factor ust be incorporated into the acceptance criterion Eq. 1 to correct for the resulting copression of phase space [4]. ii P eq x is given by the Boltzann distribution P eq x e βux iii The trial trajectories are accepted using the standard criterion P std acc Eq. 3. iv The Hailtonian Hx, p is an even function of p i.e., Hx, p = Hx, p and defined using the standard expression for kinetic energy fro classical echanics. v The function P p is unspecified, but has the following properties: a b P p = P p P p is stationary i.e., the definition of P p does not change with tie. The initial p and final p oenta are thus both treated as variates fro this distribution. Note that this does not iply P p = P p, which would suggest unifority. For HMC, the forward transition probability fro x to x is given by [4]: Γx x = dpdp P pg{x, p} {x, p }Pacc std {x, p} {x, p } 5 = dpdp P pg{x, p} {x, p }in 1, e β[hx,p Hx,p], 6 where g{x, p} {x, p } = δx x tδp p tδx x0δp p0 and δ is the Kronecker delta function. The function g{x, p} {x, p } gives the probability that an MD trajectory originating fro {x0, p0} = {x, p} will end up in state {x t, p t} = {x, p } after an elapsed tie t. Siilarly, the reverse transition probability Γx x can be written as [4]: Γx x = dpdp P pg{x, p} {x, p }Pacc std {x, p} {x, p } 7 = dpdp P pg{x, p} {x, p }in 1, e β[hx,p Hx,p]. 8 In Eqs. 7 and 8 we follow the notation convention used in Ref. [4], whereby p denotes the generic set of oenta over which integration is perfored, and hence we can use p in conjunction with x. The notation is arbitrary, provided that one defines and uses p and p consistently and does not confuse these variables in anipulating the equations or aking variable substitutions. In particular, because P refers to a generic oentu distribution associated with the initial state of the transition under consideration, one can write, for a x x transition Eqs. 7 or 8 P pg{x, p} {x, p }Pacc std {x, p} {x, p } or P p g{x, p } {x, p}pacc std {x, p } {x, p}. Changing redefining integration variables [5] p p and p p and noting that dpdp = d [ p ] d [ p] yields: Γx x = dpdp P p g{x, p } {x, p}in 1, e β[hx, p Hx, p ]. 9 For a tie-reversible and volue-preserving integrator g{x, p } {x, p} = g{x, p} {x, p }. Additionally, assuing Hx, p = Hx, p and P p = P p, we obtain: Γx x = dpdp P p g{x, p} {x, p }in 1, e β[hx,p Hx,p ]. 10
3 Invoking the identity in 1, e β[hx,p Hx,p ] = e βhp,x e βhp,x in Γx x = = = e βux e βux 1, e β[hx,p Hx,p] : dpdp P p g{x, p} {x, p } e βhp,x e in 1, e β[hx,p Hx,p] βhp,x 11 dpdp P pg{x, p} {x, p } P p e βux e βkp P p e βux e in 1, e β[hx,p Hx,p] βkp 12 dpdp P pg{x, p} {x, p } P p e βkp P pe βkp in 1, e β[hx,p Hx,p]. 13 Coparison with Eq. 6 reveals that the integrands differ by the boxed ter in Eq. 13. To satisfy detailed balance Eq. 4, the choice of P p ust therefore either ensure that the boxed ter in Eq. 13 is always unity or has no effect on the value of the integral when the integration over p and p is carried out. It is clear fro our derivation that neither of these criteria will be autoatically satisfied by an arbitrary distribution P p, even if it has the properties listed in point v above. The forer criterion is satisfied, however, by the Maxwell-Boltzann distribution P p = PMB stdp Eq. 2. Inserting the MB distribution into Eqs. 6 and 13, we find that the boxed ter in Eq. 13 is identically equal to 1. Subsequent inspection of Eqs. 6 and 13 reveals that Γx x = e βux Γx x, which is the detailed balance e condition Eq. 4. By contrast, identifying a distribution that unabiguously βux satisfies only the latter criterion is extreely challenging, as the integrals in Eq. 6 and 13 cannot usually be evaluated analytically. Other choices for P p are certainly possible, but HMC acceptance criterion in Eqs. 6 and 13 ust be redefined as prescribed by Eq. 1 to avoid violating detailed balance. 3 III. DETAILED BALANCE FOR THE ONE-DIMENSIONAL HARMONIC OSCILLATOR The analysis above deonstrates that choosing P p = PMB std p is generally required to satisfy detailed balance when HMC sapling is perfored using the standard acceptance criterion Pacc std Eq. 3. Here, we show analytically for a siple odel syste and velocity Verlet integrator that detailed balance can only be satisfied by P p = PMB stdp when P acc = Pacc std. For this case, the condition that P p = PMB std p is thus not only sufficient but also necessary. Consider the siple case of applying the HMC algorith to saple the configurations of a particle of ass in a one-diensional haronic well described by the potential energy function Ux = 1 2 kx2. Further, suppose that we define the kinetic energy using the standard expression Kp p2 2 and consider HMC oves that consist of only a single olecular dynaics integration step propagated using the tie reversible and volue-preserving velocity Verlet schee. Under these conditions, {x, p } is related to {x, p} via: x = x + p fx p = p [fx + fx], 15 where is the integration tie step and fx = dux dx = kx is the force acting on the particle. The detailed balance condition Eq. 4 ust be satisfied for all x and x. Thus we are free to choose x and x arbitrarily to siplify our analysis. Letting x = 0 and x 0, and noting that f0 = 0 and fx = kx, we obtain: x = p 16 p = p 1 2 kx = p1 1 2 k2. 17 These equations show that for a given choice of the sapling paraeter, only p = x will lead to an attepted HMC ove fro {x, p} {x, p }. Siilarly, they also reveal that p = x 1 2 kx is unique for a given. Therefore g will be given by: Additionally, we find that: g{x, p} {x, p } = δ p x δ p x kx. 18 Hx, p = x Hx, p = x k2 x
4 Using these expressions, we can now analyze the conditions under which the detailed balance Eq. 4 will be satisfied. For the selected x and x, the detailed balance condition Eq. 4 reduces to: Γx x = e βux e βux Γx x 21 = e β 2 kx 2 Γx x, 22 where we have substituted Ux = 0 and Ux = 1 2 kx 2. Inserting Eqs. 18, 19, and 20 into the expression for the forward transition probability given by Eq. 6, we obtain: Γx x = dpdp P pg{x, p} {x, p }in 1, e β[hx,p Hx,p] 23 = dpdp P pδ p x δ p x + 1 [ ] 2 kx in 1, e β x k2 x 2 2 x x = P in 1, e β 8 k2 x = P x 4 e β 8 k2 x 2 2, 26 where the last equality holds because β 8 k2 x The reverse transition probability can also be evaluated in a siilar fashion fro Eq. 11: Γx x = dpdp P p g{x, p} {x, p } e βhp,x e in 1, e β[hx,p Hx,p] 27 βhp,x = dpdp P p δ p x δ p x + 1 e βhp,x 2 kx e in 1, e β[hx,p Hx,p] 28 βhp,x x = P 1 2 kx e β 8 k2 x 2 2 in 1, e β 8 k2 x x = P 1 2 kx e β 8 k2 x 2 2 e β 8 k2 x = P x 1 2 kx Substituting Eqs. 26 and 31 into Eq. 22, we find that for detailed balance to hold, we need to have: x P e β 8 k2 x 2 2 = e β x 2 kx 2 P 1 2 kx which is clearly not satisfied for an arbitrary P. Letting φp := ln P p + βp 2 /2, we obtain: x ln P β x 8 k2 x 2 2 = ln P 1 2 kx β 2 kx 2 x φ βx β x 8 k2 x 2 2 = φ 1 2 kx β x kx β 2 kx 2 x βx φ β x 8 k2 x 2 2 = φ 1 βx 2 kx β 2 kx 2 β 8 k2 x 2 2 β 2 kx 2 x x φ = φ 1 2 kx which can be atheatically stated as φbx = φax with b = /, a = / k/2. Because φbx = φax for all x, or ore generally, φb 2 x = φ [bbx ] = φ [abx ] = φ [bax ] = φ [aax ] = φa 2 x, φb n x = φa n x 35
5 for n Z. Siilarly, φx = φ [b n x b n ] = φ [a n x b n ] [ a ] n = φ x. 36 b Without loss of generality, suppose a < b. Because Eq. 36 is valid for any n Z, it ust also hold in the liit that n. Taking this liit, we find that: [ a ] n [ φx = φ x a ] n = li b φ x = φ0 = const. n b Thus, in order for Eq. 33 to hold for any x, φp needs to be constant: or equivalently, φp = ln P p + βp 2 /2 = const., 37 p2 β P p e 2 = e βkp, 38 which is the Maxwell-Boltzann distribution. Hence, this derivation proves for the 1-D haronic oscillator odel syste that for detailed balance to hold when the standard HMC acceptance criterion Eq. 3 is used in conjunction with Kp p2 2, P p can only have the Maxwell-Boltzann distribution. 5 IV. A STATIONARY AND EVEN MOMENTUM DISTRIBUTION IS NOT SUFFICIENT FOR DETAILED BALANCE The derivations outlined in Secs. II and III assue that P p is even i.e., P p = P p and stationary i.e., the initial oenta p and final oenta p are treated as variates fro the sae distribution, P p. As deonstrated by our analysis, however, these properties are not sufficient to satisfy detailed balance. This fact can also be shown nuerically by perforing HMC siulations of the Lennard-Jones odel for argon σ Ar = n, ɛ Ar /k B = K, in which the potential is truncated at 3 σ Ar and standard long-range tail corrections are applied. We consider three HMC schees in which the target sapling teperature is specified by β 1 = k B T 1 1 : Schee I: HMC using P p e β2kp and Pacc std = in 1, e β 1 Ux e β 1 Kp, where β e β 1 Ux e β 1 Kp 2 = β 1. This schee corresponds to the standard HMC algorith, where the initial oenta are drawn fro the Maxwell-Boltzann distribution at the target sapling teperature β 1 = β 2 = k B T 1 1. As proved in Sec. II, this schee satisfies detailed balance. Schee II: HMC using P p e β2kp and Pacc od = in 1, e β 1 Ux e β 2 Kp, where β e β 1 Ux e β 2 Kp 2 < β 1. In this schee, the initial oenta are drawn fro a Maxwell-Boltzann distribution corresponding to an artificially high teperature β 2 in 1, e β 1 Ux P p e β 1 Ux P p detailed balance is satisfied for all choices of β 1 and β 2. = k B T 2 1, but the acceptance criterion has been odified such that Pacc od =, as prescribed by Eq. 1. This odified acceptance criterion is sufficient to ensure that = in 1, e β 1 Ux e β 1 Kp Schee III: HMC using P p e β2kp and Pacc std, where β 2 < β 1. In this schee, the initial oenta are drawn fro a Maxwell-Boltzann distribution corresponding to an artificially e β 1 Ux e β 1 Kp high teperature β 2 = k B T 2 1, but no odifications are ade to the standard acceptance criterion. As proved in Sec. III, this sapling schee, does not satisfy detailed balance, even though P p is an even function and stationary Fig. 1 We siulate N = 500 particles at a fixed density of g/c 3 and choose T 1 = K as our target sapling teperature. In ters of Lennard-Jones units, this corresponds to ρ = and T = 0.9. In each case, sapling is perfored for HMC oves, where each ove consists of 10 olecular dynaics integration steps using a 30 fs tie step. As expected, Schee I correctly predicts the average potential energy, as illustrated by the excellent agreeent with benchark data fro the National Institute for Science and Technology for the Lennard-Jones fluid Table I. As expected, it also satisfies the noralization condition e β H = 1 Table I. This rigorous statistical echanical relationship ust hold for all valid HMC sapling schees. Thus it provides a convenient consistency check for
6 detecting sapling errors associated with detailed balance violations. In the absence of sapling errors, we expect that e β H 1 σ 1 1, where σ is the estiated uncertainty. Accordingly, we find that Schee I satisfies this expectation. Siilarly, Schee II correctly predicts the average potential energy and obeys the noralization condition. Even though the initial velocities in Schee II are drawn fro the Maxwell-Boltzann distribution at an artificially high teperature T 2 = K, the odified HMC acceptance ensures that detailed balance is satisfied and that the correct result is recovered. By contrast, Schee III does not satisfy detailed balance and hence fails to obey the noralization condition and correctly predict the equation of state Table I. The sapling errors produced by Schee III becoe ore pronounced as the velocity generation teperature T 2 increases fro to K. TABLE I: Equation of state data for Lennard-Jones argon at a target sapling teperature T 1 = K U /N kj ol 1 e β H 1 σ 1 NIST Ref. Data a b Schee I T 1 = K Schee II T 1 = K; T 2 = K Schee III T 1 = K; T 2 = K Schee III T 1 = K; T 2 = K Schee III T 1 = K; T 2 = K a Benchark data fro the National Institute of Standards and Technology which has been converted fro Lennard-Jones units b Nuber in parentheses denotes uncertainty in the last significant digit 6 FIG. 1: The proposed red and accepted blue velocity distributions are statistically indistinguishable for all sapling schees and choices of T 2. Thus, the velocity distribution is stationary regardless of whether detailed balance is satisfied by the algorith. The green lines show the Maxwell-Boltzann distribution at the indicated teperature T 2. Schee III does not satisfy detailed balance and hence does not yield the correct equation of state for the Lennard-Jones fluid Table I, even though the velocity distribution is an even function and stationary. These two criteria are therefore not sufficient to ensure detailed balance In each sapling schee, P p is an even function of p by construction. Additionally, P p is stationary under each sapling schee, regardless of whether detailed balance is satisfied by the algorith Fig. 1. Thus, in agreeent with our proof in Sec. III, these results deonstrate that stationarity and evenness of P p are not sufficient conditions for detailed balance. If one wishes to use the standard HMC acceptance criterion and definition of kinetic energy given in Sec. I, initial velocities should be drawn fro the Maxwell-Boltzann distribution at the target sapling
7 teperature to satisfy detailed balance see Sec. III. As Schee II illustrates, detailed balance can be satisfied using other choices for P p, if the HMC acceptance criterion is odified appropriately. 7 V. HMC USING PARTIAL MOMENTUM UPDATES The HMC algorith outlined in Sec. I assues that a copletely new set of initial oenta are drawn fro P p in Step 1. It also possible to perfor HMC using correlated saples fro P p that are generated by partially refreshing the oenta fro the end of the previous MC step. Wagoner and Pande [6] have rigorously proved, however, that the sign of the oenta ust be changed negated either upon acceptance or rejection of each HMC ove for this schee to satisfy balance i.e., leave P eq x e βux stationary and thus saple fro the correct equilibriu distribution. Indeed, nuerical studies have shown that oitting this iportant step results in detectable sapling errors [7]. [1] S. Duane, A. D. Kennedy, B. J. Pendleton, and D. Roweth, Phys. Lett. B 195, [2] B. Mehlig, D. W. Heerann, and B. M. Forest, Phys. Rev. B 45, [3] M. Allen and D. Quigley, Mol. Phys. 111, [4] N. Matubayasi and M. Nakahara, J. Che. Phys. 110, [5] M. Tuckeran, Statistical Mechanics: Theory and Molecular Siulation Oxford University Press, [6] J. A. Wagoner and V. S. Pande, J. Che. Phys. 137, [7] E. Akhatskaya, N. Bou-Rabee, and S. Reich, J. Cop. Phys. 19,
Ph 20.3 Numerical Solution of Ordinary Differential Equations
Ph 20.3 Nuerical Solution of Ordinary Differential Equations Due: Week 5 -v20170314- This Assignent So far, your assignents have tried to failiarize you with the hardware and software in the Physics Coputing
More information12 Towards hydrodynamic equations J Nonlinear Dynamics II: Continuum Systems Lecture 12 Spring 2015
18.354J Nonlinear Dynaics II: Continuu Systes Lecture 12 Spring 2015 12 Towards hydrodynaic equations The previous classes focussed on the continuu description of static (tie-independent) elastic systes.
More information2 Q 10. Likewise, in case of multiple particles, the corresponding density in 2 must be averaged over all
Lecture 6 Introduction to kinetic theory of plasa waves Introduction to kinetic theory So far we have been odeling plasa dynaics using fluid equations. The assuption has been that the pressure can be either
More informationKinetic Theory of Gases: Elementary Ideas
Kinetic Theory of Gases: Eleentary Ideas 17th February 2010 1 Kinetic Theory: A Discussion Based on a Siplified iew of the Motion of Gases 1.1 Pressure: Consul Engel and Reid Ch. 33.1) for a discussion
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 11 Jan 2007
Transport and Helfand oents in the Lennard-Jones fluid. II. Theral conductivity arxiv:cond-at/7125v1 [cond-at.stat-ech] 11 Jan 27 S. Viscardy, J. Servantie, and P. Gaspard Center for Nonlinear Phenoena
More informationEXPLORING PHASE SPACES OF BIOMOLECULES WITH MONTE CARLO METHODS
Journal of Optoelectronics and Advanced Materials Vol. 7, o. 3, June 2005, p. 1563-1571 EXPLORIG PHASE SPACES OF BIOMOLECULES WITH MOTE CARLO METHODS A. Bu u * ational Institute for Research and Developent
More informationKinetic Theory of Gases: Elementary Ideas
Kinetic Theory of Gases: Eleentary Ideas 9th February 011 1 Kinetic Theory: A Discussion Based on a Siplified iew of the Motion of Gases 1.1 Pressure: Consul Engel and Reid Ch. 33.1) for a discussion of
More informationFirst of all, because the base kets evolve according to the "wrong sign" Schrödinger equation (see pp ),
HW7.nb HW #7. Free particle path integral a) Propagator To siplify the notation, we write t t t, x x x and work in D. Since x i, p j i i j, we can just construct the 3D solution. First of all, because
More informationClassical systems in equilibrium
35 Classical systes in equilibriu Ideal gas Distinguishable particles Here we assue that every particle can be labeled by an index i... and distinguished fro any other particle by its label if not by any
More informationPhysics 139B Solutions to Homework Set 3 Fall 2009
Physics 139B Solutions to Hoework Set 3 Fall 009 1. Consider a particle of ass attached to a rigid assless rod of fixed length R whose other end is fixed at the origin. The rod is free to rotate about
More informationLecture #8-3 Oscillations, Simple Harmonic Motion
Lecture #8-3 Oscillations Siple Haronic Motion So far we have considered two basic types of otion: translation and rotation. But these are not the only two types of otion we can observe in every day life.
More informationNonuniqueness of canonical ensemble theory. arising from microcanonical basis
onuniueness of canonical enseble theory arising fro icrocanonical basis arxiv:uant-ph/99097 v2 25 Oct 2000 Suiyoshi Abe and A. K. Rajagopal 2 College of Science and Technology, ihon University, Funabashi,
More informationAnnealing contour Monte Carlo algorithm for structure optimization in an off-lattice protein model
JOURNAL OF CHEMICAL PHYSICS VOLUME 120, NUMBER 14 8 APRIL 2004 Annealing contour Monte Carlo algorith for structure optiization in an off-lattice protein odel Faing Liang a) Departent of Statistics, Texas
More informationAn Approximate Model for the Theoretical Prediction of the Velocity Increase in the Intermediate Ballistics Period
An Approxiate Model for the Theoretical Prediction of the Velocity... 77 Central European Journal of Energetic Materials, 205, 2(), 77-88 ISSN 2353-843 An Approxiate Model for the Theoretical Prediction
More informationExtension of CSRSM for the Parametric Study of the Face Stability of Pressurized Tunnels
Extension of CSRSM for the Paraetric Study of the Face Stability of Pressurized Tunnels Guilhe Mollon 1, Daniel Dias 2, and Abdul-Haid Soubra 3, M.ASCE 1 LGCIE, INSA Lyon, Université de Lyon, Doaine scientifique
More informationSchool of Chemistry, University of Nottingham, University Park, Nottingham, NG7 2RD, United Kingdom
pubs.acs.org/jpcb Calculation of Partition Functions and Free Energies of a Binary Mixture Using the Energy Partitioning Method: Application to Carbon Dioxide and Methane Haina Do,* Jonathan D. Hirst,
More informationOn the approximation of Feynman-Kac path integrals
On the approxiation of Feynan-Kac path integrals Stephen D. Bond, Brian B. Laird, and Benedict J. Leikuhler University of California, San Diego, Departents of Matheatics and Cheistry, La Jolla, CA 993,
More informationChapter 1 Introduction and Kinetics of Particles
Chapter 1 Introduction and Kinetics of Particles 1.1 Introduction There are two ain approaches in siulating the transport equations (heat, ass, and oentu), continuu and discrete. In continuu approach,
More informationKinetic Molecular Theory of Ideal Gases
Lecture -3. Kinetic Molecular Theory of Ideal Gases Last Lecture. IGL is a purely epirical law - solely the consequence of experiental obserations Explains the behaior of gases oer a liited range of conditions.
More informationIllustration of transition path theory on a collection of simple examples
THE JOURNAL OF CHEMICAL PHYSICS 125, 084110 2006 Illustration of transition path theory on a collection of siple exaples Philipp Metzner a and Christof Schütte b Departent of Matheatics and Coputer Science,
More informationPhysically Based Modeling CS Notes Spring 1997 Particle Collision and Contact
Physically Based Modeling CS 15-863 Notes Spring 1997 Particle Collision and Contact 1 Collisions with Springs Suppose we wanted to ipleent a particle siulator with a floor : a solid horizontal plane which
More informationMolecular dynamics algorithm for multiple time scales: Systems with long range forces
Molecular dynaics algorith for ultiple tie scales: Systes with long range forces Mark E. Tuckeran@ Bruce J. Berne Departent of Cheistry, Colubia University, New York, New York 10027 Glenn J. Martyna Departent
More informationQuasistationary distributions of dissipative nonlinear quantum oscillators in strong periodic driving fields
PHYSICAL REVIEW E VOLUME 61, NUMBER 5 MAY 2 Quasistationary distributions of dissipative nonlinear quantu oscillators in strong periodic driving fields Heinz-Peter Breuer, 1 Wolfgang Huber, 2 and Francesco
More informationForce and dynamics with a spring, analytic approach
Force and dynaics with a spring, analytic approach It ay strie you as strange that the first force we will discuss will be that of a spring. It is not one of the four Universal forces and we don t use
More informationComparison of Charged Particle Tracking Methods for Non-Uniform Magnetic Fields. Hann-Shin Mao and Richard E. Wirz
42nd AIAA Plasadynaics and Lasers Conferencein conjunction with the8th Internati 27-30 June 20, Honolulu, Hawaii AIAA 20-3739 Coparison of Charged Particle Tracking Methods for Non-Unifor Magnetic
More information13 Harmonic oscillator revisited: Dirac s approach and introduction to Second Quantization
3 Haronic oscillator revisited: Dirac s approach and introduction to Second Quantization. Dirac cae up with a ore elegant way to solve the haronic oscillator proble. We will now study this approach. The
More informationIII.H Zeroth Order Hydrodynamics
III.H Zeroth Order Hydrodynaics As a first approxiation, we shall assue that in local equilibriu, the density f 1 at each point in space can be represented as in eq.iii.56, i.e. f 0 1 p, q, t = n q, t
More informationKinetic Molecular Theory of. IGL is a purely empirical law - solely the
Lecture -3. Kinetic Molecular Theory of Ideal Gases Last Lecture. IGL is a purely epirical law - solely the consequence of experiental obserations Explains the behaior of gases oer a liited range of conditions.
More informationCh 12: Variations on Backpropagation
Ch 2: Variations on Backpropagation The basic backpropagation algorith is too slow for ost practical applications. It ay take days or weeks of coputer tie. We deonstrate why the backpropagation algorith
More informationThe Thermal Conductivity Theory of Non-uniform Granular Flow and the Mechanism Analysis
Coun. Theor. Phys. Beijing, China) 40 00) pp. 49 498 c International Acadeic Publishers Vol. 40, No. 4, October 5, 00 The Theral Conductivity Theory of Non-unifor Granular Flow and the Mechanis Analysis
More informationScattering and bound states
Chapter Scattering and bound states In this chapter we give a review of quantu-echanical scattering theory. We focus on the relation between the scattering aplitude of a potential and its bound states
More informationSpine Fin Efficiency A Three Sided Pyramidal Fin of Equilateral Triangular Cross-Sectional Area
Proceedings of the 006 WSEAS/IASME International Conference on Heat and Mass Transfer, Miai, Florida, USA, January 18-0, 006 (pp13-18) Spine Fin Efficiency A Three Sided Pyraidal Fin of Equilateral Triangular
More informationInspection; structural health monitoring; reliability; Bayesian analysis; updating; decision analysis; value of information
Cite as: Straub D. (2014). Value of inforation analysis with structural reliability ethods. Structural Safety, 49: 75-86. Value of Inforation Analysis with Structural Reliability Methods Daniel Straub
More informationUfuk Demirci* and Feza Kerestecioglu**
1 INDIRECT ADAPTIVE CONTROL OF MISSILES Ufuk Deirci* and Feza Kerestecioglu** *Turkish Navy Guided Missile Test Station, Beykoz, Istanbul, TURKEY **Departent of Electrical and Electronics Engineering,
More informationSEISMIC FRAGILITY ANALYSIS
9 th ASCE Specialty Conference on Probabilistic Mechanics and Structural Reliability PMC24 SEISMIC FRAGILITY ANALYSIS C. Kafali, Student M. ASCE Cornell University, Ithaca, NY 483 ck22@cornell.edu M. Grigoriu,
More informationPHY 171. Lecture 14. (February 16, 2012)
PHY 171 Lecture 14 (February 16, 212) In the last lecture, we looked at a quantitative connection between acroscopic and icroscopic quantities by deriving an expression for pressure based on the assuptions
More informationChapter 12. Quantum gases Microcanonical ensemble
Chapter 2 Quantu gases In classical statistical echanics, we evaluated therodynaic relations often for an ideal gas, which approxiates a real gas in the highly diluted liit. An iportant difference between
More informationNonlinear Stabilization of a Spherical Particle Trapped in an Optical Tweezer
Nonlinear Stabilization of a Spherical Particle Trapped in an Optical Tweezer Aruna Ranaweera ranawera@engineering.ucsb.edu Bassa Baieh baieh@engineering.ucsb.edu Andrew R. Teel teel@ece.ucsb.edu Departent
More informationIN modern society that various systems have become more
Developent of Reliability Function in -Coponent Standby Redundant Syste with Priority Based on Maxiu Entropy Principle Ryosuke Hirata, Ikuo Arizono, Ryosuke Toohiro, Satoshi Oigawa, and Yasuhiko Takeoto
More informationMolecular interactions in beams
Molecular interactions in beas notable advanceent in the experiental study of interolecular forces has coe fro the developent of olecular beas, which consist of a narrow bea of particles, all having the
More informationNon-Parametric Non-Line-of-Sight Identification 1
Non-Paraetric Non-Line-of-Sight Identification Sinan Gezici, Hisashi Kobayashi and H. Vincent Poor Departent of Electrical Engineering School of Engineering and Applied Science Princeton University, Princeton,
More informationChemistry 432 Problem Set 11 Spring 2018 Solutions
1. Show that for an ideal gas Cheistry 432 Proble Set 11 Spring 2018 Solutions P V 2 3 < KE > where is the average kinetic energy of the gas olecules. P 1 3 ρ v2 KE 1 2 v2 ρ N V P V 1 3 N v2 2 3 N
More informationUse of PSO in Parameter Estimation of Robot Dynamics; Part One: No Need for Parameterization
Use of PSO in Paraeter Estiation of Robot Dynaics; Part One: No Need for Paraeterization Hossein Jahandideh, Mehrzad Navar Abstract Offline procedures for estiating paraeters of robot dynaics are practically
More informationOn the characterization of non-linear diffusion equations. An application in soil mechanics
On the characterization of non-linear diffusion equations. An application in soil echanics GARCÍA-ROS, G., ALHAMA, I., CÁNOVAS, M *. Civil Engineering Departent Universidad Politécnica de Cartagena Paseo
More informationBootstrapping Dependent Data
Bootstrapping Dependent Data One of the key issues confronting bootstrap resapling approxiations is how to deal with dependent data. Consider a sequence fx t g n t= of dependent rando variables. Clearly
More informationA simple phenomenologic model for particle transport in spaceperiodic potentials in underdamped systems
A siple phenoenologic odel for particle transport in spaceperiodic potentials in underdaped systes IG MARCHENKO 1,(a,b), II MARCHENKO 3, A ZHIGLO 1 1 NSC Kharov Institute of Physics and Technology, Aadeichesaya
More informationSupporting information for Self-assembly of multicomponent structures in and out of equilibrium
Supporting inforation for Self-assebly of ulticoponent structures in and out of equilibriu Stephen Whitela 1, Rebecca Schulan 2, Lester Hedges 1 1 Molecular Foundry, Lawrence Berkeley National Laboratory,
More informationSF Chemical Kinetics.
SF Cheical Kinetics. Lecture 5. Microscopic theory of cheical reaction inetics. Microscopic theories of cheical reaction inetics. basic ai is to calculate the rate constant for a cheical reaction fro first
More informationChapter 6 1-D Continuous Groups
Chapter 6 1-D Continuous Groups Continuous groups consist of group eleents labelled by one or ore continuous variables, say a 1, a 2,, a r, where each variable has a well- defined range. This chapter explores:
More informationQ5 We know that a mass at the end of a spring when displaced will perform simple m harmonic oscillations with a period given by T = 2!
Chapter 4.1 Q1 n oscillation is any otion in which the displaceent of a particle fro a fixed point keeps changing direction and there is a periodicity in the otion i.e. the otion repeats in soe way. In
More informationChapter 1: Basics of Vibrations for Simple Mechanical Systems
Chapter 1: Basics of Vibrations for Siple Mechanical Systes Introduction: The fundaentals of Sound and Vibrations are part of the broader field of echanics, with strong connections to classical echanics,
More informationSome consequences of a Universal Tension arising from Dark Energy for structures from Atomic Nuclei to Galaxy Clusters
unning Head: Universal Tension fro DE Article Type: Original esearch Soe consequences of a Universal Tension arising fro Dark Energy for structures fro Atoic Nuclei to Galaxy Clusters C Sivara Indian Institute
More informationEstimation of the Mean of the Exponential Distribution Using Maximum Ranked Set Sampling with Unequal Samples
Open Journal of Statistics, 4, 4, 64-649 Published Online Septeber 4 in SciRes http//wwwscirporg/ournal/os http//ddoiorg/436/os4486 Estiation of the Mean of the Eponential Distribution Using Maiu Ranked
More informationKINETIC THEORY. Contents
KINETIC THEORY This brief paper on inetic theory deals with three topics: the hypotheses on which the theory is founded, the calculation of pressure and absolute teperature of an ideal gas and the principal
More informationModel Fitting. CURM Background Material, Fall 2014 Dr. Doreen De Leon
Model Fitting CURM Background Material, Fall 014 Dr. Doreen De Leon 1 Introduction Given a set of data points, we often want to fit a selected odel or type to the data (e.g., we suspect an exponential
More informationBROWNIAN DYNAMICS Lecture notes
Göran Wahnströ BROWNIAN DYNAMICS Lecture notes Göteborg, 6 Deceber 6 Brownian dynaics Brownian otion is the observed erratic otion of particles suspended in a fluid (a liquid or a gas) resulting fro their
More informationOptical Properties of Plasmas of High-Z Elements
Forschungszentru Karlsruhe Techni und Uwelt Wissenschaftlishe Berichte FZK Optical Properties of Plasas of High-Z Eleents V.Tolach 1, G.Miloshevsy 1, H.Würz Project Kernfusion 1 Heat and Mass Transfer
More informationANALYSIS OF HALL-EFFECT THRUSTERS AND ION ENGINES FOR EARTH-TO-MOON TRANSFER
IEPC 003-0034 ANALYSIS OF HALL-EFFECT THRUSTERS AND ION ENGINES FOR EARTH-TO-MOON TRANSFER A. Bober, M. Guelan Asher Space Research Institute, Technion-Israel Institute of Technology, 3000 Haifa, Israel
More informationAn analytical relation between relaxation time spectrum and molecular weight distribution
An analytical relation between relaxation tie spectru and olecular weight distribution Wolfgang Thi, Christian Friedrich, a) Michael Marth, and Josef Honerkap b) Freiburger Materialforschungszentru, Stefan-Meier-Straße
More informationP235 Midterm Examination Prof. Cline
P235 Mier Exaination Prof. Cline THIS IS A CLOSED BOOK EXAMINATION. Do all parts of all four questions. Show all steps to get full credit. 7:00-10.00p, 30 October 2009 1:(20pts) Consider a rocket fired
More informationStatistical associating fluid theory-dimer of the solid phase of the pearl-necklace model
Indian Journal of Pure & Applied Physics Vol., epteber 00, pp 688-69 tatistical associating fluid theory-dier of the solid phase of the pearl-necklace odel Vivek Kuar ingh & K N Khanna Departent of Physics,
More informationCelal S. Konor Release 1.1 (identical to 1.0) 3/21/08. 1-Hybrid isentropic-sigma vertical coordinate and governing equations in the free atmosphere
Celal S. Konor Release. (identical to.0) 3/2/08 -Hybrid isentropic-siga vertical coordinate governing equations in the free atosphere This section describes the equations in the free atosphere of the odel.
More informationA NEW APPROACH FOR CALCULATING AVERAGE CROSS SECTIONS IN THE UNRESOLVED ENERGY REGION
Nuclear Matheatical Coputational Sciences: A Century in Review, A Century Anew Gatlinburg, Tennessee, April 6-, 2003, on CD-ROM, Aerican Nuclear Society, LaGrange Park, IL (2003) A NEW APPROACH FOR CALCULATING
More information8.1 Force Laws Hooke s Law
8.1 Force Laws There are forces that don't change appreciably fro one instant to another, which we refer to as constant in tie, and forces that don't change appreciably fro one point to another, which
More informationIntelligent Systems: Reasoning and Recognition. Perceptrons and Support Vector Machines
Intelligent Systes: Reasoning and Recognition Jaes L. Crowley osig 1 Winter Seester 2018 Lesson 6 27 February 2018 Outline Perceptrons and Support Vector achines Notation...2 Linear odels...3 Lines, Planes
More informationOn the Diffusion Coefficient: The Einstein Relation and Beyond 3
Stoch. Models, Vol. 19, No. 3, 2003, (383-405) Research Report No. 424, 2001, Dept. Theoret. Statist. Aarhus On the Diffusion Coefficient: The Einstein Relation and Beyond 3 GORAN PESKIR 33 We present
More informationQuantum Chemistry Exam 2 Take-home Solutions
Cheistry 60 Fall 07 Dr Jean M Standard Nae KEY Quantu Cheistry Exa Take-hoe Solutions 5) (0 points) In this proble, the nonlinear variation ethod will be used to deterine an approxiate solution for the
More informationSOLUTIONS. PROBLEM 1. The Hamiltonian of the particle in the gravitational field can be written as, x 0, + U(x), U(x) =
SOLUTIONS PROBLEM 1. The Hailtonian of the particle in the gravitational field can be written as { Ĥ = ˆp2, x 0, + U(x), U(x) = (1) 2 gx, x > 0. The siplest estiate coes fro the uncertainty relation. If
More informationDarcy s law describes water flux determined by a water potential gradient: q w = k dψ dx
6 Hydraulic Properties The soil hydraulic properties are of utost iportance for a large variety of soil processes and in particular for soil water flow. The two ost coon hydraulic properties are the soil
More informationSupplemental Material for Correlation between Length and Tilt of Lipid Tails
Suppleental Material for Correlation between Length and Tilt of Lipid Tails Ditry I. Kopelevich and John F. Nagle I. RESULTS FOR ALTERNATIVE DIRECTOR DEFINITIONS A. Alternative Director Definitions The
More informationKinetic Theory of Gases. Chapter 33 1/6/2017. Kinetic Theory of Gases
1/6/017 Kinetic Theory of Gases Kinetic Theory of Gases Chapter 33 Kinetic theory of gases envisions gases as a collection of atos or olecules in otion. Atos or olecules are considered as particles. This
More information1 The properties of gases The perfect gas
1 The properties of gases 1A The perfect gas Answers to discussion questions 1A. The partial pressure of a gas in a ixture of gases is the pressure the gas would exert if it occupied alone the sae container
More informationPhysics 2210 Fall smartphysics 20 Conservation of Angular Momentum 21 Simple Harmonic Motion 11/23/2015
Physics 2210 Fall 2015 sartphysics 20 Conservation of Angular Moentu 21 Siple Haronic Motion 11/23/2015 Exa 4: sartphysics units 14-20 Midter Exa 2: Day: Fri Dec. 04, 2015 Tie: regular class tie Section
More informationChapter 2: Introduction to Damping in Free and Forced Vibrations
Chapter 2: Introduction to Daping in Free and Forced Vibrations This chapter ainly deals with the effect of daping in two conditions like free and forced excitation of echanical systes. Daping plays an
More informationCHECKLIST. r r. Newton s Second Law. natural frequency ω o (rad.s -1 ) (Eq ) a03/p1/waves/waves doc 9:19 AM 29/03/05 1
PHYS12 Physics 1 FUNDAMENTALS Module 3 OSCILLATIONS & WAVES Text Physics by Hecht Chapter 1 OSCILLATIONS Sections: 1.5 1.6 Exaples: 1.6 1.7 1.8 1.9 CHECKLIST Haronic otion, periodic otion, siple haronic
More informationDimensions and Units
Civil Engineering Hydraulics Mechanics of Fluids and Modeling Diensions and Units You already know how iportant using the correct diensions can be in the analysis of a proble in fluid echanics If you don
More informationPhysics 140 D100 Midterm Exam 2 Solutions 2017 Nov 10
There are 10 ultiple choice questions. Select the correct answer for each one and ark it on the bubble for on the cover sheet. Each question has only one correct answer. (2 arks each) 1. An inertial reference
More informationAccuracy of the Scaling Law for Experimental Natural Frequencies of Rectangular Thin Plates
The 9th Conference of Mechanical Engineering Network of Thailand 9- October 005, Phuket, Thailand Accuracy of the caling Law for Experiental Natural Frequencies of Rectangular Thin Plates Anawat Na songkhla
More informationOptimizing energy potentials for success in protein tertiary structure prediction Ting-Lan Chiu 1 and Richard A Goldstein 1,2
Research Paper 223 Optiizing energy potentials for success in protein tertiary structure prediction Ting-Lan Chiu 1 and Richard A Goldstein 1,2 Background: Success in solving the protein structure prediction
More informationBayesian Approach for Fatigue Life Prediction from Field Inspection
Bayesian Approach for Fatigue Life Prediction fro Field Inspection Dawn An and Jooho Choi School of Aerospace & Mechanical Engineering, Korea Aerospace University, Goyang, Seoul, Korea Srira Pattabhiraan
More information(a) Why cannot the Carnot cycle be applied in the real world? Because it would have to run infinitely slowly, which is not useful.
PHSX 446 FINAL EXAM Spring 25 First, soe basic knowledge questions You need not show work here; just give the answer More than one answer ight apply Don t waste tie transcribing answers; just write on
More informationm potential kinetic forms of energy.
Spring, Chapter : A. near the surface of the earth. The forces of gravity and an ideal spring are conservative forces. With only the forces of an ideal spring and gravity acting on a ass, energy F F will
More informationAnalyzing Simulation Results
Analyzing Siulation Results Dr. John Mellor-Cruey Departent of Coputer Science Rice University johnc@cs.rice.edu COMP 528 Lecture 20 31 March 2005 Topics for Today Model verification Model validation Transient
More informationFigure 1: Equivalent electric (RC) circuit of a neurons membrane
Exercise: Leaky integrate and fire odel of neural spike generation This exercise investigates a siplified odel of how neurons spike in response to current inputs, one of the ost fundaental properties of
More informationMolecular dynamic simulation of the melting and solidification processes of argon
Journal of Mechanical Science and Technology 23 (200) 1563~1570 Journal of Mechanical Science and Technology www.springerlink.co/content/1738-44x DOI 10.1007/s12206-00-0418-0 Molecular dynaic siulation
More information1 Bounding the Margin
COS 511: Theoretical Machine Learning Lecturer: Rob Schapire Lecture #12 Scribe: Jian Min Si March 14, 2013 1 Bounding the Margin We are continuing the proof of a bound on the generalization error of AdaBoost
More informationBasic concept of dynamics 3 (Dynamics of a rigid body)
Vehicle Dynaics (Lecture 3-3) Basic concept of dynaics 3 (Dynaics of a rigid body) Oct. 1, 2015 김성수 Vehicle Dynaics Model q How to describe vehicle otion? Need Reference fraes and Coordinate systes 2 Equations
More informationMSEC MODELING OF DEGRADATION PROCESSES TO OBTAIN AN OPTIMAL SOLUTION FOR MAINTENANCE AND PERFORMANCE
Proceeding of the ASME 9 International Manufacturing Science and Engineering Conference MSEC9 October 4-7, 9, West Lafayette, Indiana, USA MSEC9-8466 MODELING OF DEGRADATION PROCESSES TO OBTAIN AN OPTIMAL
More informationSolidification of Porous Material under Natural Convection by Three Phases Modeling
Solidification of Porous Material under Natural Convection by Three Phases Modeling Hassan Basirat Tabrizi, Meber, IAENG and F. Sadeghpour Abstract The perforance of natural convective flow over a rectangular
More informationP (t) = P (t = 0) + F t Conclusion: If we wait long enough, the velocity of an electron will diverge, which is obviously impossible and wrong.
4 Phys520.nb 2 Drude theory ~ Chapter in textbook 2.. The relaxation tie approxiation Here we treat electrons as a free ideal gas (classical) 2... Totally ignore interactions/scatterings Under a static
More informationKinematics and dynamics, a computational approach
Kineatics and dynaics, a coputational approach We begin the discussion of nuerical approaches to echanics with the definition for the velocity r r ( t t) r ( t) v( t) li li or r( t t) r( t) v( t) t for
More informationAbout the definition of parameters and regimes of active two-port networks with variable loads on the basis of projective geometry
About the definition of paraeters and regies of active two-port networks with variable loads on the basis of projective geoetry PENN ALEXANDR nstitute of Electronic Engineering and Nanotechnologies "D
More informationDonald Fussell. October 28, Computer Science Department The University of Texas at Austin. Point Masses and Force Fields.
s Vector Moving s and Coputer Science Departent The University of Texas at Austin October 28, 2014 s Vector Moving s Siple classical dynaics - point asses oved by forces Point asses can odel particles
More informationResearch in Area of Longevity of Sylphon Scraies
IOP Conference Series: Earth and Environental Science PAPER OPEN ACCESS Research in Area of Longevity of Sylphon Scraies To cite this article: Natalia Y Golovina and Svetlana Y Krivosheeva 2018 IOP Conf.
More informationWater a) 48 o b) 53 o c) 41.5 o d) 44 o. Glass. PHYSICS 223 Exam-2 NAME II III IV
PHYSICS 3 Exa- NAME. In the figure shown, light travels fro aterial I, through three layers of other aterials with surfaces parallel to one another, and then back into another layer of aterial I. The refractions
More informationTHERMAL ENDURANCE OF UNREINFORCED UNSATURATED POLYESTERS AND VINYL ESTER RESINS = (1) ln = COMPOSITES & POLYCON 2009
Aerican Coposites Manufacturers Association January 15-17, 29 Tapa, FL USA Abstract THERMAL ENDURANCE OF UNREINFORCED UNSATURATED POLYESTERS AND VINYL ESTER RESINS by Thore M. Klaveness, Reichhold AS In
More informationExperimental Design For Model Discrimination And Precise Parameter Estimation In WDS Analysis
City University of New York (CUNY) CUNY Acadeic Works International Conference on Hydroinforatics 8-1-2014 Experiental Design For Model Discriination And Precise Paraeter Estiation In WDS Analysis Giovanna
More information16.512, Rocket Propulsion Prof. Manuel Martinez-Sanchez Lecture 30: Dynamics of Turbopump Systems: The Shuttle Engine
6.5, Rocket Propulsion Prof. Manuel Martinez-Sanchez Lecture 30: Dynaics of Turbopup Systes: The Shuttle Engine Dynaics of the Space Shuttle Main Engine Oxidizer Pressurization Subsystes Selected Sub-Model
More information2009 Academic Challenge
009 Acadeic Challenge PHYSICS TEST - REGIONAL This Test Consists of 5 Questions Physics Test Production Tea Len Stor, Eastern Illinois University Author/Tea Leader Doug Brandt, Eastern Illinois University
More informationSOLUTIONS for Homework #3
SOLUTIONS for Hoework #3 1. In the potential of given for there is no unboun states. Boun states have positive energies E n labele by an integer n. For each energy level E, two syetrically locate classical
More information