Thu June 16 Lecture Notes: Lattice Exercises I

Similar documents
USPAS Accelerator Physics 2013 Duke University

USPAS Accelerator Physics 2017 University of California, Davis

Introduction to Transverse Beam Optics. II.) Twiss Parameters & Lattice Design

LECTURE 7. insertion MATCH POINTS. Lattice design: insertions and matching

Introduction to Accelerator Physics 2011 Mexican Particle Accelerator School

E. Wilson - CERN. Components of a synchrotron. Dipole Bending Magnet. Magnetic rigidity. Bending Magnet. Weak focusing - gutter. Transverse ellipse

FODO Cell Introduction to OptiM

Particle Accelerators: Transverse Beam Dynamics

Equations of motion in an accelerator (Lecture 7)

Accelerator Physics Homework #3 P470 (Problems: 1-5)

Introduction to Accelerator Physics 2011 Mexican Particle Accelerator School

Definition: Let f(x) be a function of one variable with continuous derivatives of all orders at a the point x 0, then the series.

Lattice Design in Particle Accelerators

Beam Dynamics. D. Brandt, CERN. CAS Bruges June 2009 Beam Dynamics D. Brandt 1

Wed Jan 25 Lecture Notes: Coordinate Transformations and Nonlinear Dynamics

Accelerator Physics Final Exam pts.

Lattice Design II: Insertions Bernhard Holzer, DESY

S9: Momentum Spread Effects and Bending S9A: Formulation

RESOLUTION MSC.362(92) (Adopted on 14 June 2013) REVISED RECOMMENDATION ON A STANDARD METHOD FOR EVALUATING CROSS-FLOODING ARRANGEMENTS

Simulation of Optics Correction for ERLs. V. Sajaev, ANL

TWISS FUNCTIONS. Lecture 1 January P.J. Bryant. JUAS18_01- P.J. Bryant - Lecture 1 Twiss functions

Transverse dynamics Selected topics. Erik Adli, University of Oslo, August 2016, v2.21

The Deutsch-Jozsa Problem: De-quantization and entanglement

Practical Lattice Design

Part I: Thin Converging Lens

Bernhard Holzer, CERN-LHC

Transverse dynamics. Transverse dynamics: degrees of freedom orthogonal to the reference trajectory

Introduction to Transverse Beam Dynamics

( 1) ( 2) ( 1) nan integer, since the potential is no longer simple harmonic.

Lattice Design and Performance for PEP-X Light Source

Supplementary material for Continuous-action planning for discounted infinite-horizon nonlinear optimal control with Lipschitz values

(C) The rationals and the reals as linearly ordered sets. Contents. 1 The characterizing results

Transverse Dynamics II

Accelerator Physics. Elena Wildner. Transverse motion. Benasque. Acknowldements to Simon Baird, Rende Steerenberg, Mats Lindroos, for course material

06. Orbit Stability and the Phase Amplitude Formulation *

Low Emittance Machines

Lecture 2: Modeling Accelerators Calculation of lattice functions and parameters. X. Huang USPAS, January 2015 Hampton, Virginia

Mathematical Notation Math Calculus & Analytic Geometry III

Magnets and Lattices. - Accelerator building blocks - Transverse beam dynamics - coordinate system

TLT-5200/5206 COMMUNICATION THEORY, Exercise 3, Fall TLT-5200/5206 COMMUNICATION THEORY, Exercise 3, Fall Problem 1.

accelerator physics and ion optics summary longitudinal optics

9.1 The Square Root Function

Practical Lattice Design

Accelerator Physics Multipoles and Closed Orbits. G. A. Krafft Old Dominion University Jefferson Lab Lecture 13

Transverse Beam Dynamics II

Injection, extraction and transfer

CI Courses / Winter 2016

D. Brandt, CERN. CAS Frascati 2008 Accelerators for Newcomers D. Brandt 1

Chapter 4 Imaging. Lecture 21. d (110) Chem 793, Fall 2011, L. Ma

Beam Transfer Lines. Brennan Goddard CERN

07. The Courant Snyder Invariant and the Betatron Formulation *

