Impact of Detector Solenoid on the CLIC Luminosity Performance

In order to obtain the necessary luminosity with a reasonable amount of beam power, the Compact Linear Collider (CLIC) design includes an unprecedented collision beam size of {\sigma} = 1 nm vertically and {\sigma} = 45 nm horizontally. Given the small and very flat beams, the luminosity can be significantly degraded from the impact of the experimental solenoid field in combination with a large crossing angle. Main effects include y-x'-coupling and increase of vertical dispersion. Additionally, Incoherent Synchrotron Radiation (ISR) from the orbit deflection created by the solenoid field, increases the beam emittance. A detailed study of the impact from a realistic solenoid field and the associated correction techniques for the CLIC Final Focus is presented. In particular, the impact of techniques to compensate the beam optics distortions due to the detector solenoid main field and its overlap with the final focus magnets are shown. The unrecoverable luminosity loss due to ISR has been evaluated, and found to be in the range 4-5 % for the solenoid design under study.


I. INTRODUCTION
The Compact Linear Collider (CLIC) is an accelerator design based on normal conducting components. In order to obtain the required luminosity with reasonable power consumption, short bunch separation (0.5 ns) and small β Ã are needed. A post-collision beam line for the spent beam and the main beam-beam products is necessary, which requires a large crossing angle of around 20 mrad [1]. A large crossing angle is also required to mitigate the effects of the parasitic bunch collisions between the incoming and outgoing beam, but is not the limiting factor for CLIC [1].
Relevant parameters for the CLIC beam delivery system (BDS) are shown in Table I. A detailed overview of the CLIC BDS can be found in [[2], Chap. 3.5]. The CLIC final focus system has strict tolerances, and the BDS is optimized taking higher order terms into consideration [3][4][5][6][7]. Compensation of beam distortions in the BDS, such as static and dynamic misalignments, has proven to be quite challenging [2,8].
In Fig. 1, the final 20 m of the CLIC BDS lattice is shown. The residual field from an experimental solenoid typically extends 10-15 m away from the interaction point (IP), depending on shielding and solenoid design. L Ã is the distance from the IP to the closest focusing magnet, which is the QD0 for CLIC. Due to the short L Ã required to reach the luminosity target, the main solenoid field overlaps with the last final focus magnets, which enhances the optical distortions at the IP [9]. The QD0 is partly inside the experimental solenoid. Special care has to be taken to make sure the interplay between the solenoid field and the magnet field is minimized.
A solenoid will in general have a radial field component on any charged particle off the solenoid center line, with a maximum around the entrance of the solenoid. This is the region of maximum β-function in a linear collider, and the beams are more sensitive to small errors. With a horizontal crossing angle, the horizontal solenoid field component will be larger than the vertical one, resulting in a strong vertical orbit displacement. In CLIC this orbit offset is typically on the order of 10 μm, for a solenoid field of 4-5 T and 1.5 TeV beam energy. The displacement results in a large vertical dispersion at the interaction point (IP). Furthermore, the beams in CLIC are exceptionally flat, which means that any coupling to the vertical plane significantly deteriorates the luminosity. Particles with large angles at the IP have a large displacement from the beam orbit in the region close to the last focusing magnet, where the radial solenoid field is strongest. Hence, the experimental solenoid introduces strong y-x 0 coupling at the IP which must be corrected.
Due to the high beam energy in CLIC, there is a significant emission of synchrotron radiation as a result of the beam deflection in the solenoid region. Earlier similar studies have shown an unrecoverable luminosity loss due to ISR of up to 25%, depending on the detector solenoid design [2,10].
The unrecoverable loss is an important concern for CLIC. Optical aberrations can be corrected in several ways; using the final focus magnets, adding skew quadrupoles, using an antisolenoid [9], dipole orbit corrector integrated into the experiment [11], and longer L Ã [12].
We present a new simulation approach which evaluates the effect of the ISR alone without the knowledge of the full compensation. This approach is verified with a semianalytical approach, as well as a more time-consuming study where the full compensation is found.
For the latter study, the tuning methods described in [8] are used to compensate for optical distortions introduced by the experimental solenoid field. A realistic design of the solenoid and antisolenoid is used [13]. Similar correction schemes have been explored for. e.g., the NLC [9], but at lower beam energies which means synchrotron radiation effects are less significant.
There are two problems with the full compensation study which are addressed with the new simulation approach. First of all, it is a computationally demanding procedure, requiring on the order of weeks of CPU time to get to the final result. Second, once the result is obtained, one does not know if the remaining luminosity loss is purely due to ISR, or if there are residual optical aberrations.

