Low-emittance tuning at the Cornell Electron Storage Ring Test Accelerator

Low-emittance tuning at the Cornell Electron Storage Ring Test Accelerator

J. Shanks js583@cornell.edu    D.L. Rubin    D. Sagan CLASSE, Cornell University, Ithaca, New York 14853, USA
July 27, 2019

In 2008 the Cornell Electron/Positron Storage Ring (CESR) was reconfigured from an electron/positron collider to serve as a testbed for the International Linear Collider (ILC) damping rings. One of the primary goals of the CESR Test Accelerator (CesrTA) project is to develop a fast low-emittance tuning method which scales well to large rings such as the ILC damping rings, and routinely achieves a vertical emittance of order at 2.085 GeV. This paper discusses the tuning methods developed at CesrTA to achieve low-emittance conditions. One iteration of beam-based measurement and correction requires about 10 minutes. A minimum vertical emittance of (+3.2/-3.4) (0.2) pm has been achieved at 2.085 GeV. In various configurations and beam energies the correction technique routinely achieves vertical emittance around  pm after correction, with RMS coupling . The measured vertical dispersion is dominated by beam position monitor systematics. The propagation of uncertainties in the emittance measurement is described in detail. Simulations modeling the effects of magnet misalignments, BPM errors, and emittance correction algorithm suggest the residual vertical emittance measured at the conclusion of the tuning procedure is dominated by sources other than optics errors and misalignments.

I Introduction

In 2008 the Cornell Electron Storage Ring (CESR) was reconfigured from an electron/positron collider to the CESR Test Accelerator (CesrTA) Palmer et al. (2009); Dugan et al. (2009); Rubin et al. (2009), a testbed for the International Linear Collider (ILC) damping rings ilc (2013). Parameters for the CESR storage ring are shown in Table 1. One of the primary objectives of the CesrTA program is to develop low-emittance tuning methods for the ILC damping rings.

Parameter Value Units
Circumference 768.4 m
Energy 2.085 (1.5-5.3) GeV
Lattice Type FODO
Tunes () (14.59, 9.63)
Symmetry Mirror
H / V Steerings 55/58
Quadrupoles 105
Skew Quadrupoles 27
Damping Wigglers 12
Wiggler 1.9 T
Position Monitors 100
2.7 nm
(target) 10 pm
Table 1: Parameters of the CESR electron/positron storage ring.

By far the most common tool for linear optics correction is Orbit Response Matrix (ORM) analysis, specifically Linear Optics from Closed Orbits (LOCO) Safranek (); Safranek et al. (2002). In particular, LOCO has been used as the cornerstone for corrections at both the Swiss Light Source (SLS) and the Australian Synchrotron, where vertical emittances of order 1 pm have been reported Aiba et al. (); Dowd et al. ().

However, the time required for measuring the response matrix scales linearly with the number of correctors. The Australian Synchrotron has demonstrated an acquisition rate of order 10 seconds per corrector. Assuming the ILC damping rings will be capable of the same acquisition rate, simply measuring the response matrix for all 800 steerings would take several hours, and thus response matrix analysis is deemed prohibitively slow for the ILC damping rings.

The tuning algorithm developed for CesrTA was required to be fast, and scale well to large rings such as the ILC damping rings. The correction procedure takes less than 5 minutes to acquire a full data set, where the duration time is limited by the slew rate of the superconducting RF cavities for dispersion measurements. One correction iteration (measure, compute corrections, load corrections, and remeasure) takes around 10 minutes. Data acquisition is fully parallelized, with pre-processing done on beam position monitor (BPM) modules, one per BPM. Measurement time for the CesrTA algorithm scales independently of number of BPMs, and does not depend on the number of correctors.

The algorithm may also prove useful to other storage rings. Betatron phase measurements are significantly faster than traditional response matrix analysis, allowing for less time to be spent on optics correction. The measurements may be performed using a witness bunch, exciting and measuring only a single bunch in a fully-loaded machine. Additionally, measurements such as betatron phase and coupling which resonant excitation do not require changing the machine conditions, minimizing hysteresis.

This paper describes the optics correction procedure developed at CesrTA that meets these requirements. Experimental results, with detailed propagation of uncertainties, are presented. Also discussed are simulations of the correction procedure, which have been essential to understanding measurement systematics and recognizing that the residual vertical emittance is dominated by sources other than optics errors and misalignments.

Ii Motivation for Beam-Based Emittance Tuning

For the ILC, the quantity of interest is not the projected vertical beam size but rather the emittance of the vertical-like normal mode, called the “-mode.” In principle, the beam could be intentionally coupled in the damping rings in order to reduce collective effects, and decoupled in the extraction line, so long as the -mode emittance is preserved. The decomposition into normal modes has been discussed elsewhere Wolski (2006); Qin et al. (2013); Baumgarten (2012), and therefore will not be covered here.

The primary static contributions to in a planar ring are tilted and vertically-offset quadrupoles, and rolled dipoles. Tilted quadrupoles couple horizontal and vertical motion which couples photon emission in the horizontal plane to the b-mode. Vertical quadrupole offsets and dipole rolls introduce vertical kicks, generating vertical dispersion and thus vertical emittance. Additional sources of vertical emittance include time-dependent variations associated with line voltage, ground motion, and feedback systems, which contribute kicks to the beam in various ways, and current-dependent effects such as intra-beam scattering (IBS).

