BUDKER INSTITUTE OF NUCLEAR PHYSICS. B.V. Chirikov and V.V. Vecheslavov MULTIPLE SEPARATRIX CROSSING: CHAOS STRUCTURE. Budker INP NOVOSIBIRSK
|
|
- Hubert Gilmore
- 5 years ago
- Views:
Transcription
1 BUDKER INSTITUTE OF NUCLEAR PHYSICS B.V. Chirikov and V.V. Vecheslavov MULTIPLE SEPARATRIX CROSSING: CHAOS STRUCTURE Budker INP NOVOSIBIRSK 2000
2 chaos structure B.V. Chirikov and V.V. Vecheslavov Budker Institute of Nuclear Physics Novosibirsk, Russia Abstract Numerical experiments on the structure of the chaotic component of motion under multiple crossing of the separatrix of a nonlinear resonance with timeívarying amplitude are described with the main attention to the problem of ergodicity. The results clearly demonstrate nonergodicity of that motion due to the presence of a regular component of relatively small measure with a very complicated structure. A simple 2Dímap per crossing has been constructed which qualitatively describes the main properties of both chaotic and regular components of the motion. An empirical relation for the correlationíaæected diæusion rate has been found including a close vicinity of the chaos border where an evidence of the critical structure has been observed. Some unsolved problems and open questions are also discussed. chirikov@inp.nsk.su vecheslavov@inp.nsk.su cæbudker Institute of Nuclear Physics SB RAS
3 1 Introduction The present work continues the studies of chaotic motion under a slow separatrix crossing. This is a particular case of adiabatic processes which are very important in physics because of the adiabatic invariance, approximate though, that is of the conservation of action variables èjè under a slow parametric perturbation. The main problem here is the degree of accuracy or of violation of that invariance. Separatrix crossing produces the largest chaotic component in phase space whose size does not depend on the adiabatic parameter æ! 0 which, however, does aæect the detailed structure of the motion as well as its time scale. In our previous paper ë1ë the single separatrix crossing for a particular model was described in detail. Remarkably, a fairly simple relation for such a model in Ref.ë2ë we used turned out to be surprisingly accurate within the most part of the chaotic component. In this paper we discribe the results of numerical experiments on multiple separatrix crossing. We focus on statistical properties of the motion, including the structure and measure of regular component disseminated into the chaotic 'sea' in a rather tricky way. The existence of regular component means nonergodicity of the motion, the question which has remained unclear for a long time until recently. To our knowledge, the nonergodicity of motion in a similar model was ærst predicted theoretically and estimated numerically in Ref.ë3ë. We have conærmed this result by diæerent methods, and found many other characteristics of the motion structure. The present work, as well as the previous one ë1ë, was stimulated by a very interesting study of the corresponding quantum adiabaticity ë4ë. We use the same classical model which is brieæy described, for reader's convenience, in the next Section èfor details see Ref.ë1ëè. 3
4 The model is speciæed by the Hamiltonian: Hèx; p; tè = p2 2 + A 0 sin èætè æ cos x; è2:1è which describes a single nonlinear resonance in the pendulum approximation èsee, e.g., Ref.ë5, 6ëè with a timeívarying amplitude Aètè = A 0 sin èætè : è2:2è The dimensionless adiabaticity parameter is deæned in the usual way as the ratio of perturbationèoscillation frequencies: æ = æ p A0 ; è2:3è where p A 0 is constant frequency of the small pendulum oscillation for the maximal amplitude. Two branches of the instant, or 'frozen', separatrix at some t = const are given by the relation p s èx 0 ; tè =æ2 p ç ç ç x 0 x; Aètè é 0 jaètèjæsin ; x 0 = 2 x, ç; Aètè é 0 : è2:4è Following previous studies of separatrix crossing we restrict ourselves below to this frozen approximation. As was shown in Ref.ë1ë the latter provides quite good accuracy of fairly simple theoretical relations. In this approximation the action variable is deæned in the standard way as J = 1 2ç I pèxè dx ; è2:5è where integral is taken over the whole period for x rotation èoæ the resonanceè and over a half of that for x oscillation èinside the resonanceè. This diæerence is necessary to avoid the discontinuity of J at separatrix where the action is given by a simple expression J = J s ètè = 4 ç p jaètèj çjmax = 4 ç p A 0 : è2:6è At æt =0èmod çè the action J = jpj, and the conjugated phase ç = x. Notice that unlike p the action J ç 0isnever negative. 4
5 transformation: J=J max! J. Then, the crossing region swept by separatrix is the unit interval, and J is simply related to the crossing time t = t cr by the expression jaèt cr èj = J 2 ; 0 ç J ç 1 è2:7è while the adiabaticity parameter becomes æ = æ. Numerical integration of the motion equations for Hamiltonian è2.1è was performed in èx; pè variables using two algorithms. In most cases it was the so-called bilateral symplectic fourth-order Runge-Kutta algorithm as in Ref.ë1ë. However, in a few most long runs we applied a very simple and also symplectic ærstíorder algorithm like in Ref.ë2ë which is actually the well known standard map ë5ë with the timeívarying parameter: ~p = ~p + A ~ ç ç 0 æ sin ~æ t ~ æ sin x; x = x + ~p; è2:8è where tilde marks the new set of quantities rescaled by the transformation ~A 0 = 1 s 2 ; ~ t = st ; ~æ= æ s ; ~p = p s : è2:9è Here s is the scaling parameter, and we remind that A 0 = 1. The primary goal of the rescaling was decreasing parameter ~ A 0 which controls the computation accuracy. Usually, itwas around ~ A0 ç 0:1. As is well known, the variation of J under adiabatic perturbation consists of the two qualitatively diæerent parts: èiè the average action which is nearly constant between the crossings, up to an exponentially small correction, and which is of primary interest in our problem, and èiiè the rapid oscillation with the motion frequency. The ratio of the two time scales is ç æ= p jaètèj ç1 which allows for eæcient suppression of the second unimportant part of J variation by a simple averaging of Jètè over a long time interval ç 1=æ èsee Ref.ë1ëè. 3 Ergodicity The ergodicity is the weakest statistical property in dynamical systems èsee, e.g., Ref.ë7ëè. Nevertheless, it is an important characteristic of the motion, necessary in the statistical theory èsee, e.g., Ref.ë8ëè. The question of ergodicity of the motion under separatrix crossing remained open for a long time until recently. The upper bound for the measure èphaseíspace areaè of a separate domain with regular motion è'stability islet'è was estimated in Ref.ë9ë asç 1 é ç æ. 5
6 ærst predicted theoretically and estimated numerically in Ref.ë3ë. The authors directly calculated the number and positions of stable trajectories for two diæerent periods. Moreover, they were able to locate some of them in computation, and thus to measure their area in phase space which turned out to be surprisingly small. Here, we make use of a diæerent, statistical, approach. To this end, we ærst obtain from numerical experiments the steadyístate distribution f s èjè in the action. In case of ergodic motion it should be constant. Examples of the distribution are shown in Fig.1 with the parameters listed in the Table below. Table. Regular component under separatrix crossing n æ ç r æ 10 2 T æ N tr N b 1 0:1 0:68 æ 0:2 2 æ 10 3 æ :05 0:75 æ 0:06 4 æ 10 5 æ :033 0:70 æ 0:2 4 æ 10 5 æ :033 0:81 æ 0:08 4 æ 10 5 æ :02 0:60 æ 0:05 2 æ 10 6 æ :01 0:75 æ 0:04 4 æ 10 6 æ æ í parameter of adiabaticity ç r í total relative measure of regular component T í number of separatrix crossings for each ofn tr trajectories N b í number of histogram bins in Fig.1 n í reference number for Fig.1 The striking feature of all the distributions is clear and rather speciæc inhomogeneity, reminiscent of a burst of 'icicles' hanging down from a nearly 'ergodic roof'. This directly demonstrates the generic nonergodic character of motion under separatrix crossing. The histograms normalized in such a way that for ergodic motion the distribution f s èjè = 1 while the sum over all the bins is also unity for any distribution. As a result the dips in the distribution è'icicles'è, indicating the regular component, are compensated by an increase in the ergodic background. The latter is clearly seen in all distributions, especially for small J, and is a measure of the regular component. Namely, the relative measure èshareè is given by the approximate relation ç r ç éf s èjè, 1 é; J é J 1 ; è3:1è 6
7
NON-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 information9 Facta Universitatis ser.: Elect. and Energ. vol. 11, No.3 è1998è this paper we have considered shaping gain for two interesting quantization procedu
FACTA UNIVERSITATIS èniçsè Series: Electronics and Energetics vol. 11, No.3 è1998è, 91-99 NONLINEAR TRANSFORMATION OF ONEíDIMENSIONAL CONSTELLATION POINTS IN ORDER TO ERROR PROBABILITY DECREASING Zoran
More informationFigure 1: Startup screen for Interactive Physics Getting Started The ærst simulation will be very simple. With the mouse, select the circle picture on
Experiment Simulations í Kinematics, Collisions and Simple Harmonic Motion Introduction Each of the other experiments you perform in this laboratory involve a physical apparatus which you use to make measurements
More informationBrazilian Journal of Physics, vol. 28, no. 3, September, Transitions Between Periodic Orbits and Control of
Brazilian Journal of Physics, vol. 28, no. 3, September, 1998 203 Transitions Between Periodic Orbits and Control of Chaos Assisted by an External Force in Hamiltonian Systems R.M. Angelo æ, K. Furuya
More informationComputational Physics (6810): Session 8
Computational Physics (6810): Session 8 Dick Furnstahl Nuclear Theory Group OSU Physics Department February 24, 2014 Differential equation solving Session 7 Preview Session 8 Stuff Solving differential
More informationLinear and Nonlinear Oscillators (Lecture 2)
Linear and Nonlinear Oscillators (Lecture 2) January 25, 2016 7/441 Lecture outline A simple model of a linear oscillator lies in the foundation of many physical phenomena in accelerator dynamics. A typical
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 informationA New Invariance Property of Lyapunov Characteristic Directions S. Bharadwaj and K.D. Mease Mechanical and Aerospace Engineering University of Califor
A New Invariance Property of Lyapunov Characteristic Directions S. Bharadwaj and K.D. Mease Mechanical and Aerospace Engineering University of California, Irvine, California, 92697-3975 Email: sanjay@eng.uci.edu,
More informationLecture 1. Introduction. The importance, ubiquity, and complexity of embedded systems are growing
Lecture 1 Introduction Karl Henrik Johansson The importance, ubiquity, and complexity of embedded systems are growing tremendously thanks to the revolution in digital technology. This has created a need
More informationarxiv:cond-mat/ v1 29 Dec 1996
Chaotic enhancement of hydrogen atoms excitation in magnetic and microwave fields Giuliano Benenti, Giulio Casati Università di Milano, sede di Como, Via Lucini 3, 22100 Como, Italy arxiv:cond-mat/9612238v1
More informationThe total current injected in the system must be zero, therefore the sum of vector entries B j j is zero, è1;:::1èb j j =: Excluding j by means of è1è
1 Optimization of networks 1.1 Description of electrical networks Consider an electrical network of N nodes n 1 ;:::n N and M links that join them together. Each linkl pq = l k joints the notes numbered
More informationproblem of detection naturally arises in technical diagnostics, where one is interested in detecting cracks, corrosion, or any other defect in a sampl
In: Structural and Multidisciplinary Optimization, N. Olhoæ and G. I. N. Rozvany eds, Pergamon, 1995, 543í548. BOUNDS FOR DETECTABILITY OF MATERIAL'S DAMAGE BY NOISY ELECTRICAL MEASUREMENTS Elena CHERKAEVA
More informationWç;Ze èdiæerential or totalè cross sections can be written schematically as ç = è1 + A 0 P èç un + çèp + A 0 èç pol, where P is the electron's polariz
SLAC-PUB-8192 July 1999 POLARIZATION ASYMMETRIES IN æe COLLISIONS AND TRIPLE GAUGE BOSON COUPLINGS REVISITED a T.G. RIZZO b Stanford Linear Accelerator Center, Stanford University, Stanford, CA 94309,
More informationBragg scattering of an atomic beam by a standing laser wave with time-periodic amplitude modulation
PHYSICAL REVIEW A VOLUME 57, NUMBER 5 MAY 1998 Bragg scattering of an atomic beam by a standing laser wave with time-periodic amplitude modulation A. R. Kolovsky Kirensky Institute of Physics, 6636 Krasnoyarsk,
More informationQuantum Invariants: Topographical Map of Quantized Actions
University of Rhode Island DigitalCommons@URI Physics Faculty Publications Physics 1990 Quantum Invariants: Topographical Map of Quantized Actions Niraj Srivastava University of Rhode Island Gerhard Müller
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 informationSolutions of a PT-symmetric Dimer with Constant Gain-loss
Solutions of a PT-symmetric Dimer with Constant Gain-loss G14DIS Mathematics 4th Year Dissertation Spring 2012/2013 School of Mathematical Sciences University of Nottingham John Pickton Supervisor: Dr
More informationNumerical Computational Techniques for Nonlinear Optimal Control
Seminar at LAAS Numerical Computational echniques for Nonlinear Optimal Control Yasuaki Oishi (Nanzan University) July 2, 2018 * Joint work with N. Sakamoto and. Nakamura 1. Introduction Optimal control
More informationMethod of Averaging for Differential Equations on an Infinite Interval
Method of Averaging for Differential Equations on an Infinite Interval Theory and Applications SUB Gottingen 7 222 045 71X ;, ' Vladimir Burd Yaroslavl State University Yaroslavl, Russia 2 ' 08A14338 Contents
More informationPhysics 106a, Caltech 13 November, Lecture 13: Action, Hamilton-Jacobi Theory. Action-Angle Variables
Physics 06a, Caltech 3 November, 08 Lecture 3: Action, Hamilton-Jacobi Theory Starred sections are advanced topics for interest and future reference. The unstarred material will not be tested on the final
More informationalso has x æ as a local imizer. Of course, æ is typically not known, but an algorithm can approximate æ as it approximates x æ èas the augmented Lagra
Introduction to sequential quadratic programg Mark S. Gockenbach Introduction Sequential quadratic programg èsqpè methods attempt to solve a nonlinear program directly rather than convert it to a sequence
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 informationLECTURE Review. In this lecture we shall study the errors and stability properties for numerical solutions of initial value.
LECTURE 24 Error Analysis for Multi-step Methods 1. Review In this lecture we shall study the errors and stability properties for numerical solutions of initial value problems of the form è24.1è dx = fèt;
More informationSTUDY OF SYNCHRONIZED MOTIONS IN A ONE-DIMENSIONAL ARRAY OF COUPLED CHAOTIC CIRCUITS
International Journal of Bifurcation and Chaos, Vol 9, No 11 (1999) 19 4 c World Scientific Publishing Company STUDY OF SYNCHRONIZED MOTIONS IN A ONE-DIMENSIONAL ARRAY OF COUPLED CHAOTIC CIRCUITS ZBIGNIEW
More informationAdiabatic Approximation
Adiabatic Approximation The reaction of a system to a time-dependent perturbation depends in detail on the time scale of the perturbation. Consider, for example, an ideal pendulum, with no friction or
More informationAssessing the Calibration of Naive Bayes' Posterior. Estimates. Paul N. Bennett. September 12, 2000 CMU-CS Carnegie Mellon University
Assessing the Calibration of Naive Bayes' Posterior Estimates Paul N. Bennett September 2, 2 CMU-CS--55 School of Computer Science Carnegie Mellon University Pittsburgh, PA 523 Abstract In this paper,
More informationMechanics IV: Oscillations
Mechanics IV: Oscillations Chapter 4 of Morin covers oscillations, including damped and driven oscillators in detail. Also see chapter 10 of Kleppner and Kolenkow. For more on normal modes, see any book
More informationTechnical Report No AN APPLICATION OF PARALLEL COMPUTATION TO DYNAMICAL SYSTEMS æ Selim G. Akl and Weiguang Yao School of Computing Queen's
Technical Report No. 2003-466 AN APPLICATION OF PARALLEL COMPUTATION TO DYNAMICAL SYSTEMS æ Selim G. Akl and Weiguang Yao School of Computing Queen's University Kingston, Ontario, Canada K7L 3N6 June 27,
More information256 Facta Universitatis ser.: Elect. and Energ. vol. 11, No.2 è1998è primarily concerned with narrow range of frequencies near ærst resonance èwhere s
FACTA UNIVERSITATIS èni ç Sè Series: Electronics and Energetics vol. 11, No.2 è1998è, 255-261 EFFICIENT CALCULATION OF RADAR CROSS SECTION FOR FINITE STRIP ARRAY ON DIELECTRIC SLAB Borislav Popovski and
More informationQuantum Chaos as a Practical Tool in Many-Body Physics ESQGP Shuryak fest
Quantum Chaos as a Practical Tool in Many-Body Physics ESQGP Shuryak fest Vladimir Zelevinsky NSCL/ Michigan State University Stony Brook October 3, 2008 Budker Institute of Nuclear Physics, Novosibirsk
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 informationAutoresonant excitation and evolution of nonlinear waves:
PHYSICAL REVIEW E VOLUME 55, NUMBER 2 FEBRUARY 1997 Autoresonant excitation and evolution of nonlinear waves: The variational approach L. Friedland Racah Institute of Physics, Hebrew University of Jerusalem,
More informationKlaus Korossy: Linear-Recursive Number Sequence Tasks 44 to be a prime example of inductive reasoning. Inductive reasoning has been an integral part o
Methods of Psychological Research Online 1998, Vol.3, No.1 Internet: http:èèwww.pabst-publishers.deèmprè Solvability and Uniqueness of Linear-Recursive Number Sequence Tasks æ Klaus Korossy Psychologisches
More informationA plane autonomous system is a pair of simultaneous first-order differential equations,
Chapter 11 Phase-Plane Techniques 11.1 Plane Autonomous Systems A plane autonomous system is a pair of simultaneous first-order differential equations, ẋ = f(x, y), ẏ = g(x, y). This system has an equilibrium
More informationtokamak and stellarator geometry, regarding both its physical character and its interaction
THE INFLUENCE OF ZONAL EXB FLOWS ON EDGE TURBULENCE IN TOKAMAKS AND STELLARATORS B. SCOTT, F. JENKO, A. KENDL Max-Planck-Institut fíur Plasmaphysik, Garching, Germany We report on æuid, gyroæuid and gyrokinetic
More informationParticle Dynamics in Time-Dependent Stadium-Like Billiards
Journal of Statistical Physics, Vol. 108, Nos. 5/6, September 2002 ( 2002) Particle Dynamics in Time-Dependent Stadium-Like Billiards Alexander Loskutov 1 and Alexei Ryabov 1 Received November 7, 2001;
More informationThe Phase Transition 55. have a peculiar gravitation in which the probability of merging is
The Phase Transition 55 from to + d there is a probability c 1 c 2 d that they will merge. Components have a peculiar gravitation in which the probability of merging is proportional to their sizes. With
More informationQUANTUM CHAOS IN NUCLEAR PHYSICS
QUANTUM CHAOS IN NUCLEAR PHYSICS Investigation of quantum chaos in nuclear physics is strongly hampered by the absence of even the definition of quantum chaos, not to mention the numerical criterion of
More informationChapter 2 PARAMETRIC OSCILLATOR
CHAPTER- Chapter PARAMETRIC OSCILLATOR.1 Introduction A simple pendulum consists of a mass m suspended from a string of length L which is fixed at a pivot P. When simple pendulum is displaced to an initial
More informationPhysics 221B Spring 1997 Notes 31 Adiabatic Invariance, the Geometric Phase, and the Born-Oppenheimer Approximation
Physics 221B Spring 1997 Notes 31 Adiabatic Invariance, the Geometric Phase, and the Born-Oppenheimer Approximation A time-dependent system is said to be adiabatic if the time-dependence is slow. First
More informationChaos in the Planar Two-Body Coulomb Problem with a Uniform Magnetic Field
Annual Review of Chaos Theory, Bifurcations and Dynamical Systems Vol. 3, (2013) 23-33, www.arctbds.com. Copyright (c) 2013 (ARCTBDS). ISSN 2253 0371. All Rights Reserved. Chaos in the Planar Two-Body
More informationDirect observation of free fractional charge elementary particles would be an undisputed signature of physics beyond the Standard Model. In this paper
SLAC-PUB-8283 Oct, 999 Search for Free Fractional Electric Charge Elementary Particles V. Halyo, P. Kim, E. R. Lee, I. T. Lee, D. Loomba, and M. L. Perl Stanford Linear Accelerator Center Stanford University,
More informationLecture 10. Transition probabilities and photoelectric cross sections
Lecture 10 Transition probabilities and photoelectric cross sections TRANSITION PROBABILITIES AND PHOTOELECTRIC CROSS SECTIONS Cross section = = Transition probability per unit time of exciting a single
More informationBifurcations of phase portraits of pendulum with vibrating suspension point
Bifurcations of phase portraits of pendulum with vibrating suspension point arxiv:1605.09448v [math.ds] 9 Sep 016 A.I. Neishtadt 1,,, K. Sheng 1 1 Loughborough University, Loughborough, LE11 3TU, UK Space
More informationDavid Sagan. In a colliding beam storage ring the Twiss parameters are aæected by the quadrupolar
ëdcs.papers.dynamic BETAëDYNAMIC BETA.TEX CBN 94í06 April 25, 1995 The Dynamic Beta Eæect in CESR David Sagan Introduction In a colliding beam storage ring the Twiss parameters are aæected by the quadrupolar
More informationB2.III Revision notes: quantum physics
B.III Revision notes: quantum physics Dr D.M.Lucas, TT 0 These notes give a summary of most of the Quantum part of this course, to complement Prof. Ewart s notes on Atomic Structure, and Prof. Hooker s
More informationTwo-Body Problem. Central Potential. 1D Motion
Two-Body Problem. Central Potential. D Motion The simplest non-trivial dynamical problem is the problem of two particles. The equations of motion read. m r = F 2, () We already know that the center of
More informationMultiple scale methods
Multiple scale methods G. Pedersen MEK 3100/4100, Spring 2006 March 13, 2006 1 Background Many physical problems involve more than one temporal or spatial scale. One important example is the boundary layer
More informationPasta-Ulam problem. The birth of nonlinear physics and dynamical chaos. Paolo Valentini
The Fermi-Pasta Pasta-Ulam problem The birth of nonlinear physics and dynamical chaos Paolo Valentini SUMMARY The Fermi-Pasta Pasta-Ulam experiment on MANIAC at LANL in 1955 The ergodic hypothesis The
More informationA.V. SAVKIN AND I.R. PETERSEN uncertain systems in which the uncertainty satisæes a certain integral quadratic constraint; e.g., see ë5, 6, 7ë. The ad
Journal of Mathematical Systems, Estimation, and Control Vol. 6, No. 3, 1996, pp. 1í14 cæ 1996 Birkhíauser-Boston Robust H 1 Control of Uncertain Systems with Structured Uncertainty æ Andrey V. Savkin
More informationParity violation. no left-handed ν$ are produced
Parity violation Wu experiment: b decay of polarized nuclei of Cobalt: Co (spin 5) decays to Ni (spin 4), electron and anti-neutrino (spin ½) Parity changes the helicity (H). Ø P-conservation assumes a
More informationA Classical Approach to the Stark-Effect. Mridul Mehta Advisor: Prof. Enrique J. Galvez Physics Dept., Colgate University
A Classical Approach to the Stark-Effect Mridul Mehta Advisor: Prof. Enrique J. Galvez Physics Dept., Colgate University Abstract The state of an atom in the presence of an external electric field is known
More information2:2:1 Resonance in the Quasiperiodic Mathieu Equation
Nonlinear Dynamics 31: 367 374, 003. 003 Kluwer Academic Publishers. Printed in the Netherlands. ::1 Resonance in the Quasiperiodic Mathieu Equation RICHARD RAND Department of Theoretical and Applied Mechanics,
More informationDavydov Soliton Collisions
Davydov Soliton Collisions Benkui Tan Department of Geophysics Peking University Beijing 100871 People s Republic of China Telephone: 86-10-62755041 email:dqgchw@ibmstone.pku.edu.cn John P. Boyd Dept.
More informationForced Oscillations in a Linear System Problems
Forced Oscillations in a Linear System Problems Summary of the Principal Formulas The differential equation of forced oscillations for the kinematic excitation: ϕ + 2γ ϕ + ω 2 0ϕ = ω 2 0φ 0 sin ωt. Steady-state
More informationarxiv: v1 [nlin.cd] 20 Mar 2017
Non-integrability of the semiclassical Jaynes Cummings models without the rotating-wave approximation Andrzej Maciejewski Janusz Gil Institute of Astronomy University of Zielona Góra Lubuska PL-65-65 Zielona
More informationA Scaled Diæerence Chi-square Test Statistic. Albert Satorra. Universitat Pompeu Fabra. and. Peter M. Bentler. August 3, 1999
A Scaled Diæerence Chi-square Test Statistic for Moment Structure Analysis æ Albert Satorra Universitat Pompeu Fabra and Peter M. Bentler University of California, Los Angeles August 3, 1999 æ Research
More informationTorsion Spring Oscillator with Dry Friction
Torsion Spring Oscillator with Dry Friction Manual Eugene Butikov Annotation. The manual includes a description of the simulated physical system and a summary of the relevant theoretical material for students
More informationContinuity of Bçezier patches. Jana Pçlnikovça, Jaroslav Plaçcek, Juraj ç Sofranko. Faculty of Mathematics and Physics. Comenius University
Continuity of Bezier patches Jana Plnikova, Jaroslav Placek, Juraj Sofranko Faculty of Mathematics and Physics Comenius University Bratislava Abstract The paper is concerned about the question of smooth
More informationLaboratory for Astronomy & Solar Physics, Code 681, NASAèGSFC, Greenbelt MD
1997 HST Calibration Workshop Space Telescope Science Institute, 1997 S. Casertano, et al., eds. Correction for the STIS Echelle Blaze Function Sara R. Heap and Thomas M. Brown 1 Laboratory for Astronomy
More informationof the BaBar Silicon Vertex Tracker D. Barni, D. Giugni, F. Lanni, F. Palombo Abstract
BaBar Note è 295 May 16th, 1996 Thermal Finite Elements Analysis of the BaBar Silicon Vertex Tracker D. Barni, D. Giugni, F. Lanni, F. Palombo Dipartimento di Fisica dell'universitça and Sezione INFN,
More informationQuantum Chaos as a Practical Tool in Many-Body Physics
Quantum Chaos as a Practical Tool in Many-Body Physics Vladimir Zelevinsky NSCL/ Michigan State University Supported by NSF Statistical Nuclear Physics SNP2008 Athens, Ohio July 8, 2008 THANKS B. Alex
More informationLecture V : Oscillatory motion and spectral analysis
Lecture V : Oscillatory motion and spectral analysis I. IDEAL PENDULUM AND STABILITY ANALYSIS Let us remind ourselves of the equation of motion for the pendulum. Remembering that the external torque applied
More informationDynamics of a double pendulum with distributed mass
Dynamics of a double pendulum with distributed mass M. Z. Rafat, a M. S. Wheatl, b T. R. Bedding c School of Physics, University of Sydney, NSW 2006, Australia Received 3 April 2008; accepted 1 December
More informationLecture 11 Weak interactions, Cabbibo-angle. angle. SS2011: Introduction to Nuclear and Particle Physics, Part 2
Lecture 11 Weak interactions, Cabbibo-angle angle SS2011: Introduction to Nuclear and Particle Physics, Part 2 1 Neutrino-lepton reactions Consider the reaction of neutrino-electron scattering: Feynman
More information2πi0.j... j. j > 1. R n-1 R 2. R n-2. j > R n. 0> + e 1> n-1. j > 0> + e 1> n-1 n. j > 0> + e 2 n 1>
Technical Report No. 006í507 Coping with Decoherence: Parallelizing the Quantum Fourier Transform Marius Nagy and Selim G. Akl School of Computing Queen's University Kingston, Ontario K7L 3N6 Canada Email:
More informationStanford Linear Accelerator Center. Abstract. After an extremely brief overview of the discovery reach for Z 0 bosons,
SLAC-PUB-8335 January 2000 Z 0 Bosons and Kaluza-Klein Excitations at Muon Colliders 1 Thomas G Rizzo 2 Stanford Linear Accelerator Center Stanford CA 94309, USA Abstract. After an extremely brief overview
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 informationOn long-term dynamics of slow-fast systems with passages through resonances
On long-term dynamics of slow-fast systems with passages through resonances Anatoly Neishtadt Loughborough University, UK Space Research Institute, Moscow, Russia Based on joint works with A.Artemyev,
More informationW 1 æw 2 G + 0 e? u K y Figure 5.1: Control of uncertain system. For MIMO systems, the normbounded uncertainty description is generalized by assuming
Chapter 5 Robust stability and the H1 norm An important application of the H1 control problem arises when studying robustness against model uncertainties. It turns out that the condition that a control
More informationNovember 9, Abstract. properties change in response to the application of a magnetic æeld. These æuids typically
On the Eæective Magnetic Properties of Magnetorheological Fluids æ November 9, 998 Tammy M. Simon, F. Reitich, M.R. Jolly, K. Ito, H.T. Banks Abstract Magnetorheological èmrè æuids represent a class of
More informationMeasurement of magnetic eld line stochasticity in nonlinearly evolving, Department of Engineering Physics,
Measurement of magnetic eld line stochasticity in nonlinearly evolving, nonequilibrium plasmas Y. Nishimura, J. D. Callen, and C. C. Hegna Department of Engineering Physics, University of Wisconsin-Madison,
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 informationarxiv: v3 [nlin.cd] 19 Nov 2013
Why Instantaneous Values of the Covariant Lyapunov Exponents Depend upon the Chosen State-Space Scale Wm. G. Hoover and Carol G. Hoover Ruby Valley Research Institute arxiv:1309.2342v3 [nlin.cd] 19 Nov
More informationSiberian Branch of Russian Academy of Science BUDKER INSTITUTE OF NUCLEAR PHYSICS. B.V. Chirikov POINCAR'E RECURRENCES IN MICROTRON
Siberian Branch of Russian Academy of Science BUDKER INSTITUTE OF NUCLEAR PHYSICS BV Chirikov POINCAR'E RECURRENCES IN MICROTRON AND THE GLOBAL CRITICAL STRUCTURE Budker INP 99-7 NOVOSIBIRSK 1999 Poincare
More informationRegular & Chaotic. collective modes in nuclei. Pavel Cejnar. ipnp.troja.mff.cuni.cz
Pavel Cejnar Regular & Chaotic collective modes in nuclei Institute of Particle and Nuclear Physics Faculty of Mathematics and Physics Charles University, Prague, Czech Republic cejnar @ ipnp.troja.mff.cuni.cz
More informationRadial Basis Functions for Process Control Lyle H. Ungar Tom Johnson Richard D. De Veaux University ofpennsylvania Voice Processing Corp. Princeton Un
Radial Basis Functions for Process Control Lyle H. Ungar Tom Johnson Richard D. De Veaux University ofpennsylvania Voice Processing Corp. Princeton University Abstract Radial basis function èrbfsè neural
More informationNANOSCALE SCIENCE & TECHNOLOGY
. NANOSCALE SCIENCE & TECHNOLOGY V Two-Level Quantum Systems (Qubits) Lecture notes 5 5. Qubit description Quantum bit (qubit) is an elementary unit of a quantum computer. Similar to classical computers,
More informationAbstract To characterize color values measured by color devices èe.g. scanners, color copiers and color camerasè in a deviceíindependent fashion these
White-point preserving color correction Graham D. Finlayson Department of Computer Science, University of York, York, United Kingdom YO1 5DD Mark S. Drew School of Computing Science, Simon Fraser University,
More informationAn Exact Local Hybrid Monte Carlo Algorithm. for Gauge Theories. A. D. Kennedy and K. M. Bitar. Supercomputer Computations Research Institute
An Exact Local Hybrid Monte Carlo Algorithm for Gauge Theories A. D. Kennedy and K. M. Bitar Supercomputer Computations Research Institute Florida State University Tallahassee, FL 06{405, USA Internet:
More informationGRAVITATIONAL WAVES AND THE END OF INFLATION. Richard Easther (Yale)
GRAVITATIONAL WAVES AND THE END OF INFLATION Richard Easther (Yale) OUTLINE Inflation: a reminder Ending inflation: Parametric resonance / preheating [SKIP: technical calculation] Gravitational wave generation
More informationarxiv: v1 [nlin.cd] 25 May 2017
Synchronization, phase slips and coherent structures in area-preserving maps Swetamber Das, Sasibhusan Mahata, and Neelima Gupte Department of Physics, Indian Institute of Technology Madras, Chennai, 636,
More informationLECTURE 17. Algorithms for Polynomial Interpolation
LECTURE 7 Algorithms for Polynomial Interpolation We have thus far three algorithms for determining the polynomial interpolation of a given set of data.. Brute Force Method. Set and solve the following
More informationOscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. small angle approximation. Oscillatory solution
Oscillatory Motion Simple pendulum: linear Hooke s Law restoring force for small angular deviations d 2 θ dt 2 = g l θ small angle approximation θ l Oscillatory solution θ(t) =θ 0 sin(ωt + φ) F with characteristic
More informationTHE DYNAMICS OF CHEMICAL REACTORS WITH HEAT INTEGRATION. John Morud Sigurd Skogestad æ. Chemical Engineering, University oftrondheim-nth
THE DYNAMICS OF CHEMICAL REACTORS WITH HEAT INTEGRATION John Morud Sigurd Skogestad æ Chemical Engineering, University oftrondheim-nth N-7034 Trondheim, Norway Keywords: Chemical reactor, linear analysis,
More informationAbstract. subplans to arrive at a more eæcient implementation. Our algorithm improves upon the
An OèT 3 è algorithm for the economic lot-sizing problem with constant capacities C.P.M. van Hoesel æ A.P.M. Wagelmans y Revised June 1994 Abstract We develop an algorithm that solves the constant capacities
More informationOscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. Oscillatory solution
Oscillatory Motion Simple pendulum: linear Hooke s Law restoring force for small angular deviations d 2 θ dt 2 = g l θ θ l Oscillatory solution θ(t) =θ 0 sin(ωt + φ) F with characteristic angular frequency
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 informationBIFURCATIONS AND CHAOS IN A PERIODICALLY FORCED PROTOTYPE ADAPTIVE CONTROL SYSTEM 1
KYBERNETIKA VOLUME 3 0 (1994), NUMBER 2, PAG ES 121-128 BIFURCATIONS AND CHAOS IN A PERIODICALLY FORCED PROTOTYPE ADAPTIVE CONTROL SYSTEM 1 YURI A. KUZNETSOV AND CARLO PICCARDI An adaptive control system
More informationTheory of Adiabatic Invariants A SOCRATES Lecture Course at the Physics Department, University of Marburg, Germany, February 2004
Preprint CAMTP/03-8 August 2003 Theory of Adiabatic Invariants A SOCRATES Lecture Course at the Physics Department, University of Marburg, Germany, February 2004 Marko Robnik CAMTP - Center for Applied
More informationLagrangian Coherent Structures (LCS)
Lagrangian Coherent Structures (LCS) CDS 140b - Spring 2012 May 15, 2012 ofarrell@cds.caltech.edu A time-dependent dynamical system ẋ (t; t 0, x 0 )=v(x(t;,t 0, x 0 ),t) x(t 0 ; t 0, x 0 )=x 0 t 2 I R
More informationTwo models for the parametric forcing of a nonlinear oscillator
Nonlinear Dyn (007) 50:147 160 DOI 10.1007/s11071-006-9148-3 ORIGINAL ARTICLE Two models for the parametric forcing of a nonlinear oscillator Nazha Abouhazim Mohamed Belhaq Richard H. Rand Received: 3
More informationNonlinear Single-Particle Dynamics in High Energy Accelerators
Nonlinear Single-Particle Dynamics in High Energy Accelerators Part 4: Canonical Perturbation Theory Nonlinear Single-Particle Dynamics in High Energy Accelerators There are six lectures in this course
More informationDynamics of a double pendulum with distributed mass. Abstract
The following article has been accepted by the American Journal of Physics. After it is published, it will be found at http://scitation.aip.org/ajp. Dynamics of a double pendulum with distributed mass
More informationTheory of selective excitation in stimulated Raman scattering
Theory of selective excitation in stimulated Raman scattering S. A. Malinovskaya, P. H. Bucksbaum, and P. R. Berman Michigan Center for Theoretical Physics, FOCUS Center, and Department of Physics, University
More informationIn the previous chapters we have presented synthesis methods for optimal H 2 and
Chapter 8 Robust performance problems In the previous chapters we have presented synthesis methods for optimal H 2 and H1 control problems, and studied the robust stabilization problem with respect to
More informationQuantum chaos on graphs
Baylor University Graduate seminar 6th November 07 Outline 1 What is quantum? 2 Everything you always wanted to know about quantum but were afraid to ask. 3 The trace formula. 4 The of Bohigas, Giannoni
More informationTransitioning to Chaos in a Simple Mechanical Oscillator
Transitioning to Chaos in a Simple Mechanical Oscillator Hwan Bae Physics Department, The College of Wooster, Wooster, Ohio 69, USA (Dated: May 9, 8) We vary the magnetic damping, driver frequency, and
More informationPhysics 115/242 Comparison of methods for integrating the simple harmonic oscillator.
Physics 115/4 Comparison of methods for integrating the simple harmonic oscillator. Peter Young I. THE SIMPLE HARMONIC OSCILLATOR The energy (sometimes called the Hamiltonian ) of the simple harmonic oscillator
More informationFORM 1 Physics 240 í Fall Exam 1 CONTI COULTER GRAFE SKALSEY ZORN. Please æll in the information requested on the scantron answer sheet
FORM 1 Physics 240 í Fall 1999 Exam 1 Name: Discussion section instructor and time: ècircle ONE TIMEè CONTI COULTER GRAFE SKALSEY ZORN MW 10-11 èè23è MW 8-9 èè4è MW 2-3 èè7è MW 11-12 èè6è MW 9-10 èè12è
More information