II. EXPERIMENTAL SOLENOID FIELD
Two detectors will be running in a push-pull configuration in CLIC. In the conceptual design report it is foreseen that one detector will follow the SiD design [14], while the second detector will have The International Large Detector (ILD) design [15]. An important difference in the two magnet designs is the peak longitudinal field, which is 4 T for the ILD detector magnet, FIG. 2. The longitudinal (a) and radial (b) SiD solenoid field with (blue, solid) and without (green, dashed) antisolenoid, along a beam line with a 10 μrad inclination with respect to the solenoid axis. The QD0 entrance is at 3.5 m, and the IP is at 0 m. and 5 T in the SiD case. Nevertheless, previous studies have found that the luminosity loss from the two detector designs is fairly similar due to the relative increase in stray fields from the ILD solenoid compared to the SiD solenoid [10].
The longitudinal and radial fields along the beam line for the SiD detector magnet are shown in Fig. 2, both with (blue, solid) and without (green, dashed) the antisolenoid [14]. For all detector designs currently considered, the antisolenoid is foreseen to be integrated into the CLIC detector. This is in contrast to, e.g., the International Linear Collider (ILC) [16], where the antisolenoid is integrated into the QD0 design [17]. The antisolenoid significantly reduces the longitudinal field inside the QD0, increasing the radial field at the entrance of the QD0 (at 3.5 m). This reduces the optical aberrations originating from the combination of the quadrupolar field of the QD0 and the solenoid stray fields [9].
In this paper we discuss only the L Ã ¼ 3.5 m lattice design. An increase of L Ã to around 6-8 m has been considered in order to have the QD0 outside the detector [12,18]. This would have the benefit of an improved mechanical stability of the QD0 since it would be attached to the tunnel ground/beam line. It would further reduce the interference between the QD0 and the experiment, both in terms of taking up physical space for detector components, and in terms of magnetic field interplay. The downside of a longer L Ã is primarily reduced luminosity. A different L Ã is expected to significantly impact the result of the studies presented here.

III. SEMIANALYTICAL APPROACH
The problem of evaluating the luminosity loss due to the detector solenoid can be divided into two parts. The first part consists of the evaluation and the correction of the optical distortions, which should be possible to correct for by using the antisolenoid and tuning knobs. The second part, the ISR from the vertical orbit deflection, increases the beam emittance. This emittance increase cannot be compensated for, and can be considered a minimum luminosity loss for a given solenoid design.
From Refs. [19,20], we have the following estimate for the increase of vertical beam size due to synchrotron radiation where Here, R 36 is the transport matrix element 36 for the given slice dz to the IP, γ is the relativistic gamma, and ρðzÞ is the radius of curvature at z. r e and λ e are the electron classical and Compton radius, respectively. The beam size increase should be added in quadrature to the initial beam size. We calculate the R 36 ðzÞ by tracking backwards an offmomentum particle (þ4 GeV) from the IP under the assumption that the dispersion at the IP is 0. ρðzÞ is calculated numerically from the orbit shown in Fig. 3. The estimate from this analytical formula then gives us for the solenoid field map which has the antisolenoid included.
The initial core 1 sigma beam size is about 1 nm in CLIC, which means an increase of 6.4%. If we assume that luminosity is inversely proportional to beam size, we get a luminosity loss of 6%. If we instead use the rms vertical beam size which is around 1.3 nm, we get a relative increase to the beam size of around 3.7%. However, for the peak luminosity [21] in particular, the core beam size is usually considered to be the more relevant parameter. This is an encouraging result, considering that up to 25% luminosity loss due to ISR caused by the solenoid was expected for the nominal CLIC machine [2].
Estimating luminosity loss only via beam size growth has a considerable level of uncertainty. The tails of the beam typically increase the rms beam size, while not affecting luminosity as significantly. For this reason, we always use GUINEA-PIG [23] to simulate the luminosity in our tracking studies.

