Analytical approaches for the. dynamics of charged particle beams
|
|
- Mavis Norton
- 6 years ago
- Views:
Transcription
1 Analytical approaches for the dynamics of charged particle beams J. Struckmeier, GSI Frankfurt am Main, 21 December
2 Outline Review: Hamiltonian mechanics of N particle systems Liouville s theorem Continuous description of beam dynamics Vlasov equation Moment analysis of the Vlasov equation Beam envelope equations Comparison of analytic, simulated, and measured quantities Remarks on irreversibility in computer simulations Outlook: non-liouvillean effects 2
3 The Three Pillars of today s Physics Analytical Description ("Theory") Computer Simulation Measurements 3
4 Hamilton mechanics: basis of our description of dynamical systems. Hamiltonian H for a system of N (interacting) particles H = m 2 N vi 2 + V (x 1,..., x N, t), i=1 i = 1,..., N With H given, we obtain from Hamilton s variational principle the canonical equations for each i = 1,..., N ẋ i = H v i = v i, m v i = H x i = V (x 1,..., x N, t) x i Together with the initial condition ( x1 (t 0 ), v 1 (t 0 ),..., x N (t 0 ), v N (t 0 ) ), the system of N particles is completely determined by the generally coupled system of 6N first order differential equations. 4
5 Properties of this description deterministic: The system is fully determined, no sources exist that could cause loss of information on the system s state reversible: the equations of motion ( canonical equations ) are invariant with respect to time-inversion transformations t t x i x i, v i v i. If we know the initial condition (at t 0 ) and the Hamiltonian, we can integrate the canonical equations to any instant of time, either t > t 0, or t < t 0 The canonical equations form the basis for computer simulations. In most cases, we have N sim N (representative sample). Analytically hopeless : even the three-body-problem lacks an analytical solution. For analytical studies, simplifying assumptions are necessary. 5
6 Liouville s theorem Liouville s theorem (proper) states that the volume element dγ dγ = dx 1 dv 1... dx N dv N = dx 1dv 1... dx Ndv N of a Hamiltonian N particle system is invariant with respect to canonical transformations in particular with respect to the system s time evolution. just a one invariant out of the set of Poincaré invariants. We may easily proof Liouville s theorem by showing that the determinant of the Jacobi matrix associated with a canonical transformation is unity. Note: the projection of dγ to subspaces especially to the 6-dimensional phase space or 2-dimensional phase-space planes does not necessarily provide a conserved quantity! We must be very careful if we draw conclusions that apply to subspaces! 6
7 We summarize: The single particle approach with N sim N is appropriate for computer simulations. It does not work for analytical studies. Alternative: description in terms of statistical mechanics continuous description of beam dynamics. 7
8 Continuous description of beam dynamics For the continuous description of the ion beam, we define a 6-dimensional phase-space probability density function f = f(x, v, t) f dx dv provides the probability of finding a particle inside the volume dx dv around the phase-space point (x, v) at time t. f is a smooth function of the phase-space variables. f does not provide information on individual particles. The special Liouville theorem df dt = 0 applies for systems with non-negligible Coulomb interaction only if the system can be described to sufficient accuracy by a continuous Hamilton function. This means that we must be allowed to treat the internal space charge forces analogously to an external force field. 8
9 We observe: The statistical description smoothes out the effects reflecting the actual charge granularity. The usual equation of motion of the phase-space probability density df/dt = 0 does not cover intra-beam scattering effects, for example. df/dt = 0 does not hold for stochastic cooling! 9
10 Remark on stochastic cooling Strahl Pick-up Verstaerker Kicker 10
11 The cooling process hinges on charge density fluctuations within the beam that are sampled by a pick-up electrode. This signal is used to counteract the related momentum fluctuations within the beam at some specific point further downstream the ring. The continuous description does not apply! The special Liouville theorem for f does not hold, hence stochastic cooling does not conflict with Liouville s theorem (in its special form). 11
12 Vlasov equation Explicitly, the Liouville equation for f(x, v, t) writes df dt = f 3 t + i=1 v i f x i + 3 i=1 ẍ i f v i = 0, Inserting ẍ, this equation is usually referred to as the Vlasov equation f t + v xf + 1 m (F ext + qe ssc ) v f = 0. Herein F ext stands for the external forces applied to the beam by focusing elements (such as quadrupoles) and E ssc for the smooth part of the space-charge forces. E ssc depends on f via Poisson s equation: div E ssc (x, t) = q f(x, v, t) dv. ɛ 0 12
13 We summarize: The statistical description replaces the original problem of solving 6N coupled first-order equations by one equation of motion for the 6-dimensional probability density f(x, v, t). The coupled set of Vlasov and Poisson equations provides the closed equation of motion for f if all effects due to the actual charge granularity can be neglected. Even this equation is generally not worth-while to solve directly. We do not really want to know f(x, v, t) in detail. We must switch to even more global quantities in our analytical description of ion beams. 13
14 Moment analysis of the Vlasov equation The usual way to switch to even more global physical quantities is to consider moments of f(x, v, t) x 2 i (t) = x 2 i f(x, v, t) dxdv, i = 1, 2, 3. x 2 i is proportional to the actual beam width in x i. Idea: Instead of solving the equation of motion for f, we are going to set up and solve the equations of motion for the second beam moments. The derivatives of the moments are calculated according to d x 2 dt i = x 2 f i t dxdv, inserting f/ t from the Vlasov-Poisson equation. 14
15 Integrating by parts, we obtain for each phase-space plane i = 1, 2, 3 a coupled set of moment equations d dt x 2 i 2 xi v i = 0 m d dt xi v i m v 2 i xi F ext,i q xi E ssc,i = 0 m d dt v 2 i 2 vi F ext,i 2q vi E ssc,i = 0 As usual, we define the rms emittance ε i (t) as ε 2 i (t) = x 2 i v 2 2 i xi v i The time derivative of the rms emittance may be arranged as d dt ε2 i (t) = d dt ε2 i (t) ext + d dt ε2 i (t) ssc 15
16 RMS beam envelope equations For unbunched beams with elliptic cross section in real space, we have x i E sc,i = I 4πɛ 0 cβ x 2 i x 2 + y 2. We obtain the rms envelope equation from the first two moment equations d 2 dt 2 x2 + ω 2 x(t) x 2 d 2 dt 2 y2 + ω 2 y(t) y 2 qi 4πɛ 0 mcβ qi 4πɛ 0 mcβ 1 x2 + y 2 1 x2 + y 2 ε2 x(t) x 2 3 = 0 ε2 y(t) y 2 3 = 0. We must assume the rms emittances to be approximately constant to be able to integrate this coupled set. valid for linear external forces and if no transient effects occur. 16
17 2.5 GSI QUADRUPOLE CHANNEL, SIGMA-0=120 DEG., SIGMA=35 DEG., MATCHED kv x (cm) y (cm) Cells
18 Comparison of analytic, simulation, and measured quantities In our analytic description, the second moment in x was defined as x 2 = x 2 f(x, v, t) dxdv. x2 is proportional to the actual beam width in x. In a computer simulation based on N sim representative particles, this quantity is given by x 2 = 1 N sim How do we measure this quantity? N sim i=1 x 2 i. 18
19 19
20 The result of such a slit and detector measurement is given by the current matrix I i i 1m. I =......, i n1... i nm with m the number of slit positions, and n the number of detectors. For a fixed step size x of the slit positions, the second beam moments can now be directly calculated x 2 = I 1 t m n ) 2 i kj (j x x, j=1 k=1 with I t = m n i kj, x = m n i kj j x. j=1 k=1 j=1 k=1 20
21 N Teilchensystem Bewegungsgleichungen (kan. Gleichungen) N Teilchensystem zum Zeitpunkt t0 zum Zeitpunkt t1 Projektion auf abgeleitete Größe abgeleitete Größe zum Zeitpunkt t0 Bewegungsgleichung? abgeleitete Größe zum Zeitpunkt t1 21
22 We summarize: We must clearly distinguish quantities that have only meaning at a fixed instant of time ( figures of merit, like 90 % emittances ) from quantities that are related to the system s equations of motion. If we want to understand the dynamics of measured quantities, we must ensure that these quantities are physical in the sense that they must be consistent with the equations of motion. an equation of motion must exist for the measured quantities! The rms description provides a basis to compare measured emittances, beam profiles, etc. with results of computer simulations and analytical approximations. The rms description avoids all problems that are associated with the definitions of boundaries. Best possible accuracy! 22
23 Remarks on irreversibility Consider the transformation that reverses the direction of time flow: t t x i x i, v i v i, E E, B B. We may investigate the single particle equations of motion with respect to this transformation mẍ i F ext (x i, t) q ( E i + v i B i ) = 0, i = 1,..., Nsim. We observe that this equation is indeed invariant under time reversal. We may convince ourselves that the Vlasov equation is also invariant. Earlier states are fully restored just like a movie that is reversed. The equations describe the reversible aspects of the system evolution. But: computer simulations of dynamical systems are irreversible even if the coded equations of motion are strictly reversible and the integration algorithm maintains the symplectic nature of the canonical equations ( symplectic integrator technique ). 23
24 ε(t/τ) / ε(0) Cells (t/τ) 24
25 ε(t/τ) / ε(0) Cells (t/τ) 25
26 ε(t/τ) / ε(0) Cells (t/τ) 26
27 27
28 Non Liouvillean Effects An additional equation, hence a generalization of the special Liouville theorem is necessary (Chandrasekhar (1943)): [ ] df f dt = 0 +. t If the non-liouvillean effects are small compared to the macroscopic forces (smooth space charge and external forces), we can describe them NL by the Fokker-Planck equation: [ ] f = {F i (v, t) f} + t NL v i i i,j 2 v i v j {D ij (v, t) f} D ij are referred to as elements of the diffusion tensor and the F i as elements of the drift vector that describes the dynamical friction forces. They must be determined appropriately depending on the nature of the underlying physical process. 28
29 We summarize: The actual time evolution of the simulated system always comprises irreversible aspects even if the actually coded equations of motion are strictly reversible. A computer simulation based on individual particles can never be an exact realization of a solution of the Vlasov equation. We must include these effects in our analytical description in order to fully understand our simulation results. 29
Calculation of matched beams under space-charge conditions
Calculation of matched beams under space-charge conditions Jürgen Struckmeier j.struckmeier@gsi.de www.gsi.de/~struck Vortrag im Rahmen des Winterseminars Aktuelle Probleme der Beschleuniger- und Plasmaphysik
More information4. Statistical description of particle beams
4. Statistical description of particle beams 4.1. Beam moments 4. Statistical description of particle beams 4.1. Beam moments In charged particle beam dynamics, we are commonly not particularly interested
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 informationReview: control, feedback, etc. Today s topic: state-space models of systems; linearization
Plan of the Lecture Review: control, feedback, etc Today s topic: state-space models of systems; linearization Goal: a general framework that encompasses all examples of interest Once we have mastered
More informationLangevin Methods. Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 10 D Mainz Germany
Langevin Methods Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 1 D 55128 Mainz Germany Motivation Original idea: Fast and slow degrees of freedom Example: Brownian motion Replace
More informationPhase Space Gymnastics
Phase Space Gymnastics As accelerator technology advances, the requirements on accelerator beam quality become increasingly demanding. Phase space gymnastics becomes a new focus of accelerator physics
More informationNonlinear Single-Particle Dynamics in High Energy Accelerators
Nonlinear Single-Particle Dynamics in High Energy Accelerators Part 2: Basic tools and concepts Nonlinear Single-Particle Dynamics in High Energy Accelerators This course consists of eight lectures: 1.
More informationHypocoercivity for kinetic equations with linear relaxation terms
Hypocoercivity for kinetic equations with linear relaxation terms Jean Dolbeault dolbeaul@ceremade.dauphine.fr CEREMADE CNRS & Université Paris-Dauphine http://www.ceremade.dauphine.fr/ dolbeaul (A JOINT
More informationLecture 5. Alexey Boyarsky. October 21, Legendre transformation and the Hamilton equations of motion
Lecture 5 Alexey Boyarsky October 1, 015 1 The Hamilton equations of motion 1.1 Legendre transformation and the Hamilton equations of motion First-order equations of motion. In the Lagrangian formulation,
More informationUne approche hypocoercive L 2 pour l équation de Vlasov-Fokker-Planck
Une approche hypocoercive L 2 pour l équation de Vlasov-Fokker-Planck Jean Dolbeault dolbeaul@ceremade.dauphine.fr CEREMADE CNRS & Université Paris-Dauphine http://www.ceremade.dauphine.fr/ dolbeaul (EN
More informationCaltech Ph106 Fall 2001
Caltech h106 Fall 2001 ath for physicists: differential forms Disclaimer: this is a first draft, so a few signs might be off. 1 Basic properties Differential forms come up in various parts of theoretical
More informationMATHEMATICAL PHYSICS
MATHEMATICAL PHYSICS Third Year SEMESTER 1 015 016 Classical Mechanics MP350 Prof. S. J. Hands, Prof. D. M. Heffernan, Dr. J.-I. Skullerud and Dr. M. Fremling Time allowed: 1 1 hours Answer two questions
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 informationSelf-consistent analysis of three dimensional uniformly charged ellipsoid with zero emittance
1 SLAC-PUB-883 17 May 1 Self-consistent analysis of three dimensional uniformly charged ellipsoid with zero emittance Yuri K. Batygin Stanford Linear Accelerator Center, Stanford University, Stanford,
More informationThe dynamics of small particles whose size is roughly 1 µmt or. smaller, in a fluid at room temperature, is extremely erratic, and is
1 I. BROWNIAN MOTION The dynamics of small particles whose size is roughly 1 µmt or smaller, in a fluid at room temperature, is extremely erratic, and is called Brownian motion. The velocity of such particles
More informationBrownian Motion: Fokker-Planck Equation
Chapter 7 Brownian Motion: Fokker-Planck Equation The Fokker-Planck equation is the equation governing the time evolution of the probability density of the Brownian particla. It is a second order differential
More informationGravitational-like interactions in a cloud of cold atoms?
Gravitational-like interactions in a cloud of cold atoms? J. Barré (1), B. Marcos (1), A. Olivetti (1), D. Wilkowski (2), M. Chalony (2) (+ R. Kaiser (2) ) 1 Laboratoire JA Dieudonné, U. de Nice-Sophia
More informationLectures on basic plasma physics: Kinetic approach
Lectures on basic plasma physics: Kinetic approach Department of applied physics, Aalto University April 30, 2014 Motivation Layout 1 Motivation 2 Boltzmann equation (a nasty bastard) 3 Vlasov equation
More informationG : Statistical Mechanics
G25.2651: Statistical Mechanics Notes for Lecture 1 Defining statistical mechanics: Statistical Mechanics provies the connection between microscopic motion of individual atoms of matter and macroscopically
More informationPHYS 705: Classical Mechanics. Hamiltonian Formulation & Canonical Transformation
1 PHYS 705: Classical Mechanics Hamiltonian Formulation & Canonical Transformation Legendre Transform Let consider the simple case with ust a real value function: F x F x expresses a relationship between
More information1 Phase Spaces and the Liouville Equation
Phase Spaces and the Liouville Equation emphasize the change of language from deterministic to probablistic description. Under the dynamics: ½ m vi = F i ẋ i = v i with initial data given. What is the
More informationHamilton-Jacobi theory
Hamilton-Jacobi theory November 9, 04 We conclude with the crowning theorem of Hamiltonian dynamics: a proof that for any Hamiltonian dynamical system there exists a canonical transformation to a set of
More informationHamiltonian and Non-Hamiltonian Reductions of Charged Particle Dynamics: Diffusion and Self-Organization
NNP2017 11 th July 2017 Lawrence University Hamiltonian and Non-Hamiltonian Reductions of Charged Particle Dynamics: Diffusion and Self-Organization N. Sato and Z. Yoshida Graduate School of Frontier Sciences
More information(a) Sections 3.7 through (b) Sections 8.1 through 8.3. and Sections 8.5 through 8.6. Problem Set 4 due Monday, 10/14/02
Physics 601 Dr. Dragt Fall 2002 Reading Assignment #4: 1. Dragt (a) Sections 1.5 and 1.6 of Chapter 1 (Introductory Concepts). (b) Notes VI, Hamilton s Equations of Motion (to be found right after the
More informationIndex. Accelerator model 8 Adiabatic damping 32, 141 Air-bag model 338 Alternating explicit time scheme 112 Azimuthal modes, see Modes
Index Accelerator model 8 Adiabatic damping 32, 141 Air-bag model 338 Alternating explicit time scheme 112 Azimuthal modes, see Modes Beam breakup in linacs dipole mode 136 higher modes 160 quadrupole
More informationHydrodynamics, Thermodynamics, and Mathematics
Hydrodynamics, Thermodynamics, and Mathematics Hans Christian Öttinger Department of Mat..., ETH Zürich, Switzerland Thermodynamic admissibility and mathematical well-posedness 1. structure of equations
More informationSolutions of M3-4A16 Assessed Problems # 3 [#1] Exercises in exterior calculus operations
D. D. Holm Solutions to M3-4A16 Assessed Problems # 3 15 Dec 2010 1 Solutions of M3-4A16 Assessed Problems # 3 [#1] Exercises in exterior calculus operations Vector notation for differential basis elements:
More informationGas Dynamics: Basic Equations, Waves and Shocks
Astrophysical Dynamics, VT 010 Gas Dynamics: Basic Equations, Waves and Shocks Susanne Höfner Susanne.Hoefner@fysast.uu.se Astrophysical Dynamics, VT 010 Gas Dynamics: Basic Equations, Waves and Shocks
More informationClassical Statistical Mechanics: Part 1
Classical Statistical Mechanics: Part 1 January 16, 2013 Classical Mechanics 1-Dimensional system with 1 particle of mass m Newton s equations of motion for position x(t) and momentum p(t): ẋ(t) dx p =
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 informationAccelerator Physics. Tip World Scientific NEW JERSEY LONDON SINGAPORE BEIJING SHANGHAI HONG KONG TAIPEI BANGALORE. Second Edition. S. Y.
Accelerator Physics Second Edition S. Y. Lee Department of Physics, Indiana University Tip World Scientific NEW JERSEY LONDON SINGAPORE BEIJING SHANGHAI HONG KONG TAIPEI BANGALORE Contents Preface Preface
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 informationClassical Mechanics in Hamiltonian Form
Classical Mechanics in Hamiltonian Form We consider a point particle of mass m, position x moving in a potential V (x). It moves according to Newton s law, mẍ + V (x) = 0 (1) This is the usual and simplest
More informationPhysics 208, Spring 2016 Exam #2
Physics 208, Spring 2016 Exam #2 A Name (Last, First): ID #: Section #: You have 75 minutes to complete the exam. Formulae are provided on an attached sheet. You may NOT use any other formula sheet. You
More informationThe kinetic equation (Lecture 11)
The kinetic equation (Lecture 11) January 29, 2016 190/441 Lecture outline In the preceding lectures we focused our attention on a single particle motion. In this lecture, we will introduce formalism for
More informationSynchrotron Motion with Space-Charge
Synchrotron Motion with Space-Charge Basic Principles without space-charge RF resonant cavity providing accelerating voltage V (t). Often V = V 0 sin(φ s + ω rf t), where ω rf is the angular frequency
More informationMagnets and Lattices. - Accelerator building blocks - Transverse beam dynamics - coordinate system
Magnets and Lattices - Accelerator building blocks - Transverse beam dynamics - coordinate system Both electric field and magnetic field can be used to guide the particles path. r F = q( r E + r V r B
More informationFigure-8 Storage RingF8SRNon Neutral Plasma Confinement in Curvilinear Guiding Fields
Figure-8 Storage Ring F8SR Non Neutral Plasma Confinement in Curvilinear Guiding Fields Joschka F. Wagner Institute of Applied Physics (IAP) Workgroup Prof. Ulrich Ratzinger Non Neutral Plasma-Group (NNP)
More informationNew Physical Principle for Monte-Carlo simulations
EJTP 6, No. 21 (2009) 9 20 Electronic Journal of Theoretical Physics New Physical Principle for Monte-Carlo simulations Michail Zak Jet Propulsion Laboratory California Institute of Technology, Advance
More informationStationary mean-field games Diogo A. Gomes
Stationary mean-field games Diogo A. Gomes We consider is the periodic stationary MFG, { ɛ u + Du 2 2 + V (x) = g(m) + H ɛ m div(mdu) = 0, (1) where the unknowns are u : T d R, m : T d R, with m 0 and
More informationMetric spaces and metrizability
1 Motivation Metric spaces and metrizability By this point in the course, this section should not need much in the way of motivation. From the very beginning, we have talked about R n usual and how relatively
More informationThe Particle-Field Hamiltonian
The Particle-Field Hamiltonian For a fundamental understanding of the interaction of a particle with the electromagnetic field we need to know the total energy of the system consisting of particle and
More informationELEG 3143 Probability & Stochastic Process Ch. 6 Stochastic Process
Department of Electrical Engineering University of Arkansas ELEG 3143 Probability & Stochastic Process Ch. 6 Stochastic Process Dr. Jingxian Wu wuj@uark.edu OUTLINE 2 Definition of stochastic process (random
More informationanalysis for transport equations and applications
Multi-scale analysis for transport equations and applications Mihaï BOSTAN, Aurélie FINOT University of Aix-Marseille, FRANCE mihai.bostan@univ-amu.fr Numerical methods for kinetic equations Strasbourg
More informationFeynman s path integral approach to quantum physics and its relativistic generalization
Feynman s path integral approach to quantum physics and its relativistic generalization Jürgen Struckmeier j.struckmeier@gsi.de, www.gsi.de/ struck Vortrag im Rahmen des Winterseminars Aktuelle Probleme
More informationA Primer on Three Vectors
Michael Dine Department of Physics University of California, Santa Cruz September 2010 What makes E&M hard, more than anything else, is the problem that the electric and magnetic fields are vectors, and
More informationFluid-Particles Interaction Models Asymptotics, Theory and Numerics I
Fluid-Particles Interaction Models Asymptotics, Theory and Numerics I J. A. Carrillo collaborators: T. Goudon (Lille), P. Lafitte (Lille) and F. Vecil (UAB) (CPDE 2005), (JCP, 2008), (JSC, 2008) ICREA
More informationUtility of Correlation Functions
Utility of Correlation Functions 1. As a means for estimating power spectra (e.g. a correlator + WK theorem). 2. For establishing characteristic time scales in time series (width of the ACF or ACV). 3.
More informationChap. 1. Some Differential Geometric Tools
Chap. 1. Some Differential Geometric Tools 1. Manifold, Diffeomorphism 1.1. The Implicit Function Theorem ϕ : U R n R n p (0 p < n), of class C k (k 1) x 0 U such that ϕ(x 0 ) = 0 rank Dϕ(x) = n p x U
More informationViscosity Solutions for Dummies (including Economists)
Viscosity Solutions for Dummies (including Economists) Online Appendix to Income and Wealth Distribution in Macroeconomics: A Continuous-Time Approach written by Benjamin Moll August 13, 2017 1 Viscosity
More informationExperimental Storage Ring - ESR E max = 420 MeV/u, 10 Tm, electron-, stochastic- and laser cooling. Indian Institute of Technology Ropar
Experimental Storage Ring - ESR E max = 420 MeV/u, 10 Tm, electron-, stochastic- and laser cooling Specification of the ESR Particle detectors Re-injection to SIS Two 5 kv rf-cavities Fast Injection Schottky
More information1. Electricity and Magnetism (Fall 1995, Part 1) A metal sphere has a radius R and a charge Q.
1. Electricity and Magnetism (Fall 1995, Part 1) A metal sphere has a radius R and a charge Q. (a) Compute the electric part of the Maxwell stress tensor T ij (r) = 1 {E i E j 12 } 4π E2 δ ij both inside
More informationPhysics 351 Wednesday, February 14, 2018
Physics 351 Wednesday, February 14, 2018 HW4 due Friday. For HW help, Bill is in DRL 3N6 Wed 4 7pm. Grace is in DRL 2C2 Thu 5:30 8:30pm. Respond at pollev.com/phys351 or text PHYS351 to 37607 once to join,
More informationMechanics Physics 151
Mechanics Physics 151 Lecture 4 Continuous Systems and Fields (Chapter 13) What We Did Last Time Built Lagrangian formalism for continuous system Lagrangian L Lagrange s equation = L dxdydz Derived simple
More informationPath integrals for classical Markov processes
Path integrals for classical Markov processes Hugo Touchette National Institute for Theoretical Physics (NITheP) Stellenbosch, South Africa Chris Engelbrecht Summer School on Non-Linear Phenomena in Field
More informationGEOMETRIC QUANTIZATION
GEOMETRIC QUANTIZATION 1. The basic idea The setting of the Hamiltonian version of classical (Newtonian) mechanics is the phase space (position and momentum), which is a symplectic manifold. The typical
More informationEva Miranda. UPC-Barcelona and BGSMath. XXV International Fall Workshop on Geometry and Physics Madrid
b-symplectic manifolds: going to infinity and coming back Eva Miranda UPC-Barcelona and BGSMath XXV International Fall Workshop on Geometry and Physics Madrid Eva Miranda (UPC) b-symplectic manifolds Semptember,
More informationfor changing independent variables. Most simply for a function f(x) the Legendre transformation f(x) B(s) takes the form B(s) = xs f(x) with s = df
Physics 106a, Caltech 1 November, 2018 Lecture 10: Hamiltonian Mechanics I The Hamiltonian In the Hamiltonian formulation of dynamics each second order ODE given by the Euler- Lagrange equation in terms
More informationTemperature and Pressure Controls
Ensembles Temperature and Pressure Controls 1. (E, V, N) microcanonical (constant energy) 2. (T, V, N) canonical, constant volume 3. (T, P N) constant pressure 4. (T, V, µ) grand canonical #2, 3 or 4 are
More informationELECTRON DYNAMICS WITH SYNCHROTRON RADIATION
ELECTRON DYNAMICS WITH SYNCHROTRON RADIATION Lenny Rivkin Ecole Polythechnique Federale de Lausanne (EPFL) and Paul Scherrer Institute (PSI), Switzerland CERN Accelerator School: Introduction to Accelerator
More informationStochastic equations for thermodynamics
J. Chem. Soc., Faraday Trans. 93 (1997) 1751-1753 [arxiv 1503.09171] Stochastic equations for thermodynamics Roumen Tsekov Department of Physical Chemistry, University of Sofia, 1164 Sofia, ulgaria The
More informationslhc-project-report-0029
EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH European Laboratory for Particle Physics slhc Project slhc Project Report 29 Transverse Emittance Measurement with the Three-Monitor-Method at the Linac4 K. Hanke,
More informationComputing the Universe, Dynamics IV: Liouville, Boltzmann, Jeans [DRAFT]
Computing the Universe, Dynamics IV: Liouville, Boltzmann, Jeans [DRAFT] Martin Weinberg June 24, 1998 1 Liouville s Theorem Liouville s theorem is a fairly straightforward consequence of Hamiltonian dynamics.
More informationKinetic theory of gases
Kinetic theory of gases Toan T. Nguyen Penn State University http://toannguyen.org http://blog.toannguyen.org Graduate Student seminar, PSU Jan 19th, 2017 Fall 2017, I teach a graduate topics course: same
More informationPoisson Brackets and Lie Operators
Poisson Brackets and Lie Operators T. Satogata January 22, 2008 1 Symplecticity and Poisson Brackets 1.1 Symplecticity Consider an n-dimensional 2n-dimensional phase space) linear system. Let the canonical
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 informationGaussian processes for inference in stochastic differential equations
Gaussian processes for inference in stochastic differential equations Manfred Opper, AI group, TU Berlin November 6, 2017 Manfred Opper, AI group, TU Berlin (TU Berlin) inference in SDE November 6, 2017
More informationThe Geometry of Euler s equation. Introduction
The Geometry of Euler s equation Introduction Part 1 Mechanical systems with constraints, symmetries flexible joint fixed length In principle can be dealt with by applying F=ma, but this can become complicated
More information4 Classical Coherence Theory
This chapter is based largely on Wolf, Introduction to the theory of coherence and polarization of light [? ]. Until now, we have not been concerned with the nature of the light field itself. Instead,
More informationHamilton s principle and Symmetries
Hamilton s principle and Symmetries Sourendu Gupta TIFR, Mumbai, India Classical Mechanics 2011 August 18, 2011 The Hamiltonian The change in the Lagrangian due to a virtual change of coordinates is dl
More informationarxiv:physics/ v1 [physics.class-ph] 3 Apr 2000
Dirac monopole with Feynman brackets Alain Bérard arxiv:physics/0004008v1 [physicsclass-ph] 3 Apr 2000 LPLI-Institut de Physique, 1 blvd DArago, F-57070 Metz, France Y Grandati LPLI-Institut de Physique,
More informationMathematical modelling of collective behavior
Mathematical modelling of collective behavior Young-Pil Choi Fakultät für Mathematik Technische Universität München This talk is based on joint works with José A. Carrillo, Maxime Hauray, and Samir Salem
More informationMathematical Structures of Statistical Mechanics: from equilibrium to nonequilibrium and beyond Hao Ge
Mathematical Structures of Statistical Mechanics: from equilibrium to nonequilibrium and beyond Hao Ge Beijing International Center for Mathematical Research and Biodynamic Optical Imaging Center Peking
More informationIntroduction to electron and photon beam physics. Zhirong Huang SLAC and Stanford University
Introduction to electron and photon beam physics Zhirong Huang SLAC and Stanford University August 03, 2015 Lecture Plan Electron beams (1.5 hrs) Photon or radiation beams (1 hr) References: 1. J. D. Jackson,
More informationFluid Mechanics for Astrophysicists
Fluid Mechanics for Astrophysicists MAUCA Fundamental Course 2018 Course Notes Instructor: Oliver HAHN oliver.hahn@oca.eu Contents 1 From Microscopic Particles to the Fluid Equations 3 1.1 Describing many
More informationExtensions of the Longitudinal Envelope Equation
Extensions of the Longitudinal Envelope Equation David Neuffer A ccelerator Physics Group Fermilab, PO Box 500, Batavia IL 6050 Introduction Recently, longitudinal space charge effects have become of increased
More informationSpace Charge in Linear Machines
Space Charge in Linear Machines Massimo.Ferrario@LNF.INFN.IT Egham September 6 th 017 Relativistic equation of motion dp dt = F p = γm o v γm o dv dt + m ov dγ dt = F β = v c dγ dt = d dt " a v % m o γ
More informationVariational principles and Hamiltonian Mechanics
A Primer on Geometric Mechanics Variational principles and Hamiltonian Mechanics Alex L. Castro, PUC Rio de Janeiro Henry O. Jacobs, CMS, Caltech Christian Lessig, CMS, Caltech Alex L. Castro (PUC-Rio)
More informationThe First Principle Calculation of Green Kubo Formula with the Two-Time Ensemble Technique
Commun. Theor. Phys. (Beijing, China 35 (2 pp. 42 46 c International Academic Publishers Vol. 35, No. 4, April 5, 2 The First Principle Calculation of Green Kubo Formula with the Two-Time Ensemble Technique
More information7. Kinetics controlled by fluctuations: Kramers theory of activated processes
7. Kinetics controlled by fluctuations: Kramers theory of activated processes Macroscopic kinetic processes (time dependent concentrations) Elementary kinetic process Reaction mechanism Unimolecular processes
More informationFYS 3120: Classical Mechanics and Electrodynamics
FYS 3120: Classical Mechanics and Electrodynamics Formula Collection Spring semester 2014 1 Analytical Mechanics The Lagrangian L = L(q, q, t), (1) is a function of the generalized coordinates q = {q i
More informationPIC solvers for Intense Beams in Synchrotrons
PIC solvers for Intense Beams in Synchrotrons Oliver Boine-Franenheim TU Darmstadt and GSI Contents o Application and challenges for beam simulations with space charge o D electrostatic Particle-In-Cell
More informationFokker-Planck Equation with Detailed Balance
Appendix E Fokker-Planck Equation with Detailed Balance A stochastic process is simply a function of two variables, one is the time, the other is a stochastic variable X, defined by specifying: a: the
More informationÜbungen zur Elektrodynamik (T3)
Arnold Sommerfeld Center Ludwig Maximilians Universität München Prof. Dr. Ivo Sachs SoSe 08 Übungen zur Elektrodynamik (T3) Lösungen zum Übungsblatt 7 Lorentz Force Calculate dx µ and ds explicitly in
More informationProblem 1: Lagrangians and Conserved Quantities. Consider the following action for a particle of mass m moving in one dimension
105A Practice Final Solutions March 13, 01 William Kelly Problem 1: Lagrangians and Conserved Quantities Consider the following action for a particle of mass m moving in one dimension S = dtl = mc dt 1
More informationPhysics 208, Spring 2016 Exam #3
Physics 208, Spring 206 Exam #3 A Name (Last, First): ID #: Section #: You have 75 minutes to complete the exam. Formulae are provided on an attached sheet. You may NOT use any other formula sheet. You
More informationSimulation of transverse emittance measurements using the single slit method
Simulation of transverse emittance measurements using the single slit method Rudolf Höfler Vienna University of Technology DESY Zeuthen Summer Student Program 007 Abstract Emittance measurements using
More informationName: Fall 2014 CLOSED BOOK
Name: Fall 2014 1. Rod AB with weight W = 40 lb is pinned at A to a vertical axle which rotates with constant angular velocity ω =15 rad/s. The rod position is maintained by a horizontal wire BC. Determine
More information1 Introduction. 2 Diffusion equation and central limit theorem. The content of these notes is also covered by chapter 3 section B of [1].
1 Introduction The content of these notes is also covered by chapter 3 section B of [1]. Diffusion equation and central limit theorem Consider a sequence {ξ i } i=1 i.i.d. ξ i = d ξ with ξ : Ω { Dx, 0,
More informationControlled Diffusions and Hamilton-Jacobi Bellman Equations
Controlled Diffusions and Hamilton-Jacobi Bellman Equations Emo Todorov Applied Mathematics and Computer Science & Engineering University of Washington Winter 2014 Emo Todorov (UW) AMATH/CSE 579, Winter
More informationPenning Traps. Contents. Plasma Physics Penning Traps AJW August 16, Introduction. Clasical picture. Radiation Damping.
Penning Traps Contents Introduction Clasical picture Radiation Damping Number density B and E fields used to increase time that an electron remains within a discharge: Penning, 936. Can now trap a particle
More informationOn the numerical integration of variational equations
On the numerical integration of variational equations Haris Skokos Max Planck Institute for the Physics of Complex Systems Dresden, Germany E-mail: hskokos@pks.mpg.de, URL: http://www.pks.mpg.de/~hskokos/
More informationDynamical Systems and Mathematical Models A system is any collection of objects and their interconnections that we happen to be interested in; be it
Dynamical Systems and Mathematical Models A system is any collection of objects and their interconnections that we happen to be interested in; be it physical, engineering, economic, financial, demographic,
More informationG : Statistical Mechanics
G25.2651: Statistical Mechanics Notes for Lecture 15 Consider Hamilton s equations in the form I. CLASSICAL LINEAR RESPONSE THEORY q i = H p i ṗ i = H q i We noted early in the course that an ensemble
More informationStochastic Hamiltonian systems and reduction
Stochastic Hamiltonian systems and reduction Joan Andreu Lázaro Universidad de Zaragoza Juan Pablo Ortega CNRS, Besançon Geometric Mechanics: Continuous and discrete, nite and in nite dimensional Ban,
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 informationUnimodularity and preservation of measures in nonholonomic mechanics
Unimodularity and preservation of measures in nonholonomic mechanics Luis García-Naranjo (joint with Y. Fedorov and J.C. Marrero) Mathematics Department ITAM, Mexico City, MEXICO ẋ = f (x), x M n, f smooth
More informationStatistical Mechanics of Active Matter
Statistical Mechanics of Active Matter Umberto Marini Bettolo Marconi University of Camerino, Italy Naples, 24 May,2017 Umberto Marini Bettolo Marconi (2017) Statistical Mechanics of Active Matter 2017
More informationSimulation of Coulomb Collisions in Plasma Accelerators for Space Applications
Simulation of Coulomb Collisions in Plasma Accelerators for Space Applications D. D Andrea 1, W.Maschek 1 and R. Schneider 2 Vienna, May 6 th 2009 1 Institut for Institute for Nuclear and Energy Technologies
More informationThe 3-wave PDEs for resonantly interac7ng triads
The 3-wave PDEs for resonantly interac7ng triads Ruth Mar7n & Harvey Segur Waves and singulari7es in incompressible fluids ICERM, April 28, 2017 What are the 3-wave equa.ons? What is a resonant triad?
More information