Electromagnetic dipole moments of charged baryons with bent crystals at the LHC.

We propose a unique program of measurements of electric and magnetic dipole moments of charm, beauty and strange charged baryons at the LHC, based on the phenomenon of spin precession of channeled particles in bent crystals. Studies of crystal channeling and spin precession of positively- and negatively-charged particles are presented, along with feasibility studies and expected sensitivities for the proposed experiment using a layout based on the LHCb detector.

Introduction

The magnetic dipole moment (MDM) and the electric dipole moment (EDM) are static properties of particles that determine the spin motion in an external electromagnetic field, as described by the T-BMT equation [1-3]. Several measurements of baryon MDMs contributed to confirm the validity of the quark model [4]. Measurements of the MDM of heavy baryons, i.e. baryons containing charm or beauty quarks, have never been performed due to the difficulties imposed by the short lifetime of these particles of about . These measurements would provide important anchor points for QCD calculations, helping to discriminate between different models [5, 6], and would improve the current understanding of the internal structure of hadrons. The possibility to measure the MDM of positively-charged charm baryons at the Large Hadron Collider (LHC) using bent crystals has been proposed in Refs. [7, 8] and recently revisited [9, 10].

The EDM is the only static property of a particle that requires the violation of parity (P) and time reversal (T) symmetries and thus, relying on invariance, the violation of symmetry. The EDM of a baryon may arise from the structure of quarks and gluons, and any process involving a photon and a flavour–diagonal coupling. In the Standard Model (SM), contributions to the EDM of baryons are highly suppressed but can be largely enhanced in some of its extensions. Hence, the experimental searches for the EDM of fundamental particles provide powerful probes for physics beyond the SM.

Indirect bounds on charm (beauty) quark EDM are set from different experimental measurements and span over several orders of magnitude, i.e. charm (beauty) EDM  [11-15] ( [13-16]), depending on different models and assumptions. As an example, an indirect bound on the charm quark EDM is derived from the experimental limit on the neutron EDM to be  [11], and a charm quark EDM of comparable magnitude is possible in extensions of the SM [12]. For the beauty quark, indirect EDM limits  [13] and  [16] are derived, and a relatively large beauty quark EDM is possible in presence of new physics.

Recently, it has been proposed to search for the EDM of positively-charged charm baryons using bent crystals at the LHC  [17]. Similarly to the MDM case, the method relies on baryons produced by the interaction of protons, extracted from the LHC beam, on a fixed target. The baryons are subsequently channeled in a bent crystal. The spin precession of short-lived particles is induced by the intense electromagnetic field between the crystal atomic planes. The EDM and the MDM information can be extracted by measuring the spin polarization of the channeled baryons the end of the bent crystal. This technique can be extended to strange and beauty positively-charged baryons.

In this paper we address several key aspects of this unique experimental program. In Sects. 2 and 4, after introducing the channeling of charged particles in bent crystals, we study the deflection and the spin precession of positively- and negatively-charged baryons in a bent crystal using Geant4 simulations. We assess the experimental technique and study the possibility to extend the EDM searches and MDM measurements to negatively-charged baryons. This would allow to perform tests of the symmetry by measuring the MDM of particles and antiparticles. In Sect. 3 we prove that the spin evolution equations describing MDM and EDM effects hold for non-harmonic planar channel potential, therefore spin precession effects for both positively- and negatively-charged particles depend uniquely on the crystal curvature. In fact, the same equations also apply for axial-channeled particles, mostly relevant for negatively-charged particles, although its application to MDM and EDM physics will require further investigation. Section 5 focuses on the description of a possible fixed-target setup installed in front of the LHCb detector. The feasibility of the measurements has been evaluated relying on both parametric and Geant4 simulations along with a geometrical model of the detector. Finally, Sect. 6 presents sensitivity studies for EDM searches and MDM measurements of , charm baryons, , beauty antibaryons, and , strange antibaryons. Baryons and antibaryons will be referred hereafter generically as baryons, unless otherwise stated.

Channeling of multi-TeV charged particles

In a crystal the strong electric field experienced by a charged particle in the proximity of the ordered structure of the atoms exerts a strong confinement force onto the particle itself. The particle trajectory can be bound to stay parallel to a crystalline plane or to an atomic string, which becomes a preferential pathway in the crystal. This phenomenon is called channeling and can occur if the angle between the particle trajectory and a crystal plane (planar channeling) or a crystal axis (axial channeling) is lower than a Lindhard angle , where is the potential-well depth, p the particle momentum and its velocity [18, 19]. This process has been studied in laboratory up to the highest available energy, in particular at the LHC where the planar channeling of 6.5 protons has been observed [20].

When the crystal is bent, its planes or atomic strings are bent too. The particle pathways are adiabatically bent following the crystal curvature, resulting in a net deflection of the incoming direction by an angle equal to that of crystal bending: charged particle steering is then possible through channeling in bent crystals. Various applications as circular accelerator halo collimation [21-25] or beam extraction from an accelerator ring for fixed-target experiments [26] have been studied and proposed also for the LHC.

Beam steering of positively-charged particles (positive particles in the following) with channeling has progressed significantly over the last years, featuring silicon crystals with about 80% deflection efficiency at an energy of several hundred GeV [27, 28]. Positive particles in channeling condition are repelled by the atomic electric field and follow trajectories that tend to be far from the lattice sites. Negative particles, on the contrary, are attracted by the same field and can repeatedly oscillate across the nuclei of the crystal. For this reason negative particles are more likely to collide with the nuclei of the crystal lattice and therefore can easily escape from a channeling bound state. The average length that a channeled particle traverses before exiting from the planar or axial potential well is called dechanneling length , being much shorter for negative than for positive charges.

The first measurements of for ultra-high energy negative particles has been done at CERN using relatively short bent silicon crystal at few hundreds . It has been measured to be few [29-31] using negative hadron and electron beams, and successfully compared with simulations. Therefore, simulations can be used to extrapolate the efficiency of channeling processes to the multi-TeV energy range and for different crystal curvature radii.

For the ultra-relativistic energy range it was also demonstrated the ability to steer negative particle beams through the axial channeling regime with an efficiency above 90% [32, 33]. Indeed, in the case of axial channeling, there are particles which are above the electric field barrier and are not trapped along a single atomic string. Those particles are anyway deflected due to their stochastic scattering with different atomic strings of the crystal, avoiding the fast dechanneling occurring for negative particles. However, reaching a condition of axial alignment for beam steering is relatively more difficult than for planar channeling since the orientation of the crystal with respect to two (and not only one) rotational axes must be found.