Without beam-based corrections of dispersion and coupling, the vertical emittance would be limited by the quality of survey and alignment.

Survey and alignment is accomplished primarily using a Leica AT402 Absolute Laser Tracker for establishing the reference network and a Leica DNA03 digital level for elevation runs. Establishing the reference network is done via “free stationing,” with over 100 stations in the network. The reference network consists of a triplet of reference targets–attached to the inner tunnel wall, outer tunnel wall, and embedded in the floor–approximately every 8 to 10 meters in the tunnel arcs, and less regular, but on average equally dense wall, floor, and ceiling targets in the long-straight sections, flares, and main south areas. Typical combined (bundled) 1 uncertainties for the reference targets are on order 35 m for height (z). Magnets are surveyed into place to better than 100 m of measured versus ideal positions using fiducials mechanically aligned to their irons. The reference network and magnets are resurveyed regularly. Comparisons of reference target variations from survey to survey establish the uncertainties used in the analysis.

The measured distributions of surveyed quadrupole and dipole offsets and tilts for CESR are shown in Fig. 1. The root mean square (RMS) of the position and tilt measurements are summarized in Table 8, and include an estimated uncertainty in the displacement of the magnetic center from physical center of the magnet.

Simulations using random distributions of magnet, beam position monitor, and multipole errors consistent with measurements (summarized in the Appendix) have been used to study the effect of these errors on the vertical emittance. Repeating for 100 random sets of magnet errors, the resulting distributions of emittance, dispersion, and coupling yield statistical information about the probability of achieving the target emittance, and are shown in Fig. 2. The coupling is characterized using the coupling matrix, an extension of the Edwards and Teng formalism Edwards and Teng (1973) and defined in Sagan and Rubin (1999). In particular, only the out-of-phase coupling matrix element is considered, for reasons discussed in Section III.

Figure 1: Survey and alignment results for CesrTA as of December 2012 CesrTA run, compared to alignment in September 2008 at the start of the CesrTA program. Left to right: dipole roll, quadrupole tilt, and quadrupole vertical offset.
Figure 2: Resulting distributions of vertical emittance, vertical dispersion, and coupling when applying random distributions of errors at the amplitudes specified in Table 8, along with systematic and random multipoles specified in Table 9.

Without any beam-based corrections, simulations show that out of 100 random seeds, only three yielded the target vertical emittance of 10 pm; the mean vertical emittance of the 100 seeds is . It is evident that the survey and alignment techniques used are insufficient by themselves to reach the CesrTA emittance target. Some form of beam-based correction is clearly required in order to achieve and maintain low-emittance operating conditions.

Iii Measurement Techniques

Beam position monitors (BPMs) are used to collect data for most beam-based optics characterization techniques used in emittance tuning at CesrTA. CESR is instrumented with 100 button-style peak-detection BPMs. A cross-section of a typical CESR BPM is shown in Fig. 3. New electronics, developed in-house, were installed in 2009 Palmer et al. (2010). The new BPM system is capable of bunch-by-bunch, turn-by-turn readout for bunch spacings 4ns with a buffer of 300,000 bunch-turns. At each BPM, all four button channels are read out by separate controller cards, therefore channel-to-channel crosstalk is minimized. Bunch-to-bunch cross-talk is below 4% after 4ns, and is effectively zero after 50ns; there is no turn-to-turn cross-talk. Single-turn orbit reproducibility is measured to be for consecutive turns; as this is determined from beam-based measurements, it includes not only the sensitivity of the BPM, but also all contributions such as electronic stability and beam pipe vibration. Depending on the user’s request for data, some level of pre-processing is done onboard the BPM electronics before committing data to file, or the raw bunch-by-bunch turn-by-turn button signals are written directly to file.

BPMs are used to measure closed orbit, betatron phase and coupling, and dispersion. Turn-by-turn trajectory data is also used for BPM calibrations. For all beam-based measurements in low-emittance tuning at CesrTA, a single bunch of 0.8 mA = particles is used. Therefore, bunch-to-bunch effects do not contribute to emittance measurements.

Figure 3: Cross-section of a CESR BPM. Dimensions are in inches.

The closed orbit at each BPM is measured by averaging 1024 turns of turn-by-turn bunch trajectory data onboard the BPM modules. A closed orbit measurement takes roughly 5 seconds, with measurement reproducibility of around .

Dispersion measurements are performed in the usual way, by varying the RF frequency by a known amount, which changes the beam energy, and measuring the change in closed orbit. A standard dispersion measurement at CESR varies the 500 MHz cavities by  kHz (corresponding to ) and takes several minutes to acquire. The measurement time is determined by the slew rate of the RF frequency of the superconducting cavities. The measurement reproducibility is better than .

