P(E) = Γ
(E−Er)2+ (Γ/2)2 (2.20)
The resonance energy is given byEr.
It is now possible to plot the resonance strength and the width at half maximum ∆ of the Gamow window at different temperatures side by side with the nuclear core level scheme of the product nuclei.
2.4 Ion Optics
The general task in beam optics is to transport charged particles from point A to point B along a desired path. The collection of bending and focusing magnets installed along this ideal path are called the magnet lattice and the complete optical system including the bending and focusing parameters is called a beam transport system (see also [28]).
2.4.1 Particle Beam Guidance
To guide a charged particle along a predefined path, magnetic fields are used which deflect particles as determined by the equilibrium of the centrifugal force and Lorentz force
mγv2~k+q
~ v×B~
= 0 (2.21)
where~k= (kx, ky,0) is the local curvature vector of the trajectory, which is pointing in the direction of the centrifugal force. Assuming that the magnetic field vectorB~ is oriented normal to the velocity vector~v, the treatment of linear beam dynamics is restricted to purely transverse fields. This has no fundamental reason other than to simplify the formulation of particle beam dynamics. The transverse components of the particle velocities for relativistic beams are small compared to the particle velocity.
A curvilinear coordinate system(x, y, z)following the ideal path is used. The direction of the particle is denoted with the coordinates: The bending radius for the particle trajectory in a magnetic field is from equation2.21withp=γmv
~kx,y=±qc βE
B~y,x (2.22)
and the angular frequency of revolution of a particle on a complete orbit normal to the fieldB is ωL=
qc EB
(2.23)
which is also called the cyclotron or Larmor frequency. The sign in equation2.22has been chosen to meet the definition of curvature in analytical geometry, where the curvature is negative, if the tangent to the trajectory rotates counterclockwise. Often, the beam rigidity, defined as
|Bρ|= p0
q (2.24)
is used to normalize the magnet strength. Using more practical units the expressions for the beam rigidity and bending radius become
Bρ[Tm] = 10
2.998βE[GeV] (2.25)
Physics Department 11 Technische Universität München
12 CHAPTER 2. THEORY
and
1
ρ[m−1] = B
Bρ= 0.2998Z A
|B[T]|
βE[GeV/u] (2.26)
where the sign for the bending radius was dropped. For relativistic particles this expression is further simplified sinceβ ≈1. The deflection angle in a magnetic field is
θ= Z dz
ρ (2.27)
or for a uniform field like in a dipole magnet of arc length lm the deflection angle isθ=lm/ρ.
2.4.2 Particle Beam Focusing
Similar to the properties of light rays, particle beams also have a tendency to spread out due to an inherent beam divergence. To keep the particle beam together and to generate specifically desired beam properties at selected points along the beam transport line, focusing devices are required. Any magnetic field, that deflects a particle by an angle proportional to its distance r from the axis of the focusing device, will act in the same way as a glass lens does in the approximation of paraxial, geometric optics for visible light. Iff is the focal length, the deflection angleαis defined by
α=−r
f (2.28)
A similar focusing property can be provided for charged particle beams by the use of azimuthal magnetic fieldsBφ with the property
α=−l
ρ =−qc
βEBφl=−qc
βEgrl (2.29)
wherelis the path length of the particle trajectory in the magnetic fieldBφandgis the field gradient defined by Bφ=gr or byg= dBφ/dr.
In beam dynamics, it is customary to define an energy independent focusing strength. Similar to the definition of the bending curvature a focusing strengthK is define by
K= q pg= qc
βEg (2.30)
2.4.3 Equation of Motion
Magnetic fields are used to guide charged particles along a prescribed path or at least keep them close by. This path, or reference trajectory, is defined geometrically by straight sections and bending magnets only. Dipole magnets deflect the path and quadrupole and higher order magnets do not influence this path but provide the focusing forces necessary to keep all particles close to the reference path.
The most convenient coordinate system to describe particle motion is the curvilinear coordinate system as seen in Figure 2.3. The curvatures are functions of the coordinate z and are nonzero only where there are bending magnets. The equations of motion are derived in the horizontal plane only. The generalization to both horizontal and vertical plane is straightforward. Using the notation of Figure 2.4 the deflection angle of the ideal path is dφ0 = dz/ρ0 or utilizing the curvature to preserve the directionality of the deflection
dφ0= +k0dz (2.31)
Technische Universität München 12 Physics Department
2.4. ION OPTICS 13
y Fig. 2.3 Locally defined right-handed
coordinate system used with bending design.
Here s is the distance along the design orbit, x is the distance from this orbit along the radius of curvature, and y is the distance from the design orbit out of the bend plane.
along the trajectory are a superposition of angle and offset perturbations. From these observations, we may strongly suspect that the perturbed trajectory is a harmonic oscillation, but we must first develop the proper analysis tools to verify this suspicion.
As stated in Chapter 1, as we will typically take the distance along the design orbit to be the independent variable (in this case indicated by
s), we implicitlywish to analyze the charged particle dynamics near the design orbit. In the present case, this orbit is specified by a certain radius of curvature
R(and thus a certain momentum
p0 = qB0R), and center of curvature,(x0,
y0). With this choice of analysis geometry, we can locally define a new right-handed coordinate system
(x,y,s), as shown in Fig. 2.3. In this coordinate system, xis the distance of the orbit under consideration from the design orbit, in the direction measured along the radius and normal to
s. The distancey(formerly indicated by the coordinate
zin Section 2.1) is measured from the design orbit
to the particle orbit under consideration, in the direction out of the bend plane.
1 1This convention, in which the symboly is defined to be the distance out of the bend plane is typical of the American literature. European beam physicists more often use the symbolz instead, but we do not follow this conven-tion even though it connects more naturally to our previous discussion. This is because our adopted convention makes subsequent derivations somewhat easier to understand, and also because it allows the connection between linear accelerator and circular accel-erator coordinate systems to become more obvious.The choice of a right-handed system in this case is a function of the direction of the bend, and in simple circular accelerators, one is free to construct the curvilinear coordinates once and for all. On the other hand, when we encounter bends in the opposing direction, as in chicane systems (see Chapter 3), we will choose to consistently define the coordinate
x, so that it is positive alongthe direction away from the origin of the bend. As we will also choose to leave the vertical direction unchanged in this transformation, a left-handed coordinate system will result when the bend direction is changed.
The coordinate system shown in Fig. 2.3 is quite similar to a cylindrical coordinate system, with
xrelated to the radial variable
ρby the definition
x≡ρ−R,sreplacing the azimuthal angle
φ (ds
=Rd
φ), and
y, as previouslynoted, replacing
z. Thus we can write the equations of motion for orbits in thissystem by using the Lagrange–Euler formulation (see Problem 2.1), as
dp
ρdt
= γm0v2φρ −qvφB0
, (2.9)
where
vφ =ρφ˙is the azimuthal velocity.
Equation (2.9) can be cast as a familiar differential equation by using
xas a small variable (x
R, which is also equivalent, as will be seen below, tothe paraxial ray approximation) to linearize the relation. This is accomplished through use of a lowest order Taylor series expansion of the motion about the design orbit equilibrium
(px =pρ =0) at
ρ =R,dp
xdt
∼= −γ0m0v20R2 x.
(2.10)
Figure 2.3: Locally defined right-handed coordinate system used with bending design. Heresis the distance along the design orbit,xis the distance from this orbit along the radius of curvature, andy is the distance from the design orbit out of the bend plane. [29]
5.3 Equation of Motion 107
Fig. 5.4 Particle trajectories
The ideal curvature 0 is evaluated along the reference trajectory u D 0 for a particle with the ideal momentum. In linear approximation with respect to the coordinates the path length element for an arbitrary trajectory is
dsD.1C0u/dzCO.2/; (5.22) where u D x or y is the distance of the particle trajectory from the reference trajectory in the deflecting plane.
The magnetic fields depend onzin such a way that the fields are zero in magnet free sections and assume a constant value within the magnets. This assumption results in a step function distribution of the magnetic fields and is referred to as the hard edge model, generally used in beam dynamics. The path is therefore composed of a series of segments with constant curvatures. To obtain the equations of motion with respect to the ideal path we subtract from the curvature for an individual particle the curvature0of the ideal path at the same location.
Sinceuis the deviation of a particle from the ideal path, we get for the equation of motion in the deflecting plane with respect to the ideal path from Fig.5.4 and (5.20), (5.21) withu00D .d'=dzd'0=dz/,
u00 D .1C0u/C0; (5.23)
where the derivations are taken with respect to z. In particle beam dynamics, we generally assume paraxial beams,u02 1since the divergence of the trajectories u0 is typically of the order of 103 rad or less and terms in u02 can therefore be neglected. Where this assumption leads to intolerable inaccuracies the equation of motion must be modified accordingly.
The equation of motion for charged particles in electromagnetic fields can be derived from (5.23) and the Lorentz force. In case of horizontal deflection, the curvature is D x and expressing the general field by its components, we have
Figure 2.4: Particle trajectories in deflecting systems. Reference pathz and individual particle tra-jectoryshave in general different bending radii [28]
wherek0 is the curvature of the ideal path. The deflection angle for an arbitrary trajectory is then given by
dφ= +kds (2.32)
The ideal curvaturek0is evaluated along the reference trajectory u= 0for a particle with the ideal momentum. In linear approximation with respect to the coordinates the path length element for an arbitrary trajectory is
ds= (1 +k0u)dz+O(2), (2.33)
whereu=xoryis the distance of the particle trajectory from the reference trajectory in the deflecting plane.
The magnetic fields depend onz in such a way that the fields are zero in magnet free sections and assume a constant value within the magnets. This assumption results in a step function distribution of the magnetic fields and is referred to as the hard edge model, generally used in beam dynamics. The path is therefore composed of a series of segments with constant curvatures. To obtain the equations of motion with respect to the ideal path the curvaturek0 of the ideal path is subtracted from the
Physics Department 13 Technische Universität München
14 CHAPTER 2. THEORY curvaturekfor an individual particle at the same location. Sinceuis the deviation of a particle from the ideal path, for the equation of motion in the deflecting plane with respect to the ideal path from Figure2.4 and equation2.31,2.32withu00=−(dφ/dz−dφ0/dz),
u00=−(1 +k0u)k+k0 (2.34)
where the derivations are taken with respect to z. Paraxial beams were assumed, u02 1, since the divergence of the trajectoriesu0 is typically of the order of10−3rad or less and terms inu02 can therefore be neglected. Where this assumption leads to intolerable inaccuracies the equation of motion must be modified accordingly.
The equation of motion for charged particles in electromagnetic fields can be derived from equation 2.34 and the Lorentz force. In case of horizontal deflection, the curvature isk =kx and expressing the general field by its components with equation 2.22
kx= 1 1 +δ
k0x+Kx+1 2mx2
(2.35)
where the field is expanded into components up to second order. The three lowest order multipoles, a bending magnet, a quadrupole and a sextupole were used. A real particle beam is never monochro-matic and therefore effects due to small momentum errors must be considered. This can be done by expanding the particle momentum in the vicinity of the ideal momentump0
1
p= 1
p0(1 +δ)≈ 1
p0(1−δ+. . .) (2.36) The horizontal plane (u=xand k=kx) is now applied to equation 2.34to get with equation 2.35 and2.36the equation of motion
x00+ (K+k20x)x=k0x(δ−δ2) + (K+k20x)xδ−1
2mx2−k0Kx2+O(3) (2.37) The definitions of energy independent field strength parameters as defined in equation 2.30and2.22 was used. The equation of motion in the vertical plane can be derived in a similar way by setting u=y in equation2.34andk=ky. The equation of motion in the vertical plane thereby is
y00−(K−k02y)y=k0y(δ−δ2)−(K−k20y)yδ+1
2my2+k0Ky2+O(3) (2.38)
2.4.4 Solutions for the linear Equations of Motion
Equations 2.37 and 2.38 are the equations of motion for strong focusing beam transport systems, where the magnitude of the focusing strength is a free parameter. No general analytical solutions are available for arbitrary distributions of magnets. The best tool in the mathematical formulation of a solution to the equations of motion is the ability of magnet builders and alignment specialists to build magnets with almost ideal field properties and to place them precisely along a predefined ideal path.
In addition, the capability to produce almost monochromatic particle beams is of great importance for the determination of the properties of particle beams. As a consequence, all terms on the right-hand side of 2.37and 2.38can and will be treated as small perturbations and mathematical perturbation methods can be employed to describe the effects of these perturbations on particle motion.
The left-hand side of the equations of motion resembles that of a harmonic oscillator although with a time dependent frequency. By a proper transformation of the variables equations2.37and2.38can be expressed exactly in the form of the equation for a harmonic oscillator with constant frequency.
Technische Universität München 14 Physics Department
2.5. DIRECTION COSINES AS A BEAM PROPERTY 15