The dependence of the channeling efficiency for positive particles on the particle energy, the crystal length, and bending radius is well known for crystals with the length along the beam comparable to [19]. The dechanneling length scales almost proportionally to the particle momentum-velocity , i.e. , and can be calculated aswhere and are the mass and the classical radius of the electron, I is the mean ionization energy, is the interplanar spacing, is the Thomas-Fermi screening radius, and and are the charge number and boost of the incident particle.

The channeling efficiency under harmonic approximation scales proportionally to , with , where R is the crystal bending radius and is the critical radius for channeling, i.e. the minimum bending radius for which channeling occurs, and holds,Since the fraction of particles which remains channeled in a bent crystal has to oscillate near the potential well edge, the probability to leave the channeling state increases for small bending radii. Under harmonic approximation for a bent crystal the dechanneling length is shortened by a factor ,Differently, the energy dependence of for negative particles is not theoretically known. Contrarily to the dechanneling of positive particles, the large variation of the transverse energy by negative particles in the interaction with atomic nuclei cannot be treated as a stochastic process, and requires a new formalism to be developed.

To quantify the deflection efficiency of the channeling process in bent crystals, Monte Carlo simulations for both positive and negative particles have been carried out using the Geant4 toolkit [34, 35] in a version allowing for crystalline structures [36]. The channeling process is implemented by including Dynecharm++ [37] and Echarm [38] into the Geant4 channeling package [39]. Such model was validated against experimental data for negative pions [30], protons [40-46], electrons [31, 47–50] and positrons [31] in a range of energies spanning from (electrons at MAMI) to (protons at CERN). Dynecharm++ allows the tracking of a relativistic charged particle inside a crystalline medium via the numerical integration of the T-BMT classical equations of motion [1-3]. The continuum potential approximation proposed by Lindhard is used [18]. Echarm allows the computation if the electrical characteristics of the crystal is within this approximation.

The Geant4 application was developed on top of the 10.3 version of the toolkit, which allows for crystalline structures [36]. The Geant4 physics lists used were the G4HadronElasticPhysics and the G4HadronPhysicsFTFP_BERT and a custom G4EmStandardPhysics_option4 with single scattering instead of multiple-scattering.

Figure 1 shows the dependence of the channeling efficiency for protons and antiprotons with the particle momentum for a 1, 1, 7 long crystal bent along the (110) plane by a 0.1, 1 and 14 bending angle, respectively. The efficiency for positive particles is not spoiled by , being cm for all the momenta, but the unfavourable ratio between the critical radius and the bending radius is causing it to decrease. Indeed, as illustrated in Fig. 2, this ratio rapidly increases, lowering the deflection efficiency. For negative particles, the deflection efficiency is largely dominated by the crystal length. Therefore, the efficiency remains always lower than for positive particles for all the momenta. Such simulations show that the crystal geometric parameters have to be carefully chosen depending on the energy range in which the crystal has to be operated. In the figure, particles which are not captured under channeling at the crystal entrance are reflected to the opposite side with an angle which depends on the particle momentum [51].

Fig. 1

Dependence of the channeling efficiency of protons and antiprotons with the particle momentum for 1 , 1 and 7 long crystals bent along the (110) plane by a 0.1, 1 and 14 bending angle, respectively. The curves for the anti-proton interacting with the 1 and 7 long crystals are superimposed

Fig. 2

Outgoing angular distributions at various momenta for protons and antiprotons impinging on a 1 long crystal bent along the (110) plane by a 0.1 bending angle. The dechanneling length for positive particles in straight crystals () and in the bent crystal (), along with the critical radius for channeling (), are calculated for the different momenta following Ref. [19]

Spin precession

The spin precession of a charged particle is induced by the interaction of its electromagnetic dipole moments, e.g. MDM and EDM, with external electromagnetic fields. The time evolution of the spin-polarization vector is regulated by the T-BMT equationwhere the precession angular velocity vector is composed by three contributions corresponding to the MDM, EDM, and Thomas precession:where and are the electric and the magnetic fields in the laboratory frame, and q, and are the electric charge, boost and vector velocity of the particle, respectively. The g and d dimensionless factors, also referred to as the gyromagnetic and gyroelectric ratios, define the magnetic and electric dipole moment of a particle with spin J (in Gaussian units) as and , respectively, where is the particle magneton.1

The lifetime of baryons with heavy quark constituents is too short, e.g. the baryon lifetime is about , for a standard magnet to induce any detectable effect to the spin-polarization vector before they decay. The possibility to measure the MDM of short-lived baryons using channeling in bent crystals was firstly pointed out by V. G. Baryshevsky in 1979. The method is based on the interaction of the MDM of the channeled particles with the intense electric field between crystal atomic planes. As an example, a sketch of the deflection of the baryon trajectory and spin precession in a bent crystal is shown in Fig. 3. Charm baryons produced by interaction of protons on a fixed target, e.g. a target, are polarized perpendicularly to the production plane due to parity conservation in strong interactions. The production plane xz, also shown in the figure, is determined by the proton and baryon momenta; the latter defines the z axis. The initial polarization vector is perpendicular to the production plane, along the y axis, and at the end of the crystal it is rotated of an angle . The crystal is bent in the yz plane by an angle . The measurement of the MDM of charm baryons using bent crystals has been widely discussed since the 80’s [52-54]. Lately, the possibility of the measurement at LHC energies has been considered [7-9]. Recently, the search for charm baryon EDM using bent crystals at LHC has been proposed [17].

Fig. 3

Sketch of the deflection of the baryon trajectory and spin precession in a bent crystal. The initial polarization vector is perpendicular to the production plane, along the y axis, due to parity conservation in strong interactions. The spin precession in the yz and xy plane are induced by the MDM and the EDM, respectively. The red (dashed) arrows indicate the (magnified) spin-polarization component proportional to the particle EDM. The angle indicates the spin precession due to the MDM

The spin precession of particles channeled in bent crystals was firstly observed by the E761 Collaboration [55]. Using a 800 proton beam impinging on a Cu target, baryons with 375 average momentum were produced and channeled in two bent crystals with opposite bending angles. The MDM of the baryon was measured and proved the viability of this technique for the measurement of the MDM of short-lived particles.

Planar channeling

In the case of planar channeling, the intense electric field between the crystal planes, , which deflects charged particles, transforms into a strong electromagnetic field , in the instantaneous rest frame of the particle and induces spin precession. In the limit of large boost, the spin precession induced by the MDM in the yz plane is [56]- In order to obtain rad are required large crystal bending angles , and large Lorentz factors , which can be uniquely achieved at LHC.