IV. DETERMINISTIC APPROACH
A novel simulation approach is proposed, which separately evaluates the losses from ISR alone, before the full compensation is known. The procedure to evaluate the luminosity loss due to ISR is described in the lower part of FIG. 3. The vertical orbit of the last 7 m before the IP given in μm. The QD0 is displaced vertically in order to get the orbit to end at y ¼ 0. Fig. 4, where the tracking including the solenoid field is done using the new 4th order symplectic integrator described in the Appendix. The beam is first tracked forward without synchrotron radiation, and without the solenoid field present. This provides the optimal beam distribution at the interaction point. The ideal IP beam distribution is tracked backwards through the beam line, with the solenoid field turned on but still without synchrotron radiation. The result is a beam distribution with a perfect compensation for the coupling introduced by the solenoid field. Finally, the synchrotron radiation is turned on, and the beam is tracked forward through the beam line. The estimated luminosity is compared to a normal tracking of the beam without the solenoid field, but including ISR.
Using this approach we evaluated the simulated SiD field maps presented in Fig. 2. The loss of peak luminosity due to ISR in the detector solenoid including the antisolenoid is found to be ð4.1 AE 0.2Þ%, where the error bar is from the calculation of the luminosity in GUINEA-PIG. This result compares well to the result from the semianalytical calculation.
Without the antisolenoid, we find a luminosity loss of around 5%. 1% of additional luminosity would not alone be enough to justify the installation of an antisolenoid which significantly complicates the detector design. However, the most important purpose of the antisolenoid is to protect the permanent magnet material in the QD0 and allow it to safely reach its high gradient. The antisolenoid also makes the compensation easier, as it decouples the solenoid field from the field inside the QD0 [9]. With an alternative larger L Ã of 6 m or more [18], where the QD0 is outside the detector, the need for an antisolenoid could be reassessed.

V. FULL COMPENSATION
The main distortions responsible for the luminosity loss are vertical dispersion and y-x 0 coupling, shown in Figs. 5 and 6. From these results it is evident that the antisolenoid alone is not able to fully compensate the optical distortions caused by the main solenoid field. With the solenoid alone on the order of 1% of nominal luminosity remains, before any compensation of the beam distribution is applied. When we add the antisolenoid the luminosity increases by an order of magnitude, but is still far off acceptable performance. Other compensation methods are required in addition to the antisolenoid, in order to fully recover the luminosity.
In order to recover the residual optical distortions induced by the main solenoid field we can use knobs based on transversal sextupole displacement in addition to the antisolenoid. These linear combinations of displacements of the five sextupoles in the CLIC final focus system ideally give 5 degrees of freedom for correcting coupling and dispersion terms (vertical displacements), and 5 for correcting focusing and dispersion terms (horizontal displacements). The main couplings caused by the experimental solenoid are vertical dispersion and y-x 0 coupling, both corrected by vertical knobs. The knobs have been proven successful when applied to the tuning against magnet misalignment of the CLIC final focus, as reported in [8]. Additionally, a vertical displacement of QD0 is effective at correcting the vertical offset and dispersion at the IP. We also add horizontal displacement and roll of the QD0, for more local corrections. With the same three knobs for QF1, we have a total of 17 knobs to recover the residual luminosity loss due to the given experimental solenoid and antisolenoid design.
In our simulation, each knob is evaluated separately. We start with the QD0 knobs, then QF1, then the vertical sextupole knobs, and finally the horizontal knobs. For each knob we make a parabolic fit of the luminosity as a function of the knob value and move the magnets accordingly. This is repeated with smaller and smaller steps to make sure we are close to the optimum. We iterate over the entire algorithm a few times to make sure we have at least found a local optimum. Note that we have not taken into account any magnet imperfections or misalignments in these simulations. ISR is still activated, which means that we expect to reach a luminosity of about 96% with the antisolenoid compared to the beam line without solenoid field included, based on the result from the deterministic simulation.
In Fig. 7 we see the resulting luminosity as a function of iterations. Each dip corresponds to the iteration where the algorithm moved to a new knob. We see that most of the aberrations are corrected after the first round of QD0 and vertical sextupole knobs. Without the antisolenoid, it is not possible to obtain the same luminosity level. The number of iterations to reach optimal luminosity is about the same.
We find a luminosity loss of ð8.0 AE 1.6Þ% with the antisolenoid in these simulations. Including the error bars, this fits well with both the semianalytical estimate of 6%, and the deterministic approach which estimated 4% luminosity loss. The results give us confidence that the deterministic approach makes valid assumptions for evaluating the luminosity loss due to ISR alone.
In Fig. 8 the vertical dispersion and y-x 0 coupling can be seen with the full compensation. Only dipolar and quadrupolar terms have been used for this compensation. A check was made keeping ISR off for the entire simulation. The routine then completely canceled out the optical aberrations using only these linear elements. While one can then conclude that the solenoid itself introduces mostly linear coupling terms, the solenoid in combination with ISR can lead to nonlinear effects that may require nonlinear correctors. Indeed, the limitation of a tuning-based algorithm is the lack of knowledge about the absolute optimal luminosity.