Roberto s Notes on Differential Calculus Chapter 8: Graphical analysis Section 1. Extreme points

S5: Linear Transverse Particle Equations of Motion without Space Charge, Acceleration, and Momentum Spread S5A: Hill's Equation

S5: Linear Transverse Particle Equations of Motion without Space Charge, Acceleration, and Momentum Spread S5A: Hill's Equation

Bernhard Holzer, CERN-LHC

Accelerator Physics Homework #7 P470 (Problems: 1-4)

( x) f = where P and Q are polynomials.

Class Review. Content

Nonlinear Single-Particle Dynamics in High Energy Accelerators

Conceptual design of an accumulator ring for the Diamond II upgrade

Feasibility of a Multi-Pass Thomson Scattering System with Confocal Spherical Mirrors

Scattering of Solitons of Modified KdV Equation with Self-consistent Sources

Mathematical Notation Math Calculus & Analytic Geometry III

The A, B, C and D are determined by these 4 BCs to obtain

between electron energy levels. Using the spectrum of electron energy (1) and the law of energy-momentum conservation for photon absorption:

New Functions from Old Functions

Physics 2B Chapter 17 Notes - First Law of Thermo Spring 2018

Review D: Potential Energy and the Conservation of Mechanical Energy

Feedback Linearization

accelerator physics and ion optics summary longitudinal optics

Accelerator Physics Closed Orbits and Chromaticity. G. A. Krafft Old Dominion University Jefferson Lab Lecture 14

Parametrization of the Driven Betatron Oscillation

Bernhard Holzer, CERN-LHC

S1: Particle Equations of Motion S1A: Introduction: The Lorentz Force Equation

Point Spread Function of Symmetrical Optical System Apodised with Gaussian Filter

Nonlinear Single-Particle Dynamics in High Energy Accelerators

Chapter 2 Laser Light Amplification by Stimulated Emission of Radiation

9.3 Graphing Functions by Plotting Points, The Domain and Range of Functions

Lecture 8 Optimization

Beam-Based Measurement of Dynamical Characteristics in Nuclotron

Lattices and Emittance

Your Comments. I was not sure about the angle at which the beam reflects from a concave mirror.

. This is the Basic Chain Rule. x dt y dt z dt Chain Rule in this context.

Low Emittance Machines

Objectives. By the time the student is finished with this section of the workbook, he/she should be able

Lattice Design for the Taiwan Photon Source (TPS) at NSRRC

Christoffel symbols and Gauss Theorema Egregium

arxiv: v1 [math.gt] 20 Feb 2008

Physics 598ACC Accelerators: Theory and Applications

SEPARATED AND PROPER MORPHISMS

Lattice Design II: Insertions Bernhard Holzer, CERN

Poisson Brackets and Lie Operators

Wind-Driven Circulation: Stommel s gyre & Sverdrup s balance

Numerical Methods - Lecture 2. Numerical Methods. Lecture 2. Analysis of errors in numerical methods

Math Review and Lessons in Calculus

Pre-AP Physics Chapter 1 Notes Yockers JHS 2008

Finite Dimensional Hilbert Spaces are Complete for Dagger Compact Closed Categories (Extended Abstract)

Lecture : Feedback Linearization

Aggregate Growth: R =αn 1/ d f

Polarization Preservation and Control in a Figure-8 Ring

Transcription:

Thu June 6 ecture Notes: attice Exercises I T. Satogata: June USPAS Accelerator Physics Most o these notes ollow the treatment in the class text, Conte and MacKay, Chapter 6 on attice Exercises. The portions on orbit correction ollow Edwards and Syphers, section 3.4. By the end o these lectures, you should understand a bit about orbit correction in a circular accelerator, how to calculate optics parameters or a FODO lattice and other combinations o drits, dipoles, and quads, FODO lattice stability, and FODO lattice dispersion. Review Rcall that the parameterization o the one-dimensional periodic transer matrix is given by Conte/MacKay 5.-: where M periodic I cos µ + J sin µ e µj. αs βs J γs αs J I. αs β s/ and γs + αs /βs. Here I ve included the explicit s-dependence. βs is the square o the envelope unction ws, or amplitude scaling o particle motion, Conte/Mackay 5.69: xs Wβs cos[ψs + ψ ].3 where ψ s /βs ψs ds/βs βs βs + C.4 The beta unction has the periodicity o the lattice, as Waldo showed, but the phase advance ψs does NOT have this periodicity integrating it around the accelerator gives us an extra phase advance o µ πq where Q is called the betatron tune. See Conte/Mackay Eqns. 5.54 55. There are beta unctions, phase advances, and betatron tunes or each plane H, V. Note that some reerences use ν or tunes instead o Q and that since the transer matrix around one turn o the accelerator is periodic, we have M one turn I cos πq + J sin πq e πqj.5 In the weak ocusing betatron case where the ield index is n, Q x n and Q y n, the tunes must also be in the range Q x,y. Indeed, this is practically what is deined as weak ocusing. In contrast, accelerators with alternating quadrupole ocusing that lead to Q > are termed strong ocusing. To propagate position and angle rom one location s s to another s location, both

with known Twiss parameters, we use MacKay and Conte 5.5: βs β Ms s [cos µs + α sin µs] β βs sin µs [αs α ] cos µs+[+α αs] sin µs β.6 [cos µs αs sin µs] β βs βs βs cos µ sin µ β βs α sin µ cos µ β.7 β αs βs T s RµT.8 Even though.6 looks horrible, it can be decomposed into local transormations T s into normalized phase space, and a separate pure rotation matrix Rµ that rotates through the phase advance. You ll see this in problem 5-4 o tonight s homework. Figure shows particle motion through a FODO lattice, a series o alternating ocusing and deocusing quadrupoles. The dark lines that are plotted are ± βs, which is directly proportional to the maximum amplitude o the particle oscillation. You can see why the beta unction is also called the envelope unction! Waldo s showed something like this and you ve also seen it in the labs. et s analyze the FODO lattice in more detail and calculate some things about it ater a brie digression into a real application o what you ve learned so ar. Figure : This igure shows a plot o betatron oscillations, equation.3 and Conte/MacKay Fig. 5.3, through a regular FODO lattice. Waldo showed this picture earlier. I you ollow the trajectory, you can see that the kicks alternate between ocusing and deocusing, and the betatron oscillation is slow compared to the periodicity o the FODO lattice. This example has about 6.5 FODO cells per betatron oscillation, so this FODO lattice has a phase advance per cell o about π/6.5 55 degrees. What is the relationship between the ocal length and the FODO cell length here? Steering Errors and Orbit Correction Suppose we have a single steering error in an otherwise perect circular accelerator that gives a kick x B l. Bρ loated without loss o generality at s. What is the new closed -turn periodic orbit? Say that the orbit immediately downstream o this kick is x, x. I we propagate this orbit through one turn with a one-turn matrix M one turn and then proceed through the kick, we should get the same conditions. This the deinition o closed orbit, the ixed point o the one-turn map. So M one turn x x + x x x..

Solving or x, x gives x x I M one turn x..3 We can now use M one turn e πqj to perorm a little trick: I M one turn I e πqj [ e πqj e πqj e πqj] The closed orbit solution is then x I M one turn x x sin πq x.4 J sin πq e πqj.5 sin πq Je πqj.6 J cos πq + I sin πq sin πq.7 β cos πq sin πq α cos πq.8 where the twiss parameters β, α are at the position o the kick. We can then propagate this to any location s in the ring with the general Ms s matrix.6 to ind xs x β / sβ / sin πq cos[ψs ψ πq].9 xs x x s. This result really represents a standing wave solution through the accelerator. Note some interesting points about.9: The sensitivity to the kick is proportional to the square root o the product o the beta unctions at the kick location and position observation. Dipole errors have the worst eect at areas o high β. The orbit change is inversely proportional to sin πq. For integer Q, dipole error eects diverge. This is an example o a resonance. We will later see that quadrupole error eects are inversely proportional to sin πq and so on to higher orders. Even though the kick is only an angle kick, there is a closed orbit change x at the kick location. An example o a comparison o measurements and a calculation rom.9 or a single dipole corrector change in RHIC is shown in Fig...9 can be used as a undamental equation o orbit correction since it relates a position change anywhere in the ring xs to a change in a dipole corrector change x elsewhere in the ring. Because it is linear we can then write x x x n x x x x x n x x x x x x x m x x m x n x n x n x m x x x m. This non-square! matrix can be inverted optimally using, or example, singular value decomposition to ind a solution to the changes in all m dipole correctors x, x,, x m to eect an orbit change in all n beam position monitor locations x, x,, x n. 3