The equations describing the spin precession of planar channeled positive particles in presence of MDM and EDM are derived in Ref. [17]. In the limit of large boost, and assuming small EDM effects compared to the main MDM spin precession, a polarization component orthogonal to the bending plane is induced,The MDM driven precession taking place in the bending plane is given byInside the crystal, positive particles feel a non-zero mean electric field thanks to the centripetal force induced by the crystal bending, and the spin precession depends basically on . The planar channel potential seen by positive particles can be assumed to be approximately harmonic, as described in Sect. 2. For negative particles this assumption is no more valid since their motion is regulated by a non-harmonic potential, as shown in Fig. 4.

Fig. 4

Harmonic (red dash-dotted line) and non-harmonic (blue continuous line) electric potential versus interplanar distance for positive and negative particles in a (110) crystal. The electric potential is extracted from Geant4 simulations. For the sake of comparison the electric potential for negative particles is shifted by half of the interplanar distance

In the following, we demonstrate that in presence of a non-harmonic potential V, identical spin precession equations derived for the harmonic potential case hold. We consider the layout of Fig. 3, with the crystal bent along an atomic plane. Polar coordinates are introduced for describing the particle trajectory in the bending planewhere is the revolution frequency for the particle traversing the bent crystal, and the electric field described by the planar channel potential isNeglecting EDM contributions, the spin evolution resulting from Eqs. (4), (5), (9) and (10) isfor the initial condition and whereThe radial coordinate is constant up to , therefore the spin precession depends on . Over a complete oscillation in the channel potential the effect of this term is equivalent to that of the electric field in the particle equilibrium radial position ,which is determined solely by the centripetal force induced by the bent trajectory,This statement follows by computingfor a complete particle oscillation. By changing the integration variable to and , then is determined by the non-relativistic energy conservation for the radial motion of channeled particles [19]in which and is the total radial energy, assumed to be constant during a particle oscillation. The relation holds because the longitudinal motion is ultra-relativistic and independent from the radial one, which is non-relativistic since the potential depth is . The integration boundaries are chosen to be the particle oscillation limits, in whichFinally, the integral can be trivially computedSummarizing, spin precession effects given by the actual shape of the planar channel potential cancel out at each particle oscillation and the net spin precession depends uniquely on the crystal curvature. This result generalises the same conclusion previously obtained for harmonic potentials [52, 56]. The spin evolution equations describing MDM and EDM effects, Eqs. (6), (7) and (8), hold as for an harmonic planar channel potential, and in particular for the potential seen by negative particles.

Axial channeling

The planar channeling efficiency for negative particles is smaller than for positive. It becomes negligibly small for crystals longer than with bending angle larger than , as shown in Fig. 1. In this case, the spin depolarization effect of particles scattered by crystal axes (planes) provides a possibility for measuring the MDM of negative particles [10]. In case of axial alignment, spin rotation can also be investigated. The phenomenon of axial channeling has been observed for positive and negative particles but it has not been considered for spin precession to date. Here we discuss the possibility to induce spin precession in axial-channeled particles for potential applications in MDM and EDM measurements of charged baryons. We consider the same layout of Fig. 3, in which the crystal is now bent along a crystallographic axis. Polar coordinates are introduced for the bending plane as in Eq. (9) and the electric field described by the axial channel potential isIn analogy with planar-channeled particles described in Sect. 3.1, the longitudinal velocity is ultra-relativistic. The velocity components orthogonal to the channel are non-relativistic; their contribution to the spin precession described in Eq. (5), is negligible. The particle velocity inside the bent crystal is therefore simplified asin which and to a very good approximation for a crystal bending angle mrad. The spin evolves according to the spin precession equationwith precession vector following from Eqs. (5), (19) and (20),The main difference with respect to the planar channeling case is the presence of the electric field component, which in principle complicates the separation between the MDM and EDM induced spin rotation. Nonetheless, in the following it is shown that the contribution of the dV / dx terms can be neglected and the spin precession evolution derived for the planar channeling case applies also to axial-channeled particles.

During a particle oscillation the spin can be assumed to be constant since the typical spin precession frequency Hz is three orders of magnitude lower than the oscillation frequency of the particle trapped in the channel, Hz. The two dominant components and , describing spin precession in the yz and xz planes, respectively, can be considered to act independently of each other; namely, the spin rotation in the yz plane is not influenced by the spin rotation in the xz plane and vice versa. In this case Eq. (13) can be applied to both contributions: while the centripetal force induces a net spin precession in the yz plane identical to that of planarly-channeled particles, the effect of dV / dx mediates to zero over each particle oscillation, since no centripetal force acts in the x direction.

The limit of the employed assumption is checked estimating the typical amount of spin precession accumulated during an incomplete particle oscillation, which may lead to an imperfect cancellation of the dV / dx contribution. This amount is at the order ofin which the integration is carried on half of an oscillation. Here, is taken with and the mass [4]. The typical electric field magnitude of the axial channel is estimated as the ratio between the potential depth and the channel width for a Ge crystal, with values taken from Ref. [19]. The oscillation frequency is for , where the constant describing the potential curvature for a Ge axial channel is about eight times the value for a axial channel, according to Ref. [19].

Neglecting EDM effects (), Eqs. (21) and  (22) show that dV / dx contributes to the spin component via . Since the spin component is not constant during a complete oscillation, the dV / dx contribution is not exactly zero and can be conservatively estimated in Eq. (24), in which changes by an amount of order ,The integrated effect along the whole crystal is conservatively estimated by multiplying by the number of particle oscillations,in which a crystal length of is taken. Such a component does not affect the main MDM spin precession in the yz plane and it is negligible compared to the experimental sensitivity on the particle polarization. Indeed, according to the sensitivity studies detailed in Sect. 6, the uncertainty on the spin-polarization component, constituting the EDM signature, will be at the order of .

In summary, the spin evolution equations describing MDM and EDM effects, Eqs. (6), (7) and (8), hold also for axial-channeled particles, which is mostly relevant for negative particles with relatively high axial channeling efficiency. Nevertheless, the application of this result to the measurement of the MDM and EDM of particles has to be further studied.

Geant4 simulations of spin precession

The simulation method based on the numerical integration of the classical equations of motion allows to introduce the modification of the particle spin under the effect of the strong electric field generated by the crystalline lattice. Indeed, the step-by-step variation of the spin is tracked in Geant4 by numeric integration of the T-BMT equation [1-3].

The Geant4 application for spin precession has been validated against the solely available experimental data provided by the E761 experiment at FNAL [55]. In that experiment, two 4.5 long crystals bent along the (111) plane were exposed to a beam with 375 GeV/c momentum. The deflection angle of the two crystals were and , with measured precession angles of and , respectively. As expected, the spin in the two crystals precesses in opposite directions. The average of experimental values is consistent with the predicted value of . A uniformly bent crystal with the E761 geometrical parameters has been implemented within Geant4 and exposed to a monochromatic and perfectly collimated beam with 375 momentum. Figure 5 shows the distributions of the trajectory deflection angle and of the spin precession angle for both the up- and down-bending crystals. A precession of and was obtained for the two cases, respectively, in good agreement with the predicted values of .

