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][2][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 10 −13 −10 −12 s. 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 CP T invariance, the vioa e-mail: nicola.neri@mi.infn.it lation of CP 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.
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 7 TeV 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 CP T 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 axialchanneled 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 Λ + c , Ξ + c charm baryons, Ξ + b , Ω + b 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 θ L = √ 2U 0 /( pβc), where U 0 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 TeV 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][22][23][24][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 L d , being much shorter for negative than for positive charges.
The first measurements of L d for ultra-high energy negative particles has been done at CERN using relatively short bent silicon crystal at few hundreds GeV. It has been measured to be few mm [29][30][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 L d [19]. The dechanneling length scales almost proportionally to the particle momentum-velocity pβ, i.e. L d ∝ pβ, and can be calculated as where m e and r e are the mass and the classical radius of the electron, I is the mean ionization energy, d p is the interplanar spacing, a TF is the Thomas-Fermi screening radius, and Z i and γ are the charge number and boost of the incident particle.
The channeling efficiency under harmonic approximation scales proportionally to 1 − η c , with η c = R c /R, where R is the crystal bending radius and R c ∝ pβ 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 (1 − η c ), Differently, the energy dependence of L d 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][41][42][43][44][45][46], electrons [31,[47][48][49][50] and positrons [31] in a range of energies spanning from 855 MeV (electrons at MAMI) to 400 GeV (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][2][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 mm, 1 cm, 7 cm long Si crystal bent along the (110) plane by a 0.1, 1 and 14 mrad bending angle, respectively. The efficiency for positive particles is not spoiled by L d , being L d > 13 cm for all the momenta, but the Fig. 1 Dependence of the channeling efficiency of protons and antiprotons with the particle momentum for 1 mm, 1 cm and 7 cm long Si crystals bent along the (110) plane by a 0.1, 1 and 14 mrad bending angle, respectively. The curves for the anti-proton interacting with the 1 and 7 cm long Si crystals are superimposed unfavourable ratio between the critical radius R c 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].

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 s is regulated by the T-BMT equation where the precession angular velocity vector is composed by three contributions corresponding to the MDM, EDM, and Thomas precession: where E and B 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 Sketch of the deflection of the Λ + c baryon trajectory and spin precession in a bent crystal. The initial polarization vector s 0 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) s x spin-polarization component proportional to the particle EDM. The Φ angle indicates the spin precession due to the MDM units) as μ = J gμ B s and δ = J dμ B s, respectively, where μ B = eh/(2mc) is the particle magneton. 1 The lifetime of baryons with heavy quark constituents is too short, e.g. the Λ + c baryon lifetime is about 10 −13 s, for a standard magnet to induce any detectable effect to the spinpolarization 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 Λ + c 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 W 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 s 0 = (0, s 0 , 0) 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 θ C . The measurement of the MDM of charm baryons using bent crystals has been widely discussed since the 80's [52][53][54]. Lately, the possibility of the measurement at LHC energies has been considered [7][8][9]. Recently, the search for charm baryon EDM using bent crystals at LHC has been proposed [17].
The spin precession of particles channeled in bent crystals was firstly observed by the E761 Collaboration [55]. Using a 800 GeV/c proton beam impinging on a Cu target, Σ + baryons with 375 GeV/c 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, E, which deflects charged particles, transforms into a strong electromagnetic field E * ≈ γ E, B * ≈ −γ β × E/c 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 Φ ≈ 1 rad are required large crystal bending angles θ C ≈ 1−10 mrad, and large Lorentz factors γ ≈ 10 2 −10 3 , 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 by Inside 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 θ C . 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 nonharmonic potential, as shown in Fig. 4. In the following, we demonstrate that in presence of a nonharmonic 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 plane where Ω is the revolution frequency for the particle traversing the bent crystal, and the electric field described by the planar channel potential V (ρ) is Neglecting EDM contributions, the spin evolution resulting from Eqs. (4), (5), (9) and (10) is for the initial condition s 0 = (0, s 0 , 0) and where The radial coordinate ρ is constant up to δρ/ρ = O(Å/m) = 10 −10 , therefore the spin precession depends on t 0 dV /dρ dt . 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 ρ 0 , which is determined solely by the centripetal force f c induced by the bent trajectory, This statement follows by computing for a complete particle oscillation. By changing the integration variable to dρ and dt = dρ/ρ, thenρ is determined by the non-relativistic energy conservation for the radial motion of channeled particles [19] in which M = mγ and W r 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 O(100 eV) m. The integration boundaries ρ 1,2 are chosen to be the particle oscillation limits, in which Finally, the integral can be trivially computed m 2e Summarizing, 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 1 cm with bending angle larger than 1 mrad, 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 axialchanneled 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 V (x, ρ) is In analogy with planar-channeled particles described in Sect. 3.1, the longitudinal velocity is ultra-relativistic. The velocity components orthogonal to the channel are nonrelativistic; their contribution to the spin precession described in Eq. (5), is negligible. The particle velocity inside the bent crystal is therefore simplified as in which cos(Ωt) ≈ 1 and sin(Ωt) ≈ 0 to a very good approximation for a crystal bending angle θ C ≤ 15 mrad. The spin evolves according to the spin precession equation with precession vector following from Eqs. (5), (19) and (20), The main difference with respect to the planar channeling case is the presence of the E x 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 ω = 2μ E(ρ 0 )/h ≈ 10 10 Hz is three orders of magnitude lower than the oscillation frequency of the particle trapped in the channel, Ω k ≈ 10 13 Hz. The two dominant components Ω x ∝ dV /dρ and Ω y ∝ dV /dx, 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 planarlychanneled 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 of in which the integration is carried on half of an oscillation.
Here, μ is taken with (g−2)/2 = −0.3 and the Λ + c mass [4]. The typical electric field magnitude of the axial channel |E| ≈ 4×10 11 eV/ m is estimated as the ratio between the potential depth ≈ 200 eV and the channel width ≈ 5Å for a Ge crystal, with values taken from Ref. [19]. The oscillation frequency is Ω k = kc 2 /eW ≈ 5.4 × 10 13 for W = 1 TeV, where the constant describing the potential curvature k ≈ 3.2 × 10 22 eV/ m 2 for a Ge axial channel is about eight times the value for a Si axial channel, according to Ref. [19].
Neglecting EDM effects (d = 0), Eqs. (21) and (22) show that dV /dx contributes to the spin component s x via ds x /dt = −s z Ω y . Since the s z 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 s z (t) = s z + δs z (t) changes by an amount of order Δ, The integrated effect along the whole crystal is conservatively estimated by multiplying δs x by the number of particle oscillations, in which a crystal length of L = 10 cm 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 s x spin-polarization component, constituting the EDM signature, will be at the order of 10 −2 .
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.