Figure : attice diagram and vertical dierence orbit rom RHIC, comparing measurement and modeled dierence orbit rom.9. Figure 3: Two ways o calculating transer or a regular FODO cell o length with two thin quadrupoles. The one on the let is not symmetric and gives optics at one side o the ocusing quad. The one on the right is symmetric and gives optics in the center o the ocusing quadrupole, where we expect to locate β max even in the case o thick quadrupoles. 3 FODO attice One o the most common lattice layouts in modern separated-unction accelerators is the FODO lattice, a lattice with alternating ocusing and deocusing quadrupoles separated by drit spaces or dipoles. We have started to see how this combination gives you a net ocusing. Using the thin lens approximation or equal-strength quadrupoles i.e. assuming their ocal lengths are much larger than their lengths, separating them by / length drit spaces so the FODO cell length is, we ind: / M FODO 4 + 4 + / 3. 3. et s assume that we have a string o FODO cells. The system is periodic through every FODO cell, so the transer matrix o each FODO cell must be expressible in the orm o.: 4 + 4 + cos µ + α sin µ β sin µ γ sin µ cos µ α sin µ 3.3 4

I these are equal, their traces must be equal: 4 cos µ 4 sin µ 3.4 sin µ ± 4 3.5 using the trig identity cosµ cos µ sin µ sin µ. 3.5 relates the FODO cell phase advance µ to the parameters and that we used to construct it. For a FODO lattice, the maximum o the beta unction β β max is located at the center o the ocusing quadrupole. The transer matrix rom the center o one ocusing quadrupole with ocal length to another in this lattice is ound by splitting the ocusing quadrupole in hal ocal length as shown in Fig. 3. This even works i we use the thick lens quadrupole transer matrix. M to / 4 / / / / + 4 3.6 4 8 8 At the maximum o the beta unction in the split ocusing quadrupole, αs β s/, so this matrix can also be written as cos µ β M to max sin µ 3.7 sin µ/β max cos µ Equating irst diagonal, then upper-right o-diagonal terms gives cos µ cos µ /8 sinµ/ 4 + sinµ/ β max sin µ 3.8 The minimum o the beta unction, β min, occurs at the center o the deocusing quadrupole, so it may be ound by ollowing the above analysis replacing with, producing: sinµ/ β min 3.9 sin µ For example, the phase advance per FODO cell at RHIC is about 78 degrees and the length o the RHIC FODO cell is about 3 m, so we can calculate that β max 5 m and β min.7 m. Figure 4 shows one sixth o the RHIC injection optics lattice, and you can readily see that the ormulas or β min,max work very well. Note that we can calculate the Twiss parameters at any point in the FODO cell by changing the cut point or our lattice transport. Alternatively you can use the result o tonight s homework problem 5-5 to propagate these Twiss unctions through the FODO cell. Propagating through an entire FODO cell should give an identity matrix regardless o the cut point because the optics are periodic. We can test this or our example above using M to rom 3.6. This is thankully NOT part o your homework. But it works: β α γ M to M M M M M M M M + M M M M M M M M /8 + /4 /6 3 /4 /8 β α γ 3. 3. 5

