âčïž Skipped - page is already crawled
| Filter | Status | Condition | Details |
|---|---|---|---|
| HTTP status | PASS | download_http_code = 200 | HTTP 200 |
| Age cutoff | PASS | download_stamp > now() - 6 MONTH | 0.8 months ago |
| History drop | PASS | isNull(history_drop_reason) | No drop reason |
| Spam/ban | PASS | fh_dont_index != 1 AND ml_spam_score = 0 | ml_spam_score=0 |
| Canonical | PASS | meta_canonical IS NULL OR = '' OR = src_unparsed | Not set |
| Property | Value |
|---|---|
| URL | https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8 |
| Last Crawled | 2026-03-18 01:49:02 (24 days ago) |
| First Indexed | not set |
| HTTP Status Code | 200 |
| Meta Title | First measurements of beam-beam effects in beam-separation, luminosity-calibration scans at the LHC | The European Physical Journal C | Springer Nature Link |
| Meta Description | At the CERN Large Hadron Collider (LHC), absolute luminosity calibrations obtained by the van der Meer (vdM) method are affected by the mutual electromagne |
| Meta Canonical | null |
| Boilerpipe Text | 1
Introduction
A precise determination of the absolute luminosity scale is essential for a wide range of LHC measurements, as it directly affects the normalization of many key physics cross sections. At the LHC, this calibration relies primarily on the van der Meer (vdM) method [
1
,
2
], which uses dedicated beam-separation scans performed under specially tailored beam conditions to relate measured interaction rates to the absolute luminosity inferred from beam parameters. While sub-percent uncertainties on the integrated luminosity have been achieved for Run-2 data by ATLASÂ [
3
] and CMSÂ [
4
], further improvements and robust experimental validation of beam-dynamical effects remain essential.
For relativistic proton beams colliding with zero crossing angle, which is the typical configuration for the
vdM
calibration, the luminosity at a given interaction point (IP) of the LHC is proportional to the overlap integral of the particle-density distributions
\(\rho _{\text {1,i}}\)
,
\(\rho _{\text {2,i}}\)
in bunch-pair
i
 [
5
]:
$$\begin{aligned} {\mathcal {L}}_{inst}&= 2c f_{\text {rev}} \sum ^{N_b}_i n_{\text {1,i}}, n_{\text {2,i}} \nonumber \\&\quad \times \iiiint ^{+\infty }_{-\infty }\rho _{\text {1,i}}(x,y,z-ct) \nonumber \\&\quad \times \rho _{\text {2,i}}(x,y,z+ct)\,dxdydzdt, \end{aligned}$$
(1)
where
c
is the speed of light,
\(f_{\text {rev}}\)
is the LHC revolution frequency,
\(N_b\)
is the number of colliding bunches, and
\(n_{\text {1,i}}, n_{\text {2,i}}\)
are the corresponding total charges per bunch for each colliding pair. Although only a very small fraction of the beam particles actually collide, the two opposing beams interact with each other electromagnetically, a dynamical process known as the beam-beam interaction [
5
,
6
]. In the LHC, as in any other synchrotron, the beams are not continuous but are divided into discrete âbunchesâ, each with a high proton density, containing approximately
\(10^{11}\)
protons over a bunch length of 7.5â9Â cm (in terms of RMS). This electromagnetic interaction occurs when the two beams share a common beam pipe, and is governed by a non-linear force that depends on the radial distance
r
of a test particle to the center of the opposing âsourceâ bunch with
n
protons Gaussian-distributed in the transverse planes (valid approximation for the LHC beams):
$$\begin{aligned} F = - \frac{ne^2}{2\pi \epsilon _0r}\biggl (1-\exp \biggl [ -\frac{r^2}{2\sigma ^2}\biggr ] \biggr ), \end{aligned}$$
(2)
where
e
is the elementary proton charge,
\(\epsilon _0\)
is the vacuum permittivity, and
\(\sigma \)
is the RMS radius of the transverse charge distribution within the bunch, in a simplified case of a round shape
\((\sigma = \sigma _x = \sigma _y)\)
. At very small amplitudes
\((r\rightarrow 0)\)
the force can be linearized giving the expression for the so called beam-beam parameter, which is often used to assess the strength of the force:
$$\begin{aligned} \xi = \frac{nr_p\beta ^*}{4\pi \gamma \sigma ^2} = \frac{nr_p}{4\pi \epsilon _n}, \end{aligned}$$
(3)
where
\(r_p\)
is the classical proton radius,
\(\beta ^*\)
is the optical beta-function at the interaction point (IP) â one of the CourantâSnyder parameters [
7
],
\(\gamma \)
is the Lorentz factor, and
\(\epsilon _n\)
is the normalized emittance.
Beam-beam effects were first observed almost 40 years ago [
8
] at the lepton colliders, and have been used extensively for accelerator diagnostics and optimization at both SLC and LEPÂ [
9
,
10
,
11
]. At the LHC, they were first observed in late Run 1Â [
12
], but only in terms of the beam-beam induced coherent deflections. Single particle effects, in contrast, were presumed too small in hadron colliders that until recently, there was no motivation to measure them with a precision that would be meaningful in the context evoked here.
Nonetheless, the full complexity of this interaction must be accounted for when performing the luminosity calibration via the vdM method, using the separation scans. These are designed to measure the detector-specific constant
\(\sigma _{vis}\)
that relates the observed rate
\(\mu ^{vis}\)
to the absolute instantaneous luminosity
\({\mathcal {L}}_{inst}\)
:
$$\begin{aligned} \mu ^{vis} = \frac{{\mathcal {L}}_{inst}\sigma _{vis}}{f_{rev}}. \end{aligned}$$
(4)
Under the assumption of uncorrelated particle densities in
x
and
y
planes, the transverse convolved bunch widths
\(\Sigma _x, \Sigma _y\)
can be extracted from the measured beam-separation dependence of the collision rate in the corresponding direction. The combined information from these scans, and bunch intensities can be used to calculate the instantaneous luminosity at the head-on position. The proportionality between the measured rate at the centered head-on position
\(\mu _{pk}\)
and the reconstructed luminosity from measured bunch parameters during the vdM calibration defines the luminometer-specific visible cross-section:
$$\begin{aligned} \sigma _{vis} = \frac{2\pi \Sigma _x\Sigma _y}{n_1 n_2} \mu _{pk}. \end{aligned}$$
(5)
If left uncorrected, the
\(\sigma _{vis}\)
measurement by the
vdM
method is biased by the beam-separation dependence of the mutual electromagnetic interaction of the two beams: the colliding bunches experience deflection-induced orbit shifts, as well as optical distortions akin to the dynamic-
\(\beta \)
effect, that both depend on the transverse beam separation and must therefore be accounted for when deriving the absolute luminosity scale. The correction strategy discussed extensively in Ref. [
13
] relies on simulation studies only, and thus requires experimental validation to establish its reliability.
The present report summarizes a campaign of dedicated measurements carried out in the spring of 2022 at the LHC, aimed at confirming the accuracy of the simulation studies detailed in Ref. [
13
], and at validating the beam-beam correction strategy used in the luminosity-calibration analyses of the ALICE, ATLAS, CMS and LHCb experiments. More specifically, the goal is to quantify the impact of the beam-beam interaction on the tune spectra, orbit, the transverse-density distribution of the colliding bunches, as well as, for the first time, on the luminosity, by systematically varying the strength or number of the beam-beam interactions. The measurements are repeated for different values of the beam-beam parameter to verify the scaling proposed in [
13
], for various choices of the scanning IP and for several multi-IP configurations. The beam conditions are chosen to be representative of
vdM
calibration sessions at the LHC, but optimized so as to maximize both the sensitivity of the measurements, and the operational efficiency.
This paper is organized as follows. The measurement strategy is detailed in Sect.Â
2
, followed by the main ingredients of this accelerator experiment: beam conditions (Sect.Â
3
), luminometers and accelerator instrumentation (Sect.Â
4
), and optimization of the ring lattice to maximize the sensitivity of the measurements to beamâbeam-induced optical distortions of the colliding bunches (Sect.Â
5
). SectionÂ
6
is devoted to the experimental characterization of the impact of beam-beam effects on the tune spectra, on the luminosity at different IPs, as well as on other observables such as transverse single-bunch widths and IP orbits. At each step, the results are confronted with the predictions of the COherent Multibunch Beam-beam Interactions (COMBI) tracking code [
14
], that produced most of the results detailed in Ref. [
13
]. A collection of macroparticles is simulated turn by turn, based on a simplified model of the accelerator lattice using linear transfer and the self-consistent computation of the beam-beam forces at the IP. The conclusions are presented in Sect.Â
7
.
2
Measurement strategy
The beam-beam interaction at an interaction point, referred to here as the
scanning
IP, induces measurable changes in key beam parameters such as the tune spectrum, closed orbit, and transverse single-beam profiles. These effects can be observed, albeit with limited precision, using standard accelerator diagnostics by comparing measurements obtained under varying transverse beam separation at the scanning IP. However, the resulting modifications to the transverse density distribution of the colliding bunches remain below the sensitivity threshold of conventional instrumentation, including synchrotron-radiation-based beam-profile monitors. Beam-beam-induced variations typically lie in the 0.5â1% range, small enough that only a highly precise measurement of collision-rate changes can provide adequate sensitivity.
Fig. 1
Scheme of the LHC and its interaction points (IP) with indicated IP1 as witness IP while performing a scan at IP5
Full size image
Importantly, the beam-beam interaction induces changes in beam properties that are intrinsically entangled with the
pp
collisions when observing luminosity at the scanning IP. Consequently, under realistic collider conditions, it is not possible to directly access a reference luminosity signal that is free of beam-beam effects â unlike in simulation studies, where such effects can be explicitly disabled. In a multi-collision configuration, however, luminosity shifts at one or more other IPs (Fig.Â
1
), where the beams are kept in continuous head-on collision and therefore act as non-scanning, or
witness
, IPs, can be used to monitor the separation-dependent beamâbeam effects. These effects are induced at the scanning IP, where the beams are deliberately brought in and out of collision, and propagate around the rings to the witness IPs, where they manifest as changes in the measured luminosity. In principle, any of the four LHC IPs can serve as a scanning IP, and any non-scanning IP can be designated as the witness IP for a given measurement. In practice, instrumental and operational constraints at the time the experiment was carried out restrict the choice of witness IP to the ATLAS (IP1) and CMS (IP5) collision points; IPs 1, 2 or 5, or combinations thereof, are used as scanning IP(s). Collisions at IP8 were deliberately avoided because the large longitudinal offset of the LHCb collision point, that breaks the eight-fold symmetry of the LHC rings, precludes the possibility of both members of a colliding-bunch pair to collide at all four IPs. In addition, the large crossing angle at IP8 would complicate the interpretation of the beam-beam effects that would occur at that IP if some of the bunches collided there.
The propagation of the beamâbeam-induced, amplitude-dependent
\(\beta \)
-beating from the scanning to a witness IP is controlled by the betatron phase advance between these two IPs. Since the targeted signatures are delicate to measure at best, it is natural to try and enhance them by adjusting this phase advance so as to maximize the sensitivity of luminosity shifts at the witness IP to beam-beam effects at the scanning IP; in doing so, however, the overall tunes must be preserved. Given the central and symmetric roles played by IP1 and IP5 in this experiment, the adjustment of their relative phase advance drove the optimization procedure that is detailed in Sect.Â
5
, and that achieved a threefold improvement in measurement sensitivity.
Central to the quality of the measurements is the controlled variation of the beam-beam parameter at the scanning IP. This can be achieved either:
by using
step scans
,
i.e.
by fully separating the beams at the scanning IP in either the horizontal or the vertical plane, and then bringing them back into head-on collision; or
by using
separation scans
,
i.e.
by scanning the beams transversely with respect to each other in either the horizontal or the vertical plane; or
by taking advantage of the natural beam-intensity decay and emittance growth to progressively reduce the beam-beam parameter.
Equally important is to ensure that the beams remain in head-on collision at the witness IP(s), such that reproducible luminosity shifts measured at these locations can be unambiguously correlated with controlled changes in the strength of the beam-beam interaction at the scanning IP. It was verified that the beam-beam-induced orbit shift and the potential non-closure of the orbit bumps used to control the beam separation at the scanning IP, did not significantly affect the actual separation, and thereby the measured collision rate, at the witness IP. In a few cases however, orbit drifts of uncontrolled origin ended up degrading the quality of some of the measurements.
Simulations demonstrate [
13
] that multi-IP effects strongly influence the magnitude of beam-beam biases to vdM calibrations. Their characterization, therefore, constitutes an essential component of this experiment. The scan protocol is a generalization of that in the two-IP case. One IP (say, IP1) is designated as the witness IP, with beams colliding head-on there and at IP2, and IP5; the beams are then taken out of collision at IP2 only, then at both IP2 and IP5, and finally returned to a three-collision configuration in the reverse order. Since the phase-advance combinations are different in each configuration, the corresponding luminosity shifts predicted at the witness IP are also different and can be meaningfully confronted with the data.
3
Beam conditions
In the round-beam approximation, the beam-beam parameter
\(\xi \)
does not depend on the beam energy (Eq. (
3
)). Therefore, and in order to minimize the overhead associated with commissioning and operating the LHC in a non-standard optical configuration, as well as to allow for frequent refilling with different bunch patterns, the measurements were carried out at injection energy (
\(450\ \text {GeV}\)
per beam), and under non-standard machine-protection conditions. The improvement in operational efficiency came at a cost:
the counting rate of the luminometers (or equivalently their visible cross-section) is about an order of magnitude smaller at
\(\sqrt{s} = 900\ \text {GeV}\)
than at
\(13\ \text {TeV}\)
, reflecting the combination of a 40% drop in the inelastic
pp
cross-section, a factor of two to three drop in the multiplicity of the final-state particles, and a significant softening of their momentum spectrum;
because of the total-intensity constraints dictated by machine-protection requirements, the injected beam could not exceed 4 bunches per beam (compared to 150 during a routine
vdM
-calibration session), with a maximum allowed population of
\(1.25\times 10^{11}~p\)
/bunch;
the combination of the low beam energy, that suppresses synchrotron-radiation damping, and of the large bunch intensity, that enhances intra-beam scattering, resulted in relatively rapid emittance growth and rather short single-beam lifetimes. To partially mitigate these effects, all collision-rate measurements are expressed in terms of specific luminosity, thereby automatically accounting for the beam-intensity decay.
The bunch intensity was deliberately increased by about 25% beyond the
vdM
-scan values typical of LHC Run 2, to maximize both the luminosity and the beam-beam parameter
\(\xi \)
. The latter reached 0.010 per IP (TableÂ
1
), compared to 0.003â0.006 during normal
pp
vdM
sessions. With a target injected emittance of
\(1.5 \ \upmu \text {m} \cdot \text {rad}\)
, the
\(\beta \)
function at the IP set to
\(\beta ^* = 11\ \text {m}\)
, and zero crossing angle at IP1 and IP5, the statistical uncertainty affecting a typically 60-second long per-bunch luminosity measurement in the presence of head-on collisions, lay around 0.5%.
Table 1 Range of beam-beam parameter values during the various stages of the experiment
Full size table
The bunch patterns were chosen such that all bunches in each beam collided either at IP1 and IP5 only, or at all three of IP1, 2 and 5. The orbit-stabilization feedback system was turned off during the data-taking periods to prevent it from interfering with the beamâbeam-induced orbit shift. The chromaticity was set to its standard value in physics fills, of +10 units, to guarantee coherent stability against the machine impedance, and in view of the very small number of bunches and of their large longitudinal spacing, the beam-stabilizing Landau octupoles were set to their minimum current (1A). The damping time of the bunch-by-bunch transverse feedback was set to 1000 turns, a rather loose setting intended to preserve longer the natural beam oscillation and therefore improve the precision of the tune measurements.
4
Beam instrumentation
4.1
Luminometers
The luminometer systems in use by the ATLAS and CMS collaborations at IP1 and IP5 are described in Refs. [
15
,
16
] respectively; no luminosity measurements were available at IP2 during the collider experiment described in this paper. Maximizing the statistical sensitivity leads to choosing the luminosity algorithm with the highest possible acceptance, namely:
the hit rate per bunch crossing in the ATLAS Minimum Bias Trigger Scintillators (MBTS)Â [
15
], and
the occupancy in the CMS Hadron Forward Calorimeter (HFOC)Â [
16
].
In both cases, raw luminometer counts are converted to collision rates using the Poisson formalism [
15
,
16
], and the instantaneous luminosity is averaged over typically 60-second time bins, that during the scans are synchronized with the scan steps. The statistical uncertainties are estimated from either the total number of raw luminometer counts per time bin (ATLAS), or from the RMS of the approximately 40 luminosity samples recorded in a given time bin (CMS). Since the absolute luminosity scale is irrelevant, and in order to simplify the interpretation of the results, all measurements in this paper are presented in terms of fractional shifts in the bunch-averaged specific luminosity, relative to a reference time that depends on the type of measurement considered.
4.2
Bunch-charge measurement
The total intensity of each beam is measured by a direct-current current transformer (DCCT), the absolute scale of which is calibrated against a very high precision pulse generator [
17
]. The bunch-by-bunch charge fractions, in turn, normalized to the total stored intensity, are measured, separately for the two beams, by Fast Beam Current Transformers (FBCT)Â [
18
]. The latter devices also provide the fill pattern,
i.e.
the relative location, at a given instant and around the two rings, of the nominally filled bunches. Since the absolute scale of the DCCTs is known to much better accuracy than that of the FBCTs, the bunch charges are typically computed as the product of the FBCT bunch-charge fractions and the total circulating beam intensity reported by the corresponding DCCT.
4.3
Emittance measurement
The synchrotron-light beam-profile monitors [
19
], dubbed BSRTs (âBeam Synchrotron Radiation Telescopeâ), are mainly used to track emittance evolution over time, but they can also measure the relative changes in transverse single-beam size that result from beamâbeam-induced
\(\beta \)
-beating. BSRT data are recorded every second; in what follows, they are presented averaged over one-minute intervals for easier comparison with other measurements.
The absolute length scale of the BSRT profiles suffers from significant uncertainties, if only because the online optical corrections to the images can be updated only a few times per year, and are unable to track the evolution of the efficiency of the light sensors, that depends both on time and on the position of the light spot on the sensor array. For each beam therefore, the BSRT is complemented by a wire-scanner (WS) profile monitor that is located in the same straight section. The advantage of the WS is that its accuracy is significantly better than that of the BSRT; its down side, in the present context, is that the WS cannot acquire data continuously and must be triggered manually. For the results presented in this paper, the wire scanners were flown through the beams at the start of each group of measurements to provide single-beam emittance measurements that are as accurate as possible, and then at regular intervals thereafter. The BSRT is used to interpolate the time evolution of the emittance between two sets of WS measurements.
To interpret measured RMS bunch widths in terms of emittance requires the knowledge of the
\(\beta \)
functions at the locations of the BSRT and the WS. These are determined by the phase-advance method, and compared in TableÂ
2
to their model value computed using MAD-X; the agreement is typically better than 5%, and the worst disagreement amounts to 9%. These discrepancies are attributed to imperfections in the magnetic lattice. Uncertainties in the measurement mainly arise from the interpolation of the lattice functions between the two closest beam-position monitors (BPMs) where the
\(\beta \)
-functions are measured [
20
,
21
], and either the BSRT or the WS location. Based on the MAD-X LHC lattice model, these uncertainties are estimated not to exceed 3%, and this value is assigned as the systematic uncertainty on the
\(\beta \)
functions used in BSRT- and WS-based emittance measurements.
Table 2
\(\beta \)
functions at the BSRT and WS locations, for the optical configuration detailed in Sect.Â
5.3
Full size table
The beam-averaged emittance,
i.e.
the average of the emittances of the beam-1 bunch and of the corresponding beam-2 bunch, can be obtained directly from emittance scans [
22
] at IP1 and/or IP5. In this approach, the convolved transverse bunch sizes measured using beam-separation scans are translated into emittances using the
\(\beta ^*\)
values determined by
k
-modulation [
23
] at the relevant IP (see Sect.Â
4.3
). The associated uncertainty (TableÂ
3
) is combined with the statistical uncertainty in convolved transverse width to estimate the error affecting the measured emittance. Comparison with beam-averaged emittances extracted from the single-beam profile monitors reveals excellent agreement between WS and emittance-scan results (Fig.Â
2
). The relative time-evolution of the BSRT emittances is consistent with that observed using the WS or emittance scans, but the absolute magnitudes differ significantly, especially in the horizontal plane.
Table 3
\(\beta ^*\)
functions at the ATLAS and CMS IPs measured by
k
-modulation. The target value in all cases is
\(\beta ^*=\)
11Â m
Full size table
Fig. 2
Time evolution of the horizontal (top) and vertical (bottom) beam-averaged normalized emittances during LHC fill 8037, as reported by the BSRT (blue), the WS (orange), and using beam-separation scans at IP1 (red) and IP5 (purple). Different shades of the same color correspond to two different bunches present in this fill. For ATLAS emittance scans, only the bunch-averaged value is shown
Full size image
4.4
Beam-beam parameter determination
The systematic and quantitative comparison of the measured and of the predicted beam-beam impact on the observables detailed in Sect.Â
2
requires continuous monitoring of the actual beam-beam parameter throughout the duration of the experiment. Only the BSRT provides uninterrupted emittance determination throughout the fill. In view of the scale biases apparent in Fig.Â
2
, however, WS-based emittances are used, whenever possible, as input to the determination of the beam-beam parameter; BSRT emittances re-scaled to close-in-time, absolute WS measurements provide time-interpolated emittance values whenever WS data are unavailable.
The beam-beam parameter evolution during the first fill of the experiment is illustrated in Fig.Â
3
. The error bands are dominated by the 3% systematic uncertainty in the measured optical functions. The strength of the beam-beam interaction drops by a factor of two over a couple of hours, from the combined effect of beam-intensity decay and of emittance growth.
Fig. 3
Time evolution of the beam- and plane-averaged beam-beam parameter
\(\xi \)
inferred from the measured bunch charges and emittances, separately for the two colliding-bunch pairs present in the fill pattern. The color bands indicate the systematic uncertainty (see text)
Full size image
4.5
Tune measurements
Coherent spectra can be monitored using either the LHC Transverse Damper (ADT)Â [
24
], or the Base-Band Tune (BBQ)Â [
25
] system; their comparison typically yields consistent results. Spectrograms with a tune resolution of 0.0001 are computed using bunch positions recorded at every turn, and averaged over the one-minute time bins mentioned in Sect.Â
4.1
. To mitigate the influence of noise on the tune measurement, a median filter is applied at the pre-processing stage with a self-defined local window size. Additionally, the 50Â Hz noise lines present in the spectra are masked. The filtered data are fitted by the sum of a Gaussian function and a constant baseline, and the frequency at the peak of the Gaussian is interpreted as the measured mean tune. The systematic uncertainty on measured tune shifts
\(\Delta Q\)
is estimated empirically to be
\(\sigma _{\Delta Q}=0.001\)
.
4.6
IP-orbit monitoring
Orbit displacements at each IP are measured using the Diode ORbit and OScillation (DOROS) BPM system, that provides sub-micrometer beam-position resolution [
26
]. These strip-line BPMs are located in the two quadrupoles on either side of and closest to the IP, allowing the position of both beams to be measured simultaneously. The position and the angle of each beam at the IP are inferred from, respectively, the average and the difference of the positions measured in the two final-triplet quadrupoles.
5
Optimization of the phase advance between interaction points
In the two-IP configuration, the basic theory of linearized beam-beam
\(\beta \)
-beating provides an analytical determination of the optimum phase advance between the two IPs,
i.e.
of the setting that maximizes the sensitivity, at the witness IP, to beam-beam effects at the scanning IP (Sect.Â
5.1
). A generalization to the three- and four-IP configurations is discussed in Sect.Â
5.2
. The required modifications to the baseline LHC optics, as well as their implementation and their experimental validation, are detailed in Sect.Â
5.3
.
5.1
Linear beam-beam
\(\beta \)
-beating
In circular colliders, the
\(\beta \)
-beating induced by beam-beam effects at the IP(s) leads to changes in transverse beam size [
5
,
27
], thereby impacting the luminosity. For a single collision point, the change in
\(\beta ^*\)
of a single small-amplitude particle in (for instance) the horizontal plane
x
, is given by:
$$\begin{aligned} \frac{\beta _x^*}{\beta _{0,x}^*} = \frac{1}{\sqrt{1-4\pi \xi \cot {(2\pi Q_x)} - 4\pi ^2\xi ^2}}. \end{aligned}$$
(6)
The magnitude of the effect depends on the absolute value
\(\xi \)
of the beam-beam parameter and on the nominal betatron tune
\(Q_x\)
; its periodicity isÂ
\( [\pi ]\)
. It can be shown that approaching the half-integer betatron tune from below minimizes the
\(\beta \)
function at the IP, thereby maximizing the luminosity; this phenomenon has been exploited with great success in (among others) the KEKB and PEP-II
B
factories. In the LHC, where during collisions the nominal fractional betatron tunes are fixed at
\(q_x=0.31,\,q_y=0.32\)
, this dynamic-
\(\beta \)
effect depends only on the beam-beam parameter.
This conclusion no longer strictly holds in the presence of more than one collision point. In the case of two IPs with identical values of
\(\xi \)
and
\(\beta ^*\)
, the dynamic-
\(\beta \)
effect can be described analytically [
28
]:
$$\begin{aligned} \frac{\beta _x^*}{\beta ^*_{0,x}} = \frac{\sin {2\pi Q_x} + 4\pi \xi (\cos (2\pi Q_x-2\mu _{1,x})-\cos 2\pi Q_x )}{\pm \sqrt{1-\cos ^2{2\pi (Q_x+\Delta Q_x)}}}, \end{aligned}$$
(7)
where
\(\mu _{1,x}\)
is the phase advance from the scanning IP to the witness IP. The sign in the denominator is linked to that of the numerator so as to ensure that the ratio remains positive; the denominator must be positive (resp. negative) below (resp. above) the half-integer tune,
i.e.
when
\(m< Q_x <m+\frac{1}{2}\)
(resp.
\(m+\frac{1}{2}< Q_x < m+1\)
), where
m
is a positive integer. The tune shift
\(\Delta Q_x\)
is related to
\(Q_x\)
,
\(\xi \)
and
\(\mu _{1,x}\)
by:
$$\begin{aligned} \cos {2\pi (Q_x+\Delta Q_x)}&=(1-16\pi ^2\xi ^2)\cos {2\pi Q_x} + 8\pi \xi \sin {2\pi Q_x} \nonumber \\&\quad + 16\pi ^2\xi ^2\cos {(2\pi Q_x-2\mu _{1,x})}. \end{aligned}$$
(8)
The phase-advance dependence of the resulting
\(\beta \)
-beating is illustrated in Fig.Â
4
. The minimum of each curve corresponds to the phase advance that yields the largest mutual beamâbeam-induced luminosity enhancement between the two IPs. The optimal settings are
\(\mu _{1,x}^{min}/2\pi =0.405,\,\mu _{1,y}^{min}/2\pi =0.410\)
; their difference reflects that between the nominal horizontal and vertical fractional tunes
\(q_x\)
and
\(q_y\)
.
Fig. 4
Analytically computed beam-beam induced
\(\beta \)
-beating in a two-IP configuration, as a function of the phase advance
\( \mu _1\)
between the 2 IPs (
\(q_x=0.31\)
,
\(q_y=0.32\)
, head-on collisions at both IPs,
\(\xi =7\times 10^{-3}\)
per IP)
Full size image
The extrema of the curves in Fig.Â
4
can be determined analytically by differentiating Eq. (
7
) with respect to the phase advance
\(\mu _{1,x}\)
. The general expression for this derivative
\(\frac{d}{d\mu _{1,x}} \biggl ( \frac{\beta _x^*}{\beta ^*_{0,x}} \biggr )\)
, detailed in Appendix A, is the product of two factors, each of which can be zero:
imposing that
\(\sin {(2\pi Q_x - 2\mu _{1,x})} = 0\)
yields the solution
$$\begin{aligned} \frac{\mu _{1,x}}{2\pi } =\frac{Q_x}{2}-\frac{m}{4} \, \end{aligned}$$
(9)
where
m
is an integer of either sign, or zero;
the other factor is zero for either
$$\begin{aligned} \frac{\mu _{1,x}}{2\pi } = \frac{Q_x}{2} - A(\xi , Q_x)\, \end{aligned}$$
(10)
or
$$\begin{aligned} \frac{\mu _{1,x}}{2\pi } = \frac{Q_x - 1}{2} + A(\xi , Q_x) \, \end{aligned}$$
(11)
with
\(A(\xi , Q_x)\)
defined in Eq. (
12
)
$$\begin{aligned} A(\xi , Q_x) = \arccos \left( \frac{ 32\pi ^3\xi ^3 \sin (4\pi Q_x) + 48\pi ^2\xi ^2 \cos ^2(2\pi Q_x) - 32\pi ^2\xi ^2 - 6\pi \xi \sin (4\pi Q_x) + \cos ^2(2\pi Q_x) - 1}{16\pi ^2\xi ^2 \left( 4\pi \xi \sin (2\pi Q_x) + \cos (2\pi Q_x) \right) } \right) / (4\pi ) \end{aligned}$$
(12)
The two sets of solutions do not overlap. The subset of solutions that minimize (rather than maximize)
\(\beta ^*\)
is identified by requiring that the second derivative be positive. This restricts the main solution to:
$$\begin{aligned} \frac{\mu _{1,x}}{2\pi } =\frac{Q_x}{2}+\frac{m+1}{4} . \end{aligned}$$
(13)
The optimal phase
\(\mu _{x,1}\)
that minimizes the
\(\beta ^*\)
for range of
\((\xi , Q_x)\)
values is shown in Fig.Â
5
.
Fig. 5
Optimal phase advance value
\(\mu _{x,1}\)
for various values of
\((\xi , Q_x)\)
in the two-IPs configuration
Full size image
The phase-advance dependence of the dynamic-
\(\beta \)
effect can also be quantified in terms of the luminosity shift, at the witness IP, that is associated with the electromagnetic interaction of the two bunches at the scanning IP (Fig.Â
6
). The optimal setting,
i.e.
that which maximizes the sensitivity, at the witness IP, to beam-beam effects at the scanning IP, is indicated by the green dotted vertical line at
\(\Delta \mu _{\text {IP1-IP5}}=0.41\ [2\pi ]\)
; the full suppression of the phase-related luminosity enhancement is indicated by the red dotted line. These results are consistent with the predictions of the analytical model (Fig.Â
4
). In the two-IP configuration considered here, in which IP1 and IP5 can both play the role of either the scanning or the witness IP, and due to the periodicity of this curve (
\( \pi \)
) and to the fractional-tune values used during collisions (
\(q_x / q_y = 0.31 / 0.32)\)
, the phaseâadvance-related luminosity enhancement is the same at the two IPs; at the optimum setting, it is three times larger than it would be using the nominal LHC lattice. The maximum effect depends mainly on the transverse tunes, and thus could have been further enhanced if it had been possible to move these away from their nominal values.
A more comprehensive simulation study, that evaluates the impact of the beam-beam interaction on the luminosity for more realistic Gaussian proton-density distributions, multiple collision points and as a function of the beam separation, is reported in Ref. [
29
].
Fig. 6
Beamâbeam-induced luminosity shift at IP5 predicted by COMBI simulations, as a function of the phase advance
\( \mu _1=\Delta \mu _{\text {IP1-IP5}}\)
between the two IPs (
\(q_x=0.31\)
,
\(q_y=0.32\)
, head-on collisions at both IPs,
\(\xi =7\times 10^{-3}\)
per IP)
Full size image
5.2
Optimized phases for multi-IP configurations
The general expression to obtain the first-order change of the
\(\beta \)
function at an arbitrarily chosen reference IP (at
\(s=\mu _{0,x}=0\)
) can be approximated as multiple quadrupole errors [
30
]:
$$\begin{aligned} \frac{\beta ^*_x}{\beta ^*_{0,x}} = \frac{2\pi \xi }{\sin {2\pi Q_x}} \sum _{i\in IPs}^{N-1} \cos (2\pi Q_x-2\mu _{i,x})\, \end{aligned}$$
(14)
where
\(\mu _{i,x}\)
is the horizontal phase advance between the reference and the
\(i^{th}\)
IP for a total of
N
IPs. Therefore, the optimal configuration, which minimizes the
\(\beta ^*\)
change at all IPs, must satisfy the condition:
$$\begin{aligned} \sum _{i=1}^{N-1} \cos (2\pi Q_x-2\mu _{i,x}) = -(N-1)\, \end{aligned}$$
(15)
which yields the same solution for all IPs as in the 2 IPs configuration (Eq.Â
13
):
$$\begin{aligned} \frac{\mu _{i,x}}{2\pi } =\frac{Q_x}{2}+\frac{n+1}{4} . \end{aligned}$$
(16)
where
n
is an integer of either sign, or zero. To maximize the effect on luminosity, the phase advance between the reference interaction point at
\( i = 0 \)
and the nearest IP, at
\( i = 1 \)
must satisfy
\( \mu _{1,x}/2\pi = (Q_x + 0.5)/2 \)
, as given by the first solution (for
\(n=0\)
) of Eq. (
13
). This condition implies that the phase advance to any subsequent adjacent IPs (
\( i \ge 2 \)
) must be an integer multiple ofÂ
\( \pi \)
. Thus, this requirement exposes a fundamental constraint: it is not possible to minimize
\( \beta ^* \)
simultaneously at all IPs. This would only be possible for (
\( Q_x = 0.5 \)
) which is practically not achievable, as it would require operating the storage ring precisely at the half-integer resonance.
5.3
Optics validation
The optimum phase advance from IP1 to IP5 determined in Sect.Â
5.1
was implemented using tuning trim quadrupoles in the LHC arcs. The targeted fractional phase advance was
\(0.9 \, [2 \pi ]\)
for both beams and both planes (TableÂ
4
). This setting, which differs from that in Fig.Â
6
by one period (
\( \pi \)
), is equivalent in terms of sensitivity optimization but has the advantage that it minimizes the absolute magnitude of the perturbations inflicted upon the nominal lattice.
Several measurements were carried out to validate this adjustment: comparison of the targeted and achieved shifts in phase advance (Sect.Â
5.3.1
),
\(\beta \)
-function measurements across the full LHC circumference (Sect.Â
5.3.2
) as well as at the IPs (Sect.Â
5.3.3
), and dispersion measurements to estimate the dispersive contribution to the transverse beam size at the IP (Sect.Â
5.3.4
).
5.3.1
Phase advance
The phase advance around each LHC ring is determined using turn-by-turn orbit measurements [
20
,
21
]. The targeted and achieved phase-advance shifts agree within
\(0.01\, [2\pi ]\)
or better (TableÂ
4
). The phase advance change from IP1 and IP5 was compensated on the opposite side of the ring, ensuring that the transverse tunes remain unchanged. The target fractional local phase advance is set to
\(0.9,[2\pi ]\)
. TableÂ
4
lists the corresponding phase shifts applied to the nominal lattice to reach this target on one side of the ring; this is consistent with the analytically derived value of
\(0.41\,[2\pi ]\)
for the compensating shift on the opposite side. The change in phase advance was also measured at the BSRT location, where the discrepancy between the model and measurements is more pronounced, reaching up to
\(0.05\,[2\pi ]\)
.
Table 4 Nominal and optimized IP1 to IP5 phase-advance settings. Column 4 (resp. 5) displays the targeted (resp. measured) shift in phase advance after the optimized settings have been applied
Full size table
5.3.2
Residual
\(\beta \)
-beating
The
\(\beta \)
-beating was quantified and found to be well within acceptable tolerances for various configurations of the trim strength. The measurement results, relative to the reference measurements without any phase advance adjustments, are illustrated in Fig.Â
7
. The maximal observed deviation is approximately 25% around IP5 for Beam 1 in the vertical plane. Such significant deviations in proximity to the IPs are anticipated and are frequently attributed to erroneous BPM readings. Overall, the measurements indicated that
\(\beta \)
-beating remained within acceptable bounds, predominantly within a 10% margin.
Fig. 7
Measured
\(\beta \)
-function change from the full strength of the phase advance knob with reference to nominal lattice along the LHC ring, for Beam 1 (top two figures for
x
and
y
) and Beam 2 (bottom two figures for
x
and
y
)
Full size image
The measured
\(\beta \)
functions are also compared to the MADX LHC lattice model predictions, which include the phase advance adjustment. The results, illustrated along the LHC lattice in Fig.Â
8
, indicate that the measurements align with the model within an acceptable margin of 10%.
Fig. 8
Measured deviations in the
\(\beta \)
-function along the LHC ring relative to the MADX model incorporating the phase advance knob, for Beam 1 (top two figures for
x
and
y
) and Beam 2 (bottom two figures for
x
and
y
). The error bars represent solely the statistical uncertainties as determined in the measurements
Full size image
5.3.3
\(\beta \)
functions at the interaction points
The DOROS BPMs are strategically located at the end of the inner triplets on both sides of each of the experiments. Due to the phase advance of
\(\pi \)
between the DOROS on the left and right sides, these are inadequate for the precise determination of the measurement of
\(\beta ^*\)
at the IP. This can only be achieved with a special
k
-modulation technique [
31
], which was performed at the interaction points used in the experiment. The results are summarized in TableÂ
3
. While the achieved precision of the
\(\beta ^*\)
determination varies between measurements, uncertainties at the level of
\({\mathcal {O}}(1\%)\)
are not uncommon and are not systematic to a specific beam or interaction point. Such uncertainties would propagate to only a few-percent effect on derived quantities such as the emittance or the beamâbeam parameter.
5.3.4
Dispersion
It was essential to perform dispersion
\(D^*\)
measurements at the interaction points, under the new optics configurations, as it can contribute to the measured beam size. Assuming the LHC design momentum spread value of
\(\frac{\Delta p}{p}=10^{-4}\)
, the average dispersion measurement around IP1 and IP5 is estimated to contribute up to 2% to the observed bunch sizes (worst case measured for B2Y). This effect is considered static throughout the experiment.
6
Experimental results
This section presents experimental measurements of beamâbeam effects in configurations involving multiple collision points. The impact of the optics modifications introduced in Sect.Â
5
is first assessed through tune spectra measurements, as detailed in Sect.Â
6.1
. To quantify the cumulative effect of one or more additional collisions and assess the reproducibility of the results, a stepwise collision scan was employed. In this procedure, beams were sequentially brought into head-on collision at selected interaction points (IPs), followed by full transverse separation (
\(\Delta > 6\sigma \)
) in either the horizontal or vertical plane at all but one witness IP. Two distinct filling schemes were used to isolate the influence of a single additional collision (2-IP configuration: ATLAS and CMS, see Sect.Â
6.2
) and to investigate more complex multi-collision scenarios (3-IP configuration: including ALICE, see Sect.Â
6.3
). The first measurement of separation-dependent beamâbeam effects, analogous to those encountered during luminosity calibration scans, is presented in Sect.Â
6.4
, alongside comparisons with numerical predictions. An additional measurement conducted during the 2022 van der Meer calibration is detailed in Sect.Â
6.5
. A consolidated overview of all findings is provided in Sect.
7
.
6.1
Measured effect on the tune spectra of the phase advance changes
The coherent transverse tune spectra serve as a diagnostic for evaluating the effectiveness of the phase advance modifications. COMBI simulations were carried out for two extreme lattice configurations, one designed to maximize and the other to suppress beamâbeam effects. In the head-on scenario (Fig.Â
9
), the coherent modes are suppressed for the maximizing configuration, as the phase advance between the two collision points is close to
\(\pi \)
, as was previously shown in simulation in Ref. [
32
]. In this configuration, the effect of the first collision is counteracted by the second. This setting is particularly suitable for studying incoherent spectra, as the diminished coherent signal allows for a clearer extraction of the mean tune shift from fitted spectra. During separation scans, the symmetry between horizontal (
x
) and vertical (
y
) planes is broken. In the scanning plane (horizontal), coherent modes re-emerge at
\(1.5\,\sigma \)
separation (Fig.Â
10
), with their frequency shifting in response to the separation distance (additional example for
\(3\,\sigma \)
is shown in Fig.Â
11
).
Fig. 9
Comparison of COMBI simulated tune spectra in
x
-plane (left) and
y
-plane (right) for the suppressing (in blue) and maximizing (in green) phase configurations, for head-on collision
Full size image
Fig. 10
Comparison of COMBI simulated tune spectra in
x
-plane (left) and
y
-plane (right) for the suppressing (in blue) and maximizing (in green) phase configurations, at
\(1.5\,\sigma \)
separation in the horizontal plane
Full size image
Fig. 11
Comparison of COMBI simulated tune spectra in
x
-plane (left) and
y
-plane (right) for the suppressing (in blue) and maximizing (in green) phase configurations, at
\(3\,\sigma \)
separation in the horizontal plane
Full size image
These simulated trends of the coherent spectra are also observed in experimental data. Transverse spectra recorded for head-on conditions confirm the suppression of coherent modes in the maximizing phase configuration (Fig.Â
12
). Additional comparisons for separated beams at the IP (see Figs.Â
13
andÂ
14
), shown here for the case of vertical separation, reveal similar features. In this configuration, the non-scanning (horizontal) plane exhibits spectra that qualitatively agree with simulation predictions (see Figs.Â
10
andÂ
11
), while in the scanning (vertical) plane the coherent modes do not reappear as expected. An equivalent behavior is observed when the separation is applied in the horizontal plane, with the roles of the two transverse planes interchanged; this case is not shown for brevity. This discrepancy is attributed to the effect of the transverse damper, which acts independently on the two beam centroids and suppresses coherent oscillations arising from beamâbeam coupling, consistent with earlier observations [
33
].
Fig. 12
Measured tune spectra before and after phase change, corresponding to the suppressing (in blue) and maximizing (in green) phase configurations, at the head-on position, for each of the transverse planes (Beam 1 shown as an example). Dash-dotted lines show Gaussian fits to the data
Full size image
Fig. 13
Comparison of the measured tune spectra at suppressing (blue) and maximizing (green) phase configurations with
\(1.3\,\sigma \)
separation in the vertical plane
Full size image
Fig. 14
Comparison of the measured tune spectra at suppressing (blue) and maximizing (green) phase configurations with
\(2.6\,\sigma \)
separation in the vertical plane
Full size image
6.2
Head-on collisions in a two-IP configuration at
\(\sqrt{s} = 900\ \text {GeV}\)
Each test consisted of four measurements designed to quantify the beamâbeam effect induced by a single additional head-on collision. The sequence was executed with CMS as the witness IP, followed by ATLAS, and concluded with a repeat at CMS. Over the course of the fill, the beamâbeam parameter decreased from approximately 0.01 to 0.005. The outcome of each scan is analyzed using three independent observables: luminosity, beam size, and tune spectra, discussed separately below.
6.2.1
Impact of beam-beam effects on head-on luminosity
An overview of the luminosity measurements during the first test with CMS as the witness IP is shown in Fig.Â
15
, where the specific luminosity at CMS is normalized to the initial head-on step to highlight relative changes. The uncertainty at each step is given by the standard error of the mean; shaded bands represent the standard deviation. A comparable dataset for ATLAS is presented in Fig.Â
16
, confirming the expected luminosity change due to beamâbeam effects originating from another IP. Relative luminosity changes for each step in all tests are summarized in TableÂ
5
. These values are computed relative to the average of adjacent head-on points, which removes the linear luminosity decay over time. Assuming slow change of beam parameters, the results within each test are expected to be consistent.
Fig. 15
CMS specific luminosity, normalized to the first measured both IPs head-on point (blue), during the first step-function scan at IP1. The measurements with ATLAS fully separated are shown in red. The error bars indicate error on the mean
Full size image
Fig. 16
ATLAS specific luminosity, normalized to the first measured both IPs head-on point (blue), during the step-function scan at IP5. The points for which IP5 was separated are shown in red. The error bars indicate the statistical error
Full size image
Table 5 Measured luminosity shifts in each of full separation steps â summary of all three step-function scans performed for IP1 and IP5 as the witness IP. The beam-beam parameter at the start of each scan
\(\xi _\text {start}\)
is included for reference
Full size table
All measurements are compiled in Fig.Â
17
and compared to COMBI simulation predictions. The results are plotted as a function of the beamâbeam parameter (estimated as described in Sect.Â
4.4
), following a reverse chronological order. Agreement with simulation is excellent across the entire range.
Fig. 17
COMBI simulated luminosity enhancement induced by a single head-on collision as a function of the beam-beam parameter (dashed line), compared to the test results at both ATLAS (red points) and CMS (green points)
Full size image
Fig. 18
Measurements of beam-beam effects on beam widths shown with red points for each beam and plane, during the experiment with IP5 as the witness IP at
\(\xi =0.0056\)
. COMBI predictions are shown in blue with calculated systematic error from the phase advance error between the IP and BSRT locations
Full size image
6.2.2
Impact of beam-beam effects on single-beam sizes
Complementary measurements of the bunch size from the BSRT are shown in Fig.Â
18
. Data from all steps and bunches are combined, and a linear trend corresponding to continuous emittance growth is subtracted using the difference between the first and last head-on steps. For comparison, COMBI simulation results are also included in the figure. The dominant source of uncertainty in COMBI predictions arises from discrepancies between the measured phase advance and the MADX model, with the largest observed deviation being
\(0.05\,[2\pi ]\)
(as discussed in Sect.Â
5.3.1
). This uncertainty is propagated using a modified form of Eq. (
6
) adapted to the BSRT location. The magnitude of this uncertainty depends on the specific phase advance configuration, as some are more sensitive to small deviations. This is reflected in the varying sizes of the error bands associated with each prediction. The most accurate predictions are obtained in the vertical plane of Beam 2, where the beamâbeamâinduced
\(\beta \)
-beating propagates to the BSRT with a particularly favorable phase advance. In this configuration, the measured effect is also the most pronounced, and excellent agreement with COMBI simulations is observed. In other cases, the expected absolute effect is smaller, typically below 1%, and the influence of phase advance uncertainty becomes more significant, limiting the predictive reliability of the simulations. Experimentally, it is also possible that in these configurations the beamâbeamâinduced effect is too subtle to resolve, potentially masked by instrumental limitations.
6.2.3
Impact of beam-beam effects on tune spectra
FigureÂ
19
shows representative ADT spectra recorded during the step-function scans. These spectra are affected by coupling, noise, and various nonlinearities, so they do not follow an idealized analytical distribution. Nonetheless, the peak position can be extracted using a Gaussian fit (indicated by dashed lines), yielding an estimate of the mean tune shift. In the single-collision configuration (blue lines), unsuppressed coherent modes dominate the spectrum and partially obscure the underlying incoherent component. The Gaussian fit is therefore applied to the incoherent contribution between coherent lines, rather than to the global spectral maximum, which results in a larger uncertainty on the extracted peak position. Despite this, in each of the performed step-scans at both collision points, the single-collision tune shift, defined as twice the fitted mean shift, is found to be consistent within
\(\pm 5\%\)
with the expected coherent beamâbeam tune shift corresponding to the conditions of the respective scan. For head-on collisions with Gaussian transverse beam distributions, the coherent tune shift of the dipole mode is expected to be
\(Y\xi \)
, where
Y
is the Yokoya factor (
\(Y = 1.21\)
for round beams)Â [
34
]. As expected, the two-collision configuration yields a tune shift approximately double that of the single-collision case [
13
].
Fig. 19
Measured vertical tune spectra for different collision configurations: reference nominal tune (red), without the presence of any collisions, compared to single collision (blue) and two collisions (green) at
\(\xi \simeq 0.0086\)
. The dash-dotted lines in corresponding colors indicate Gaussian fits to the data used to extract mean tunes
Full size image
Fig. 20
Specific luminosity, normalized to the first head-on point, measured at ATLAS (top) and CMS (bottom) during beam separation at the other IPs, with horizontal separation applied in steps 2â4 and vertical separation in steps 6â8. The tested configurations are indicated in the legend. The blue line shows an interpolation of the head-on measurements, serving as a reference for relative luminosity changes
Full size image
6.3
Head-on collisions in a three-IP configuration at
\(\sqrt{s} = 900\ \text {GeV}\)
An additional head-on collision was included in the configuration using IP2 (ALICE), as it was the only IP where this was feasible. A specially designed filling scheme ensured uniform collision exposure across all bunches, allowing measurements to be combined and thereby reducing statistical uncertainties. However, no reliable luminosity data were available from IP2 during this experiment.
Table 6 Luminosity shift per full separation step, calculated with respect to an interpolated reference obtained from the two closest head-on points. Results for horizontal and vertical separations are shown as Shift part 1 and Shift part 2, respectively
Full size table
6.3.1
Impact of beam-beam effects on head-on luminosity
A step-function scan was conducted by sequentially collapsing or separating one IP at a time, with the separation applied first in the horizontal plane and, in a second stage of the experiment, in the vertical plane. A comprehensive summary of the measured specific luminosity at ATLAS and CMS is presented in Fig.Â
20
. Exceptionally stable beam conditions during this test enabled measurements at both experiments under nearly identical beamâbeam parameters. The beamâbeam-induced luminosity changes observed during each step are quantified in TableÂ
6
, referenced to the time-interpolated specific luminosity from three head-on points.
Orbit-drift corrections are applied only to these head-on reference points and result in changes of the extracted luminosity shifts of at most 0.2%. Good agreement with COMBI simulations is observed across the separation steps, with measured shifts closely matching the simulated values for configurations involving IP1 and IP5. For single-IP separation, the level of agreement depends on the witness IP and the separation step. When IP1 is used as the witness, the luminosity shift in the first separation step is reproduced almost exactly by COMBI, while a smaller shift is observed in the second step. This behaviour is expected, as the beamâbeam parameter decreases during the scan and the measurement becomes increasingly sensitive to residual orbit perturbations. For configurations involving IP2, the comparison is less conclusive. As no reliable luminosity data were available from IP2 during this experiment, it was not possible to precisely optimize the collision conditions at this IP. As a result, the actual beamâbeam conditions at IP2 may differ from those assumed in the simulations, affecting the comparison. Residual orbit drifts affecting the separation steps cannot be fully corrected. As illustrated by the DOROS measurements in Fig.Â
21
, the orbit data contain both slow drifts and beamâbeamâinduced orbit offsets arising from changes in the collision configuration when other IPs are separated. While the beamâbeamâinduced offsets are included in the COMBI simulations, the residual orbit drifts are not. Since these contributions cannot be disentangled in the measurements, a full correction is not possible. Consequently, the effective beamâbeam impact measured during the experiment may be weaker than that assumed in the simulations.
Fig. 21
Residual separation between the two beams at all three IPs in horizontal (top) and vertical (bottom) planes, for test with observations at ATLAS, where the separation was performed first in the horizontal and then in the vertical plane at the other IPs
Full size image
6.3.2
Tune shift scaling with the number of beamâbeam interactions
Analysis of the transverse tune spectra in this configuration is further complicated by the reappearance of coherent modes when beams are separated at one or more IPs, as predicted by simulation. In such cases, the maximum tune shift is determined from the frequency difference between the
\(\pi \)
and
\(\sigma \)
modes, rather than from fitting the incoherent peak. For configurations where all three IPs were in collision, coherent modes remained suppressed due to the retained symmetry between IP1 and IP5. FigureÂ
22
presents the measured tune shifts for all tested collision configurations. The unperturbed tune measurement, shown at â0â collisions, serves as the reference from which all tune shifts are calculated. An empirical uncertainty of 0.001 is assumed for each independent measurement. Both transverse planes and separation directions are included to demonstrate consistency. As expected, the measured tune shift scales linearly with the number of head-on collisions.
Fig. 22
Measured mean tune shift as a function of the number of head-on collisions, shown for both transverse planes and for both separation directions (horizontal and vertical) used in the experiment. The separation refers to how collisions were removed: it was applied in only one transverse direction at a time (H sep. or V sep.) at the non-colliding interaction points, which are not explicitly labeled on the x-axis
Full size image
6.4
Beam-beam signatures during transverse-separation scans at
\(\sqrt{s} = 900\ \text {GeV}\)
In the van der Meer (vdM) calibration method, beamâbeam corrections are applied separately at each separation step. Therefore, it is essential to validate the simulated dependence of luminosity shifts on nominal beam separation. Full vdM-like scans were performed in both transverse planes at ATLAS and CMS, using a configuration where only these two IPs were in collision.
6.4.1
Impact of beamâbeam effects on the luminosity scan-curves
FigureÂ
23
shows the luminosity curves measured at ATLAS, acting as the witness IP, during a separation scan performed at CMS. The curves are fitted with Gaussian profiles for visual guidance. The beam separation at CMS is expressed in units of the measured transverse beam size,
\(\sigma _{\textrm{meas}}\)
, which was found to be smaller than the nominal beam size,
\(\sigma _{\textrm{nom}}\)
, used to define the scan steps. As a result, scan points corresponding to nominal separations beyond
\(6\,\sigma _{\textrm{nom}}\)
extend to larger values when expressed in units of
\(\sigma _{\textrm{meas}}\)
.
To account for emittance evolution during the scan, both linear and exponential fits are applied to the head-on luminosity points before, after, and in the middle of the scan. These trends are subtracted from the data in Fig.Â
23
, and all measurements are presented relative to this fitted head-on reference (normalized to unity).
Fig. 23
ATLAS specific luminosity relative changes (blue points) during separation scans in
x
-plane (top) and
y
-plane (bottom) at CMS, corrected for exponential (exp. corr.) luminosity decay. Gaussian function with a constant fitted to these points is shown with blue dashed line and compared to fit performed on the data corrected for linear luminosity change (lin. corr.) in time (dotted line). These are compared to COMBI simulation predictions (purple solid line)
Full size image
At CMS, similar scans were performed; however, degraded beam stability was observed during this period of the fill. The scan results exhibit systematic deviations, especially at large separation, which could not be symmetrically corrected using simple linear or exponential models. Despite these issues, all four scan cases (two transverse planes at two IPs) yield a total specific luminosity reduction of approximately 2.4%, consistent with expectations based on the measured beamâbeam parameter
\(\xi \)
. Comparison with COMBI simulations reveal up to 30% discrepancy for IP5 as the witness IP, while agreement improves to within 10% for measurements at IP1. These differences are comparable to those observed in TableÂ
6
, suggesting a common underlying origin. While a finer scan granularity can improve the description of the steep-response region below
\(3\,\sigma _{\textrm{meas}}\)
, it would not address the dominant discrepancy between head-on and fully separated luminosity, which is likely driven by simulation assumptions and scan-related systematics.
6.4.2
Beamâbeam-induced orbit deflections
Beamâbeam deflections are clearly observed during separation scans. FigureÂ
24
shows the fully corrected residual horizontal beam displacement at IP1, measured with the DOROS beam position monitors, as a function of the beam separation. This representation allows a direct comparison with the analytical beamâbeam deflection prediction based on the BassettiâErskine formalism [
35
].
The raw orbit data were corrected for slow orbit drifts by aligning the beam positions at the luminosity optimization points before and after each scan. Deviations of the applied beam separation from the nominal settings were accounted for by fitting a common length scale simultaneously with the beamâbeam deflection. This approach avoids an overestimation of the separation calibration. The measured orbit displacement at the DOROS location was corrected for the contribution arising from the beamâbeam deflection angle by applying a geometric amplification factor, proportional to the distance between the IP and the BPM in the small-angle approximation. After all corrections, the measured beamâbeam deflection as a function of beam separation closely follows the analytical prediction. The maximum horizontal deflection observed at IP1 is 10.2Â
\(\upmu \)
m, compared to an expected value of 8.4Â
\(\upmu \)
m, with agreement within the total measurement uncertainty. The dominant uncertainty contribution (3.5Â
\(\upmu \)
m, corresponding to about 90% of the total) arises from a possible residual half crossing angle at the IP of up to 0.5Â
\(\upmu \)
rad. Additional systematic uncertainties include the fitted separation length scale, tune uncertainty of 0.005, and an extended 2Â % uncertainty on
\(\beta ^*\)
. Statistical uncertainties are negligible for most separation steps.
Fig. 24
Orbit offset caused by the beamâbeam interaction measured by DOROS BPMs in the horizontal plane at IP1 during the horizontal separation scan. Analytical prediction is shown for comparison [
35
]
Full size image
6.4.3
Impact of the beam-beam effects on beam sizes and tune shift
Separation-dependent changes in beam size were measured using the synchrotron light monitor. FigureÂ
25
shows good consistency across individual bunch measurements, as well as strong agreement with COMBI simulations. This reinforces the utility of the BSRT as a reliable diagnostic for beamâbeam effects. The uncertainties shown are statistical only and are treated as uncorrelated between scan points. The normalization to the zero-separation value is applied for visualization purposes and does not impose a constraint on that point.
Fig. 25
BSRT beam-size measurement in the vertical plane for each bunch (points) during horizontal (top) and vertical (bottom) separation scans at IP1. The exponential evolution correction is applied based on adjacent head-on points (before, in the middle and after the scan) to highlight the relative changes. COMBI predictions are shown (gray dashed line) for comparison
Full size image
Transverse tune shifts were extracted from ADT-measured spectra. As illustrated in Fig.Â
26
, the mean tune shift evolution is reconstructed for both the scanning and non-scanning planes. The measured tune shifts closely follow the expected COMBI trends, further validating the simulation models used for beamâbeam corrections.
Fig. 26
Measured tune shift in units of the beam-beam parameter
\(\xi _{BB}\)
for the horizontal (blue) and vertical (orange) planes, during a horizontal separation scan at IP1. The error bars indicate typically assumed conservative empirical systematic error on the tune measurement of 0.001. The solid curves represent the COMBI prediction
Full size image
6.5
Beam-beam signatures during vdM scans at
\(\sqrt{s} = 13.6\ \text {TeV}\)
During the first vdM calibration of LHC Run 3, an unusually high beamâbeam parameter was recorded, offering a unique opportunity to study beamâbeam effects under enhanced interaction conditions. These effects are examined using online ATLAS luminosity data, with representative results shown in Fig.Â
27
.
The example focuses on a standard vdM scan pair conducted at CMS, while monitoring luminosity at the non-scanning IP (ATLAS). The plotted luminosity changes are normalized to the central scan point, where both IPs are in head-on collision. Variations in this normalized luminosity across the scan directly reflect the impact of beamâbeam interaction at the scanning IP on the luminosity at the non-scanning IP. These patterns are reproducible across multiple scan pairs acquired within the same fill. At the largest nominal separations (approximately
\(\pm 0.6\,\)
mm), the beams are sufficiently separated to eliminate beamâbeam effects entirely. Notably, a qualitative difference is observed between the horizontal and vertical scan planes, an asymmetry that is not present in the above presented dedicated experiment. This behavior is attributed to differences in the phase advance between IPs in the standard vdM optics configuration. For vertical (
y
) scans, experimental data agree well with simulation, confirming the robustness of the modeling under these conditions. In contrast, the horizontal (
x
) scans exhibit a systematic discrepancy: the observed beamâbeam impact is broader than that predicted by simulation. The model assumes that the separation steps are accurately defined using the beam size measured during the vdM scans, but this mismatch may suggest an unmodeled phase advance deviation or another optics imperfection. The total difference between data and simulation is within 0.2%. This dataset originates from a standard vdM calibration, for which only online luminosity measurement was available, and many contributions to the overall vdM uncertainty, such as bunch-to-bunch variations, averaging over different bunch families, and emittance evolution during the scan, are not disentangled. Therefore, the 0.2% difference should be interpreted as a consistency check under typical operational conditions rather than a precision validation of the beamâbeam model. Nevertheless, this indicates that the model remains sufficiently accurate for beamâbeam corrections in vdM analyses.
Fig. 27
Changes from the separation-dependent beam-beam effects at the non-scanning IP â ATLAS online luminosity shown with points. Data comes from a standard vdM scan pair during the 2022
pp
luminosity calibration. Separation steps and direction are indicated with gray markers. The last CMS scan is shown, with the highest beam-beam parameter of
\( \xi =5.3 \times 10^{-3} \)
. The corresponding COMBI simulation results are also shown (triangles)
Full size image
7
Conclusions
Beamâbeam effects were directly observed and quantitatively validated at the LHC. The impact of beamâbeam interactions on head-on luminosity, bunch sizes, tune spectra, and beam orbit positions was measured across various collision configurations. The observed effects align closely with COMBI simulations. In particular, the luminosity enhancement due to amplitude-dependent
\(\beta \)
-beating was measured for the first time and is found to agree with expectations over a wide range of the beamâbeam parameter
\(\xi \)
. Together with the observed linear scaling with respect to the beamâbeam parameter, this provides experimental confirmation of the numerical studies presented in [
13
]. Consistent signatures were found across multiple scan configurations, including two- and three-IP setups, and the characteristic tune shifts were well reproduced by both measurement and simulation.
Accurate instrumentation and experimental design were essential. Key to the results was the stepwise collision scheme, isolating individual beamâbeam contributions, with the implementation of the witness IP. Synchrotron light measurements from the BSRT, especially in the vertical plane of Beam 2, showed excellent agreement with COMBI predictions and bunch-by-bunch consistency. Tune spectra from the ADT system provided direct access to mean tune shifts, while orbit data from the DOROS system enabled measurement of beamâbeam deflection down to 1Â
\(\upmu \)
m, matching analytical estimates. These tools, combined with precise control of IP configurations, enabled a comprehensive cross-check of beamâbeam models.
Some systematic limitations were identified. During one of the tests, beam stability during the scans was insufficient to allow for symmetric modeling of luminosity decay in both transverse planes. In several scans, the nominal separation steps were based on overestimated beam sizes, resulting in larger-than-expected scan intervals and less sensitivity at low separation. Beamâbeam deflection measurements showed sign reversals at high separation, attributed to nonlinear orbit distortions from hysteresis effects in superconducting magnets. Phase advance errors between IPs and diagnostic systems further reduced the accuracy of simulations in some configurations.
Targeted improvements would allow for higher precision and new insights. Finer separation steps, especially below
\(3\,\sigma \)
, would improve sensitivity to nonlinear beamâbeam effects and allow more accurate comparison with models. Suppressing coherent beamâbeam modes via phase advance tuning of not only two but all three IPs would improve tune shift measurements. Moreover, such optimal phase advance configuration between the IPs can be exploited to enhance luminosity. This effect could be directly propagated into at least a few percent increase in the total collected integrated luminosity, as discussed in Ref. [
36
]. Aligning the phase advance between IPs and diagnostic devices, particularly the BSRT, would reduce sensitivity to optics errors, enabling better beam size measurements for both beams and transverse planes. Moreover, the suppression of the coherent beamâbeam modes facilitates the observation and analysis of the central, incoherent part of the tune spectrum, which is otherwise obscured by the dominant coherent peaks.
Further experimental strategies should also focus on directly probing the beamâbeam interaction at its point of occurrence. This would require access to the full transverse beam distributions with high spatial resolution and sensitivity, going beyond projected sizes or centroids. Such measurements could directly reveal the nonlinear distortions and tails induced by the beamâbeam force. In parallel, validating correction schemes that parametrize multi-IP beamâbeam effects on luminosity via equivalent tune shifts remains essential. These models underpin practical luminosity corrections in complex fill configurations and would benefit from systematic benchmarking.
References
S. Meer, Calibration of the effective beam height in the ISR. Technical report, CERN, Geneva (1968).
https://cds.cern.ch/record/296752
C. Rubbia, Measurement of the luminosity of
\(p{{\overline{p}}}\)
collider with a (generalized) Van der Meer Method. Technical report, CERN, Geneva (1977).
https://cds.cern.ch/record/1025746
ATLAS collaboration, Luminosity determination in
\(pp\)
collisions at
\(\sqrt{s}=13\)
TeV using the ATLAS detector at the LHC. Eur. Phys. J. C
83
(10), 982 (2023).
https://doi.org/10.1140/epjc/s10052-023-11747-w
.
arXiv:2212.09379
CMS Collaboration, Precision luminosity measurement in proton-proton collisions at sqrts = 13 TeV with the CMS detector. Technical report, CERN, Geneva (2025).
https://cds.cern.ch/record/2940794
W. Herr, T. Pieloni, Beam-Beam Effects. Technical report, CERN (2014). Contribution to the CASâCERN Accelerator School: Advanced Accelerator Physics Course, Trondheim, Norway, 18â29 Aug 2013, p. 29.
https://doi.org/10.5170/CERN-2014-009.431
.
https://cds.cern.ch/record/1982430
A.W. Chao, Beam-beam instability. AIP Conf. Proc.
127
, 201â242 (1985).
https://doi.org/10.1063/1.35187
Article
Â
ADS
Â
Google Scholar
Â
E.D. Courant, H.S. Snyder, Theory of the alternating-gradient synchrotron. Ann. Phys.
3
, 1â48 (1958).
https://doi.org/10.1016/0003-4916(58)90012-5
Article
Â
ADS
Â
Google Scholar
Â
P. Bambade, R. Erickson, W.A. Koska, W. Kozanecki, N. Phinney, S.R. Wagner, Observation of beam-beam deflections at the interaction point of the SLAC Linear Collider. Phys. Rev. Lett.
62
(25), 2949â2952 (1989).
https://doi.org/10.1103/PhysRevLett.62.2949
Article
Â
ADS
Â
Google Scholar
Â
W.A. Koska, P. Bambade, W. Kozanecki, N. Phinney, S.R. Wagner, Beam-beam deflection as a beam tuning tool at the slac linear collider. Nucl. Instrum. Methods A
286
, 32 (1990).
https://doi.org/10.1016/0168-9002(90)90203-I
Article
Â
ADS
Â
Google Scholar
Â
C. Bovet, M.D. Hildreth, M. Lamont, H. Schmickler, J. Wenninger, Luminosity optimisation using beam-beam deflections at LEP, in
Conf. Proc. C
(1996).
https://cds.cern.ch/record/306910
D. Brandt, W. Herr, M. Meddahi, A. Verdier, Is lep beam-beam limited at its highest energy? in
Proceedings of the 1999 Particle Accelerator Conference (PAC â99), New York
, pp. 3005â3007. IEEE/JACoW. THP25; CERN report CERN-SL-99-030-AP (1999).
https://accelconf.web.cern.ch/p99/PAPERS/THP25.pdf
W. Kozanecki, T. Pieloni, J. Wenninger, Observation of Beam-beam Deflections with LHC Orbit Data. Technical report, CERN (2013).
https://cds.cern.ch/record/1581723?ln=en
A. Babaev, T. Barklow, O. Karacheban, W. Kozanecki, I. Kralik, A. Mehta, G. Pasztor, T. Pieloni, D. Stickland, C. Tambasco, R. Tomas, J. WaĆczyk, Impact of beam-beam effects on absolute luminosity calibrations at the CERN large hadron collider. Eur. Phys. J. C
84
, 17 (2024).
https://doi.org/10.48550/arXiv.2306.10394
Article
Â
ADS
Â
Google Scholar
Â
T. Pieloni, A study of beam-beam effects in hadron colliders with a large number of bunches. PhD Thesis, Ăcole Polytechnique FĂ©dĂ©rale de Lausanne (2008).
https://doi.org/10.5075/epfl-thesis-4211
ATLAS Collaboration, Improved luminosity determination in pp collisions at
\(\sqrt{s} = 7\)
TeV using the ATLAS detector at the LHC. Eur. Phys. J. C
73
, 2518 (2013).
https://doi.org/10.1140/epjc/s10052-013-2518-3
C.M.S. Collaboration, Development of the cms detector for the cern lhc run 3. J. Instrum.
19
(05), 05064 (2024).
https://doi.org/10.1088/1748-0221/19/05/p05064
Article
Â
Google Scholar
Â
C. Barschel, M. Ferro-Luzzi, J.-J. Gras, M. Ludwig, P. Odier, S. Thoulet, Results of the LHC DCCT calibration studies. Technical Report CERN-ATS-Note-2012-026 PERF (2012).
https://cds.cern.ch/record/1425904?ln=en
G. Anders et al., Study of the relative LHC bunch populations for luminosity calibration. Technical Report. CERN-ATS-Note-2012-028 PERF (2012).
https://cds.cern.ch/record/1427726
G. Trad, Development and Optimisation of the SPS and LHC Beam Diagnostics Based on Synchrotron Radiation Monitors Presented 22 Jan 2015. Presented 22 Jan 2015.
https://cds.cern.ch/record/2266055
R. TomĂĄs, O. BrĂŒning, M. Giovannozzi, P. Hagen, M. Lamont, F. Schmidt, G. Vanbavinckhove, M. Aiba, R. Calaga, R. Miyamoto, Cern large hadron collider optics model, measurements, and corrections. Phys. Rev. ST Accel. Beams
13
, 121004 (2010).
https://doi.org/10.1103/PhysRevSTAB.13.121004
R. TomĂĄs, M. Aiba, A. Franchi, U. Iriso, Review of linear optics measurement and correction for charged particle accelerators. Phys. Rev. Accel. Beams
20
, 054801 (2017).
https://doi.org/10.1103/PhysRevAccelBeams.20.054801
Article
Â
ADS
Â
Google Scholar
Â
O. Karacheban, P. Tsrunchev, on behalf of CMS, Emittance scans for cms luminosity calibration. EPJ Web Conf.
201
, 04001 (2019).
https://doi.org/10.1051/epjconf/201920104001
F. Carlier, R. TomĂĄs, Accuracy and feasibility of the
\({\beta }^{*}\)
measurement for lhc and high luminosity lhc using
\(k\)
modulation. Phys. Rev. Accel. Beams
20
, 011005 (2017).
https://doi.org/10.1103/PhysRevAccelBeams.20.011005
Article
Â
ADS
Â
Google Scholar
Â
M. SöderĂ©n, J. Komppula, G. Kotzian, S. Rains, D. Valuch, ADT and Obsbox in LHC Run 2, plans for LS2, pp. 165â171 (2019)
A. Boccardi, M. Gasior, R. Jones, P. Karlsson, R. Steinhagen, First results from the lhc bbq tune and chromaticity systems. Technical report, CERN, Geneva (2009).
https://cds.cern.ch/record/1156349
M. Gasior, G. Baud, J. Olexa, G. Valentino, First Operational Experience with the LHC Diode ORbit and OScillation (DOROS) System, p. 07 (2017)
https://doi.org/10.18429/JACoW-IBIC2016-MOPG07
A.W. Chao, Coherent beam-beam effects, in
Frontiers of Particle Beams: Intensity Limitations
(Springer, Berlin, Heidelberg, 1992), pp. 363â414.
https://doi.org/10.1007/3-540-55250-2_36
J. WaĆczyk, A. Dabrowski, D. Stickland, T. Pieloni, W. Kozanecki, X. Buffat, Impact of multiple beam-beam encounters on LHC absolute-luminosity calibrations by the van der Meer method. JACoW IPAC
2023
, 053 (2023).
https://doi.org/10.18429/JACoW-IPAC2023-WEPA053
Article
Â
Google Scholar
Â
J. WaĆczyk, Precision luminosity measurement at hadron colliders. PhD Thesis, Ăcole Polytechnique FĂ©dĂ©rale de Lausanne (2024).
https://doi.org/10.5075/epfl-thesis-10500
T. Pieloni, J. Barranco, X. Buffat, P.C. Jorge, L.M. Medrano, C. Tambasco, R. TomĂĄs, Dynamic beta and beta-beating effects in the presence of the beam-beam interactions in
57th ICFA Advanced Beam Dynamics Workshop on High-Intensity and High-Brightness Hadron Beams
, p. 027 (2016).
https://doi.org/10.18429/JACoW-HB2016-MOPR027
F. Carlier, et al. Challenges of K-modulation measurements in the LHC Run 3, in
Proceedings of IPACâ23. IPACâ23â14th International Particle Accelerator Conference
(JACoW Publishing, Geneva, Venice, 2023), pp. 531â534.
https://doi.org/10.18429/JACoW-IPAC2023-MOPL014
Y. Alexahin, A study of the coherent beam-beam effect in the framework of the Vlasov perturbation theory. Nucl. Instrum. Methods Phys. Res., Sect. A
480
(1â3), 253â288 (2002).
https://doi.org/10.1016/S0168-9002(01)01219-0
Article
Â
ADS
Â
Google Scholar
Â
R. Calaga, W. Herr, G. Papotti, T. Pieloni, X. Buffat, S. White, R. Giachino, Coherent beam-beam mode in the LHC. CERN (2014).
https://doi.org/10.5170/CERN-2014-004.227
.
http://cds.cern.ch/record/1957039
K. Yokoya, Y. Funakoshi, E. Kikutani, H. Koiso, J. Urakawa, Tune shift of coherent beam-beam oscillations. Part. Accel.
27
, 181â186 (1990)
Google Scholar
Â
M. Bassetti, G.A. Erskine, Closed expression for the electrical field of a two-dimensional gaussian charge. Technical report, CERN, Geneva (1980).
https://cds.cern.ch/record/122227
J. WaĆczyk, T. Pieloni, Optimizing beam-beam beta-beating for luminosity enhancement at the LHC. JACoW IPAC
2025
, 018 (2025)
Google Scholar
Â
Download references |
| Markdown | [Skip to main content](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#main)
Advertisement
[](https://pubads.g.doubleclick.net/gampad/jump?iu=/270604982/springerlink/10052/article&sz=728x90&pos=top&articleid=s10052-026-15476-8)
[](https://link.springer.com/)
[Account](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8)
[Menu](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#eds-c-header-nav)
[Find a journal](https://link.springer.com/journals/) [Publish with us](https://www.springernature.com/gp/authors) [Track your research](https://link.springernature.com/home/)
[Search](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#eds-c-header-popup-search)
[Saved research](https://link.springer.com/saved-research)
[Cart](https://order.springer.com/public/cart)
1. [Home](https://link.springer.com/)
2. [The European Physical Journal C](https://link.springer.com/journal/10052)
3. Article
# First measurements of beam-beam effects in beam-separation, luminosity-calibration scans at the LHC
- Regular Article - Experimental Physics
- Experimental Physics
- [Open access](https://www.springernature.com/gp/open-science/about/the-fundamentals-of-open-access-and-open-research)
- Published:
17 March 2026
- Volume 86, article number 282, (2026)
- [Cite this article](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#citeas)
You have full access to this [open access](https://www.springernature.com/gp/open-science/about/the-fundamentals-of-open-access-and-open-research) article
[Download PDF](https://link.springer.com/content/pdf/10.1140/epjc/s10052-026-15476-8.pdf)
[Save article](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/save-research?_csrf=skr420D_Tsy443125TKQwosqfCsJrZRn)
[View saved research](https://link.springer.com/saved-research)
[ The European Physical Journal C](https://link.springer.com/journal/10052)
[Aims and scope](https://link.springer.com/journal/10052/aims-and-scope)
[Submit manuscript](https://mc.manuscriptcentral.com/epjc)
First measurements of beam-beam effects in beam-separation, luminosity-calibration scans at the LHC
[Download PDF](https://link.springer.com/content/pdf/10.1140/epjc/s10052-026-15476-8.pdf)
- [Joanna WaĆczyk](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#auth-Joanna-Wa_czyk-Aff1-Aff2)
[ORCID: orcid.org/0000-0002-8562-1863](https://orcid.org/0000-0002-8562-1863)[1](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Aff1),[2](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Aff2),
- [Tatiana Pieloni](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#auth-Tatiana-Pieloni-Aff2) [ORCID: orcid.org/0000-0002-3218-0048](https://orcid.org/0000-0002-3218-0048)[2](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Aff2),
- [Xavier Buffat](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#auth-Xavier-Buffat-Aff1) [ORCID: orcid.org/0000-0002-2338-8808](https://orcid.org/0000-0002-2338-8808)[1](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Aff1),
- [Anne Dabrowski](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#auth-Anne-Dabrowski-Aff1) [ORCID: orcid.org/0000-0003-2570-9676](https://orcid.org/0000-0003-2570-9676)[1](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Aff1),
- [David Stickland](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#auth-David-Stickland-Aff3) [ORCID: orcid.org/0000-0003-4702-8820](https://orcid.org/0000-0003-4702-8820)[3](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Aff3) &
- âŠ
- [Rogelio Tomas](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#auth-Rogelio-Tomas-Aff1) [ORCID: orcid.org/0000-0002-9857-1703](https://orcid.org/0000-0002-9857-1703)[1](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Aff1)
Show authors
## Abstract
At the CERN Large Hadron Collider (LHC), absolute luminosity calibrations obtained by the van der Meer (*vdM*) method are affected by the mutual electromagnetic interaction of the two beams. In the last few years, increasingly stringent, physics-driven requirements on the accuracy of the absolute luminosity scale, combined with the unparalleled precision of the luminosity measurements achieved by the LHC experiments, motivated comprehensive simulation studies to accurately model beamâbeam effects and provide simulation-based corrections to *vdM*\-based luminosity calibrations. This paper reports the first attempt at an experimental validation of that model, in dedicated accelerator measurements under controlled conditions, of the predicted impact of beam-beam effects on the luminosity measured at the LHC under luminosity-calibration conditions. This is the first measurement of beam-beam interaction effects on luminosity at the LHC. The results show good agreement with the predictions of multiparticle simulations, and justify the recently proposed strategy for correcting beam-beam biases on absolute luminosity calibrations at hadron colliders. This marks a critical step forward in precise luminosity calibration for past, current, and future datasets, including the High-Luminosity Large Hadron Collider era.
### Explore related subjects
Discover the latest articles, books and news in related subjects, suggested using machine learning.
- [Accelerator Physics](https://link.springer.com/subjects/accelerator-physics)
- [Experimental Particle Physics](https://link.springer.com/subjects/experimental-particle-physics)
- [Measurement Science and Instrumentation](https://link.springer.com/subjects/measurement-science-and-instrumentation)
- [Metrology](https://link.springer.com/subjects/metrology)
- [Optical Metrology](https://link.springer.com/subjects/optical-metrology)
- [Particle Physics](https://link.springer.com/subjects/particle-physics)
## 1 Introduction
A precise determination of the absolute luminosity scale is essential for a wide range of LHC measurements, as it directly affects the normalization of many key physics cross sections. At the LHC, this calibration relies primarily on the van der Meer (vdM) method \[[1](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR1 "S. Meer, Calibration of the effective beam height in the ISR. Technical report, CERN, Geneva (1968).
https://cds.cern.ch/record/296752
"), [2](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR2 "C. Rubbia, Measurement of the luminosity of
$$p{{\overline{p}}}$$
p
p
ÂŻ
collider with a (generalized) Van der Meer Method. Technical report, CERN, Geneva (1977).
https://cds.cern.ch/record/1025746
")\], which uses dedicated beam-separation scans performed under specially tailored beam conditions to relate measured interaction rates to the absolute luminosity inferred from beam parameters. While sub-percent uncertainties on the integrated luminosity have been achieved for Run-2 data by ATLAS \[[3](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR3 "ATLAS collaboration, Luminosity determination in
$$pp$$
pp
collisions at
$$\sqrt{s}=13$$
s
=
13
TeV using the ATLAS detector at the LHC. Eur. Phys. J. C 83(10), 982 (2023).
https://doi.org/10.1140/epjc/s10052-023-11747-w
.
arXiv:2212.09379
")\] and CMS \[[4](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR4 "CMS Collaboration, Precision luminosity measurement in proton-proton collisions at sqrts = 13 TeV with the CMS detector. Technical report, CERN, Geneva (2025).
https://cds.cern.ch/record/2940794
")\], further improvements and robust experimental validation of beam-dynamical effects remain essential.
For relativistic proton beams colliding with zero crossing angle, which is the typical configuration for the *vdM* calibration, the luminosity at a given interaction point (IP) of the LHC is proportional to the overlap integral of the particle-density distributions \\(\\rho \_{\\text {1,i}}\\), \\(\\rho \_{\\text {2,i}}\\) in bunch-pair *i* \[[5](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR5 "W. Herr, T. Pieloni, Beam-Beam Effects. Technical report, CERN (2014). Contribution to the CASâCERN Accelerator School: Advanced Accelerator Physics Course, Trondheim, Norway, 18â29 Aug 2013, p. 29.
https://doi.org/10.5170/CERN-2014-009.431
.
https://cds.cern.ch/record/1982430
")\]:
\$\$\\begin{aligned} {\\mathcal {L}}\_{inst}&= 2c f\_{\\text {rev}} \\sum ^{N\_b}\_i n\_{\\text {1,i}}, n\_{\\text {2,i}} \\nonumber \\\\&\\quad \\times \\iiiint ^{+\\infty }\_{-\\infty }\\rho \_{\\text {1,i}}(x,y,z-ct) \\nonumber \\\\&\\quad \\times \\rho \_{\\text {2,i}}(x,y,z+ct)\\,dxdydzdt, \\end{aligned}\$\$
(1)
where *c* is the speed of light, \\(f\_{\\text {rev}}\\) is the LHC revolution frequency, \\(N\_b\\) is the number of colliding bunches, and \\(n\_{\\text {1,i}}, n\_{\\text {2,i}}\\) are the corresponding total charges per bunch for each colliding pair. Although only a very small fraction of the beam particles actually collide, the two opposing beams interact with each other electromagnetically, a dynamical process known as the beam-beam interaction \[[5](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR5 "W. Herr, T. Pieloni, Beam-Beam Effects. Technical report, CERN (2014). Contribution to the CASâCERN Accelerator School: Advanced Accelerator Physics Course, Trondheim, Norway, 18â29 Aug 2013, p. 29.
https://doi.org/10.5170/CERN-2014-009.431
.
https://cds.cern.ch/record/1982430
"), [6](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR6 "A.W. Chao, Beam-beam instability. AIP Conf. Proc. 127, 201â242 (1985).
https://doi.org/10.1063/1.35187
")\]. In the LHC, as in any other synchrotron, the beams are not continuous but are divided into discrete âbunchesâ, each with a high proton density, containing approximately \\(10^{11}\\) protons over a bunch length of 7.5â9 cm (in terms of RMS). This electromagnetic interaction occurs when the two beams share a common beam pipe, and is governed by a non-linear force that depends on the radial distance *r* of a test particle to the center of the opposing âsourceâ bunch with *n* protons Gaussian-distributed in the transverse planes (valid approximation for the LHC beams):
\$\$\\begin{aligned} F = - \\frac{ne^2}{2\\pi \\epsilon \_0r}\\biggl (1-\\exp \\biggl \[ -\\frac{r^2}{2\\sigma ^2}\\biggr \] \\biggr ), \\end{aligned}\$\$
(2)
where *e* is the elementary proton charge, \\(\\epsilon \_0\\) is the vacuum permittivity, and \\(\\sigma \\) is the RMS radius of the transverse charge distribution within the bunch, in a simplified case of a round shape \\((\\sigma = \\sigma \_x = \\sigma \_y)\\). At very small amplitudes \\((r\\rightarrow 0)\\) the force can be linearized giving the expression for the so called beam-beam parameter, which is often used to assess the strength of the force:
\$\$\\begin{aligned} \\xi = \\frac{nr\_p\\beta ^\*}{4\\pi \\gamma \\sigma ^2} = \\frac{nr\_p}{4\\pi \\epsilon \_n}, \\end{aligned}\$\$
(3)
where \\(r\_p\\) is the classical proton radius, \\(\\beta ^\*\\) is the optical beta-function at the interaction point (IP) â one of the CourantâSnyder parameters \[[7](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR7 "E.D. Courant, H.S. Snyder, Theory of the alternating-gradient synchrotron. Ann. Phys. 3, 1â48 (1958).
https://doi.org/10.1016/0003-4916(58)90012-5
")\], \\(\\gamma \\) is the Lorentz factor, and \\(\\epsilon \_n\\) is the normalized emittance.
Beam-beam effects were first observed almost 40 years ago \[[8](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR8 "P. Bambade, R. Erickson, W.A. Koska, W. Kozanecki, N. Phinney, S.R. Wagner, Observation of beam-beam deflections at the interaction point of the SLAC Linear Collider. Phys. Rev. Lett. 62(25), 2949â2952 (1989).
https://doi.org/10.1103/PhysRevLett.62.2949
")\] at the lepton colliders, and have been used extensively for accelerator diagnostics and optimization at both SLC and LEP \[[9](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR9 "W.A. Koska, P. Bambade, W. Kozanecki, N. Phinney, S.R. Wagner, Beam-beam deflection as a beam tuning tool at the slac linear collider. Nucl. Instrum. Methods A 286, 32 (1990).
https://doi.org/10.1016/0168-9002(90)90203-I
"),[10](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR10 "C. Bovet, M.D. Hildreth, M. Lamont, H. Schmickler, J. Wenninger, Luminosity optimisation using beam-beam deflections at LEP, in Conf. Proc. C (1996).
https://cds.cern.ch/record/306910
"),[11](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR11 "D. Brandt, W. Herr, M. Meddahi, A. Verdier, Is lep beam-beam limited at its highest energy? in Proceedings of the 1999 Particle Accelerator Conference (PAC â99), New York, pp. 3005â3007. IEEE/JACoW. THP25; CERN report CERN-SL-99-030-AP (1999).
https://accelconf.web.cern.ch/p99/PAPERS/THP25.pdf
")\]. At the LHC, they were first observed in late Run 1 \[[12](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR12 "W. Kozanecki, T. Pieloni, J. Wenninger, Observation of Beam-beam Deflections with LHC Orbit Data. Technical report, CERN (2013).
https://cds.cern.ch/record/1581723?ln=en
")\], but only in terms of the beam-beam induced coherent deflections. Single particle effects, in contrast, were presumed too small in hadron colliders that until recently, there was no motivation to measure them with a precision that would be meaningful in the context evoked here.
Nonetheless, the full complexity of this interaction must be accounted for when performing the luminosity calibration via the vdM method, using the separation scans. These are designed to measure the detector-specific constant \\(\\sigma \_{vis}\\) that relates the observed rate \\(\\mu ^{vis}\\) to the absolute instantaneous luminosity \\({\\mathcal {L}}\_{inst}\\):
\$\$\\begin{aligned} \\mu ^{vis} = \\frac{{\\mathcal {L}}\_{inst}\\sigma \_{vis}}{f\_{rev}}. \\end{aligned}\$\$
(4)
Under the assumption of uncorrelated particle densities in *x* and *y* planes, the transverse convolved bunch widths \\(\\Sigma \_x, \\Sigma \_y\\) can be extracted from the measured beam-separation dependence of the collision rate in the corresponding direction. The combined information from these scans, and bunch intensities can be used to calculate the instantaneous luminosity at the head-on position. The proportionality between the measured rate at the centered head-on position \\(\\mu \_{pk}\\) and the reconstructed luminosity from measured bunch parameters during the vdM calibration defines the luminometer-specific visible cross-section:
\$\$\\begin{aligned} \\sigma \_{vis} = \\frac{2\\pi \\Sigma \_x\\Sigma \_y}{n\_1 n\_2} \\mu \_{pk}. \\end{aligned}\$\$
(5)
If left uncorrected, the \\(\\sigma \_{vis}\\) measurement by the *vdM* method is biased by the beam-separation dependence of the mutual electromagnetic interaction of the two beams: the colliding bunches experience deflection-induced orbit shifts, as well as optical distortions akin to the dynamic-\\(\\beta \\) effect, that both depend on the transverse beam separation and must therefore be accounted for when deriving the absolute luminosity scale. The correction strategy discussed extensively in Ref. \[[13](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR13 "A. Babaev, T. Barklow, O. Karacheban, W. Kozanecki, I. Kralik, A. Mehta, G. Pasztor, T. Pieloni, D. Stickland, C. Tambasco, R. Tomas, J. WaĆczyk, Impact of beam-beam effects on absolute luminosity calibrations at the CERN large hadron collider. Eur. Phys. J. C 84, 17 (2024).
https://doi.org/10.48550/arXiv.2306.10394
")\] relies on simulation studies only, and thus requires experimental validation to establish its reliability.
The present report summarizes a campaign of dedicated measurements carried out in the spring of 2022 at the LHC, aimed at confirming the accuracy of the simulation studies detailed in Ref. \[[13](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR13 "A. Babaev, T. Barklow, O. Karacheban, W. Kozanecki, I. Kralik, A. Mehta, G. Pasztor, T. Pieloni, D. Stickland, C. Tambasco, R. Tomas, J. WaĆczyk, Impact of beam-beam effects on absolute luminosity calibrations at the CERN large hadron collider. Eur. Phys. J. C 84, 17 (2024).
https://doi.org/10.48550/arXiv.2306.10394
")\], and at validating the beam-beam correction strategy used in the luminosity-calibration analyses of the ALICE, ATLAS, CMS and LHCb experiments. More specifically, the goal is to quantify the impact of the beam-beam interaction on the tune spectra, orbit, the transverse-density distribution of the colliding bunches, as well as, for the first time, on the luminosity, by systematically varying the strength or number of the beam-beam interactions. The measurements are repeated for different values of the beam-beam parameter to verify the scaling proposed in \[[13](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR13 "A. Babaev, T. Barklow, O. Karacheban, W. Kozanecki, I. Kralik, A. Mehta, G. Pasztor, T. Pieloni, D. Stickland, C. Tambasco, R. Tomas, J. WaĆczyk, Impact of beam-beam effects on absolute luminosity calibrations at the CERN large hadron collider. Eur. Phys. J. C 84, 17 (2024).
https://doi.org/10.48550/arXiv.2306.10394
")\], for various choices of the scanning IP and for several multi-IP configurations. The beam conditions are chosen to be representative of *vdM* calibration sessions at the LHC, but optimized so as to maximize both the sensitivity of the measurements, and the operational efficiency.
This paper is organized as follows. The measurement strategy is detailed in Sect. [2](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec2), followed by the main ingredients of this accelerator experiment: beam conditions (Sect. [3](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec3)), luminometers and accelerator instrumentation (Sect. [4](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec4)), and optimization of the ring lattice to maximize the sensitivity of the measurements to beamâbeam-induced optical distortions of the colliding bunches (Sect. [5](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec11)). Section [6](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec19) is devoted to the experimental characterization of the impact of beam-beam effects on the tune spectra, on the luminosity at different IPs, as well as on other observables such as transverse single-bunch widths and IP orbits. At each step, the results are confronted with the predictions of the COherent Multibunch Beam-beam Interactions (COMBI) tracking code \[[14](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR14 "T. Pieloni, A study of beam-beam effects in hadron colliders with a large number of bunches. PhD Thesis, Ăcole Polytechnique FĂ©dĂ©rale de Lausanne (2008).
https://doi.org/10.5075/epfl-thesis-4211
")\], that produced most of the results detailed in Ref. \[[13](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR13 "A. Babaev, T. Barklow, O. Karacheban, W. Kozanecki, I. Kralik, A. Mehta, G. Pasztor, T. Pieloni, D. Stickland, C. Tambasco, R. Tomas, J. WaĆczyk, Impact of beam-beam effects on absolute luminosity calibrations at the CERN large hadron collider. Eur. Phys. J. C 84, 17 (2024).
https://doi.org/10.48550/arXiv.2306.10394
")\]. A collection of macroparticles is simulated turn by turn, based on a simplified model of the accelerator lattice using linear transfer and the self-consistent computation of the beam-beam forces at the IP. The conclusions are presented in Sect. [7](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec33).
## 2 Measurement strategy
The beam-beam interaction at an interaction point, referred to here as the *scanning* IP, induces measurable changes in key beam parameters such as the tune spectrum, closed orbit, and transverse single-beam profiles. These effects can be observed, albeit with limited precision, using standard accelerator diagnostics by comparing measurements obtained under varying transverse beam separation at the scanning IP. However, the resulting modifications to the transverse density distribution of the colliding bunches remain below the sensitivity threshold of conventional instrumentation, including synchrotron-radiation-based beam-profile monitors. Beam-beam-induced variations typically lie in the 0.5â1% range, small enough that only a highly precise measurement of collision-rate changes can provide adequate sensitivity.
**Fig. 1**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/1)
Scheme of the LHC and its interaction points (IP) with indicated IP1 as witness IP while performing a scan at IP5
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/1)
Importantly, the beam-beam interaction induces changes in beam properties that are intrinsically entangled with the *pp* collisions when observing luminosity at the scanning IP. Consequently, under realistic collider conditions, it is not possible to directly access a reference luminosity signal that is free of beam-beam effects â unlike in simulation studies, where such effects can be explicitly disabled. In a multi-collision configuration, however, luminosity shifts at one or more other IPs (Fig. [1](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig1)), where the beams are kept in continuous head-on collision and therefore act as non-scanning, or *witness*, IPs, can be used to monitor the separation-dependent beamâbeam effects. These effects are induced at the scanning IP, where the beams are deliberately brought in and out of collision, and propagate around the rings to the witness IPs, where they manifest as changes in the measured luminosity. In principle, any of the four LHC IPs can serve as a scanning IP, and any non-scanning IP can be designated as the witness IP for a given measurement. In practice, instrumental and operational constraints at the time the experiment was carried out restrict the choice of witness IP to the ATLAS (IP1) and CMS (IP5) collision points; IPs 1, 2 or 5, or combinations thereof, are used as scanning IP(s). Collisions at IP8 were deliberately avoided because the large longitudinal offset of the LHCb collision point, that breaks the eight-fold symmetry of the LHC rings, precludes the possibility of both members of a colliding-bunch pair to collide at all four IPs. In addition, the large crossing angle at IP8 would complicate the interpretation of the beam-beam effects that would occur at that IP if some of the bunches collided there.
The propagation of the beamâbeam-induced, amplitude-dependent \\(\\beta \\)\-beating from the scanning to a witness IP is controlled by the betatron phase advance between these two IPs. Since the targeted signatures are delicate to measure at best, it is natural to try and enhance them by adjusting this phase advance so as to maximize the sensitivity of luminosity shifts at the witness IP to beam-beam effects at the scanning IP; in doing so, however, the overall tunes must be preserved. Given the central and symmetric roles played by IP1 and IP5 in this experiment, the adjustment of their relative phase advance drove the optimization procedure that is detailed in Sect. [5](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec11), and that achieved a threefold improvement in measurement sensitivity.
Central to the quality of the measurements is the controlled variation of the beam-beam parameter at the scanning IP. This can be achieved either:
- by using *step scans*, *i.e.* by fully separating the beams at the scanning IP in either the horizontal or the vertical plane, and then bringing them back into head-on collision; or
- by using *separation scans*, *i.e.* by scanning the beams transversely with respect to each other in either the horizontal or the vertical plane; or
- by taking advantage of the natural beam-intensity decay and emittance growth to progressively reduce the beam-beam parameter.
Equally important is to ensure that the beams remain in head-on collision at the witness IP(s), such that reproducible luminosity shifts measured at these locations can be unambiguously correlated with controlled changes in the strength of the beam-beam interaction at the scanning IP. It was verified that the beam-beam-induced orbit shift and the potential non-closure of the orbit bumps used to control the beam separation at the scanning IP, did not significantly affect the actual separation, and thereby the measured collision rate, at the witness IP. In a few cases however, orbit drifts of uncontrolled origin ended up degrading the quality of some of the measurements.
Simulations demonstrate \[[13](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR13 "A. Babaev, T. Barklow, O. Karacheban, W. Kozanecki, I. Kralik, A. Mehta, G. Pasztor, T. Pieloni, D. Stickland, C. Tambasco, R. Tomas, J. WaĆczyk, Impact of beam-beam effects on absolute luminosity calibrations at the CERN large hadron collider. Eur. Phys. J. C 84, 17 (2024).
https://doi.org/10.48550/arXiv.2306.10394
")\] that multi-IP effects strongly influence the magnitude of beam-beam biases to vdM calibrations. Their characterization, therefore, constitutes an essential component of this experiment. The scan protocol is a generalization of that in the two-IP case. One IP (say, IP1) is designated as the witness IP, with beams colliding head-on there and at IP2, and IP5; the beams are then taken out of collision at IP2 only, then at both IP2 and IP5, and finally returned to a three-collision configuration in the reverse order. Since the phase-advance combinations are different in each configuration, the corresponding luminosity shifts predicted at the witness IP are also different and can be meaningfully confronted with the data.
## 3 Beam conditions
In the round-beam approximation, the beam-beam parameter \\(\\xi \\) does not depend on the beam energy (Eq. ([3](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Equ3))). Therefore, and in order to minimize the overhead associated with commissioning and operating the LHC in a non-standard optical configuration, as well as to allow for frequent refilling with different bunch patterns, the measurements were carried out at injection energy (\\(450\\ \\text {GeV}\\) per beam), and under non-standard machine-protection conditions. The improvement in operational efficiency came at a cost:
- the counting rate of the luminometers (or equivalently their visible cross-section) is about an order of magnitude smaller at \\(\\sqrt{s} = 900\\ \\text {GeV}\\) than at \\(13\\ \\text {TeV}\\), reflecting the combination of a 40% drop in the inelastic *pp* cross-section, a factor of two to three drop in the multiplicity of the final-state particles, and a significant softening of their momentum spectrum;
- because of the total-intensity constraints dictated by machine-protection requirements, the injected beam could not exceed 4 bunches per beam (compared to 150 during a routine *vdM*\-calibration session), with a maximum allowed population of \\(1.25\\times 10^{11}~p\\)/bunch;
- the combination of the low beam energy, that suppresses synchrotron-radiation damping, and of the large bunch intensity, that enhances intra-beam scattering, resulted in relatively rapid emittance growth and rather short single-beam lifetimes. To partially mitigate these effects, all collision-rate measurements are expressed in terms of specific luminosity, thereby automatically accounting for the beam-intensity decay.
The bunch intensity was deliberately increased by about 25% beyond the *vdM*\-scan values typical of LHC Run 2, to maximize both the luminosity and the beam-beam parameter \\(\\xi \\). The latter reached 0.010 per IP (Table [1](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Tab1)), compared to 0.003â0.006 during normal *pp* *vdM* sessions. With a target injected emittance of \\(1.5 \\ \\upmu \\text {m} \\cdot \\text {rad}\\), the \\(\\beta \\) function at the IP set to \\(\\beta ^\* = 11\\ \\text {m}\\), and zero crossing angle at IP1 and IP5, the statistical uncertainty affecting a typically 60-second long per-bunch luminosity measurement in the presence of head-on collisions, lay around 0.5%.
**Table 1 Range of beam-beam parameter values during the various stages of the experiment**
[Full size table](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/tables/1)
The bunch patterns were chosen such that all bunches in each beam collided either at IP1 and IP5 only, or at all three of IP1, 2 and 5. The orbit-stabilization feedback system was turned off during the data-taking periods to prevent it from interfering with the beamâbeam-induced orbit shift. The chromaticity was set to its standard value in physics fills, of +10 units, to guarantee coherent stability against the machine impedance, and in view of the very small number of bunches and of their large longitudinal spacing, the beam-stabilizing Landau octupoles were set to their minimum current (1A). The damping time of the bunch-by-bunch transverse feedback was set to 1000 turns, a rather loose setting intended to preserve longer the natural beam oscillation and therefore improve the precision of the tune measurements.
## 4 Beam instrumentation
### 4\.1 Luminometers
The luminometer systems in use by the ATLAS and CMS collaborations at IP1 and IP5 are described in Refs. \[[15](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR15 "ATLAS Collaboration, Improved luminosity determination in pp collisions at
$$\sqrt{s} = 7$$
s
=
7
TeV using the ATLAS detector at the LHC. Eur. Phys. J. C 73, 2518 (2013).
https://doi.org/10.1140/epjc/s10052-013-2518-3
"), [16](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR16 "C.M.S. Collaboration, Development of the cms detector for the cern lhc run 3. J. Instrum. 19(05), 05064 (2024).
https://doi.org/10.1088/1748-0221/19/05/p05064
")\] respectively; no luminosity measurements were available at IP2 during the collider experiment described in this paper. Maximizing the statistical sensitivity leads to choosing the luminosity algorithm with the highest possible acceptance, namely:
- the hit rate per bunch crossing in the ATLAS Minimum Bias Trigger Scintillators (MBTS) \[[15](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR15 "ATLAS Collaboration, Improved luminosity determination in pp collisions at
$$\sqrt{s} = 7$$
s
=
7
TeV using the ATLAS detector at the LHC. Eur. Phys. J. C 73, 2518 (2013).
https://doi.org/10.1140/epjc/s10052-013-2518-3
")\], and
- the occupancy in the CMS Hadron Forward Calorimeter (HFOC) \[[16](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR16 "C.M.S. Collaboration, Development of the cms detector for the cern lhc run 3. J. Instrum. 19(05), 05064 (2024).
https://doi.org/10.1088/1748-0221/19/05/p05064
")\].
In both cases, raw luminometer counts are converted to collision rates using the Poisson formalism \[[15](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR15 "ATLAS Collaboration, Improved luminosity determination in pp collisions at
$$\sqrt{s} = 7$$
s
=
7
TeV using the ATLAS detector at the LHC. Eur. Phys. J. C 73, 2518 (2013).
https://doi.org/10.1140/epjc/s10052-013-2518-3
"), [16](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR16 "C.M.S. Collaboration, Development of the cms detector for the cern lhc run 3. J. Instrum. 19(05), 05064 (2024).
https://doi.org/10.1088/1748-0221/19/05/p05064
")\], and the instantaneous luminosity is averaged over typically 60-second time bins, that during the scans are synchronized with the scan steps. The statistical uncertainties are estimated from either the total number of raw luminometer counts per time bin (ATLAS), or from the RMS of the approximately 40 luminosity samples recorded in a given time bin (CMS). Since the absolute luminosity scale is irrelevant, and in order to simplify the interpretation of the results, all measurements in this paper are presented in terms of fractional shifts in the bunch-averaged specific luminosity, relative to a reference time that depends on the type of measurement considered.
### 4\.2 Bunch-charge measurement
The total intensity of each beam is measured by a direct-current current transformer (DCCT), the absolute scale of which is calibrated against a very high precision pulse generator \[[17](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR17 "C. Barschel, M. Ferro-Luzzi, J.-J. Gras, M. Ludwig, P. Odier, S. Thoulet, Results of the LHC DCCT calibration studies. Technical Report CERN-ATS-Note-2012-026 PERF (2012).
https://cds.cern.ch/record/1425904?ln=en
")\]. The bunch-by-bunch charge fractions, in turn, normalized to the total stored intensity, are measured, separately for the two beams, by Fast Beam Current Transformers (FBCT) \[[18](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR18 "G. Anders et al., Study of the relative LHC bunch populations for luminosity calibration. Technical Report. CERN-ATS-Note-2012-028 PERF (2012).
https://cds.cern.ch/record/1427726
")\]. The latter devices also provide the fill pattern, *i.e.* the relative location, at a given instant and around the two rings, of the nominally filled bunches. Since the absolute scale of the DCCTs is known to much better accuracy than that of the FBCTs, the bunch charges are typically computed as the product of the FBCT bunch-charge fractions and the total circulating beam intensity reported by the corresponding DCCT.
### 4\.3 Emittance measurement
The synchrotron-light beam-profile monitors \[[19](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR19 "G. Trad, Development and Optimisation of the SPS and LHC Beam Diagnostics Based on Synchrotron Radiation Monitors Presented 22 Jan 2015. Presented 22 Jan 2015.
https://cds.cern.ch/record/2266055
")\], dubbed BSRTs (âBeam Synchrotron Radiation Telescopeâ), are mainly used to track emittance evolution over time, but they can also measure the relative changes in transverse single-beam size that result from beamâbeam-induced \\(\\beta \\)\-beating. BSRT data are recorded every second; in what follows, they are presented averaged over one-minute intervals for easier comparison with other measurements.
The absolute length scale of the BSRT profiles suffers from significant uncertainties, if only because the online optical corrections to the images can be updated only a few times per year, and are unable to track the evolution of the efficiency of the light sensors, that depends both on time and on the position of the light spot on the sensor array. For each beam therefore, the BSRT is complemented by a wire-scanner (WS) profile monitor that is located in the same straight section. The advantage of the WS is that its accuracy is significantly better than that of the BSRT; its down side, in the present context, is that the WS cannot acquire data continuously and must be triggered manually. For the results presented in this paper, the wire scanners were flown through the beams at the start of each group of measurements to provide single-beam emittance measurements that are as accurate as possible, and then at regular intervals thereafter. The BSRT is used to interpolate the time evolution of the emittance between two sets of WS measurements.
To interpret measured RMS bunch widths in terms of emittance requires the knowledge of the \\(\\beta \\) functions at the locations of the BSRT and the WS. These are determined by the phase-advance method, and compared in Table [2](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Tab2) to their model value computed using MAD-X; the agreement is typically better than 5%, and the worst disagreement amounts to 9%. These discrepancies are attributed to imperfections in the magnetic lattice. Uncertainties in the measurement mainly arise from the interpolation of the lattice functions between the two closest beam-position monitors (BPMs) where the \\(\\beta \\)\-functions are measured \[[20](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR20 "R. TomĂĄs, O. BrĂŒning, M. Giovannozzi, P. Hagen, M. Lamont, F. Schmidt, G. Vanbavinckhove, M. Aiba, R. Calaga, R. Miyamoto, Cern large hadron collider optics model, measurements, and corrections. Phys. Rev. ST Accel. Beams 13, 121004 (2010).
https://doi.org/10.1103/PhysRevSTAB.13.121004
"), [21](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR21 "R. TomĂĄs, M. Aiba, A. Franchi, U. Iriso, Review of linear optics measurement and correction for charged particle accelerators. Phys. Rev. Accel. Beams 20, 054801 (2017).
https://doi.org/10.1103/PhysRevAccelBeams.20.054801
")\], and either the BSRT or the WS location. Based on the MAD-X LHC lattice model, these uncertainties are estimated not to exceed 3%, and this value is assigned as the systematic uncertainty on the \\(\\beta \\) functions used in BSRT- and WS-based emittance measurements.
**Table 2 \\(\\beta \\) functions at the BSRT and WS locations, for the optical configuration detailed in Sect. [5\.3](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec14)**
[Full size table](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/tables/2)
The beam-averaged emittance, *i.e.* the average of the emittances of the beam-1 bunch and of the corresponding beam-2 bunch, can be obtained directly from emittance scans \[[22](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR22 "O. Karacheban, P. Tsrunchev, on behalf of CMS, Emittance scans for cms luminosity calibration. EPJ Web Conf. 201, 04001 (2019).
https://doi.org/10.1051/epjconf/201920104001
")\] at IP1 and/or IP5. In this approach, the convolved transverse bunch sizes measured using beam-separation scans are translated into emittances using the \\(\\beta ^\*\\) values determined by *k*\-modulation \[[23](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR23 "F. Carlier, R. TomĂĄs, Accuracy and feasibility of the
$${\beta }^{*}$$
ÎČ
â
measurement for lhc and high luminosity lhc using
$$k$$
k
modulation. Phys. Rev. Accel. Beams 20, 011005 (2017).
https://doi.org/10.1103/PhysRevAccelBeams.20.011005
")\] at the relevant IP (see Sect. [4\.3](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec7)). The associated uncertainty (Table [3](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Tab3)) is combined with the statistical uncertainty in convolved transverse width to estimate the error affecting the measured emittance. Comparison with beam-averaged emittances extracted from the single-beam profile monitors reveals excellent agreement between WS and emittance-scan results (Fig. [2](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig2)). The relative time-evolution of the BSRT emittances is consistent with that observed using the WS or emittance scans, but the absolute magnitudes differ significantly, especially in the horizontal plane.
**Table 3 \\(\\beta ^\*\\) functions at the ATLAS and CMS IPs measured by *k*\-modulation. The target value in all cases is \\(\\beta ^\*=\\) 11 m**
[Full size table](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/tables/3)
**Fig. 2**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/2)
Time evolution of the horizontal (top) and vertical (bottom) beam-averaged normalized emittances during LHC fill 8037, as reported by the BSRT (blue), the WS (orange), and using beam-separation scans at IP1 (red) and IP5 (purple). Different shades of the same color correspond to two different bunches present in this fill. For ATLAS emittance scans, only the bunch-averaged value is shown
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/2)
### 4\.4 Beam-beam parameter determination
The systematic and quantitative comparison of the measured and of the predicted beam-beam impact on the observables detailed in Sect. [2](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec2) requires continuous monitoring of the actual beam-beam parameter throughout the duration of the experiment. Only the BSRT provides uninterrupted emittance determination throughout the fill. In view of the scale biases apparent in Fig. [2](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig2), however, WS-based emittances are used, whenever possible, as input to the determination of the beam-beam parameter; BSRT emittances re-scaled to close-in-time, absolute WS measurements provide time-interpolated emittance values whenever WS data are unavailable.
The beam-beam parameter evolution during the first fill of the experiment is illustrated in Fig. [3](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig3). The error bands are dominated by the 3% systematic uncertainty in the measured optical functions. The strength of the beam-beam interaction drops by a factor of two over a couple of hours, from the combined effect of beam-intensity decay and of emittance growth.
**Fig. 3**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/3)
Time evolution of the beam- and plane-averaged beam-beam parameter \\(\\xi \\) inferred from the measured bunch charges and emittances, separately for the two colliding-bunch pairs present in the fill pattern. The color bands indicate the systematic uncertainty (see text)
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/3)
### 4\.5 Tune measurements
Coherent spectra can be monitored using either the LHC Transverse Damper (ADT) \[[24](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR24 "M. SöderĂ©n, J. Komppula, G. Kotzian, S. Rains, D. Valuch, ADT and Obsbox in LHC Run 2, plans for LS2, pp. 165â171 (2019)")\], or the Base-Band Tune (BBQ) \[[25](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR25 "A. Boccardi, M. Gasior, R. Jones, P. Karlsson, R. Steinhagen, First results from the lhc bbq tune and chromaticity systems. Technical report, CERN, Geneva (2009).
https://cds.cern.ch/record/1156349
")\] system; their comparison typically yields consistent results. Spectrograms with a tune resolution of 0.0001 are computed using bunch positions recorded at every turn, and averaged over the one-minute time bins mentioned in Sect. [4\.1](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec5). To mitigate the influence of noise on the tune measurement, a median filter is applied at the pre-processing stage with a self-defined local window size. Additionally, the 50 Hz noise lines present in the spectra are masked. The filtered data are fitted by the sum of a Gaussian function and a constant baseline, and the frequency at the peak of the Gaussian is interpreted as the measured mean tune. The systematic uncertainty on measured tune shifts \\(\\Delta Q\\) is estimated empirically to be \\(\\sigma \_{\\Delta Q}=0.001\\).
### 4\.6 IP-orbit monitoring
Orbit displacements at each IP are measured using the Diode ORbit and OScillation (DOROS) BPM system, that provides sub-micrometer beam-position resolution \[[26](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR26 "M. Gasior, G. Baud, J. Olexa, G. Valentino, First Operational Experience with the LHC Diode ORbit and OScillation (DOROS) System, p. 07 (2017)
https://doi.org/10.18429/JACoW-IBIC2016-MOPG07
")\]. These strip-line BPMs are located in the two quadrupoles on either side of and closest to the IP, allowing the position of both beams to be measured simultaneously. The position and the angle of each beam at the IP are inferred from, respectively, the average and the difference of the positions measured in the two final-triplet quadrupoles.
## 5 Optimization of the phase advance between interaction points
In the two-IP configuration, the basic theory of linearized beam-beam \\(\\beta \\)\-beating provides an analytical determination of the optimum phase advance between the two IPs, *i.e.* of the setting that maximizes the sensitivity, at the witness IP, to beam-beam effects at the scanning IP (Sect. [5\.1](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec12)). A generalization to the three- and four-IP configurations is discussed in Sect. [5\.2](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec13). The required modifications to the baseline LHC optics, as well as their implementation and their experimental validation, are detailed in Sect. [5\.3](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec14).
### 5\.1 Linear beam-beam \\(\\beta \\)\-beating
In circular colliders, the \\(\\beta \\)\-beating induced by beam-beam effects at the IP(s) leads to changes in transverse beam size \[[5](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR5 "W. Herr, T. Pieloni, Beam-Beam Effects. Technical report, CERN (2014). Contribution to the CASâCERN Accelerator School: Advanced Accelerator Physics Course, Trondheim, Norway, 18â29 Aug 2013, p. 29.
https://doi.org/10.5170/CERN-2014-009.431
.
https://cds.cern.ch/record/1982430
"), [27](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR27 "A.W. Chao, Coherent beam-beam effects, in Frontiers of Particle Beams: Intensity Limitations (Springer, Berlin, Heidelberg, 1992), pp. 363â414.
https://doi.org/10.1007/3-540-55250-2_36
")\], thereby impacting the luminosity. For a single collision point, the change in \\(\\beta ^\*\\) of a single small-amplitude particle in (for instance) the horizontal plane *x*, is given by:
\$\$\\begin{aligned} \\frac{\\beta \_x^\*}{\\beta \_{0,x}^\*} = \\frac{1}{\\sqrt{1-4\\pi \\xi \\cot {(2\\pi Q\_x)} - 4\\pi ^2\\xi ^2}}. \\end{aligned}\$\$
(6)
The magnitude of the effect depends on the absolute value \\(\\xi \\) of the beam-beam parameter and on the nominal betatron tune \\(Q\_x\\); its periodicity is \\( \[\\pi \]\\). It can be shown that approaching the half-integer betatron tune from below minimizes the \\(\\beta \\) function at the IP, thereby maximizing the luminosity; this phenomenon has been exploited with great success in (among others) the KEKB and PEP-II *B* factories. In the LHC, where during collisions the nominal fractional betatron tunes are fixed at \\(q\_x=0.31,\\,q\_y=0.32\\), this dynamic-\\(\\beta \\) effect depends only on the beam-beam parameter.
This conclusion no longer strictly holds in the presence of more than one collision point. In the case of two IPs with identical values of \\(\\xi \\) and \\(\\beta ^\*\\), the dynamic-\\(\\beta \\) effect can be described analytically \[[28](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR28 "J. WaĆczyk, A. Dabrowski, D. Stickland, T. Pieloni, W. Kozanecki, X. Buffat, Impact of multiple beam-beam encounters on LHC absolute-luminosity calibrations by the van der Meer method. JACoW IPAC 2023, 053 (2023).
https://doi.org/10.18429/JACoW-IPAC2023-WEPA053
")\]:
\$\$\\begin{aligned} \\frac{\\beta \_x^\*}{\\beta ^\*\_{0,x}} = \\frac{\\sin {2\\pi Q\_x} + 4\\pi \\xi (\\cos (2\\pi Q\_x-2\\mu \_{1,x})-\\cos 2\\pi Q\_x )}{\\pm \\sqrt{1-\\cos ^2{2\\pi (Q\_x+\\Delta Q\_x)}}}, \\end{aligned}\$\$
(7)
where \\(\\mu \_{1,x}\\) is the phase advance from the scanning IP to the witness IP. The sign in the denominator is linked to that of the numerator so as to ensure that the ratio remains positive; the denominator must be positive (resp. negative) below (resp. above) the half-integer tune, *i.e.* when \\(m\< Q\_x \<m+\\frac{1}{2}\\) (resp. \\(m+\\frac{1}{2}\< Q\_x \< m+1\\)), where *m* is a positive integer. The tune shift \\(\\Delta Q\_x\\) is related to \\(Q\_x\\), \\(\\xi \\) and \\(\\mu \_{1,x}\\) by:
\$\$\\begin{aligned} \\cos {2\\pi (Q\_x+\\Delta Q\_x)}&=(1-16\\pi ^2\\xi ^2)\\cos {2\\pi Q\_x} + 8\\pi \\xi \\sin {2\\pi Q\_x} \\nonumber \\\\&\\quad + 16\\pi ^2\\xi ^2\\cos {(2\\pi Q\_x-2\\mu \_{1,x})}. \\end{aligned}\$\$
(8)
The phase-advance dependence of the resulting \\(\\beta \\)\-beating is illustrated in Fig. [4](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig4). The minimum of each curve corresponds to the phase advance that yields the largest mutual beamâbeam-induced luminosity enhancement between the two IPs. The optimal settings are \\(\\mu \_{1,x}^{min}/2\\pi =0.405,\\,\\mu \_{1,y}^{min}/2\\pi =0.410\\); their difference reflects that between the nominal horizontal and vertical fractional tunes \\(q\_x\\) and \\(q\_y\\).
**Fig. 4**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/4)
Analytically computed beam-beam induced \\(\\beta \\)\-beating in a two-IP configuration, as a function of the phase advance \\( \\mu \_1\\) between the 2 IPs (\\(q\_x=0.31\\), \\(q\_y=0.32\\), head-on collisions at both IPs, \\(\\xi =7\\times 10^{-3}\\) per IP)
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/4)
The extrema of the curves in Fig. [4](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig4) can be determined analytically by differentiating Eq. ([7](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Equ7)) with respect to the phase advance \\(\\mu \_{1,x}\\). The general expression for this derivative \\(\\frac{d}{d\\mu \_{1,x}} \\biggl ( \\frac{\\beta \_x^\*}{\\beta ^\*\_{0,x}} \\biggr )\\), detailed in Appendix A, is the product of two factors, each of which can be zero:
- imposing that \\(\\sin {(2\\pi Q\_x - 2\\mu \_{1,x})} = 0\\) yields the solution
\$\$\\begin{aligned} \\frac{\\mu \_{1,x}}{2\\pi } =\\frac{Q\_x}{2}-\\frac{m}{4} \\, \\end{aligned}\$\$
(9)
where *m* is an integer of either sign, or zero;
- the other factor is zero for either
\$\$\\begin{aligned} \\frac{\\mu \_{1,x}}{2\\pi } = \\frac{Q\_x}{2} - A(\\xi , Q\_x)\\, \\end{aligned}\$\$
(10)
or
\$\$\\begin{aligned} \\frac{\\mu \_{1,x}}{2\\pi } = \\frac{Q\_x - 1}{2} + A(\\xi , Q\_x) \\, \\end{aligned}\$\$
(11)
with \\(A(\\xi , Q\_x)\\) defined in Eq. ([12](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Equ12))
\$\$\\begin{aligned} A(\\xi , Q\_x) = \\arccos \\left( \\frac{ 32\\pi ^3\\xi ^3 \\sin (4\\pi Q\_x) + 48\\pi ^2\\xi ^2 \\cos ^2(2\\pi Q\_x) - 32\\pi ^2\\xi ^2 - 6\\pi \\xi \\sin (4\\pi Q\_x) + \\cos ^2(2\\pi Q\_x) - 1}{16\\pi ^2\\xi ^2 \\left( 4\\pi \\xi \\sin (2\\pi Q\_x) + \\cos (2\\pi Q\_x) \\right) } \\right) / (4\\pi ) \\end{aligned}\$\$
(12)
The two sets of solutions do not overlap. The subset of solutions that minimize (rather than maximize) \\(\\beta ^\*\\) is identified by requiring that the second derivative be positive. This restricts the main solution to:
\$\$\\begin{aligned} \\frac{\\mu \_{1,x}}{2\\pi } =\\frac{Q\_x}{2}+\\frac{m+1}{4} . \\end{aligned}\$\$
(13)
The optimal phase \\(\\mu \_{x,1}\\) that minimizes the \\(\\beta ^\*\\) for range of \\((\\xi , Q\_x)\\) values is shown in Fig. [5](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig5).
**Fig. 5**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/5)
Optimal phase advance value \\(\\mu \_{x,1}\\) for various values of \\((\\xi , Q\_x)\\) in the two-IPs configuration
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/5)
The phase-advance dependence of the dynamic-\\(\\beta \\) effect can also be quantified in terms of the luminosity shift, at the witness IP, that is associated with the electromagnetic interaction of the two bunches at the scanning IP (Fig. [6](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig6)). The optimal setting, *i.e.* that which maximizes the sensitivity, at the witness IP, to beam-beam effects at the scanning IP, is indicated by the green dotted vertical line at \\(\\Delta \\mu \_{\\text {IP1-IP5}}=0.41\\ \[2\\pi \]\\); the full suppression of the phase-related luminosity enhancement is indicated by the red dotted line. These results are consistent with the predictions of the analytical model (Fig. [4](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig4)). In the two-IP configuration considered here, in which IP1 and IP5 can both play the role of either the scanning or the witness IP, and due to the periodicity of this curve (\\( \\pi \\)) and to the fractional-tune values used during collisions (\\(q\_x / q\_y = 0.31 / 0.32)\\), the phaseâadvance-related luminosity enhancement is the same at the two IPs; at the optimum setting, it is three times larger than it would be using the nominal LHC lattice. The maximum effect depends mainly on the transverse tunes, and thus could have been further enhanced if it had been possible to move these away from their nominal values.
A more comprehensive simulation study, that evaluates the impact of the beam-beam interaction on the luminosity for more realistic Gaussian proton-density distributions, multiple collision points and as a function of the beam separation, is reported in Ref. \[[29](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR29 "J. WaĆczyk, Precision luminosity measurement at hadron colliders. PhD Thesis, Ăcole Polytechnique FĂ©dĂ©rale de Lausanne (2024).
https://doi.org/10.5075/epfl-thesis-10500
")\].
**Fig. 6**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/6)
Beamâbeam-induced luminosity shift at IP5 predicted by COMBI simulations, as a function of the phase advance \\( \\mu \_1=\\Delta \\mu \_{\\text {IP1-IP5}}\\) between the two IPs (\\(q\_x=0.31\\), \\(q\_y=0.32\\), head-on collisions at both IPs, \\(\\xi =7\\times 10^{-3}\\) per IP)
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/6)
### 5\.2 Optimized phases for multi-IP configurations
The general expression to obtain the first-order change of the \\(\\beta \\) function at an arbitrarily chosen reference IP (at \\(s=\\mu \_{0,x}=0\\)) can be approximated as multiple quadrupole errors \[[30](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR30 "T. Pieloni, J. Barranco, X. Buffat, P.C. Jorge, L.M. Medrano, C. Tambasco, R. TomĂĄs, Dynamic beta and beta-beating effects in the presence of the beam-beam interactions in 57th ICFA Advanced Beam Dynamics Workshop on High-Intensity and High-Brightness Hadron Beams, p. 027 (2016).
https://doi.org/10.18429/JACoW-HB2016-MOPR027
")\]:
\$\$\\begin{aligned} \\frac{\\beta ^\*\_x}{\\beta ^\*\_{0,x}} = \\frac{2\\pi \\xi }{\\sin {2\\pi Q\_x}} \\sum \_{i\\in IPs}^{N-1} \\cos (2\\pi Q\_x-2\\mu \_{i,x})\\, \\end{aligned}\$\$
(14)
where \\(\\mu \_{i,x}\\) is the horizontal phase advance between the reference and the \\(i^{th}\\) IP for a total of *N* IPs. Therefore, the optimal configuration, which minimizes the \\(\\beta ^\*\\) change at all IPs, must satisfy the condition:
\$\$\\begin{aligned} \\sum \_{i=1}^{N-1} \\cos (2\\pi Q\_x-2\\mu \_{i,x}) = -(N-1)\\, \\end{aligned}\$\$
(15)
which yields the same solution for all IPs as in the 2 IPs configuration (Eq. [13](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Equ13)):
\$\$\\begin{aligned} \\frac{\\mu \_{i,x}}{2\\pi } =\\frac{Q\_x}{2}+\\frac{n+1}{4} . \\end{aligned}\$\$
(16)
where *n* is an integer of either sign, or zero. To maximize the effect on luminosity, the phase advance between the reference interaction point at \\( i = 0 \\) and the nearest IP, at \\( i = 1 \\) must satisfy \\( \\mu \_{1,x}/2\\pi = (Q\_x + 0.5)/2 \\), as given by the first solution (for \\(n=0\\)) of Eq. ([13](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Equ13)). This condition implies that the phase advance to any subsequent adjacent IPs (\\( i \\ge 2 \\)) must be an integer multiple of \\( \\pi \\). Thus, this requirement exposes a fundamental constraint: it is not possible to minimize \\( \\beta ^\* \\) simultaneously at all IPs. This would only be possible for (\\( Q\_x = 0.5 \\)) which is practically not achievable, as it would require operating the storage ring precisely at the half-integer resonance.
### 5\.3 Optics validation
The optimum phase advance from IP1 to IP5 determined in Sect. [5\.1](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec12) was implemented using tuning trim quadrupoles in the LHC arcs. The targeted fractional phase advance was \\(0.9 \\, \[2 \\pi \]\\) for both beams and both planes (Table [4](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Tab4)). This setting, which differs from that in Fig. [6](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig6) by one period (\\( \\pi \\)), is equivalent in terms of sensitivity optimization but has the advantage that it minimizes the absolute magnitude of the perturbations inflicted upon the nominal lattice.
Several measurements were carried out to validate this adjustment: comparison of the targeted and achieved shifts in phase advance (Sect. [5\.3.1](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec15)), \\(\\beta \\)\-function measurements across the full LHC circumference (Sect. [5\.3.2](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec16)) as well as at the IPs (Sect. [5\.3.3](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec17)), and dispersion measurements to estimate the dispersive contribution to the transverse beam size at the IP (Sect. [5\.3.4](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec18)).
#### 5\.3.1 Phase advance
The phase advance around each LHC ring is determined using turn-by-turn orbit measurements \[[20](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR20 "R. TomĂĄs, O. BrĂŒning, M. Giovannozzi, P. Hagen, M. Lamont, F. Schmidt, G. Vanbavinckhove, M. Aiba, R. Calaga, R. Miyamoto, Cern large hadron collider optics model, measurements, and corrections. Phys. Rev. ST Accel. Beams 13, 121004 (2010).
https://doi.org/10.1103/PhysRevSTAB.13.121004
"), [21](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR21 "R. TomĂĄs, M. Aiba, A. Franchi, U. Iriso, Review of linear optics measurement and correction for charged particle accelerators. Phys. Rev. Accel. Beams 20, 054801 (2017).
https://doi.org/10.1103/PhysRevAccelBeams.20.054801
")\]. The targeted and achieved phase-advance shifts agree within \\(0.01\\, \[2\\pi \]\\) or better (Table [4](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Tab4)). The phase advance change from IP1 and IP5 was compensated on the opposite side of the ring, ensuring that the transverse tunes remain unchanged. The target fractional local phase advance is set to \\(0.9,\[2\\pi \]\\). Table [4](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Tab4) lists the corresponding phase shifts applied to the nominal lattice to reach this target on one side of the ring; this is consistent with the analytically derived value of \\(0.41\\,\[2\\pi \]\\) for the compensating shift on the opposite side. The change in phase advance was also measured at the BSRT location, where the discrepancy between the model and measurements is more pronounced, reaching up to \\(0.05\\,\[2\\pi \]\\).
**Table 4 Nominal and optimized IP1 to IP5 phase-advance settings. Column 4 (resp. 5) displays the targeted (resp. measured) shift in phase advance after the optimized settings have been applied**
[Full size table](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/tables/4)
#### 5\.3.2 Residual \\(\\beta \\)\-beating
The \\(\\beta \\)\-beating was quantified and found to be well within acceptable tolerances for various configurations of the trim strength. The measurement results, relative to the reference measurements without any phase advance adjustments, are illustrated in Fig. [7](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig7). The maximal observed deviation is approximately 25% around IP5 for Beam 1 in the vertical plane. Such significant deviations in proximity to the IPs are anticipated and are frequently attributed to erroneous BPM readings. Overall, the measurements indicated that \\(\\beta \\)\-beating remained within acceptable bounds, predominantly within a 10% margin.
**Fig. 7**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/7)
Measured \\(\\beta \\)\-function change from the full strength of the phase advance knob with reference to nominal lattice along the LHC ring, for Beam 1 (top two figures for *x* and *y*) and Beam 2 (bottom two figures for *x* and *y*)
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/7)
The measured \\(\\beta \\) functions are also compared to the MADX LHC lattice model predictions, which include the phase advance adjustment. The results, illustrated along the LHC lattice in Fig. [8](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig8), indicate that the measurements align with the model within an acceptable margin of 10%.
**Fig. 8**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/8)
Measured deviations in the \\(\\beta \\)\-function along the LHC ring relative to the MADX model incorporating the phase advance knob, for Beam 1 (top two figures for *x* and *y*) and Beam 2 (bottom two figures for *x* and *y*). The error bars represent solely the statistical uncertainties as determined in the measurements
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/8)
#### 5\.3.3 \\(\\beta \\) functions at the interaction points
The DOROS BPMs are strategically located at the end of the inner triplets on both sides of each of the experiments. Due to the phase advance of \\(\\pi \\) between the DOROS on the left and right sides, these are inadequate for the precise determination of the measurement of \\(\\beta ^\*\\) at the IP. This can only be achieved with a special *k*\-modulation technique \[[31](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR31 "F. Carlier, et al. Challenges of K-modulation measurements in the LHC Run 3, in Proceedings of IPACâ23. IPACâ23â14th International Particle Accelerator Conference (JACoW Publishing, Geneva, Venice, 2023), pp. 531â534.
https://doi.org/10.18429/JACoW-IPAC2023-MOPL014
")\], which was performed at the interaction points used in the experiment. The results are summarized in Table [3](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Tab3). While the achieved precision of the \\(\\beta ^\*\\) determination varies between measurements, uncertainties at the level of \\({\\mathcal {O}}(1\\%)\\) are not uncommon and are not systematic to a specific beam or interaction point. Such uncertainties would propagate to only a few-percent effect on derived quantities such as the emittance or the beamâbeam parameter.
#### 5\.3.4 Dispersion
It was essential to perform dispersion \\(D^\*\\) measurements at the interaction points, under the new optics configurations, as it can contribute to the measured beam size. Assuming the LHC design momentum spread value of \\(\\frac{\\Delta p}{p}=10^{-4}\\), the average dispersion measurement around IP1 and IP5 is estimated to contribute up to 2% to the observed bunch sizes (worst case measured for B2Y). This effect is considered static throughout the experiment.
## 6 Experimental results
This section presents experimental measurements of beamâbeam effects in configurations involving multiple collision points. The impact of the optics modifications introduced in Sect. [5](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec11) is first assessed through tune spectra measurements, as detailed in Sect. [6\.1](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec20). To quantify the cumulative effect of one or more additional collisions and assess the reproducibility of the results, a stepwise collision scan was employed. In this procedure, beams were sequentially brought into head-on collision at selected interaction points (IPs), followed by full transverse separation (\\(\\Delta \> 6\\sigma \\)) in either the horizontal or vertical plane at all but one witness IP. Two distinct filling schemes were used to isolate the influence of a single additional collision (2-IP configuration: ATLAS and CMS, see Sect. [6\.2](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec21)) and to investigate more complex multi-collision scenarios (3-IP configuration: including ALICE, see Sect. [6\.3](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec25)). The first measurement of separation-dependent beamâbeam effects, analogous to those encountered during luminosity calibration scans, is presented in Sect. [6\.4](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec28), alongside comparisons with numerical predictions. An additional measurement conducted during the 2022 van der Meer calibration is detailed in Sect. [6\.5](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec32). A consolidated overview of all findings is provided in Sect. [7](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec33).
### 6\.1 Measured effect on the tune spectra of the phase advance changes
The coherent transverse tune spectra serve as a diagnostic for evaluating the effectiveness of the phase advance modifications. COMBI simulations were carried out for two extreme lattice configurations, one designed to maximize and the other to suppress beamâbeam effects. In the head-on scenario (Fig. [9](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig9)), the coherent modes are suppressed for the maximizing configuration, as the phase advance between the two collision points is close to \\(\\pi \\), as was previously shown in simulation in Ref. \[[32](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR32 "Y. Alexahin, A study of the coherent beam-beam effect in the framework of the Vlasov perturbation theory. Nucl. Instrum. Methods Phys. Res., Sect. A 480(1â3), 253â288 (2002).
https://doi.org/10.1016/S0168-9002(01)01219-0
")\]. In this configuration, the effect of the first collision is counteracted by the second. This setting is particularly suitable for studying incoherent spectra, as the diminished coherent signal allows for a clearer extraction of the mean tune shift from fitted spectra. During separation scans, the symmetry between horizontal (*x*) and vertical (*y*) planes is broken. In the scanning plane (horizontal), coherent modes re-emerge at \\(1.5\\,\\sigma \\) separation (Fig. [10](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig10)), with their frequency shifting in response to the separation distance (additional example for \\(3\\,\\sigma \\) is shown in Fig. [11](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig11)).
**Fig. 9**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/9)
Comparison of COMBI simulated tune spectra in *x*\-plane (left) and *y*\-plane (right) for the suppressing (in blue) and maximizing (in green) phase configurations, for head-on collision
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/9)
**Fig. 10**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/10)
Comparison of COMBI simulated tune spectra in *x*\-plane (left) and *y*\-plane (right) for the suppressing (in blue) and maximizing (in green) phase configurations, at \\(1.5\\,\\sigma \\) separation in the horizontal plane
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/10)
**Fig. 11**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/11)
Comparison of COMBI simulated tune spectra in *x*\-plane (left) and *y*\-plane (right) for the suppressing (in blue) and maximizing (in green) phase configurations, at \\(3\\,\\sigma \\) separation in the horizontal plane
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/11)
These simulated trends of the coherent spectra are also observed in experimental data. Transverse spectra recorded for head-on conditions confirm the suppression of coherent modes in the maximizing phase configuration (Fig. [12](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig12)). Additional comparisons for separated beams at the IP (see Figs. [13](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig13) and [14](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig14)), shown here for the case of vertical separation, reveal similar features. In this configuration, the non-scanning (horizontal) plane exhibits spectra that qualitatively agree with simulation predictions (see Figs. [10](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig10) and [11](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig11)), while in the scanning (vertical) plane the coherent modes do not reappear as expected. An equivalent behavior is observed when the separation is applied in the horizontal plane, with the roles of the two transverse planes interchanged; this case is not shown for brevity. This discrepancy is attributed to the effect of the transverse damper, which acts independently on the two beam centroids and suppresses coherent oscillations arising from beamâbeam coupling, consistent with earlier observations \[[33](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR33 "R. Calaga, W. Herr, G. Papotti, T. Pieloni, X. Buffat, S. White, R. Giachino, Coherent beam-beam mode in the LHC. CERN (2014).
https://doi.org/10.5170/CERN-2014-004.227
.
http://cds.cern.ch/record/1957039
")\].
**Fig. 12**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/12)
Measured tune spectra before and after phase change, corresponding to the suppressing (in blue) and maximizing (in green) phase configurations, at the head-on position, for each of the transverse planes (Beam 1 shown as an example). Dash-dotted lines show Gaussian fits to the data
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/12)
**Fig. 13**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/13)
Comparison of the measured tune spectra at suppressing (blue) and maximizing (green) phase configurations with \\(1.3\\,\\sigma \\) separation in the vertical plane
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/13)
**Fig. 14**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/14)
Comparison of the measured tune spectra at suppressing (blue) and maximizing (green) phase configurations with \\(2.6\\,\\sigma \\) separation in the vertical plane
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/14)
### 6\.2 Head-on collisions in a two-IP configuration at \\(\\sqrt{s} = 900\\ \\text {GeV}\\)
Each test consisted of four measurements designed to quantify the beamâbeam effect induced by a single additional head-on collision. The sequence was executed with CMS as the witness IP, followed by ATLAS, and concluded with a repeat at CMS. Over the course of the fill, the beamâbeam parameter decreased from approximately 0.01 to 0.005. The outcome of each scan is analyzed using three independent observables: luminosity, beam size, and tune spectra, discussed separately below.
#### 6\.2.1 Impact of beam-beam effects on head-on luminosity
An overview of the luminosity measurements during the first test with CMS as the witness IP is shown in Fig. [15](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig15), where the specific luminosity at CMS is normalized to the initial head-on step to highlight relative changes. The uncertainty at each step is given by the standard error of the mean; shaded bands represent the standard deviation. A comparable dataset for ATLAS is presented in Fig. [16](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig16), confirming the expected luminosity change due to beamâbeam effects originating from another IP. Relative luminosity changes for each step in all tests are summarized in Table [5](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Tab5). These values are computed relative to the average of adjacent head-on points, which removes the linear luminosity decay over time. Assuming slow change of beam parameters, the results within each test are expected to be consistent.
**Fig. 15**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/15)
CMS specific luminosity, normalized to the first measured both IPs head-on point (blue), during the first step-function scan at IP1. The measurements with ATLAS fully separated are shown in red. The error bars indicate error on the mean
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/15)
**Fig. 16**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/16)
ATLAS specific luminosity, normalized to the first measured both IPs head-on point (blue), during the step-function scan at IP5. The points for which IP5 was separated are shown in red. The error bars indicate the statistical error
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/16)
**Table 5 Measured luminosity shifts in each of full separation steps â summary of all three step-function scans performed for IP1 and IP5 as the witness IP. The beam-beam parameter at the start of each scan \\(\\xi \_\\text {start}\\) is included for reference**
[Full size table](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/tables/5)
All measurements are compiled in Fig. [17](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig17) and compared to COMBI simulation predictions. The results are plotted as a function of the beamâbeam parameter (estimated as described in Sect. [4\.4](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec8)), following a reverse chronological order. Agreement with simulation is excellent across the entire range.
**Fig. 17**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/17)
COMBI simulated luminosity enhancement induced by a single head-on collision as a function of the beam-beam parameter (dashed line), compared to the test results at both ATLAS (red points) and CMS (green points)
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/17)
**Fig. 18**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/18)
Measurements of beam-beam effects on beam widths shown with red points for each beam and plane, during the experiment with IP5 as the witness IP at \\(\\xi =0.0056\\). COMBI predictions are shown in blue with calculated systematic error from the phase advance error between the IP and BSRT locations
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/18)
#### 6\.2.2 Impact of beam-beam effects on single-beam sizes
Complementary measurements of the bunch size from the BSRT are shown in Fig. [18](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig18). Data from all steps and bunches are combined, and a linear trend corresponding to continuous emittance growth is subtracted using the difference between the first and last head-on steps. For comparison, COMBI simulation results are also included in the figure. The dominant source of uncertainty in COMBI predictions arises from discrepancies between the measured phase advance and the MADX model, with the largest observed deviation being \\(0.05\\,\[2\\pi \]\\) (as discussed in Sect. [5\.3.1](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec15)). This uncertainty is propagated using a modified form of Eq. ([6](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Equ6)) adapted to the BSRT location. The magnitude of this uncertainty depends on the specific phase advance configuration, as some are more sensitive to small deviations. This is reflected in the varying sizes of the error bands associated with each prediction. The most accurate predictions are obtained in the vertical plane of Beam 2, where the beamâbeamâinduced \\(\\beta \\)\-beating propagates to the BSRT with a particularly favorable phase advance. In this configuration, the measured effect is also the most pronounced, and excellent agreement with COMBI simulations is observed. In other cases, the expected absolute effect is smaller, typically below 1%, and the influence of phase advance uncertainty becomes more significant, limiting the predictive reliability of the simulations. Experimentally, it is also possible that in these configurations the beamâbeamâinduced effect is too subtle to resolve, potentially masked by instrumental limitations.
#### 6\.2.3 Impact of beam-beam effects on tune spectra
Figure [19](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig19) shows representative ADT spectra recorded during the step-function scans. These spectra are affected by coupling, noise, and various nonlinearities, so they do not follow an idealized analytical distribution. Nonetheless, the peak position can be extracted using a Gaussian fit (indicated by dashed lines), yielding an estimate of the mean tune shift. In the single-collision configuration (blue lines), unsuppressed coherent modes dominate the spectrum and partially obscure the underlying incoherent component. The Gaussian fit is therefore applied to the incoherent contribution between coherent lines, rather than to the global spectral maximum, which results in a larger uncertainty on the extracted peak position. Despite this, in each of the performed step-scans at both collision points, the single-collision tune shift, defined as twice the fitted mean shift, is found to be consistent within \\(\\pm 5\\%\\) with the expected coherent beamâbeam tune shift corresponding to the conditions of the respective scan. For head-on collisions with Gaussian transverse beam distributions, the coherent tune shift of the dipole mode is expected to be \\(Y\\xi \\), where *Y* is the Yokoya factor (\\(Y = 1.21\\) for round beams) \[[34](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR34 "K. Yokoya, Y. Funakoshi, E. Kikutani, H. Koiso, J. Urakawa, Tune shift of coherent beam-beam oscillations. Part. Accel. 27, 181â186 (1990)")\]. As expected, the two-collision configuration yields a tune shift approximately double that of the single-collision case \[[13](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR13 "A. Babaev, T. Barklow, O. Karacheban, W. Kozanecki, I. Kralik, A. Mehta, G. Pasztor, T. Pieloni, D. Stickland, C. Tambasco, R. Tomas, J. WaĆczyk, Impact of beam-beam effects on absolute luminosity calibrations at the CERN large hadron collider. Eur. Phys. J. C 84, 17 (2024).
https://doi.org/10.48550/arXiv.2306.10394
")\].
**Fig. 19**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/19)
Measured vertical tune spectra for different collision configurations: reference nominal tune (red), without the presence of any collisions, compared to single collision (blue) and two collisions (green) at \\(\\xi \\simeq 0.0086\\). The dash-dotted lines in corresponding colors indicate Gaussian fits to the data used to extract mean tunes
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/19)
**Fig. 20**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/20)
Specific luminosity, normalized to the first head-on point, measured at ATLAS (top) and CMS (bottom) during beam separation at the other IPs, with horizontal separation applied in steps 2â4 and vertical separation in steps 6â8. The tested configurations are indicated in the legend. The blue line shows an interpolation of the head-on measurements, serving as a reference for relative luminosity changes
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/20)
### 6\.3 Head-on collisions in a three-IP configuration at \\(\\sqrt{s} = 900\\ \\text {GeV}\\)
An additional head-on collision was included in the configuration using IP2 (ALICE), as it was the only IP where this was feasible. A specially designed filling scheme ensured uniform collision exposure across all bunches, allowing measurements to be combined and thereby reducing statistical uncertainties. However, no reliable luminosity data were available from IP2 during this experiment.
**Table 6 Luminosity shift per full separation step, calculated with respect to an interpolated reference obtained from the two closest head-on points. Results for horizontal and vertical separations are shown as Shift part 1 and Shift part 2, respectively**
[Full size table](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/tables/6)
#### 6\.3.1 Impact of beam-beam effects on head-on luminosity
A step-function scan was conducted by sequentially collapsing or separating one IP at a time, with the separation applied first in the horizontal plane and, in a second stage of the experiment, in the vertical plane. A comprehensive summary of the measured specific luminosity at ATLAS and CMS is presented in Fig. [20](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig20). Exceptionally stable beam conditions during this test enabled measurements at both experiments under nearly identical beamâbeam parameters. The beamâbeam-induced luminosity changes observed during each step are quantified in Table [6](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Tab6), referenced to the time-interpolated specific luminosity from three head-on points.
Orbit-drift corrections are applied only to these head-on reference points and result in changes of the extracted luminosity shifts of at most 0.2%. Good agreement with COMBI simulations is observed across the separation steps, with measured shifts closely matching the simulated values for configurations involving IP1 and IP5. For single-IP separation, the level of agreement depends on the witness IP and the separation step. When IP1 is used as the witness, the luminosity shift in the first separation step is reproduced almost exactly by COMBI, while a smaller shift is observed in the second step. This behaviour is expected, as the beamâbeam parameter decreases during the scan and the measurement becomes increasingly sensitive to residual orbit perturbations. For configurations involving IP2, the comparison is less conclusive. As no reliable luminosity data were available from IP2 during this experiment, it was not possible to precisely optimize the collision conditions at this IP. As a result, the actual beamâbeam conditions at IP2 may differ from those assumed in the simulations, affecting the comparison. Residual orbit drifts affecting the separation steps cannot be fully corrected. As illustrated by the DOROS measurements in Fig. [21](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig21), the orbit data contain both slow drifts and beamâbeamâinduced orbit offsets arising from changes in the collision configuration when other IPs are separated. While the beamâbeamâinduced offsets are included in the COMBI simulations, the residual orbit drifts are not. Since these contributions cannot be disentangled in the measurements, a full correction is not possible. Consequently, the effective beamâbeam impact measured during the experiment may be weaker than that assumed in the simulations.
**Fig. 21**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/21)
Residual separation between the two beams at all three IPs in horizontal (top) and vertical (bottom) planes, for test with observations at ATLAS, where the separation was performed first in the horizontal and then in the vertical plane at the other IPs
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/21)
#### 6\.3.2 Tune shift scaling with the number of beamâbeam interactions
Analysis of the transverse tune spectra in this configuration is further complicated by the reappearance of coherent modes when beams are separated at one or more IPs, as predicted by simulation. In such cases, the maximum tune shift is determined from the frequency difference between the \\(\\pi \\) and \\(\\sigma \\) modes, rather than from fitting the incoherent peak. For configurations where all three IPs were in collision, coherent modes remained suppressed due to the retained symmetry between IP1 and IP5. Figure [22](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig22) presents the measured tune shifts for all tested collision configurations. The unperturbed tune measurement, shown at â0â collisions, serves as the reference from which all tune shifts are calculated. An empirical uncertainty of 0.001 is assumed for each independent measurement. Both transverse planes and separation directions are included to demonstrate consistency. As expected, the measured tune shift scales linearly with the number of head-on collisions.
**Fig. 22**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/22)
Measured mean tune shift as a function of the number of head-on collisions, shown for both transverse planes and for both separation directions (horizontal and vertical) used in the experiment. The separation refers to how collisions were removed: it was applied in only one transverse direction at a time (H sep. or V sep.) at the non-colliding interaction points, which are not explicitly labeled on the x-axis
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/22)
### 6\.4 Beam-beam signatures during transverse-separation scans at \\(\\sqrt{s} = 900\\ \\text {GeV}\\)
In the van der Meer (vdM) calibration method, beamâbeam corrections are applied separately at each separation step. Therefore, it is essential to validate the simulated dependence of luminosity shifts on nominal beam separation. Full vdM-like scans were performed in both transverse planes at ATLAS and CMS, using a configuration where only these two IPs were in collision.
#### 6\.4.1 Impact of beamâbeam effects on the luminosity scan-curves
Figure [23](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig23) shows the luminosity curves measured at ATLAS, acting as the witness IP, during a separation scan performed at CMS. The curves are fitted with Gaussian profiles for visual guidance. The beam separation at CMS is expressed in units of the measured transverse beam size, \\(\\sigma \_{\\textrm{meas}}\\), which was found to be smaller than the nominal beam size, \\(\\sigma \_{\\textrm{nom}}\\), used to define the scan steps. As a result, scan points corresponding to nominal separations beyond \\(6\\,\\sigma \_{\\textrm{nom}}\\) extend to larger values when expressed in units of \\(\\sigma \_{\\textrm{meas}}\\).
To account for emittance evolution during the scan, both linear and exponential fits are applied to the head-on luminosity points before, after, and in the middle of the scan. These trends are subtracted from the data in Fig. [23](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig23), and all measurements are presented relative to this fitted head-on reference (normalized to unity).
**Fig. 23**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/23)
ATLAS specific luminosity relative changes (blue points) during separation scans in *x*\-plane (top) and *y*\-plane (bottom) at CMS, corrected for exponential (exp. corr.) luminosity decay. Gaussian function with a constant fitted to these points is shown with blue dashed line and compared to fit performed on the data corrected for linear luminosity change (lin. corr.) in time (dotted line). These are compared to COMBI simulation predictions (purple solid line)
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/23)
At CMS, similar scans were performed; however, degraded beam stability was observed during this period of the fill. The scan results exhibit systematic deviations, especially at large separation, which could not be symmetrically corrected using simple linear or exponential models. Despite these issues, all four scan cases (two transverse planes at two IPs) yield a total specific luminosity reduction of approximately 2.4%, consistent with expectations based on the measured beamâbeam parameter \\(\\xi \\). Comparison with COMBI simulations reveal up to 30% discrepancy for IP5 as the witness IP, while agreement improves to within 10% for measurements at IP1. These differences are comparable to those observed in Table [6](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Tab6), suggesting a common underlying origin. While a finer scan granularity can improve the description of the steep-response region below \\(3\\,\\sigma \_{\\textrm{meas}}\\), it would not address the dominant discrepancy between head-on and fully separated luminosity, which is likely driven by simulation assumptions and scan-related systematics.
#### 6\.4.2 Beamâbeam-induced orbit deflections
Beamâbeam deflections are clearly observed during separation scans. Figure [24](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig24) shows the fully corrected residual horizontal beam displacement at IP1, measured with the DOROS beam position monitors, as a function of the beam separation. This representation allows a direct comparison with the analytical beamâbeam deflection prediction based on the BassettiâErskine formalism \[[35](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR35 "M. Bassetti, G.A. Erskine, Closed expression for the electrical field of a two-dimensional gaussian charge. Technical report, CERN, Geneva (1980).
https://cds.cern.ch/record/122227
")\].
The raw orbit data were corrected for slow orbit drifts by aligning the beam positions at the luminosity optimization points before and after each scan. Deviations of the applied beam separation from the nominal settings were accounted for by fitting a common length scale simultaneously with the beamâbeam deflection. This approach avoids an overestimation of the separation calibration. The measured orbit displacement at the DOROS location was corrected for the contribution arising from the beamâbeam deflection angle by applying a geometric amplification factor, proportional to the distance between the IP and the BPM in the small-angle approximation. After all corrections, the measured beamâbeam deflection as a function of beam separation closely follows the analytical prediction. The maximum horizontal deflection observed at IP1 is 10.2 \\(\\upmu \\)m, compared to an expected value of 8.4 \\(\\upmu \\)m, with agreement within the total measurement uncertainty. The dominant uncertainty contribution (3.5 \\(\\upmu \\)m, corresponding to about 90% of the total) arises from a possible residual half crossing angle at the IP of up to 0.5 \\(\\upmu \\)rad. Additional systematic uncertainties include the fitted separation length scale, tune uncertainty of 0.005, and an extended 2 % uncertainty on \\(\\beta ^\*\\). Statistical uncertainties are negligible for most separation steps.
**Fig. 24**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/24)
Orbit offset caused by the beamâbeam interaction measured by DOROS BPMs in the horizontal plane at IP1 during the horizontal separation scan. Analytical prediction is shown for comparison \[[35](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR35 "M. Bassetti, G.A. Erskine, Closed expression for the electrical field of a two-dimensional gaussian charge. Technical report, CERN, Geneva (1980).
https://cds.cern.ch/record/122227
")\]
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/24)
#### 6\.4.3 Impact of the beam-beam effects on beam sizes and tune shift
Separation-dependent changes in beam size were measured using the synchrotron light monitor. Figure [25](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig25) shows good consistency across individual bunch measurements, as well as strong agreement with COMBI simulations. This reinforces the utility of the BSRT as a reliable diagnostic for beamâbeam effects. The uncertainties shown are statistical only and are treated as uncorrelated between scan points. The normalization to the zero-separation value is applied for visualization purposes and does not impose a constraint on that point.
**Fig. 25**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/25)
BSRT beam-size measurement in the vertical plane for each bunch (points) during horizontal (top) and vertical (bottom) separation scans at IP1. The exponential evolution correction is applied based on adjacent head-on points (before, in the middle and after the scan) to highlight the relative changes. COMBI predictions are shown (gray dashed line) for comparison
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/25)
Transverse tune shifts were extracted from ADT-measured spectra. As illustrated in Fig. [26](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig26), the mean tune shift evolution is reconstructed for both the scanning and non-scanning planes. The measured tune shifts closely follow the expected COMBI trends, further validating the simulation models used for beamâbeam corrections.
**Fig. 26**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/26)
Measured tune shift in units of the beam-beam parameter \\(\\xi \_{BB}\\) for the horizontal (blue) and vertical (orange) planes, during a horizontal separation scan at IP1. The error bars indicate typically assumed conservative empirical systematic error on the tune measurement of 0.001. The solid curves represent the COMBI prediction
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/26)
### 6\.5 Beam-beam signatures during vdM scans at \\(\\sqrt{s} = 13.6\\ \\text {TeV}\\)
During the first vdM calibration of LHC Run 3, an unusually high beamâbeam parameter was recorded, offering a unique opportunity to study beamâbeam effects under enhanced interaction conditions. These effects are examined using online ATLAS luminosity data, with representative results shown in Fig. [27](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig27).
The example focuses on a standard vdM scan pair conducted at CMS, while monitoring luminosity at the non-scanning IP (ATLAS). The plotted luminosity changes are normalized to the central scan point, where both IPs are in head-on collision. Variations in this normalized luminosity across the scan directly reflect the impact of beamâbeam interaction at the scanning IP on the luminosity at the non-scanning IP. These patterns are reproducible across multiple scan pairs acquired within the same fill. At the largest nominal separations (approximately \\(\\pm 0.6\\,\\)mm), the beams are sufficiently separated to eliminate beamâbeam effects entirely. Notably, a qualitative difference is observed between the horizontal and vertical scan planes, an asymmetry that is not present in the above presented dedicated experiment. This behavior is attributed to differences in the phase advance between IPs in the standard vdM optics configuration. For vertical (*y*) scans, experimental data agree well with simulation, confirming the robustness of the modeling under these conditions. In contrast, the horizontal (*x*) scans exhibit a systematic discrepancy: the observed beamâbeam impact is broader than that predicted by simulation. The model assumes that the separation steps are accurately defined using the beam size measured during the vdM scans, but this mismatch may suggest an unmodeled phase advance deviation or another optics imperfection. The total difference between data and simulation is within 0.2%. This dataset originates from a standard vdM calibration, for which only online luminosity measurement was available, and many contributions to the overall vdM uncertainty, such as bunch-to-bunch variations, averaging over different bunch families, and emittance evolution during the scan, are not disentangled. Therefore, the 0.2% difference should be interpreted as a consistency check under typical operational conditions rather than a precision validation of the beamâbeam model. Nevertheless, this indicates that the model remains sufficiently accurate for beamâbeam corrections in vdM analyses.
**Fig. 27**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/27)
Changes from the separation-dependent beam-beam effects at the non-scanning IP â ATLAS online luminosity shown with points. Data comes from a standard vdM scan pair during the 2022 *pp* luminosity calibration. Separation steps and direction are indicated with gray markers. The last CMS scan is shown, with the highest beam-beam parameter of \\( \\xi =5.3 \\times 10^{-3} \\). The corresponding COMBI simulation results are also shown (triangles)
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/27)
## 7 Conclusions
Beamâbeam effects were directly observed and quantitatively validated at the LHC. The impact of beamâbeam interactions on head-on luminosity, bunch sizes, tune spectra, and beam orbit positions was measured across various collision configurations. The observed effects align closely with COMBI simulations. In particular, the luminosity enhancement due to amplitude-dependent \\(\\beta \\)\-beating was measured for the first time and is found to agree with expectations over a wide range of the beamâbeam parameter \\(\\xi \\). Together with the observed linear scaling with respect to the beamâbeam parameter, this provides experimental confirmation of the numerical studies presented in \[[13](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR13 "A. Babaev, T. Barklow, O. Karacheban, W. Kozanecki, I. Kralik, A. Mehta, G. Pasztor, T. Pieloni, D. Stickland, C. Tambasco, R. Tomas, J. WaĆczyk, Impact of beam-beam effects on absolute luminosity calibrations at the CERN large hadron collider. Eur. Phys. J. C 84, 17 (2024).
https://doi.org/10.48550/arXiv.2306.10394
")\]. Consistent signatures were found across multiple scan configurations, including two- and three-IP setups, and the characteristic tune shifts were well reproduced by both measurement and simulation.
Accurate instrumentation and experimental design were essential. Key to the results was the stepwise collision scheme, isolating individual beamâbeam contributions, with the implementation of the witness IP. Synchrotron light measurements from the BSRT, especially in the vertical plane of Beam 2, showed excellent agreement with COMBI predictions and bunch-by-bunch consistency. Tune spectra from the ADT system provided direct access to mean tune shifts, while orbit data from the DOROS system enabled measurement of beamâbeam deflection down to 1 \\(\\upmu \\)m, matching analytical estimates. These tools, combined with precise control of IP configurations, enabled a comprehensive cross-check of beamâbeam models.
Some systematic limitations were identified. During one of the tests, beam stability during the scans was insufficient to allow for symmetric modeling of luminosity decay in both transverse planes. In several scans, the nominal separation steps were based on overestimated beam sizes, resulting in larger-than-expected scan intervals and less sensitivity at low separation. Beamâbeam deflection measurements showed sign reversals at high separation, attributed to nonlinear orbit distortions from hysteresis effects in superconducting magnets. Phase advance errors between IPs and diagnostic systems further reduced the accuracy of simulations in some configurations.
Targeted improvements would allow for higher precision and new insights. Finer separation steps, especially below \\(3\\,\\sigma \\), would improve sensitivity to nonlinear beamâbeam effects and allow more accurate comparison with models. Suppressing coherent beamâbeam modes via phase advance tuning of not only two but all three IPs would improve tune shift measurements. Moreover, such optimal phase advance configuration between the IPs can be exploited to enhance luminosity. This effect could be directly propagated into at least a few percent increase in the total collected integrated luminosity, as discussed in Ref. \[[36](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR36 "J. WaĆczyk, T. Pieloni, Optimizing beam-beam beta-beating for luminosity enhancement at the LHC. JACoW IPAC 2025, 018 (2025)")\]. Aligning the phase advance between IPs and diagnostic devices, particularly the BSRT, would reduce sensitivity to optics errors, enabling better beam size measurements for both beams and transverse planes. Moreover, the suppression of the coherent beamâbeam modes facilitates the observation and analysis of the central, incoherent part of the tune spectrum, which is otherwise obscured by the dominant coherent peaks.
Further experimental strategies should also focus on directly probing the beamâbeam interaction at its point of occurrence. This would require access to the full transverse beam distributions with high spatial resolution and sensitivity, going beyond projected sizes or centroids. Such measurements could directly reveal the nonlinear distortions and tails induced by the beamâbeam force. In parallel, validating correction schemes that parametrize multi-IP beamâbeam effects on luminosity via equivalent tune shifts remains essential. These models underpin practical luminosity corrections in complex fill configurations and would benefit from systematic benchmarking.
## Data Availability Statement
Data will be made available on reasonable request. \[Authorsâ comment: The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.\]
## Code Availability Statement
Code/software will be made available on reasonable request. \[Authorsâ comment: The code/software generated during and/or analysed during the current study is available from the corresponding author on reasonable request.\]
## References
1. S. Meer, Calibration of the effective beam height in the ISR. Technical report, CERN, Geneva (1968). <https://cds.cern.ch/record/296752>
2. C. Rubbia, Measurement of the luminosity of \\(p{{\\overline{p}}}\\) collider with a (generalized) Van der Meer Method. Technical report, CERN, Geneva (1977). <https://cds.cern.ch/record/1025746>
3. ATLAS collaboration, Luminosity determination in \\(pp\\) collisions at \\(\\sqrt{s}=13\\) TeV using the ATLAS detector at the LHC. Eur. Phys. J. C **83**(10), 982 (2023). <https://doi.org/10.1140/epjc/s10052-023-11747-w>. [arXiv:2212.09379](http://arxiv.org/abs/2212.09379)
4. CMS Collaboration, Precision luminosity measurement in proton-proton collisions at sqrts = 13 TeV with the CMS detector. Technical report, CERN, Geneva (2025). <https://cds.cern.ch/record/2940794>
5. W. Herr, T. Pieloni, Beam-Beam Effects. Technical report, CERN (2014). Contribution to the CASâCERN Accelerator School: Advanced Accelerator Physics Course, Trondheim, Norway, 18â29 Aug 2013, p. 29. <https://doi.org/10.5170/CERN-2014-009.431>. <https://cds.cern.ch/record/1982430>
6. A.W. Chao, Beam-beam instability. AIP Conf. Proc. **127**, 201â242 (1985). <https://doi.org/10.1063/1.35187>
[Article](https://doi.org/10.1063%2F1.35187) [ADS](http://adsabs.harvard.edu/cgi-bin/nph-data_query?link_type=ABSTRACT&bibcode=1985AIPC..127..201C) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Beam-beam%20instability&journal=AIP%20Conf.%20Proc.&doi=10.1063%2F1.35187&volume=127&pages=201-242&publication_year=1985&author=Chao%2CAW)
7. E.D. Courant, H.S. Snyder, Theory of the alternating-gradient synchrotron. Ann. Phys. **3**, 1â48 (1958). <https://doi.org/10.1016/0003-4916(58)90012-5>
[Article](https://doi.org/10.1016%2F0003-4916%2858%2990012-5) [ADS](http://adsabs.harvard.edu/cgi-bin/nph-data_query?link_type=ABSTRACT&bibcode=1958AnPhy...3....1C) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Theory%20of%20the%20alternating-gradient%20synchrotron&journal=Ann.%20Phys.&doi=10.1016%2F0003-4916%2858%2990012-5&volume=3&pages=1-48&publication_year=1958&author=Courant%2CED&author=Snyder%2CHS)
8. P. Bambade, R. Erickson, W.A. Koska, W. Kozanecki, N. Phinney, S.R. Wagner, Observation of beam-beam deflections at the interaction point of the SLAC Linear Collider. Phys. Rev. Lett. **62**(25), 2949â2952 (1989). <https://doi.org/10.1103/PhysRevLett.62.2949>
[Article](https://doi.org/10.1103%2FPhysRevLett.62.2949) [ADS](http://adsabs.harvard.edu/cgi-bin/nph-data_query?link_type=ABSTRACT&bibcode=1989PhRvL..62.2949B) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Observation%20of%20beam-beam%20deflections%20at%20the%20interaction%20point%20of%20the%20SLAC%20Linear%20Collider&journal=Phys.%20Rev.%20Lett.&doi=10.1103%2FPhysRevLett.62.2949&volume=62&issue=25&pages=2949-2952&publication_year=1989&author=Bambade%2CP&author=Erickson%2CR&author=Koska%2CWA&author=Kozanecki%2CW&author=Phinney%2CN&author=Wagner%2CSR)
9. W.A. Koska, P. Bambade, W. Kozanecki, N. Phinney, S.R. Wagner, Beam-beam deflection as a beam tuning tool at the slac linear collider. Nucl. Instrum. Methods A **286**, 32 (1990). <https://doi.org/10.1016/0168-9002(90)90203-I>
[Article](https://doi.org/10.1016%2F0168-9002%2890%2990203-I) [ADS](http://adsabs.harvard.edu/cgi-bin/nph-data_query?link_type=ABSTRACT&bibcode=1990NIMPA.286...32K) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Beam-beam%20deflection%20as%20a%20beam%20tuning%20tool%20at%20the%20slac%20linear%20collider&journal=Nucl.%20Instrum.%20Methods%20A&doi=10.1016%2F0168-9002%2890%2990203-I&volume=286&publication_year=1990&author=Koska%2CWA&author=Bambade%2CP&author=Kozanecki%2CW&author=Phinney%2CN&author=Wagner%2CSR)
10. C. Bovet, M.D. Hildreth, M. Lamont, H. Schmickler, J. Wenninger, Luminosity optimisation using beam-beam deflections at LEP, in *Conf. Proc. C* (1996). <https://cds.cern.ch/record/306910>
11. D. Brandt, W. Herr, M. Meddahi, A. Verdier, Is lep beam-beam limited at its highest energy? in *Proceedings of the 1999 Particle Accelerator Conference (PAC â99), New York*, pp. 3005â3007. IEEE/JACoW. THP25; CERN report CERN-SL-99-030-AP (1999). <https://accelconf.web.cern.ch/p99/PAPERS/THP25.pdf>
12. W. Kozanecki, T. Pieloni, J. Wenninger, Observation of Beam-beam Deflections with LHC Orbit Data. Technical report, CERN (2013). <https://cds.cern.ch/record/1581723?ln=en>
13. A. Babaev, T. Barklow, O. Karacheban, W. Kozanecki, I. Kralik, A. Mehta, G. Pasztor, T. Pieloni, D. Stickland, C. Tambasco, R. Tomas, J. WaĆczyk, Impact of beam-beam effects on absolute luminosity calibrations at the CERN large hadron collider. Eur. Phys. J. C **84**, 17 (2024). <https://doi.org/10.48550/arXiv.2306.10394>
[Article](https://doi.org/10.48550%2FarXiv.2306.10394) [ADS](http://adsabs.harvard.edu/cgi-bin/nph-data_query?link_type=ABSTRACT&bibcode=2024EPJC...84...17B) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Impact%20of%20beam-beam%20effects%20on%20absolute%20luminosity%20calibrations%20at%20the%20CERN%20large%20hadron%20collider&journal=Eur.%20Phys.%20J.%20C&doi=10.48550%2FarXiv.2306.10394&volume=84&publication_year=2024&author=Babaev%2CA&author=Barklow%2CT&author=Karacheban%2CO&author=Kozanecki%2CW&author=Kralik%2CI&author=Mehta%2CA&author=Pasztor%2CG&author=Pieloni%2CT&author=Stickland%2CD&author=Tambasco%2CC&author=Tomas%2CR&author=Wa%C5%84czyk%2CJ)
14. T. Pieloni, A study of beam-beam effects in hadron colliders with a large number of bunches. PhD Thesis, Ăcole Polytechnique FĂ©dĂ©rale de Lausanne (2008). <https://doi.org/10.5075/epfl-thesis-4211>
15. ATLAS Collaboration, Improved luminosity determination in pp collisions at \\(\\sqrt{s} = 7\\) TeV using the ATLAS detector at the LHC. Eur. Phys. J. C **73**, 2518 (2013). <https://doi.org/10.1140/epjc/s10052-013-2518-3>
16. C.M.S. Collaboration, Development of the cms detector for the cern lhc run 3. J. Instrum. **19**(05), 05064 (2024). <https://doi.org/10.1088/1748-0221/19/05/p05064>
[Article](https://doi.org/10.1088%2F1748-0221%2F19%2F05%2Fp05064) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Development%20of%20the%20cms%20detector%20for%20the%20cern%20lhc%20run%203&journal=J.%20Instrum.&doi=10.1088%2F1748-0221%2F19%2F05%2Fp05064&volume=19&issue=05&publication_year=2024&author=Collaboration%2CCMS)
17. C. Barschel, M. Ferro-Luzzi, J.-J. Gras, M. Ludwig, P. Odier, S. Thoulet, Results of the LHC DCCT calibration studies. Technical Report CERN-ATS-Note-2012-026 PERF (2012). <https://cds.cern.ch/record/1425904?ln=en>
18. G. Anders et al., Study of the relative LHC bunch populations for luminosity calibration. Technical Report. CERN-ATS-Note-2012-028 PERF (2012). <https://cds.cern.ch/record/1427726>
19. G. Trad, Development and Optimisation of the SPS and LHC Beam Diagnostics Based on Synchrotron Radiation Monitors Presented 22 Jan 2015. Presented 22 Jan 2015. <https://cds.cern.ch/record/2266055>
20. R. TomĂĄs, O. BrĂŒning, M. Giovannozzi, P. Hagen, M. Lamont, F. Schmidt, G. Vanbavinckhove, M. Aiba, R. Calaga, R. Miyamoto, Cern large hadron collider optics model, measurements, and corrections. Phys. Rev. ST Accel. Beams **13**, 121004 (2010). <https://doi.org/10.1103/PhysRevSTAB.13.121004>
21. R. TomĂĄs, M. Aiba, A. Franchi, U. Iriso, Review of linear optics measurement and correction for charged particle accelerators. Phys. Rev. Accel. Beams **20**, 054801 (2017). <https://doi.org/10.1103/PhysRevAccelBeams.20.054801>
[Article](https://doi.org/10.1103%2FPhysRevAccelBeams.20.054801) [ADS](http://adsabs.harvard.edu/cgi-bin/nph-data_query?link_type=ABSTRACT&bibcode=2017PhRvS..20e4801T) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Review%20of%20linear%20optics%20measurement%20and%20correction%20for%20charged%20particle%20accelerators&journal=Phys.%20Rev.%20Accel.%20Beams&doi=10.1103%2FPhysRevAccelBeams.20.054801&volume=20&publication_year=2017&author=Tom%C3%A1s%2CR&author=Aiba%2CM&author=Franchi%2CA&author=Iriso%2CU)
22. O. Karacheban, P. Tsrunchev, on behalf of CMS, Emittance scans for cms luminosity calibration. EPJ Web Conf. **201**, 04001 (2019). <https://doi.org/10.1051/epjconf/201920104001>
23. F. Carlier, R. TomĂĄs, Accuracy and feasibility of the \\({\\beta }^{\*}\\) measurement for lhc and high luminosity lhc using \\(k\\) modulation. Phys. Rev. Accel. Beams **20**, 011005 (2017). <https://doi.org/10.1103/PhysRevAccelBeams.20.011005>
[Article](https://doi.org/10.1103%2FPhysRevAccelBeams.20.011005) [ADS](http://adsabs.harvard.edu/cgi-bin/nph-data_query?link_type=ABSTRACT&bibcode=2017PhRvS..20a1005C) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Accuracy%20and%20feasibility%20of%20the%20%24%24%7B%5Cbeta%20%7D%5E%7B%2A%7D%24%24%20%CE%B2%20%E2%88%97%20measurement%20for%20lhc%20and%20high%20luminosity%20lhc%20using%20%24%24k%24%24%20k%20modulation&journal=Phys.%20Rev.%20Accel.%20Beams&doi=10.1103%2FPhysRevAccelBeams.20.011005&volume=20&publication_year=2017&author=Carlier%2CF&author=Tom%C3%A1s%2CR)
24. M. SöderĂ©n, J. Komppula, G. Kotzian, S. Rains, D. Valuch, ADT and Obsbox in LHC Run 2, plans for LS2, pp. 165â171 (2019)
25. A. Boccardi, M. Gasior, R. Jones, P. Karlsson, R. Steinhagen, First results from the lhc bbq tune and chromaticity systems. Technical report, CERN, Geneva (2009). <https://cds.cern.ch/record/1156349>
26. M. Gasior, G. Baud, J. Olexa, G. Valentino, First Operational Experience with the LHC Diode ORbit and OScillation (DOROS) System, p. 07 (2017) <https://doi.org/10.18429/JACoW-IBIC2016-MOPG07>
27. A.W. Chao, Coherent beam-beam effects, in *Frontiers of Particle Beams: Intensity Limitations* (Springer, Berlin, Heidelberg, 1992), pp. 363â414. <https://doi.org/10.1007/3-540-55250-2_36>
28. J. WaĆczyk, A. Dabrowski, D. Stickland, T. Pieloni, W. Kozanecki, X. Buffat, Impact of multiple beam-beam encounters on LHC absolute-luminosity calibrations by the van der Meer method. JACoW IPAC **2023**, 053 (2023). <https://doi.org/10.18429/JACoW-IPAC2023-WEPA053>
[Article](https://doi.org/10.18429%2FJACoW-IPAC2023-WEPA053) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Impact%20of%20multiple%20beam-beam%20encounters%20on%20LHC%20absolute-luminosity%20calibrations%20by%20the%20van%20der%20Meer%20method&journal=JACoW%20IPAC&doi=10.18429%2FJACoW-IPAC2023-WEPA053&volume=2023&publication_year=2023&author=Wa%C5%84czyk%2CJ&author=Dabrowski%2CA&author=Stickland%2CD&author=Pieloni%2CT&author=Kozanecki%2CW&author=Buffat%2CX)
29. J. WaĆczyk, Precision luminosity measurement at hadron colliders. PhD Thesis, Ăcole Polytechnique FĂ©dĂ©rale de Lausanne (2024). <https://doi.org/10.5075/epfl-thesis-10500>
30. T. Pieloni, J. Barranco, X. Buffat, P.C. Jorge, L.M. Medrano, C. Tambasco, R. TomĂĄs, Dynamic beta and beta-beating effects in the presence of the beam-beam interactions in *57th ICFA Advanced Beam Dynamics Workshop on High-Intensity and High-Brightness Hadron Beams*, p. 027 (2016). <https://doi.org/10.18429/JACoW-HB2016-MOPR027>
31. F. Carlier, et al. Challenges of K-modulation measurements in the LHC Run 3, in *Proceedings of IPACâ23. IPACâ23â14th International Particle Accelerator Conference* (JACoW Publishing, Geneva, Venice, 2023), pp. 531â534. <https://doi.org/10.18429/JACoW-IPAC2023-MOPL014>
32. Y. Alexahin, A study of the coherent beam-beam effect in the framework of the Vlasov perturbation theory. Nucl. Instrum. Methods Phys. Res., Sect. A **480**(1â3), 253â288 (2002). <https://doi.org/10.1016/S0168-9002(01)01219-0>
[Article](https://doi.org/10.1016%2FS0168-9002%2801%2901219-0) [ADS](http://adsabs.harvard.edu/cgi-bin/nph-data_query?link_type=ABSTRACT&bibcode=2002NIMPA.480..253A) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=A%20study%20of%20the%20coherent%20beam-beam%20effect%20in%20the%20framework%20of%20the%20Vlasov%20perturbation%20theory&journal=Nucl.%20Instrum.%20Methods%20Phys.%20Res.%2C%20Sect.%20A&doi=10.1016%2FS0168-9002%2801%2901219-0&volume=480&issue=1%E2%80%933&pages=253-288&publication_year=2002&author=Alexahin%2CY)
33. R. Calaga, W. Herr, G. Papotti, T. Pieloni, X. Buffat, S. White, R. Giachino, Coherent beam-beam mode in the LHC. CERN (2014). <https://doi.org/10.5170/CERN-2014-004.227>. <http://cds.cern.ch/record/1957039>
34. K. Yokoya, Y. Funakoshi, E. Kikutani, H. Koiso, J. Urakawa, Tune shift of coherent beam-beam oscillations. Part. Accel. **27**, 181â186 (1990)
[Google Scholar](http://scholar.google.com/scholar_lookup?&title=Tune%20shift%20of%20coherent%20beam-beam%20oscillations&journal=Part.%20Accel.&volume=27&pages=181-186&publication_year=1990&author=Yokoya%2CK&author=Funakoshi%2CY&author=Kikutani%2CE&author=Koiso%2CH&author=Urakawa%2CJ)
35. M. Bassetti, G.A. Erskine, Closed expression for the electrical field of a two-dimensional gaussian charge. Technical report, CERN, Geneva (1980). <https://cds.cern.ch/record/122227>
36. J. WaĆczyk, T. Pieloni, Optimizing beam-beam beta-beating for luminosity enhancement at the LHC. JACoW IPAC **2025**, 018 (2025)
[Google Scholar](http://scholar.google.com/scholar_lookup?&title=Optimizing%20beam-beam%20beta-beating%20for%20luminosity%20enhancement%20at%20the%20LHC&journal=JACoW%20IPAC&volume=2025&publication_year=2025&author=Wa%C5%84czyk%2CJ&author=Pieloni%2CT)
[Download references](https://citation-needed.springer.com/v2/references/10.1140/epjc/s10052-026-15476-8?format=refman&flavour=references)
## Acknowledgements
We would especially like to acknowledge W. Kozanecki for his substantial contributions to the development of this work and for carefully reviewing the manuscript, and Y. Wu for participation in the experiment and online data checks. We are also indebted to the ATLAS and CMS luminosity experts, in particular R. Hawkings, E. Torrence and A. Shevelev. We thank LHC experts for help with setting up and carrying out the experiment, in particular M. Hostettler, T. Persson, M. Solfaroli Camillocci, G. Trad, J. Wenninger. As well as optics team, in particular F. Carlier and M. Le Garrec. Some of the authors have received support from the Swiss Accelerator Research and Technology Institute (CHART).
## Author information
### Authors and Affiliations
1. CERN, Geneva, Switzerland
Joanna WaĆczyk, Xavier Buffat, Anne Dabrowski & Rogelio Tomas
2. Laboratory of Particle Accelerator Physics, EPFL, Lausanne, Switzerland
Joanna WaĆczyk & Tatiana Pieloni
3. Princeton University, Princeton, NJ, USA
David Stickland
Authors
1. Joanna WaĆczyk
[View author publications](https://link.springer.com/search?sortBy=newestFirst&contributor=Joanna%20Wa%C5%84czyk)
Search author on:[PubMed](https://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=search&term=Joanna%20Wa%C5%84czyk) [Google Scholar](https://scholar.google.co.uk/scholar?as_q=&num=10&btnG=Search+Scholar&as_epq=&as_oq=&as_eq=&as_occt=any&as_sauthors=%22Joanna%20Wa%C5%84czyk%22&as_publication=&as_ylo=&as_yhi=&as_allsubj=all&hl=en)
2. Tatiana Pieloni
[View author publications](https://link.springer.com/search?sortBy=newestFirst&contributor=Tatiana%20Pieloni)
Search author on:[PubMed](https://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=search&term=Tatiana%20Pieloni) [Google Scholar](https://scholar.google.co.uk/scholar?as_q=&num=10&btnG=Search+Scholar&as_epq=&as_oq=&as_eq=&as_occt=any&as_sauthors=%22Tatiana%20Pieloni%22&as_publication=&as_ylo=&as_yhi=&as_allsubj=all&hl=en)
3. Xavier Buffat
[View author publications](https://link.springer.com/search?sortBy=newestFirst&contributor=Xavier%20Buffat)
Search author on:[PubMed](https://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=search&term=Xavier%20Buffat) [Google Scholar](https://scholar.google.co.uk/scholar?as_q=&num=10&btnG=Search+Scholar&as_epq=&as_oq=&as_eq=&as_occt=any&as_sauthors=%22Xavier%20Buffat%22&as_publication=&as_ylo=&as_yhi=&as_allsubj=all&hl=en)
4. Anne Dabrowski
[View author publications](https://link.springer.com/search?sortBy=newestFirst&contributor=Anne%20Dabrowski)
Search author on:[PubMed](https://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=search&term=Anne%20Dabrowski) [Google Scholar](https://scholar.google.co.uk/scholar?as_q=&num=10&btnG=Search+Scholar&as_epq=&as_oq=&as_eq=&as_occt=any&as_sauthors=%22Anne%20Dabrowski%22&as_publication=&as_ylo=&as_yhi=&as_allsubj=all&hl=en)
5. David Stickland
[View author publications](https://link.springer.com/search?sortBy=newestFirst&contributor=David%20Stickland)
Search author on:[PubMed](https://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=search&term=David%20Stickland) [Google Scholar](https://scholar.google.co.uk/scholar?as_q=&num=10&btnG=Search+Scholar&as_epq=&as_oq=&as_eq=&as_occt=any&as_sauthors=%22David%20Stickland%22&as_publication=&as_ylo=&as_yhi=&as_allsubj=all&hl=en)
6. Rogelio Tomas
[View author publications](https://link.springer.com/search?sortBy=newestFirst&contributor=Rogelio%20Tomas)
Search author on:[PubMed](https://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=search&term=Rogelio%20Tomas) [Google Scholar](https://scholar.google.co.uk/scholar?as_q=&num=10&btnG=Search+Scholar&as_epq=&as_oq=&as_eq=&as_occt=any&as_sauthors=%22Rogelio%20Tomas%22&as_publication=&as_ylo=&as_yhi=&as_allsubj=all&hl=en)
### Corresponding author
Correspondence to [Joanna WaĆczyk](mailto:jwanczyk@cern.ch).
## Appendix A: \\(\\beta \\)-beating derivative
### Appendix A: \\(\\beta \\)\-beating derivative
Derivative of Eq. ([7](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Equ7)) with respect to the phase advance \\(\\mu \_{1,x}\\), which characterizes the change in the linear beam-beam \\(\\beta \\)\-beating due to variations in phase advance:

## Rights and permissions
**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit <http://creativecommons.org/licenses/by/4.0/>.
Funded by SCOAP3.
[Reprints and permissions](https://s100.copyright.com/AppDispatchServlet?title=First%20measurements%20of%20beam-beam%20effects%20in%20beam-separation%2C%20luminosity-calibration%20scans%20at%20the%20LHC&author=Joanna%20Wa%C5%84czyk%20et%20al&contentID=10.1140%2Fepjc%2Fs10052-026-15476-8©right=The%20Author%28s%29&publication=1434-6052&publicationDate=2026-03-17&publisherName=SpringerNature&orderBeanReset=true&oa=CC%20BY)
## About this article
[](https://crossmark.crossref.org/dialog/?doi=10.1140/epjc/s10052-026-15476-8)
### Cite this article
WaĆczyk, J., Pieloni, T., Buffat, X. *et al.* First measurements of beam-beam effects in beam-separation, luminosity-calibration scans at the LHC. *Eur. Phys. J. C* **86**, 282 (2026). https://doi.org/10.1140/epjc/s10052-026-15476-8
[Download citation](https://citation-needed.springer.com/v2/references/10.1140/epjc/s10052-026-15476-8?format=refman&flavour=citation)
- Received: 13 October 2025
- Accepted: 21 February 2026
- Published: 17 March 2026
- Version of record: 17 March 2026
- DOI: https://doi.org/10.1140/epjc/s10052-026-15476-8
### Share this article
Anyone you share the following link with will be able to read this content:
Get shareable link
Sorry, a shareable link is not currently available for this article.
Copy shareable link to clipboard
Provided by the Springer Nature SharedIt content-sharing initiative
Advertisement
## Search
## Navigation
- [Find a journal](https://link.springer.com/journals/)
- [Publish with us](https://www.springernature.com/gp/authors)
- [Track your research](https://link.springernature.com/home/)
### Discover content
- [Journals A-Z](https://link.springer.com/journals/a/1)
- [Books A-Z](https://link.springer.com/books/a/1)
### Publish with us
- [Journal finder](https://link.springer.com/journals)
- [Publish your research](https://www.springernature.com/gp/authors)
- [Language editing](https://authorservices.springernature.com/go/sn/?utm_source=SNLinkfooter&utm_medium=Web&utm_campaign=SNReferral)
- [Open access publishing](https://www.springernature.com/gp/open-science/about/the-fundamentals-of-open-access-and-open-research)
### Products and services
- [Our products](https://www.springernature.com/gp/products)
- [Librarians](https://www.springernature.com/gp/librarians)
- [Societies](https://www.springernature.com/gp/societies)
- [Partners and advertisers](https://www.springernature.com/gp/partners)
### Our brands
- [Springer](https://link.springer.com/brands/springer)
- [Nature Portfolio](https://www.nature.com/)
- [BMC](https://link.springer.com/brands/bmc)
- [Palgrave Macmillan](https://link.springer.com/brands/palgrave)
- [Apress](https://link.springer.com/brands/apress)
- [Discover](https://link.springer.com/brands/discover)
- Your privacy choices/Manage cookies
- [Your US state privacy rights](https://www.springernature.com/gp/legal/ccpa)
- [Accessibility statement](https://link.springer.com/accessibility)
- [Terms and conditions](https://link.springer.com/termsandconditions)
- [Privacy policy](https://link.springer.com/privacystatement)
- [Help and support](https://support.springernature.com/en/support/home)
- [Legal notice](https://link.springer.com/legal-notice)
- [Cancel contracts here](https://support.springernature.com/en/support/solutions/articles/6000255911-subscription-cancellations)
Not affiliated
[](https://www.springernature.com/)
© 2026 Springer Nature |
| Readable Markdown | ## 1 Introduction
A precise determination of the absolute luminosity scale is essential for a wide range of LHC measurements, as it directly affects the normalization of many key physics cross sections. At the LHC, this calibration relies primarily on the van der Meer (vdM) method \[[1](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR1 "S. Meer, Calibration of the effective beam height in the ISR. Technical report, CERN, Geneva (1968).
https://cds.cern.ch/record/296752
"), [2](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR2 "C. Rubbia, Measurement of the luminosity of
$$p{{\overline{p}}}$$
p
p
ÂŻ
collider with a (generalized) Van der Meer Method. Technical report, CERN, Geneva (1977).
https://cds.cern.ch/record/1025746
")\], which uses dedicated beam-separation scans performed under specially tailored beam conditions to relate measured interaction rates to the absolute luminosity inferred from beam parameters. While sub-percent uncertainties on the integrated luminosity have been achieved for Run-2 data by ATLAS \[[3](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR3 "ATLAS collaboration, Luminosity determination in
$$pp$$
pp
collisions at
$$\sqrt{s}=13$$
s
=
13
TeV using the ATLAS detector at the LHC. Eur. Phys. J. C 83(10), 982 (2023).
https://doi.org/10.1140/epjc/s10052-023-11747-w
.
arXiv:2212.09379
")\] and CMS \[[4](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR4 "CMS Collaboration, Precision luminosity measurement in proton-proton collisions at sqrts = 13 TeV with the CMS detector. Technical report, CERN, Geneva (2025).
https://cds.cern.ch/record/2940794
")\], further improvements and robust experimental validation of beam-dynamical effects remain essential.
For relativistic proton beams colliding with zero crossing angle, which is the typical configuration for the *vdM* calibration, the luminosity at a given interaction point (IP) of the LHC is proportional to the overlap integral of the particle-density distributions \\(\\rho \_{\\text {1,i}}\\), \\(\\rho \_{\\text {2,i}}\\) in bunch-pair *i* \[[5](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR5 "W. Herr, T. Pieloni, Beam-Beam Effects. Technical report, CERN (2014). Contribution to the CASâCERN Accelerator School: Advanced Accelerator Physics Course, Trondheim, Norway, 18â29 Aug 2013, p. 29.
https://doi.org/10.5170/CERN-2014-009.431
.
https://cds.cern.ch/record/1982430
")\]:
\$\$\\begin{aligned} {\\mathcal {L}}\_{inst}&= 2c f\_{\\text {rev}} \\sum ^{N\_b}\_i n\_{\\text {1,i}}, n\_{\\text {2,i}} \\nonumber \\\\&\\quad \\times \\iiiint ^{+\\infty }\_{-\\infty }\\rho \_{\\text {1,i}}(x,y,z-ct) \\nonumber \\\\&\\quad \\times \\rho \_{\\text {2,i}}(x,y,z+ct)\\,dxdydzdt, \\end{aligned}\$\$
(1)
where *c* is the speed of light, \\(f\_{\\text {rev}}\\) is the LHC revolution frequency, \\(N\_b\\) is the number of colliding bunches, and \\(n\_{\\text {1,i}}, n\_{\\text {2,i}}\\) are the corresponding total charges per bunch for each colliding pair. Although only a very small fraction of the beam particles actually collide, the two opposing beams interact with each other electromagnetically, a dynamical process known as the beam-beam interaction \[[5](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR5 "W. Herr, T. Pieloni, Beam-Beam Effects. Technical report, CERN (2014). Contribution to the CASâCERN Accelerator School: Advanced Accelerator Physics Course, Trondheim, Norway, 18â29 Aug 2013, p. 29.
https://doi.org/10.5170/CERN-2014-009.431
.
https://cds.cern.ch/record/1982430
"), [6](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR6 "A.W. Chao, Beam-beam instability. AIP Conf. Proc. 127, 201â242 (1985).
https://doi.org/10.1063/1.35187
")\]. In the LHC, as in any other synchrotron, the beams are not continuous but are divided into discrete âbunchesâ, each with a high proton density, containing approximately \\(10^{11}\\) protons over a bunch length of 7.5â9 cm (in terms of RMS). This electromagnetic interaction occurs when the two beams share a common beam pipe, and is governed by a non-linear force that depends on the radial distance *r* of a test particle to the center of the opposing âsourceâ bunch with *n* protons Gaussian-distributed in the transverse planes (valid approximation for the LHC beams):
\$\$\\begin{aligned} F = - \\frac{ne^2}{2\\pi \\epsilon \_0r}\\biggl (1-\\exp \\biggl \[ -\\frac{r^2}{2\\sigma ^2}\\biggr \] \\biggr ), \\end{aligned}\$\$
(2)
where *e* is the elementary proton charge, \\(\\epsilon \_0\\) is the vacuum permittivity, and \\(\\sigma \\) is the RMS radius of the transverse charge distribution within the bunch, in a simplified case of a round shape \\((\\sigma = \\sigma \_x = \\sigma \_y)\\). At very small amplitudes \\((r\\rightarrow 0)\\) the force can be linearized giving the expression for the so called beam-beam parameter, which is often used to assess the strength of the force:
\$\$\\begin{aligned} \\xi = \\frac{nr\_p\\beta ^\*}{4\\pi \\gamma \\sigma ^2} = \\frac{nr\_p}{4\\pi \\epsilon \_n}, \\end{aligned}\$\$
(3)
where \\(r\_p\\) is the classical proton radius, \\(\\beta ^\*\\) is the optical beta-function at the interaction point (IP) â one of the CourantâSnyder parameters \[[7](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR7 "E.D. Courant, H.S. Snyder, Theory of the alternating-gradient synchrotron. Ann. Phys. 3, 1â48 (1958).
https://doi.org/10.1016/0003-4916(58)90012-5
")\], \\(\\gamma \\) is the Lorentz factor, and \\(\\epsilon \_n\\) is the normalized emittance.
Beam-beam effects were first observed almost 40 years ago \[[8](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR8 "P. Bambade, R. Erickson, W.A. Koska, W. Kozanecki, N. Phinney, S.R. Wagner, Observation of beam-beam deflections at the interaction point of the SLAC Linear Collider. Phys. Rev. Lett. 62(25), 2949â2952 (1989).
https://doi.org/10.1103/PhysRevLett.62.2949
")\] at the lepton colliders, and have been used extensively for accelerator diagnostics and optimization at both SLC and LEP \[[9](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR9 "W.A. Koska, P. Bambade, W. Kozanecki, N. Phinney, S.R. Wagner, Beam-beam deflection as a beam tuning tool at the slac linear collider. Nucl. Instrum. Methods A 286, 32 (1990).
https://doi.org/10.1016/0168-9002(90)90203-I
"),[10](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR10 "C. Bovet, M.D. Hildreth, M. Lamont, H. Schmickler, J. Wenninger, Luminosity optimisation using beam-beam deflections at LEP, in Conf. Proc. C (1996).
https://cds.cern.ch/record/306910
"),[11](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR11 "D. Brandt, W. Herr, M. Meddahi, A. Verdier, Is lep beam-beam limited at its highest energy? in Proceedings of the 1999 Particle Accelerator Conference (PAC â99), New York, pp. 3005â3007. IEEE/JACoW. THP25; CERN report CERN-SL-99-030-AP (1999).
https://accelconf.web.cern.ch/p99/PAPERS/THP25.pdf
")\]. At the LHC, they were first observed in late Run 1 \[[12](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR12 "W. Kozanecki, T. Pieloni, J. Wenninger, Observation of Beam-beam Deflections with LHC Orbit Data. Technical report, CERN (2013).
https://cds.cern.ch/record/1581723?ln=en
")\], but only in terms of the beam-beam induced coherent deflections. Single particle effects, in contrast, were presumed too small in hadron colliders that until recently, there was no motivation to measure them with a precision that would be meaningful in the context evoked here.
Nonetheless, the full complexity of this interaction must be accounted for when performing the luminosity calibration via the vdM method, using the separation scans. These are designed to measure the detector-specific constant \\(\\sigma \_{vis}\\) that relates the observed rate \\(\\mu ^{vis}\\) to the absolute instantaneous luminosity \\({\\mathcal {L}}\_{inst}\\):
\$\$\\begin{aligned} \\mu ^{vis} = \\frac{{\\mathcal {L}}\_{inst}\\sigma \_{vis}}{f\_{rev}}. \\end{aligned}\$\$
(4)
Under the assumption of uncorrelated particle densities in *x* and *y* planes, the transverse convolved bunch widths \\(\\Sigma \_x, \\Sigma \_y\\) can be extracted from the measured beam-separation dependence of the collision rate in the corresponding direction. The combined information from these scans, and bunch intensities can be used to calculate the instantaneous luminosity at the head-on position. The proportionality between the measured rate at the centered head-on position \\(\\mu \_{pk}\\) and the reconstructed luminosity from measured bunch parameters during the vdM calibration defines the luminometer-specific visible cross-section:
\$\$\\begin{aligned} \\sigma \_{vis} = \\frac{2\\pi \\Sigma \_x\\Sigma \_y}{n\_1 n\_2} \\mu \_{pk}. \\end{aligned}\$\$
(5)
If left uncorrected, the \\(\\sigma \_{vis}\\) measurement by the *vdM* method is biased by the beam-separation dependence of the mutual electromagnetic interaction of the two beams: the colliding bunches experience deflection-induced orbit shifts, as well as optical distortions akin to the dynamic-\\(\\beta \\) effect, that both depend on the transverse beam separation and must therefore be accounted for when deriving the absolute luminosity scale. The correction strategy discussed extensively in Ref. \[[13](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR13 "A. Babaev, T. Barklow, O. Karacheban, W. Kozanecki, I. Kralik, A. Mehta, G. Pasztor, T. Pieloni, D. Stickland, C. Tambasco, R. Tomas, J. WaĆczyk, Impact of beam-beam effects on absolute luminosity calibrations at the CERN large hadron collider. Eur. Phys. J. C 84, 17 (2024).
https://doi.org/10.48550/arXiv.2306.10394
")\] relies on simulation studies only, and thus requires experimental validation to establish its reliability.
The present report summarizes a campaign of dedicated measurements carried out in the spring of 2022 at the LHC, aimed at confirming the accuracy of the simulation studies detailed in Ref. \[[13](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR13 "A. Babaev, T. Barklow, O. Karacheban, W. Kozanecki, I. Kralik, A. Mehta, G. Pasztor, T. Pieloni, D. Stickland, C. Tambasco, R. Tomas, J. WaĆczyk, Impact of beam-beam effects on absolute luminosity calibrations at the CERN large hadron collider. Eur. Phys. J. C 84, 17 (2024).
https://doi.org/10.48550/arXiv.2306.10394
")\], and at validating the beam-beam correction strategy used in the luminosity-calibration analyses of the ALICE, ATLAS, CMS and LHCb experiments. More specifically, the goal is to quantify the impact of the beam-beam interaction on the tune spectra, orbit, the transverse-density distribution of the colliding bunches, as well as, for the first time, on the luminosity, by systematically varying the strength or number of the beam-beam interactions. The measurements are repeated for different values of the beam-beam parameter to verify the scaling proposed in \[[13](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR13 "A. Babaev, T. Barklow, O. Karacheban, W. Kozanecki, I. Kralik, A. Mehta, G. Pasztor, T. Pieloni, D. Stickland, C. Tambasco, R. Tomas, J. WaĆczyk, Impact of beam-beam effects on absolute luminosity calibrations at the CERN large hadron collider. Eur. Phys. J. C 84, 17 (2024).
https://doi.org/10.48550/arXiv.2306.10394
")\], for various choices of the scanning IP and for several multi-IP configurations. The beam conditions are chosen to be representative of *vdM* calibration sessions at the LHC, but optimized so as to maximize both the sensitivity of the measurements, and the operational efficiency.
This paper is organized as follows. The measurement strategy is detailed in Sect. [2](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec2), followed by the main ingredients of this accelerator experiment: beam conditions (Sect. [3](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec3)), luminometers and accelerator instrumentation (Sect. [4](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec4)), and optimization of the ring lattice to maximize the sensitivity of the measurements to beamâbeam-induced optical distortions of the colliding bunches (Sect. [5](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec11)). Section [6](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec19) is devoted to the experimental characterization of the impact of beam-beam effects on the tune spectra, on the luminosity at different IPs, as well as on other observables such as transverse single-bunch widths and IP orbits. At each step, the results are confronted with the predictions of the COherent Multibunch Beam-beam Interactions (COMBI) tracking code \[[14](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR14 "T. Pieloni, A study of beam-beam effects in hadron colliders with a large number of bunches. PhD Thesis, Ăcole Polytechnique FĂ©dĂ©rale de Lausanne (2008).
https://doi.org/10.5075/epfl-thesis-4211
")\], that produced most of the results detailed in Ref. \[[13](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR13 "A. Babaev, T. Barklow, O. Karacheban, W. Kozanecki, I. Kralik, A. Mehta, G. Pasztor, T. Pieloni, D. Stickland, C. Tambasco, R. Tomas, J. WaĆczyk, Impact of beam-beam effects on absolute luminosity calibrations at the CERN large hadron collider. Eur. Phys. J. C 84, 17 (2024).
https://doi.org/10.48550/arXiv.2306.10394
")\]. A collection of macroparticles is simulated turn by turn, based on a simplified model of the accelerator lattice using linear transfer and the self-consistent computation of the beam-beam forces at the IP. The conclusions are presented in Sect. [7](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec33).
## 2 Measurement strategy
The beam-beam interaction at an interaction point, referred to here as the *scanning* IP, induces measurable changes in key beam parameters such as the tune spectrum, closed orbit, and transverse single-beam profiles. These effects can be observed, albeit with limited precision, using standard accelerator diagnostics by comparing measurements obtained under varying transverse beam separation at the scanning IP. However, the resulting modifications to the transverse density distribution of the colliding bunches remain below the sensitivity threshold of conventional instrumentation, including synchrotron-radiation-based beam-profile monitors. Beam-beam-induced variations typically lie in the 0.5â1% range, small enough that only a highly precise measurement of collision-rate changes can provide adequate sensitivity.
**Fig. 1**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/1)
Scheme of the LHC and its interaction points (IP) with indicated IP1 as witness IP while performing a scan at IP5
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/1)
Importantly, the beam-beam interaction induces changes in beam properties that are intrinsically entangled with the *pp* collisions when observing luminosity at the scanning IP. Consequently, under realistic collider conditions, it is not possible to directly access a reference luminosity signal that is free of beam-beam effects â unlike in simulation studies, where such effects can be explicitly disabled. In a multi-collision configuration, however, luminosity shifts at one or more other IPs (Fig. [1](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig1)), where the beams are kept in continuous head-on collision and therefore act as non-scanning, or *witness*, IPs, can be used to monitor the separation-dependent beamâbeam effects. These effects are induced at the scanning IP, where the beams are deliberately brought in and out of collision, and propagate around the rings to the witness IPs, where they manifest as changes in the measured luminosity. In principle, any of the four LHC IPs can serve as a scanning IP, and any non-scanning IP can be designated as the witness IP for a given measurement. In practice, instrumental and operational constraints at the time the experiment was carried out restrict the choice of witness IP to the ATLAS (IP1) and CMS (IP5) collision points; IPs 1, 2 or 5, or combinations thereof, are used as scanning IP(s). Collisions at IP8 were deliberately avoided because the large longitudinal offset of the LHCb collision point, that breaks the eight-fold symmetry of the LHC rings, precludes the possibility of both members of a colliding-bunch pair to collide at all four IPs. In addition, the large crossing angle at IP8 would complicate the interpretation of the beam-beam effects that would occur at that IP if some of the bunches collided there.
The propagation of the beamâbeam-induced, amplitude-dependent \\(\\beta \\)\-beating from the scanning to a witness IP is controlled by the betatron phase advance between these two IPs. Since the targeted signatures are delicate to measure at best, it is natural to try and enhance them by adjusting this phase advance so as to maximize the sensitivity of luminosity shifts at the witness IP to beam-beam effects at the scanning IP; in doing so, however, the overall tunes must be preserved. Given the central and symmetric roles played by IP1 and IP5 in this experiment, the adjustment of their relative phase advance drove the optimization procedure that is detailed in Sect. [5](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec11), and that achieved a threefold improvement in measurement sensitivity.
Central to the quality of the measurements is the controlled variation of the beam-beam parameter at the scanning IP. This can be achieved either:
- by using *step scans*, *i.e.* by fully separating the beams at the scanning IP in either the horizontal or the vertical plane, and then bringing them back into head-on collision; or
- by using *separation scans*, *i.e.* by scanning the beams transversely with respect to each other in either the horizontal or the vertical plane; or
- by taking advantage of the natural beam-intensity decay and emittance growth to progressively reduce the beam-beam parameter.
Equally important is to ensure that the beams remain in head-on collision at the witness IP(s), such that reproducible luminosity shifts measured at these locations can be unambiguously correlated with controlled changes in the strength of the beam-beam interaction at the scanning IP. It was verified that the beam-beam-induced orbit shift and the potential non-closure of the orbit bumps used to control the beam separation at the scanning IP, did not significantly affect the actual separation, and thereby the measured collision rate, at the witness IP. In a few cases however, orbit drifts of uncontrolled origin ended up degrading the quality of some of the measurements.
Simulations demonstrate \[[13](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR13 "A. Babaev, T. Barklow, O. Karacheban, W. Kozanecki, I. Kralik, A. Mehta, G. Pasztor, T. Pieloni, D. Stickland, C. Tambasco, R. Tomas, J. WaĆczyk, Impact of beam-beam effects on absolute luminosity calibrations at the CERN large hadron collider. Eur. Phys. J. C 84, 17 (2024).
https://doi.org/10.48550/arXiv.2306.10394
")\] that multi-IP effects strongly influence the magnitude of beam-beam biases to vdM calibrations. Their characterization, therefore, constitutes an essential component of this experiment. The scan protocol is a generalization of that in the two-IP case. One IP (say, IP1) is designated as the witness IP, with beams colliding head-on there and at IP2, and IP5; the beams are then taken out of collision at IP2 only, then at both IP2 and IP5, and finally returned to a three-collision configuration in the reverse order. Since the phase-advance combinations are different in each configuration, the corresponding luminosity shifts predicted at the witness IP are also different and can be meaningfully confronted with the data.
## 3 Beam conditions
In the round-beam approximation, the beam-beam parameter \\(\\xi \\) does not depend on the beam energy (Eq. ([3](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Equ3))). Therefore, and in order to minimize the overhead associated with commissioning and operating the LHC in a non-standard optical configuration, as well as to allow for frequent refilling with different bunch patterns, the measurements were carried out at injection energy (\\(450\\ \\text {GeV}\\) per beam), and under non-standard machine-protection conditions. The improvement in operational efficiency came at a cost:
- the counting rate of the luminometers (or equivalently their visible cross-section) is about an order of magnitude smaller at \\(\\sqrt{s} = 900\\ \\text {GeV}\\) than at \\(13\\ \\text {TeV}\\), reflecting the combination of a 40% drop in the inelastic *pp* cross-section, a factor of two to three drop in the multiplicity of the final-state particles, and a significant softening of their momentum spectrum;
- because of the total-intensity constraints dictated by machine-protection requirements, the injected beam could not exceed 4 bunches per beam (compared to 150 during a routine *vdM*\-calibration session), with a maximum allowed population of \\(1.25\\times 10^{11}~p\\)/bunch;
- the combination of the low beam energy, that suppresses synchrotron-radiation damping, and of the large bunch intensity, that enhances intra-beam scattering, resulted in relatively rapid emittance growth and rather short single-beam lifetimes. To partially mitigate these effects, all collision-rate measurements are expressed in terms of specific luminosity, thereby automatically accounting for the beam-intensity decay.
The bunch intensity was deliberately increased by about 25% beyond the *vdM*\-scan values typical of LHC Run 2, to maximize both the luminosity and the beam-beam parameter \\(\\xi \\). The latter reached 0.010 per IP (Table [1](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Tab1)), compared to 0.003â0.006 during normal *pp* *vdM* sessions. With a target injected emittance of \\(1.5 \\ \\upmu \\text {m} \\cdot \\text {rad}\\), the \\(\\beta \\) function at the IP set to \\(\\beta ^\* = 11\\ \\text {m}\\), and zero crossing angle at IP1 and IP5, the statistical uncertainty affecting a typically 60-second long per-bunch luminosity measurement in the presence of head-on collisions, lay around 0.5%.
**Table 1 Range of beam-beam parameter values during the various stages of the experiment**
[Full size table](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/tables/1)
The bunch patterns were chosen such that all bunches in each beam collided either at IP1 and IP5 only, or at all three of IP1, 2 and 5. The orbit-stabilization feedback system was turned off during the data-taking periods to prevent it from interfering with the beamâbeam-induced orbit shift. The chromaticity was set to its standard value in physics fills, of +10 units, to guarantee coherent stability against the machine impedance, and in view of the very small number of bunches and of their large longitudinal spacing, the beam-stabilizing Landau octupoles were set to their minimum current (1A). The damping time of the bunch-by-bunch transverse feedback was set to 1000 turns, a rather loose setting intended to preserve longer the natural beam oscillation and therefore improve the precision of the tune measurements.
## 4 Beam instrumentation
### 4\.1 Luminometers
The luminometer systems in use by the ATLAS and CMS collaborations at IP1 and IP5 are described in Refs. \[[15](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR15 "ATLAS Collaboration, Improved luminosity determination in pp collisions at
$$\sqrt{s} = 7$$
s
=
7
TeV using the ATLAS detector at the LHC. Eur. Phys. J. C 73, 2518 (2013).
https://doi.org/10.1140/epjc/s10052-013-2518-3
"), [16](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR16 "C.M.S. Collaboration, Development of the cms detector for the cern lhc run 3. J. Instrum. 19(05), 05064 (2024).
https://doi.org/10.1088/1748-0221/19/05/p05064
")\] respectively; no luminosity measurements were available at IP2 during the collider experiment described in this paper. Maximizing the statistical sensitivity leads to choosing the luminosity algorithm with the highest possible acceptance, namely:
- the hit rate per bunch crossing in the ATLAS Minimum Bias Trigger Scintillators (MBTS) \[[15](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR15 "ATLAS Collaboration, Improved luminosity determination in pp collisions at
$$\sqrt{s} = 7$$
s
=
7
TeV using the ATLAS detector at the LHC. Eur. Phys. J. C 73, 2518 (2013).
https://doi.org/10.1140/epjc/s10052-013-2518-3
")\], and
- the occupancy in the CMS Hadron Forward Calorimeter (HFOC) \[[16](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR16 "C.M.S. Collaboration, Development of the cms detector for the cern lhc run 3. J. Instrum. 19(05), 05064 (2024).
https://doi.org/10.1088/1748-0221/19/05/p05064
")\].
In both cases, raw luminometer counts are converted to collision rates using the Poisson formalism \[[15](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR15 "ATLAS Collaboration, Improved luminosity determination in pp collisions at
$$\sqrt{s} = 7$$
s
=
7
TeV using the ATLAS detector at the LHC. Eur. Phys. J. C 73, 2518 (2013).
https://doi.org/10.1140/epjc/s10052-013-2518-3
"), [16](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR16 "C.M.S. Collaboration, Development of the cms detector for the cern lhc run 3. J. Instrum. 19(05), 05064 (2024).
https://doi.org/10.1088/1748-0221/19/05/p05064
")\], and the instantaneous luminosity is averaged over typically 60-second time bins, that during the scans are synchronized with the scan steps. The statistical uncertainties are estimated from either the total number of raw luminometer counts per time bin (ATLAS), or from the RMS of the approximately 40 luminosity samples recorded in a given time bin (CMS). Since the absolute luminosity scale is irrelevant, and in order to simplify the interpretation of the results, all measurements in this paper are presented in terms of fractional shifts in the bunch-averaged specific luminosity, relative to a reference time that depends on the type of measurement considered.
### 4\.2 Bunch-charge measurement
The total intensity of each beam is measured by a direct-current current transformer (DCCT), the absolute scale of which is calibrated against a very high precision pulse generator \[[17](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR17 "C. Barschel, M. Ferro-Luzzi, J.-J. Gras, M. Ludwig, P. Odier, S. Thoulet, Results of the LHC DCCT calibration studies. Technical Report CERN-ATS-Note-2012-026 PERF (2012).
https://cds.cern.ch/record/1425904?ln=en
")\]. The bunch-by-bunch charge fractions, in turn, normalized to the total stored intensity, are measured, separately for the two beams, by Fast Beam Current Transformers (FBCT) \[[18](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR18 "G. Anders et al., Study of the relative LHC bunch populations for luminosity calibration. Technical Report. CERN-ATS-Note-2012-028 PERF (2012).
https://cds.cern.ch/record/1427726
")\]. The latter devices also provide the fill pattern, *i.e.* the relative location, at a given instant and around the two rings, of the nominally filled bunches. Since the absolute scale of the DCCTs is known to much better accuracy than that of the FBCTs, the bunch charges are typically computed as the product of the FBCT bunch-charge fractions and the total circulating beam intensity reported by the corresponding DCCT.
### 4\.3 Emittance measurement
The synchrotron-light beam-profile monitors \[[19](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR19 "G. Trad, Development and Optimisation of the SPS and LHC Beam Diagnostics Based on Synchrotron Radiation Monitors Presented 22 Jan 2015. Presented 22 Jan 2015.
https://cds.cern.ch/record/2266055
")\], dubbed BSRTs (âBeam Synchrotron Radiation Telescopeâ), are mainly used to track emittance evolution over time, but they can also measure the relative changes in transverse single-beam size that result from beamâbeam-induced \\(\\beta \\)\-beating. BSRT data are recorded every second; in what follows, they are presented averaged over one-minute intervals for easier comparison with other measurements.
The absolute length scale of the BSRT profiles suffers from significant uncertainties, if only because the online optical corrections to the images can be updated only a few times per year, and are unable to track the evolution of the efficiency of the light sensors, that depends both on time and on the position of the light spot on the sensor array. For each beam therefore, the BSRT is complemented by a wire-scanner (WS) profile monitor that is located in the same straight section. The advantage of the WS is that its accuracy is significantly better than that of the BSRT; its down side, in the present context, is that the WS cannot acquire data continuously and must be triggered manually. For the results presented in this paper, the wire scanners were flown through the beams at the start of each group of measurements to provide single-beam emittance measurements that are as accurate as possible, and then at regular intervals thereafter. The BSRT is used to interpolate the time evolution of the emittance between two sets of WS measurements.
To interpret measured RMS bunch widths in terms of emittance requires the knowledge of the \\(\\beta \\) functions at the locations of the BSRT and the WS. These are determined by the phase-advance method, and compared in Table [2](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Tab2) to their model value computed using MAD-X; the agreement is typically better than 5%, and the worst disagreement amounts to 9%. These discrepancies are attributed to imperfections in the magnetic lattice. Uncertainties in the measurement mainly arise from the interpolation of the lattice functions between the two closest beam-position monitors (BPMs) where the \\(\\beta \\)\-functions are measured \[[20](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR20 "R. TomĂĄs, O. BrĂŒning, M. Giovannozzi, P. Hagen, M. Lamont, F. Schmidt, G. Vanbavinckhove, M. Aiba, R. Calaga, R. Miyamoto, Cern large hadron collider optics model, measurements, and corrections. Phys. Rev. ST Accel. Beams 13, 121004 (2010).
https://doi.org/10.1103/PhysRevSTAB.13.121004
"), [21](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR21 "R. TomĂĄs, M. Aiba, A. Franchi, U. Iriso, Review of linear optics measurement and correction for charged particle accelerators. Phys. Rev. Accel. Beams 20, 054801 (2017).
https://doi.org/10.1103/PhysRevAccelBeams.20.054801
")\], and either the BSRT or the WS location. Based on the MAD-X LHC lattice model, these uncertainties are estimated not to exceed 3%, and this value is assigned as the systematic uncertainty on the \\(\\beta \\) functions used in BSRT- and WS-based emittance measurements.
**Table 2 \\(\\beta \\) functions at the BSRT and WS locations, for the optical configuration detailed in Sect. [5\.3](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec14)**
[Full size table](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/tables/2)
The beam-averaged emittance, *i.e.* the average of the emittances of the beam-1 bunch and of the corresponding beam-2 bunch, can be obtained directly from emittance scans \[[22](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR22 "O. Karacheban, P. Tsrunchev, on behalf of CMS, Emittance scans for cms luminosity calibration. EPJ Web Conf. 201, 04001 (2019).
https://doi.org/10.1051/epjconf/201920104001
")\] at IP1 and/or IP5. In this approach, the convolved transverse bunch sizes measured using beam-separation scans are translated into emittances using the \\(\\beta ^\*\\) values determined by *k*\-modulation \[[23](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR23 "F. Carlier, R. TomĂĄs, Accuracy and feasibility of the
$${\beta }^{*}$$
ÎČ
â
measurement for lhc and high luminosity lhc using
$$k$$
k
modulation. Phys. Rev. Accel. Beams 20, 011005 (2017).
https://doi.org/10.1103/PhysRevAccelBeams.20.011005
")\] at the relevant IP (see Sect. [4\.3](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec7)). The associated uncertainty (Table [3](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Tab3)) is combined with the statistical uncertainty in convolved transverse width to estimate the error affecting the measured emittance. Comparison with beam-averaged emittances extracted from the single-beam profile monitors reveals excellent agreement between WS and emittance-scan results (Fig. [2](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig2)). The relative time-evolution of the BSRT emittances is consistent with that observed using the WS or emittance scans, but the absolute magnitudes differ significantly, especially in the horizontal plane.
**Table 3 \\(\\beta ^\*\\) functions at the ATLAS and CMS IPs measured by *k*\-modulation. The target value in all cases is \\(\\beta ^\*=\\) 11 m**
[Full size table](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/tables/3)
**Fig. 2**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/2)
Time evolution of the horizontal (top) and vertical (bottom) beam-averaged normalized emittances during LHC fill 8037, as reported by the BSRT (blue), the WS (orange), and using beam-separation scans at IP1 (red) and IP5 (purple). Different shades of the same color correspond to two different bunches present in this fill. For ATLAS emittance scans, only the bunch-averaged value is shown
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/2)
### 4\.4 Beam-beam parameter determination
The systematic and quantitative comparison of the measured and of the predicted beam-beam impact on the observables detailed in Sect. [2](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec2) requires continuous monitoring of the actual beam-beam parameter throughout the duration of the experiment. Only the BSRT provides uninterrupted emittance determination throughout the fill. In view of the scale biases apparent in Fig. [2](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig2), however, WS-based emittances are used, whenever possible, as input to the determination of the beam-beam parameter; BSRT emittances re-scaled to close-in-time, absolute WS measurements provide time-interpolated emittance values whenever WS data are unavailable.
The beam-beam parameter evolution during the first fill of the experiment is illustrated in Fig. [3](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig3). The error bands are dominated by the 3% systematic uncertainty in the measured optical functions. The strength of the beam-beam interaction drops by a factor of two over a couple of hours, from the combined effect of beam-intensity decay and of emittance growth.
**Fig. 3**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/3)
Time evolution of the beam- and plane-averaged beam-beam parameter \\(\\xi \\) inferred from the measured bunch charges and emittances, separately for the two colliding-bunch pairs present in the fill pattern. The color bands indicate the systematic uncertainty (see text)
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/3)
### 4\.5 Tune measurements
Coherent spectra can be monitored using either the LHC Transverse Damper (ADT) \[[24](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR24 "M. SöderĂ©n, J. Komppula, G. Kotzian, S. Rains, D. Valuch, ADT and Obsbox in LHC Run 2, plans for LS2, pp. 165â171 (2019)")\], or the Base-Band Tune (BBQ) \[[25](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR25 "A. Boccardi, M. Gasior, R. Jones, P. Karlsson, R. Steinhagen, First results from the lhc bbq tune and chromaticity systems. Technical report, CERN, Geneva (2009).
https://cds.cern.ch/record/1156349
")\] system; their comparison typically yields consistent results. Spectrograms with a tune resolution of 0.0001 are computed using bunch positions recorded at every turn, and averaged over the one-minute time bins mentioned in Sect. [4\.1](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec5). To mitigate the influence of noise on the tune measurement, a median filter is applied at the pre-processing stage with a self-defined local window size. Additionally, the 50 Hz noise lines present in the spectra are masked. The filtered data are fitted by the sum of a Gaussian function and a constant baseline, and the frequency at the peak of the Gaussian is interpreted as the measured mean tune. The systematic uncertainty on measured tune shifts \\(\\Delta Q\\) is estimated empirically to be \\(\\sigma \_{\\Delta Q}=0.001\\).
### 4\.6 IP-orbit monitoring
Orbit displacements at each IP are measured using the Diode ORbit and OScillation (DOROS) BPM system, that provides sub-micrometer beam-position resolution \[[26](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR26 "M. Gasior, G. Baud, J. Olexa, G. Valentino, First Operational Experience with the LHC Diode ORbit and OScillation (DOROS) System, p. 07 (2017)
https://doi.org/10.18429/JACoW-IBIC2016-MOPG07
")\]. These strip-line BPMs are located in the two quadrupoles on either side of and closest to the IP, allowing the position of both beams to be measured simultaneously. The position and the angle of each beam at the IP are inferred from, respectively, the average and the difference of the positions measured in the two final-triplet quadrupoles.
## 5 Optimization of the phase advance between interaction points
In the two-IP configuration, the basic theory of linearized beam-beam \\(\\beta \\)\-beating provides an analytical determination of the optimum phase advance between the two IPs, *i.e.* of the setting that maximizes the sensitivity, at the witness IP, to beam-beam effects at the scanning IP (Sect. [5\.1](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec12)). A generalization to the three- and four-IP configurations is discussed in Sect. [5\.2](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec13). The required modifications to the baseline LHC optics, as well as their implementation and their experimental validation, are detailed in Sect. [5\.3](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec14).
### 5\.1 Linear beam-beam \\(\\beta \\)\-beating
In circular colliders, the \\(\\beta \\)\-beating induced by beam-beam effects at the IP(s) leads to changes in transverse beam size \[[5](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR5 "W. Herr, T. Pieloni, Beam-Beam Effects. Technical report, CERN (2014). Contribution to the CASâCERN Accelerator School: Advanced Accelerator Physics Course, Trondheim, Norway, 18â29 Aug 2013, p. 29.
https://doi.org/10.5170/CERN-2014-009.431
.
https://cds.cern.ch/record/1982430
"), [27](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR27 "A.W. Chao, Coherent beam-beam effects, in Frontiers of Particle Beams: Intensity Limitations (Springer, Berlin, Heidelberg, 1992), pp. 363â414.
https://doi.org/10.1007/3-540-55250-2_36
")\], thereby impacting the luminosity. For a single collision point, the change in \\(\\beta ^\*\\) of a single small-amplitude particle in (for instance) the horizontal plane *x*, is given by:
\$\$\\begin{aligned} \\frac{\\beta \_x^\*}{\\beta \_{0,x}^\*} = \\frac{1}{\\sqrt{1-4\\pi \\xi \\cot {(2\\pi Q\_x)} - 4\\pi ^2\\xi ^2}}. \\end{aligned}\$\$
(6)
The magnitude of the effect depends on the absolute value \\(\\xi \\) of the beam-beam parameter and on the nominal betatron tune \\(Q\_x\\); its periodicity is \\( \[\\pi \]\\). It can be shown that approaching the half-integer betatron tune from below minimizes the \\(\\beta \\) function at the IP, thereby maximizing the luminosity; this phenomenon has been exploited with great success in (among others) the KEKB and PEP-II *B* factories. In the LHC, where during collisions the nominal fractional betatron tunes are fixed at \\(q\_x=0.31,\\,q\_y=0.32\\), this dynamic-\\(\\beta \\) effect depends only on the beam-beam parameter.
This conclusion no longer strictly holds in the presence of more than one collision point. In the case of two IPs with identical values of \\(\\xi \\) and \\(\\beta ^\*\\), the dynamic-\\(\\beta \\) effect can be described analytically \[[28](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR28 "J. WaĆczyk, A. Dabrowski, D. Stickland, T. Pieloni, W. Kozanecki, X. Buffat, Impact of multiple beam-beam encounters on LHC absolute-luminosity calibrations by the van der Meer method. JACoW IPAC 2023, 053 (2023).
https://doi.org/10.18429/JACoW-IPAC2023-WEPA053
")\]:
\$\$\\begin{aligned} \\frac{\\beta \_x^\*}{\\beta ^\*\_{0,x}} = \\frac{\\sin {2\\pi Q\_x} + 4\\pi \\xi (\\cos (2\\pi Q\_x-2\\mu \_{1,x})-\\cos 2\\pi Q\_x )}{\\pm \\sqrt{1-\\cos ^2{2\\pi (Q\_x+\\Delta Q\_x)}}}, \\end{aligned}\$\$
(7)
where \\(\\mu \_{1,x}\\) is the phase advance from the scanning IP to the witness IP. The sign in the denominator is linked to that of the numerator so as to ensure that the ratio remains positive; the denominator must be positive (resp. negative) below (resp. above) the half-integer tune, *i.e.* when \\(m\< Q\_x \<m+\\frac{1}{2}\\) (resp. \\(m+\\frac{1}{2}\< Q\_x \< m+1\\)), where *m* is a positive integer. The tune shift \\(\\Delta Q\_x\\) is related to \\(Q\_x\\), \\(\\xi \\) and \\(\\mu \_{1,x}\\) by:
\$\$\\begin{aligned} \\cos {2\\pi (Q\_x+\\Delta Q\_x)}&=(1-16\\pi ^2\\xi ^2)\\cos {2\\pi Q\_x} + 8\\pi \\xi \\sin {2\\pi Q\_x} \\nonumber \\\\&\\quad + 16\\pi ^2\\xi ^2\\cos {(2\\pi Q\_x-2\\mu \_{1,x})}. \\end{aligned}\$\$
(8)
The phase-advance dependence of the resulting \\(\\beta \\)\-beating is illustrated in Fig. [4](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig4). The minimum of each curve corresponds to the phase advance that yields the largest mutual beamâbeam-induced luminosity enhancement between the two IPs. The optimal settings are \\(\\mu \_{1,x}^{min}/2\\pi =0.405,\\,\\mu \_{1,y}^{min}/2\\pi =0.410\\); their difference reflects that between the nominal horizontal and vertical fractional tunes \\(q\_x\\) and \\(q\_y\\).
**Fig. 4**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/4)
Analytically computed beam-beam induced \\(\\beta \\)\-beating in a two-IP configuration, as a function of the phase advance \\( \\mu \_1\\) between the 2 IPs (\\(q\_x=0.31\\), \\(q\_y=0.32\\), head-on collisions at both IPs, \\(\\xi =7\\times 10^{-3}\\) per IP)
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/4)
The extrema of the curves in Fig. [4](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig4) can be determined analytically by differentiating Eq. ([7](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Equ7)) with respect to the phase advance \\(\\mu \_{1,x}\\). The general expression for this derivative \\(\\frac{d}{d\\mu \_{1,x}} \\biggl ( \\frac{\\beta \_x^\*}{\\beta ^\*\_{0,x}} \\biggr )\\), detailed in Appendix A, is the product of two factors, each of which can be zero:
- imposing that \\(\\sin {(2\\pi Q\_x - 2\\mu \_{1,x})} = 0\\) yields the solution
\$\$\\begin{aligned} \\frac{\\mu \_{1,x}}{2\\pi } =\\frac{Q\_x}{2}-\\frac{m}{4} \\, \\end{aligned}\$\$
(9)
where *m* is an integer of either sign, or zero;
- the other factor is zero for either
\$\$\\begin{aligned} \\frac{\\mu \_{1,x}}{2\\pi } = \\frac{Q\_x}{2} - A(\\xi , Q\_x)\\, \\end{aligned}\$\$
(10)
or
\$\$\\begin{aligned} \\frac{\\mu \_{1,x}}{2\\pi } = \\frac{Q\_x - 1}{2} + A(\\xi , Q\_x) \\, \\end{aligned}\$\$
(11)
with \\(A(\\xi , Q\_x)\\) defined in Eq. ([12](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Equ12))
\$\$\\begin{aligned} A(\\xi , Q\_x) = \\arccos \\left( \\frac{ 32\\pi ^3\\xi ^3 \\sin (4\\pi Q\_x) + 48\\pi ^2\\xi ^2 \\cos ^2(2\\pi Q\_x) - 32\\pi ^2\\xi ^2 - 6\\pi \\xi \\sin (4\\pi Q\_x) + \\cos ^2(2\\pi Q\_x) - 1}{16\\pi ^2\\xi ^2 \\left( 4\\pi \\xi \\sin (2\\pi Q\_x) + \\cos (2\\pi Q\_x) \\right) } \\right) / (4\\pi ) \\end{aligned}\$\$
(12)
The two sets of solutions do not overlap. The subset of solutions that minimize (rather than maximize) \\(\\beta ^\*\\) is identified by requiring that the second derivative be positive. This restricts the main solution to:
\$\$\\begin{aligned} \\frac{\\mu \_{1,x}}{2\\pi } =\\frac{Q\_x}{2}+\\frac{m+1}{4} . \\end{aligned}\$\$
(13)
The optimal phase \\(\\mu \_{x,1}\\) that minimizes the \\(\\beta ^\*\\) for range of \\((\\xi , Q\_x)\\) values is shown in Fig. [5](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig5).
**Fig. 5**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/5)
Optimal phase advance value \\(\\mu \_{x,1}\\) for various values of \\((\\xi , Q\_x)\\) in the two-IPs configuration
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/5)
The phase-advance dependence of the dynamic-\\(\\beta \\) effect can also be quantified in terms of the luminosity shift, at the witness IP, that is associated with the electromagnetic interaction of the two bunches at the scanning IP (Fig. [6](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig6)). The optimal setting, *i.e.* that which maximizes the sensitivity, at the witness IP, to beam-beam effects at the scanning IP, is indicated by the green dotted vertical line at \\(\\Delta \\mu \_{\\text {IP1-IP5}}=0.41\\ \[2\\pi \]\\); the full suppression of the phase-related luminosity enhancement is indicated by the red dotted line. These results are consistent with the predictions of the analytical model (Fig. [4](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig4)). In the two-IP configuration considered here, in which IP1 and IP5 can both play the role of either the scanning or the witness IP, and due to the periodicity of this curve (\\( \\pi \\)) and to the fractional-tune values used during collisions (\\(q\_x / q\_y = 0.31 / 0.32)\\), the phaseâadvance-related luminosity enhancement is the same at the two IPs; at the optimum setting, it is three times larger than it would be using the nominal LHC lattice. The maximum effect depends mainly on the transverse tunes, and thus could have been further enhanced if it had been possible to move these away from their nominal values.
A more comprehensive simulation study, that evaluates the impact of the beam-beam interaction on the luminosity for more realistic Gaussian proton-density distributions, multiple collision points and as a function of the beam separation, is reported in Ref. \[[29](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR29 "J. WaĆczyk, Precision luminosity measurement at hadron colliders. PhD Thesis, Ăcole Polytechnique FĂ©dĂ©rale de Lausanne (2024).
https://doi.org/10.5075/epfl-thesis-10500
")\].
**Fig. 6**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/6)
Beamâbeam-induced luminosity shift at IP5 predicted by COMBI simulations, as a function of the phase advance \\( \\mu \_1=\\Delta \\mu \_{\\text {IP1-IP5}}\\) between the two IPs (\\(q\_x=0.31\\), \\(q\_y=0.32\\), head-on collisions at both IPs, \\(\\xi =7\\times 10^{-3}\\) per IP)
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/6)
### 5\.2 Optimized phases for multi-IP configurations
The general expression to obtain the first-order change of the \\(\\beta \\) function at an arbitrarily chosen reference IP (at \\(s=\\mu \_{0,x}=0\\)) can be approximated as multiple quadrupole errors \[[30](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR30 "T. Pieloni, J. Barranco, X. Buffat, P.C. Jorge, L.M. Medrano, C. Tambasco, R. TomĂĄs, Dynamic beta and beta-beating effects in the presence of the beam-beam interactions in 57th ICFA Advanced Beam Dynamics Workshop on High-Intensity and High-Brightness Hadron Beams, p. 027 (2016).
https://doi.org/10.18429/JACoW-HB2016-MOPR027
")\]:
\$\$\\begin{aligned} \\frac{\\beta ^\*\_x}{\\beta ^\*\_{0,x}} = \\frac{2\\pi \\xi }{\\sin {2\\pi Q\_x}} \\sum \_{i\\in IPs}^{N-1} \\cos (2\\pi Q\_x-2\\mu \_{i,x})\\, \\end{aligned}\$\$
(14)
where \\(\\mu \_{i,x}\\) is the horizontal phase advance between the reference and the \\(i^{th}\\) IP for a total of *N* IPs. Therefore, the optimal configuration, which minimizes the \\(\\beta ^\*\\) change at all IPs, must satisfy the condition:
\$\$\\begin{aligned} \\sum \_{i=1}^{N-1} \\cos (2\\pi Q\_x-2\\mu \_{i,x}) = -(N-1)\\, \\end{aligned}\$\$
(15)
which yields the same solution for all IPs as in the 2 IPs configuration (Eq. [13](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Equ13)):
\$\$\\begin{aligned} \\frac{\\mu \_{i,x}}{2\\pi } =\\frac{Q\_x}{2}+\\frac{n+1}{4} . \\end{aligned}\$\$
(16)
where *n* is an integer of either sign, or zero. To maximize the effect on luminosity, the phase advance between the reference interaction point at \\( i = 0 \\) and the nearest IP, at \\( i = 1 \\) must satisfy \\( \\mu \_{1,x}/2\\pi = (Q\_x + 0.5)/2 \\), as given by the first solution (for \\(n=0\\)) of Eq. ([13](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Equ13)). This condition implies that the phase advance to any subsequent adjacent IPs (\\( i \\ge 2 \\)) must be an integer multiple of \\( \\pi \\). Thus, this requirement exposes a fundamental constraint: it is not possible to minimize \\( \\beta ^\* \\) simultaneously at all IPs. This would only be possible for (\\( Q\_x = 0.5 \\)) which is practically not achievable, as it would require operating the storage ring precisely at the half-integer resonance.
### 5\.3 Optics validation
The optimum phase advance from IP1 to IP5 determined in Sect. [5\.1](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec12) was implemented using tuning trim quadrupoles in the LHC arcs. The targeted fractional phase advance was \\(0.9 \\, \[2 \\pi \]\\) for both beams and both planes (Table [4](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Tab4)). This setting, which differs from that in Fig. [6](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig6) by one period (\\( \\pi \\)), is equivalent in terms of sensitivity optimization but has the advantage that it minimizes the absolute magnitude of the perturbations inflicted upon the nominal lattice.
Several measurements were carried out to validate this adjustment: comparison of the targeted and achieved shifts in phase advance (Sect. [5\.3.1](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec15)), \\(\\beta \\)\-function measurements across the full LHC circumference (Sect. [5\.3.2](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec16)) as well as at the IPs (Sect. [5\.3.3](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec17)), and dispersion measurements to estimate the dispersive contribution to the transverse beam size at the IP (Sect. [5\.3.4](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec18)).
#### 5\.3.1 Phase advance
The phase advance around each LHC ring is determined using turn-by-turn orbit measurements \[[20](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR20 "R. TomĂĄs, O. BrĂŒning, M. Giovannozzi, P. Hagen, M. Lamont, F. Schmidt, G. Vanbavinckhove, M. Aiba, R. Calaga, R. Miyamoto, Cern large hadron collider optics model, measurements, and corrections. Phys. Rev. ST Accel. Beams 13, 121004 (2010).
https://doi.org/10.1103/PhysRevSTAB.13.121004
"), [21](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR21 "R. TomĂĄs, M. Aiba, A. Franchi, U. Iriso, Review of linear optics measurement and correction for charged particle accelerators. Phys. Rev. Accel. Beams 20, 054801 (2017).
https://doi.org/10.1103/PhysRevAccelBeams.20.054801
")\]. The targeted and achieved phase-advance shifts agree within \\(0.01\\, \[2\\pi \]\\) or better (Table [4](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Tab4)). The phase advance change from IP1 and IP5 was compensated on the opposite side of the ring, ensuring that the transverse tunes remain unchanged. The target fractional local phase advance is set to \\(0.9,\[2\\pi \]\\). Table [4](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Tab4) lists the corresponding phase shifts applied to the nominal lattice to reach this target on one side of the ring; this is consistent with the analytically derived value of \\(0.41\\,\[2\\pi \]\\) for the compensating shift on the opposite side. The change in phase advance was also measured at the BSRT location, where the discrepancy between the model and measurements is more pronounced, reaching up to \\(0.05\\,\[2\\pi \]\\).
**Table 4 Nominal and optimized IP1 to IP5 phase-advance settings. Column 4 (resp. 5) displays the targeted (resp. measured) shift in phase advance after the optimized settings have been applied**
[Full size table](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/tables/4)
#### 5\.3.2 Residual \\(\\beta \\)\-beating
The \\(\\beta \\)\-beating was quantified and found to be well within acceptable tolerances for various configurations of the trim strength. The measurement results, relative to the reference measurements without any phase advance adjustments, are illustrated in Fig. [7](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig7). The maximal observed deviation is approximately 25% around IP5 for Beam 1 in the vertical plane. Such significant deviations in proximity to the IPs are anticipated and are frequently attributed to erroneous BPM readings. Overall, the measurements indicated that \\(\\beta \\)\-beating remained within acceptable bounds, predominantly within a 10% margin.
**Fig. 7**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/7)
Measured \\(\\beta \\)\-function change from the full strength of the phase advance knob with reference to nominal lattice along the LHC ring, for Beam 1 (top two figures for *x* and *y*) and Beam 2 (bottom two figures for *x* and *y*)
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/7)
The measured \\(\\beta \\) functions are also compared to the MADX LHC lattice model predictions, which include the phase advance adjustment. The results, illustrated along the LHC lattice in Fig. [8](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig8), indicate that the measurements align with the model within an acceptable margin of 10%.
**Fig. 8**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/8)
Measured deviations in the \\(\\beta \\)\-function along the LHC ring relative to the MADX model incorporating the phase advance knob, for Beam 1 (top two figures for *x* and *y*) and Beam 2 (bottom two figures for *x* and *y*). The error bars represent solely the statistical uncertainties as determined in the measurements
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/8)
#### 5\.3.3 \\(\\beta \\) functions at the interaction points
The DOROS BPMs are strategically located at the end of the inner triplets on both sides of each of the experiments. Due to the phase advance of \\(\\pi \\) between the DOROS on the left and right sides, these are inadequate for the precise determination of the measurement of \\(\\beta ^\*\\) at the IP. This can only be achieved with a special *k*\-modulation technique \[[31](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR31 "F. Carlier, et al. Challenges of K-modulation measurements in the LHC Run 3, in Proceedings of IPACâ23. IPACâ23â14th International Particle Accelerator Conference (JACoW Publishing, Geneva, Venice, 2023), pp. 531â534.
https://doi.org/10.18429/JACoW-IPAC2023-MOPL014
")\], which was performed at the interaction points used in the experiment. The results are summarized in Table [3](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Tab3). While the achieved precision of the \\(\\beta ^\*\\) determination varies between measurements, uncertainties at the level of \\({\\mathcal {O}}(1\\%)\\) are not uncommon and are not systematic to a specific beam or interaction point. Such uncertainties would propagate to only a few-percent effect on derived quantities such as the emittance or the beamâbeam parameter.
#### 5\.3.4 Dispersion
It was essential to perform dispersion \\(D^\*\\) measurements at the interaction points, under the new optics configurations, as it can contribute to the measured beam size. Assuming the LHC design momentum spread value of \\(\\frac{\\Delta p}{p}=10^{-4}\\), the average dispersion measurement around IP1 and IP5 is estimated to contribute up to 2% to the observed bunch sizes (worst case measured for B2Y). This effect is considered static throughout the experiment.
## 6 Experimental results
This section presents experimental measurements of beamâbeam effects in configurations involving multiple collision points. The impact of the optics modifications introduced in Sect. [5](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec11) is first assessed through tune spectra measurements, as detailed in Sect. [6\.1](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec20). To quantify the cumulative effect of one or more additional collisions and assess the reproducibility of the results, a stepwise collision scan was employed. In this procedure, beams were sequentially brought into head-on collision at selected interaction points (IPs), followed by full transverse separation (\\(\\Delta \> 6\\sigma \\)) in either the horizontal or vertical plane at all but one witness IP. Two distinct filling schemes were used to isolate the influence of a single additional collision (2-IP configuration: ATLAS and CMS, see Sect. [6\.2](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec21)) and to investigate more complex multi-collision scenarios (3-IP configuration: including ALICE, see Sect. [6\.3](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec25)). The first measurement of separation-dependent beamâbeam effects, analogous to those encountered during luminosity calibration scans, is presented in Sect. [6\.4](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec28), alongside comparisons with numerical predictions. An additional measurement conducted during the 2022 van der Meer calibration is detailed in Sect. [6\.5](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec32). A consolidated overview of all findings is provided in Sect. [7](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec33).
### 6\.1 Measured effect on the tune spectra of the phase advance changes
The coherent transverse tune spectra serve as a diagnostic for evaluating the effectiveness of the phase advance modifications. COMBI simulations were carried out for two extreme lattice configurations, one designed to maximize and the other to suppress beamâbeam effects. In the head-on scenario (Fig. [9](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig9)), the coherent modes are suppressed for the maximizing configuration, as the phase advance between the two collision points is close to \\(\\pi \\), as was previously shown in simulation in Ref. \[[32](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR32 "Y. Alexahin, A study of the coherent beam-beam effect in the framework of the Vlasov perturbation theory. Nucl. Instrum. Methods Phys. Res., Sect. A 480(1â3), 253â288 (2002).
https://doi.org/10.1016/S0168-9002(01)01219-0
")\]. In this configuration, the effect of the first collision is counteracted by the second. This setting is particularly suitable for studying incoherent spectra, as the diminished coherent signal allows for a clearer extraction of the mean tune shift from fitted spectra. During separation scans, the symmetry between horizontal (*x*) and vertical (*y*) planes is broken. In the scanning plane (horizontal), coherent modes re-emerge at \\(1.5\\,\\sigma \\) separation (Fig. [10](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig10)), with their frequency shifting in response to the separation distance (additional example for \\(3\\,\\sigma \\) is shown in Fig. [11](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig11)).
**Fig. 9**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/9)
Comparison of COMBI simulated tune spectra in *x*\-plane (left) and *y*\-plane (right) for the suppressing (in blue) and maximizing (in green) phase configurations, for head-on collision
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/9)
**Fig. 10**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/10)
Comparison of COMBI simulated tune spectra in *x*\-plane (left) and *y*\-plane (right) for the suppressing (in blue) and maximizing (in green) phase configurations, at \\(1.5\\,\\sigma \\) separation in the horizontal plane
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/10)
**Fig. 11**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/11)
Comparison of COMBI simulated tune spectra in *x*\-plane (left) and *y*\-plane (right) for the suppressing (in blue) and maximizing (in green) phase configurations, at \\(3\\,\\sigma \\) separation in the horizontal plane
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/11)
These simulated trends of the coherent spectra are also observed in experimental data. Transverse spectra recorded for head-on conditions confirm the suppression of coherent modes in the maximizing phase configuration (Fig. [12](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig12)). Additional comparisons for separated beams at the IP (see Figs. [13](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig13) and [14](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig14)), shown here for the case of vertical separation, reveal similar features. In this configuration, the non-scanning (horizontal) plane exhibits spectra that qualitatively agree with simulation predictions (see Figs. [10](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig10) and [11](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig11)), while in the scanning (vertical) plane the coherent modes do not reappear as expected. An equivalent behavior is observed when the separation is applied in the horizontal plane, with the roles of the two transverse planes interchanged; this case is not shown for brevity. This discrepancy is attributed to the effect of the transverse damper, which acts independently on the two beam centroids and suppresses coherent oscillations arising from beamâbeam coupling, consistent with earlier observations \[[33](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR33 "R. Calaga, W. Herr, G. Papotti, T. Pieloni, X. Buffat, S. White, R. Giachino, Coherent beam-beam mode in the LHC. CERN (2014).
https://doi.org/10.5170/CERN-2014-004.227
.
http://cds.cern.ch/record/1957039
")\].
**Fig. 12**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/12)
Measured tune spectra before and after phase change, corresponding to the suppressing (in blue) and maximizing (in green) phase configurations, at the head-on position, for each of the transverse planes (Beam 1 shown as an example). Dash-dotted lines show Gaussian fits to the data
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/12)
**Fig. 13**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/13)
Comparison of the measured tune spectra at suppressing (blue) and maximizing (green) phase configurations with \\(1.3\\,\\sigma \\) separation in the vertical plane
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/13)
**Fig. 14**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/14)
Comparison of the measured tune spectra at suppressing (blue) and maximizing (green) phase configurations with \\(2.6\\,\\sigma \\) separation in the vertical plane
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/14)
### 6\.2 Head-on collisions in a two-IP configuration at \\(\\sqrt{s} = 900\\ \\text {GeV}\\)
Each test consisted of four measurements designed to quantify the beamâbeam effect induced by a single additional head-on collision. The sequence was executed with CMS as the witness IP, followed by ATLAS, and concluded with a repeat at CMS. Over the course of the fill, the beamâbeam parameter decreased from approximately 0.01 to 0.005. The outcome of each scan is analyzed using three independent observables: luminosity, beam size, and tune spectra, discussed separately below.
#### 6\.2.1 Impact of beam-beam effects on head-on luminosity
An overview of the luminosity measurements during the first test with CMS as the witness IP is shown in Fig. [15](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig15), where the specific luminosity at CMS is normalized to the initial head-on step to highlight relative changes. The uncertainty at each step is given by the standard error of the mean; shaded bands represent the standard deviation. A comparable dataset for ATLAS is presented in Fig. [16](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig16), confirming the expected luminosity change due to beamâbeam effects originating from another IP. Relative luminosity changes for each step in all tests are summarized in Table [5](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Tab5). These values are computed relative to the average of adjacent head-on points, which removes the linear luminosity decay over time. Assuming slow change of beam parameters, the results within each test are expected to be consistent.
**Fig. 15**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/15)
CMS specific luminosity, normalized to the first measured both IPs head-on point (blue), during the first step-function scan at IP1. The measurements with ATLAS fully separated are shown in red. The error bars indicate error on the mean
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/15)
**Fig. 16**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/16)
ATLAS specific luminosity, normalized to the first measured both IPs head-on point (blue), during the step-function scan at IP5. The points for which IP5 was separated are shown in red. The error bars indicate the statistical error
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/16)
**Table 5 Measured luminosity shifts in each of full separation steps â summary of all three step-function scans performed for IP1 and IP5 as the witness IP. The beam-beam parameter at the start of each scan \\(\\xi \_\\text {start}\\) is included for reference**
[Full size table](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/tables/5)
All measurements are compiled in Fig. [17](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig17) and compared to COMBI simulation predictions. The results are plotted as a function of the beamâbeam parameter (estimated as described in Sect. [4\.4](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec8)), following a reverse chronological order. Agreement with simulation is excellent across the entire range.
**Fig. 17**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/17)
COMBI simulated luminosity enhancement induced by a single head-on collision as a function of the beam-beam parameter (dashed line), compared to the test results at both ATLAS (red points) and CMS (green points)
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/17)
**Fig. 18**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/18)
Measurements of beam-beam effects on beam widths shown with red points for each beam and plane, during the experiment with IP5 as the witness IP at \\(\\xi =0.0056\\). COMBI predictions are shown in blue with calculated systematic error from the phase advance error between the IP and BSRT locations
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/18)
#### 6\.2.2 Impact of beam-beam effects on single-beam sizes
Complementary measurements of the bunch size from the BSRT are shown in Fig. [18](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig18). Data from all steps and bunches are combined, and a linear trend corresponding to continuous emittance growth is subtracted using the difference between the first and last head-on steps. For comparison, COMBI simulation results are also included in the figure. The dominant source of uncertainty in COMBI predictions arises from discrepancies between the measured phase advance and the MADX model, with the largest observed deviation being \\(0.05\\,\[2\\pi \]\\) (as discussed in Sect. [5\.3.1](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Sec15)). This uncertainty is propagated using a modified form of Eq. ([6](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Equ6)) adapted to the BSRT location. The magnitude of this uncertainty depends on the specific phase advance configuration, as some are more sensitive to small deviations. This is reflected in the varying sizes of the error bands associated with each prediction. The most accurate predictions are obtained in the vertical plane of Beam 2, where the beamâbeamâinduced \\(\\beta \\)\-beating propagates to the BSRT with a particularly favorable phase advance. In this configuration, the measured effect is also the most pronounced, and excellent agreement with COMBI simulations is observed. In other cases, the expected absolute effect is smaller, typically below 1%, and the influence of phase advance uncertainty becomes more significant, limiting the predictive reliability of the simulations. Experimentally, it is also possible that in these configurations the beamâbeamâinduced effect is too subtle to resolve, potentially masked by instrumental limitations.
#### 6\.2.3 Impact of beam-beam effects on tune spectra
Figure [19](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig19) shows representative ADT spectra recorded during the step-function scans. These spectra are affected by coupling, noise, and various nonlinearities, so they do not follow an idealized analytical distribution. Nonetheless, the peak position can be extracted using a Gaussian fit (indicated by dashed lines), yielding an estimate of the mean tune shift. In the single-collision configuration (blue lines), unsuppressed coherent modes dominate the spectrum and partially obscure the underlying incoherent component. The Gaussian fit is therefore applied to the incoherent contribution between coherent lines, rather than to the global spectral maximum, which results in a larger uncertainty on the extracted peak position. Despite this, in each of the performed step-scans at both collision points, the single-collision tune shift, defined as twice the fitted mean shift, is found to be consistent within \\(\\pm 5\\%\\) with the expected coherent beamâbeam tune shift corresponding to the conditions of the respective scan. For head-on collisions with Gaussian transverse beam distributions, the coherent tune shift of the dipole mode is expected to be \\(Y\\xi \\), where *Y* is the Yokoya factor (\\(Y = 1.21\\) for round beams) \[[34](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR34 "K. Yokoya, Y. Funakoshi, E. Kikutani, H. Koiso, J. Urakawa, Tune shift of coherent beam-beam oscillations. Part. Accel. 27, 181â186 (1990)")\]. As expected, the two-collision configuration yields a tune shift approximately double that of the single-collision case \[[13](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR13 "A. Babaev, T. Barklow, O. Karacheban, W. Kozanecki, I. Kralik, A. Mehta, G. Pasztor, T. Pieloni, D. Stickland, C. Tambasco, R. Tomas, J. WaĆczyk, Impact of beam-beam effects on absolute luminosity calibrations at the CERN large hadron collider. Eur. Phys. J. C 84, 17 (2024).
https://doi.org/10.48550/arXiv.2306.10394
")\].
**Fig. 19**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/19)
Measured vertical tune spectra for different collision configurations: reference nominal tune (red), without the presence of any collisions, compared to single collision (blue) and two collisions (green) at \\(\\xi \\simeq 0.0086\\). The dash-dotted lines in corresponding colors indicate Gaussian fits to the data used to extract mean tunes
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/19)
**Fig. 20**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/20)
Specific luminosity, normalized to the first head-on point, measured at ATLAS (top) and CMS (bottom) during beam separation at the other IPs, with horizontal separation applied in steps 2â4 and vertical separation in steps 6â8. The tested configurations are indicated in the legend. The blue line shows an interpolation of the head-on measurements, serving as a reference for relative luminosity changes
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/20)
### 6\.3 Head-on collisions in a three-IP configuration at \\(\\sqrt{s} = 900\\ \\text {GeV}\\)
An additional head-on collision was included in the configuration using IP2 (ALICE), as it was the only IP where this was feasible. A specially designed filling scheme ensured uniform collision exposure across all bunches, allowing measurements to be combined and thereby reducing statistical uncertainties. However, no reliable luminosity data were available from IP2 during this experiment.
**Table 6 Luminosity shift per full separation step, calculated with respect to an interpolated reference obtained from the two closest head-on points. Results for horizontal and vertical separations are shown as Shift part 1 and Shift part 2, respectively**
[Full size table](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/tables/6)
#### 6\.3.1 Impact of beam-beam effects on head-on luminosity
A step-function scan was conducted by sequentially collapsing or separating one IP at a time, with the separation applied first in the horizontal plane and, in a second stage of the experiment, in the vertical plane. A comprehensive summary of the measured specific luminosity at ATLAS and CMS is presented in Fig. [20](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig20). Exceptionally stable beam conditions during this test enabled measurements at both experiments under nearly identical beamâbeam parameters. The beamâbeam-induced luminosity changes observed during each step are quantified in Table [6](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Tab6), referenced to the time-interpolated specific luminosity from three head-on points.
Orbit-drift corrections are applied only to these head-on reference points and result in changes of the extracted luminosity shifts of at most 0.2%. Good agreement with COMBI simulations is observed across the separation steps, with measured shifts closely matching the simulated values for configurations involving IP1 and IP5. For single-IP separation, the level of agreement depends on the witness IP and the separation step. When IP1 is used as the witness, the luminosity shift in the first separation step is reproduced almost exactly by COMBI, while a smaller shift is observed in the second step. This behaviour is expected, as the beamâbeam parameter decreases during the scan and the measurement becomes increasingly sensitive to residual orbit perturbations. For configurations involving IP2, the comparison is less conclusive. As no reliable luminosity data were available from IP2 during this experiment, it was not possible to precisely optimize the collision conditions at this IP. As a result, the actual beamâbeam conditions at IP2 may differ from those assumed in the simulations, affecting the comparison. Residual orbit drifts affecting the separation steps cannot be fully corrected. As illustrated by the DOROS measurements in Fig. [21](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig21), the orbit data contain both slow drifts and beamâbeamâinduced orbit offsets arising from changes in the collision configuration when other IPs are separated. While the beamâbeamâinduced offsets are included in the COMBI simulations, the residual orbit drifts are not. Since these contributions cannot be disentangled in the measurements, a full correction is not possible. Consequently, the effective beamâbeam impact measured during the experiment may be weaker than that assumed in the simulations.
**Fig. 21**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/21)
Residual separation between the two beams at all three IPs in horizontal (top) and vertical (bottom) planes, for test with observations at ATLAS, where the separation was performed first in the horizontal and then in the vertical plane at the other IPs
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/21)
#### 6\.3.2 Tune shift scaling with the number of beamâbeam interactions
Analysis of the transverse tune spectra in this configuration is further complicated by the reappearance of coherent modes when beams are separated at one or more IPs, as predicted by simulation. In such cases, the maximum tune shift is determined from the frequency difference between the \\(\\pi \\) and \\(\\sigma \\) modes, rather than from fitting the incoherent peak. For configurations where all three IPs were in collision, coherent modes remained suppressed due to the retained symmetry between IP1 and IP5. Figure [22](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig22) presents the measured tune shifts for all tested collision configurations. The unperturbed tune measurement, shown at â0â collisions, serves as the reference from which all tune shifts are calculated. An empirical uncertainty of 0.001 is assumed for each independent measurement. Both transverse planes and separation directions are included to demonstrate consistency. As expected, the measured tune shift scales linearly with the number of head-on collisions.
**Fig. 22**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/22)
Measured mean tune shift as a function of the number of head-on collisions, shown for both transverse planes and for both separation directions (horizontal and vertical) used in the experiment. The separation refers to how collisions were removed: it was applied in only one transverse direction at a time (H sep. or V sep.) at the non-colliding interaction points, which are not explicitly labeled on the x-axis
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/22)
### 6\.4 Beam-beam signatures during transverse-separation scans at \\(\\sqrt{s} = 900\\ \\text {GeV}\\)
In the van der Meer (vdM) calibration method, beamâbeam corrections are applied separately at each separation step. Therefore, it is essential to validate the simulated dependence of luminosity shifts on nominal beam separation. Full vdM-like scans were performed in both transverse planes at ATLAS and CMS, using a configuration where only these two IPs were in collision.
#### 6\.4.1 Impact of beamâbeam effects on the luminosity scan-curves
Figure [23](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig23) shows the luminosity curves measured at ATLAS, acting as the witness IP, during a separation scan performed at CMS. The curves are fitted with Gaussian profiles for visual guidance. The beam separation at CMS is expressed in units of the measured transverse beam size, \\(\\sigma \_{\\textrm{meas}}\\), which was found to be smaller than the nominal beam size, \\(\\sigma \_{\\textrm{nom}}\\), used to define the scan steps. As a result, scan points corresponding to nominal separations beyond \\(6\\,\\sigma \_{\\textrm{nom}}\\) extend to larger values when expressed in units of \\(\\sigma \_{\\textrm{meas}}\\).
To account for emittance evolution during the scan, both linear and exponential fits are applied to the head-on luminosity points before, after, and in the middle of the scan. These trends are subtracted from the data in Fig. [23](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig23), and all measurements are presented relative to this fitted head-on reference (normalized to unity).
**Fig. 23**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/23)
ATLAS specific luminosity relative changes (blue points) during separation scans in *x*\-plane (top) and *y*\-plane (bottom) at CMS, corrected for exponential (exp. corr.) luminosity decay. Gaussian function with a constant fitted to these points is shown with blue dashed line and compared to fit performed on the data corrected for linear luminosity change (lin. corr.) in time (dotted line). These are compared to COMBI simulation predictions (purple solid line)
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/23)
At CMS, similar scans were performed; however, degraded beam stability was observed during this period of the fill. The scan results exhibit systematic deviations, especially at large separation, which could not be symmetrically corrected using simple linear or exponential models. Despite these issues, all four scan cases (two transverse planes at two IPs) yield a total specific luminosity reduction of approximately 2.4%, consistent with expectations based on the measured beamâbeam parameter \\(\\xi \\). Comparison with COMBI simulations reveal up to 30% discrepancy for IP5 as the witness IP, while agreement improves to within 10% for measurements at IP1. These differences are comparable to those observed in Table [6](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Tab6), suggesting a common underlying origin. While a finer scan granularity can improve the description of the steep-response region below \\(3\\,\\sigma \_{\\textrm{meas}}\\), it would not address the dominant discrepancy between head-on and fully separated luminosity, which is likely driven by simulation assumptions and scan-related systematics.
#### 6\.4.2 Beamâbeam-induced orbit deflections
Beamâbeam deflections are clearly observed during separation scans. Figure [24](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig24) shows the fully corrected residual horizontal beam displacement at IP1, measured with the DOROS beam position monitors, as a function of the beam separation. This representation allows a direct comparison with the analytical beamâbeam deflection prediction based on the BassettiâErskine formalism \[[35](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR35 "M. Bassetti, G.A. Erskine, Closed expression for the electrical field of a two-dimensional gaussian charge. Technical report, CERN, Geneva (1980).
https://cds.cern.ch/record/122227
")\].
The raw orbit data were corrected for slow orbit drifts by aligning the beam positions at the luminosity optimization points before and after each scan. Deviations of the applied beam separation from the nominal settings were accounted for by fitting a common length scale simultaneously with the beamâbeam deflection. This approach avoids an overestimation of the separation calibration. The measured orbit displacement at the DOROS location was corrected for the contribution arising from the beamâbeam deflection angle by applying a geometric amplification factor, proportional to the distance between the IP and the BPM in the small-angle approximation. After all corrections, the measured beamâbeam deflection as a function of beam separation closely follows the analytical prediction. The maximum horizontal deflection observed at IP1 is 10.2 \\(\\upmu \\)m, compared to an expected value of 8.4 \\(\\upmu \\)m, with agreement within the total measurement uncertainty. The dominant uncertainty contribution (3.5 \\(\\upmu \\)m, corresponding to about 90% of the total) arises from a possible residual half crossing angle at the IP of up to 0.5 \\(\\upmu \\)rad. Additional systematic uncertainties include the fitted separation length scale, tune uncertainty of 0.005, and an extended 2 % uncertainty on \\(\\beta ^\*\\). Statistical uncertainties are negligible for most separation steps.
**Fig. 24**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/24)
Orbit offset caused by the beamâbeam interaction measured by DOROS BPMs in the horizontal plane at IP1 during the horizontal separation scan. Analytical prediction is shown for comparison \[[35](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR35 "M. Bassetti, G.A. Erskine, Closed expression for the electrical field of a two-dimensional gaussian charge. Technical report, CERN, Geneva (1980).
https://cds.cern.ch/record/122227
")\]
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/24)
#### 6\.4.3 Impact of the beam-beam effects on beam sizes and tune shift
Separation-dependent changes in beam size were measured using the synchrotron light monitor. Figure [25](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig25) shows good consistency across individual bunch measurements, as well as strong agreement with COMBI simulations. This reinforces the utility of the BSRT as a reliable diagnostic for beamâbeam effects. The uncertainties shown are statistical only and are treated as uncorrelated between scan points. The normalization to the zero-separation value is applied for visualization purposes and does not impose a constraint on that point.
**Fig. 25**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/25)
BSRT beam-size measurement in the vertical plane for each bunch (points) during horizontal (top) and vertical (bottom) separation scans at IP1. The exponential evolution correction is applied based on adjacent head-on points (before, in the middle and after the scan) to highlight the relative changes. COMBI predictions are shown (gray dashed line) for comparison
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/25)
Transverse tune shifts were extracted from ADT-measured spectra. As illustrated in Fig. [26](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig26), the mean tune shift evolution is reconstructed for both the scanning and non-scanning planes. The measured tune shifts closely follow the expected COMBI trends, further validating the simulation models used for beamâbeam corrections.
**Fig. 26**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/26)
Measured tune shift in units of the beam-beam parameter \\(\\xi \_{BB}\\) for the horizontal (blue) and vertical (orange) planes, during a horizontal separation scan at IP1. The error bars indicate typically assumed conservative empirical systematic error on the tune measurement of 0.001. The solid curves represent the COMBI prediction
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/26)
### 6\.5 Beam-beam signatures during vdM scans at \\(\\sqrt{s} = 13.6\\ \\text {TeV}\\)
During the first vdM calibration of LHC Run 3, an unusually high beamâbeam parameter was recorded, offering a unique opportunity to study beamâbeam effects under enhanced interaction conditions. These effects are examined using online ATLAS luminosity data, with representative results shown in Fig. [27](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#Fig27).
The example focuses on a standard vdM scan pair conducted at CMS, while monitoring luminosity at the non-scanning IP (ATLAS). The plotted luminosity changes are normalized to the central scan point, where both IPs are in head-on collision. Variations in this normalized luminosity across the scan directly reflect the impact of beamâbeam interaction at the scanning IP on the luminosity at the non-scanning IP. These patterns are reproducible across multiple scan pairs acquired within the same fill. At the largest nominal separations (approximately \\(\\pm 0.6\\,\\)mm), the beams are sufficiently separated to eliminate beamâbeam effects entirely. Notably, a qualitative difference is observed between the horizontal and vertical scan planes, an asymmetry that is not present in the above presented dedicated experiment. This behavior is attributed to differences in the phase advance between IPs in the standard vdM optics configuration. For vertical (*y*) scans, experimental data agree well with simulation, confirming the robustness of the modeling under these conditions. In contrast, the horizontal (*x*) scans exhibit a systematic discrepancy: the observed beamâbeam impact is broader than that predicted by simulation. The model assumes that the separation steps are accurately defined using the beam size measured during the vdM scans, but this mismatch may suggest an unmodeled phase advance deviation or another optics imperfection. The total difference between data and simulation is within 0.2%. This dataset originates from a standard vdM calibration, for which only online luminosity measurement was available, and many contributions to the overall vdM uncertainty, such as bunch-to-bunch variations, averaging over different bunch families, and emittance evolution during the scan, are not disentangled. Therefore, the 0.2% difference should be interpreted as a consistency check under typical operational conditions rather than a precision validation of the beamâbeam model. Nevertheless, this indicates that the model remains sufficiently accurate for beamâbeam corrections in vdM analyses.
**Fig. 27**
[](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/27)
Changes from the separation-dependent beam-beam effects at the non-scanning IP â ATLAS online luminosity shown with points. Data comes from a standard vdM scan pair during the 2022 *pp* luminosity calibration. Separation steps and direction are indicated with gray markers. The last CMS scan is shown, with the highest beam-beam parameter of \\( \\xi =5.3 \\times 10^{-3} \\). The corresponding COMBI simulation results are also shown (triangles)
[Full size image](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8/figures/27)
## 7 Conclusions
Beamâbeam effects were directly observed and quantitatively validated at the LHC. The impact of beamâbeam interactions on head-on luminosity, bunch sizes, tune spectra, and beam orbit positions was measured across various collision configurations. The observed effects align closely with COMBI simulations. In particular, the luminosity enhancement due to amplitude-dependent \\(\\beta \\)\-beating was measured for the first time and is found to agree with expectations over a wide range of the beamâbeam parameter \\(\\xi \\). Together with the observed linear scaling with respect to the beamâbeam parameter, this provides experimental confirmation of the numerical studies presented in \[[13](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR13 "A. Babaev, T. Barklow, O. Karacheban, W. Kozanecki, I. Kralik, A. Mehta, G. Pasztor, T. Pieloni, D. Stickland, C. Tambasco, R. Tomas, J. WaĆczyk, Impact of beam-beam effects on absolute luminosity calibrations at the CERN large hadron collider. Eur. Phys. J. C 84, 17 (2024).
https://doi.org/10.48550/arXiv.2306.10394
")\]. Consistent signatures were found across multiple scan configurations, including two- and three-IP setups, and the characteristic tune shifts were well reproduced by both measurement and simulation.
Accurate instrumentation and experimental design were essential. Key to the results was the stepwise collision scheme, isolating individual beamâbeam contributions, with the implementation of the witness IP. Synchrotron light measurements from the BSRT, especially in the vertical plane of Beam 2, showed excellent agreement with COMBI predictions and bunch-by-bunch consistency. Tune spectra from the ADT system provided direct access to mean tune shifts, while orbit data from the DOROS system enabled measurement of beamâbeam deflection down to 1 \\(\\upmu \\)m, matching analytical estimates. These tools, combined with precise control of IP configurations, enabled a comprehensive cross-check of beamâbeam models.
Some systematic limitations were identified. During one of the tests, beam stability during the scans was insufficient to allow for symmetric modeling of luminosity decay in both transverse planes. In several scans, the nominal separation steps were based on overestimated beam sizes, resulting in larger-than-expected scan intervals and less sensitivity at low separation. Beamâbeam deflection measurements showed sign reversals at high separation, attributed to nonlinear orbit distortions from hysteresis effects in superconducting magnets. Phase advance errors between IPs and diagnostic systems further reduced the accuracy of simulations in some configurations.
Targeted improvements would allow for higher precision and new insights. Finer separation steps, especially below \\(3\\,\\sigma \\), would improve sensitivity to nonlinear beamâbeam effects and allow more accurate comparison with models. Suppressing coherent beamâbeam modes via phase advance tuning of not only two but all three IPs would improve tune shift measurements. Moreover, such optimal phase advance configuration between the IPs can be exploited to enhance luminosity. This effect could be directly propagated into at least a few percent increase in the total collected integrated luminosity, as discussed in Ref. \[[36](https://link.springer.com/article/10.1140/epjc/s10052-026-15476-8#ref-CR36 "J. WaĆczyk, T. Pieloni, Optimizing beam-beam beta-beating for luminosity enhancement at the LHC. JACoW IPAC 2025, 018 (2025)")\]. Aligning the phase advance between IPs and diagnostic devices, particularly the BSRT, would reduce sensitivity to optics errors, enabling better beam size measurements for both beams and transverse planes. Moreover, the suppression of the coherent beamâbeam modes facilitates the observation and analysis of the central, incoherent part of the tune spectrum, which is otherwise obscured by the dominant coherent peaks.
Further experimental strategies should also focus on directly probing the beamâbeam interaction at its point of occurrence. This would require access to the full transverse beam distributions with high spatial resolution and sensitivity, going beyond projected sizes or centroids. Such measurements could directly reveal the nonlinear distortions and tails induced by the beamâbeam force. In parallel, validating correction schemes that parametrize multi-IP beamâbeam effects on luminosity via equivalent tune shifts remains essential. These models underpin practical luminosity corrections in complex fill configurations and would benefit from systematic benchmarking.
## References
1. S. Meer, Calibration of the effective beam height in the ISR. Technical report, CERN, Geneva (1968). <https://cds.cern.ch/record/296752>
2. C. Rubbia, Measurement of the luminosity of \\(p{{\\overline{p}}}\\) collider with a (generalized) Van der Meer Method. Technical report, CERN, Geneva (1977). <https://cds.cern.ch/record/1025746>
3. ATLAS collaboration, Luminosity determination in \\(pp\\) collisions at \\(\\sqrt{s}=13\\) TeV using the ATLAS detector at the LHC. Eur. Phys. J. C **83**(10), 982 (2023). <https://doi.org/10.1140/epjc/s10052-023-11747-w>. [arXiv:2212.09379](http://arxiv.org/abs/2212.09379)
4. CMS Collaboration, Precision luminosity measurement in proton-proton collisions at sqrts = 13 TeV with the CMS detector. Technical report, CERN, Geneva (2025). <https://cds.cern.ch/record/2940794>
5. W. Herr, T. Pieloni, Beam-Beam Effects. Technical report, CERN (2014). Contribution to the CASâCERN Accelerator School: Advanced Accelerator Physics Course, Trondheim, Norway, 18â29 Aug 2013, p. 29. <https://doi.org/10.5170/CERN-2014-009.431>. <https://cds.cern.ch/record/1982430>
6. A.W. Chao, Beam-beam instability. AIP Conf. Proc. **127**, 201â242 (1985). <https://doi.org/10.1063/1.35187>
[Article](https://doi.org/10.1063%2F1.35187) [ADS](http://adsabs.harvard.edu/cgi-bin/nph-data_query?link_type=ABSTRACT&bibcode=1985AIPC..127..201C) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Beam-beam%20instability&journal=AIP%20Conf.%20Proc.&doi=10.1063%2F1.35187&volume=127&pages=201-242&publication_year=1985&author=Chao%2CAW)
7. E.D. Courant, H.S. Snyder, Theory of the alternating-gradient synchrotron. Ann. Phys. **3**, 1â48 (1958). <https://doi.org/10.1016/0003-4916(58)90012-5>
[Article](https://doi.org/10.1016%2F0003-4916%2858%2990012-5) [ADS](http://adsabs.harvard.edu/cgi-bin/nph-data_query?link_type=ABSTRACT&bibcode=1958AnPhy...3....1C) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Theory%20of%20the%20alternating-gradient%20synchrotron&journal=Ann.%20Phys.&doi=10.1016%2F0003-4916%2858%2990012-5&volume=3&pages=1-48&publication_year=1958&author=Courant%2CED&author=Snyder%2CHS)
8. P. Bambade, R. Erickson, W.A. Koska, W. Kozanecki, N. Phinney, S.R. Wagner, Observation of beam-beam deflections at the interaction point of the SLAC Linear Collider. Phys. Rev. Lett. **62**(25), 2949â2952 (1989). <https://doi.org/10.1103/PhysRevLett.62.2949>
[Article](https://doi.org/10.1103%2FPhysRevLett.62.2949) [ADS](http://adsabs.harvard.edu/cgi-bin/nph-data_query?link_type=ABSTRACT&bibcode=1989PhRvL..62.2949B) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Observation%20of%20beam-beam%20deflections%20at%20the%20interaction%20point%20of%20the%20SLAC%20Linear%20Collider&journal=Phys.%20Rev.%20Lett.&doi=10.1103%2FPhysRevLett.62.2949&volume=62&issue=25&pages=2949-2952&publication_year=1989&author=Bambade%2CP&author=Erickson%2CR&author=Koska%2CWA&author=Kozanecki%2CW&author=Phinney%2CN&author=Wagner%2CSR)
9. W.A. Koska, P. Bambade, W. Kozanecki, N. Phinney, S.R. Wagner, Beam-beam deflection as a beam tuning tool at the slac linear collider. Nucl. Instrum. Methods A **286**, 32 (1990). <https://doi.org/10.1016/0168-9002(90)90203-I>
[Article](https://doi.org/10.1016%2F0168-9002%2890%2990203-I) [ADS](http://adsabs.harvard.edu/cgi-bin/nph-data_query?link_type=ABSTRACT&bibcode=1990NIMPA.286...32K) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Beam-beam%20deflection%20as%20a%20beam%20tuning%20tool%20at%20the%20slac%20linear%20collider&journal=Nucl.%20Instrum.%20Methods%20A&doi=10.1016%2F0168-9002%2890%2990203-I&volume=286&publication_year=1990&author=Koska%2CWA&author=Bambade%2CP&author=Kozanecki%2CW&author=Phinney%2CN&author=Wagner%2CSR)
10. C. Bovet, M.D. Hildreth, M. Lamont, H. Schmickler, J. Wenninger, Luminosity optimisation using beam-beam deflections at LEP, in *Conf. Proc. C* (1996). <https://cds.cern.ch/record/306910>
11. D. Brandt, W. Herr, M. Meddahi, A. Verdier, Is lep beam-beam limited at its highest energy? in *Proceedings of the 1999 Particle Accelerator Conference (PAC â99), New York*, pp. 3005â3007. IEEE/JACoW. THP25; CERN report CERN-SL-99-030-AP (1999). <https://accelconf.web.cern.ch/p99/PAPERS/THP25.pdf>
12. W. Kozanecki, T. Pieloni, J. Wenninger, Observation of Beam-beam Deflections with LHC Orbit Data. Technical report, CERN (2013). <https://cds.cern.ch/record/1581723?ln=en>
13. A. Babaev, T. Barklow, O. Karacheban, W. Kozanecki, I. Kralik, A. Mehta, G. Pasztor, T. Pieloni, D. Stickland, C. Tambasco, R. Tomas, J. WaĆczyk, Impact of beam-beam effects on absolute luminosity calibrations at the CERN large hadron collider. Eur. Phys. J. C **84**, 17 (2024). <https://doi.org/10.48550/arXiv.2306.10394>
[Article](https://doi.org/10.48550%2FarXiv.2306.10394) [ADS](http://adsabs.harvard.edu/cgi-bin/nph-data_query?link_type=ABSTRACT&bibcode=2024EPJC...84...17B) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Impact%20of%20beam-beam%20effects%20on%20absolute%20luminosity%20calibrations%20at%20the%20CERN%20large%20hadron%20collider&journal=Eur.%20Phys.%20J.%20C&doi=10.48550%2FarXiv.2306.10394&volume=84&publication_year=2024&author=Babaev%2CA&author=Barklow%2CT&author=Karacheban%2CO&author=Kozanecki%2CW&author=Kralik%2CI&author=Mehta%2CA&author=Pasztor%2CG&author=Pieloni%2CT&author=Stickland%2CD&author=Tambasco%2CC&author=Tomas%2CR&author=Wa%C5%84czyk%2CJ)
14. T. Pieloni, A study of beam-beam effects in hadron colliders with a large number of bunches. PhD Thesis, Ăcole Polytechnique FĂ©dĂ©rale de Lausanne (2008). <https://doi.org/10.5075/epfl-thesis-4211>
15. ATLAS Collaboration, Improved luminosity determination in pp collisions at \\(\\sqrt{s} = 7\\) TeV using the ATLAS detector at the LHC. Eur. Phys. J. C **73**, 2518 (2013). <https://doi.org/10.1140/epjc/s10052-013-2518-3>
16. C.M.S. Collaboration, Development of the cms detector for the cern lhc run 3. J. Instrum. **19**(05), 05064 (2024). <https://doi.org/10.1088/1748-0221/19/05/p05064>
[Article](https://doi.org/10.1088%2F1748-0221%2F19%2F05%2Fp05064) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Development%20of%20the%20cms%20detector%20for%20the%20cern%20lhc%20run%203&journal=J.%20Instrum.&doi=10.1088%2F1748-0221%2F19%2F05%2Fp05064&volume=19&issue=05&publication_year=2024&author=Collaboration%2CCMS)
17. C. Barschel, M. Ferro-Luzzi, J.-J. Gras, M. Ludwig, P. Odier, S. Thoulet, Results of the LHC DCCT calibration studies. Technical Report CERN-ATS-Note-2012-026 PERF (2012). <https://cds.cern.ch/record/1425904?ln=en>
18. G. Anders et al., Study of the relative LHC bunch populations for luminosity calibration. Technical Report. CERN-ATS-Note-2012-028 PERF (2012). <https://cds.cern.ch/record/1427726>
19. G. Trad, Development and Optimisation of the SPS and LHC Beam Diagnostics Based on Synchrotron Radiation Monitors Presented 22 Jan 2015. Presented 22 Jan 2015. <https://cds.cern.ch/record/2266055>
20. R. TomĂĄs, O. BrĂŒning, M. Giovannozzi, P. Hagen, M. Lamont, F. Schmidt, G. Vanbavinckhove, M. Aiba, R. Calaga, R. Miyamoto, Cern large hadron collider optics model, measurements, and corrections. Phys. Rev. ST Accel. Beams **13**, 121004 (2010). <https://doi.org/10.1103/PhysRevSTAB.13.121004>
21. R. TomĂĄs, M. Aiba, A. Franchi, U. Iriso, Review of linear optics measurement and correction for charged particle accelerators. Phys. Rev. Accel. Beams **20**, 054801 (2017). <https://doi.org/10.1103/PhysRevAccelBeams.20.054801>
[Article](https://doi.org/10.1103%2FPhysRevAccelBeams.20.054801) [ADS](http://adsabs.harvard.edu/cgi-bin/nph-data_query?link_type=ABSTRACT&bibcode=2017PhRvS..20e4801T) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Review%20of%20linear%20optics%20measurement%20and%20correction%20for%20charged%20particle%20accelerators&journal=Phys.%20Rev.%20Accel.%20Beams&doi=10.1103%2FPhysRevAccelBeams.20.054801&volume=20&publication_year=2017&author=Tom%C3%A1s%2CR&author=Aiba%2CM&author=Franchi%2CA&author=Iriso%2CU)
22. O. Karacheban, P. Tsrunchev, on behalf of CMS, Emittance scans for cms luminosity calibration. EPJ Web Conf. **201**, 04001 (2019). <https://doi.org/10.1051/epjconf/201920104001>
23. F. Carlier, R. TomĂĄs, Accuracy and feasibility of the \\({\\beta }^{\*}\\) measurement for lhc and high luminosity lhc using \\(k\\) modulation. Phys. Rev. Accel. Beams **20**, 011005 (2017). <https://doi.org/10.1103/PhysRevAccelBeams.20.011005>
[Article](https://doi.org/10.1103%2FPhysRevAccelBeams.20.011005) [ADS](http://adsabs.harvard.edu/cgi-bin/nph-data_query?link_type=ABSTRACT&bibcode=2017PhRvS..20a1005C) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Accuracy%20and%20feasibility%20of%20the%20%24%24%7B%5Cbeta%20%7D%5E%7B%2A%7D%24%24%20%CE%B2%20%E2%88%97%20measurement%20for%20lhc%20and%20high%20luminosity%20lhc%20using%20%24%24k%24%24%20k%20modulation&journal=Phys.%20Rev.%20Accel.%20Beams&doi=10.1103%2FPhysRevAccelBeams.20.011005&volume=20&publication_year=2017&author=Carlier%2CF&author=Tom%C3%A1s%2CR)
24. M. SöderĂ©n, J. Komppula, G. Kotzian, S. Rains, D. Valuch, ADT and Obsbox in LHC Run 2, plans for LS2, pp. 165â171 (2019)
25. A. Boccardi, M. Gasior, R. Jones, P. Karlsson, R. Steinhagen, First results from the lhc bbq tune and chromaticity systems. Technical report, CERN, Geneva (2009). <https://cds.cern.ch/record/1156349>
26. M. Gasior, G. Baud, J. Olexa, G. Valentino, First Operational Experience with the LHC Diode ORbit and OScillation (DOROS) System, p. 07 (2017) <https://doi.org/10.18429/JACoW-IBIC2016-MOPG07>
27. A.W. Chao, Coherent beam-beam effects, in *Frontiers of Particle Beams: Intensity Limitations* (Springer, Berlin, Heidelberg, 1992), pp. 363â414. <https://doi.org/10.1007/3-540-55250-2_36>
28. J. WaĆczyk, A. Dabrowski, D. Stickland, T. Pieloni, W. Kozanecki, X. Buffat, Impact of multiple beam-beam encounters on LHC absolute-luminosity calibrations by the van der Meer method. JACoW IPAC **2023**, 053 (2023). <https://doi.org/10.18429/JACoW-IPAC2023-WEPA053>
[Article](https://doi.org/10.18429%2FJACoW-IPAC2023-WEPA053) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=Impact%20of%20multiple%20beam-beam%20encounters%20on%20LHC%20absolute-luminosity%20calibrations%20by%20the%20van%20der%20Meer%20method&journal=JACoW%20IPAC&doi=10.18429%2FJACoW-IPAC2023-WEPA053&volume=2023&publication_year=2023&author=Wa%C5%84czyk%2CJ&author=Dabrowski%2CA&author=Stickland%2CD&author=Pieloni%2CT&author=Kozanecki%2CW&author=Buffat%2CX)
29. J. WaĆczyk, Precision luminosity measurement at hadron colliders. PhD Thesis, Ăcole Polytechnique FĂ©dĂ©rale de Lausanne (2024). <https://doi.org/10.5075/epfl-thesis-10500>
30. T. Pieloni, J. Barranco, X. Buffat, P.C. Jorge, L.M. Medrano, C. Tambasco, R. TomĂĄs, Dynamic beta and beta-beating effects in the presence of the beam-beam interactions in *57th ICFA Advanced Beam Dynamics Workshop on High-Intensity and High-Brightness Hadron Beams*, p. 027 (2016). <https://doi.org/10.18429/JACoW-HB2016-MOPR027>
31. F. Carlier, et al. Challenges of K-modulation measurements in the LHC Run 3, in *Proceedings of IPACâ23. IPACâ23â14th International Particle Accelerator Conference* (JACoW Publishing, Geneva, Venice, 2023), pp. 531â534. <https://doi.org/10.18429/JACoW-IPAC2023-MOPL014>
32. Y. Alexahin, A study of the coherent beam-beam effect in the framework of the Vlasov perturbation theory. Nucl. Instrum. Methods Phys. Res., Sect. A **480**(1â3), 253â288 (2002). <https://doi.org/10.1016/S0168-9002(01)01219-0>
[Article](https://doi.org/10.1016%2FS0168-9002%2801%2901219-0) [ADS](http://adsabs.harvard.edu/cgi-bin/nph-data_query?link_type=ABSTRACT&bibcode=2002NIMPA.480..253A) [Google Scholar](http://scholar.google.com/scholar_lookup?&title=A%20study%20of%20the%20coherent%20beam-beam%20effect%20in%20the%20framework%20of%20the%20Vlasov%20perturbation%20theory&journal=Nucl.%20Instrum.%20Methods%20Phys.%20Res.%2C%20Sect.%20A&doi=10.1016%2FS0168-9002%2801%2901219-0&volume=480&issue=1%E2%80%933&pages=253-288&publication_year=2002&author=Alexahin%2CY)
33. R. Calaga, W. Herr, G. Papotti, T. Pieloni, X. Buffat, S. White, R. Giachino, Coherent beam-beam mode in the LHC. CERN (2014). <https://doi.org/10.5170/CERN-2014-004.227>. <http://cds.cern.ch/record/1957039>
34. K. Yokoya, Y. Funakoshi, E. Kikutani, H. Koiso, J. Urakawa, Tune shift of coherent beam-beam oscillations. Part. Accel. **27**, 181â186 (1990)
[Google Scholar](http://scholar.google.com/scholar_lookup?&title=Tune%20shift%20of%20coherent%20beam-beam%20oscillations&journal=Part.%20Accel.&volume=27&pages=181-186&publication_year=1990&author=Yokoya%2CK&author=Funakoshi%2CY&author=Kikutani%2CE&author=Koiso%2CH&author=Urakawa%2CJ)
35. M. Bassetti, G.A. Erskine, Closed expression for the electrical field of a two-dimensional gaussian charge. Technical report, CERN, Geneva (1980). <https://cds.cern.ch/record/122227>
36. J. WaĆczyk, T. Pieloni, Optimizing beam-beam beta-beating for luminosity enhancement at the LHC. JACoW IPAC **2025**, 018 (2025)
[Google Scholar](http://scholar.google.com/scholar_lookup?&title=Optimizing%20beam-beam%20beta-beating%20for%20luminosity%20enhancement%20at%20the%20LHC&journal=JACoW%20IPAC&volume=2025&publication_year=2025&author=Wa%C5%84czyk%2CJ&author=Pieloni%2CT)
[Download references](https://citation-needed.springer.com/v2/references/10.1140/epjc/s10052-026-15476-8?format=refman&flavour=references) |
| Shard | 129 (laksa) |
| Root Hash | 17645177711233004329 |
| Unparsed URL | com,springer!link,/article/10.1140/epjc/s10052-026-15476-8 s443 |