Fig. 5

Distributions of the trajectory deflection angle and spin precession angle for baryons of 375 momentum interacting with 4.5 long (top) up-bent and (bottom) down-bent crystals, uniformly bent along the (111) plane at angle. Similar crystals were used for the E761 experiment at FNAL [55]

The same Geant4 toolkit can also be used to simulate the spin precession of other positive and negative particles in a bent crystal, to compare with the expected analytical values. For this purpose, a 1 long crystal bent along the (110) plane by a 1 bending angle has been used. A beam of short-lived particles with no angular divergence is generated in the simulation immediately before the crystal and the precession angle at the end of the crystal is evaluated. Table 1 shows the simulation results for pairs of short-lived particles and their antiparticles, i.e.  , , , and , in presence of MDM and EDM, respectively. Figure 6 shows the distribution of the measured trajectory deflection angles and spin precession angles for the / case.

Table 1

Average spin precession angle () and EDM polarization component () obtained from Geant4 simulation compared to the expected values ( and , respectively), due to the gyromagnetic factor of the particle expressed as and the gyroelectric factor , along with the mean channeling deflection efficiency (), for different 1 particles impinging on a 1 long crystal bent along the (110) plane at 1 angle. The normalization of the polarization vector has been taken unity, i.e.

Particle	g^′	Φ_exp [^∘]	Φ_sim [^∘]	s_x,exp	s_x,sim	ε_c [%]
Λc+	-0.30	-7.518	-7.474±0.015	7.17×10^-4	(7.19±0.03)×10^-4	71.0±0.08
Λ¯c-	0.30	7.518	7.59±0.07	-7.17×10^-4	(-7.20±0.13)×10^-4	0.51±0.07
Ξ^-	-1.92	-83.09	-83.0±0.9	1.132×10^-2	(1.145±0.020)×10^-2	0.47±0.07
Ξ¯^+	1.92	83.09	83.21±0.23	-1.132×10^-2	(-1.149±0.005)×10^-2	70.6±0.08
Ω^-	-2.20	-75.38	-75.4±0.6	8.50×10^-3	(8.51±0.12)×10^-3	0.39±0.6
Ω¯^+	2.20	75.38	75.53±0.15	-8.50×10^-3	(-8.51±0.03)×10^-3	70.9±0.8
Ξb-	-1.38	-13.65	-13.64±0.14	5.154×10^-4	(5.15±0.10)×10^-4	0.51±0.07
Ξ¯b+	1.38	13.65	13.78±0.03	-5.154×10^-4	(-5.167±0.021)×10^-4	71.0±0.8

Fig. 6

Distributions of the trajectory deflection angle and spin precession angle for (top) and (bottom) baryons of 1 momentum interacting with a 1 long crystal uniformly bent along the (110) plane at 1 angle

The experiment

The MDM and EDM information can be extracted using Eqs. (6), (7) and (8), from the measurement of the spin polarization of channeled baryons at the exit of the crystal, via the study of the angular distribution of final state particles. For decaying to two-body final states such as , , and , the angular distribution is described bywhere is a parity-violating coefficient depending on the final state f, the direction of the final state baryon in the helicity frame, the corresponding solid angle, and the polarization vector. Equation (26) holds for all spin-1/2 decays into two-body final states with spins , and .

The initial polarization would require in principle the measurement of the angular distribution for unchanneled baryons. In practice, however, this is not required since the measurement of the three components of the final polarization vector for channeled baryons allows a simultaneous determination of g, d and , up to discrete ambiguities, as discussed in Appendix A. These, in turn, can be resolved exploiting the dependence of the angular distribution with the boost .

Possible experimental layout

The possible experimental layout is based on the double crystal scheme [8] sketched in Fig. 7. It consists of four main elements:To protect against radiation damage and minimise detector occupancies, the design has to guarantee that non-interacting protons and unchanneled particles follow the beam pipe towards the absorbers.

Fig. 7

Conceptual layout of the fixed-target setup shown in side view with down- and up-bending crystals. The zoom in shows the spin precession in the down-bending crystal for channeled baryons

two crystal kickers positioned about upstream of the target for up and down deflection at angle of the 7 protons from the LHC beam halo. This technique has been demonstrated to be feasible without affecting the LHC beam lifetime [20];

two amorphous targets about thick intercepting the deflected proton beam where charm, beauty and strange baryons are produced. The fixed target has to be installed in front of the detector, as close as possible to obtain good vertex resolution;

up- and down-bending crystals to induce opposite spin precession to channeled baryons. The use of two crystals with opposite bendings is crucial to prove the robustness of the results and control systematic uncertainties. The target should be attached to the crystal to maximize the yield of channeled baryons;

two absorbers positioned downstream of the detector to stop the deflected proton beam and background particles induced by the interactions with the target and crystal materials [57].

Despite its challenges, the setup is based on two key elements already existing and tested successfully at the LHC: high-purity bent crystals and high-accuracy positioning systems (goniometers). Two bent crystal types with different characteristics are required. The crystal kicker is very similar to that one tested at the LHC  [20], of 4 length and bending angle; the second crystal should have larger angle, order , as discussed later. The remotely controlled goniometers equipped with the bent crystals are mounted on standard collimation supports and make use of fast plug-in technology, which ensures fast handling of the object. They are based on a piezoelectric actuator and feature angular resolution of and linear resolution of . This is necessary to align the crystal kicker with respect to the beam halo, and position the long-bending crystal to intercept the deflected beam.

Two possible configurations have been considered for the fixed-target and detector setup, referred to hereafter as and . The former is based on the upgraded LHCb detector [58], which will become operational in 2021 and will run for the rest of the decade (LHC Run 3 and Run 4), whereas the latter is an hypothetical dedicated detector considered to function at even higher luminosities and providing an angular coverage to minimise the crystal bending angle.

LHCb is a single-arm forward spectrometer [59, 60] dedicated to the study of particles containing or quarks at LHC. The detector tracking system consists of three main devices: a vertex locator (VELO) surrounding the interaction region and a large-area detector (TT) upstream of the dipole magnet inducing an integrated field of about 4 T, and three stations (T1–T3) downstream of the magnet. Particle identification is provided by two ring-imaging Cherenkov detectors, a calorimeter system and muon chambers. With the upgrade, most of the sub-detectors will be replaced [61-63] and a full software based trigger will become operational [64], providing significantly increased efficiencies in hadronic final states and allowing the experiment to operate at higher luminosities. The fixed-target like geometry combined with the capability to reconstruct with good efficiency highly-boosted baryons makes of LHCb the most suitable detector for this proposal.

