Position Dependence Of The K Parameter In The Delta Undulator
|
|
- Abel Benson
- 5 years ago
- Views:
Transcription
1 LCLS-TN-16-1 Position Dependence Of The K Parameter In The Delta Undulator Zachary Wolf SLAC January 19, 016 Abstract In order to understand the alignment tolerances of the Delta undulator, we must know how the K value depends on the transverse position of the beam. In this note we calculate the K value as a function of beam position for di erent polarization modes and row shift settings. The e ect of transverse eld rollo of the quadrants is included, and we will see that it plays a large role in the position dependence of K. In particular, it gives a gradient to K at the magnetic center in the elliptical modes, potentially making the undulator alignment requirements di cult to meet. 1 Introduction 1 The position dependence of the K parameter in the Delta undulator determines how accurately the undulator must be aligned in order to be resonant with other undulators. In an EPU the K parameter depends on position di erently for di erent quadrant row settings. In addition, the K parameter depends on the shape of the quadrant elds. In this note we explore these dependencies. The organization of the note is as follows. We rst calculate the magnetic elds in the Delta undulator. This is done in a rotated coordinate frame relative to the laboratory frame in which the undulator is used. Rotating the coordinates does not a ect the calculated K value and it simpli es the calculations. From the elds, we calculate the trajectory slopes which we use to calculate the slippage. The slope of the linear part of the slippage change with longitudinal position lets us calculate the K value of the undulator. The K value is calculated at di erent transverse positions and at di erent undulator row settings in order to study the position dependence of K. Magnetic Field In The Undulator.1 The Scalar Potential The scalar potential from a single magnet array satis es Laplace s equation. It varies sinusoidally down the array, it decreases away from the array, and it falls o as one moves away from the center line of the array in the transverse direction. In the coordinate system shown in gure 1, the scalar potential can be expressed as a Fourier sum of terms with a fundamental term of the form = 0 cos( s) exp ( r) cos ( (z z 0 )) (1) 1 Worupported in part by the DOE Contract DE-AC0-76SF This work was performed in support of the LCLS project at SLAC. 1
2 where 0 is a constant, = = u where u is the undulator period, z is the coordinate down the undulator, and z 0 gives the quadrant position along z. and determine the behavior of the potential in the transverse directions. The potential decreases as one moves radially away from the magnet and this is expressed by exp ( r). The potential also decreases as one moves to the side and this is given by the cos( s) dependence over a limited range of s. Since the Laplacian of the potential is zero, we have the constraint ks + kr ku = 0 () Figure 1: Coordinate system for the scalar potential of a single quadrant. Although the scalar potential is a sum of terms, we assume that the fundamental term dominates the K value. In a Fourier expansion of this type, typically the amplitude of the fundamental term is the largest. Furthermore, harmonics with a cos (n (z z 0 )) dependence give contributions to the transverse velocity of the beam that go as 1=n, and therefore the contributions to K from these harmonics go as 1=n. The smaller amplitudes of the higher harmonics combined with the 1=n dependence lead us to expect their contribution to K will be small. The scalar potential also has terms with di erent values (and di erent values as given by the constraint). Satisfactory ts of B r as a function of s above a magnet array can be done with only a few cos( s) terms. The fundamental term is the largest and has a value of determined by the e ective width of the magnet block. Additional terms have values which are not integer multiples of the fundamental, but are determined by the detailed magnet shape. These higher order terms are typically much smaller than the fundamental. The higher order terms may be useful to construct a precise eld behavior by shaping the magnet. But for now, we consider only the fundamental term which dominates K. We use the form of the potential given in equation 1 to calculate the magnetic scalar potential in the undulator by rotating the four quadrants and summing their rotated scalar potentials. The Delta undulator is oriented as shown on the left side of gure. In the laboratory, z is along the beam direction, y L is up, and x L makes a right handed system. For our calculations, it is more convenient to use the rotated system x, y, z, where x is in the direction from quadrant 3 to quadrant 1, y is in the direction from quadrant 4 to quadrant, and z is in the beam direction. Using equation 1, the scalar potential for each of the quadrants in the x, y, z system is 1 (x; y; z) = 0Q cos( y) exp ( x) cos ( (z z 01 )) (3) (x; y; z) = 0Q cos( x) exp ( y) cos ( (z z 0 )) (4) 3 (x; y; z) = 0Q cos( y) exp ( x) cos ( (z z 03 )) (5) 4 (x; y; z) = 0Q cos( x) exp ( y) cos ( (z z 04 )) (6)
3 Figure : The left side of the gure shows the Delta undulator in its con guration in the tunnel where y L is up, z is in the beam direction, and x L makes a right handed system. For our calculations, we use the rotated coordinate system on the right, where x is along the line pointing from quadrant 3 to quadrant 1, y is along the line pointing from quadrant 4 to quadrant, and z makes a right handed system. where z 0i is the longitudinal shift of quadrant i, and 0Q is the amplitude of the scalar potential of all the identical quadrants on the axis of the undulator where x = 0 and y = 0. Quadrants 3 and 4 are loaded with opposite polarity magnets as quadrants 1 and in order to make a vertical eld planar undulator in the laboratory frame when all the rows are aligned. This accounts for the minus signs, 0Q, in the potentials for quadrants 3 and 4. The scalar potential for the undulator is the sum of the scalar potentials for the quadrants. Quadrants 1 and 3 both depend on x, and quadrants and 4 both depend on y. We rst add the scalar potentials for quadrants 1 and 3, and then add the scalar potentials for quadrants and 4, and then add the sums to get the scalar potential for the whole undulator. We will interpret this as forming the entire undulator from two crossed planar adjustable phase undulators. The scalar potential for the combination of quadrants 1 and 3 is given by Let = 0Q cos( y) exp ( x) cos ( (z z 01 )) 0Q cos( y) exp ( x) cos ( (z z 03 )) (7) z 01 = Z + z 03 = Z where Z = z 01 + z 03 is the average z-position of the quadrants, and is the z-shift between the quadrants. quadrants becomes = 0Q cos( y) sinh ( x) cos (8) (9) (10) = z 01 z 03 (11) With these de nitions, the scalar potential for the pair of + 0Q cos( y) cosh ( x) sin 3 cos ( (z Z )) sin ( (z Z )) (1)
4 This is the scalar potential for a planar adjustable phase undulator. The full range of amplitudes and phases are covered if the range of includes u ; u and the range of Z includes u ; u. Similarly, the scalar potential for the combination of quadrants and 4 is given by Let 4 = 0Q cos( x) exp ( y) cos ( (z z 0 )) 0Q cos( x) exp ( y) cos ( (z z 04 )) () z 0 = Z z 04 = Z 4 4 where Z 4 = z 0 + z 04 is the average z-position of the quadrants, and (14) (15) (16) 4 = z 0 z 04 (17) is the z-shift between the quadrants. With these de nitions, the scalar potential for the pair of quadrants becomes 4 4 = 0Q cos( x) sinh ( y) cos cos ( (z Z 4 )) 4 + 0Q cos( x) cosh ( y) sin sin ( (z Z 4 )) (18) This is again the potential for a planar adjustable phase undulator. The scalar potential for the undulator is the sum of the scalar potentials for the quadrant pairs. = + 4 (19) Performing the sum, we nd = 0Q cos( y) sinh ( x) cos cos ( (z Z )) + 0Q cos( y) cosh ( x) sin sin ( (z Z )) 4 + 0Q cos( x) sinh ( y) cos cos ( (z Z 4 )) 4 + 0Q cos( x) cosh ( y) sin sin ( (z Z 4 )) (0) By putting in the various values for Z, Z 4,, and 4, we get the scalar potential for the various undulator modes at di erent K values. Z. Wolf, "Variable Phase PPM Undulator Study", LCLS-TN-11-1, May,
5 . Magnetic Field In The Undulator Coordinate System The magnetic eld in the undulator is given by B = r. Taking the gradient we nd B x (x; y; z) = 0Q [ k r cos cos( y) cosh ( x) cos ( (z Z )) + k r sin cos( y) sinh ( x) sin ( (z Z )) 4 cos sin( x) sinh ( y) cos ( (z Z 4 )) sin 4 sin( x) cosh ( y) sin ( (z Z 4 ))] (1) B y (x; y; z) = 0Q [ cos sin( y) sinh ( x) cos ( (z Z )) sin sin( y) cosh ( x) sin ( (z Z )) + k r 4 cos cos( x) cosh ( y) cos ( (z Z 4 )) + k r 4 sin cos( x) sinh ( y) sin ( (z Z 4 ))] () B z (x; y; z) = 0Q [ cos + sin cos + sin 4 4 In order to simplify these formulas further, let cos( y) sinh ( x) sin ( (z Z )) cos( y) cosh ( x) cos ( (z Z )) cos( x) sinh ( y) sin ( (z Z 4 )) cos( x) cosh ( y) cos ( (z Z 4 ))] (3) B 0 = 0Q (4) B xc = B 0 cos (5) B xs = B 0 sin (6) 4 B yc = B 0 cos (7) 4 B ys = B 0 sin (8) 5
6 With these substitutions, the elds become B x (x; y; z) = B xc cos( y) cosh ( x) cos ( (z Z )) + B xs cos( y) sinh ( x) sin ( (z Z )) B yc sin( x) sinh ( y) cos ( (z Z 4 )) B ys sin( x) cosh ( y) sin ( (z Z 4 )) (9) B y (x; y; z) = B xc sin( y) sinh ( x) cos ( (z Z )) B xs sin( y) cosh ( x) sin ( (z Z )) + B yc cos( x) cosh ( y) cos ( (z Z 4 )) + B ys cos( x) sinh ( y) sin ( (z Z 4 )) (30) B z (x; y; z) = B xc cos( y) sinh ( x) sin ( (z Z )) We can change the origin of z such that If we let +B xs cos( y) cosh ( x) cos ( (z Z )) B yc cos( x) sinh ( y) sin ( (z Z 4 )) and drop the prime, we write the elds as a function of position as +B ys cos( x) cosh ( y) cos ( (z Z 4 )) (31) z = z 0 + Z (3) = (Z Z 4 ) (33) B x (x; y; z) = B xc cos( y) cosh ( x) cos ( z) +B xs cos( y) sinh ( x) sin ( z) B yc sin( x) sinh ( y) cos ( z + ) B ys sin( x) cosh ( y) sin ( z + ) (34) B y (x; y; z) = B xc sin( y) sinh ( x) cos ( z) B xs sin( y) cosh ( x) sin ( z) +B yc cos( x) cosh ( y) cos ( z + ) +B ys cos( x) sinh ( y) sin ( z + ) (35) 6
7 B z (x; y; z) = B xc cos( y) sinh ( x) sin ( z) 3 Trajectories In The Undulator 3.1 Equations Of Motion +B xs cos( y) cosh ( x) cos ( z) B yc cos( x) sinh ( y) sin ( z + ) +B ys cos( x) cosh ( y) cos ( z + ) (36) The trajectory of a charged particle beam in the undulator is determined by the Lorentz force law dp dt = q (v B) (37) where p = mv, = 1 1, and m is the particle rest mass. The energy of the particle is constant in the magnetic eld of the undulator, except for radiation losses which we neglect. With constant energy, is constant. The Lorentz force law becomes v = q (v B) (38) m With the substitution d=dt = v z d=, the changes in the individual velocity components with z are given by dv x dv y dv z = = = q (v y B z mv z v z B y ) (39) q (v z B x mv z v x B z ) (40) q (v x B y mv z v y B x ) (41) 3. Iterative Solution We solve these equations iteratively by expanding in powers of a small parameter. Let = qb 0 m c (4) be the small dimensionless expansion parameter. motion are dv x dv y dv z In terms of this parameter, the equations of = c B 0 v z (v y B z v z B y ) (43) = c B 0 v z (v z B x v x B z ) (44) = c B 0 v z (v x B y v y B x ) (45) At this point we could divide all velocities by c, divide all positions by 1=, and divide all magnetic elds by B 0 in order to eliminate the factor kuc B 0v z. We think of small quantities in terms of these dimensionless variables. In the interest of clarity, however, we keep the factor and continue to use mks units, so there is no confusion about the meaning of a quantity. 7
8 Expanding the particle position in terms of the expansion parameter, we have x = x (0) + x (1) + x () + (46) y = y (0) + y (1) + y () + (47) z = z (48) where the last line indicates that the z-position of the particle is the independent variable. numbers in parenthesis indicate the order of the expansion. The velocities are where The magnetic elds are given by The v x = v x(0) + v x(1) + v x() + (49) v y = v y(0) + v y(1) + v y() + (50) v z = v z(0) + v z(1) + v z() + (51) B x = B x j 0 x B x j 0 x (1) y B x j 0 y (1) + (5) B y = B y j 0 x B y j 0 x (1) y B y j 0 y (1) + (53) B z = B z j 0 x B z j 0 x (1) y B z j 0 y (1) + (54) B i j 0 = B i (x (0) ; y (0) ; z) j B i j 0 j (x (0) ; y (0) ; z) (56) where i = x; y; z, j = x; y. The transverse positions are found from the velocities as follows: v x = dx dt x(z) = x 0 + = x 0 + Z z 0 Z z = x 0 + v x(0) v z(0) z + 0 v x (z 0 1 ) v z (z 0 ) 0 v x(0) + v x(1) + 1 Z z 0 v x(1) v z(0) (57) 1 v z(1) + 0 v z(0) v z(0)! v x(0) v z(1) vz(0) 0 + (58) In these equations x 0 is the constant initial position where the beam enters the undulator. Similarly y(z) = y 0 + v Z! z y(0) v y(1) v y(0) v z(1) z + v z(0) 0 v z(0) vz(0) 0 + (59) We see that and x (1) = y (1) = x (0) = x 0 + v x(0) v z(0) z (60) y (0) = y 0 + v y(0) v z(0) z (61) Z z 0 Z z 0 v x(1) v z(0) v y(1) v z(0) v x(0) v z(1) v z(0) v y(0) v z(1) v z(0)!! 0 (6) 0 (63) 8
9 With the expansions in the small parameter, the equations of motion become dv x(0) + dv x(1) + dv x() + c = 1 v z(1) + B 0 v z(0) v z(0) f(v y(0) + v y(1) + )(B z j 0 x B z j 0 x (1) y B z j 0 y (1) + ) (v z(0) + v z(1) + )(B y j 0 x B y j 0 x (1) y B y j 0 y (1) + )g (64) dv y(0) + dv y(1) + dv y() + c = 1 v z(1) + B 0 v z(0) v z(0) f(v z(0) + v z(1) + )(B x j 0 x B x j 0 x (1) y B x j 0 y (1) + ) (v x(0) + v x(1) + )(B z j 0 x B z j 0 x (1) y B z j 0 y (1) + )g (65) dv z(0) + dv z(1) + dv z() + c = 1 v z(1) + B 0 v z(0) v z(0) f(v x(0) + v x(1) + )(B y j 0 x B y j 0 x (1) y B y j 0 y (1) + ) (v y(0) + v y(1) + )(B x j 0 x B x j 0 x (1) y B x j 0 y (1) + )g (66) We now proceed to solve these equations for zeroth, and rst orders Zeroth Order To zeroth order in, the equations of motion are The solutions can be written by inspection dv x(0) dv y(0) dv z(0) = 0 (67) = 0 (68) = 0 (69) v x(0) = v x0 (70) v y(0) = v y0 (71) v z(0) = v z0 (7) where the v x0, v y0, v z0 are the constant initial velocities. The corresponding zeroth order transverse positions are x (0) = x 0 + v x0 v z0 z (73) y (0) = y 0 + v y0 v z0 z (74) 9
10 where the x 0 and y 0 have been previously de ned to be the initial positions. At this point, we make a simplifying assumption. We assume that all particles have zero initial transverse velocity. In this case and 3.. First Order To rst order in, the equations of motion are dv x(1) dv y(1) dv z(1) With zero initial transverse velocity, these equations become = = = v x(0) = 0 (75) v y(0) = 0 (76) v z(0) = v z0 (77) x (0) = x 0 (78) y (0) = y 0 (79) c B 0 v z(0) fv y(0) B z j 0 v z(0) B y j 0 g (80) c B 0 v z(0) fv z(0) B x j 0 v x(0) B z j 0 g (81) c B 0 v z(0) fv x(0) B y j 0 v y(0) B x j 0 g (8) dv x(1) dv y(1) dv z(1) = c B 0 B y j 0 (83) = c B 0 B x j 0 (84) = 0 (85) Inserting the expressions for the elds in the undulator coordinate system, we have dv x(1) = c B 0 [ B xc sin( y 0 ) sinh ( x 0 ) cos ( z) B xs sin( y 0 ) cosh ( x 0 ) sin ( z) +B yc cos( x 0 ) cosh ( y 0 ) cos ( z + ) +B ys cos( x 0 ) sinh ( y 0 ) sin ( z + )] (86) dv y(1) = c B 0 [B xc cos( y 0 ) cosh ( x 0 ) cos ( z) +B xs cos( y 0 ) sinh ( x 0 ) sin ( z) B yc sin( x 0 ) sinh ( y 0 ) cos ( z + ) B ys sin( x 0 ) cosh ( y 0 ) sin ( z + ) (87) 10
11 dv z(1) = 0 (88) We integrate to nd the velocity components to rst order: v x(1) = c B 0 [ B xc sin( y 0 ) sinh ( x 0 ) sin ( z) +B xs sin( y 0 ) cosh ( x 0 ) cos ( z) +B yc cos( x 0 ) cosh ( y 0 ) sin ( z + ) B ys cos( x 0 ) sinh ( y 0 ) cos ( z + )] (89) v y(1) = c B 0 [B xc cos( y 0 ) cosh ( x 0 ) sin ( z) B xs cos( y 0 ) sinh ( x 0 ) cos ( z) B yc sin( x 0 ) sinh ( y 0 ) sin ( z + ) +B ys sin( x 0 ) cosh ( y 0 ) cos ( z + ) (90) v z(1) = 0 (91) 3.3 Solution To The Equations Of Motion The velocity of a particle in the undulator can now be given to rst order. Using v x = v x(0) + v x(1) + (9) v y = v y(0) + v y(1) + (93) v z = v z(0) + v z(1) + (94) with we get = qb 0 m c (95) v x = q m [ B xc sin( y 0 ) sinh ( x 0 ) sin ( z) +B xs sin( y 0 ) cosh ( x 0 ) cos ( z) +B yc cos( x 0 ) cosh ( y 0 ) sin ( z + ) B ys cos( x 0 ) sinh ( y 0 ) cos ( z + )] (96) 11
12 v y = q m [B xc cos( y 0 ) cosh ( x 0 ) sin ( z) B xs cos( y 0 ) sinh ( x 0 ) cos ( z) B yc sin( x 0 ) sinh ( y 0 ) sin ( z + ) +B ys sin( x 0 ) cosh ( y 0 ) cos ( z + )] (97) In order to simplify future equations, let v z = v z0 (98) K 0 = qb 0 m c (99) Then c v x = K 0 [ cos sin( y 0 ) sinh ( x 0 ) sin ( z) + k s sin sin( y 0 ) cosh ( x 0 ) cos ( z) + k r 4 cos cos( x 0 ) cosh ( y 0 ) sin ( z + ) sin 4 cos( x 0 ) sinh ( y 0 ) cos ( z + )] (100) 4 K Value c v y = K 0 [ k r cos cos( y 0 ) cosh ( x 0 ) sin ( z) sin cos( y 0 ) sinh ( x 0 ) cos ( z) 4 cos sin( x 0 ) sinh ( y 0 ) sin ( z + ) + sin 4 sin( x 0 ) cosh ( y 0 ) cos ( z + )] (101) 4.1 Slippage Slippage is the distance between the light wave and the electrons. The slippage is given by ds = (c v z )dt (10) = ( c v z 1) (103) We take the electron energy to be constant, so the Lorentz factor is constant. 1 = 1 v x c v y c v z c (104) 1
13 So vz c = 1 1 Taking the square root, we nd for large and small v x =c and small v y =c v x c v y c (105) v z c = 1 1 v x c v y c (106) Taking the inverse, we nd So the slippage changes with z according to 4. De ne K c v z = v x c + v y c (107) ds = 1 + v x c + v y c (108) The slippage changes linearly with z with additional oscillatory terms. K value. ds = K So ds K = = c v x + v y We know c v x and c v y. We need to square them and nd the values averaged in z. 4.3 Calculate D E v c x From equation 100, the formula for c v x is c v x = K 0 [ 1 cos sin( y 0 ) sinh ( x 0 ) sin ( z) + k s sin sin( y 0 ) cosh ( x 0 ) cos ( z) + k r 4 cos cos( x 0 ) cosh ( y 0 ) sin ( z + ) sin 4 The linear term de nes the (109) (110) (111) cos( x 0 ) sinh ( y 0 ) cos ( z + )] (11) In order to simplify the calculations, we make the following de nitions for c 1 to c 4. c v x = c 1 sin ( z) + c cos ( z) + c 3 sin ( z + ) + c 4 cos ( z + ) (1)
14 Explicitly, we have Squaring c v x, we have c 1 = K 0 cos sin( y 0 ) sinh ( x 0 ) (114) c = K 0 sin sin( y 0 ) cosh ( x 0 ) (115) 4 c 3 = K 0 cos cos( x 0 ) cosh ( y 0 ) (116) c 4 = K 0 sin 4 cos( x 0 ) sinh ( y 0 ) (117) c v x = c 1 sin ( z) + c cos ( z) + c 3 sin ( z + ) + c 4 cos ( z + ) +c 1 c sin ( z) cos ( z) +c 1 c 3 sin ( z) sin ( z + ) +c 1 c 4 sin ( z) cos ( z + ) +c c 3 cos ( z) sin ( z + ) +c c 4 cos ( z) cos ( z + ) +c 3 c 4 sin ( z + ) cos ( z + ) (118) Taking the average value along z, we have 4.4 Calculate D E v c y v c x From equation 101, the formula for c v y is c v y = K 0 [ k r cos = 1 c 1 + c + c 3 + c c 1 c 3 cos() c 1 c 4 sin () +c c 3 sin () +c c 4 cos () +0 (119) cos( y 0 ) cosh ( x 0 ) sin ( z) sin cos( y 0 ) sinh ( x 0 ) cos ( z) 4 cos sin( x 0 ) sinh ( y 0 ) sin ( z + ) + sin 4 We make the following de nitions for d 1 to d 4. sin( x 0 ) cosh ( y 0 ) cos ( z + )] (10) 14
15 Explicitly, we have Squaring c v y, we have v y = d 1 sin ( z) + d cos ( z) + d 3 sin ( z + ) + d 4 cos ( z + ) d 1 = K 0 cos cos( y 0 ) cosh ( x 0 ) (11) d = K 0 sin cos( y 0 ) sinh ( x 0 ) (1) 4 d 3 = K 0 cos sin( x 0 ) sinh ( y 0 ) () d 4 = K 0 sin 4 sin( x 0 ) cosh ( y 0 ) (14) c v y = d 1 sin ( z) + d cos ( z) + d 3 sin ( z + ) + d 4 cos ( z + ) +d 1 d sin ( z) cos ( z) +d 1 d 3 sin ( z) sin ( z + ) +d 1 d 4 sin ( z) cos ( z + ) +d d 3 cos ( z) sin ( z + ) +d d cos ( z) cos ( z + ) +d 3 d 4 sin ( z + ) cos ( z + ) (15) Taking the average value along z, we have v c y = 1 d 1 + d + d 3 + d d 1 d 3 cos() d 1 d 4 sin () +d d 3 sin () +d d 4 cos () +0 (16) 4.5 Calculate K General Expression Since K = c v x + v y we can use equations 119 and 16 to write an expression for K K = c 1 + c + c 3 + c 4 + d 1 + d + d 3 + d 4 +(c 1 c 3 + c c 4 + d 1 d 3 + d d 4 ) cos() +( c 1 c 4 + c c 3 d 1 d 4 + d d 3 ) sin () (17) 15
16 Recalling from equations 114 to 117 and 11 to 14, the c i and d i in this expression for K have values c 1 = K 0 cos sin( y 0 ) sinh ( x 0 ) (18) c = K 0 sin sin( y 0 ) cosh ( x 0 ) (19) 4 c 3 = K 0 cos cos( x 0 ) cosh ( y 0 ) (0) c 4 = K 0 sin 4 cos( x 0 ) sinh ( y 0 ) (1) d 1 = K 0 cos cos( y 0 ) cosh ( x 0 ) () d = K 0 sin cos( y 0 ) sinh ( x 0 ) (3) 4 d 3 = K 0 cos sin( x 0 ) sinh ( y 0 ) (4) d 4 = K 0 sin 4 sin( x 0 ) cosh ( y 0 ) (5) These are the relations we seek. With these relations, we know K as a function of position for all polarization modes and row settings of the undulator Special Case 1, = 0 Consider the special case of two crossed planar undulators without eld rollo. In this case = 0 and =. The c and d parameters in equation 17 then become c 1 = 0 (6) c = 0 (7) 4 c 3 = K 0 cos cosh ( y 0 ) (8) c 4 = K 0 sin 4 Inserting these values in the expression for K we nd sinh ( y 0 ) (9) d 1 = K 0 cos cosh ( x 0 ) (140) d = K 0 sin sinh ( x 0 ) (141) d 3 = 0 (14) d 4 = 0 (143) K = c 3 + c 4 + d 1 + d (144) 16
17 K = K0[cos 4 + cos cosh ( y 0 ) + sin 4 sinh ( y 0 ) cosh ( x 0 ) + sin sinh ( x 0 )] (145) Note that making the eld rollo small greatly simpli es the formula for K. The position dependence is second order near the magnetic center for all quadrant row positions. The polarization mode does not enter the equation. It is very desirable to shape the magnets so that is small Special Case, x 0 = 0; y 0 = 0 Consider the special case of being on the magnetic axis where x 0 = 0 and y 0 = 0. equations 18 to 5 become Inserting these values in equation 17 for K, we nd K = c 3 + d 1 = K 0 In this case, c 1 = 0 (146) c = 0 (147) 4 c 3 = K 0 cos (148) c 4 = 0 (149) d 1 = K 0 cos (150) d = 0 (151) d 3 = 0 (15) d 4 = 0 (153) kr cos + cos 4 (154) This equation shows that regardless of the value of, on the magnetic axis, the undulator has the ideal behavior where the K value is set by shifting the rows of the two crossed undulators. This is true for all polarization modes Special Case 3, = u, 4 = u This is the case when the undulator is "turned o ", that is, on the magnetic axis K = 0. case, In this = 4 = At general positions we have using equations 18 to 5 (155) (156) c 1 = 0 (157) c = K 0 sin( y 0 ) cosh ( x 0 ) (158) c 3 = 0 (159) c 4 = K 0 cos( x 0 ) sinh ( y 0 ) (160) 17
18 Inserting these values in equation 17 for K, we have d 1 = 0 (161) d = K 0 cos( y 0 ) sinh ( x 0 ) (16) d 3 = 0 (163) d 4 = K 0 sin( x 0 ) cosh ( y 0 ) (164) K = c + c 4 + d + d 4 +(c c 4 + d d 4 ) cos() (165) K = K 0f + + kr " ks sin ( y 0 ) cosh ( x 0 ) + ku ku kr cos ( x 0 ) sinh ( y 0 ) cos ( y 0 ) sinh ks ( x 0 ) + sin ( x 0 ) cosh ( y 0 ) # sin( y 0 ) cosh ( x 0 ) cos( x 0 ) sinh ( y 0 ) cos()g (166) cos( y 0 ) sinh ( x 0 ) sin( x 0 ) cosh ( y 0 ) Note that away from the magnetic axis, the undulator does not "turn o ". value. For the special case of = 0, It still has a nite K The undulator does not have K = 0 o axis even if = 0. K = K0 sinh ( y 0 ) + sinh ( x 0 ) (167) 5 K Values Near The Magnetic Center Typically the undulator is aligned so that the beam is close to the magnetic center. In this case we expand the c and d parameters of equations 18 to 5 to second order in x 0, x 0, y 0 and y 0. We get the following results. k c 1 = K s 0 x 0 y 0 cos (168) k s c = K 0 y 0 sin (169) 1 c 3 = K 0 1 x y0 4 cos (170) k r 4 c 4 = K 0 y 0 sin (171) 18
19 d 1 = K 0 1 d = K 0 k r x 0 sin k d 3 = K sk r 0 x 0 y 0 cos d 4 = K 0 k s x 0 sin 1 y0 + 1 x cos We can now nd K to second order in x 0, x 0, y 0 and y 0. Equation 17 gives K = c 1 + c + c 3 + c 4 + d 1 + d + d 3 + d 4 +(c 1 c 3 + c c 4 + d 1 d 3 + d d 4 ) cos() (17) (173) (174) (175) +( c 1 c 4 + c c 3 d 1 d 4 + d d 3 ) sin () (176) Inserting the second order position dependence, we get K k s 1 = K 0 y 0 sin + K 0 1 x y0 4 cos k r K 0 y 0 sin + K 0 1 y0 + 1 x 0 cos k r k s 4 + K 0 x 0 sin + K 0 x 0 sin 8 h k K i h i 9 s kr 0 x 0 y 0 cos k K r 0 cos k 4 u h >< k + K i h s 0 k + u y 0 sin k k u K i r 0 y 0 sin k 4 u >= h k + K r i h 0 cos k k u K i cos() s kr 0 x 0 y 0 cos k 4 u h >: k + K r 0 x 0 sin k i h k u K s 0 x 0 sin k 4 i u >; 8 h < k K s 0 k + u y 0 sin k i h K r 0 cos k 4 i 9 u h : k K r 0 cos k i h k u K s 0 x 0 sin k 4 i = sin () (177) u ; We proceed to simplify this expression. K = K 0 ku ksy 4 0 sin + K 0 ku kry 4 0 sin 4 + K 0 ku 8 >< + >: 8 < + : k 4 rx 0 sin K 0 ku K 0 k u K 0 ku K 0 k u K 0 ku K 0 ku + K 0 ku kr 1 ksx 0 + kry 0 cos 4 + K 0 ku kr 1 ksy 0 + krx 0 + K 0 ku ksx 4 0 sin 4 cos ku 4 9 ksk rx 0 y 0 cos k u ksk ry 0 sin k u sin ku 4 ksk rx 0 y 0 cos k u cos ku 4 ksk rx 0 sin k u sin ku 4 k s y 0 sin k s x 0 cos cos ku 4 sin ku 4 >= cos() cos >; 9 = sin () (178) ; 19
20 Ordering the terms in powers of x 0, x 0, y 0 and y 0, we nd K = K 0 ku f kr cos + cos 4 +ksk y0 sin k u cos ku 4 r sin ku 4 +k 4 sy 0 sin +k 4 ry 0 cos 4 +krk sy 0 cos +krx 4 0 sin k sk r 6 4 x 0 cos k sin () u + x 0 cos 4 + k 4ry 0 sin 4 + k 4rx 0 cos + k 4sx 0 sin 4 x 0 y 0 cos k u cos ku 4 +y0 sin k u sin ku 4 +x 0 y 0 cos k u cos ku 4 +x 0 sin k u sin ku 4 Rearranging the second order terms and simplifying further, we nd K = K 0 ku f kr cos + cos 4 +ksk y0 sin k u cos ku 4 r x 0 cos k u sin ku 4 + ksy krx 4 0 sin + ksx kry 4 0 sin k sk r sin () + kr k rx 0 + ksy kr k sx 0 + kry 0 cos cos ku 4 x 0 y 0 cos k u + x 0 + y0 sin ku sin ku cos()g (179) cos 4 cos()g (180) This is our desired result. Note that when 6= 0 and 6=, there is a rst order gradient in the K value. The gradient of K near the magnetic center is problematic for aligning the undulator. To rst order in x 0, x 0, y 0 and y 0, K can be expressed as s K = K 0 kr cos k + cos u 4 " 1 + k r cos k s + cos 4 y0 sin k u cos ku 4 x 0 cos k u sin ku 4 # sin () (181) 0
21 The gradient in K at the magnetic center is 0 = 0 = K 0 q cos k s q cos + cos 4 k s + cos 4 We will calculate the size of these terms below. 6 Numerical Calculations 6.1 Parameter Values cos In the calculations below, the following parameter values are used: 4 sin sin () (18) sin 4 cos sin () (183) u = :03 m (184) = u = 196 1/m (185) = 70 1/m (186) = 186 1/m (187) B 0 = :61 T (188) K 0 = 1:8 (189) The value for was measured in the Delta undulator 3. The value for came from the constraint kr = ks + ku. The peak quadrant eld B 0 was calculated from the peak measured eld of 1: T in linear polarization vertical eld mode. The peak eld in this mode B peak is given by 4 0Q = p k and B 0 is given by 0Q, so B 0 = B u peak pkr. The value for K 0 was calculated from B 0 using equation Primary Modes From equation 180 we see that the linear terms near the magnetic center are largest when = = and they are smallest when = 0 or. We now consider these cases. = 0 gives the linear polarization vertical magnetic eld mode, = gives the linear polarization horizontal magnetic eld mode, = = gives the circular polarization right hand magnetic eld mode, and = = gives the circular polarization left hand magnetic eld mode. In these primary modes = 4 and the K value is set by the and 4 row shifts Linear Polarization Vertical Magnetic Field Figure 3 shows how the K value changes with x 0, the horizontal position in the undulator frame, when y 0 = 0. For this plot, the rows were shifted so that = 4 = :008 m, and = 0. This is the linear polarization vertical eld mode. On the magnetic axis, K is given by s kr K = K 0 cos k + cos u 4 (190) 3 Yurii Levashov, private communication. 4 Z. Wolf, H.-D. Nuhn, "Setting The K Value And Polarization Mode Of The Delta Undulator", LCLS-TN-14-, September,
22 The maximum value is obtained when cos k u = 1 and cos ku 4 = 1, and is K max = K 0 kr p (191) When = 4 = :008 m, cos p p = 1= and cos ku 4 = 1=, so kr K = K 0 (19) and K K max = 1 p (193) Figure 3 has several curves. The blue curve is calculated from equation 17. The green dashed curve is an independent calculation starting with the expressions for the magnetic eld components and numerically integrating to get the transverse velocities and from them to calculate K. This provides a check on the calculations in this note. The red curve is from equation 180. It gives an approximate expression for K near the magnetic center. The black line is a linear t to the curves at the magnetic center. It uses equation 190 to nd the K value at the magnetic center and equation 18 to nd the slope. The values of K at the magnetic center, K(0; 0), and the slope at the magnetic center, dk=dx(0; 0), are shown in the gure. The red curve deviates from the blue curve away from x 0 = 0, y 0 = 0. We expect the second order approximation to be valid for x 0, x 0, y 0 and y 0 much less than 1. Since <, we expect the approximation to be valid for or x 0 1 (194) x 0 5 mm (195) For x 0 = 0:5 mm, x 0 = 0:1, ( x 0 ) = 0:01, and the missing third order terms should contribute about K 0:001 = :003 to K. This is roughly in line with the deviations of the red line that we see. We can zoom in on the region near x 0 = 0 in order to estimate the alignment tolerance of the undulator. This is done in gure 4. From the blue curve, K changes by about 10 4 relative to its value at the magnetic center when x 0 is approximately 70 microns from the magnetic center. This is in line with the hyperbolic cosine position dependence of a planar undulator, and is in fact very close to the value in the LCLS undulators which have a similar period (30 mm). The plot of K vs y 0 at x 0 = 0 is shown in gure 5. The gure is identical to the plot of K vs x 0 at y 0 = 0, as we expect from the symmetry of the undulator. 6.. Linear Polarization Horizontal Magnetic Field Figure 6 shows how the K value changes with horizontal position x 0 when y 0 = 0. For this plot, the rows were shifted such that = 4 = :008 m, and =. This is the linear polarization horizontal eld mode. K changes with position in this mode in the same way as in the linear polarization vertical eld mode, as one would expect from symmetry Circular Polarization Right Hand Magnetic Field Figure 7 shows how the K value changes with horizontal position x 0 when y 0 = 0. For this plot, the rows were shifted such that = 4 = :008 m, and = =. This is the circular polarization right hand eld mode. Note the gradient at the magnetic center.
23 Figure 3: This gure shows how K varies with x 0 when y 0 = 0 in linear polarization vertical magnetic eld mode. Figure 8 shows how K varies near the magnetic center. The slope of the K change with x 0 is very large, 161 1/m. A misalignment of the undulator by 100 microns changes K by 0:016, or 0:6%. This is a very large change for a modest misalignment. Figure 9 shows how K varies with y 0 when x 0 = 0. The sign of the change of K with y 0 is opposite to the sign of the change of K with x 0, but the magnitude is the same, as expected from equation 180 with = Circular Polarization Left Hand Magnetic Field Figure 10 shows how the K value changes with horizontal position x 0 when y 0 = 0. For this plot, the rows were shifted such that = 4 = :008 m, and = =. This is the circular polarization left hand eld mode. Again note the gradient at the magnetic center. The variation of K with position in this mode is very similar to the case of right hand circular polarization, as one would expect from symmetry. 3
24 Figure 4: Relative deviation of K from its value at the magnetic center in a small region near the magnetic center. Figure 5: K variation with y 0 when x 0 = 0 in linear polarization vertical magnetic eld mode. 4
25 Figure 6: This gure shows how K varies with x 0 when y 0 = 0 in linear polarization horizontal magnetic eld mode. Figure 7: This gure shows how K varies with x 0 when y 0 = 0 in circular polarization right hand magnetic eld mode. 5
26 Figure 8: K variation with x 0 when y 0 = 0 in circular polarization right hand magnetic eld mode. It is a zoom view near the magnetic center. Figure 9: This gure shows how K varies with y 0 when x 0 = 0 in circular polarization right hand magnetic eld mode. 6
27 Figure 10: Circular polarization left hand magnetic eld mode K vs x 0 with y 0 = 0. 7
28 7 Laboratory Coordinate System The calculations for the Delta undulator are most easily done in the undulator reference frame. This was used in all the calculations up to this point. For comparing the calculations to measurements, however, it is useful to also work in the laboratory reference frame. The relation between frames is x 0L = 1 p (x 0 y 0 ) (196) y 0L = 1 p (x 0 + y 0 ) (197) or equivalently x 0 = 1 p (x 0L + y 0L ) (198) y 0 = 1 p ( x 0L + y 0L ) (199) The subscript L indicates the laboratory frame, and the unsubscripted quantities refer to the undulator frame. Consider equation 181 giving to rst order in x 0, x 0, y 0 and y 0 the expression for the K value near the magnetic center and reproduced below. " K = K 0 s k r cos + cos 4 k 1 + s y0 sin k u cos ku 4 kr cos k u + cos k 4 u x 0 cos k u sin ku 4 Expressing this formula in laboratory coordinates, we nd " 1 + k r cos K = K 0 k s which can be expanded to give s k r + cos 4 K = K 0 s k r cos + cos 4 " p 1 ( x 0L + y 0L ) sin k u cos ku 4 1 p (x 0L + y 0L ) cos cos + cos 4 sin ku 4 # sin () # # sin () (00) x 0L f K 0 q k r sin ksk r cos k u + cos k 4 u cos 4 + cos 4 sin sin ()g +y 0L f K 0 q k r sin ksk r cos k u + cos k 4 u cos cos sin sin ()g (01)
29 The gradients in K at the magnetic center in the laboratory 0L = K 0 sin k s q cos k u + cos k 4 u 4 cos + cos 4 sin sin 0L = K 0 sin k s q cos k u + cos k 4 u 4 cos cos In the primary modes, = 4 =. In 0L = K 0 k s sin 4 sin sin () (03) sin 0L = 0 (05) In the linear modes, sin () = 0 and there are no gradients. In the circular modes, however, sin () = 1 for right and left handed magnetic elds, respectively. In this case for the circular = K 0 k s 0L (06) Inserting the parameter values used previously, namely 0L = 0 (07) = 196 1/m (08) = 186 1/m (09) K 0 = 1:8 (10) sin = p 0L = 7 0L = 0 () where the plus sign is for right handed magnetic elds and the minus sign is for left handed magnetic elds. This is a very large gradient in x 0L. A misalignment of 100 microns in x L causes K to change by :07 for this setting. Note that smaller values of give a smaller gradient. A misalignment in y L does not produce a change in K. 9
30 8 Conclusion The general expressions for the position dependence of K were derived in this note. We found that it is highly desirable to have as small as possible. With > 0, there are rst order gradients in K with position, and the K value will be di cult to set by conventional magnetic measurement and alignment techniques. Acknowledgements I am grateful to Heinz-Dieter Nuhn for many discussions about this work. 30
E ect Of Quadrant Taper On Delta Undulator Phase Errors
LCLS-TN-5- E ect Of Quadrant Taper On Delta Undulator Phase Errors Zachary Wolf SLAC February 8, 05 Abstract The elds in the Delta undulator are very sensitive to the position of the magnet arrays. A taper
More informationCalculation Of Fields On The Beam Axis Of The Delta Undulator
LCLS-TN-3- Calculation Of Fields On The Beam Axis Of The Delta Undulator Zachary Wolf SLAC April, 20 Abstract It is not possible to do Hall probe measurements on the beam axis of the Delta undulator since
More informationCalculation Of Fields On The Beam Axis Of The Delta Undulator
LCLS-TN-3- Calculation Of Fields On The Beam Axis Of The Delta Undulator Zachary Wolf SLAC December 8, 203 Abstract It is not possible to do Hall probe measurements on the beam axis of the Delta undulator
More informationA PPM Phase Shifter Design
LCLS-TN-- A PPM Phase Shifter Design Zachary Wolf SLAC July 30, 0 Abstract This note describes the design of a pure permanent magnet (PPM) phase shifter for LCLS-II. A PPM phase shifter has the advantage
More informationDelta Undulator End Design
LCLS-TN-1- Delta Undulator End Design Zachary Wolf SLAC January 9, 01 Abstract This note describes the magnetic end design of the Delta undulator. Both three block ends and two block ends were studied.
More informationUndulator Field Integral Measurements
LCLS-TN-05-22 Undulator Field Integral Measurements Zachary Wolf SLAC August 5, 2005 Abstract The LCLS undulator eld integrals must be very small so that the beam trajectory slope and o set stay within
More informationUndulator K Value Temperature Dependence
LCLS-TN-09-4 Undulator K Value Temperature Dependence Zachary Wolf, Yurii Levashov SLAC October 6, 2009 Abstract LCLS undulator SN27 was used to measure the temperature dependence of its K value. The results
More informationBeam Pipe Corrector Study
LCLS-TN-12-1 Beam Pipe Corrector Study Zachary Wolf, Yurii Levashov SLAC March 19, 2012 Abstract This note describes the test of a corrector magnet which can be built into the beam pipe of the LCLS-II
More informationIntroduction To LCLS Undulator Tuning
LCLS-TN-04-7 Introduction To LCLS Undulator Tuning Zachary Wolf June 3, 2004 Abstract This note gives a general introduction to undulator tuning for the LCLS. It starts with a theoretical discussion in
More informationA Magnetic Measurement Plan For The Delta Undulator
LCLS-TN-3-4 A Magnetic Measurement Plan For The Delta Undulator Zachary Wolf SLAC March 28, 23 Abstract This note describes a magnetic measurement plan for the Delta undulator. The note starts by discussing
More informationE ect Of SXR Beam Pipe Temperature On Undulator Magnet Temperature
LCLS-TN-18-7 E ect Of SXR Beam Pipe Temperature On Undulator Magnet Temperature Zachary Wolf SLAC September 10, 2018 Abstract This note documents a test that was made with an LCLS-II SXR undulator with
More informationE ect Of Quadrant Bow On Delta Undulator Phase Errors
LCLS-TN-15-1 E ect Of Qadrant Bow On Delta Undlator Phase Errors Zachary Wolf SLAC Febrary 18, 015 Abstract The Delta ndlator qadrants are tned individally and are then assembled to make the tned ndlator.
More information4.3 - Linear Combinations and Independence of Vectors
- Linear Combinations and Independence of Vectors De nitions, Theorems, and Examples De nition 1 A vector v in a vector space V is called a linear combination of the vectors u 1, u,,u k in V if v can be
More informationSetting The K Value And Polarization Mode Of The Delta Undulator
LCLS-TN-4- Setting The Vale And Polarization Mode Of The Delta Undlator Zachary Wolf, Heinz-Dieter Nhn SLAC September 4, 04 Abstract This note provides the details for setting the longitdinal positions
More informationQuadrupole Induced Resonant Particle Transport in a Pure Electron Plasma
Quadrupole Induced Resonant Particle Transport in a Pure Electron Plasma E. Gilson 1 and J. Fajans 2 Department of Physics University of California, Berkeley Berkeley, California, 94720-7300 Abstract.
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 information2 Formulation. = arg = 2 (1)
Acoustic di raction by an impedance wedge Aladin H. Kamel (alaahassan.kamel@yahoo.com) PO Box 433 Heliopolis Center 11757, Cairo, Egypt Abstract. We consider the boundary-value problem for the Helmholtz
More informationWakefield Effects of Collimators in LCLS-II
Wakefield Effects of Collimators in LCLS-II LCLS-II TN-15-38 10/7/2015 K. Bane SLAC, Menlo Park, CA 94025, USA October 7, 2015 LCLSII-TN-15-38 L C L S - I I T E C H N I C A L N O T E LCLS-TN-15-38 October
More informationparameter symbol value beam energy E 15 GeV transverse rms beam size x;y 25 m rms bunch length z 20 m charge per bunch Q b 1nC electrons per bunch N b
LCLS{TN{98{2 March 1998 SLAC/AP{109 November 1997 Ion Eects in the LCLS Undulator 1 Frank Zimmermann Stanford Linear Accelerator Center, Stanford University, Stanford, CA 94309 Abstract I calculate the
More informationfree- electron lasers I (FEL)
free- electron lasers I (FEL) produc'on of copious radia'on using charged par'cle beam relies on coherence coherence is achieved by insuring the electron bunch length is that desired radia'on wavelength
More informationThree Loose Ends: Edge Focusing; Chromaticity; Beam Rigidity.
Linear Dynamics, Lecture 5 Three Loose Ends: Edge Focusing; Chromaticity; Beam Rigidity. Andy Wolski University of Liverpool, and the Cockcroft Institute, Daresbury, UK. November, 2012 What we Learned
More informationLecture 3: Modeling Accelerators Fringe fields and Insertion devices. X. Huang USPAS, January 2015 Hampton, Virginia
Lecture 3: Modeling Accelerators Fringe fields and Insertion devices X. Huang USPAS, January 05 Hampton, Virginia Fringe field effects Dipole Quadrupole Outline Modeling of insertion devices Radiation
More informationSingle Particle Motion
Single Particle Motion C ontents Uniform E and B E = - guiding centers Definition of guiding center E gravitation Non Uniform B 'grad B' drift, B B Curvature drift Grad -B drift, B B invariance of µ. Magnetic
More informationTraveling Wave Undulators for FELs and Synchrotron Radiation Sources
LCLS-TN-05-8 Traveling Wave Undulators for FELs and Synchrotron Radiation Sources 1. Introduction C. Pellegrini, Department of Physics and Astronomy, UCLA 1 February 4, 2005 We study the use of a traveling
More informationSome Explicit Solutions of the Cable Equation
Some Explicit Solutions of the Cable Equation Marco Herrera-Valdéz and Sergei K. Suslov Mathematical, Computational and Modeling Sciences Center, Arizona State University, Tempe, AZ 85287 1904, U.S.A.
More informationACCURATE FREE VIBRATION ANALYSIS OF POINT SUPPORTED MINDLIN PLATES BY THE SUPERPOSITION METHOD
Journal of Sound and Vibration (1999) 219(2), 265 277 Article No. jsvi.1998.1874, available online at http://www.idealibrary.com.on ACCURATE FREE VIBRATION ANALYSIS OF POINT SUPPORTED MINDLIN PLATES BY
More information94 CHAPTER 3. VECTORS AND THE GEOMETRY OF SPACE
94 CHAPTER 3. VECTORS AND THE GEOMETRY OF SPACE 3.3 Dot Product We haven t yet de ned a multiplication between vectors. It turns out there are di erent ways this can be done. In this section, we present
More informationWarped Products. by Peter Petersen. We shall de ne as few concepts as possible. A tangent vector always has the local coordinate expansion
Warped Products by Peter Petersen De nitions We shall de ne as few concepts as possible. A tangent vector always has the local coordinate expansion a function the di erential v = dx i (v) df = f dxi We
More informationRADIATION OF ELECTROMAGNETIC WAVES
Chapter RADIATION OF ELECTROMAGNETIC WAVES. Introduction In the preceding Chapters, some general properties of TEM waves have been discussed. A question arises as to how to excite electromagnetic waves.
More informationAccelerator Physics Homework #3 P470 (Problems: 1-5)
Accelerator Physics Homework #3 P470 (Problems: -5). Particle motion in the presence of magnetic field errors is (Sect. II.2) y + K(s)y = B Bρ, where y stands for either x or z. Here B = B z for x motion,
More information23.1 Chapter 8 Two-Body Central Force Problem (con)
23 Lecture 11-20 23.1 Chapter 8 Two-Body Central Force Problem (con) 23.1.1 Changes of Orbit Before we leave our discussion of orbits we shall discuss how to change from one orbit to another. Consider
More informationX-Band RF Harmonic Compensation for Linear Bunch Compression in the LCLS
SLAC-TN-5- LCLS-TN-1-1 November 1,1 X-Band RF Harmonic Compensation for Linear Bunch Compression in the LCLS Paul Emma SLAC November 1, 1 ABSTRACT An X-band th harmonic RF section is used to linearize
More informationPhysics 6303 Lecture 15 October 10, Reminder of general solution in 3-dimensional cylindrical coordinates. sinh. sin
Physics 6303 Lecture 15 October 10, 2018 LAST TIME: Spherical harmonics and Bessel functions Reminder of general solution in 3-dimensional cylindrical coordinates,, sin sinh cos cosh, sin sin cos cos,
More informationTaylor series - Solutions
Taylor series - Solutions. f(x) sin(x) sin(0) + x cos(0) + x x ( sin(0)) +!! ( cos(0)) + + x4 x5 (sin(0)) + 4! 5! 0 + x + 0 x x! + x5 5! x! + 0 + x5 (cos(0)) + x6 6! ( sin(0)) + x 7 7! + x9 9! 5! + 0 +
More informationLCLS Undulators Present Status and Future Upgrades
LCLS Undulators Present Status and Future Upgrades Heinz-Dieter Nuhn LCLS Undulator Group Leader 1 1 Heinz-Dieter Nuhn Linac Coherent Light Source INJECTOR LINAC BEAM TRANSPORT UNDULATOR HALL 2 2 Heinz-Dieter
More informationStudent Solutions Manual. for. Web Sections. Elementary Linear Algebra. 4 th Edition. Stephen Andrilli. David Hecker
Student Solutions Manual for Web Sections Elementary Linear Algebra th Edition Stephen Andrilli David Hecker Copyright, Elsevier Inc. All rights reserved. Table of Contents Lines and Planes and the Cross
More informationElectromagnetic Radiation from Relativistic Electron Beams
Electromagnetic Radiation from Relativistic Electron Beams Helmut Wiedemann Stanford University April 2000 Abstract These notes represent an introduction to the physics of synchrotron radiation. Such radiation
More informationPart II Effect of Insertion Devices on the Electron Beam
Part II Effect of Insertion Devices on the Electron Beam Pascal ELLEAUME European Synchrotron Radiation Facility, Grenoble II, 1/14, P. Elleaume, CAS, Brunnen July -9, 3. Effect of an Insertion Device
More informationDynamic Green Function Solution of Beams Under a Moving Load with Dierent Boundary Conditions
Transaction B: Mechanical Engineering Vol. 16, No. 3, pp. 273{279 c Sharif University of Technology, June 2009 Research Note Dynamic Green Function Solution of Beams Under a Moving Load with Dierent Boundary
More informationTransverse Dynamics II
Transverse Dynamics II JAI Accelerator Physics Course Michaelmas Term 217 Dr. Suzie Sheehy Royal Society University Research Fellow University of Oxford Acknowledgements These lectures have been produced
More informationLecture notes 5: Diffraction
Lecture notes 5: Diffraction Let us now consider how light reacts to being confined to a given aperture. The resolution of an aperture is restricted due to the wave nature of light: as light passes through
More informationGreen s function for the wave equation
Green s function for the wave equation Non-relativistic case January 2018 1 The wave equations In the Lorentz Gauge, the wave equations for the potentials are (Notes 1 eqns 44 and 43): 1 2 2 2 2 0 (1)
More information1 The Linearized Einstein Equations
The Linearized Einstein Equations. The Assumption.. Simplest Version The simplest version of the linearized theory begins with at Minkowski spacetime with basis vectors x and metric tensor components 8
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 informationSecond-Order Linear ODEs
Chap. 2 Second-Order Linear ODEs Sec. 2.1 Homogeneous Linear ODEs of Second Order On pp. 45-46 we extend concepts defined in Chap. 1, notably solution and homogeneous and nonhomogeneous, to second-order
More informationMath 443/543 Graph Theory Notes 5: Graphs as matrices, spectral graph theory, and PageRank
Math 443/543 Graph Theory Notes 5: Graphs as matrices, spectral graph theory, and PageRank David Glickenstein November 3, 4 Representing graphs as matrices It will sometimes be useful to represent graphs
More information1 Solutions in cylindrical coordinates: Bessel functions
1 Solutions in cylindrical coordinates: Bessel functions 1.1 Bessel functions Bessel functions arise as solutions of potential problems in cylindrical coordinates. Laplace s equation in cylindrical coordinates
More informationLinearized methods for ordinary di erential equations
Applied Mathematics and Computation 104 (1999) 109±19 www.elsevier.nl/locate/amc Linearized methods for ordinary di erential equations J.I. Ramos 1 Departamento de Lenguajes y Ciencias de la Computacion,
More informationLow Emittance Machines
Advanced Accelerator Physics Course RHUL, Egham, UK September 2017 Low Emittance Machines Part 1: Beam Dynamics with Synchrotron Radiation Andy Wolski The Cockcroft Institute, and the University of Liverpool,
More informationType II hidden symmetries of the Laplace equation
Type II hidden symmetries of the Laplace equation Andronikos Paliathanasis Larnaka, 2014 A. Paliathanasis (Univ. of Naples) Type II symmetries Larnaka, 2014 1 / 25 Plan of the Talk 1 The generic symmetry
More informationÆThe Betatron. Works like a tranformer. Primary winding : coils. Secondary winding : beam. Focusing from beveled gap.
Weak Focusing Not to be confused with weak folk cussing. Lawrence originally thought that the cyclotron needed to have a uniform (vertical) field. Actually unstable: protons with p vert 0 would crash into
More informationCompound Damped Pendulum: An Example
Compound Damped Pendulum: An Example Temple H. Fay Department of Mathematical Technology 1 Tshwane University of Technology Pretoria, South Africa thf ay@hotmail:com Abstract: In this article, we use an
More informationResistive wall impedance and tune shift for a chamber with a finite thickness
PHYSICAL REVIEW E 66 05650 00 Resistive wall impedance and tune shift for a chamber with a finite thickness Yoshihiro Shobuda and Kaoru Yokoya KEK High Energy Accelerator Research Organization - Oho Tsukuba
More informationJanuary 13, 2005 CBN 05-1 IDEAL WIGGLER 1. A.Mikhailichenko, Cornell University, LEPP, Ithaca NY 14853
January 1, 005 CBN 05-1 IDEAL WIGGLER 1 A.Mikhailichenko, Cornell University, LEPP, Ithaca NY 1485 Abstract. Described is the wiggler with reduced nonlinear components for usage in the damping ring of
More informationThe A, B, C and D are determined by these 4 BCs to obtain
Solution:. Floquet transformation: (a) Defining a new coordinate η = y/ β and φ = (/ν) s 0 ds/β, we find ds/dφ = νβ, and dη dφ = ds dη dφ d 2 η dφ 2 = ν2 β ( β y ) ( 2 β 3/2 β y = ν β /2 y ) 2 β /2 β y,
More informationa11 a A = : a 21 a 22
Matrices The study of linear systems is facilitated by introducing matrices. Matrix theory provides a convenient language and notation to express many of the ideas concisely, and complicated formulas are
More informationSLAC-TN-93-2 SLAC/SSRL-0002 January 1993 (SSRL/ACD) COMPARISON OF PLANAR HELICAL UNDULATOR DESIGNS FOR SPEAR BEAMLINE FIVE*
SLAC-TN-93-2 SLAC/SSRL-0002 January 1993 (SSRL/ACD) COMPARISON OF PLANAR HELICAL UNDULATOR DESIGNS FOR SPEAR BEAMLINE FIVE* Roger Car-r Stanford Linear Accelerator Center Stanford Synchrotron Radiation
More informationnotes 5d. Details of the Construction of the Smale horseshoe map.
by 3940-07 notes 5d Details of the Construction of the Smale horseshoe map. Let f (x; y) = ( 4 x + 8 ; 5y ). Then f maps D onto the rectangle R de ned 8 x 3 8 y 9 : In particular, it squeezes D horizontally
More informationEffect of Insertion Devices. Effect of IDs on beam dynamics
Effect of Insertion Devices The IDs are normally made of dipole magnets ith alternating dipole fields so that the orbit outside the device is un-altered. A simple planer undulator ith vertical sinusoidal
More informationQuasi-static analysis of asymmetric cylindrical coplanar waveguides with nite-extent ground
INT. J. ELECTRONICS, 2000, VOL. 87, NO. 2, 185± 198 Quasi-static analysis of asymmetric cylindrical coplanar waveguides with nite-extent ground NIHAD DIB*{ and ASEM AL-ZOUBI{ Quasi-static characteristics
More information2 Feynman rules, decay widths and cross sections
2 Feynman rules, decay widths and cross sections 2.1 Feynman rules Normalization In non-relativistic quantum mechanics, wave functions of free particles are normalized so that there is one particle in
More informationFigure 1: Grad, Div, Curl, Laplacian in Cartesian, cylindrical, and spherical coordinates. Here ψ is a scalar function and A is a vector field.
Figure 1: Grad, Div, Curl, Laplacian in Cartesian, cylindrical, and spherical coordinates. Here ψ is a scalar function and A is a vector field. Figure 2: Vector and integral identities. Here ψ is a scalar
More informationThin airfoil theory. Chapter Compressible potential flow The full potential equation
hapter 4 Thin airfoil theory 4. ompressible potential flow 4.. The full potential equation In compressible flow, both the lift and drag of a thin airfoil can be determined to a reasonable level of accuracy
More informationApproximate Solutions of the Grad-Schlüter-Shafranov Equation
Approximate Solutions of the Grad-Schlüter-Shafranov Equation Gerson Otto Ludwig Associated Plasma Laboratory, National Space Research Institute 17-010, São José dos Campos, SP, Brazil ludwig@plasma.inpe.br
More informationLCLS-II Undulator Tolerance Budget*
LCLS-TN-13-5 rev2 April 16, 2013 LCLS-II Undulator Tolerance Budget* Heinz-Dieter Nuhn, Stanford Linear Accelerator Center, MS-18, 2575 Sandhill Road, Menlo Park, California 94025. 1 Introduction This
More informationFree electron lasers
Preparation of the concerned sectors for educational and R&D activities related to the Hungarian ELI project Free electron lasers Lecture 2.: Insertion devices Zoltán Tibai János Hebling 1 Outline Introduction
More informationStochastic solutions of nonlinear pde s: McKean versus superprocesses
Stochastic solutions of nonlinear pde s: McKean versus superprocesses R. Vilela Mendes CMAF - Complexo Interdisciplinar, Universidade de Lisboa (Av. Gama Pinto 2, 1649-3, Lisbon) Instituto de Plasmas e
More informationEquations of linear stellar oscillations
Chapter 4 Equations of linear stellar oscillations In the present chapter the equations governing small oscillations around a spherical equilibrium state are derived. The general equations were presented
More informationProperties of the Gradient
Properties of the Gradient Gradients and Level Curves In this section, we use the gradient and the chain rule to investigate horizontal and vertical slices of a surface of the form z = g (x; y) : To begin
More informationSecond Printing Errata for Advanced Engineering Mathematics, by Dennis G. Zill and Michael R. Cullen
59X_errataSecondPrint.qxd /7/7 9:9 AM Page Second Printing Errata for Advanced Engineering Mathematics, by Dennis G. Zill and Michael R. Cullen (Yellow highlighting indicates corrected material) Page 7
More informationErrata Instructor s Solutions Manual Introduction to Electrodynamics, 3rd ed Author: David Griffiths Date: September 1, 2004
Errata Instructor s Solutions Manual Introduction to Electrodynamics, 3rd ed Author: David Griffiths Date: September, 004 Page 4, Prob..5 (b): last expression should read y +z +3x. Page 4, Prob..6: at
More informationDust devils, water spouts, tornados
Balanced flow Things we know Primitive equations are very comprehensive, but there may be a number of vast simplifications that may be relevant (e.g., geostrophic balance). Seems that there are things
More informationGenerators For Nonlinear Three Term Recurrences Involving The Greatest Integer Function
Applied Mathematics E-Notes, 11(2011), 73-81 c ISSN 1607-2510 Available free at mirror sites of http://www.math.nthu.edu.tw/amen/ Generators For Nonlinear Three Term Recurrences Involving The Greatest
More informationINTRODUCTION TO FINITE ELEMENT METHODS ON ELLIPTIC EQUATIONS LONG CHEN
INTROUCTION TO FINITE ELEMENT METHOS ON ELLIPTIC EQUATIONS LONG CHEN CONTENTS 1. Poisson Equation 1 2. Outline of Topics 3 2.1. Finite ifference Method 3 2.2. Finite Element Method 3 2.3. Finite Volume
More informationOrdinary Differential Equations (ODEs)
c01.tex 8/10/2010 22: 55 Page 1 PART A Ordinary Differential Equations (ODEs) Chap. 1 First-Order ODEs Sec. 1.1 Basic Concepts. Modeling To get a good start into this chapter and this section, quickly
More informationMATH 2203 Final Exam Solutions December 14, 2005 S. F. Ellermeyer Name
MATH 223 Final Exam Solutions ecember 14, 25 S. F. Ellermeyer Name Instructions. Your work on this exam will be graded according to two criteria: mathematical correctness and clarity of presentation. In
More informationFixed Point Theorem and Sequences in One or Two Dimensions
Fied Point Theorem and Sequences in One or Two Dimensions Dr. Wei-Chi Yang Let us consider a recursive sequence of n+ = n + sin n and the initial value can be an real number. Then we would like to ask
More informationECON0702: Mathematical Methods in Economics
ECON0702: Mathematical Methods in Economics Yulei Luo SEF of HKU January 12, 2009 Luo, Y. (SEF of HKU) MME January 12, 2009 1 / 35 Course Outline Economics: The study of the choices people (consumers,
More informationMagnetostatics and the vector potential
Magnetostatics and the vector potential December 8, 2015 1 The divergence of the magnetic field Starting with the general form of the Biot-Savart law, B (x 0 ) we take the divergence of both sides with
More informationHamiltonian Mechanics
Alain J. Brizard Saint Michael's College Hamiltonian Mechanics 1 Hamiltonian The k second-order Euler-Lagrange equations on con guration space q =(q 1 ; :::; q k ) d @ _q j = @q j ; (1) can be written
More informationEstimate of Undulator Magnet Damage Due to Beam Finder Wire Measurements
LCLS-TN-06-6 Estimate of Undulator Magnet Damage Due to Beam Finder Wire Measurements J. Welch April 5, 2006 Beam Finder Wire (BFW) devices will be installed at each break in the Undulator magnet line.
More informationMassachusetts Institute of Technology Physics Department. Physics 8.20 IAP 2005 Special Relativity January 28, 2005 FINAL EXAM
Massachusetts Institute of Technology Physics Department Physics 8.20 IAP 2005 Special Relativity January 28, 2005 FINAL EXAM Instructions You have 2.5 hours for this test. Papers will be picked up promptly
More informationThe Kuhn-Tucker Problem
Natalia Lazzati Mathematics for Economics (Part I) Note 8: Nonlinear Programming - The Kuhn-Tucker Problem Note 8 is based on de la Fuente (2000, Ch. 7) and Simon and Blume (1994, Ch. 18 and 19). The Kuhn-Tucker
More informationMAE210C: Fluid Mechanics III Spring Quarter sgls/mae210c 2013/ Solution II
MAE210C: Fluid Mechanics III Spring Quarter 2013 http://web.eng.ucsd.edu/ sgls/mae210c 2013/ Solution II D 4.1 The equations are exactly the same as before, with the difference that the pressure in the
More informationChapter 5. Sound Waves and Vortices. 5.1 Sound waves
Chapter 5 Sound Waves and Vortices In this chapter we explore a set of characteristic solutions to the uid equations with the goal of familiarizing the reader with typical behaviors in uid dynamics. Sound
More informationMagnetic Counterforce for Insertion Devices*
SLAC-PUB-9594 October 2002 Magnetic Counterforce for Insertion Devices* Roger Carr Stanford Linear Accelerator Center, Stanford University, Stanford, CA 94309 Abstract In a standard insertion device, such
More informationRadiation from a Moving Charge
Radiation from a Moving Charge 1 The Lienard-Weichert radiation field For details on this theory see the accompanying background notes on electromagnetic theory (EM_theory.pdf) Much of the theory in this
More informationValidity domain and limitation of non-retarded Green's tensor for electromagnetic scattering at surfaces
1 October 2 Optics Communications 184 (2) 37±47 www.elsevier.com/locate/optcom Validity domain and limitation of non-retarded Green's tensor for electromagnetic scattering at surfaces Philippe Gay-Balmaz
More informationStructural Dynamics. Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma).
Structural Dynamics Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma). We will now look at free vibrations. Considering the free
More informationEmittance Limitation of a Conditioned Beam in a Strong Focusing FEL Undulator. Abstract
SLAC PUB 11781 March 26 Emittance Limitation of a Conditioned Beam in a Strong Focusing FEL Undulator Z. Huang, G. Stupakov Stanford Linear Accelerator Center, Stanford, CA 9439 S. Reiche University of
More informationStochastic Processes
Stochastic Processes A very simple introduction Péter Medvegyev 2009, January Medvegyev (CEU) Stochastic Processes 2009, January 1 / 54 Summary from measure theory De nition (X, A) is a measurable space
More informationSome Sample Calculations for the Far Field Harmonic Power and Angular Pattern in LCLS-1 and LCLS-2
Some Sample Calculations for the Far Field Harmonic Power and Angular Pattern in LCLS-1 and LCLS-2 W.M. Fawley February 2013 SLAC-PUB-15359 ABSTRACT Calculations with the GINGER FEL simulation code are
More informationElectromagnetic Implosion Using an Array
Sensor and Simulation Notes Note 57 July 2006 Electromagnetic Implosion Using an Array Carl E. Baum University of New Mexico Department of Electrical and Computer Engineering Albuquerque New Mexico 873
More informationMOTT-RUTHERFORD SCATTERING AND BEYOND. Abstract. The electron charge is considered to be distributed or extended in
MOTT-RUTHERFORD SCATTERING AND BEYOND Abstract The electron charge is considered to be distributed or extended in space. The differential of the electron charge is set equal to a function of the electron
More information1. (a) +EA; (b) EA; (c) 0; (d) 0 2. (a) 2; (b) 3; (c) 1 3. (a) equal; (b) equal; (c) equal e; (b) 150e 5. 3 and 4 tie, then 2, 1
CHAPTER 24 GAUSS LAW 659 CHAPTER 24 Answer to Checkpoint Questions 1. (a) +EA; (b) EA; (c) ; (d) 2. (a) 2; (b) 3; (c) 1 3. (a) eual; (b) eual; (c) eual 4. +5e; (b) 15e 5. 3 and 4 tie, then 2, 1 Answer
More informationE. Wilson - CERN. Components of a synchrotron. Dipole Bending Magnet. Magnetic rigidity. Bending Magnet. Weak focusing - gutter. Transverse ellipse
Transverse Dynamics E. Wilson - CERN Components of a synchrotron Dipole Bending Magnet Magnetic rigidity Bending Magnet Weak focusing - gutter Transverse ellipse Fields and force in a quadrupole Strong
More informationCERN Accelerator School. Intermediate Accelerator Physics Course Chios, Greece, September Low Emittance Rings
CERN Accelerator School Intermediate Accelerator Physics Course Chios, Greece, September 2011 Low Emittance Rings Part 1: Beam Dynamics with Synchrotron Radiation Andy Wolski The Cockcroft Institute, and
More informationRings, Integral Domains, and Fields
Rings, Integral Domains, and Fields S. F. Ellermeyer September 26, 2006 Suppose that A is a set of objects endowed with two binary operations called addition (and denoted by + ) and multiplication (denoted
More informationENGIN 211, Engineering Math. Laplace Transforms
ENGIN 211, Engineering Math Laplace Transforms 1 Why Laplace Transform? Laplace transform converts a function in the time domain to its frequency domain. It is a powerful, systematic method in solving
More informationMOTION IN TWO OR THREE DIMENSIONS
MOTION IN TWO OR THREE DIMENSIONS 3 Sections Covered 3.1 : Position & velocity vectors 3.2 : The acceleration vector 3.3 : Projectile motion 3.4 : Motion in a circle 3.5 : Relative velocity 3.1 Position
More information