%written by B. Palmer
% edited by R. Fernow
\section{IONIZATION COOLING}
\label{subsec-compcool}
\subsection{Introduction}
The muon beam at the end of the decay channel is very intense, with approximately $N_{\mu}\approx 7.5\times 10^{12}$~muons/bunch but with a rather large normalized transverse emittance $\epsilon_{x,N}(rms)\approx 15\times 10^3$~mm-mrad and a normalized longitudinal emittance of $\epsilon_{z,N}(rms)\approx 612$~mm.
For a high luminosity collider, the phase-space volume of the muon beam must be reduced to be able to accelerate the beam and inject into the collider and that must be done within a time of the order of the
$\mu$ lifetime $\approx 2~\mu s$. This requires the development of a new method for beam cooling.
Cooling by synchrotron radiation, conventional stochastic
cooling and conventional electron cooling are all too slow. Optical stochastic
cooling\cite{ref21}, electron cooling in a plasma discharge\cite{ref22} and
cooling in a crystal lattice\cite{ref23,ref23a} are being studied, but appear
technologically difficult. Ionization cooling\cite{ref2b,ref24b} of muons seems relatively straightforward, and it is uniquely applicable to muons because of their minimal interaction with matter.
Ionization cooling involves passing the beam through some
material in which the muons lose both transverse and longitudinal momentum
by ionization loss (dE/dx). The longitudinal muon
momentum is then restored by reacceleration, leaving a net loss of
transverse momentum (transverse cooling). The process is repeated
many times to achieve a large cooling factor.
The energy spread can also be reduced (longitudinal emittance exchange)
by introducing a transverse variation in the absorber density or
thickness (e.g.\ a wedge)
at a location where there is dispersion (the position is energy dependent). This method results in an a corresponding increase of transverse momentum and is thus not true cooling, but an exchange of longitudinal and transverse emittances that, with transverse cooling, allows cooling in all dimensions.
We define the rms normalized emittance as
\begin{equation}
\epsilon_{i,N}(rms)=\sqrt{-^2}/m_{\mu}c
\end{equation}
where $x_i$ and $p_i$ are the beam canonical conjugate variables with $i\equiv {x,y,z}$ and $<...>$ stands for the statistical average over the particle distribution.
The appropriate figure of merit for cooling is the final value of the 6-D relativistically invariant emittance $\epsilon_6$ -- the area in the 6-dimensional phase space $(x,y,z,p_x,p_y,p_z)$-- since, to a fairly good approximation, it is preserved during
acceleration and storage in the collider ring. In practical calculations we neglect all correlations and use the following simplified expresssion for 6-D normalized emittance,
\begin{equation}
\epsilon_{6,N}\approx \epsilon_{x,N}\times \epsilon_{x,N}\times \left(\sigma_z^{\textit{rms}}\sigma_E^{\textit{rms}}\gamma \beta\right)
\end{equation}
Theoretical studies have shown that, assuming realistic parameters for the
cooling hardware, ionization cooling can be expected to
reduce the phase-space volume
occupied by the initial muon beam by a factor of $10^5$ -- $10^6$.
A complete
cooling channel would
consist of 20 -- 30 cooling stages, each stage yielding about a factor of
two in phase-space reduction.
It is
recognized that the feasibility of constructing a muon
ionization cooling channel is on the critical path to understanding the
viability of the whole muon collider concept.
The muon cooling channel is the most novel part of a muon collider complex.
Steady progress has been made both in improving the design of the channel
and in adding detail to the computer simulations. A vigorous experimental program is needed to verify and benchmark the computer simulations.
In the following parts of this section we will briefly describe the physics
underlying the process of ionization cooling, we will show results of simulations for some chosen examples and mention a six year R\&D program to demonstrate the feasibility of using ionization cooling techniques to reduce the 6-dimensional phase space emittance.
\subsection{Cooling Theory}
In ionization cooling, the beam loses both transverse and longitudinal momentum
as it passes through a material medium. At the same time its emittance is increased due to stochastic multiple scattering and Landau straggling. Subsequently, the longitudinal
momentum can be restored by reacceleration, leaving a net loss of
transverse momentum.
The approximate equation for transverse cooling (with energies in GeV) is
\begin{equation}
\frac{d\epsilon_n}{ds} = -{1\over \beta^2}\frac{dE_{\mu}}{ds}\ \frac{\epsilon_n}{E_{\mu}} +
\frac{\beta_{\perp} (0.014)^2}{2\beta^3 E_{\mu}m_{\mu}\ L_R}, \label{eq1}
\end{equation}
where $\beta$ is the normalized velocity, $E_{\mu}$ is the total energy,
$m_{\mu}$ is the muon mass, $\epsilon_n$ is the normalized transverse emittance, $\beta_{\perp}$ is the betatron
function at the absorber, $dE_{\mu}/ds$ is the energy loss, and $L_R$ is the
radiation length of the material. The first term in this equation is the
cooling term, and the second is the heating due to multiple scattering.
This heating term is minimized if $\beta_{\perp}$ is small (strong-focusing)
and $L_R$ is large (a low-Z absorber).
The equation for energy spread (longitudinal emittance) is:
\begin{equation}
{\frac{d(\Delta E_\mu)^2}{ds}}\ =
-2\ {\frac{d\left( {\frac{dE_\mu}{ds}} \right)} {dE_\mu}}
<(\Delta E_{\mu})^2 >\ +
{\frac{d(\Delta E_{\mu})^2_{{\rm straggling}}}{ds}}\label{eq2}
\end{equation}
where the first term is the cooling (or heating) due to energy loss,
and the second term is the heating due to straggling.
Energy spread can be reduced by artificially increasing
${d(dE_\mu/ds)\over dE_{\mu}}$ by placing a transverse variation in absorber
density or thickness at a location where position is energy dependent, i.e.\ where there is
dispersion. The use of such wedges can reduce energy spread, but it
simultaneously increases transverse emittance in the direction of the
dispersion. Six dimensional phase space is not reduced, we only exchange longitudinal phase space with transverse phase space.
Ionization cooling of muons seems relatively straightforward in theory, but
will require extensive simulation studies and hardware development for its
optimization. There are practical problems in designing lattices that can
transport, and focus the large emittance beam. There will also be effects
from space charge and wake fields.
We have developed a number of tools for studying the ionization cooling
process. First, the basic theory was used to identify the most promising
beam properties, material type and focusing arrangements for cooling. Given
the practical limits on magnetic field strengths, this gives an estimate of
the minimum achievable emittance for a given configuration.
Next, the differential equations for cooling and heating described above have
been incorporated into a computer code. Allowance for the shifts in the
betatron phase advance due to space charge and aberrations was included.
This code was used to develop an overall cooling scenario, which broke the
cooling system into a number of stages, and determined the properties of the
beam, rf, and focusing lattice at each stage.
Finally, several tracking codes were either written or modified to study
the cooling process in detail. Two new codes (SIMUCOOL\cite{van}, and ICOOL\cite{ref25})
use Monte Carlo techniques to track particles one at a time through the
cooling system. All the codes attempt to include all relevant physical
processes to some degree, e.g.\ energy loss, straggling, multiple scattering,
and to use Maxwellian models of the focusing fields. They do not yet
take into account any space charge or wake field effects. In addition, we
have also used a modified version of PARMELA\cite{parmela} for tracking, which does
include space charge effects, and a double precision version of GEANT\cite{ref40,paul}.
We have recently developed\cite{envelope} a model of beam cooling based on a second
order moment expansion. A computer code solving the equations for
transverse cooling gives results that agree with tracking codes.
The code is being extended to include energy spread and bends. It
is very fast and is appropriate for preliminary design and optimization of
the cooling channel.
All of these codes are continually being
updated and optimized for studying the cooling problem.
\subsection{Cooling System}
The cooling is obtained in a series of cooling stages. Each stage
consists of a succession of two components:
\begin{enumerate}
\item Transverse cooling sections using materials in a strong focusing (low$\beta_\perp$) environment alternated with linac accelerators.
\item Emittance exchange in lattices that generate dispersion, with absorbing material wedges to reduce momentum spread.
\end{enumerate}
In the examples that follow it is seen that each such stage lowers the 6 dimensional emittance by a factor of about 2.
Since the required total 6 dimensional cooling is about 10$^6$, about 20 such stages are required.
The total length of the system would be of the order of 600 m, and the total acceleration would be approximately
6 GeV. The fraction of muons remaining at the end of the cooling system
is estimated to be $\approx 60 \%$.
The baseline solution for transverse cooling involves the use of liquid hydrogen absorbers in strong solenoid focusing fields, interleaved with short linac sections. The solenoidal fields in successive absorbers are reversed to avoid build up of the canonical angular momentum.
The focusing magnetic fields are small ($\approx$ 1 T) in the early stages where the emittances are large, but must increase as the emittance falls.
Two transverse cooling examples have been designed and simulated. The first would lie toward the end of a full cooling sequence and uses 15 T solenoids. The second, using 31 T solenoids, meets the requirements for the Higgs factory and could be the final
cooling stage for this machine.
The baseline solution for emittance exchange involves the use of bent solenoids to generate dispersion and wedges of hydrogen or LiH to reduce the energy spread. A simulated example is given for exchange that would be needed after the 15 T transverse cooling case.
A lithium lens solution may prove more economical for the final stages, and might allow even lower emittances to be obtained.
In this case, the lithium lens serves simultaneously to maintain the low
$\beta_{\perp}$, and attenuate the beam momenta. Similar lenses, with
surface fields of $10\,$T, were developed at Novosibirsk (BINP) and have been used, at low repetition rates, as focusing elements at FNAL and CERN\cite{reviewtev,ref26,ref26a,ref26b,ref27}. Lenses for the cooling application, which would operate at 15 Hz, would need to employ flowing liquid lithium to provide adequate cooling. Higher surface fields would also be desirable.
Studies have simulated cooling in multiple lithium lenses, and have shown cooling through several orders of magnitude. But these studies have, so far, used ideal matching and acceleration. Further studies are underway.
Cooling in beam recirculators could lead to reduction of costs of the cooling section\cite{balbekov}, but again, full simulations with all higher order effects have not yet been successfully achieved.
%We require a reduction of the normalized transverse emittance by almost three
%orders of magnitude (from $1\times 10^{-2}$ to $5\times 10^{-5}\,$m-rad for %the high energy collider), and a
%reduction of the longitudinal emittance by one order of magnitude.
\subsection{15 T Solenoid Transverse Cooling Example}
The lattice consists of 11 identical 2~m long \textit{cells}. In each cell there is a liquid hydrogen absorber (64~cm long, 10~cm diameter) in the 15~T solenoid focusing magnet (64~cm long, 12~cm diameter ). The direction of the fields in the magnets alternates
from one cell to the next.
Between the magnets there are magnetic matching sections (1.3~m long, 32~cm inside diameter) where the field is lowered and then reversed. Inside the matching sections are short 805~MHz high gradient (36~MeV/m) linacs.
Figure~\ref{altsol} shows the cross section of one cell of such a system, together with the beta function, and the magnetic field along the axis. For
convenience in modelling, the cell shown in Fig.~\ref{altsol}(a) starts and ends
symmetrically in the middle of hydrogen absorber regions at the location of
the peaks in the axial magnetic field. In reality each cell would start at
the beginning of the hydrogen region and extend to the end of the rf module.
\begin{figure*}[bht!]
\centerline{\epsfig{file=radii_beta_axial.eps,height=5.0in,width=5.0in}}
\caption[Cross section of one cell of an alternating solenoid cooling system ]{(a) Cross section of one cell of an alternating solenoid cooling system; (b) axial magnetic field {\it vs.} z; (c) $\beta $~function {\it vs.} z. }
\label{altsol}
\end{figure*}
A GEANT simulation of muons traversing a section of the cooling section is shown in Fig.~\ref{dpgeant}.
\begin{figure*}[hbt!]
\centerline{\epsfig{file=dpgeant.eps,height=3.0in,width=6.0in}}
\caption[GEANT simulation of muons traversing a section ]{GEANT simulation of muons traversing a section of the alternating solenoid cooling channel. The variation of the magnetic field B$_z$ is shown for $1{1\over 2}$ cells of the figure.}
\label{dpgeant}
\end{figure*}
Additional simulations were performed using the program ICOOL\cite{ref25}; the only effects which are not yet included
are space charge and wake fields. Analytic vacuum calculations indicate that these effects should, for the later stages be significant, but not overwhelming, but a correct simulation must be done before we are assured that no significant problems exist. Particles are introduced with transverse and longitudinal emittance (186~MeV/c, 1400~$\pi$~mm~mrad transverse, and 1100~$\pi$~mm~longitudinal), together with a number of naturally occuring correlations: 1) the angular momentum appropriate
for the starting axial magnetic field, 2) a correlation between momentum and amplitude squared to give forward velocities independent of amplitude, and 3) a distortion of the longitudinal bunch distribution to reflect the asymmetric nature of the alpha rf bucket.
Figure~\ref{simu8}(a) shows the average momentum of the beam as a function of distance along the channel. The momentum drops as the beam crosses the liquid hydrogen absorbers. The gradient and phase
of the rf cavities have been adjusted so that the reacceleration given to
the reference particle just
equals the mean energy loss. This causes the average momentum of the beam to
remain in a
narrow band around 195~MeV/c. Figure~\ref{simu8}(b) shows the mechanical and canonical angular
momenta as a function of distance along the channel. The mechanical angular
momentum shows
the rotational motion of the beam around the axial solenoidal field. It
periodically reverses sign
when the solenoids alternate direction. The canonical angular momentum is
defined such that it
removes the axial field dependence\cite{aac_paper}. Without the
absorbers the
beam would have a constant (0) value for the canonical angular momentum.
However, the
presence of absorber causes the canonical angular momentum to grow and would
lead to severe
emittance growth at the end of a long channel. This growth is stopped by
alternating the direction
of the solenoid field, as shown in the figure. Simulations have shown that 2~m is a reasonable
(half) period for the field since the net growth in canonical angular
momentum is small. In
addition synchrobetatron resonances are avoided since it forces the average
betatron wavelength
to be 2~m, whereas the synchrotron oscillation wavelength is $\approx 14$~m.
\begin{figure*}[bht!]
\centerline{\epsfig{file=simu8.ps,height=4.0in,width=6.0in}}
\caption[Average momentum \textit{vs.} length ]{a) Average momentum \textit{vs.} length; b) Average angular momentum: mechanical (solid curve) and canonical (dashed curve), \textit{vs.} length.}
\label{simu8}
\end{figure*}
Figure~\ref{simu6}(a) shows the rms and maximum radius of any particle in the beam
distribution as a
function of distance along the channel. The rms radius shows that most of
the beam is confined
to within 2~cm of the axis. The peak rms radius decreases towards the end of
the channel as a
result of the cooling. The maximum particle radius is about 8~cm, which
determines the radius of
the windows required in the rf cavities.
Figure~\ref{simu6}(b) shows the rms momentum
spread corrected for
the correlation between P$_z$ and transverse amplitude imposed on the initial
particle distribution.
The momentum spread grows as a function of distance since the alternating
solenoid system only
cools the transverse emittance.
Figure~\ref{simu6}(c) shows the rms bunch length as a
function of distance
along the channel. Again this grows with distance since this channel does
not cool longitudinally.
\begin{figure*}[bht!]
\centerline{\epsfig{file=simu6.ps,height=4.0in,width=6.0in}}
\caption[rms and maximum beam radii ]{a) rms and maximum beam radii, b) rms corrected momentum, c) rms bunch length; all \textit{vs.} length.}
\label{simu6}
\end{figure*}
Figure~\ref{coolingeg}(a) shows the decrease in transverse normalized emittance as a
function of distance along
the channel. The system provides cooling by a factor of $\approx $~2 in both the x and y transverse phase
spaces. From the slope of the line we note that the rate of cooling is
dropping. This sets $\approx$~22~m as
the maximum useful length for this type of system. It must be followed by a
longitudinal
emittance exchange region to reduce the momentum spread and bunch length
approximately
back to their starting values. Figure~\ref{coolingeg}(b) shows the increase in longitudinal
normalized emittance
in the channel due to the increase in momentum spread and bunch length. Finally, Fig.~\ref{coolingeg}(c) shows the decrease in the 6-dimensional normalized emittance as a function
of distance along the
channel. There is a net decrease in 6-dimensional emittance by a factor of
$\approx$~2 in the channel. Table~\ref{cooltab} gives the initial and final parameters.
\begin{figure*}[bht!]
\centerline{\epsfig{file=simu7.ps,height=4.0in,width=6.0in}}
\caption[Emittance \textit{vs.} length ]{Emittance \textit{vs.} length: a) transverse emittance; b) longitudinal emittance; and c) 6-dimensional emittance.}
\label{coolingeg}
\end{figure*}
\begin{table*}[thb]
\caption{Initial and final beam parameters in a 15 T transverse cooling section.}
\label{cooltab}
\begin{tabular}{llccc}
& & initial & final & factor\\
\hline
Particles tracked & & 1000 & 980 & 98\% \\
Reference momentum & MeV/c &186 & 186 & \\
Transverse Emittance & $\pi$ mm mrad & 1400 & 600 & 0.43 \\
Longitudinal Emittance & $\pi$ mm mrad & 1100 & 2300 & 2.09 \\
6D emittance & $10^{-12}\ \times\ (\pi\ m\ rad)^3$ & 2000 & 800 & 0.40 \\
rms beam size in hydrogen & cm & 0.8 & 0.55 & 0.69 \\
rms beam size in linac & cm & 2.0 & 1.4 & 0.70 \\
max beam rad in linac & cm & 7.0 & 7.0 & 1.0 \\
rms bunch length & cm & 1.5 & 2.2 & 1.5 \\
max bunch full width & cm & 13 & 19 & 1.5 \\
rms dp/p & \% & 3.8 & 5.6 & 1.5 \\
\end{tabular}
\end{table*}
This simulation has been confirmed, with minor differences, by the codes double precision GEANT\cite{paul} and PARMELA\cite{parmela}.
\subsection{31 T Solenoid Transverse Cooling Example}
As above, the lattice consists of 11 identical 2~m long \textit{cells}. The maximum fields are higher (31~T) but the bore is smaller (8~cm). The
liquid hydrogen absorber is slightly shorter then in the previous case (64~cm long, 6~cm diameter). As before,
the direction of the fields in the magnets alternates from one cell to the next.
Between the magnets there are again magnetic matching sections (1.3~m long, 32~cm inside diameter) and the rf is also the same (805~MHz, 36~MeV/m).
Table~\ref{cooltab2} gives the initial and final parameters for
the 31~T example, together with the required emittances for a Higgs factory. In setting these requirements a dilution of $20\,\%$ is assumed in each of the three emittances.
\begin{table*}[thb]
\caption{Initial and final beam parameters in a 31 T transverse cooling section.}
\label{cooltab2}
\begin{tabular}{ll|ccc|c}
& & initial & final & factor& Req.\\
\hline
Particles tracked & & 4000&3984 &99.6\%& \\
Reference momentum & MeV/c & 186 & 186 & \\
Transverse Emittance & $\pi$ mm mrad & 460 & 240 & 0.52 & 240 \\
Longitudinal Emittance & $\pi$ mm mrad & 850 & 1600 & 1.9& \\
6D emittance & $10^{-12}\ \times\ (\pi\ m\ rad)^3$ & 150 & 95 & 0.63&98 \\
rms beam size in hydrogen & cm & 0.44 & 0.33 & 0.75& \\
rms beam size in linac & cm & 1.4 & 1.1 & 0.80& \\
max beam rad in linac & cm & 6.0 & 6.0 & 1.0& \\
rms bunch length & cm & 1.5 & 1.8 & 1.2& \\
max bunch full width & cm & 11 & 19 & 1.7& \\
rms dp/p & \% & 3.5 & 5.0 & 1.4& \\
\end{tabular}
\end{table*}
\subsection{Bent Solenoid Emittance Exchange Example}
The exchange of longitudinal and transverse emittance requires dispersion in a large acceptance channel, followed by low Z wedge absorbers. In this example, the lattice consists of a continuous (3.5~T) solenoid and dispersion is generated by bending the
solenoid.
In a bent solenoid, in the absence of any dipole field, there is a drift perpendicular to the bend that is proportional to the particle's momentum. In this lattice, a uniform dipole field is introduced over the bend to exactly cancel this drift for particles with the reference momentum. Particles above and below this reference momentum then drift up or down, depending on their momentum, giving the required dispersion (0.4 m). The dispersion is removed, and the momentum spread reduced, by introducing liquid hydrogen wedges. The hydrogen wedges
would be contained by thin beryllium or aluminum foils, but these were not included in this simulation.
After one bend and one set of wedges, the beam is asymmetric in cross section. Symmetry is restored by a following bend and wedge system rotated by 90 degrees with respect to the first. Figure~\ref{cell_exch} shows a representation of the two bends and wedges. The total solenoid length was 8.5 m. The beam tube outside diameters are 20 cm, and the minimum bend radii 34 cm.
\begin{figure*}[tbh!]
\centerline{\epsfig{file=double_bent_sol.ps,height=5.5in,width=5.5in}}
\vspace{0.5cm}
\caption{Representation of a bent solenoid longitudinal emittance exchange section}
\label{cell_exch}
\end{figure*}
Figure~\ref{simu4}(a) shows the magnetic fields (B$_z$, B$_y$, and B$_x$) as a function of the position along the cell. The solenoid bend curvature is exactly that given by the trajectory of a reference particle (equal in momentum to the average momenta given in Fig.~\ref{simu4}(b))
in the given transverse fields. The actual shape of the bend turns out to be very important.
Discontinuities in the bend radius can excite Larmor oscillations
which increase the transverse emittance. We have shown, for example,
that the transverse emittance growth in a bent solenoid depends on
discontinuities of the bend radius as a function of distance, and
its first and second derivatives, the size and tilt of the solenoidal
coils, auxiliary fields and the six dimensional phase space of the beam,
thus optimization is not straightforward. One solution to this problem is to have long, adiabatic bends,
however this adds unnecessary length to the emittance exchange section.
We are studying options with coupling sections to tight bends roughly
half a Larmor length long, which seems to minimize transverse
emittance growth while also minimizing the length of the section.
Due to similar problems, the length and longitudinal distribution of
the wedge material has also been found to effect emittance growth and
similar solutions have been successfully considered. For example
emittance growth can be minimized when the vector sum of the Larmor
phases at the absorber elements is small or zero.
\begin{figure*}[bht!]
\centerline{\epsfig{file=simu4.ps,height=4.in,width=6.in}}
\caption[Axial B$_z$ and dipole B$_y$, B$_x$, magnetic fields ]{a) Axial B$_z$ and dipole B$_y$, B$_x$, magnetic fields; b)
average momentum; both as a function of the position along the cell.}
\label{simu4}
\end{figure*}
The simulations were performed using the program ICOOL\cite{ref25}. The maximum beam radius is 10~cm. Transmission was 100 \%.
Figure~\ref{simu2}(a) shows the rms longitudinal momentum spread relative to the reference momentum as a function of the position along the cell. The functional spread decreases from an initial value of approximately 5\%, to a final value of approximately
2.2\%. At the same time, since this is an emittance exchange, the transverse beam area grows, as shown in Fig.~\ref{simu2}(b). One notes that the area increases not only in the regions of bends (region 1 and 8), but also in the regions of wedges (2-6 and 9-11). This is probably due to failures in matching that have yet to be understood.
\begin{figure*}[bht!]
\centerline{\epsfig{file=simu2.ps,height=4.in,width=6.in}}
\caption[rms longitudinal $\delta p_z$ with respect to the reference momentum ]{a) rms longitudinal $\delta p_z$ with respect to the reference momentum and b) transverse beam area, both as a function of length. }
\label{simu2}
\end{figure*}
Figure~\ref{xfig5} shows scatter plots of the transverse particle positions against their momenta: a) y-p at the start, b) y-p after the first bend c) y-p after the first set of wedges, d) x-p after the first wedges, e) x-p after the second bend, and f) at the end.
The dispersion is clearly observed in Fig.~\ref{xfig5}(b) (after the first bend) and in Fig.~\ref{xfig5}(e) (after the second). It is seen to be removed, with a corresponding decrease in momentum spread, in Fig.~\ref{xfig5}(c) (after the first set of wedges) and Fig.~\ref{xfig5}(f) (after the second set of wedges).
\begin{figure*}[thb!]
\centerline{\epsfig{file=simu5.ps,height=5.5in,width=5.5in}}
\caption[y-p$_z$ plots and x-p$_z$ plots before and after a wedge ]{y-p$_z$ plots: a) at the start, b) after the first bend, c) after the first set of wedges. x-p$_z$ plots: d) after the first wedges, e) after the second bend, and f) at the end}
\label{xfig5}
\end{figure*}
Figure~\ref{xfig1} shows a scatterplot of the particle squared radii vs.\ their longitudinal momentum, (a) at the start, and (b) at the end. The decrease in momentum spread and rise in beam area are clearly evident.
The initial and final beam parameters are given in Tb.~\ref{xtab}. The example was not well optimized. There are significant final beam fluctuations, presumably due to mismatching, and some of the 37\% rise in 5 dimensional phase space may be from this
cause. RF must be included to allow the 6 dimensional emittance to be studied. Significantly more work is needed to further optimize the emittance exchange section.
\begin{figure}[tbh!]
\centerline{\epsfig{file=simu1.ps,height=2.5in,width=3.5in}}
\caption[Scatterplot of squared radii vs.\ longitudinal momentum ]{Scatterplot of squared radii vs.\ longitudinal momentum a) at the start, and b) at the end.}
\label{xfig1}
\end{figure}
\begin{table*}[bth!]
\caption{Initial and final beam parameters in a longitudinal emittance exchange section.}
\label{xtab}
\begin{tabular}{llccc}
& & initial & final & factor\\
\hline
Longitudinal Momentum spread & MeV/c & 9.26 & 3.35 & 0.36 \\
Ave. Momentum & MeV/c & 180 & 150 & 0.83 \\
Transverse size & cm & 1.33 & 2.26 & 1.70 \\
Transverse Momentum spread & MeV/c & 6.84 & 7.84 & 1.15 \\
Transverse Emittance & $\pi$ mm mrad & 870 & 1694 & 1.95 \\
Emit$_{trans}^2$ $\times$ $\Delta$p$_{long}$& ($\pi$ m mrad)$^2$ MeV/c & 7.0 & 9.6 & 1.37 \\
\end{tabular}
\end{table*}
Emittance exchange in solid LiH wedges, with ideal dispersion and matching, has also been successfully simulated using SIMUCOOL\cite{dave}. Dispersion generation by weak focusing spectrometers\cite{balbekov} and dipoles with solenoids\cite{wan} have also been studied.
\subsection{RF for the Cooling Systems}
RF acceleration sections are used to
restore the longitudinal momentum of the muon beam
after losses in the cooling media.
These rf structures are embedded in solenoidal fields which
reverse direction within each section.
% and vary in strength
%from 0 to 5~T.
In the two transverse cooling examples, above, the rf frequency is 805 MHz and the peak gradient is 36~MeV/m. The focusing fields, which extend over the cavities, vary from 0 to 5~T, reversing in the center. It should be pointed out that in the earlier and later stages, the bunches are longer, and lower frequencies will be required.
In order to realize maximum acceleration gradients within the acceleration
cavities, we take advantage of the penetrating properties of a muon
beam by placing thin windows between each rf cell, thereby creating an
accelerating structure closely approximating the classic pill-box cavity. This permits operating conditions in which the axial accelerating field is equal to the maximum wall field and gives a high shunt impedance.
The windows in the 15~T example are
16~cm diameter, $125~\mu$m thick Be foils. In the 31~T case, they are 10~cm diameter and $50~\mu$m thick.
%\begin{center}
\begin{figure*}[thb!]
\centering
\epsfig{file=rf_a.eps,height=3.5in,angle=-90}
\vspace{-2.5cm}
\caption[Two full cell sections plus two half cell sections ]{Two full cell sections plus two half cell sections
of the interleaved $\pi$/2 mode accelerating cavities. The
volumes labeled \textbf{C} are powered separately from the volumes labeled \textbf{D}.}
\label{rf_inter}
\end{figure*}
%\end{center}
\begin{figure*}[thb!]
\centering
\epsfig{file=rf_1.3m.eps,height=12cm,angle=-90}
\caption[A 1.3 m acceleration section ]{A 1.3 m acceleration section with quarter section
cut away for viewing the inter cavity windows.}
\label{rf_section}
\end{figure*}
For these rf structures, we will use an interleaved cavity design in which
two parts are independently powered (Fig.~\ref{rf_inter}).
The mode of the system will be referred to
as $\pi$/2 interleaved.
Each section supports a standing wave $\pi$ mode, with each acceleration
cell $\pi$/2 long, giving a good transit time factor.
To reduce the peak rf power requirements (by a factor of 2) the cells may be operated
at liquid nitrogen temperatures.
%This also facilitates maximizing the shunt impedance.
\begin{table}[htb!]
\centering{\caption{Characteristics of the rf system.}
\begin{tabular}{lc}
\multicolumn{2}{c}{RF Cavity Parameters} \\
\noalign{\vspace{2pt}}
\hline
\noalign{\vspace{2pt}}
& \\
RF frequency [MHz] & 805 \\
Cavity Length [cm] & 8.1 \\
Cavity Inner Radius [cm] & 14.65 \\
Cavity Outer Radius [cm] & 19 \\
Q/1000 & 2 $\times$ 20 \\
Peak Axial Gradient [MV/m] & 30 \\
Shunt Impedance [M$\Omega$/m] & 2 $\times$ 54 \\
Zt$^2$ [M$\Omega$/m] & 2 $\times$ 44 \\
Fill Time (3 $\tau$) [$\mu$s] & 26 \\
RF Peak Power [MW/m] & 16/2 \\
Ave. Power (15Hz) [KW/m] & 7/2 \\
Be window aperture [cm] & 16 (10 for 31 T case) \\
Be window thickness [$\mu$m] & 127 (50 for 31 T case) \\
\end{tabular}
\label{rf_table}
}
\end{table}
The characteristics of the rf systems currently being studied are summarized
in Table~\ref{rf_table}. Figure~\ref {rf_section} shows a full 1.3~m section with interleaved cavities.
Each cell is 8.1~cm in length and the 1.3~m section consists of 16~cells.
\subsection{The Liquid Lithium Lens}
The final cooling element ultimately determines the luminosity of the
collider. In order to obtain smaller transverse emittance as the muon
beam travels down the cooling channel, focusing strength must increase
and $\beta_{\perp}$'s must
decrease. A current within a conductor produces an active
lens absorber, which can maintain the beam at small $\beta_{\perp}$
throughout
an extended absorber length, while simultaneously attenuating
the beam momentum. An active lens absorber, such as a
lithium lens, may prove to be the most efficient cooling element
for the final stages.
The cooling power of a Li lens is illustrated
in Figure~\ref{lilesim}, where the $x-p_x$ phase space distributions at
the beginning and at the end of the absorber are shown. This example
corresponds to a 1~m long lens, with 1~cm radius, and a surface field of
10~T. The beam momentum entering the lens was 267~MeV/c, with Gaussian
transverse spatial and momentum distributions:
$\sigma_x = \sigma_y = 2.89$ mm,
$\sigma_{p_x}=\sigma_{p_y} = 0.0267$ GeV/c, and a normalized emittance of
$\epsilon_{x,N}^{\textrm{in}} = 710$~mm~mrad. The normalized emittance at the end of the absorber was $\epsilon_{x,N}^{\textrm{out}} = 450$~mm~mrad (cooling factor
$\sim 1.57$), and the beam momentum
159 MeV/c.
The results were obtained using a detailed GEANT simulation of a single
stage. In order to construct a Li lens cooling channel, a series of
lenses has to be used; the parameters of the lens would have to vary to
match the changing beam emittance, and
acceleration of the beam in between the lenses has to be included.
\begin{figure}[hth!]
\epsfysize=3.5in
\vspace{-1.0cm}
\centerline{\epsffile{lilenssim.eps}}
\caption{$x-p_x$ phase space distribution at the beginning and at the
end of the absorber described in the text.}
\label{lilesim}
\end{figure}
Lithium lenses have been used with high reliability
as focusing elements at FNAL and CERN [70-74]. Although these lenses
have
many similar properties to those required for ionization cooling, there
are
some very crucial differences which will require significant advances in
liquid lens technology: ionization cooling requires longer lenses
($\sim$ 1 meter), higher fields ($\sim$ 10 $T$), and higher operation
rates (15 Hz). The last requirement calls for operating the lenses with
lithium in the liquid phase.
A liquid Li lens consists of a small diameter rod-like chamber filled
with
liquid Li through which a large current is drawn. Each lens is assumed
to be
about a meter long and we plan to use a number of lenses
mounted sequentially. The azimuthal magnetic field
focuses the beam to give the minimum achievable emittance
$\epsilon_{x,N} \approx C \beta_{\perp}$ where the constant $C$ depends on
the properties of the material, for example $C_{Li}=79$ mm~mrad/cm.
The focusing term can be written as
$\beta_{\perp} \sim 0.08$[cm]$ \sqrt{p/J}$
with $p$ is the muon momentum in MeV/c, and $J$
is the current density in MA/cm$^2$. Increasing $J$ is desirable,
and decreasing $p$ can be useful, however below about 250 MeV/c
the slope of ${dE\over dx}(E)$ tends to increase longitudinal emittance.
The requirement for the highest current density guarantees
large ohmic power deposition, and the current density will be
limited by the maximum tolerable deposited energy, which
will produce instantaneous heating, expansion, and pressure effects.
Understanding the above effects is part of the ongoing liquid Li lens
development R\&D.
The structural design of the lithium lens is determined by how the
pressure
pulse and heat deposition are handled. We assume that the Li will be
flowing rapidly under high pressure, confined by electrical
insulators radially and fairly thick Be windows longitudinally.
Operation at 15 Hz for long periods poses severe challenges. Shock,
fatigue and other failure modes are being evaluated, in addition
to studies of material compatibilities, corrosion and degradation to
insure
safe operation over long periods. It seems that the
minimum required radius of the lens may be the most important parameter
to
determine,
since mechanical problems rise while losses decrease as a function of
radius.
Transferring the beam from one lens to another, with linacs to
reaccelerate and provide longitudinal focussing,
is also a challenging problem, because of the multiple scattering
introduced
in the windows,straggling and the large divergence of the beams.
We are in the process of evaluating a number of designs for this
transfer
channel, using solenoids, quadrupoles and other focusing elements
together with
Li lenses into detailed tracking simulations.
A group from BINP, ANL and FNAL have begun to design and develop a
liquid lens
prototype. The goal is to study
the possible solutions to the problems of confining and
minimizing the pressure pulse, and the heat deposition effects, and
to incorporate the best idea into the design of the prototype.
The group plans a study program at BINP and
ANL, using simulations and experimental tests of small prototypes.
The design of two lenses, whose behavior will be tested at first on a
bench
and then with muon beams at the
{\it Ionization Cooling Demonstration Facility}, Fermilab P904, will
then
follow.
\subsection{Ionization Cooling Experiment}
An R\&D program has been proposed to
design and prototype the critical sections of a muon ionization
cooling channel. The goal of this experimental R\&D program
is to develop the muon ionization cooling hardware to the point
where a complete
ionization cooling channel can be confidently designed for the First Muon
Collider. Details can be found in Fermilab proposal
P904~\cite{steveproposal}. Following is a brief description of the main elements of the experimental program.
The proposed R\&D program consists of:
\begin{itemize}
\item Developing the appropriate rf re-acceleration structure.
It is proposed to construct a 3-cell prototype rf cavity with
thin beryllium windows, which will be tested at high power at
liquid nitrogen temperature and within a high-field solenoid.
\item Prototyping initially a 2~m section, and eventually a 10~m section,
of an alternating solenoid
transverse cooling stage. It is proposed to test the performance of
these sections in a muon beam of the appropriate momentum.
\item Prototyping an emittamce exchange (wedge) section and measuring its
performance in a muon beam of the appropriate momentum.
\item Prototyping and bench testing
$\sim 1$~m long liquid lithium lenses, and
developing lenses with the highest achievable surface fields, and hence
the maximum radial focusing.
\item Prototyping a lithium lens--rf--lens system and measuring its
performance in a muon beam of the appropriate momentum.
\item Developing, prototyping, and testing a hybrid lithium lens/wedge
cooling system.
\end{itemize}
The measurements that are needed to demonstrate the cooling capability and
optimize the design of the alternating solenoid, wedge,
and lithium lens cooling stages will
require the construction and operation of an ionization cooling test facility.
This facility will need
\begin{enumerate}
\item a muon beam with a central momentum that can be
chosen in the range 100 -- 300~MeV/c,
\item an experimental area that can
accommodate a cooling and instrumentation setup of initially $\sim 30$~m
in length, and eventually up to $\sim 50$~m in length, and
\item instrumentation to precisely measure the positions of the incoming
and outgoing particles in six-dimensional phase-space and confirm that they
are muons.
\end{enumerate} In the initial design shown in Fig.~\ref{cooling_expt},
the instrumentation consists of
identical measuring systems before and after the cooling apparatus\cite{lumcdpre}.
Each measuring system consists of
(a) an upstream time measuring device to determine the arrival time of the
particles to one quarter of an RF cycle ($\sim\pm300$~ps),
(b) an upstream momentum spectrometer in
which the track trajectories are measured by low pressure TPC's on either
side of a bent solenoid,
(c) an accelerating RF cavity to change the
particles momentum by an amount that depends on its arrival time,
(d) a downstream momentum spectrometer, which is identical to the upstream
spectrometer, and together with the RF cavity and the upstream spectrometer
forms a precise time measurement system with a precision of a few ps. The
measuring systems are 8~m long, and are contained within a high-field
solenoidal channel
to keep the beam particles within the acceptance of the cooling apparatus.
It is proposed to accomplish this ionization cooling R\&D program
in a period
of about 6~years. At the end of this period we believe that it will be
possible to assess the feasibility and cost of constructing an ionization
cooling channel for the First Muon Collider, and if it proves feasible,
begin a detailed design
of the complete cooling channel.
\begin{figure}[thb!]
%\noindent
\begin{minipage}{.85\linewidth} % fig
\centering\epsfig{figure=muonexp7a_c.eps,width=0.85\linewidth}
\end{minipage}\hfill
\begin{minipage}{.85\linewidth} % fig
\centering\epsfig{figure=muonexp7_c.eps,width=0.85\linewidth}
\end{minipage}
\caption{Schematic of the cooling test apparatus arrangement.}
\label{cooling_expt}
\end{figure}
%%%% TH LINES BELOW ARE GOOD FOR TWOCOLUMN FORMAT
%\begin{figure*}[thb!]
%\noindent
%\begin{minipage}{.85\linewidth} % fig
%\begin{minipage}{2.05\linewidth} % fig
%\centering\epsfig{figure=muonexp7a_c.eps,width=2.05\linewidth}
%\end{minipage}\hfill
%\begin{minipage}{2.05\linewidth} % fig
%\centering\epsfig{figure=muonexp7_c.eps,width=2.05\linewidth}
%\end{minipage}
%\caption{Schematic of the cooling test apparatus arrangement.}
%\label{cooling_expt}
%\end{figure*}