GEANTsimulations 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][2][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 cm long Si crystals bent along the (111) plane were exposed to a Σ + beam with 375 GeV/c momentum. The deflection angle of the two crystals were + 1.649 mrad and − 1.649 mrad, with measured precession angles of − 72 • ± 26 • and + 51 • ± 23 • , respectively. As expected, the spin in the two crystals precesses in opposite directions. The average of experimental values 60 • ± 17 • is consistent with the predicted value of 62 • ± 2 • . 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 GeV/c 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 + 63.3 • ± 0.2 • and − 63.3 • ± 0.2 • was obtained for the two cases, respectively, in good agreement with the predicted values of ±63.0 • . 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 cm long Si crystal bent along the (110) plane by a 1 mrad 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.
, in presence of MDM and EDM, respectively. Figure 6 shows the distribution of the measured trajectory

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 Λ + c decaying to two-body final states such as f = Δ ++ K − , pK * 0 , Λ(1520)π + and Λπ + , the angular distribution is described by Table 1 Average spin precession angle (Φ sim ) and EDM polarization component (s x,sim ) obtained from Geant4 simulation compared to the expected values (Φ exp and s x,exp , respectively), due to the gyromagnetic factor of the particle expressed as g = (g − 2)/2 and the gyroelectric factor d = 5 × 10 −2 , along with the mean channeling deflection efficiency (ε c ), for different 1 TeV/c particles impinging on a 1 cm long Si crystal bent along the (110) plane at 1 mrad angle. The normalization of the polarization vector s 0 has been taken unity, i.e. s 0 = 1 where α f is a parity-violating coefficient depending on the final state f ,k the direction of the final state baryon in the Λ + c helicity frame, Ω the corresponding solid angle, and s the Λ + c polarization vector. Equation (26) holds for all spin-1/2 decays into two-body final states with spins 1/2 + 0, 1/2 + 1 and 3/2 + 0.
The initial polarization s 0 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 s 0 , up to discrete ambiguities, as discussed in Appendix A. These, in turn, can be resolved exploiting the dependence of the angular distribution with the Λ + c 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.
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 mm length and ≈ 100 µrad bending angle; the second crystal should have larger angle, order 10 mrad, 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 0.1 µrad and linear resolution of 5 µm. 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 S1 and S2. 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 b or c quarks at LHC. The detector tracking system consists of three main devices: a vertex locator (VELO) surrounding the p p interaction region and a large-area detector (TT) upstream of the dipole magnet inducing an integrated field of about 4 T m, 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][62][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 S1 scenario the fixed-target setup is positioned outside of the VELO vessel container [61], ≈ 1.16 m upstream of the nominal p p collision point and ≈ 0.87 m upstream of the first VELO module. At this position, a minimal crystal bending angle of about 14 mrad 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 S1 scenario, a beam intensity of 5×10 8 p/ s 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 10 15 protons on target ( PoT). 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.
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 S2 scenario to integrate ∼ 10 17 PoT. 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 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 to provide enhanced vertex resolution. Nevertheless, a minimal bending angle of about 5 mrad is needed to guarantee a good separation between channeled and unchanneled particles. This value can be inferred from the angular distribution of Λ + c baryons produced by 7 TeV 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 ∝ γ −1 ≈ 1 mrad. The corresponding Pythia distribution of polar angle versus momentum, illustrated in Fig. 9, shows that for |θ y | > 5 mrad and momentum higher than ≈ 1 TeV there are practically no unchanneled Λ + c baryons.

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 θ C for Si and Ge 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 Λ + c → pK − π + events, produced in 7 TeV proton beam collisions on protons at rest using Pythia. The Λ + c 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 θ y lies within the (−θ L , θ L ) interval, where θ L is the Lindhard angle, introduced in Sect. 2. For p ≈ 1 TeV in Si (Ge) this angle is about 6 (7) µrad, about three orders of magnitude smaller than the angular divergence of Λ + c baryons shown in Fig. 9. Therefore, the trapping efficiency ε t is largely dominated by the angular opening of the baryons produced in the target. This imposes the crystal to be directly attached to the W target. The per-event deflection efficiency is parameterised as where dechanneling losses inside the bent crystal are described by the factor (1 − η c ) introduced in Sect. 2, which accounts for the shortening of the dechanneling length with respect to a straight crystal, and θ d , which is the ratio of L d in a straight crystal to R c . The per-event critical radius, R c = pc/U (x c ), must be below the crystal curvature radius R = L/θ C , where U (x c ) is the interplanar electric field at the critical transverse coordinate, below which the particle is lost from channeling mode. For Si 110 (Ge 110), x c ≈ 0.885 (0.915) Å and U (x c ) ≈ 5.7 (10) GeV/ cm [19]. A scan in the (L , θ C ) 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 Λ + c 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 Λ + c baryons to be higher than 800 GeV/c. The dependence of the spin polarization is obtained from Eqs. (7) and (8)   shows regions whose uncertainties on d (g) are increased by 20% with respect to the minimum, for Si and Ge in both S1 and S2 scenarios. The gyromagnetic ratio prefers higher L values, as σ g ∝ 1/γ . These wide regions provide the optimal parameters summarized in Table 2, chosen around the minimum of the gyroelectric factor. The parametric approach adopted above to account for dechanneling losses has been validated with Geant4 simulations, discussed in Sects. 2 and 4, using 1 TeVΛ + c baryons with no angular divergence. For a 7 (5) cm long Si (Ge) crystal bent along the (110) plane by a 14 (15) mrad bending angle, the channeling deflection efficiency ε c is found to be 37.8 ± 0.6 (31.7 ± 0.6)%. In this case the ratio of the optimal crystal bending radius to the critical bending radius, R/R c , is 2.85 (3.33) for the Si (Ge) crystal. The enhanced channeling efficiency of Ge compared to Si is experimentally demonstrated [67][68][69], and it explains the possibility for Ge 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 TeV 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 cm long Si crystal with a 14 mrad bending. For a 5 cm long Si crystal with 5 mrad 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 ν = F p/ f , 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, f = 11245 Hz × 2400 = 27 MHz assuming 2400 bunches. The probability p can be estimated itself as p = 1−e −T /λ , where T is the thickness of the target material and λ its nuclear interaction length. For the W target, T = 0.5 cm and λ = 9.95 cm results in p = 0.05, which in turn gives ν = 0.91 for a flux F = 5 × 10 8 protons/s. Similarly, for the Si (Ge) crystal with T = L as given in Table 2 and using λ = 46.52 (26.86) cm we obtain ν = 2.59 (3.15). Summing together the W and Si (Ge) contributions gives ν = 3.49 (4.05). Similar results are obtained evaluating the number of primary interactions as ν = F N A σ pp ρT A part /A T , where N A is the Avogadro number, ρ the target density, A T the atomic mass, σ pp = 48 mb the total p p cross section at √ s ≈ 115 TeV, and A part the number of participant nucleons, estimated using the Glauber model for pPb collisions [70] and rescaling to W and Si (Ge) assuming a spherical geometry for nuclei. These values are about a factor two smaller than for nominal p p collisions for the LHCb upgrade, ν pp = Lσ pp / f = 7.6, with L = 2 × 10 33 cm −2 s −1 and a total p p cross section at √ s = 14 TeV of σ pp = 102.5 mb. Secondary interactions have been studied implementing the geometry of the W target and Ge crystal of S1 scenario in the Geant4 framework. Interactions of protons with the W target and the Ge crystal are generated using Epos [71], tuned to the corresponding average number of primary interactions, whereas Λ + c 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  Fig. 11 Radiography of the W target and Ge 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 Λ + c decay vertex after the crystal is clearly visible Fig. 12 Sketch of a signal event: the Λ + c baryon, produced in the W target at the primary vertex (PV ) of the event, is channeled and deflected at an angle determined by the crystal curvature θ C . The secondary vertex (SV ) of the Λ + c baryon decay is outside the crystal, and its position is determined by the vertex of the decay products similar to the corresponding average in the case of nominal p p collisions for the LHCb upgrade.

Characterization of signal events
A sketch of a Λ + c signal event is shown in Fig. 12. Particles from fixed-target collisions are emitted in a cone with opening angle ≈ 1 mrad around the proton direction. Only particles entering the crystal with momenta parallel to the atomic planes within the Lindhard critical angle of few µrad, θ y ∈ (−θ L , θ L ), are channeled and bent inside the detector acceptance. From the distribution of the polar angle θ y versus the momentum, shown in Fig. 9, the fraction of channeled particles is estimated ≈ 10 −4 , 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.
Baryon decay products are also highly collimated, and might be difficult to reconstruct by the detector. For Λ + c → pK − π + decays, Pythia simulations indicate that the average angular separation between all pair of tracks range from 1.9 to 2.7 mrad, depending on their masses. For the S1 sce-nario, this results in an average radial separation at the first detection layer of about 1.7 mm, well above the 55 µrad pixel pitch [61]. This would limit the target position a few cm upstream the first sensor in S2 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 θ y (invariant mass) to be within a few σ of θ C (nominal baryon mass) and the momentum p 800 GeV/c, 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 µrad, combined with good momentum and mass resolutions, ≈ 1% [60] and ≈ 20 MeV/c 2 [73], respectively. Particle identification is very limited at momentum regime of several hundreds GeV/c, and is neglected in this study. We assign the particle mass hypothesis based on momentum hierarchy, e.g. for Λ + c → pK − π + , the highest momentum track is assigned to be the p, the second the K − and the third the π + .
With such criteria the vast majority of the signal candidates are produced in the W 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 D + → K − π + π + and D + s → K + K − π + decays for Λ + c → pK − π + 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 s 0 , see Fig. 3. This, in turn, requires a non-zero polar angle θ x , 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 θ x .

Expected sensitivities
The main contribution to the statistical uncertainty on the d and g factors of Λ + c baryons, and similarly for all other baryons, can be estimated in the limit of γ 1 as where N reco Λ + c is the number of channeled and reconstructed Λ + c baryons. These estimates assume negligibly small uncertainties on θ C , γ and s 0 , 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 Λ + c baryons, and similarly for all other baryons, channeled in the bent crystal and reconstructed by the detector can be estimated as where B(Λ + c → f ) is the branching fraction of the Λ + c decay into the final state f , ε CH is the fraction of channeled baryons in the crystal, ε DF is a "decay flight" efficiency that accounts for the fraction of channeled Λ + c baryons decaying after the crystal and within the detector fiducial volume, and ε det is the detector reconstruction efficiency for Λ + c → f decays. The number of Λ + c baryons produced with 7 TeV protons on a fixed target can be estimated as where F is the proton rate, t the data taking time, S the beam transverse area, N T the number of target nucleons, σ ( pp → Λ + c X ) the cross-section for Λ + c production in p p interactions at √ s = 114.6 GeV center-of-mass energy, and R q/q is the antibaryon-to-baryon ratio for the case of antibaryon production. The number of target nucleons is where N A is the Avogadro number, ρ (T ) is the target density (thickness), A T (A N ) is the atomic mass (mass number), and R pT is a nuclear modification factor taking into account that the number of participant nucleons of the target nuclei, A part , differs from A N due to nuclear matter effects. For hard processes in the absence of strong final-state interaction A part scales with A N [70], thus we take R pT ≈ 1. For the tungsten target ρ = 19.25 g/cm 3 , A T = 183.84 g/mol, A N = 183.84, and T = 0.5 cm.
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 s 0 , nuclear modification factors R pT , and antibaryon-to-baryon ratios R q/q for pW collisions at √ s ≈ 115 GeV, along with the anomalous magnetic moment g = (g − 2)/2, decay channels, branching ratios B and decay asymmetry parameters α f , for the different charm, beauty and strange charged baryons. For comparison purposes, the Λ + c case has been con-sidered in the Δ ++ K − and Λπ + final states. Other quantities like particle masses, spins and lifetimes are taken from Ref. [4]. For Ξ + b and Ω + b antibaryons, B includes the fragmentation fraction from b quarks, and σ is the total pp → bb beauty cross section   For the Λ + c case the fragmentation fraction f Λ + c ≈ 5.6% is derived from [74], consistent with theoretical predictions [75]. The Λ + c baryon branching fractions to Δ ++ K − and Λπ + final states are taken from Ref. [4]. The Ξ + c fragmentation fraction is estimated considering that all known c-hadron fractions, which amount to about 92%, leave room for the unknown Ξ + c , Ξ 0 c and Ω 0 c fractions [76,77]. Assum- , which is used to rescale by a factor f Ξ + c / f Λ + c ≈ 0.71 the Λ + c cross section. The absolute Ξ + c → Δ ++ K − branching fraction is estimated from B(Ξ + c → pK − π + ), measured relative to that of Ξ + c → Ξ − π + π + , 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 Ξ + c and Λ + c decays. No other quasi-two body Λ + c or Ξ + c decays to the final state pK − π + are considered for this study. However, there are additional contributions, e.g. Λ + c → K * 0 p and Λ + c → Λ(1520)π + , with similar branching fractions [4,78] that can be exploited to improve the sensitivity.
For Ξ + b and Ω + b baryons produced from 7 TeV protons impinging on fixed target, the total beauty cross section can be estimated by rescaling the pp → bb cross section measured at √ s = 7 TeV [79]. As a working hypothesis the ratios R q/q for bottom baryons are assumed to be ≈ 0.5, on the basis of the results for charm hadron production at lower energies [80][81][82]. Branching fractions for Ξ + b baryons are known for very few final states. Two suitable two-body decays, requiring a simple two-body angular analysis, are considered. Firstly, the Ξ + b → Ξ + J/ψ decay, where the J/ψ 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 Ξ + b fragmentation function is ≈ 6 × 10 −7 [4]. Secondly, the Ξ + b → Ξ 0 c π + decay, where the charm antibaryon can be reconstructed in the Ξ + π − π + π − , Ξ + π − or pK + K + π − final states. This decay has not been observed but its branching fraction can be estimated by comparing the efficiency corrected signal yields for [83]. The average fragmentation fraction f Λ 0 b ≈ 7% [84], the mea- [4], and the Ξ 0 c branching ratios, estimated similarly to the previous Ξ + c case, are used for this calculation. Summing together the contributions of the Ξ + b → Ξ + J/ψ and the Ξ + b → Ξ 0 c π + decays, we obtain a global branching fraction times the fragmentation function as shown in Table 3. Similar decays can be considered for the Ω + b baryon. In this case, the Ω + b → Ω 0 c π + and Ω 0 c branching ratios are unknown, and we assume the latter to be the same as for Ξ 0 c decays and scale [4]. For the Ξ + and Ω + antibaryons, which contain two and three s valence antiquarks respectively, the cross sections are estimated by scaling the Λ production cross section using the universal strangeness suppression factor at high energies, λ s ≈ 0.32 [85]. In turn, the Λ production cross section is estimated from the inclusive pp → ΛX cross section mea-sured at beam momenta of 158 GeV [86] and 405 GeV [87] ( √ s ≈ 17.2 and 27.6 GeV, respectively). The ratios R q/q are taken to be 0.8 and 0.9 for Ξ + /Ξ − , Ω + /Ω − , respectively, as inferred from Au+Au collisions at √ s N N = 130 GeV [88]. All branching ratios are in this case known [4].

Efficiencies
The channeling efficiency ε CH in Si and Ge crystals includes both the trapping efficiency ε t and deflection efficiency ε c , 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 800 GeV/c. Crystal parameters, optimized for charm baryons, are taken to be common for all baryon species.
The decay flight efficiency ε DF has two contributions: the survival efficiency, ε s , which accounts for the fraction of channeled baryons decaying after the crystal, and the probability ε l for long-lived baryons to decay within the VELO region, ≈ 80 cm downstream of the nominal p p collision point. When one of the baryon decay products is a long-lived Λ, ε l also accommodates the probability of the Λ to decay before the large-area tracking system upstream the magnet, ≈ 2 m 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 S1 and S2 scenarios.
The detector efficiency ε det 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 S1 scenario, is expected to have efficiency for charm hadrons comparable to the current high level trigger [60], ≈ 80%, 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 ≈ 100%. The tracking efficiency is estimated to be 70% per track. Following the discussion in Sect. 5, the geometrical efficiency is taken ≈ 50%. 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 ε DF , as discussed before.

Spin polarization of baryons
The asymmetry parameter of the Λ + c → Λπ + decay has been measured to be −0.91 ± 0.15 [4]. For other Λ + c decays no measurements are available but an effective α f parameter can be calculated from a Dalitz plot analysis of Λ + c → pK − π + decays [78], e.g. α Λ + c →Δ ++ K − = −0.67 ± 0.30 [17]. Eventually, a Dalitz plot analysis of Λ + c → pK − π + decays would provide the ultimate sensitivity to dipole moments. For the Ξ + c → Δ ++ K − decay the asymmetry parameter is taken to be similar to the Λ + c decay to the same final state, whereas for all beauty antibaryon decays it is assumed to be about the same as for the Λ + c → Λπ + 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 s Ω + = s Λ [90,91]. The initial polarization of Λ + c particles produced from the interaction of 7 TeV protons on a fixed target is unknown. However, a measurement from interaction of 230 GeV/c π − on copper target yields s 0 = −0.65 +0. 22 −0.18 for Λ + c transverse momentum larger than 1.1 GeV/c [92]. Moreover, from data produced in the interactions of 500 GeV/c π − on five thin target foils (one platinum, four diamond) [78], the polarization of the Λ + c is measured as a function of the Λ + c transverse momentum. The average polarization is about −0.1, reaching −0.7 for p 2 ⊥ between about 1.24 and 5.2 GeV 2 /c 2 . Considering these measurements and the average transverse momentum of channeled Λ + c baryons given in Table 4, we assume |s 0 | = 0.6 for both Λ + c and Ξ + c baryons [54]. The same polarization is assumed for the Ξ + b and Ω + b antibaryons. Similarly, the initial polarization of Ξ + and Ω + antibaryons produced from the interaction of 7 TeV protons on fixed target is unknown. From proton production below the TeV region [90,91,[93][94][95][96], the Ξ + are found to be polarized with the same sign and magnitude as the Ξ − , increasing about linearly with momentum and reaching ≈ −0.2 at 250 GeV/c, whereas the Ω − is consistent with no polarization. As a working hypothesis initial polarizations of |s 0 | = 0.5 and 0.2 are assumed for Ξ + and Ω + , respectively, considering the large momentum of channeled antibaryons ≈ 1 TeV/c. Theoretical predictions of g − 2 for the Λ + c and Ξ + c baryons range between −0.64 and 0.22 [5,54], thus a central value g = (g − 2)/2 = −0.3 is considered. For Ξ + b and Ω + b 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, Fig. 13 (Top) Estimated yields of channeled and reconstructed signal baryons per incident proton on target, and (middle) EDM and (bottom) MDM sensitivities, for Si and Ge with crystal parameters optimized for S1 and S2 scenarios. A total of 10 15 and 10 17 protons on target have been considered for S1 and S2 respectively. For comparison purposes, the Λ + c case has been studied in the Δ ++ K − and Λπ + final states. Blue lines show the sensitivity of the current Ξ − and Ω − MDM measurements shown in (top) Fig. 13. These rates procure the expected EDM and MDM sensitivities reported in (middle and bottom) Fig. 13, for both Si and Ge crystals and the two considered experimental scenarios, S1 with 10 15 PoT and S2 with 10 17 PoT. Germanium crystals provide in all cases sig-nificantly better EDM (MDM) sensitivities, which are for S1 scenario of order 10 −17 , 10 −14 and 10 −16 e cm (10 −3 , 10 −1 and 10 −3 μ N ) for charm, beauty and strange baryons, respectively. Here μ N = e /2m p c is the nuclear magneton, and m p the proton mass. Sensivities for S2 scenario would improve by about one order of magnitude.

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-TeV 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 CP T 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 Λ + c , Ξ + c charm baryons, Ξ + b , Ω + b 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 10 15 PoT could be reached within a few weeks of dedicated detector operations spanned over several years at a flux of 5×10 8 p/ s. 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 cm (5 cm) length and 14 mrad (15 mrad) bending angle for Si (Ge), whereas for the dedicated experiment are found to be 12 cm (7 cm) and 7 mrad (8 mrad). In all cases, Ge 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.