Quadrupole focusing errors are determined by measuring betatron phase advance at each BPM, using turn-by-turn data acquired while resonantly exciting the beam Sagan et al. (2000). Resonant excitation is achieved through a pair of “tune trackers,” which are stripline kickers phase-locked to the horizontal and vertical betatron tunes Meller and Palmer (2011). The tune trackers excite the beam to amplitudes of several millimeters. Phase and amplitude data are extracted from 40,960 consecutive turns by a processor onboard the BPM module for each button. The button-by-button phase and amplitude are post-processed into horizontal and vertical phase, the out-of-phase component of the coupling matrix , and the two in-phase components of the coupling matrix . All of the above information is processed from one measurement of the machine. Betatron phase data for all 100 BPMs is collected and analyzed in 10 seconds. Reproducibility of betatron phase measurements is of order .

When characterizing coupling, only the coupling matrix element is used, and the other two measured components ( and ) are neglected. It is not possible to measure . represents the out-of-phase propagation of coupling (a sine-like vertical motion at the horizontal tune, when the horizontal is excited with a cosine-like signal), and it can be shown that is insensitive to BPM rotations. Second, two independent measurements of can be made simultaneously, from excitations of the horizontal and vertical modes, adding redundancy in the measurement. Third, because the (normalized) coupling matrix elements mix from one BPM to the next, measuring and correcting globally is sufficient to correct the entire coupling matrix.

For clarity it is perhaps worth comparing the betatron phase and coupling technique with ORM measurement and analysis. The response matrix is established by measurement of the closed orbit (position) at each BPM in response to excitation of the distributed steering corrector magnets. From the set of measured closed orbits (two for each steering), the linear transfer matrix (betatron phase advance and coupling) from one BPM to the next may be reconstructed.

In the phase and coupling measurement as implemented for CESR, the transfer matrices are similarly reconstructed from measured trajectories. The turn by turn capability of the beam position monitors is exploited to significantly reduce the measurement time. Rather than drive DC correctors to generate a distribution of trajectories, the beam is resonantly excited at the betatron tunes by a single source at a fixed location. As the tunes are non-integer, the 40,960 consecutive trajectories smoothly sample phase space, and are acquired in a fraction of a second. The trajectories could in principle be analyzed using ORM techniques. Alternatively (and equivalently), we extract betatron phase and amplitude and coupling information with the help of pre-processing in each of the BPM modules. The objective for the emittance correction method at CesrTA is a technique suitable for a ring with very large circumference, like the ILC damping ring, therefore the betatron phase and coupling technique is favored.

Iv BPM Calibrations

In order to ensure that measurements reflect actual machine conditions, BPMs must be well-calibrated. The primary characteristics to consider are: button-by-button timing, button-to-button relative gains, BPM tilts, and BPM-to-quadrupole transverse offsets.

Many modern lightsource BPMs take four signals into one controller that pre-processes the raw signals into horizontal and vertical data. CESR BPMs have four separate controller cards, one for each button, which read out independently. This allows for greater flexibility in measurements and post-processing, however some characteristics such as timing and gains must be calibrated on each of the four button channels rather than once per BPM.

Each of the required calibrations are now discussed in the order of implementation.

iv.1 BPM Timing

Each controller card has independent timing, therefore every button on every BPM must be timed in separately. A mistimed channel results in sampling the bunch passage off-peak, which reduces the observed signal amplitude for that button.

The time-in procedure consists of sampling the temporal profile of a bunch passage at a resolution of 10 ps and fitting to determine the peak. The process converges in than one minute for all four buttons on all 100 BPMs, with less than 10 ps drift over a period of four hours.

iv.2 Button-to-Button Relative Gains

Differential response of the four BPM buttons due to variations in relative electronic gain will introduce a systematic measurement error. Measurements that depend mostly on position, such as orbit, dispersion, and the in-phase components of the coupling matrix , are sensitive to relative button gains. Measurements using relative phase, such as betatron phase advance and the out-of-phase coupling matrix element are largely insensitive to gain errors.

The method of gain mapping used at CesrTA was developed by Rubin et al. at Cornell Rubin et al. (2010), and is based on a second-order expansion of the button signal response. The method utilizes turn-by-turn data, therefore data acquisition is fast, on the order of several seconds to collect data for characterizing all 100 BPMs.

The analysis relies on the fact that a linear relation exists between two combinations of the four button signals. For turns of turn-by-turn trajectory data there are button measurements at each BPM. There are only four unknowns, namely the button gains, and the system is overconstrained for orbits; typically 1024 turns are used. Data acquisition takes about 10 seconds, and the fitting process takes less than a minute to determine all four button gains on all 100 BPMs.

All gain calibration techniques for peak-detection-style BPMs are sensitive to timing errors. This method is insensitive to detector rotation or offset, as the method uses raw button signals across a large cross-section of the BPM, and does not rely on distinguishing between horizontal and vertical modes.

Typical BPM gain variations before correction are of order 5%, and are calibrated with a reproducibility of a few tenths of a percent.

iv.3 BPM Electronic Centering

A relative offset between the electronic center of a BPM and the magnetic center of the nearest quadrupole will appear in measurements as an offset in the quad. If the relative offset is not calibrated, steering the beam to the electronic center of the BPM will result in kicks from the quadrupole, generating dispersion. To minimize vertical dispersion (and thus the emittance) generated during orbit correction, the relative offset between the electronic center of a BPM and the magnetic center of the nearest quadrupole must be measured.

