Formulation of statistical mechanics for chaotic systems
|
|
- Edgar Greene
- 5 years ago
- Views:
Transcription
1 PRAMANA c Inian Acaemy of Sciences Vol. 72, No. 2 journal of February 29 physics pp Formulation of statistical mechanics for chaotic systems VISHNU M BANNUR 1, an RAMESH BABU THAYYULLATHIL 2 1 Department of Physics, University of Calicut, Kozhikoe , Inia 2 Department of Physics, Cochin University of Science an Technology, Kochi , Inia Corresponing author. vmbannur@yahoo.co.in MS receive 14 July 28; revise 21 August 28; accepte 26 August 28 Abstract. We formulate the statistical mechanics of chaotic system with few egrees of freeom an investigate the quartic oscillator system using microcanonical an canonical ensembles. Results of statistical mechanics are numerically verifie by consiering the ynamical evolution of quartic oscillator system with two egrees of freeom. Keywors. Statistical mechanics; chaotic system; quartic oscillators. PACS Nos a; 5.2.Gg; 5.7.-a 1. Introuction The stuy of statistical mechanics an thermoynamics of chaotic systems with few egrees of freeom is very important in unerstaning its various formal aspects from a ynamical point of view [1] an for the stuy of chaotic system using the well-evelope concepts of statistical mechanics [2,3]. Since the trajectory of a chaotic system is almost ergoic in phase-space, it can be approximately escribe by the principles of statistical mechanics. The chaotic systems like quartic oscillator (QO), Henon Heiles oscillator (HHO) etc., with two egrees of freeom exhibit chaos epening on the value of parameters of the system. Hence it is possible to stuy the thermoynamics of the system from the ynamical point of view. Issues relate to the efinition of various thermoynamical quantities such as entropy, temperature etc., are iscusse in the literature [2 4]. There are also other highly ebate issues like, the stuy of Fourier heat law, thermalization of oscillator chain etc. [5 12] using chaotic systems. We know that in a thermoynamical system, as the number of microsystems N the system is riven to ergoicity an the statistical property comes from the implicit assumption of collisions. On the contrary, in a chaotic system, N is finite an the nonlinear interactions present in the Hamiltonian rive the system to an almost ergoic region, provie the system is in a chaotic region. There may be few non-chaotic islans epening on 315
2 Vishnu M Bannur an Ramesh Babu Thayyullathil the values of the parameters, where the system may not show any ergoicity an the statistical properties an thermoynamics of the system may not have much meaning in that region. An initial stuy of statistical mechanics of chaotic system was carrie out by Berichevsky an Alberti [2] using Henon Heiles oscillator an later it was extene to quartic oscillator [3] also. Various issues like concepts of temperature, entropy an istribution function an equipartition of energy were stuie in etail an also verifie numerically using the formulation of microcanonical ensemble. Here, we again strengthen the iea of statistical mechanics of chaotic systems using Kinchin s formulation base on microcanonical ensemble [13]. Further, we exten the stuy to canonical ensemble of such a system an as an example, we consier QO an obtain various thermoynamic quantities an the results are verifie numerically. 2. Henon Heiles oscillator an quartic oscillator First we briefly iscuss QO an HHO moels, stuie extensively in the context of chaos an statistical mechanics of chaos [2,3]. Both are Hamiltonian systems with two egrees of freeom an the Hamiltonians are given by an H = (p2 1 + p 2 2) 2 H = (p2 1 + p 2 2) 2 + (1 α) (q1 4 + q ) q2 1q2, 2 (1) + q q q2 1q q3 2, (2) for QO an HHO respectively. Here q s an p s are generalize coorinates an momenta respectively an α is a parameter. QO is chaotic except for the values of α =, 2 an, an is highly chaotic for α close to 1 as shown in refs [14 16]. HHO is chaotic for energy E = 1/6 an evelops non-chaotic islans as the energy is ecrease. In orer to stuy the statistical mechanics using microcanonical ensemble first we nee to evaluate the boune phase-space volume for a given energy E an it is given by Σ(E) = p 1 p 2 q 1 q 2 = CE 3/2, (3) H E for QO with α 1 an for HHO ( Σ(E) = πe E ) 32 E2 +, (4) as iscusse in etail in ref. [3]. Note that for α = 1 the phase-space volume of QO is infinite an hence unefine. In our numerical work we chose α =.99. In ref. [3] we also iscusse the generalize QO with N egrees of freeom for the Hamiltonian 316 Pramana J. Phys., Vol. 72, No. 2, February 29
3 Statistical mechanics for chaotic systems H = 1 2 N p 2 i i=1 N α ij qi 2 qj 2, (5) i,j=1 where α ij are parameters. This system reuces to our earlier N = 2 quartic oscillator if α 11 = α 22 = (1 α)/6 an α 12 = α 21 = 1/2. Phase-space volume for this system is given by Σ(E) = p 1...p N q 1...q N = C 1 E 3N/4, (6) H E where C 1 is a constant. 3. Microcanonical ensemble In the microcanonical ensemble, one consier a system with constant energy an in our case we consier a Hamiltonian system like HHO or QO with constant energy. Depening on the value of the parameters, the system may exhibit chaos an their phase-space trajectory is almost ergoic in the whole or some region of the phasespace. In those chaotic regions we may have averages of any observable O given as O t = O ph = O en, (7) which follows from ergoic theorem. Here O t, O ph an O en refers to time, phase-space an ensemble averages respectively. For example, the average of momentum square of ith egrees of freeom p 2 i approaches the same value on both time average an phase-space average, as verifie in ref. [3]. The phase-space average is efine as O ph 1 Σ O, (8) Γ E H E+ where Σ, Γ an Ω ( Σ/) are the number of microstates in phase-space hypersphere, hypershell an on the hypersphere surface respectively. Using the above efinition of average we now have the equipartition theorem, [ ] 1 ln Σ p 2 1 = p 2 2 = = T B. (9) Here T B may be ientifie as a temperature by consiering a thought experiment of equilibration of two systems in thermal contact. Suppose we have two systems 1 an 2 an they are brought into thermal contact. After equilibration, the average momentum square of ith particle in system 1 is, [ ] 1 ln p 2 Σ1 i (1) = = T B1, 1 Pramana J. Phys., Vol. 72, No. 2, February
4 Vishnu M Bannur an Ramesh Babu Thayyullathil by equipartition theorem (EPT) on system 1. It may be also evaluate by consiering the combine system (1 + 2) as [ ] 1 ln Σ p 2 i (1+2) = = T B, by EPT on system (1+2) together. Similarly, for the other system [ ] 1 ln p 2 Σ2 j (2) = = T B2, 2 by EPT on system 2 which is also equal to [ ] 1 ln Σ p 2 j (1+2) = = T B, by EPT on system (1+2) together. T B1 = T B2 an hence we may call Therefore, once the system is equilibrate T B = ( ln Σ/) 1 (1) as the temperature of the system. Note that it is ifferent from normal efinition of temperature, [ ] 1 ln( Σ/) T S = = Φ 1, where Φ = ( H/ H 2 ), which was foun recently [4] an reformulate an verifie using quartic oscillators as well as Henon Heiles oscillators in ref. [3]. In the limit N both are the same. Let us now iscuss the entropy of our system which is an extensive variable. Following Khinchin [13] for a system with energy E, which consists of two subsystems 1 an 2 we have, Ω(E) = E 1 Ω 1 (E 1 )Ω 2 (E E 1 ), which can be written as Ω(E) = E 1 Ω 1 (E 1 )Ω 2 (E E 1 ). (11) Since the relation given in eq. (11) is a convolution integral, it is convenient to work with Laplace transform of various Ω(E). Defining φ(α) = Ee αe Ω(E) we write φ(α) = φ 1 (α)φ 2 (α) where φ 1 an φ 2 are Laplace transforms of Ω 1 an Ω 2 respectively. Again, following Khinchin, we choose α = θ as the simple root of the equation ln φ(α) α = E. (12) α=θ 318 Pramana J. Phys., Vol. 72, No. 2, February 29
5 Statistical mechanics for chaotic systems In simple systems like N-quartic oscillators [3] where Σ(E) = CE 3N/4 or ieal gas [17] with Σ(E) = ( ) N V (2πm) 3N/2 h 3 E 3N/2, (3N/2)! etc., this θ is relate to T B via EPT. To see this connection consier systems with Σ(E) = CE l an the corresponing Ω(E) = ClE l 1. Then EPT gives us Since T 1 B φ(α) = = ln Σ = Ω Σ = l E. (13) E e αe ClE l 1 = Cl! α l (14) an using eq. (14) in eq. (12) gives us the relation θ = l E. (15) Comparing eqs (13) an (15) we get the important relation θ = T 1 B. Now the extensive property reas as φ(α) = φ 1 (α)φ 2 (α) φ(t B ) = φ 1 (T B )φ 2 (T B ). In this formalism [13] the entropy is efine as an it becomes S = E θ + ln φ(θ) E T B + ln φ(t B ), (16) S = E 1 + E 2 T B + ln(φ 1 (T B )φ 2 (T B )) = S 1 + S 2, (17) where S 1 an S 2 are entropies of the subsystems. Hence the extensive property of the entropy efine by eq. (16) is very well-establishe. For our examples, an this gives us Finally, it leas to Σ = CE l ln Σ = ln C + l ln E S = ln Σ + ln l! (l ln l l) + constant. S = ln Σ = T 1 B, which is the well-known thermoynamic relation. Pramana J. Phys., Vol. 72, No. 2, February
6 Vishnu M Bannur an Ramesh Babu Thayyullathil 4. Canonical ensemble formulation of quartic oscillators In the classical theory we have the canonical partition function Q = p 1 p 2 q 1 q 2 e β[ p 2 1 +p ((1 α)/6 (q4 1 +q4 2 )+ q2 1 q2 2 )], (18) which reuces to Q = (2 π T ) 3/2 2 1 α 2(4 α) K ( ) 2 + α, (19) 4 α where T is the temperature an β = 1/T an K(z) is the complete elliptical integral of the first kin. From the partition function Q one can obtain the ensity of states which may be useful for the stuy of transition from chaos to integrability of a system [18]. For a more general quartic oscillator with N egrees of freeom we have Q N = N p N qe βh, where H = i p 2 i 2 + i,j α ij 2 q2 i q 2 j, which on simplification reuces to Q N = T 3N/4 C N, (2) where C N is inepenent of T. Following the stanar proceure, various thermoynamic quantities may be erive from the Helmholtz free energy, A = T ln Q N = 3N 4 T ln T T ln C N. (21) The entropy S an the average energy U are an S = 3N 4 ln T + 3N 4 + ln C N U = 3 4 NT. Further, the specific heat C V an square of the energy fluctuation E 2 are obtaine as C V = 3 4 N, E 2 U 2 = 1 3N/4. 32 Pramana J. Phys., Vol. 72, No. 2, February 29
7 Statistical mechanics for chaotic systems It is interesting to note that the energy temperature relation here is the same as that of the microcanonical ensemble formulation, i.e. T 1 B = (3N/4)/E, provie we take the temperature T to be the same as T B. Also note that the temperature T S Φ 1 efine in ref. [4] is obtaine from the efinition Φ = 1 N p N qe β H H Q N H 2, an after some algebra it gives us Φ = β = 1 T, which is also the same as T. This means that T B an T S are the same an equal to T, since T is the temperature of the reservoir with the number of egrees of freeom N an at this limit T B = T S. Thus, for the ynamical evaluation of temperature in canonical ensemble formulation T B may be more preferable, contrary to the observation in ref. [4], since the expression for Φ is more complicate than p 2 i for very large N. 5. Numerical results Numerical results for microcanonical ensemble of QO were reporte earlier [3] where we have showe that the temperatures, which are obtaine by time averaging, agrees with that of the phase-space averaging whenever the system an the initial points are chosen to be in the chaotic region. Here, in the case of canonical ensemble also, the time averages almost agree with that of the phase-space averages for the chaotic QO system. To establish this we have performe a ynamical calculation of various quantities like temperatures (T B an T S ), average energy (U) an specific heat (C V ). In oing this we couple one of the egrees of freeom q 1 of our QO to a thermal bath. An the thermal bath is moelle by Nose Hoover thermostats given by η = q2 1 T 1, (22) where T is the temperature. The final equations of motion which follows from the Hamiltonian in eq. (1) are an q 1 = q 2 = (1 α) q1 3 q 1 q2 2 η q 1 (23) 3 (1 α) q2 3 q 2 q (24) An we nee the solutions to eqs (22) (24). To check for the accuracy of the numerical metho we may use the following constant of motion of the above set of equations: Pramana J. Phys., Vol. 72, No. 2, February
8 Vishnu M Bannur an Ramesh Babu Thayyullathil Cv < T B2 T B1 U T S <.75 < e+6 t Figure 1. Plots of TB of q 1 (TB1) an q 2 (TB2) egrees of freeom, TS, U an CV of chaotic QO system as a function of time. E (p2 1 + p 2 2) 2 + (1 α) q2 1q T η2 2 + T (q q 4 2) t τη(τ). (25) Results are plotte in figure 1 for the thermostat temperature T =.5. Here TB, T S, U an CV are obtaine by time averaging of ynamical quantities. These are plotte as a function of time in figure 1 an asymptotically they all approach the values obtaine from the phase-space averages, i.e., from eq. (21) for N = 2, which are given by T B =.5, T S =.5, U =.75 an C V = 1.5. The QO system is chosen to be highly chaotic with α close to 1. For reasons of numerical stability we took α =.99. As can be seen from the plot, TB1 immeiately approaches.5, the temperature of the thermostats, since it is in irect contact with the reservoir. But the other ynamical quantities TB2, T S, U an CV slowly approach the values.5,.5,.75 an 1.5 respectively in accorance with the ergoic theorem. 6. Conclusions The formulation of statistical mechanics of chaotic systems base on Khinchin s formalism [13] of statistical mechanics for a finite N system is presente. We investigate both microcanonical an canonical ensemble systems of quartic oscillators. The microcanonical ensemble of chaotic systems were stuie an numerically verifie earlier by Berichevsky an Alberti [2], an also by one of us [3]. Here we strengthene the earlier observations an extene the stuy to canonical ensemble 322 Pramana J. Phys., Vol. 72, No. 2, February 29
9 Statistical mechanics for chaotic systems an numerical verification is one using a chaotic QO system. Various quantities like temperatures, energy an specific heat are numerically evaluate by time averaging an they are almost in agreement with our formalism. References [1] M Bianucci, R Manella, B J West an P Grigolini, Phys. Rev. E51, 32 (1995) [2] V L Berichevsky an M V Alberti, Phys. Rev. A44, 858 (1991) [3] V M Bannur, P K Kaw an J C Parikh, Phys. Rev. E55, 2525 (1997) V M Bannur, Phys. Rev. E58, 47 (1998) [4] H H Rugh, Phys. Rev. Lett. 78, 772 (1997) [5] B Hu, L Yang an Y Zhang, Phys. Rev. Lett. 97, (26) [6] L Yang, P Grassberger an B Hu, Phys. Rev. E74, 6211 (26) [7] G R Lee-Daswell, B G Nickel an C G Gray, Phys. Rev. E72, 3122 (25) [8] S Corezzi, M Bianucci an P Grigolini, Phys. Rev. E52, 6881 (1995) [9] G Casati, J For, F Vivali an W M Visscher, Phys. Rev. Lett. 52, 1861 (1984) [1] S Lepri, R Livi an A Politi, Phys. Rev. Lett. 78, 1896 (1997) [11] B Müller an T Trayanov, Phys. Rev. Lett. 68, 3387 (1992) C Gong, Phys. Lett. B298, 257 (1993) T S Biro, C Gong, B Müller an T Trayanov, Int. J. Mo. Phys. C5, 113 (1994) [12] V M Bannur, Phys. Rev. C72, 2494 (25) [13] A I Khinchin, Mathematical founations of statistical mechanics (Dover Pub., 1949) [14] M Lakshmanan an R Sahaevan, Phys. Rep. 224, 1 (1993) [15] M P Joy an M Sabir, Pramana J. Phys. 4, 17 (1993) [16] M S Santhanam an J N Banyopahay, Phys. Rev. Lett. 95, (25) [17] R K Pathria, Statistical mechanics (Butterworth-Heinemann, Oxfor, 1997) [18] H Aiba an T Suzuki, Phys. Rev. E72, (25) Pramana J. Phys., Vol. 72, No. 2, February
Lecture 2 Lagrangian formulation of classical mechanics Mechanics
Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,
More informationChapter 2 Lagrangian Modeling
Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie
More informationensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y
Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay
More informationLie symmetry and Mei conservation law of continuum system
Chin. Phys. B Vol. 20, No. 2 20 020 Lie symmetry an Mei conservation law of continuum system Shi Shen-Yang an Fu Jing-Li Department of Physics, Zhejiang Sci-Tech University, Hangzhou 3008, China Receive
More informationarxiv:cond-mat/ v1 31 Jan 2003
Nosé-Hoover sampling of quantum entangle istribution functions D. Mentrup an J. Schnack 1 arxiv:con-mat/030165 v1 31 Jan 003 Abstract Universität Osnabrück, Fachbereich Physik Barbarastr. 7, D-49069 Osnabrück
More informationComputing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions
Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5
More informationAPPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France
APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation
More informationMath Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors
Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+
More informationELEC3114 Control Systems 1
ELEC34 Control Systems Linear Systems - Moelling - Some Issues Session 2, 2007 Introuction Linear systems may be represente in a number of ifferent ways. Figure shows the relationship between various representations.
More informationStrength Analysis of CFRP Composite Material Considering Multiple Fracture Modes
5--XXXX Strength Analysis of CFRP Composite Material Consiering Multiple Fracture Moes Author, co-author (Do NOT enter this information. It will be pulle from participant tab in MyTechZone) Affiliation
More informationDust Acoustic Compressive Waves in a Warm Dusty Plasma Having Non-Thermal Ions and Non-Isothermal Electrons
Plasma Science an Technology, Vol.17, No.9, Sep. 015 Dust Acoustic Compressive Waves in a Warm Dusty Plasma Having Non-Thermal Ions an Non-Isothermal Electrons Apul N. DEV 1, Manoj K. DEKA, Rajesh SUBEDI
More informationPhysics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1
Physics 5153 Classical Mechanics The Virial Theorem an The Poisson Bracket 1 Introuction In this lecture we will consier two applications of the Hamiltonian. The first, the Virial Theorem, applies to systems
More informationThe derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)
Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)
More informationThe Ehrenfest Theorems
The Ehrenfest Theorems Robert Gilmore Classical Preliminaries A classical system with n egrees of freeom is escribe by n secon orer orinary ifferential equations on the configuration space (n inepenent
More informationAverage value of position for the anharmonic oscillator: Classical versus quantum results
verage value of position for the anharmonic oscillator: Classical versus quantum results R. W. Robinett Department of Physics, The Pennsylvania State University, University Park, Pennsylvania 682 Receive
More informationA Review of Multiple Try MCMC algorithms for Signal Processing
A Review of Multiple Try MCMC algorithms for Signal Processing Luca Martino Image Processing Lab., Universitat e València (Spain) Universia Carlos III e Mari, Leganes (Spain) Abstract Many applications
More informationThermal conductivity of graded composites: Numerical simulations and an effective medium approximation
JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University
More informationA Modification of the Jarque-Bera Test. for Normality
Int. J. Contemp. Math. Sciences, Vol. 8, 01, no. 17, 84-85 HIKARI Lt, www.m-hikari.com http://x.oi.org/10.1988/ijcms.01.9106 A Moification of the Jarque-Bera Test for Normality Moawa El-Fallah Ab El-Salam
More informationChapter 4. Electrostatics of Macroscopic Media
Chapter 4. Electrostatics of Macroscopic Meia 4.1 Multipole Expansion Approximate potentials at large istances 3 x' x' (x') x x' x x Fig 4.1 We consier the potential in the far-fiel region (see Fig. 4.1
More informationSemiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom
PHYSICAL REVIEW A 69, 063409 (2004) Semiclassical analysis of long-wavelength multiphoton processes: The Ryberg atom Luz V. Vela-Arevalo* an Ronal F. Fox Center for Nonlinear Sciences an School of Physics,
More informationIntroduction to the Vlasov-Poisson system
Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its
More informationarxiv:nlin/ v1 [nlin.cd] 21 Mar 2002
Entropy prouction of iffusion in spatially perioic eterministic systems arxiv:nlin/0203046v [nlin.cd] 2 Mar 2002 J. R. Dorfman, P. Gaspar, 2 an T. Gilbert 3 Department of Physics an Institute for Physical
More informationVariational principle for limit cycles of the Rayleigh van der Pol equation
PHYICAL REVIEW E VOLUME 59, NUMBER 5 MAY 999 Variational principle for limit cycles of the Rayleigh van er Pol equation R. D. Benguria an M. C. Depassier Faculta e Física, Pontificia Universia Católica
More informationHyperbolic Moment Equations Using Quadrature-Based Projection Methods
Hyperbolic Moment Equations Using Quarature-Base Projection Methos J. Koellermeier an M. Torrilhon Department of Mathematics, RWTH Aachen University, Aachen, Germany Abstract. Kinetic equations like the
More informationInterpolated Rigid-Body Motions and Robotics
Interpolate Rigi-Boy Motions an Robotics J.M. Selig Faculty of Business, Computing an Info. Management. Lonon South Bank University, Lonon SE AA, U.K. seligjm@lsbu.ac.uk Yaunquing Wu Dept. Mechanical Engineering.
More informationCHM 532 Notes on Creation and Annihilation Operators
CHM 53 Notes on Creation an Annihilation Operators These notes provie the etails concerning the solution to the quantum harmonic oscillator problem using the algebraic metho iscusse in class. The operators
More informationSwitching Time Optimization in Discretized Hybrid Dynamical Systems
Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set
More informationarxiv:physics/ v2 [physics.ed-ph] 23 Sep 2003
Mass reistribution in variable mass systems Célia A. e Sousa an Vítor H. Rorigues Departamento e Física a Universiae e Coimbra, P-3004-516 Coimbra, Portugal arxiv:physics/0211075v2 [physics.e-ph] 23 Sep
More informationDynamic Equations and Nonlinear Dynamics of Cascade Two-Photon Laser
Commun. Theor. Phys. (Beiing, China) 45 (6) pp. 4 48 c International Acaemic Publishers Vol. 45, No. 6, June 5, 6 Dynamic Equations an Nonlinear Dynamics of Cascae Two-Photon Laser XIE Xia,,, HUANG Hong-Bin,
More information(# = %(& )(* +,(- Closed system, well-defined energy (or e.g. E± E/2): Microcanonical ensemble
Recall from before: Internal energy (or Entropy): &, *, - (# = %(& )(* +,(- Closed system, well-defined energy (or e.g. E± E/2): Microcanonical ensemble & = /01Ω maximized Ω: fundamental statistical quantity
More informationarxiv: v1 [physics.class-ph] 20 Dec 2017
arxiv:1712.07328v1 [physics.class-ph] 20 Dec 2017 Demystifying the constancy of the Ermakov-Lewis invariant for a time epenent oscillator T. Pamanabhan IUCAA, Post Bag 4, Ganeshkhin, Pune - 411 007, Inia.
More informationKonrad-Zuse-Zentrum für Informationstechnik Berlin
Konra-use-entrum für Informationstechnik Berlin Christof Schutte Smoothe Molecular Dynamics For Thermally Embee Systems SC 95-4 (May 995) Smoothe Molecular Dynamics For Thermally Embee Systems Christof
More information1 Heisenberg Representation
1 Heisenberg Representation What we have been ealing with so far is calle the Schröinger representation. In this representation, operators are constants an all the time epenence is carrie by the states.
More informationRole of parameters in the stochastic dynamics of a stick-slip oscillator
Proceeing Series of the Brazilian Society of Applie an Computational Mathematics, v. 6, n. 1, 218. Trabalho apresentao no XXXVII CNMAC, S.J. os Campos - SP, 217. Proceeing Series of the Brazilian Society
More informationA Simple Model for the Calculation of Plasma Impedance in Atmospheric Radio Frequency Discharges
Plasma Science an Technology, Vol.16, No.1, Oct. 214 A Simple Moel for the Calculation of Plasma Impeance in Atmospheric Raio Frequency Discharges GE Lei ( ) an ZHANG Yuantao ( ) Shanong Provincial Key
More informationQubit channels that achieve capacity with two states
Qubit channels that achieve capacity with two states Dominic W. Berry Department of Physics, The University of Queenslan, Brisbane, Queenslan 4072, Australia Receive 22 December 2004; publishe 22 March
More informationApplication of the homotopy perturbation method to a magneto-elastico-viscous fluid along a semi-infinite plate
Freun Publishing House Lt., International Journal of Nonlinear Sciences & Numerical Simulation, (9), -, 9 Application of the homotopy perturbation metho to a magneto-elastico-viscous flui along a semi-infinite
More informationLagrangian and Hamiltonian Mechanics
Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical
More informationSome Examples. Uniform motion. Poisson processes on the real line
Some Examples Our immeiate goal is to see some examples of Lévy processes, an/or infinitely-ivisible laws on. Uniform motion Choose an fix a nonranom an efine X := for all (1) Then, {X } is a [nonranom]
More informationEVOLUTION OF PARTICLE SIZE DISTRIBUTION IN AIR IN THE RAINFALL PROCESS VIA THE MOMENT METHOD
137 THERMAL SCIENCE, Year 1, Vol. 16, No. 5, pp. 137-1376 EVOLUTION OF PARTICLE SIZE DISTRIBUTION IN AIR IN THE RAINFALL PROCESS VIA THE MOMENT METHOD by Fu-Jun GAN a an Jian-Zhong LIN a,b * a Department
More informationSchrödinger s equation.
Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of
More informationarxiv: v1 [cond-mat.stat-mech] 9 Jan 2012
arxiv:1201.1836v1 [con-mat.stat-mech] 9 Jan 2012 Externally riven one-imensional Ising moel Amir Aghamohammai a 1, Cina Aghamohammai b 2, & Mohamma Khorrami a 3 a Department of Physics, Alzahra University,
More informationSeparation of Variables
Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical
More informationANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS
ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS MICHAEL HOLST, EVELYN LUNASIN, AND GANTUMUR TSOGTGEREL ABSTRACT. We consier a general family of regularize Navier-Stokes an Magnetohyroynamics
More informationImplicit Differentiation
Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,
More informationImplicit Lyapunov control of closed quantum systems
Joint 48th IEEE Conference on Decision an Control an 28th Chinese Control Conference Shanghai, P.R. China, December 16-18, 29 ThAIn1.4 Implicit Lyapunov control of close quantum systems Shouwei Zhao, Hai
More informationComputing Derivatives
Chapter 2 Computing Derivatives 2.1 Elementary erivative rules Motivating Questions In this section, we strive to unerstan the ieas generate by the following important questions: What are alternate notations
More informationarxiv: v1 [nucl-th] 21 Jan 2019
Paramagnetic squeezing of a uniformly expaning quark-gluon plasma in an out of equilibrium N. Saooghi an S. M. A. Tabatabaee Department of Physics, Sharif University of Technology, P.O. Box 11155-9161,
More informationLagrangian and Hamiltonian Dynamics
Lagrangian an Hamiltonian Dynamics Volker Perlick (Lancaster University) Lecture 1 The Passage from Newtonian to Lagrangian Dynamics (Cockcroft Institute, 22 February 2010) Subjects covere Lecture 2: Discussion
More informationOptimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations
Optimize Schwarz Methos with the Yin-Yang Gri for Shallow Water Equations Abessama Qaouri Recherche en prévision numérique, Atmospheric Science an Technology Directorate, Environment Canaa, Dorval, Québec,
More informationTotal Energy Shaping of a Class of Underactuated Port-Hamiltonian Systems using a New Set of Closed-Loop Potential Shape Variables*
51st IEEE Conference on Decision an Control December 1-13 212. Maui Hawaii USA Total Energy Shaping of a Class of Uneractuate Port-Hamiltonian Systems using a New Set of Close-Loop Potential Shape Variables*
More informationθ x = f ( x,t) could be written as
9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)
More informationTime-of-Arrival Estimation in Non-Line-Of-Sight Environments
2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor
More informationA SIMPLE ENGINEERING MODEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PRODUCTS
International Journal on Engineering Performance-Base Fire Coes, Volume 4, Number 3, p.95-3, A SIMPLE ENGINEERING MOEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PROCTS V. Novozhilov School of Mechanical
More informationHarmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method
1 Harmonic Moelling of Thyristor Briges using a Simplifie Time Domain Metho P. W. Lehn, Senior Member IEEE, an G. Ebner Abstract The paper presents time omain methos for harmonic analysis of a 6-pulse
More informationOn Characterizing the Delay-Performance of Wireless Scheduling Algorithms
On Characterizing the Delay-Performance of Wireless Scheuling Algorithms Xiaojun Lin Center for Wireless Systems an Applications School of Electrical an Computer Engineering, Purue University West Lafayette,
More informationAssignment 1. g i (x 1,..., x n ) dx i = 0. i=1
Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition
More informationarxiv: v1 [physics.flu-dyn] 8 May 2014
Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa
More informationDiscrete Hamilton Jacobi Theory and Discrete Optimal Control
49th IEEE Conference on Decision an Control December 15-17, 2010 Hilton Atlanta Hotel, Atlanta, GA, USA Discrete Hamilton Jacobi Theory an Discrete Optimal Control Tomoi Ohsawa, Anthony M. Bloch, an Melvin
More informationLecture 6: Calculus. In Song Kim. September 7, 2011
Lecture 6: Calculus In Song Kim September 7, 20 Introuction to Differential Calculus In our previous lecture we came up with several ways to analyze functions. We saw previously that the slope of a linear
More informationThe total derivative. Chapter Lagrangian and Eulerian approaches
Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function
More informationEnergy behaviour of the Boris method for charged-particle dynamics
Version of 25 April 218 Energy behaviour of the Boris metho for charge-particle ynamics Ernst Hairer 1, Christian Lubich 2 Abstract The Boris algorithm is a wiely use numerical integrator for the motion
More informationI. INTRODUCTION In this work we analize the slow roll approximation in cosmological moels of nonminimally couple theories of gravity (NMC) starting fr
Slow Roll Ination in Non-Minimally Couple Theories: Hyperextene Gravity Approach Diego F. Torres Departamento e Fsica, Universia Nacional e La Plata C.C. 67, 900, La Plata, Buenos Aires, Argentina Abstract
More informationLower Bounds for the Smoothed Number of Pareto optimal Solutions
Lower Bouns for the Smoothe Number of Pareto optimal Solutions Tobias Brunsch an Heiko Röglin Department of Computer Science, University of Bonn, Germany brunsch@cs.uni-bonn.e, heiko@roeglin.org Abstract.
More informationLinear First-Order Equations
5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)
More informationCompletely passive natural convection
Early View publication on wileyonlinelibrary.com (issue an page numbers not yet assigne; citable using Digital Object Ientifier DOI) ZAMM Z. Angew. Math. Mech., 1 6 (2011) / DOI 10.1002/zamm.201000030
More informationOptimization of Geometries by Energy Minimization
Optimization of Geometries by Energy Minimization by Tracy P. Hamilton Department of Chemistry University of Alabama at Birmingham Birmingham, AL 3594-140 hamilton@uab.eu Copyright Tracy P. Hamilton, 1997.
More informationOptimum design of tuned mass damper systems for seismic structures
Earthquake Resistant Engineering Structures VII 175 Optimum esign of tune mass amper systems for seismic structures I. Abulsalam, M. Al-Janabi & M. G. Al-Taweel Department of Civil Engineering, Faculty
More informationApproximate reduction of dynamic systems
Systems & Control Letters 57 2008 538 545 www.elsevier.com/locate/sysconle Approximate reuction of ynamic systems Paulo Tabuaa a,, Aaron D. Ames b, Agung Julius c, George J. Pappas c a Department of Electrical
More information6 General properties of an autonomous system of two first order ODE
6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x
More information4. Important theorems in quantum mechanics
TFY4215 Kjemisk fysikk og kvantemekanikk - Tillegg 4 1 TILLEGG 4 4. Important theorems in quantum mechanics Before attacking three-imensional potentials in the next chapter, we shall in chapter 4 of this
More informationTHE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE
Journal of Soun an Vibration (1996) 191(3), 397 414 THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE E. M. WEINSTEIN Galaxy Scientific Corporation, 2500 English Creek
More informationChapter 6: Energy-Momentum Tensors
49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.
More informationDelocalization of boundary states in disordered topological insulators
Journal of Physics A: Mathematical an Theoretical J. Phys. A: Math. Theor. 48 (05) FT0 (pp) oi:0.088/75-83/48//ft0 Fast Track Communication Delocalization of bounary states in isorere topological insulators
More informationGeneralized-Type Synchronization of Hyperchaotic Oscillators Using a Vector Signal
Commun. Theor. Phys. (Beijing, China) 44 (25) pp. 72 78 c International Acaemic Publishers Vol. 44, No. 1, July 15, 25 Generalize-Type Synchronization of Hyperchaotic Oscillators Using a Vector Signal
More informationAVERAGE MONTHLY PATTERNS OF SURFACE CURRENTS MEASURED IN THE SOUTHEASTERN MEDITERRANEAN OFF THE EGYPTIAN COAST
77 AVERAGE MONTHLY PATTERNS OF SURFACE CURRENTS MEASURED N THE SOUTHEASTERN MEDTERRANEAN OFF THE EGYPTAN COAST Makram A. Gerges nstitute of Oceanography an Fisheries, Alexanria, Egypt. ABSTRACT The paper
More informationLecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations
Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:
More informationDepartment of Physics, Visva-Bharati University, Santiniketan , India 2
On the non-star Lagrangian equations Aparna Saha Benoy Talukar Department of Physics, Visva-Bharati University, Santiniketan 7335, Inia Nutan Palli, Bolpur 7304, Inia Email. aparna_phyvb@yahoo.co.in binoy3@bsnl.in
More informationA nonlinear inverse problem of the Korteweg-de Vries equation
Bull. Math. Sci. https://oi.org/0.007/s3373-08-025- A nonlinear inverse problem of the Korteweg-e Vries equation Shengqi Lu Miaochao Chen 2 Qilin Liu 3 Receive: 0 March 207 / Revise: 30 April 208 / Accepte:
More informationSome vector algebra and the generalized chain rule Ross Bannister Data Assimilation Research Centre, University of Reading, UK Last updated 10/06/10
Some vector algebra an the generalize chain rule Ross Bannister Data Assimilation Research Centre University of Reaing UK Last upate 10/06/10 1. Introuction an notation As we shall see in these notes the
More informationAgmon Kolmogorov Inequalities on l 2 (Z d )
Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Publishe by Canaian Center of Science an Eucation Agmon Kolmogorov Inequalities on l (Z ) Arman Sahovic Mathematics Department,
More informationSemicompressible Ocean Thermodynamics and Boussinesq Energy Conservation
Article Semicompressible Ocean Thermoynamics an Boussinesq Energy Conservation William K. Dewar 1, *, Joseph Schoonover 1, Trevor McDougall 2 an Rupert Klein 3 1 Department of Earth, Ocean, an Atmospheric
More information12.11 Laplace s Equation in Cylindrical and
SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering
More informationQuantum mechanical approaches to the virial
Quantum mechanical approaches to the virial S.LeBohec Department of Physics an Astronomy, University of Utah, Salt Lae City, UT 84112, USA Date: June 30 th 2015 In this note, we approach the virial from
More informationTrack Initialization from Incomplete Measurements
Track Initialiation from Incomplete Measurements Christian R. Berger, Martina Daun an Wolfgang Koch Department of Electrical an Computer Engineering, University of Connecticut, Storrs, Connecticut 6269,
More informationarxiv: v4 [math.pr] 27 Jul 2016
The Asymptotic Distribution of the Determinant of a Ranom Correlation Matrix arxiv:309768v4 mathpr] 7 Jul 06 AM Hanea a, & GF Nane b a Centre of xcellence for Biosecurity Risk Analysis, University of Melbourne,
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
More informationSpectral Flow, the Magnus Force, and the. Josephson-Anderson Relation
Spectral Flow, the Magnus Force, an the arxiv:con-mat/9602094v1 16 Feb 1996 Josephson-Anerson Relation P. Ao Department of Theoretical Physics Umeå University, S-901 87, Umeå, SWEDEN October 18, 2018 Abstract
More informationThe Press-Schechter mass function
The Press-Schechter mass function To state the obvious: It is important to relate our theories to what we can observe. We have looke at linear perturbation theory, an we have consiere a simple moel for
More informationA Path Planning Method Using Cubic Spiral with Curvature Constraint
A Path Planning Metho Using Cubic Spiral with Curvature Constraint Tzu-Chen Liang an Jing-Sin Liu Institute of Information Science 0, Acaemia Sinica, Nankang, Taipei 5, Taiwan, R.O.C., Email: hartree@iis.sinica.eu.tw
More informationPoS(RAD COR 2007)030. Three-jet DIS final states from k -dependent parton showers. F. Hautmann University of Oxford
OUTP-8--P Three-jet DIS final states from k -epenent parton showers University of Oxfor E-mail: hautmann@thphys.ox.ac.uk H. Jung DESY E-mail: jung@mail.esy.e Experimental measurements of angular correlations
More informationAn inductance lookup table application for analysis of reluctance stepper motor model
ARCHIVES OF ELECTRICAL ENGINEERING VOL. 60(), pp. 5- (0) DOI 0.478/ v07-0-000-y An inuctance lookup table application for analysis of reluctance stepper motor moel JAKUB BERNAT, JAKUB KOŁOTA, SŁAWOMIR
More informationTable of Common Derivatives By David Abraham
Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec
More informationTRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS
TRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS Francesco Bullo Richar M. Murray Control an Dynamical Systems California Institute of Technology Pasaena, CA 91125 Fax : + 1-818-796-8914 email
More informationGlobal Optimization for Algebraic Geometry Computing Runge Kutta Methods
Global Optimization for Algebraic Geometry Computing Runge Kutta Methos Ivan Martino 1 an Giuseppe Nicosia 2 1 Department of Mathematics, Stockholm University, Sween martino@math.su.se 2 Department of
More informationNested Saturation with Guaranteed Real Poles 1
Neste Saturation with Guarantee Real Poles Eric N Johnson 2 an Suresh K Kannan 3 School of Aerospace Engineering Georgia Institute of Technology, Atlanta, GA 3332 Abstract The global stabilization of asymptotically
More informationMath 1271 Solutions for Fall 2005 Final Exam
Math 7 Solutions for Fall 5 Final Eam ) Since the equation + y = e y cannot be rearrange algebraically in orer to write y as an eplicit function of, we must instea ifferentiate this relation implicitly
More informationSection 7.1: Integration by Parts
Section 7.1: Integration by Parts 1. Introuction to Integration Techniques Unlike ifferentiation where there are a large number of rules which allow you (in principle) to ifferentiate any function, the
More informationDifferentiability, Computing Derivatives, Trig Review. Goals:
Secants vs. Derivatives - Unit #3 : Goals: Differentiability, Computing Derivatives, Trig Review Determine when a function is ifferentiable at a point Relate the erivative graph to the the graph of an
More informationNOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy,
NOTES ON EULER-BOOLE SUMMATION JONATHAN M BORWEIN, NEIL J CALKIN, AND DANTE MANNA Abstract We stuy a connection between Euler-MacLaurin Summation an Boole Summation suggeste in an AMM note from 196, which
More information