VI. CONCLUSIONS
A novel simulation approach for estimating the irreversible luminosity loss from incoherent synchrotron radiation produced by the experimental solenoid in a high energy lepton collider has been developed. The results are compatible with the slower and more complicated simulations to find the full compensation, and consistent with a semianalytical estimate of the beam size growth. This method obtains in a deterministic way the optimal luminosity that can be achieved if the correction is perfect.
For the current SiD design for CLIC, we find that we can expect a luminosity loss due to incoherent synchrotron radiation of ð4.1 AE 0.2Þ%. This is at the optimistic end of the scale given in the conceptual design report, and is a promising result for the CLIC design effort. The antisolenoid reduces the losses due to ISR by approximately 1%, and strongly reduces the optical distortions. We have shown The results when including the antisolenoid (blue) are significantly better than the results without (green). 100% is defined as the luminosity without solenoid field. through a full simulation that the beam delivery system provides enough flexibility to correct for the optical distortions introduced by the solenoid and its overlap with the last focusing magnet.

ACKNOWLEDGMENTS
The authors would like to express their thanks to Michele Modena, Hubert Gerwig, Antonio Bartalesi, and Alexander Aloev for all their help with the field map simulations and fruitful discussions. We would also like to thank Andrea Latina and Jochem Snuverink for valuable help implementing the code into PLACET. Finally, we would like to thank Paul Anton Letnes for carefully reading through and correcting the English grammar.

APPENDIX: TRACKING ROUTINE
We have implemented a 4th order symplectic integrator in the particle tracking code PLACET [24,25], with a user defined step length. The sum of magnetic fields from beam line elements and solenoid field map is used to calculate the Lorentz force at each location of a kick. With the appropriate choices for drifts and kicks, it can be shown numerically that this integrator is in fact of 4th order [26]. This integrator allows us to track the beam through a combination of beam line elements and added field map (solenoid field), something that was not possible previously.
The new integrator has been compared to the other independent tracking routines already available in PLACET (excluding the solenoid field map) and was found to be in good agreement. The 4th order integrator has also been independently compared to a Lie tracking routine [27].
In Fig. 9 the error in the vertical position at the interaction point is shown as a function of the step length used. Each set of initial coordinates is tracked with multiple step lengths. The error is estimated as the final position at the interaction point with the given step length compared to using a much shorter step length. For each step length we then get an average error for N initial coordinates as avgðjΔyjÞ ¼ The vertical beam size at the interaction point is approximately 1 nm, so the error should be well below this value. Hence, step lengths lower than 1 cm are acceptable. The results presented in this paper are obtained using a step length of 1 mm.