The method used at CESR for determining the BPM-to-quadrupole transverse offset is based on a common technique where the beam is steered such that a change in the quadrupole strength results in no change in orbit Chung et al. (1995); Tecker et al. (1997). The method employed at CESR has the additional benefit that the change in quadrupole strength is determined using betatron phase measurements Shanks et al. (2010). By measuring the difference in phase before and after the quadrupole change, the change in quadrupole strength can be more accurately determined than using the change in quadrupole current, which is susceptible to hysteresis. Therefore, fewer iterations are required on each BPM/quadrupole pair to achieve convergence.

Typical BPM-to-quadrupole offset measurements are around RMS in both horizontal and vertical, with a short-term reproducibility of order m and long-term drift of about over the course of a three-week machine studies period.

BPM-to-quadrupole relative centering will only affect orbit measurements and turn-by-turn trajectory data. Dispersion measurements are a difference of two closed orbits, therefore absolute offsets do not affect the measurement. Betatron phase and coupling measurements are computed button-by-button, therefore transverse offsets will not affect the measurement.

iv.4 BPM Tilt Calibration

If a BPM is rotated, a horizontal orbit perturbation will indicate a vertical offset. This becomes particularly significant when measuring the dispersion, as the average horizontal dispersion in CESR is large, on the order of a meter. The lowest vertical dispersion measured without BPM tilt corrections is around 12 mm RMS. The measurements cannot be fit with a model dispersion function generated by corrector magnets or magnet misalignments. Furthermore, simulations have shown that an RMS of 12 mm of actual vertical dispersion corresponds to vertical emittance, significantly larger than the emittance determined from vertical beam size measurements. This implies the measurement is dominated by systematic errors, such as uncalibrated BPM tilts. If the BPM tilts are uncorrelated with the horizontal dispersion, this constrains the distribution of BPM tilts to have a maximum RMS of 12 mrad.

BPM tilts can in principle be determined by linearly fitting turn-by-turn trajectory data for a well-decoupled beam which is resonantly excited in the horizontal mode. Residual in-phase coupling will also rotate the beam in x-y space, which introduces a lower bound on the ability to resolve BPM tilts using this method to around 5 mrad. To date, applying the BPM tilts to dispersion data does not improve the ability to fit the vertical dispersion function. Several alternative methods for measuring BPM tilts have been explored, all yielding different calibrations, and none improving the ability to fit the vertical dispersion. As such, the tilt calibrations have not been utilized during any of the corrections described in this paper, and remain the most significant known systematic in optics corrections.

V Beamsize Instrumentation - xBSM

The primary method of determining the effectiveness of vertical emittance tuning is direct observation of the vertical beam size, from which the emittance can be inferred. CESR is instrumented with two x-ray beam size monitors (xBSM), one for each species Rider et al. (2011, 2012).

The xBSMs are one-dimensional 32-diode arrays with pixel pitch. The instruments are capable of bunch-by-bunch, turn-by-turn measurements with a buffer of 250,000 bunch-turns. Dynamic range for the instruments span beam currents at the standard CesrTA operating energy of 2.085GeV.

When characterizing low-emittance conditions, the beam is typically imaged using a horizontal slit formed by two tungsten blades, which act as a one-dimensional pinhole. Beam size is determined by fitting to the beam profile over 1024 turns on a turn-by-turn basis, then averaging. In this way any effect of turn-by-turn beam centroid or x-ray optics motion is removed from the measured beam size. The fitting procedure takes into account the point-response function (prf) of the imaging device (in this case, the 1-dimensional pinhole), including the effects of the finite opening angle, depth-of-field, energy spectrum, diffraction off the tungsten blades, surface roughness of the tungsten blades, and detector response. A detailed analysis of the fitting procedure is available in Alexander et al. (2014). In practice, the resolution limit when using the pinhole optics is around . The vertical beta function at the xBSM source point is approximately 40 m, and the xBSM optics provide a magnification of approximately 2.2. Therefore, the xBSMs are able to resolve the vertical emittance down to .

Vi Low-Emittance Tuning

The low-emittance tuning procedure developed at CesrTA takes advantage of the fact that all magnets are independently powered, and all BPMs are capable of betatron phase and coupling measurements through turn-by-turn acquisition. The procedure is as follows:

  1. Measure the closed orbit and correct to a reference orbit (which aligns the beam with the xBSM beamline) using all 55 horizontal and 58 vertical steering correctors.

  2. Measure the betatron phase, transverse coupling (), and horizontal dispersion. Fit the model lattice to the measurement using all 100 quadrupoles and 27 skew quadrupole correctors, and load the computed corrections.

  3. Remeasure the closed orbit, transverse coupling, and vertical dispersion. Fit the model lattice to all machine data simultaneously using all vertical steerings and skew quadrupoles, and load the fitted corrector changes.

The turnaround time for one full set of corrections is roughly ten minutes. It is standard procedure when first recovering conditions to save magnet settings after achieving low emittance, run the machine through a well-defined hysteresis loop, re-load the previously saved conditions, and repeat the emittance tuning procedure to apply minor corrections and ensure the desired conditions are reproducible.

Lattice corrections are determined by a minimization where a machine model is fit to measurements of the lattice functions, with a merit function defined as Sagan et al. (2000):


where is the datum (for example, the vertical orbit at a BPM), is the variable (such as a corrector strength), and are user-defined weights. The merit function is minimized by adjusting corrector magnets in the model such that the model reproduces the measurements. The negative of the machine model corrector strengths are then loaded into the machine to compensate for optics errors.