To minimise the impact on the interaction region, in the scenario the fixed-target setup is positioned outside of the VELO vessel container [61], upstream of the nominal collision point and upstream of the first VELO module. At this position, a minimal crystal bending angle of about 14 is required for channeled baryons to be deflected inside the detector fiducial volume. This follows from a Monte Carlo simulation of fixed-target events using Pythia  [65] together with a simplified geometrical model of the upgraded LHCb detector, as sketched in Fig. 8. For the scenario, a beam intensity of impinging the target, and an overall data taking efficiency of 50% is assumed. In this case in six weeks of dedicated detector operations, spanned over several years during the next decade, it would be possible to achieve a statistics of about protons on target (). Accurate studies of the attainable proton flux are in progress [57], which has to be compliant with the LHC beam lifetime, machine operations and protection. Detector occupancies at such rates are expected to be manageable, as discussed later.

Fig. 8

Sketch of a fixed-target event generated with Pythia at position with a crystal bending , as seen by the simplified LHCb detector geometry model based on Refs. [61, 62] with a conservative beam clearance of the downstream stations of 30. The first VELO module is located at , upstream of the nominal LHC collision point at . The points represent the hits of the proton (green), pion (violet) and kaon (red) tracks overlaid in the (top left) side view of the VELO, (top right) front of the last VELO module, (bottom left) schematic side view of the whole LHCb detector, and (bottom right) central area of the T1 station. Events are considered within acceptance when they cross at least three VELO modules and the three T stations. The track bending due to the LHCb dipole magnet is taken into account

An increase in proton fluxes, e.g. at High Luminosity LHC (HL-LHC) [66], combined with the design of a dedicated detector capable to afford the higher occupancy levels and longer data taking periods, could potentially offer the opportunity for scenario to integrate . Such detector could also extend the angular coverage at larger pseudorapidity to minimise the crystal bending angle, increasing the channeling efficiency, and operate closer to the fixed-target to provide enhanced vertex resolution. Nevertheless, a minimal bending angle of about 5 is needed to guarantee a good separation between channeled and unchanneled particles. This value can be inferred from the angular distribution of baryons produced by 7 protons impinging on the fixed-target, which are highly collimated along the incident proton beam direction, and are isotropically distributed over the azimuthal angle. The polar angle that defines the emission cone is . The corresponding Pythia distribution of polar angle versus momentum, illustrated in Fig. 9, shows that for and momentum higher than there are practically no unchanneled baryons.

Fig. 9

Distribution of polar angle versus momentum for baryons produced in 7 proton beam collisions on protons at rest using Pythia

Crystal parameters

In Sect. 2 is discussed how the channeling efficiency depends on the crystal parameters and on the momentum range of the particles. In the following, the length L and bending angle for and crystals are estimated maximizing the sensitivity of the experiment to the electromagnetic moments, while taking into account detector acceptance. The optimization has been performed using fixed-target events, produced in 7 proton beam collisions on protons at rest using Pythia. The channeling has been simulated using a parameterisation based on current theoretical description and channeling measurements, following Ref. [19].

A particle entering the crystal is channeled when its polar angle lies within the interval, where is the Lindhard angle, introduced in Sect. 2. For in Si (Ge) this angle is about 6 (7), about three orders of magnitude smaller than the angular divergence of baryons shown in Fig. 9. Therefore, the trapping efficiency is largely dominated by the angular opening of the baryons produced in the target. This imposes the crystal to be directly attached to the target. The per-event deflection efficiency is parameterised aswhere dechanneling losses inside the bent crystal are described by the factor introduced in Sect. 2, which accounts for the shortening of the dechanneling length with respect to a straight crystal, and , which is the ratio of in a straight crystal to . The per-event critical radius, , must be below the crystal curvature radius , where is the interplanar electric field at the critical transverse coordinate, below which the particle is lost from channeling mode. For Si 110 (Ge 110),  Å and /  [19].

A scan in the plane is performed to determine the minimal error on the d and g factors separately. Since there is a wide momentum distribution for the channeled baryons, we generate and fit pseudo-experiments using a conditional probability density function constructed with the angular distribution in Eq. (26). Following the discussion at the end of Sec. 5.1, we require the momentum of the baryons to be higher than 800 . The dependence of the spin polarization is obtained from Eqs. (7) and (8), assuming , , and , as described in more details in Sect. 6. Figure 10 shows regions whose uncertainties on d (g) are increased by 20% with respect to the minimum, for Si and Ge in both and scenarios. The gyromagnetic ratio prefers higher L values, as . These wide regions provide the optimal parameters summarized in Table 2, chosen around the minimum of the gyroelectric factor.

Fig. 10

Regions of minimal uncertainty of d- and g-factors as a function of the crystal parameters L and , for (left) Si and (right) Ge in (top) and (bottom) scenarios. The markers and continuous (dotted) lines represent the minimum uncertainty and regions whose uncertainties on d (g) are increased by 20% with respect to the minimum, respectively

Table 2

Bent crystal parameters for and optimized using charm baryon decays, for the two possible experimental scenarios under consideration. The intervals give approximate regions whose uncertainties on the d factor are increased by 20% with respect to the minimum, whereas the central values are chosen for the sensitivity studies discussed in Sect. 6

	S1		S2
	Si	Ge	Si	Ge
L [cm]	[6, 12]	[4, 8]	[8, 15]	[5, 12]
	7	5	12	7
θ_C [mrad]	[13, 16]	[13, 16]	[5, 9]	[6, 12]
	14	15	7	8
R/R_c	2.85	3.33	9.77	8.75

The parametric approach adopted above to account for dechanneling losses has been validated with Geant4 simulations, discussed in Sects. 2 and 4, using 1 baryons with no angular divergence. For a 7 (5) long () crystal bent along the (110) plane by a 14 (15) bending angle, the channeling deflection efficiency is found to be ()%. In this case the ratio of the optimal crystal bending radius to the critical bending radius, , is 2.85 (3.33) for the () crystal. The enhanced channeling efficiency of compared to is experimentally demonstrated [67-69], and it explains the possibility for to use shorter crystals with larger bending angles.

The Geant4 toolkit has also been used to evaluate the deflection efficiency in long crystals with large bending angle for negative particles in the energy regime. The efficiency is spoiled by the crystal length, resulting in negligibly small efficiencies even at low momentum, as shown in Fig. 1 for the case of antiprotons traversing a 7 long crystal with a 14 bending. For a 5 long crystal with 5 bending the efficiency amounts to less than 0.1%.

Detector occupancy