Figure 4: The RHIC injection lattice and beta unctions. This shows one sixth o the RHIC lattice as Waldo has shown earlier; apply mirror antisymmetry and repeat the lattice three times to get the ull lattice. Solid lines are β x s; dotted lines are β y s. The magnet lattice is shown at the top o the graph; short blocks are bending dipoles; taller blocks are quadrupoles. ow-beta insertions are clearly visible. Is the quadrupole at s3m horizontally ocusing or deocusing? How many FODO cells are in the regular section? RHIC has a total horizontal tune o about 8.3; does this number o FODO cells make sense? 4 FODO attice Stability Sometimes we do not want horizontal and vertical quadrupoles in a FODO lattice to have the same strengths, such as when we are matching the beam size into dierent optics. I we re-parameterize the quadrupole strengths as F F D D 4. Calculating the trace o the FODO lattice matrix now gives something more general than 3.5: cos µ + D F F D sin µ F D 4 + F D 4. Recall rom Waldo s lecture that the stability requirement is that the phase advance µ is real since the eigenvalues o the transport matrix are e ±iµ, or < cos µ < or < sin µ < 4.3 What are the boundaries o real µ in terms o F and D? One is where sin µ, and we have F D 4.4 + D The second case, sin µ, requires that F < or stability or positive D. The reverse conditions where we reverse the roles o F and D also provide boundaries, and the inal stable area in the F, D parameter space is shown in Fig. 5. 5 FODO Cell Dispersion A FODO cell is a linear system in more than one way the design orbit is straight! We want to steer the beam and eventually wrap the beamline around back on itsel so we can 6

Figure 5: Stability or necktie diagram or an alternate ocusing lattice, including the FODO lattice where the dimensionless parameterization is F F and D D. The shaded area is the region o stability. Figure 6: The π lattice insertion. make a storage ring, so we also must consider cases where the space between the FODO quadrupoles is not a simple drit, but illed with the bending ield o a dipole. To make the math simpler and to save my sanity, we ll assume that the length o the FODO cell is ρθ c and there are no drits. The total bending angle o the cell is θ c. We will use the symmetric layout o the FODO cell in Fig. to make some math simpler or dispersion suppressors, which we ll discuss later. The transer matrix or the single FODO cell with the dispersive component included is M to θ c 8 θ c + 8 4 + 8 θ c 4 4 θ 8 3 c 8 θ c 8 θ c 5. 5. where we have kept only irst order in θ c since usually θ c <<. Waldo talked a bit about the dispersion unction ηs earlier. The periodic dispersion unction o the FODO cell can be obtained rom η η M M M 3 M M M 3 η η 5.3 M 3 M 3 M 33 The bottom row o M in 5. is trivial, so we can also write the dispersion as a nonhomogeneous equation: η η η M M3 η + 5.4 M 3 η η I M M3 5.5 M 3 7

where M is the upper let x block o 5.. This is calculable and gives the η max, just like we had β max : η max θ c + µ 4 sin µ 5.6 η 5.7 This is rather unsurprisingly a lot like the ormula or β max, except now the maximum dispersion also scales linearly with the total length o the FODO cell, θ c. The minimum o the dispersion is also easy to ind: η min θ c µ 4 sin µ 5.8 η 5.9 6 A Simple attice Insertion In any FODO cell one or both quadrupole magnets can be moved away leaving room or actual straight sections, slightly modiying the lattice dispersion without altering the optics o the lattice, as thoroughly discussed in the previous sections. Obviously, such ree spaces may not be suicient to host long items like septum magnets or injection and extraction, chains o r cavities, etc. onger straight sections have to be provided, leaving the optics o the machine unaected. How do we make the longest straight section possible? An elegant solution is the so-called π insertion. This insertion has a long straight section o length, encompassed by two quadrupoles, one ocusing and one deocusing, and urther symmetric straight sections o length. This is schematically shown in igure 6. Assuming the quadrupoles have equal ocal length, we calculate the transer matrix again: + M + 6. We want this insertion to match into the periodicity o the lattice, so it must match to our general periodic lattice description o.: cos µ + α sin µ β sin µ M I cos µ + J sin µ 6. γ sin µ cos µ α sin µ Taking the dierence o the diagonal terms gives α sin µ 6.3 and is maximized when µ π. It s now obvious why this is called a π insertion. Then cos µ and we can get other parameters in terms o the lengths and : α γ β + 6.4 and we can design our insertion to match into our FODO lattice. One surprising observation is that the transer matrix now reduces to M J 6.5 and recall that M I! So we can now design an interesting insertion by putting two o these insertions back to back. This is a simple case o a low-β insertion. 8