Beam-based characterization of the machine after a typical low-emittance correction is shown in Table 2. The discrepancy between the model that best fits those measurements and the design demonstrates the effectiveness of the correction.

Measurement RMS (Data) RMS (Model) Units
253 110 m
0.3 0.3 deg
13 5 mm
0.004 0.003
Table 2: Typical levels of correction for optics measurement after the full emittance tuning procedure. Measurements were taken at 0.8 mA (). RMS deviations from the design are reported for both the machine measurement and a model of the machine which is fit to the measurements. Beta beat is computed from fitting phase data.

The -mode emittance is determined from measurements of the beam size and local optics at the beam size monitor source point:


where is a parameter related to the coupling matrix, and is approximately unity when coupling is small. is the contribution from the -mode to the vertical beam size, which is computed from the total vertical beam size :


where , , and are the contributions to the vertical beam size from the horizontal emittance, vertical emittance, and vertical dispersion. is the horizontal-like normal mode. and are in-phase and out-of-phase components of the coupling matrix, respectively, and are directly measured at BPMs adjacent to the source point. Equations 3-5 are used to determine the component of the vertical beam size due to -mode emittance, , which is then used in Eqn. 2.

, , , and are determined by fitting a model of the accelerator to measurements. The beam size at the source point is calculated from the measured image at the xBSM , accounting for the magnification, energy spectrum, and point-source response of the optic/detector system.

Statistical and systematic errors associated with measurements of vertical emittance with the xBSM are outlined in Shanks et al. (2011), and include contributions from: turn-by-turn beamsize fitting uncertainty; turn-by-turn beamsize fluctuation; uncertainty in pinhole size; uncertainty in functions; uncertainty in longitudinal location of the x-ray source point; and uncertainty in dispersion at the source point. The uncertainties propagate as follows:




and and refer to the systematic and statistical uncertainties, respectively. The individual terms are computed by varying the terms in the emittance calculation by their uncertainties .

Using the above tuning method, and propagating errors according to Equations 78, the vertical emittances achieved at CesrTA are reported in Table 3.

Energy [GeV] Species  (pm)  (pm)  (pm)
2.085 (03/2013) 10.3
2.085 (03/2013) 13.3
2.305 (12/2012) 10.0
2.553 (03/2013) 10.2
Table 3: Lowest-achieved vertical emittance at CesrTA in a variety of energies. Electron conditions are only reported for the April 2013 CesrTA run at 2.085 GeV.

These results warrant a few comments. First, all four reported measurements are within of the target  pm. It is also interesting to note that the minimum-achieved emittances are independent of energy. All four measurements are within of each other, with only the single electron measurement standing out.

The component of the observed vertical beam size due to local coupling () only introduces 0.5 m in quadrature. The observed beam size is of order 20 m, and as such, the contribution of local coupling is insignificant.

It is also worth reminding the reader that the minimum measured RMS of 12 mm corresponds to a vertical emittance of 20-30 pm, much larger than what has been measured. It is clear that the dispersion measurement is dominated by systematics, and in particular, BPM tilts must be better understood. However, correcting the vertical orbit and transverse coupling indirectly reduces the vertical dispersion. The procedure therefore results in a vertical dispersion below the present resolution of the measurement.

Several alternative LET tuning methods have been explored, including Orbit Response Matrix (ORM) analysis Shanks et al. (2009) and normal-mode analysis Wolski et al. (2011). To date, no method has proven to be faster or yield consistently better results than the three-stage correction algorithm based on betatron phase and coupling measurements discussed here.

Vii LET Simulations

To better characterize what factors are limiting emittance corrections, software has been developed to evaluate the contributions of misalignments, BPM measurement errors, and choice of correction procedure. The program, ring_ma2, uses the Bmad accelerator code library Sagan (2006), and does the following:

  1. Assigns random misalignments and BPM errors with user-defined amplitudes to the ideal lattice to create a realistic machine model.

  2. Simulates beam based measurements of optics functions including the effects of BPM measurement errors.

  3. Computes and applies corrections for each iteration based on the simulated measurements.

  4. After each correction iteration, it records the effectiveness of the correction in terms of emittances and optics functions.

The entire procedure is repeated typically 100 times in order to generate statistics for analysis. The simulation is approached from a statistical perspective for three reasons. First, magnet positions continually drift, making it difficult to know the exact set of misalignments in the ring on any given day. Second, the precise distributions of magnetic centering or BPM measurement errors are not known, mandating that their distributions be approached from a statistical perspective. Third, by framing the analysis in terms of statistical probability of achieving the required emittance, the characterization process may be extrapolated to new machines which are not yet built using only the knowledge of survey and alignment tolerances.

When discussing the results of statistical analysis the 95% confidence levels (CL) are presented. That is, after applying the full optics correction procedure 95% of simulated lattices, each with a randomly chosen distribution of misalignments and measurement errors, achieve a vertical emittance below the 95%CL. The simulation is believed to be sufficiently complete such that it is very unlikely that the contribution of the static optics to the vertical emittance is greater than this number.