The interaction of the impinging protons on the fixed target might represent a challenge for the detector operations, radiation hardness, and event reconstruction.

The fluence depends mainly on the average number of primary and secondary interactions taking place in the fixed target. The average number of primary interactions can be determined as , where p is the probability for a proton to interact in the target material, F the rate of impinging protons, and f the LHC bunch collision frequency,  MHz assuming 2400 bunches. The probability p can be estimated itself as , where T is the thickness of the target material and its nuclear interaction length. For the target, and results in , which in turn gives for a flux protons/s. Similarly, for the () crystal with as given in Table 2 and using we obtain . Summing together the and () contributions gives . Similar results are obtained evaluating the number of primary interactions as , where is the Avogadro number, the target density, the atomic mass, the total cross section at , and the number of participant nucleons, estimated using the Glauber model for Pb collisions [70] and rescaling to and () assuming a spherical geometry for nuclei. These values are about a factor two smaller than for nominal collisions for the LHCb upgrade, , with   and a total cross section at of .

Secondary interactions have been studied implementing the geometry of the target and crystal of scenario in the Geant4 framework. Interactions of protons with the target and the crystal are generated using Epos  [71], tuned to the corresponding average number of primary interactions, whereas events use Pythia and EvtGen  [72] to describe their production and decay. Figure 11 illustrates the radiography of the device based on the distribution of stable charged particles. A figure of merit of the occupancy can be obtained from the average number of charged particles within the detector acceptance, using the simplified geometrical model discussed in Sect. 5.1. This value is found to be similar to the corresponding average in the case of nominal collisions for the LHCb upgrade.

Fig. 11

Radiography of the target and crystal geometries in Geant4, shown in the zy view. The distribution represents the origin vertex for stable charged particles from different physics processes (Compton, rays, hadronic interactions and pair production). The presence of the decay vertex after the crystal is clearly visible

Characterization of signal events

A sketch of a signal event is shown in Fig. 12. Particles from fixed-target collisions are emitted in a cone with opening angle around the proton direction. Only particles entering the crystal with momenta parallel to the atomic planes within the Lindhard critical angle of few , , are channeled and bent inside the detector acceptance. From the distribution of the polar angle versus the momentum, shown in Fig. 9, the fraction of channeled particles is estimated , and represents the hard component of the momentum spectrum. Most of the particles produced in the target are not channeled, remaining confined inside the beam pipe. A small tilt of few Lindhard angles of the crystal atomic planes with respect to the incoming proton direction will avoid channeling of non-interacting protons with no cost in trapping efficiency.

Fig. 12

Sketch of a signal event: the baryon, produced in the target at the primary vertex () of the event, is channeled and deflected at an angle determined by the crystal curvature . The secondary vertex () of the baryon decay is outside the crystal, and its position is determined by the vertex of the decay products

Baryon decay products are also highly collimated, and might be difficult to reconstruct by the detector. For decays, Pythia simulations indicate that the average angular separation between all pair of tracks range from 1.9 to 2.7, depending on their masses. For the scenario, this results in an average radial separation at the first detection layer of about 1.7, well above the 55 pixel pitch [61]. This would limit the target position a few upstream the first sensor in scenario.

The high momentum, the polar angle, and the invariant mass of the outcoming baryons define a distinct signature of signal events. By applying a few selection criteria based on these kinematic variables, e.g. the angle (invariant mass) to be within a few of (nominal baryon mass) and the momentum , it would be possible to obtain a high background rejection while retaining most of signal events. For this purpose, the LHCb detector performance are crucial: in particular the estimated track angle resolution of 25, combined with good momentum and mass resolutions,  [60] and  [73], respectively. Particle identification is very limited at momentum regime of several hundreds , and is neglected in this study. We assign the particle mass hypothesis based on momentum hierarchy, e.g. for , the highest momentum track is assigned to be the , the second the and the third the .

With such criteria the vast majority of the signal candidates are produced in the target and decay after the crystal, featuring maximal spin precession. This is because particles produced inside the crystal have lower probability to be channeled, and particles decaying inside the crystal are less bent and have smaller polar angles. The reconstructed invariant mass helps to reject backgrounds from other channeled hadrons having a similar topology, e.g. the more abundant and decays for signal events. Information on the primary and decay vertices could also be exploited to further reject baryons either produced at the beginning or decayed towards the end of the crystal, which might induce a small bias on the spin precession.

Before channeling, a net baryon transverse momentum is needed to define the production plane and have non-zero initial transverse polarization , see Fig. 3. This, in turn, requires a non-zero polar angle , for which there are no restrictions and has a distribution for the hard momentum component similar to that shown in Fig. 9. To prevent initial polarization dilution, the baryon spin rotation has to be determined in bins of .

Expected sensitivities

The main contribution to the statistical uncertainty on the d and g factors of baryons, and similarly for all other baryons, can be estimated in the limit of aswhere is the number of channeled and reconstructed baryons. These estimates assume negligibly small uncertainties on , and , and follow directly from Eqs. (7) and (8). An alternative approach to assess the sensitivity is to generate and fit pseudo-experiments using a probability density function based on the spin precession motion and angular distribution, Eqs. (7), (8) and (26). The two methods provide consistent results, although the former tends to underestimate the uncertainties by about a factor two compared to the latter.

The number of baryons, and similarly for all other baryons, channeled in the bent crystal and reconstructed by the detector can be estimated aswhere is the branching fraction of the decay into the final state f, is the fraction of channeled baryons in the crystal, is a “decay flight” efficiency that accounts for the fraction of channeled baryons decaying after the crystal and within the detector fiducial volume, and is the detector reconstruction efficiency for decays. The number of baryons produced with 7 protons on a fixed target can be estimated aswhere F is the proton rate, t the data taking time, S the beam transverse area, the number of target nucleons, the cross-section for production in interactions at center-of-mass energy, and is the antibaryon-to-baryon ratio for the case of antibaryon production. The number of target nucleons is , where is the Avogadro number, (T) is the target density (thickness), () is the atomic mass (mass number), and is a nuclear modification factor taking into account that the number of participant nucleons of the target nuclei, , differs from due to nuclear matter effects. For hard processes in the absence of strong final-state interaction scales with  [70], thus we take . For the tungsten target  g/cm, g/mol, , and .

All the necessary inputs and their values as used for the sensitivity study, summarised in Tables 3 and 4, are taken from a combination of measurements, estimates and Monte Carlo simulations, and are discussed in detail in the following, along with the final results.

Table 3

Production cross sections , initial polarizations , nuclear modification factors , and antibaryon-to-baryon ratios for collisions at , along with the anomalous magnetic moment , decay channels, branching ratios and decay asymmetry parameters , for the different charm, beauty and strange charged baryons. For comparison purposes, the case has been considered in the and final states. Other quantities like particle masses, spins and lifetimes are taken from Ref. [4]. For and antibaryons, includes the fragmentation fraction from quarks, and is the total beauty cross section

