Lowemittance tuning at the Cornell Electron Storage Ring Test Accelerator
Abstract
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 lowemittance 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 lowemittance conditions. One iteration of beambased 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 lowemittance tuning methods for the ILC damping rings.
Parameter  Value  Units 

Circumference  768.4  m 
Energy  2.085 (1.55.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 
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 preprocessing 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 fullyloaded 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 BeamBased Emittance Tuning
For the ILC, the quantity of interest is not the projected vertical beam size but rather the emittance of the verticallike 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 verticallyoffset quadrupoles, and rolled dipoles. Tilted quadrupoles couple horizontal and vertical motion which couples photon emission in the horizontal plane to the bmode. Vertical quadrupole offsets and dipole rolls introduce vertical kicks, generating vertical dispersion and thus vertical emittance. Additional sources of vertical emittance include timedependent variations associated with line voltage, ground motion, and feedback systems, which contribute kicks to the beam in various ways, and currentdependent effects such as intrabeam scattering (IBS).
Without beambased 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 longstraight 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 outofphase coupling matrix element is considered, for reasons discussed in Section III.
Without any beambased 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 beambased correction is clearly required in order to achieve and maintain lowemittance operating conditions.
Iii Measurement Techniques
Beam position monitors (BPMs) are used to collect data for most beambased optics characterization techniques used in emittance tuning at CesrTA. CESR is instrumented with 100 buttonstyle peakdetection BPMs. A crosssection of a typical CESR BPM is shown in Fig. 3. New electronics, developed inhouse, were installed in 2009 Palmer et al. (2010). The new BPM system is capable of bunchbybunch, turnbyturn readout for bunch spacings 4ns with a buffer of 300,000 bunchturns. At each BPM, all four button channels are read out by separate controller cards, therefore channeltochannel crosstalk is minimized. Bunchtobunch crosstalk is below 4% after 4ns, and is effectively zero after 50ns; there is no turntoturn crosstalk. Singleturn orbit reproducibility is measured to be for consecutive turns; as this is determined from beambased 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 preprocessing is done onboard the BPM electronics before committing data to file, or the raw bunchbybunch turnbyturn button signals are written directly to file.
BPMs are used to measure closed orbit, betatron phase and coupling, and dispersion. Turnbyturn trajectory data is also used for BPM calibrations. For all beambased measurements in lowemittance tuning at CesrTA, a single bunch of 0.8 mA = particles is used. Therefore, bunchtobunch effects do not contribute to emittance measurements.
The closed orbit at each BPM is measured by averaging 1024 turns of turnbyturn 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 turnbyturn 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 phaselocked 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 buttonbybutton phase and amplitude are postprocessed into horizontal and vertical phase, the outofphase component of the coupling matrix , and the two inphase 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 outofphase propagation of coupling (a sinelike vertical motion at the horizontal tune, when the horizontal is excited with a cosinelike 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 noninteger, 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 preprocessing 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 wellcalibrated. The primary characteristics to consider are: buttonbybutton timing, buttontobutton relative gains, BPM tilts, and BPMtoquadrupole transverse offsets.
Many modern lightsource BPMs take four signals into one controller that preprocesses 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 postprocessing, 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 offpeak, which reduces the observed signal amplitude for that button.
The timein 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 ButtontoButton 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 inphase components of the coupling matrix , are sensitive to relative button gains. Measurements using relative phase, such as betatron phase advance and the outofphase 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 secondorder expansion of the button signal response. The method utilizes turnbyturn 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 turnbyturn 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 peakdetectionstyle 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 crosssection 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 BPMtoquadrupole 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 BPMtoquadrupole offset measurements are around RMS in both horizontal and vertical, with a shortterm reproducibility of order m and longterm drift of about over the course of a threeweek machine studies period.
BPMtoquadrupole relative centering will only affect orbit measurements and turnbyturn 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 buttonbybutton, 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 turnbyturn trajectory data for a welldecoupled beam which is resonantly excited in the horizontal mode. Residual inphase coupling will also rotate the beam in xy 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 xray beam size monitors (xBSM), one for each species Rider et al. (2011, 2012).
The xBSMs are onedimensional 32diode arrays with pixel pitch. The instruments are capable of bunchbybunch, turnbyturn measurements with a buffer of 250,000 bunchturns. Dynamic range for the instruments span beam currents at the standard CesrTA operating energy of 2.085GeV.
When characterizing lowemittance conditions, the beam is typically imaged using a horizontal slit formed by two tungsten blades, which act as a onedimensional pinhole. Beam size is determined by fitting to the beam profile over 1024 turns on a turnbyturn basis, then averaging. In this way any effect of turnbyturn beam centroid or xray optics motion is removed from the measured beam size. The fitting procedure takes into account the pointresponse function (prf) of the imaging device (in this case, the 1dimensional pinhole), including the effects of the finite opening angle, depthoffield, 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 LowEmittance Tuning
The lowemittance 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 turnbyturn acquisition. The procedure is as follows:

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.

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.

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 welldefined hysteresis loop, reload 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):
(1) 
where is the datum (for example, the vertical orbit at a BPM), is the variable (such as a corrector strength), and are userdefined 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.
Beambased characterization of the machine after a typical lowemittance 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  – 
The mode emittance is determined from measurements of the beam size and local optics at the beam size monitor source point:
(2) 
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 :
(3)  
(4)  
(5) 
where , , and are the contributions to the vertical beam size from the horizontal emittance, vertical emittance, and vertical dispersion. is the horizontallike normal mode. and are inphase and outofphase components of the coupling matrix, respectively, and are directly measured at BPMs adjacent to the source point. Equations 35 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 pointsource 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: turnbyturn beamsize fitting uncertainty; turnbyturn beamsize fluctuation; uncertainty in pinhole size; uncertainty in functions; uncertainty in longitudinal location of the xray source point; and uncertainty in dispersion at the source point. The uncertainties propagate as follows:
(6)  
(7)  
where
(8)  
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 7–8, 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 
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 minimumachieved 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 2030 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 normalmode analysis Wolski et al. (2011). To date, no method has proven to be faster or yield consistently better results than the threestage 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:

Assigns random misalignments and BPM errors with userdefined amplitudes to the ideal lattice to create a realistic machine model.

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

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

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 inbetween.
For phase and coupling measurements, a particle is resonantly excited using a simulated phaselocked 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  
– 
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 buttonbybutton 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:
(9) 
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 turnbyturn 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 buttonbybutton errors of these classes are handled through use of a nonlinear interpolation grid which converts coordinates to button signals. Buttonbybutton 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 buttontobutton gain errors, timing errors, and measurement reproducibility, the four observed button signals at each BPM are:
(10) 
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:
(11) 
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 turnbyturn 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 buttonbybutton 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 directlymeasured 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.
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 
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 underestimates 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 horizontaltovertical 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. Turnbyturn 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 
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 lowestmeasured 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 “zerocurrent” beam; that is, neglecting any collective effects. Effects considered here include: electron cloud, fastion instability, intrabeam scattering, and wakefields.
Electron cloud and fastion 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 blowup takes place around bunch 1015 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 intrabeam 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 transversetolongitudinal scattering in regions with dispersion or transversetotransverse scattering in regions with coupling, such that the verticalmode action of the particle changes. Vertical dispersion and coupling are measured to be globally wellcorrected, 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 turnbyturn fitting procedure is no longer reliable. Instead, the turnbyturn images must be averaged first in order to improve signaltonoise, then fit as a single image. This has the disadvantages of incorporating a small amount of turnbyturn 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.
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 currentdependent effects diluting the vertical emittance at lowcurrent.
Ix Conclusions
A lowemittance tuning procedure has been developed at CesrTA, based on betatron phase and coupling measurements using resonant excitation and turnbyturn capable BPM. The tuning procedure has a fast turnaround, where one round of optics correction takes about ten minutes, and has yielded a singlebunch 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.0852.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 minimumachieved 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.
Acknowledgements.
The authors wish to thank the CESR operations and instrumentation groups, whose support was indispensable in our efforts to achieve wellcorrected conditions. This work was supported by the National Science Foundation grant PHY1002467 and the Department of Energy grant DESC0006505.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)):
(12) 
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 
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 
Element Class  Multipole  Value 