In this section the method for simulating optics measurements is discussed, including how BPM measurement errors and guide field magnet errors are modeled. Results of simulations based on input parameters representing the physical accelerator are given.

vii.1 Model Lattice with Errors

Bmad allows for introducing strength errors (including systematic and random multipole errors) and alignment errors (such as offset, roll, and pitch) to any lattice element. Magnet strength errors scale with the absolute strength of the magnet. Alignment errors are treated as additive errors, and are applied directly without scaling.

ring_ma2 also models BPM measurement errors, which are discussed in detail in Section VII.3.

vii.2 Simulated Measurements

All simulated measurements are modeled as realistically as possible. For closed orbit measurements this involves recording 1024 turns of trajectory data, including the effects of BPM measurement errors on every turn, and averaging the results. Dispersion is simulated as a difference of two closed orbits, varying the RF frequency in-between.

For phase and coupling measurements, a particle is resonantly excited using a simulated phase-locked tune tracker and allowed to equilibrate by tracking for several damping times ( turns). After the particle trajectory has equilibrated, 40,960 turns of raw BPM button data are recorded at every BPM. The data is then processed with the same code used for processing CESR phase and coupling data.

A comparison of lattice parameters derived from simulated measurements in an ideal lattice and those computed directly are summarized in Table 4 for each measurement type, and presumably represent a fundamental lower limit to the resolution of each measurement technique given no errors in the BPM measurements. Simulated measurements have differing levels of agreement for horizontal and vertical, which can be attributed to energy loss from stochastic radiation emission, leading to a “sawtooth” horizontal dispersive orbit between the two pairs of RF cavities on opposite sides of the ring. This effect is not seen in the vertical as there is no vertical dispersion in the design lattice.

Measurement RMS (Simulated - Bmad) Units
Closed Orbit 0.70, m
0.75, mm
0.1, 0.05 deg
Table 4: RMS difference between simulated measurements and Bmad-calculated values, neglecting any BPM measurement errors.

vii.3 BPM Errors

To generate simulated measurements as realistically as possible, BPM measurement errors must be taken into account. The two classes of BPM errors modeled in ring_ma2 are BPM misalignments (offsets, tilts, and shear) and button-by-button effects (button gain, timing, and electronic noise). Each class of errors will affect the measurement differently. All simulated measurements presented include the effects of all listed BPM measurement errors.

vii.3.1 BPM Misalignments

Errors in BPM misalignments (offsets and tilts) are applied in the following way:


where are the coordinates with BPM misalignments applied, is the rotation matrix for angle , and are the horizontal and vertical offset between the BPM and nearest quadrupole.

vii.3.2 Button Effects: Gain, Timing, and Reproducibility

Timing errors, gain variations, and turn-by-turn resolution affect individual button signals. Modeling their effects requires an accurate method for converting from coordinates to button signals , applying errors, and converting back to coordinates.

All button-by-button errors of these classes are handled through use of a nonlinear interpolation grid which converts coordinates to button signals. Button-by-button errors are applied to the individual channels, and the final “measured” coordinates are determined by the best fit to the set of new button signals using the same interpolation grid Helms and Hoffstaetter (2005). The nonlinear map used in these studies is for a BPM with a “CESR geometry” (see Fig. 3).

Including effects from button-to-button gain errors, timing errors, and measurement reproducibility, the four observed button signals at each BPM are:


In Equation 10 are defined to be the button signals determined through the interpolation grid for the coordinates from Equation 9. is the gain error on button , and is an effective gain error for button arising from the timing error:


where the constants are empirically determined. Note that because CESR BPMs are timed to the peak signal of a bunch passage, any timing error will decrease the button signal. This method of modeling the timing error also allows the BPM model to account for synchrotron motion, thus modulating the timings on all four buttons on a turn-by-turn basis.

BPM position measurement reproducibility is dominated by electronic noise arising from the digitization and amplification of an analog signal on each of the four controllers, and is modeled in Equation 10 as an additive error on each of the four button signals. The amplitude of the button-by-button reproducibility is set by determining the change in a single button signal consistent with changing the observed orbit by the desired reproducibility (for example m).

The top and bottom CESR BPM button blocks are welded separately to the vacuum chamber. There is then the possibility of a relative misalignment of the two blocks. In order to estimate the effect of this misalignment in simulation, upper and lower button signals are determined by offsetting the BPM in opposite directions.

vii.4 Simulation Results

Amplitudes for misalignments and BPM errors in the simulation are summarized in Tables 7, 8, and 9, and are determined either from directly-measured values or inferred from machine measurements. Offsets of quadrupoles and sextupoles include measured alignment levels along with 100 m added in quadrature to account for the estimated uncertainty in the offset of magnetic center with respect to geometric center of these elements.

The emittance correction procedure used in the simulation is identical to that used on the actual machine, outlined in Section VI. Results from ring_ma2 are shown in Fig. 4, and summarized in Table 5.

Figure 4: Results from ring_ma2, using misalignments, BPM measurement errors, and multipoles stated in Tables 7, 8, and 9, plotted before correction (red), and after first, second, and third stage of emittance correction (blue, green, and black, respectively). The dashed black line indicates typical measured values in CESR after low-emittance correction.
Measurement Initial Iter 1 Iter 2 Iter 3 Units
7.7 1.6 0.1 0.1 deg
42.6 18.7 18.7 15.4 mm
40.1 13.9 12.2 5.0 mm
6.3 3.2 0.34 0.24
255.8 33.0 27.5 4.1 pm
Table 5: 95% confidence level (CL) correction levels after each correction iteration. All values except include observational effects from BPM measurement errors. Details of the correction iterations are discussed in Section VI.

