p" N p" 2 p" 1 -ij" 1 =0 (W + ) ... q" 1 q" 2 q'' N =0 (W - )
|
|
- Leslie Pierce
- 6 years ago
- Views:
Transcription
1 Few-Body Systems 0, 6 (5) Few- Body Systems cfl y Springer-Verlag 5 PrintedinAustria Dynamical Hierarchy in Transition States of Reactions C.B. Li ; Λ, A. Shojiguchi ; ΛΛ,M.Toda ; ΛΛΛ,T. Komatsuzaki ;3; y Nonlinear Science Laoratory, Department of Earth and Planetary Sciences, Faculty of Science, Koe University, JST/CREST, Nada, Koe Japan DepartmentofPhysics, Faculty of Science, Nara Women's University, Nara , Japan 3 Department of Theoretical Studies, Institute for Molecular Science, Myodaiji, Okazaki , Japan Astract. We present a partial normalization procedure of Lie canonical perturation theory to elucidate the phase space geometry of the transition state in the multidimensional phase space for a wide range of energy aove the threshold. State selectivity and dynamical correlation along the evolution of reactions will also e discussed. Introduction All reactions take place through the `edge of staility,' which chemists have long envisioned as transition state (TS). The concept of TS has een much more widely utilized than just for the prolems near instaility in chemistry, e.g., atomic physics, cluster physics, celestial mechanics, and iology []. In 930s, the TS is defined y Wigner as a no-return dividing hypersurface in the phase space through which reacting species passes only once on the way from the reactant to the product. In other terms, the no-recrossing condition requires that the `one dimensional' motion normal to the TS (reaction mode) is separale from the motions tangent to the TS (ath modes). While the existence and the definaility of such a no-return TS were well estalished in two degrees of freedom (dof) Hamiltonian systems, those dynamical prolems in more than two dof has een an open prolem until very recently [, 3, 4, 5, 6]. For a N-dof Hamiltonian system, we can expand it near a saddle with index one (an equilirium point with one negative Hessian eigenvalue) as H = E 0 + NX i= p i +! i q i Λ address: cli@koe-u.ac.jp ΛΛ address: shojiguchi@ki-rin.phys.nara-wu.ac.jp ΛΛΛ address: toda@ki-rin.phys.nara-wu.ac.jp y address: tamiki@koe-u.ac.jp + X n=3 H n ; ()
2 Dynamical Hierarchy in Transition States of Reactions where E 0 is the energy of the saddle point and H n consists of terms of n-th degree in (p ; ;p N ;q ; ;q N ). (p ;q ) are the momentum and coordinate of the reactive normal mode, and (p ; q )(= (p ; ;p N ;q ; ;q N )) are those of the ath normal modes, respectively. The unpertured frequency of the reactive mode! is pure imaginary and those of the ath modes! i (i ) are real. In an energy regime slightly aove the saddle point energy (i.e., P quasi-regular N regime []) where! i (i ) do not satisfy the resonance condition n i= i! i ο = 0 (for any integers fn i g unless all n i 's are zero), H can e fully normalized into the Birkhoff normal form H 0 (a classical analog of Dunham Hamiltonian in quantum mechanics) y a canonical transformation (p; q)! (p 0 ; q 0 ) in terms of the Lie canonical perturation theory (LCPT). H 0 (p 0 ; q 0 )=H 0 (J 0 )=E 0 + NX i=! i J 0 i + X i;j a ij J 0 ij 0 j + X i;j;k a ijk J 0 ij 0 jj 0 k + ; () where J 0 = i(p0 =j! j j! jq 0 )= and J 0 i =(p0 i =! i +! q 0 i )= are,respectively, the `action' of the transformed reactive normal form coordinate (p 0 ;q0 ) and those of the ath normal form coordinates (p 0 ; q0 ). a ij;a ijk ;::: are coupling constants among J 0.Notethat the fully normalized Hamiltonian H 0 does not involve any `angle' variale conjugated to J 0. Thus, all J 0 are invariants of motion (i.e., H 0 is integrale) and the crossing dynamics can e solved analytically in full. As the energy of the system increases, nonlinear resonances among the ath modes ecome significant and the overlapping of these resonances gives rise to the transition to gloal chaos. Nevertheless, in the region of index one saddles, any resonance condition can never e satisfied among the pure imaginary frequency (associated with the reactive mode) and the real frequencies (with the ath modes) [0]. This leads to the existence of the following partial normal form H for a wide range of energy aove the threshold (i.e., semi-chaotic regime [, 5]); H (J ; p ; q )=E 0 +! J + NX i= p i +! i q i + f (J ; p ; q )+f (p ; q ); (3) where J = i(p =j! j j! jq )= and the functions f and f contain anharmonic terms in power of the partial normal form variales (p ; q ). In particular, f is defined so that f = 0 when J = 0. Note that Eq. (3) does not involve the `angle' variale conjugated to J, resulting in that J is an invariant of motion. However, in f and f, resonance terms can exist among the ath modes (p ; q) (i.e., H is non-integrale). Since solely the reactive dof is normalized, only (p ;q ) should e called normal form coordinate while the ath modes (p ; q ) remains to e normal mode" coordinates. Note, however, that (p ; q) are functions of the original normal mode momenta and coordinates (p; q) and so they are different from the original set. Most of the theoretical studies [4, 6, 7, 8] have focused only on the quasiregular regime. In this paper, y using a partial LCPT to normalize solely the reactive mode, we discuss how nonlinear resonances among the ath modes, which appear as energy increases aove the threshold, affect the invariance of the normal form action of the reactive mode and the state selectivity of reactions.
3 C.B. Li et al. 3 The partial normal form in the region of index one saddles First, let us riefly look into the equation of motion of the partial normal form H in Eq. (3): dp = q =!! q ; dq = dp dq = @p p ; p! ; =; ;N where! (J ; p; q )=@H =@J =! is the normalized frequency of the reactive dof. Note here that while! 0 in the full normal form H0 is an invariant of motion,! in H varies in time through its dependence on the ath modes (p ). Nevertheless, the equation of motion of (p ) is formally the same as ; q ;q that of (p 0 ;q0) in H0 [8], and preserves the hyperolic structure on the (q ;p) plane (ecause J = constant). For instance, the condition p = q = 0 (yielding dp =dt = dq =dt =0)definesa(N 3)-dimension normally hyperolic invariant manifold (NHIM) on the (N ) dimension equienergy surface in the phase space R N. Normal hyperolicity means that the stretching and contraction rates of the motion normal to the manifold (i.e., the `reactive' dof) dominate those of the motions tangent to the manifold (i.e., the ath dof). Fig. shows the geometrical structure in the multi-dimensional phase space in the region of an index one saddle of the partial normal form Hamiltonian. Note that Eq. (4) tells us that the velocity dq =dt does not change its sign does not change during the crossing. This implies that the = 0 (except ackward reactions). Therefore, the local no-return property ofs in H is guaranteed if! does not change its sign while the system crosses the saddle. The change of sign if the sign of! system never turns ack toahypersurface S defined y q of! requires the passage of! reakdown of normal hyperolicity, that is, ( + (4) =0,which corresponds to an instantaneous" (J ; p; q)) J = =0 0: It was found [9] that, as the energy increases near to the reakdown of normal hyperolicity, the system recrosses the hypersurface S much earlier than finding this condition. However, one of the most important feature of this condition is the dependency on (p ; q ), implying that some portions of the NHIM" start to reak the normal hyperolicity, not gloally at once ut locally, as energy increases. We apply the partial LCPT to the following 3-dimensional Hamiltonian, which is regarded as a prototype of isomerization reactions, H = (p + p + p 3)+a q + a q 4 + (! q +! q ) + X i=0 e ff i(q q i ) (fi i q q 3 + fl i q (q + q 3)): (5)
4 4 Dynamical Hierarchy in Transition States of Reactions -ij" >0 p" -ij" =0 (W + ) p" N -ij" <0 -ij" <0 q"... q'' N -ij" >0 -ij" =0 (W - ) Figure. The geometrical structure of the multidimensional phase space in the region of index one saddle (Eq. 4) and NHIM M, its stale/unstale invariant manifolds W /W + and the `local' no-return TS S. M and W =W +, respectively, correspond to the origin, and the stale/unstale directions represented y p = q =0,andp = j! jq on the ) plane with the same `internal' structure composed of the ath dof represented (q ;p y H (J = 0; p ; q ) = E (the dash ox in the figure). S for forward/ackward reactions are denoted y the two old lines along q = 0 on the (q ;p ) plane (0 <p p ;max =p ;min» p < 0) with the same internal structure of the ath dof represented y H (J = p =(! ); p ; q ) = E. As in the case of the full normal form, one can classify reactive/nonreactive trajectories in terms of the sign of ij [8].» The parameters were chosen as follows: a = 35=75;a ==875;! =;! 3 = 0:809;q 0 = ;q = q = 4:4;ff 0 = =6;ff = ff = ;fi 0 = 8;fi = fi =;fl 0 =0:75;fl = fl =. The ratio etween! and! 3, approximately the golden mean, was chosen as eing avoided from linear resonance. The imaginary frequency associated with q at the saddle is estimated as! ' 0:94i. The nonzero value of q 0 aims at avoiding specific symmetry of the potential energy function in q. The original Hamiltonian Eq. (5) is transformed to the full normal form H 0 and the partial normal form H up to the fifteenth order. In Fig., we display the full and partial normal form actions of the reactive dof J 0 (p; q) and J (p; q) evaluated y the full and partial LCPTs along three distinct classes of trajectories (namely, on a torus, in a stochastic layer and in a gloal chaotic sea) on the NHIM oeying the original Hamiltonian H(p; q) at E = +0: aove the threshold, at which the gloal chaotic regime appears. It is found that while the full normal form action J 0 (p; q) only fairly persists as invariant, the partial normal form action J (p; q) strongly persists as invariant, especially at the gloal chaotic regime (see the difference of the scale of the vertical axis and compare J with J 0 at t!±large in Fig. ). For the full normal form H 0 at slightly aove the threshold, the motions of the ath dof are integrale with constants of actions J 0 corresponding to good quantum numers in the TS. It is expected that this yields state selectivity to give rise to reaction channels linking to specific product states. However, for the partial normal form H eing nonintegrale (except only J )atthe higher
5 C.B. Li et al. 5 (a) -ij' (-ij" ) () -ij' (-ij" ) (c) -ij' (-ij" ) t t t Figure. J(p; 0 q) (dashed line) and J (p; q) (solid line) at E =+0: aove the threshold (a) on a torus, () in a stochastic layer and (c) in a gloal chaotic sea on the NHIM. The insets are the Poincaré surface of section defined y q 3 = 0 and p 3 > 0. The initial conditions were prepared y q (p; q) =p (p; q) =0attimet = 0 and the trajectories were propagated ackward and forward in time y the original Hamiltonian H (eventually at t!±large the system leaves from the (approximate) NHIM ecause of the finiteness of the order of LCPTs). Note that at E ' +0:5 the normal hyperolicity of the NHIM starts to e ruined due to the nonlinear resonances etween the reactive mode q and some of the hyperolic orits densely distriuted on the NHIM, whose imaginary frequencies are (near-)commensurale with! [9]. energy regime, one may not expect the state selectivity. Fig. 3 exemplifies this: at an energy regime slightly aove the threshold where most of tori remain on the NHIM, the system still `rememers' the initial distriution of the normal mode action J in the reactant and results in the `structured' J -distriution in the product state, where the different values of J at the initial state correspond to the well-defined distinct regions of J at the final state. As the energy increases more, the system having different J at the initial state can access any region of J at the final state, yielding a `non-structured' J -distriution. Note however that, as far as the original Hamiltonian can e transformed into the partial normal form (where normal hyperolicity holds), the dynamical correlation along the reaction coordinate q 3 Outlook and Conclusions (p; q) still persist with the invariant of action J (p; q). The hierarchical structure of the phase space was presented in the region of index one saddle y introducing a partial LCPT that normalizes solely the reactive dof. Quite recently, the quantum counterpart of this generic feature for Hamiltonian systems has also een explored [, ]. Acknowledgement. Parts of this work were supported y JSPS, Grant-in-Aid for Research on Priority Area `Control of Molecules in Intense Laser Fields' of the Ministry of Education, Science, Sports and Culture of Japan, Japan Science and Technology and st century COE (Center Of Excellence) of `Origin and Evolution of Planetary Systems (Koe Univ.)', MEXT. References. Toda, M. et al. (eds.): Adv. Chem. Phys. 30 (5) and references therein.
6 6 Dynamical Hierarchy in Transition States of Reactions P (a) 0. P () J J Figure 3. The product distriution of the normal mode action J at E = +0:05 (a) and +0: () aove the threshold. The initial conditions are prepared as follows: we turned off the coupling among the modes in the reactant y putting ff = in Eq. (5) in order to prepare awell-defined initial condition in the reactant state. At q = 4:4, we fixed p to confirm the crossings P to the product. Then, we set an initial distriution of the action J such that 3 P ini(j ) = i= ffi(j Ji) with three initial values of J indicated y three arrows (J3 is automatically determined with the fixed J and J at a given E) with the uniform sampling for the angle variale (and 3). Then, the ensemle of trajectories of the initial conditions were evolved using the original Hamiltonian untiltheyreach the product at q =4:4. The inset are the Poincaré surface of section defined y q 3 = 0 and p 3 > 0 on the NHIM (See the figure caption of Fig. in detail).. Komatsuzaki, T., Berry, R. S.: J. Chem. Phys. 0, 960 (999); Phys. Chem. Chem. Phys., 387 (999); Proc. Natl. Acad. Sci. USA 78, 7666 (); J. Chem. Phys. 5, 405 (); J. Chem. Phys. 6, 86 (); J. Phys. Chem. A 06, 0945 (); Adv. Chem. Phys. 3, 79 (); Adv. Chem. Phys. 30A, 43 (5) 3. Toda, M.: Phys. Rev. Lett. 74, 670 (995); Phys. Lett. A 7, 3 (997); Adv. Chem. Phys. 3, 53 (); Adv. Chem. Phys. 30A, 337 (5) 4. Koon, W.S. et al.: Chaos 0, 47 (0) 5. Wiggins, S. et al.: Phys. Rev. Lett. 86, 5478 () 6. Uzer, T. et al.: Nonlinearity 5, 957 () 7. Waalkens, H. et al.: J. Chem. Phys., 607 (4) 8. Li, C.B. et al.: J. Chem. Phys. 3, 8430 (5) 9. Li, C.B. et al.: (unpulished) 0. Moser, J.: Comm. Pure Appl. Math., 57 (958). Creagh, S.C.: Nonlinearity 7, 6 (4); Nonlinearity 8, 089 (5). Cargo, M. et al.: J. Phys. A 38, 977(5)
Fractional behavior in multidimensional Hamiltonian systems describing reactions
PHYSICAL REVIEW E 76, 565 7 Fractional behavior in multidimensional Hamiltonian systems describing reactions Akira Shojiguchi,, * Chun-Biu Li,, Tamiki Komatsuzaki,, and Mikito Toda Physics Department,
More informationDepartment of Chemistry, The University of Chicago, Chicago, IL 60637
Appendix: Overview of our methodology in "Regularity Non-recrossing Paths in Transition States of Chemical Reactions" Tamiki Komatsuzaki R. Stephen Berry Department of Chemistry, The University of Chicago,
More informationPhysics 106b: Lecture 7 25 January, 2018
Physics 106b: Lecture 7 25 January, 2018 Hamiltonian Chaos: Introduction Integrable Systems We start with systems that do not exhibit chaos, but instead have simple periodic motion (like the SHO) with
More informationMelnikov s Method Applied to a Multi-DOF Ship Model
Proceedings of the th International Ship Staility Workshop Melnikov s Method Applied to a Multi-DOF Ship Model Wan Wu, Leigh S. McCue, Department of Aerospace and Ocean Engineering, Virginia Polytechnic
More information= 0. = q i., q i = E
Summary of the Above Newton s second law: d 2 r dt 2 = Φ( r) Complicated vector arithmetic & coordinate system dependence Lagrangian Formalism: L q i d dt ( L q i ) = 0 n second-order differential equations
More information3.14. The model of Haldane on a honeycomb lattice
4 Phys60.n..7. Marginal case: 4 t Dirac points at k=(,). Not an insulator. No topological index...8. case IV: 4 t All the four special points has z 0. We just use u I for the whole BZ. No singularity.
More informationPADEU. The 1 : 2 resonance in the Sitnikov problem. 1 Introduction. T. Kovács 1
PADEU PADEU 7, (6) ISBN 963 463 557, ISSN 38-43 c Pulished y the Astron. Dept. of the Eötvös Univ. The : resonance in the Sitnikov prolem T. Kovács Eötvös University, Department of Astronomy, H-58 Budapest,
More informationHamiltonian Dynamics
Hamiltonian Dynamics CDS 140b Joris Vankerschaver jv@caltech.edu CDS Feb. 10, 2009 Joris Vankerschaver (CDS) Hamiltonian Dynamics Feb. 10, 2009 1 / 31 Outline 1. Introductory concepts; 2. Poisson brackets;
More informationMass-Related Dynamical Barriers in Triatomic Reactions
Few-Body Systems, 5 (5) Few- Body Systems c by Springer-Verlag 5 Printed in Austria Mass-Related Dynamical Barriers in Triatomic Reactions T. Yanao, W. S. Koon, and J. E. Marsden Control and Dynamical
More informationA Glance at the Standard Map
A Glance at the Standard Map by Ryan Tobin Abstract The Standard (Chirikov) Map is studied and various aspects of its intricate dynamics are discussed. Also, a brief discussion of the famous KAM theory
More informationHow Does the Smaller Alignment Index (SALI) Distinguish Order from Chaos?
Progress of Theoretical Physics Supplement No. 150, 2003 439 How Does the Smaller Alignment Index (SALI) Distinguish Order from Chaos? Charalampos Skokos, 1,2, ) Chris Antonopoulos, 1, ) Tassos C. Bountis
More informationThe Transition to Chaos
Linda E. Reichl The Transition to Chaos Conservative Classical Systems and Quantum Manifestations Second Edition With 180 Illustrations v I.,,-,,t,...,* ', Springer Dedication Acknowledgements v vii 1
More informationA Two-dimensional Mapping with a Strange Attractor
Commun. math. Phys. 50, 69 77 (1976) Communications in Mathematical Physics by Springer-Verlag 1976 A Two-dimensional Mapping with a Strange Attractor M. Henon Observatoire de Nice, F-06300 Nice, France
More informationInvariant Manifolds and Transport in the Three-Body Problem
Dynamical S C C A L T E C H Invariant Manifolds and Transport in the Three-Body Problem Shane D. Ross Control and Dynamical Systems California Institute of Technology Classical N-Body Systems and Applications
More informationAveraging II: Adiabatic Invariance for Integrable Systems (argued via the Averaging Principle)
Averaging II: Adiabatic Invariance for Integrable Systems (argued via the Averaging Principle In classical mechanics an adiabatic invariant is defined as follows[1]. Consider the Hamiltonian system with
More information898 IEEE TRANSACTIONS ON ROBOTICS AND AUTOMATION, VOL. 17, NO. 6, DECEMBER X/01$ IEEE
898 IEEE TRANSACTIONS ON ROBOTICS AND AUTOMATION, VOL. 17, NO. 6, DECEMBER 2001 Short Papers The Chaotic Mobile Robot Yoshihiko Nakamura and Akinori Sekiguchi Abstract In this paper, we develop a method
More informationMaslov indices and monodromy
Loughborough University Institutional Repository Maslov indices and monodromy This item was submitted to Loughborough University's Institutional Repository by the/an author. Additional Information: This
More informationTheory and Computation of Non-RRKM Reaction Rates in Chemical Systems with 3 or More D.O.F. Wang Sang Koon Control and Dynamical Systems, Caltech
Theory and Computation of Non-RRKM Reaction Rates in Chemical Systems with 3 or More D.O.F. Gabern, Koon, Marsden, and Ross Wang Sang Koon Control and Dynamical Systems, Caltech koon@cds.caltech.edu Outline
More informationLocal regularity and non-recrossing path in transition state a new strategy in chemical reaction theories
Journal of Molecular Structure (Theochem) 506 (2000) 55 70 www.elsevier.nl/locate/theochem Local regularity and non-recrossing path in transition state a new strategy in chemical reaction theories T. Komatsuzaki*,
More informationGLOBAL STABILITY FOR A CLASS OF DISCRETE SIR EPIDEMIC MODELS. Yoichi Enatsu and Yukihiko Nakata. Yoshiaki Muroya. (Communicated by Yasuhiro Takeuchi)
MATHEMATICAL BIOSCIENCES doi:10.3934/me.2010.7.347 AND ENGINEERING Volume 7, Numer 2, April 2010 pp. 347 361 GLOBAL STABILITY FOR A CLASS OF DISCRETE SIR EPIDEMIC MODELS Yoichi Enatsu and Yukihiko Nakata
More informationOn the convergence of normal forms for analytic control systems
Chapter 1 On the convergence of normal forms for analytic control systems Wei Kang Department of Mathematics, Naval Postgraduate School Monterey, CA 93943 wkang@nps.navy.mil Arthur J. Krener Department
More informationChaotic motion. Phys 750 Lecture 9
Chaotic motion Phys 750 Lecture 9 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t =0to
More informationBarcelona, Spain. RTBP, collinear points, periodic orbits, homoclinic orbits. Resumen
XX Congreso de Ecuaciones Diferenciales y Aplicaciones X Congreso de Matemática Aplicada Sevilla, 24-28 septiembre 27 (pp. 1 8) The dynamics around the collinear point L 3 of the RTBP E. Barrabés 1, J.M.
More informationA formula to compute the microcanonical volume of reactive initial conditions in transition state theory Waalkens, H.; Burbanks, A.; Wiggins, S.
University of Groningen A formula to compute the microcanonical volume of reactive initial conditions in transition state theory Waalkens, H.; Burbanks, A.; Wiggins, S. Published in: Journal of Physics
More informationControl of chaos in Hamiltonian systems
Control of chaos in Hamiltonian systems G. Ciraolo, C. Chandre, R. Lima, M. Vittot Centre de Physique Théorique CNRS, Marseille M. Pettini Osservatorio Astrofisico di Arcetri, Università di Firenze Ph.
More informationHamiltonian Dynamics In The Theory of Abstraction
Hamiltonian Dynamics In The Theory of Abstraction Subhajit Ganguly. email: gangulysubhajit@indiatimes.com Abstract: This paper deals with fluid flow dynamics which may be Hamiltonian in nature and yet
More informationOn the Third Order Rational Difference Equation
Int. J. Contemp. Math. Sciences, Vol. 4, 2009, no. 27, 32-334 On the Third Order Rational Difference Equation x n x n 2 x (a+x n x n 2 ) T. F. Irahim Department of Mathematics, Faculty of Science Mansoura
More informationSome Collision solutions of the rectilinear periodically forced Kepler problem
Advanced Nonlinear Studies 1 (2001), xxx xxx Some Collision solutions of the rectilinear periodically forced Kepler problem Lei Zhao Johann Bernoulli Institute for Mathematics and Computer Science University
More informationCoherent state semiclassical initial value representation for the Boltzmann operator in thermal correlation functions
JOURNAL OF CHEMICAL PHYSICS VOLUME 116, NUMBER 21 1 JUNE 2002 Coherent state semiclassical initial value representation for the Boltzmann operator in thermal correlation functions Nancy Makri Departments
More informationSymplectic maps. James D. Meiss. March 4, 2008
Symplectic maps James D. Meiss March 4, 2008 First used mathematically by Hermann Weyl, the term symplectic arises from a Greek word that means twining or plaiting together. This is apt, as symplectic
More informationarxiv:physics/ v1 [physics.plasm-ph] 7 Apr 2006
Europhysics Letters PREPRINT arxiv:physics/6455v [physics.plasm-ph] 7 Apr 26 Olique electromagnetic instailities for an ultra relativistic electron eam passing through a plasma A. Bret ETSI Industriales,
More informationMajor Concepts Lecture #11 Rigoberto Hernandez. TST & Transport 1
Major Concepts Onsager s Regression Hypothesis Relaxation of a perturbation Regression of fluctuations Fluctuation-Dissipation Theorem Proof of FDT & relation to Onsager s Regression Hypothesis Response
More informationIntroduction to Applied Nonlinear Dynamical Systems and Chaos
Stephen Wiggins Introduction to Applied Nonlinear Dynamical Systems and Chaos Second Edition With 250 Figures 4jj Springer I Series Preface v L I Preface to the Second Edition vii Introduction 1 1 Equilibrium
More informationEffect of Uniform Horizontal Magnetic Field on Thermal Instability in A Rotating Micropolar Fluid Saturating A Porous Medium
IOSR Journal of Mathematics (IOSR-JM) e-issn: 78-578, p-issn: 39-765X. Volume, Issue Ver. III (Jan. - Fe. 06), 5-65 www.iosrjournals.org Effect of Uniform Horizontal Magnetic Field on Thermal Instaility
More informationIMA Preprint Series # 2269
GEOMETRICAL MODELS OF THE PHASE SPACE STRUCTURES GOVERNING REACTION DYNAMICS By Holger Waalkens and Stephen Wiggins IMA Preprint Series # 2269 ( August 2009 ) INSTITUTE FOR MATHEMATICS AND ITS APPLICATIONS
More informationRATIONAL EXPECTATIONS AND THE COURNOT-THEOCHARIS PROBLEM
RATIONAL EXPECTATIONS AND THE COURNOT-THEOCHARIS PROBLEM TÖNU PUU Received 18 April 006; Accepted 1 May 006 In dynamic models in economics, often rational expectations are assumed. These are meant to show
More informationDynamical Systems Solutions to Exercises
Dynamical Systems Part 5-6 Dr G Bowtell Dynamical Systems Solutions to Exercises. Figure : Phase diagrams for i, ii and iii respectively. Only fixed point is at the origin since the equations are linear
More informationQUALITY CONTROL OF WINDS FROM METEOSAT 8 AT METEO FRANCE : SOME RESULTS
QUALITY CONTROL OF WINDS FROM METEOSAT 8 AT METEO FRANCE : SOME RESULTS Christophe Payan Météo France, Centre National de Recherches Météorologiques, Toulouse, France Astract The quality of a 30-days sample
More informationChaos and Dynamical Systems
Chaos and Dynamical Systems y Megan Richards Astract: In this paper, we will discuss the notion of chaos. We will start y introducing certain mathematical concepts needed in the understanding of chaos,
More informationSelf-Replication, Self-Destruction, and Spatio-Temporal Chaos in the Gray-Scott Model
Letter Forma, 15, 281 289, 2000 Self-Replication, Self-Destruction, and Spatio-Temporal Chaos in the Gray-Scott Model Yasumasa NISHIURA 1 * and Daishin UEYAMA 2 1 Laboratory of Nonlinear Studies and Computations,
More informationARTICLES. Normal-mode analysis without the Hessian: A driven molecular-dynamics approach
JOURNAL OF CHEMICAL PHYSICS VOLUME 119, NUMBER 2 8 JULY 2003 ARTICLES Normal-mode analysis without the Hessian: A driven molecular-dynamics approach Joel M. Bowman, a) Xiubin Zhang, and Alex Brown Cherry
More informationChaotic motion. Phys 420/580 Lecture 10
Chaotic motion Phys 420/580 Lecture 10 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t
More informationGeometrical Models of the Phase Space Structures Governing Reaction Dynamics Waalkens, Holger; Wiggins, S.
University of Groningen Geometrical Models of the Phase Space Structures Governing Reaction Dynamics Waalkens, Holger; Wiggins, S. Published in: Regular & chaotic dynamics DOI: 10.1134/S1560354710010016
More informationREVIEW. Hamilton s principle. based on FW-18. Variational statement of mechanics: (for conservative forces) action Equivalent to Newton s laws!
Hamilton s principle Variational statement of mechanics: (for conservative forces) action Equivalent to Newton s laws! based on FW-18 REVIEW the particle takes the path that minimizes the integrated difference
More informationNON-STATIONARY RESONANCE DYNAMICS OF THE HARMONICALLY FORCED PENDULUM
CYBERNETICS AND PHYSICS, VOL. 5, NO. 3, 016, 91 95 NON-STATIONARY RESONANCE DYNAMICS OF THE HARMONICALLY FORCED PENDULUM Leonid I. Manevitch Polymer and Composite Materials Department N. N. Semenov Institute
More informationCanonical transformations (Lecture 4)
Canonical transformations (Lecture 4) January 26, 2016 61/441 Lecture outline We will introduce and discuss canonical transformations that conserve the Hamiltonian structure of equations of motion. Poisson
More informationTime-Dependent Statistical Mechanics 5. The classical atomic fluid, classical mechanics, and classical equilibrium statistical mechanics
Time-Dependent Statistical Mechanics 5. The classical atomic fluid, classical mechanics, and classical equilibrium statistical mechanics c Hans C. Andersen October 1, 2009 While we know that in principle
More informationarxiv:nlin/ v1 [nlin.cd] 17 Jan 2002
Symbolic Dynamics of Homoclinic Orbits in a Symmetric Map arxiv:nlin/0201031v1 [nlin.cd] 17 Jan 2002 Zai-Qiao Bai and Wei-Mou Zheng Institute of Theoretical Physics, Academia Sinica, Beijing 100080, China
More informationForward-Backward Semiclassical Dynamics with Linear Scaling
J. Phys. Chem. A 1999, 103, 9479-9486 9479 Forward-Backward Semiclassical Dynamics with Linear Scaling Jiushu Shao School of Chemical Sciences, UniVersity of Illinois, 601 S. Goodwin AVenue, Urbana, Illinois
More informationMolecular Vibrations
Molecular Vibrations K. Srihari Department of Chemistry IIT Kanpur 6th March 2007 Chemical Reactions Make/break chemical bonds Rates: How fast? Calculate? Mechanism: Why specific bond(s) break? Control?
More informationSPLITTING OF SEPARATRICES FOR (FAST) QUASIPERIODIC FORCING. splitting, which now seems to be the main cause of the stochastic behavior in
SPLITTING OF SEPARATRICES FOR (FAST) QUASIPERIODIC FORCING A. DELSHAMS, V. GELFREICH, A. JORBA AND T.M. SEARA At the end of the last century, H. Poincare [7] discovered the phenomenon of separatrices splitting,
More informationInfluence of moving breathers on vacancies migration
Influence of moving reathers on vacancies migration J Cuevas a,1, C Katerji a, JFR Archilla a, JC Eileck, FM Russell a Grupo de Física No Lineal. Departamento de Física Aplicada I. ETSI Informática. Universidad
More informationTowards stability results for planetary problems with more than three bodies
Towards stability results for planetary problems with more than three bodies Ugo Locatelli [a] and Marco Sansottera [b] [a] Math. Dep. of Università degli Studi di Roma Tor Vergata [b] Math. Dep. of Università
More informationWHITE NOISE APPROACH TO FEYNMAN INTEGRALS. Takeyuki Hida
J. Korean Math. Soc. 38 (21), No. 2, pp. 275 281 WHITE NOISE APPROACH TO FEYNMAN INTEGRALS Takeyuki Hida Abstract. The trajectory of a classical dynamics is detrmined by the least action principle. As
More informationSecular and oscillatory motions in dynamical systems. Henk Broer Johann Bernoulli Instituut voor Wiskunde en Informatica Rijksuniversiteit Groningen
Secular and oscillatory motions in dynamical systems Henk Broer Johann Bernoulli Instituut voor Wiskunde en Informatica Rijksuniversiteit Groningen Contents 1. Toroidal symmetry 2. Secular (slow) versus
More informationIgor A. Khovanov,Vadim S. Anishchenko, Astrakhanskaya str. 83, Saratov, Russia
Noise Induced Escape from Dierent Types of Chaotic Attractor Igor A. Khovanov,Vadim S. Anishchenko, Dmitri G. Luchinsky y,andpeter V.E. McClintock y Department of Physics, Saratov State University, Astrakhanskaya
More informationSelective quasienergies from short time cross-correlation probability amplitudes by the filter-diagonalization method
PHYSICAL REVIEW E VOLUME 58, NUMBER 1 JULY 1998 Selective quasienergies from short time cross-correlation probability amplitudes by the filter-diagonalization method Markus Glück, 1 H. Jürgen Korsch, 1
More informationAn introduction to Birkhoff normal form
An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an
More informationAlignment indices: a new, simple method for determining the ordered or chaotic nature of orbits
INSTITUTE OF PHYSICSPUBLISHING JOURNAL OFPHYSICSA: MATHEMATICAL AND GENERAL J. Phys. A: Math. Gen. 34 (2001) 10029 10043 PII: S0305-4470(01)25097-5 Alignment indices: a new, simple method for determining
More informationMixing Quantum and Classical Mechanics: A Partially Miscible Solution
Mixing Quantum and Classical Mechanics: A Partially Miscible Solution R. Kapral S. Nielsen A. Sergi D. Mac Kernan G. Ciccotti quantum dynamics in a classical condensed phase environment how to simulate
More informationRepresentation theory of SU(2), density operators, purification Michael Walter, University of Amsterdam
Symmetry and Quantum Information Feruary 6, 018 Representation theory of S(), density operators, purification Lecture 7 Michael Walter, niversity of Amsterdam Last week, we learned the asic concepts of
More informationAnn. Acad. Rom. Sci. Ser. Math. Appl. Vol. 8, No. 1/2016. The space l (X) Namita Das
ISSN 2066-6594 Ann. Acad. Rom. Sci. Ser. Math. Appl. Vol. 8, No. 1/2016 The space l (X) Namita Das Astract In this paper we have shown that the space c n n 0 cannot e complemented in l n n and c 0 (H)
More informationROUNDOFF ERRORS; BACKWARD STABILITY
SECTION.5 ROUNDOFF ERRORS; BACKWARD STABILITY ROUNDOFF ERROR -- error due to the finite representation (usually in floatingpoint form) of real (and complex) numers in digital computers. FLOATING-POINT
More informationTitleQuantum Chaos in Generic Systems.
TitleQuantum Chaos in Generic Systems Author(s) Robnik, Marko Citation 物性研究 (2004), 82(5): 662-665 Issue Date 2004-08-20 URL http://hdl.handle.net/2433/97885 Right Type Departmental Bulletin Paper Textversion
More informationPentahedral Volume, Chaos, and Quantum Gravity
Pentahedral Volume, Chaos, and Quantum Gravity Hal Haggard May 30, 2012 Volume Polyhedral Volume (Bianchi, Doná and Speziale): ˆV Pol = The volume of a quantum polyhedron Outline 1 Pentahedral Volume 2
More informationAMADEU DELSHAMS AND RAFAEL RAMíREZ-ROS
POINCARÉ-MELNIKOV-ARNOLD METHOD FOR TWIST MAPS AMADEU DELSHAMS AND RAFAEL RAMíREZ-ROS 1. Introduction A general theory for perturbations of an integrable planar map with a separatrix to a hyperbolic fixed
More informationCurves in the configuration space Q or in the velocity phase space Ω satisfying the Euler-Lagrange (EL) equations,
Physics 6010, Fall 2010 Hamiltonian Formalism: Hamilton s equations. Conservation laws. Reduction. Poisson Brackets. Relevant Sections in Text: 8.1 8.3, 9.5 The Hamiltonian Formalism We now return to formal
More informationTIME-FREQUENCY ANALYSIS OF CLASSICAL TRAJECTORIES OF POLYATOMIC MOLECULES
International Journal of Bifurcation and Chaos, Vol. 11, No. 5 (1) 159 18 c World Scientific Publishing Company TIME-FREQUENCY ANALYSIS OF CLASSICAL TRAJECTORIES OF POLYATOMIC MOLECULES LUZ V. VELA-AREVALO
More informationSemi-Classical Theory for Non-separable Systems :
Semi-Classical Theory for Non-separable Systems : Construction of Good Action-Angle Variables for Reaction Rate Constants* BY WILLIAM H. MILLER-^ Department of Chemistry, and Materials and Molecular Research
More information138. Last Latexed: April 25, 2017 at 9:45 Joel A. Shapiro. so 2. iψ α j x j
138. Last Latexed: April 25, 2017 at 9:45 Joel A. Shapiro The Hamiltonian for a noninteracting fermion 1 is H = d 3 x H 1, with H 1 = iψ j α ψ j + mψ j βψ j j Chapter 13 Local Symmetry So far, we have
More informationLecture 11 : Overview
Lecture 11 : Overview Error in Assignment 3 : In Eq. 1, Hamiltonian should be H = p2 r 2m + p2 ϕ 2mr + (p z ea z ) 2 2 2m + eφ (1) Error in lecture 10, slide 7, Eq. (21). Should be S(q, α, t) m Q = β =
More informationStructure of near-threshold s-wave resonances
Structure of near-threshold s-wave resonances Tetsuo Hyodo Yukawa Institute for Theoretical Physics, Kyoto 203, Sep. 0th Introduction Structure of hadron excited states Various excitations of baryons M
More informationAkinori Sekiguchi and Yoshihiko Nakamura. Dept. of Mechano-Informatics, University of Tokyo Hongo, Bunkyo-Ku, Tokyo , Japan
The Chaotic Mobile Robot Akinori Sekiguchi and Yoshihiko Nakamura Dept. of Mechano-Informatics, University of Tokyo 7-- Hongo, Bunkyo-Ku, Tokyo -866, Japan ABSTRACT In this paper, we develop a method to
More informationBINARY ASTEROID PAIRS A Systematic Investigation of the Full Two-Body Problem
BINARY ASTEROID PAIRS A Systematic Investigation of the Full Two-Body Problem Michael Priolo Jerry Marsden, Ph.D., Mentor Shane Ross, Graduate Student, Co-Mentor Control and Dynamical Systems 107-81 Caltech,
More informationDynamics on the Double Morse Potential: A Paradigm for Roaming Reactions with no Saddle Points
Dynamics on the Double Morse Potential: A Paradigm for Roaming Reactions with no Saddle Points Barry K. Carpenter,, Gregory S. Ezra, Stavros C. Farantos, Zeb C. Kramer, Stephen Wiggins* School of Chemistry,
More informationON THE BREAK-UP OF INVARIANT TORI WITH THREE FREQUENCIES
ON THE BREAK-UP OF INVARIANT TORI WITH THREE FREQUENCIES J.D. MEISS Program in Applied Mathematics University of Colorado Boulder, CO Abstract We construct an approximate renormalization operator for a
More informationTravel Grouping of Evaporating Polydisperse Droplets in Oscillating Flow- Theoretical Analysis
Travel Grouping of Evaporating Polydisperse Droplets in Oscillating Flow- Theoretical Analysis DAVID KATOSHEVSKI Department of Biotechnology and Environmental Engineering Ben-Gurion niversity of the Negev
More informationNoncollision Singularities in the n-body Problem
Noncollision Singularities in the n-body Problem Danni Tu Department of Mathematics, Princeton University January 14, 2013 The n-body Problem Introduction Suppose we have n points in R 3 interacting through
More informationThe Ermakov Lewis invariants for coupled linear oscillators
J Phys A: Math Gen 31 1998) L279 L285 Printed in the UK PII: S0305-447098)90563-7 LETTER TO THE EDITOR The Ermakov Lewis invariants for coupled linear oscillators Karl-Erik Thylwe and H Jürgen Korsch Department
More informationIm + α α. β + I 1 I 1< 0 I 1= 0 I 1 > 0
ON THE HAMILTONIAN ANDRONOV-HOPF BIFURCATION M. Olle, J. Villanueva 2 and J. R. Pacha 3 2 3 Departament de Matematica Aplicada I (UPC), Barcelona, Spain In this contribution, we consider an specic type
More informationThe Breakdown of KAM Trajectories
The Breakdown of KAM Trajectories D. BENSIMONO ~t L. P. KADANOFF~ OAT& T Bell Laboraiories Murray Hill,New Jersey 07974 bthe James Franck Insrirure Chicago, Illinois 60637 INTRODUCl'ION Hamiltonian systems
More informationA Symmetric Treatment of Damped Harmonic Oscillator in Extended Phase Space
Proceedings of Institute of Mathematics of NAS of Ukraine 00, Vol. 43, Part, 645 65 A Symmetric Treatment of Damped Harmonic Oscillator in Extended Phase Space S. NASIRI and H. SAFARI Institute for Advanced
More informationThe Structure of Hyperbolic Sets
The Structure of Hyperbolic Sets p. 1/35 The Structure of Hyperbolic Sets Todd Fisher tfisher@math.umd.edu Department of Mathematics University of Maryland, College Park The Structure of Hyperbolic Sets
More informationHamiltonian flow in phase space and Liouville s theorem (Lecture 5)
Hamiltonian flow in phase space and Liouville s theorem (Lecture 5) January 26, 2016 90/441 Lecture outline We will discuss the Hamiltonian flow in the phase space. This flow represents a time dependent
More information1 Hamiltonian formalism
1 Hamiltonian formalism 1.1 Hamilton s principle of stationary action A dynamical system with a finite number n degrees of freedom can be described by real functions of time q i (t) (i =1, 2,..., n) which,
More informationHamiltonian Systems and Chaos Overview Liz Lane-Harvard, Melissa Swager
Hamiltonian Systems and Chaos Overview Liz Lane-Harvard, Melissa Swager Abstract: In this paper we will give an overview of Hamiltonian systems with specific examples, including the classical pendulum
More informationPath Integral for Spin
Path Integral for Spin Altland-Simons have a good discussion in 3.3 Applications of the Feynman Path Integral to the quantization of spin, which is defined by the commutation relations [Ŝj, Ŝk = iɛ jk
More informationChaos in the Hénon-Heiles system
Chaos in the Hénon-Heiles system University of Karlstad Christian Emanuelsson Analytical Mechanics FYGC04 Abstract This paper briefly describes how the Hénon-Helies system exhibits chaos. First some subjects
More informationCanonical Forms for BiHamiltonian Systems
Canonical Forms for BiHamiltonian Systems Peter J. Olver Dedicated to the Memory of Jean-Louis Verdier BiHamiltonian systems were first defined in the fundamental paper of Magri, [5], which deduced the
More informationA Toy Model of Consciousness and the Connectome, First Approximation, Second Part
A Toy Model of Consciousness and the Connectome, First Approximation, Second Part Aility Potential of A Neural Network As A Prognostic Connective Indicator Stan Bumle, Ph. D. Astract: The paper vixra:1305.0128
More informationOn complexified quantum mechanics and space-time
On complexified quantum mechanics and space-time Dorje C. Brody Mathematical Sciences Brunel University, Uxbridge UB8 3PH dorje.brody@brunel.ac.uk Quantum Physics with Non-Hermitian Operators Dresden:
More information1.3 Translational Invariance
1.3. TRANSLATIONAL INVARIANCE 7 Version of January 28, 2005 Thus the required rotational invariance statement is verified: [J, H] = [L + 1 Σ, H] = iα p iα p = 0. (1.49) 2 1.3 Translational Invariance One
More informationHandout 10. Applications to Solids
ME346A Introduction to Statistical Mechanics Wei Cai Stanford University Win 2011 Handout 10. Applications to Solids February 23, 2011 Contents 1 Average kinetic and potential energy 2 2 Virial theorem
More informationAuthor(s) Kawashima, Maoki; Miyashita, Seiji;
Title Quantum Phase Transition of Heisenberg Antiferromagnet Two-Dim Author(s) Todo, Synge; Yasuda, Chitoshi; Kato Kawashima, Maoki; Miyashita, Seiji; Citation Progress of Theoretical Physics Sup 512 Issue
More information= w 2. w 1. B j. A j. C + j1j2
Local Minima and Plateaus in Multilayer Neural Networks Kenji Fukumizu and Shun-ichi Amari Brain Science Institute, RIKEN Hirosawa 2-, Wako, Saitama 35-098, Japan E-mail: ffuku, amarig@brain.riken.go.jp
More information1D spirals: is multi stability essential?
1D spirals: is multi staility essential? A. Bhattacharyay Dipartimento di Fisika G. Galilei Universitá di Padova Via Marzolo 8, 35131 Padova Italy arxiv:nlin/0502024v2 [nlin.ps] 23 Sep 2005 Feruary 8,
More informationSYMPLECTIC MANIFOLDS, GEOMETRIC QUANTIZATION, AND UNITARY REPRESENTATIONS OF LIE GROUPS. 1. Introduction
SYMPLECTIC MANIFOLDS, GEOMETRIC QUANTIZATION, AND UNITARY REPRESENTATIONS OF LIE GROUPS CRAIG JACKSON 1. Introduction Generally speaking, geometric quantization is a scheme for associating Hilbert spaces
More informationNew Exact Solutions for a Class of High-order Dispersive Cubic-quintic Nonlinear Schrödinger Equation
Journal of Mathematics Research; Vol. 6, No. 4; 2014 ISSN 1916-9795 E-ISSN 1916-9809 Pulished y Canadian Center of Science and Education New Exact Solutions for a Class of High-order Dispersive Cuic-quintic
More informationVector Fields. It is standard to define F µν = µ ϕ ν ν ϕ µ, so that the action may be written compactly as
Vector Fields The most general Poincaré-invariant local quadratic action for a vector field with no more than first derivatives on the fields (ensuring that classical evolution is determined based on the
More informationQuadratic non-condon effect in optical spectra of impurity paramagnetic centers in dielectric crystals
Journal of Physics: Conference Series Quadratic non-condon effect in optical spectra of impurity paramagnetic centers in dielectric crystals To cite this article: R Yu Yunusov and O V Solovyev J Phys:
More information