Sextupole with  a3  
Vert. Steering Trim  a5  
a7  
a9  
a11  
Sextupole with  a4  
Skew Quad Trim  a6  
a8  
a10  
a12  
Wiggler  a1 
References
 Palmer et al. (2009) M. A. Palmer et al., in Proceedings of the 2009 Particle Accelerator Conference, Vancouver, BC (2009) pp. 4200–4204.
 Dugan et al. (2009) G. F. Dugan, M. A. Palmer, and D. L. Rubin, in ICFA Beam Dynamics Newsletter, No. 50, edited by J. Urakawa (International Committee on Future Accelerators, 2009) pp. 11–33.
 Rubin et al. (2009) D. Rubin, D. Sagan, and J. Shanks, in Proceedings of the 2009 Particle Accelerator Conference, Vancouver, BC (2009) pp. 2751–2753.
 ilc (2013) The International Linear Collider Technical Design Report, Tech. Rep. (ILC Global Design Effort, 2013).
 (5) J. Safranek, Nucl. Instrum. Methods Phys. Res. A338.
 Safranek et al. (2002) J. Safranek, G. Portmann, A. Terebilo, and C. Steier, MATLABBased LOCO, Tech. Rep. (SSRL/SLAC and LBNL, 2002).
 (7) M. Aiba, M. Boge, N. Milas, and A. Streun, Nucl. Instrum. Methods Phys. Res. A694.
 (8) R. Dowd, M. Boland, G. LeBlanc, and Y.R. E. Tan, Phys. Rev. ST Accel. Beams 14.
 Wolski (2006) A. Wolski, Phys. Rev. ST Accel. Beams 9, 024001 (2006).
 Qin et al. (2013) H. Qin, R. C. Davidson, M. Chung, and J. W. Burby, Phys. Rev. Lett. 111, 104801 (2013).
 Baumgarten (2012) C. Baumgarten, Phys. Rev. ST Accel. Beams 15, 124001 (2012).
 Edwards and Teng (1973) D. Edwards and L. Teng, IEEE Trans. Nucl. Sci 20 (1973).
 Sagan and Rubin (1999) D. Sagan and D. Rubin, Phys. Rev. ST Accel. Beams 2, 074001 (1999), 10.1103/PhysRevSTAB.2.074001.
 Palmer et al. (2010) M. A. Palmer et al., in Proceedings of the 2010 International Particle Accelerator Conference, Kyoto, Japan (2010) pp. 1191–1193.
 Sagan et al. (2000) D. Sagan, R. Meller, R. Littauer, and D. Rubin, Phys. Rev. ST Accel. Beams 3, 092801 (2000), 10.1103/PhysRevSTAB.3.092801.
 Meller and Palmer (2011) R. E. Meller and M. A. Palmer, in Proceedings of the 2011 Particle Accelerator Conference, New York, NY (2011) pp. 504–506.
 Rubin et al. (2010) D. L. Rubin et al., Phys. Rev. ST Accel. Beams 13, 092802 (2010), 10.1103/PhysRevSTAB.13.092802.
 Chung et al. (1995) Y. Chung, D. Barr, G. Decker, J. Evans, K., and E. Kahana, in Particle Accelerator Conference, 1995., Proceedings of the 1995, Vol. 4 (1995) pp. 2452–2454 vol.4.
 Tecker et al. (1997) F. Tecker, B. Dehning, P. Galbraith, K. Henrichsen, M. Placidi, and R. Schmidt, in Particle Accelerator Conference, 1997. Proceedings of the 1997, Vol. 3 (1997) pp. 3648–3650 vol.3.
 Shanks et al. (2010) J. P. Shanks, D. Rubin, and D. Sagan, in Proceedings of the 2010 International Particle Accelerator Conference, Kyoto, Japan (2010) pp. 4620–4622.
 Rider et al. (2011) N. T. Rider, J. P. Alexander, J. A. Dobbins, M. G. Billing, R. E. Meller, M. A. Palmer, D. P. Peterson, C. R. Strohman, and J. W. Flanagan, in Proceedings of the 2011 Particle Accelerator Conference, New York, NY (2011) pp. 687–689.
 Rider et al. (2012) N. T. Rider, M. G. Billing, M. P. Ehrlichman, D. P. Peterson, D. Rubin, J. P. Shanks, K. G. Sonnad, M. A. Palmer, and J. W. Flanagan, in Proceedings of IBIC 2012: International Beam Instrumentation Conference, Tsukuba, Japan, Paper WECD01 (2012) pp. 585–589.
 Alexander et al. (2014) J. P. Alexander, A. Chatterjee, C. Conolly, E. Edwards, M. P. Ehrlichman, E. Fontes, B. K. Heltsley, W. Hopkins, A. Lyndaker, D. P. Peterson, N. T. Rider, D. L. Rubin, J. Savino, R. Seeley, J. Shanks, and J. W. Flanagan, Nucl. Instrum. Methods Phys. Res. A748C (2014), 10.1016/j.nima.2014.02.040.
 Shanks et al. (2011) J. Shanks et al., in Proceedings of the 2011 Particle Accelerator Conference, New York, NY (2011) pp. 1540–1542.
 Shanks et al. (2009) J. P. Shanks et al., in Proceedings of the 2009 Particle Accelerator Conference, Vancouver, BC (2009) pp. 2754–2756.
 Wolski et al. (2011) A. Wolski, D. Rubin, D. Sagan, and J. Shanks, Phys. Rev. ST Accel. Beams 14, 072804 (2011), 10.1103/PhysRevSTAB.14.072804.
 Sagan (2006) D. Sagan, Nucl. Instrum. Methods Phys. Res. A558, 356 (2006).
 Helms and Hoffstaetter (2005) R. W. Helms and G. H. Hoffstaetter, Phys. Rev. ST Accel. Beams 8, 062802 (2005), 10.1103/PhysRevSTAB.8.062802.
 Billing et al. (2013) M. Billing et al., in Proceedings of the 2013 International Particle Accelerator Conference (2013).
 Ehrlichman (2013) M. Ehrlichman, Normal Mode Analysis of Single Bunch, Charge Density Dependent Behavior in Electron/Positron Beams, Ph.D. thesis, Cornell University, Ithaca, New York (2013).
 Ehrlichman et al. (2012) M. P. Ehrlichman et al., in Proceedings of the 2012 International Particle Accelerator Conference, New Orleans, LA (2012) pp. 2970–2972.
 Crittenden et al. (2005) J. A. Crittenden, A. Mikhailichenko, E. Smith, K. Smolenski, and A. Temnykh, in Proceedings of the 2005 Particle Accelerator Conference, Knoxville, TN (2005) pp. 2336–2338.
 Sagan (2013) D. Sagan, The Bmad Reference Manual, Tech. Rep. (CLASSE, 2013).