After correction, 95% of seeds achieved a vertical emittance below 4.1 pm, which is significantly smaller than the minimum measured vertical emittance of 10.3 (+3.2/-3.4) (0.2) pm at 2.085 GeV. It is clear that the simulation does not account for more than half of the measured vertical emittance.

BPM tilts are the single most significant contribution to the vertical emittance in the simulation, and dominate the simulated measurement. Considering the simulation under-estimates the measured vertical emittance, one could envision adjusting the simulation to reduce the amplitude of BPM tilts and increase magnet misalignments in order to increase the vertical emittance to levels measured in the actual machine while holding the RMS constant. However, the required change in alignment to generate the measured vertical emittance is much larger than the measured uncertainty in the alignment procedure.

Viii Diagnosis of Emittance Dilution

The two primary mechanisms for the static optics to contribute to vertical emittance are vertical dispersion and horizontal-to-vertical coupling. The measured vertical dispersion in Table 2 and the minimum measured at CesrTA () are within the distributions from the simulation, indicating that the model is realistic. Increasing the coupling in simulated lattices such that the RMS is consistent with the measurement in Table 2 introduces less than 1 pm of vertical emittance.

Additionally, significant efforts have been made to ensure that all sources of uncertainty in the emittance measurement are accounted for. The discrepancy is therefore not attributed to emittance measurement errors. This implies there are significant sources of vertical emittance that are not included in the model or ring_ma2 simulation. Potential sources of vertical emittance are now considered.

viii.1 Emittance Dilution from RF

Random RF voltage and phase jitter may contribute to emittance dilution. There are four superconducting RF cavities in CESR, split into two pairs. Each pair is powered by a single power supply. Turn-by-turn beam size was recorded while varying the total RF voltage and number of RF cavities powered. The results are summarized in Table 6. It should be noted that the studies summarized in this section were taken while one of the two West RF cavities was disabled, therefore only three RF cavities were used (one in the West, and two in the East). Nominal total RF voltage was 4.8 MV, distributed approximately evenly among the three cavities.

Total RF (MV) East RF West RF  (pm)
4.8 On On 11.5
1.7 On On 11.2
1.7 Off On 12.5
1.7 On Off 10.8
Table 6: Summary of beam stability tests at CesrTA. The measurements were conducted in April 2013, for a single bunch of positrons at 0.7-0.85 mA.

A small reduction in beam size was observed when reducing the total RF voltage from 4.8 MV to 1.7 MV, corresponding to a reduction in observed vertical emittance of 0.3 pm. The statistical uncertainty for the lowest-measured emittance is 0.2 pm. Note that although the systematic uncertainty is an order of magnitude larger, it represents a global uncertainty where all measurements would be affected uniformly by any change in the understanding of the beamsize measurement system.

A further reduction is seen when the single West RF cavity is powered down and detuned, such that only the two East RF cavities are running; the emittance increased slightly when running only on the W1 RF cavity. This indicates that the RF system is contributing to the vertical emittance, although the mechanism is not known at this time. The East and West RF cavity pairs run on separate power supplies; one hypothesis is that the West RF power supply is less stable than the East, thereby introducing vertical emittance through modulation of the RF voltage. Alternatively, by running a single cavity at a higher voltage, the amplitude of voltage jitter is also increased, potentially increasing the contribution to the emittance. The RF system in CESR is superconducting, therefore a direct examination of the alignment requires the nontrivial process of warming the cavities and opening the cryostats.

viii.2 Collective Effects

The CesrTA emittance target of 10 pm is for a “zero-current” beam; that is, neglecting any collective effects. Effects considered here include: electron cloud, fast-ion instability, intra-beam scattering, and wakefields.

Electron cloud and fast-ion instability typically require a train of 30 bunches with 0.5 mA/bunch or more in order for emittance dilution to occur, and the emittance blow-up takes place around bunch 10-15 in the train Billing et al. (2013). A single bunch is not capable of generating sufficient cloud or ion density to cause emittance dilution.

Extensive measurements and simulations of intra-beam scattering (IBS) at CesrTA indicate that the vertical emittance is largely insensitive to IBS effects at currents  mA/bunch, where the measurements reported here were taken Ehrlichman (2013); Ehrlichman et al. (2012). The mechanism through which IBS increases vertical emittance depends on either transverse-to-longitudinal scattering in regions with dispersion or transverse-to-transverse scattering in regions with coupling, such that the vertical-mode action of the particle changes. Vertical dispersion and coupling are measured to be globally well-corrected, and are well below levels required for IBS to contribute to vertical emittance dilution.

Wakefields would tend to increase the emittance linearly with current. By measuring the dependence of the vertical beam size on current, it may be possible to determine whether wakefields are contributing to the emittance at the nominal 0.8 mA/bunch used for the emittance measurements presented in Table 3. However, at such low current, photons are sparse and the turn-by-turn fitting procedure is no longer reliable. Instead, the turn-by-turn images must be averaged first in order to improve signal-to-noise, then fit as a single image. This has the disadvantages of incorporating a small amount of turn-by-turn beam motion and increasing the statistical uncertainty in the vertical emittance measurement.

