In accelerator physics, the Courant–Snyder parameters (frequently referred to as Twiss parameters or CS parameters) are a set of quantities used to describe the distribution of positions and velocities of the particles in a beam.[1] When the positions along a single dimension and velocities (or momenta) along that dimension of every particle in a beam are plotted on a phase space diagram, an ellipse enclosing the particles can be given by the equation:
\gammax2+2\alphaxx'+\betax'2=\epsilon
where
x
x'
\alpha
\beta
\gamma
\epsilon
The use of these parameters to describe the phase space properties of particle beams was popularized in the accelerator physics community by Ernest Courant and Hartland Snyder in their 1953 paper, "Theory of the Alternating-Gradient Synchrotron".[3] They are also widely referred to in accelerator physics literature as "Twiss parameters" after British astronomer Richard Q. Twiss, although it is unclear how his name became associated with the formulation.[4]
When simulating the motion of particles through an accelerator or beam transport line, it is often desirable to describe the overall properties of an ensemble of particles, rather than track the motion of each particle individually. By Liouville's Theorem it can be shown that the density occupied on a position and momentum phase space plot is constant when the beam is only affected by conservative forces. The area occupied by the beam on this plot is known as the beam emittance, although there are a number of competing definitions for the exact mathematical definition of this property.[5]
In accelerator physics, coordinate positions are usually defined with respect to an idealized reference particle, which follows the ideal design trajectory for the accelerator. The direction aligned with this trajectory is designated "z", (sometimes "s") and is also referred to as the longitudinal coordinate. Two transverse coordinate axes, x and y, are defined perpendicular to the z axis and to each other.[6]
In addition to describing the positions of each particle relative to the reference particle along the x, y, and z axes, it is also necessary to consider the rate of change of each of these values. This is typically given as a rate of change with respect to the longitudinal coordinate (x' = dx/dz) rather than with respect to time. In most cases, x' and y' are both much less than 1, as particles will be moving along the beam path much faster than transverse to it. Given this assumption, it is possible to use the small angle approximation to express x' and y' as angles rather than simple ratios. As such, x' and y' are most commonly expressed in milliradians.
When an ellipse is drawn around the particle distribution in phase space, the equation for the ellipse is given as:
\gammax2+2\alphaxx'+\betax'2=area
\epsilon
\epsilon/\pi
The other three coefficients,
\alpha
\beta
\gamma
\beta\gamma-\alpha2=1
In addition to treating the CS parameters as an empirical description of a collection of particles in phase space, it is possible to derive them based on the equations of motion of particles in electromagnetic fields.
In a strong focusing accelerator, transverse focusing is primarily provided by quadrupole magnets. The linear equation of motion for transverse motion parallel to an axis of the magnet is:
| d2x | |
| dz2 |
=-k(z)x
where
k(z)
Assuming
k(z)
x(z)=A(z)\cos(\phi(z)+\phio)
where A(z) is the amplitude of oscillation,
\phi(z)
k(z)
\phio
\beta
Ao
x(z)=Ao\sqrt{\beta}\cos(\phi(z)+\phio)
x'(z)=Ao
| \beta' | |
| 2\sqrt{\beta |
(It is important to remember that ' continues to indicated a derivative with respect to position along the direction of travel, not time.)
Given these equations of motion, taking the average values for particles in a beam yields:
\langlex2\rangle
| 2 | |
| =A | |
| o |
\beta\langle
| 2(\phi(z)+\phi | |
| \cos | |
| o)\rangle |
=
| 1 | |
| 2 |
| 2\beta | |
| A | |
| o |
\langle{x'}2\rangle=
| 2 | |
| A | |
| o |
| (-\beta'/2)2 | |
| 2\beta |
+
| 2 | |
| A | |
| o |
| 1 | |
| 2\beta |
=(1+(-\beta'/2)2)
| |||||||
| 2\beta |
\langlexx'\rangle=-
| |||||||
| 2 |
(-\beta'/2)
These can be simplified with the following definitions:
\epsilon=
| 1 | |
| 2 |
| 2 | |
| A | |
| o |
\alpha=-
| \beta' | |
| 2 |
\gamma=
| 1+\alpha2 | |
| \beta |
giving:
\langlex2\rangle=\epsilon\beta
\langle{x'}2\rangle=\epsilon\gamma
\langlexx'\rangle=-\epsilon\alpha
These are the CS parameters and emittance in another form. Combined with the relationship between the parameters, this also leads to a definition of emittance for an arbitrary (not necessarily Gaussian) particle distribution:
\epsilon2=\langlex2\rangle\langlex'2\rangle-\langlexx'\rangle2
The advantage of describing a particle distribution parametrically using the CS parameters is that the evolution of the overall distribution can be calculated using matrix optics more easily than tracking each individual particle and then combining the locations at multiple points along the accelerator path. For example, if a particle distribution with parameters
\alpha
\beta
\gamma
\alpha(L)
\beta(L)
\gamma(L)
\begin{pmatrix} \beta(L)\\ \alpha(L)\\ \gamma(L) \end{pmatrix}= \begin{pmatrix} 1&-2L&L2\\ 0&1&-L\\ 0&0&1 \end{pmatrix} \begin{pmatrix} \beta\ \alpha\ \gamma \end{pmatrix}