Particle	Λc+		Ξc+	Ξ¯b+	Ω¯b+	Ξ¯^+	Ω¯^+
Decay channel	Δ^++K^-	Λπ^+	Δ^++K^-	Ξ¯^+J/ψ	Ω¯^+J/ψ	Λ¯π^+	Λ¯K^+
				Ξ¯c0π^+	Ω¯c0π^+
Cross section, σ [mb]	0.0182		0.0129	4.67×10^-3	4.67×10^-3	3.4	1.03
|s_0|	0.6		0.6	0.6	0.6	0.5	0.2
g^′	-0.3		-0.3	1.4	5.8	1.9	2.2
B	1.09%	0.83%	0.31%	2.9×10^-6	8.3×10^-7	63.83%	43.32%
α_f	-0.67	-0.91	-0.67	0.91	0.91	0.458	-0.642
R_pT	≈1
R_q¯/q	1		1	0.5	0.5	0.8	0.9

Table 4

Channeling, survival, decay flight and detector efficiencies, along with the average energy and squared transverse momentum of channeled baryons (before decay flight requirements), for a target with a or bent crystal in the scenario. Note that already includes . The sensitivity study based on pseudo-experiments makes use of the complete energy spectrum after channeling, from which and reported here are obtained. All estimates, except (see text), are obtained from samples of charm and strange baryons generated separately for each baryon type from 7 proton beam collisions on protons at rest using Pythia. For beauty baryons, we scale the energy of other simulated baryons to obtain an average energy shift estimated assuming a linear dependence with the baryon mass difference

Particle	Λc+		Ξc+	Ξ¯b+	Ω¯b+	Ξ¯^+	Ω¯^+
Decay channel	Δ^++K^-	Λπ^+	Δ^++K^-	Ξ¯^+J/ψ	Ω¯^+J/ψ	Λ¯π^+	Λ¯K^+
				Ξ¯c0π^+	Ω¯c0π^+
	Si
ε_CH [×10^-4]	1.24 [4.14]		1.04 [3.90]	2.09 [8.91]	2.11 [9.10]	1.75 [5.57]	1.44 [3.84]
E¯ [TeV]	1.36 [2.70]		1.24 [2.40]	1.24 [2.44]	1.24 [2.48]	1.12 [1.54]	1.09 [1.33]
p⊥2¯ [GeV^2/c^2]	1.22 [0.75]		1.09 [1.19]	1.55 [1.25]	1.49 [1.25]	0.20 [0.21]	0.34 [0.32]
ε_s [%]	9.9 [6.9]		31.7 [24.9]	46.3 [41.0]	45.1 [40.1]	99.8 [99.8]	99.5 [99.3]
ε_DF [%]	9.9 [6.9]	0.42 [0.16]	31.7 [24.7]	46.3 [39.5]	45.0 [38.7]	0.08 [0.05]	0.20 [0.15]
	Ge
ε_CH [×10^-4]	2.32 [5.57]		2.06 [5.18]	3.92 [11.34]	3.98 [11.63]	3.18 [7.34]	2.57 [5.17]
E¯ [TeV]	1.37 [2.26]		1.30 [2.07]	1.31 [2.16]	1.32 [2.18]	1.19 [1.51]	1.14 [1.30]
p⊥2¯ [GeV^2/c^2]	1.16 [1.05]		1.47 [1.09]	1.51 [1.32]	1.52 [1.33]	0.22 [0.22]	0.35 [0.33]
ε_s [%]	20.0 [17.4]		44.9 [40.9]	59.0 [57.4]	57.9 [56.5]	99.9 [99.9]	99.7 [99.6]
ε_DF [%]	20.0 [17.4]	0.85 [0.52]	44.9 [40.7]	58.9 [56.0]	57.8 [55.2]	0.08 [0.06]	0.20 [0.16]
ε_det [%]	20	10	20	12	12	10	10

Baryon and antibaryon production yields

The and baryon cross sections can be estimated from the total charm production cross section measured by the PHENIX experiment in collisions at  [74] rescaled to , assuming a linear dependence on , and the corresponding fragmentation fractions. For the case the fragmentation fraction is derived from [74], consistent with theoretical predictions [75]. The baryon branching fractions to and final states are taken from Ref. [4]. The fragmentation fraction is estimated considering that all known -hadron fractions, which amount to about 92%, leave room for the unknown , and fractions [76, 77]. Assuming , we obtain , which is used to rescale by a factor the cross section. The absolute branching fraction is estimated from , measured relative to that of , considering that all known decay modes sum to the total width and assuming that the relative resonant contribution to the final state is the same in and decays. No other quasi-two body or decays to the final state are considered for this study. However, there are additional contributions, e.g. and , with similar branching fractions [4, 78] that can be exploited to improve the sensitivity.

For and baryons produced from 7 protons impinging on fixed target, the total beauty cross section can be estimated by rescaling the cross section measured at  [79]. As a working hypothesis the ratios for bottom baryons are assumed to be , on the basis of the results for charm hadron production at lower energies [80-82]. Branching fractions for baryons are known for very few final states. Two suitable two-body decays, requiring a simple two-body angular analysis, are considered. Firstly, the decay, where the and can be detected in the dimuon final state and as a positive track, respectively. This decay has been measured and its branching fraction times the fragmentation function is  [4]. Secondly, the decay, where the charm antibaryon can be reconstructed in the , or final states. This decay has not been observed but its branching fraction can be estimated by comparing the efficiency corrected signal yields for and decays [83]. The average fragmentation fraction  [84], the measured  [4], and the branching ratios, estimated similarly to the previous case, are used for this calculation. Summing together the contributions of the and the decays, we obtain a global branching fraction times the fragmentation function as shown in Table 3. Similar decays can be considered for the baryon. In this case, the and branching ratios are unknown, and we assume the latter to be the same as for decays and scale by , from the ratio between and  [4].

For the and antibaryons, which contain two and three valence antiquarks respectively, the cross sections are estimated by scaling the production cross section using the universal strangeness suppression factor at high energies,  [85]. In turn, the production cross section is estimated from the inclusive cross section measured at beam momenta of 158  [86] and 405  [87] ( and 27.6, respectively). The ratios are taken to be 0.8 and 0.9 for /, /, respectively, as inferred from Au+Au collisions at  [88]. All branching ratios are in this case known [4].

Efficiencies