Figure 5 shows the emittance calculated from a series of vertical beam size measurements from the xBSM, taken sequentially as the current was decreased from 1.1 mA to around 0.05 mA, and processed as described above.

Figure 5: Vertical emittance for a single bunch of positrons as a function of bunch current, from April 2013 CesrTA machine studies. Plotted error bars are systematic (red) and statistical (blue). The dashed horizontal line indicates the 10 pm zero-current vertical emittance target for CesrTA.

The rolloff of beam size at very low current ( mA) is likely due to lack of sufficient flux on the beam size monitor. Moreover, the emittance does not depend linearly on beam current. As such, there is no support for current-dependent effects diluting the vertical emittance at low-current.

Ix Conclusions

A low-emittance tuning procedure has been developed at CesrTA, based on betatron phase and coupling measurements using resonant excitation and turn-by-turn capable BPM. The tuning procedure has a fast turnaround, where one round of optics correction takes about ten minutes, and has yielded a single-bunch vertical emittance of (+3.2/-3.4) (0.2) pm with a single bunch of positrons with at 2.085 GeV. The correction procedure routinely achieves  pm in a variety of machine conditions at energies ranging from 2.085-2.5 GeV.

The tuning procedure developed at CesrTA is significantly faster than response matrix analysis. The method scales independently of number of BPMs or correctors, thus for large machines the CesrTA procedure will be proportionally faster than response matrix analysis.

The lack of energy dependence for the minimum-achieved vertical emittance may yield information regarding sources of emittance dilution. Further studies, including measuring the emittance at several energies during a single machine studies period, will be necessary before a conclusive statement may be made.

Collective effects do not appear to contribute to emittance dilution for a single bunch at 0.8 mA. The RF system on the other hand clearly does affect the emittance, and further investigations are necessary.

Although misalignments do not appear to be the most significant contribution to the emittance, any improvement in alignment or optics correction will likely result in a small reduction in the emittance, as contributions to the emittance add linearly. In particular, BPM tilts remain a significant outstanding issue which limits the understanding of . Simulations suggest that a reduction in RMS BPM tilt from 12 mrad to 5 mrad could reduce the contribution of the static optics to the vertical emittance by 50%. Alternative BPM tilt fitting techniques are under development.

The authors wish to thank the CESR operations and instrumentation groups, whose support was indispensable in our efforts to achieve well-corrected conditions. This work was supported by the National Science Foundation grant PHY-1002467 and the Department of Energy grant DE-SC0006505.

Appendix A Appendix: Errors for ring_ma2 Simulations

Table 8 shows the misalignments and errors used in CesrTA ring_ma2 studies. Offsets of quadrupoles and sextupoles include measured alignment levels along with 100 m added in quadrature to account for the estimated uncertainty in the offset of magnetic center with respect to geometric center of these elements.

Systematic multipoles are included for sextupoles which have vertical steering or skew quadrupole trim windings. These multipoles are computed using field modeling software, and are scaled to a measurement radius of 20 mm. There is a known random skew quadrupole component to the damping wiggler fields Crittenden et al. (2005), due to manufacturing tolerances in the radii of the pole windings, which is also included. Multipoles used in this study are summarized in Table 9, and use the following convention (summarized in the Bmad manual Sagan (2013)):


where and are the normal and skew multipoles, respectively. The values in the table are normalized by , where is the order of the primary field ( for quadrupole, for sextupole, etc.).

Error Applied RMS Units
Reproducibility 10 m
Tilt 12 mrad
Gains 0.5% %
Timing 10 ps
Offset () 170 m
Horizontal Shear m
Table 7: BPM errors introduced into model CesrTA lattice for ring_ma2 studies.
Element Class Error RMS Units
Dipole Offset 0.9 mm
Offset 2.0 mm
Offset 2.3 mm
Roll 144 rad
Pitch 600 rad
Pitch 300 rad
Quadrupole Offset 335 m
Offset 40.3 m
Magnetic Offset 100 m
Offset 5.2 mm
Tilt 148 rad
Pitch 1100 rad
Pitch 62 rad
k1 0.1% %
Sextupole Offset 280 m
Offset 280 m
Magnetic Offset 100 m
Offset 5.2 mm
Tilt 200 rad
Pitch 1200 rad
Pitch 800 rad
k2 0.1% %
Wiggler Offset 1 mm
Offset 250 m
Offset 500 m
Tilt 300 rad
Pitch 200 rad
Pitch 250 rad
Table 8: Misalignments and errors introduced into model CesrTA lattice for ring_ma2 studies. All parameters are determined either from machine measurements or survey.
Element Class Multipole Value
Sextupole with a3
Vert. Steering Trim a5
Sextupole with a4
Skew Quad Trim a6
Wiggler a1
Table 9: Multipoles used in ring_ma2 studies of CesrTA lattice. Sextupole multipoles are systematic and therefore identical at all sextupoles, whereas the wiggler multipole is random; the number quoted for wiggler is therefore the RMS of the applied distribution.


Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description