The channeling efficiency in and crystals includes both the trapping efficiency and deflection efficiency , and has been estimated separately for each baryon type following the procedure described in Sect. 5.2. The trapping efficiency itself accounts for the angular and momentum divergence of the baryons produced in the target, and is evaluated from the fraction of baryons within the Lindhard angle and momentun . Crystal parameters, optimized for charm baryons, are taken to be common for all baryon species.

The decay flight efficiency has two contributions: the survival efficiency, , which accounts for the fraction of channeled baryons decaying after the crystal, and the probability for long-lived baryons to decay within the VELO region, downstream of the nominal collision point. When one of the baryon decay products is a long-lived , also accommodates the probability of the to decay before the large-area tracking system upstream the magnet, downstream of the collision point, assuming it takes on average half of the initial baryon momentum. For simplicity, the same requirements are applied for both and scenarios.

The detector efficiency can be estimated from the product of the geometrical, trigger and tracking efficiencies, the latter including combinatorics and selection efficiencies. The software-based trigger for the LHCb upgrade detector [64], our scenario, is expected to have efficiency for charm hadrons comparable to the current high level trigger [60], , and similarly for other baryons. A specific trigger scheme for the fixed-target experiment based on the distinct signature of the signal events can enhance the trigger efficiency to . The tracking efficiency is estimated to be 70% per track. Following the discussion in Sect. 5, the geometrical efficiency is taken . For decays to final states including baryons we further apply a penalty factor 1 / 2 to account for the additional inefficiencies to reconstruct highly displaced vertices. Note that the inefficiency due to the long lifetime of the baryon is separately taken into account in , as discussed before.

Spin polarization of baryons

The asymmetry parameter of the decay has been measured to be  [4]. For other decays no measurements are available but an effective parameter can be calculated from a Dalitz plot analysis of decays [78], e.g.  [17]. Eventually, a Dalitz plot analysis of decays would provide the ultimate sensitivity to dipole moments. For the decay the asymmetry parameter is taken to be similar to the decay to the same final state, whereas for all beauty antibaryon decays it is assumed to be about the same as for the decay. For the decay the asymmetry parameter is taken from [4]. The decay is predominantly parity-conserving and thus has a negligibly small asymmetry parameter [89]. The polarization can be determined in this case from the angular distribution of the antiproton from the decay, as the polarization can be related to the polarization of its child baryon such that  [90, 91].

The initial polarization of particles produced from the interaction of 7 protons on a fixed target is unknown. However, a measurement from interaction of 230 GeV/c on copper target yields for transverse momentum larger than 1.1 GeV/c [92]. Moreover, from data produced in the interactions of 500  on five thin target foils (one platinum, four diamond) [78], the polarization of the is measured as a function of the transverse momentum. The average polarization is about , reaching for between about 1.24 and 5.2. Considering these measurements and the average transverse momentum of channeled baryons given in Table 4, we assume for both and baryons [54]. The same polarization is assumed for the and antibaryons. Similarly, the initial polarization of and antibaryons produced from the interaction of 7 protons on fixed target is unknown. From proton production below the region [90, 91, 93–96], the are found to be polarized with the same sign and magnitude as the , increasing about linearly with momentum and reaching at 250, whereas the is consistent with no polarization. As a working hypothesis initial polarizations of and 0.2 are assumed for and , respectively, considering the large momentum of channeled antibaryons /c.

Theoretical predictions of for the and baryons range between and 0.22 [5, 54], thus a central value is considered. For and antibaryons we take effective quark mass MDM calculations [6]. For all strange baryons under consideration there exist measurements [4].

Results

Combining all parameters, measurements and estimates discussed above and summarized in Tables 3 and 4, we obtain the signal yields, normalized to the incident proton flux F, shown in (top) Fig. 13. These rates procure the expected EDM and MDM sensitivities reported in (middle and bottom) Fig. 13, for both and crystals and the two considered experimental scenarios, with and with . Germanium crystals provide in all cases significantly better EDM (MDM) sensitivities, which are for scenario of order , and   (, and  ) for charm, beauty and strange baryons, respectively. Here is the nuclear magneton, and the proton mass. Sensivities for scenario would improve by about one order of magnitude.

Fig. 13

(Top) Estimated yields of channeled and reconstructed signal baryons per incident proton on target, and (middle) EDM and (bottom) MDM sensitivities, for and with crystal parameters optimized for and scenarios. A total of and protons on target have been considered for and respectively. For comparison purposes, the case has been studied in the and final states. Blue lines show the sensitivity of the current and MDM measurements

Conclusions

Electric and magnetic dipole moments of short-lived baryons are powerful probes for physics within and beyond the SM. However, EDM and MDM for charm and beauty baryons have not been accessible to date. A unique opportunity to measure at LHC the EDM and MDM of charm, beauty and strange charged baryons has been discussed here. The experimental setup is based on a fixed-target to be installed in the LHC where protons from the beam halo are deflected using a bent crystal, producing transversally polarized baryons from their interactions with the target. A second bent crystal is positioned after the target where charged baryons that are channeled deflect their trajectory and enter the detector acceptance while rotating their spin. The MDM and EDM information can be inferred from the measurement of the spin-polarization vector after the crystal by analysing the angular distribution of the baryon decay products.

The planar channeling efficiency for multi- particles and the spin precession in bent crystals is studied using Geant4 simulations. The main result is that both positive and negative particles feature spin precession and the results agree with predictions based on analytical calculations, also discussed in Ref. [17]. However, planar channeling efficiency for negative particles is consistently lower than for positive particles, thus much higher statistics is required to perform useful measurements. In that case tests based on the MDM for baryons and antibaryons could be performed. The possibility of exploiting axial channeling of negative particles has been briefly discussed but more studies are needed, including Monte Carlo Geant4 simulations, before drawing any conclusion on the possibility to measure electromagnetic dipole moments.

A program of EDM and MDM measurements for , charm baryons, , beauty antibaryons, and , strange antibaryons, has been discussed. The feasibility of the experiment based on the LHCb detector has been assessed relying on both parameterised and Geant4 simulations along with a geometrical model of the detector. Sensitivities for could be reached within a few weeks of dedicated detector operations spanned over several years at a flux of . The possibility of a dedicated experiment a covering larger pseudorapidity region, able to afford higher proton fluxes and longer data taking periods, has also been discussed. For the LHCb layout, optimal bent crystal parameters are determined to be 7 (5) length and 14 (15) bending angle for (), whereas for the dedicated experiment are found to be 12 (7) and 7 (8). In all cases, crystals provide enhanced sensitivity.

This unique physics program would provide important experimental anchor points for QCD calculations and searches for physics beyond the SM. In the case of charm and beauty baryon EDM the limits would be better than current indirect bounds based on the neutron EDM [11, 13, 16], extending the new physics discovery potential of the LHC.