1 Introduction

The experimental searches for the electric dipole moment (EDM) of fundamental particles provide powerful probes for physics beyond the Standard Model (SM). The existence of permanent EDMs requires the violation of parity (P) and time reversal (T) symmetries and thus, relying on the validity of the \(C\!PT\) theorem, the violation of \(C\!P\) symmetry. Since EDM searches started in the fifties [1, 2], there has been an intense experimental program, leading to limits on the EDM of leptons [3,4,5], neutron [6], heavy atoms [7], proton (indirect from \(^{199}\)Hg) [8], and \(\varLambda \) baryon [9]. New experiments are ongoing and others are planned, including those based on storage rings for muon [10, 11], proton and light nuclei [12,13,14]. Comprehensive reviews on EDM experiments can be found in Refs. [15,16,17,18,19,20,21].

The amount of \(C\!P\) violation in the weak interactions of quarks is not sufficient to explain the observed imbalance between matter and antimatter in the Universe. The SM Lagrangian of strong interactions contains a \(C\!P\)-violating term proportional to the QCD vacuum angle \(\theta \); however, no \(C\!P\) violation has been observed in the strong interactions. A stringent upper bound, , is derived from the experimental limit on the EDM of the neutron, \({<} 3.0 \times 10^{-26} e {\mathrm { \,cm}}\) (90% C.L.) [6]. This degree of tuning in the value of \(\theta \) is known as the “strong \(C\!P\) ” problem. Several solutions have been proposed, among which is the Peccei-Quinn mechanism [22,23,24] that predicts the axion as a candidate for dark matter.

EDM searches of fundamental particles rely on the measurement of the spin precession angle induced by the interaction with the electromagnetic field. For unstable particles this is challenging since the precession has to take place before the decay. A solution to this problem requires large samples of high energy polarized particles traversing an intense electromagnetic field.

In this work, we discuss the unique possibility to search for the EDM of the strange \(\varLambda \) baryon and of the charm \({\varLambda } ^+_{c} \) and \({\varXi } ^+_{c} \) baryons at LHC. Using the experimental upper limit of the neutron EDM, the absolute value of the \(\varLambda \) EDM is predicted to be \({<}4.4\times 10^{-26}~e\mathrm { \,cm} \) [25,26,27,28], while the indirect constraints on the charm EDM are weaker,  [29]. Any experimental observation of an EDM would indicate a new source of \(C\!P\) violation from physics beyond the SM. The EDM of the long-lived \(\varLambda \) baryon was measured to be \({<}1.5 \times 10^{-16}~e\mathrm { \,cm} \) (95% C.L.) in a fixed-target experiment at Fermilab [9]. No experimental measurements exist for short-lived charm baryons since negligibly small spin precession would be induced by magnetic fields used in current particle detectors.

By studying the spin precession of polarized \(\varLambda \) baryons, originated from weak charm baryon decays, it is possible to extract the EDM. We show that an improvement of the present limit of about two orders of magnitude is within reach of the LHCb experiment. The measurement of the magnetic dipole moment (MDM) of \(\varLambda \) and \(\overline{\varLambda }\) baryons would allow a test of \(C\!PT\) symmetry at per mille level. A similar test has been performed for the proton [30], electron [31], and muon [32], and a new experiment for the proton is planned [33].

We propose to search for the EDM of short-lived charm baryons produced by interaction of the 7\(\mathrm {\,TeV}\) LHC proton beam on a fixed target and channeled in a bent crystal in front of the LHCb detector. A sizeable spin precession angle for the short-lived \({\varLambda } ^+_{c} \) and \({\varXi } ^+_{c} \) baryons would be possible by exploiting the intense electromagnetic field between crystal atomic planes. The charm baryon decays can be reconstructed using the LHCb detector. From one month dedicated runs, sensitivities at the level of \(10^{-17}~e\mathrm { \,cm} \) can be reached. This research would extend the physics program of the proposed experiment [34, 35] for the measurement of charm baryon MDMs.

2 EDM experiment concept

The magnetic and electric dipole moment of a spin-1/2 particle is given (in Gaussian units) by \(\varvec{\mu } = g \mu _B {\mathbf s}/2\) and \(\varvec{ \delta } = d \mu _B {\mathbf s}/2\), respectively, where \(\mathbf {s}\) is the spin-polarization vectorFootnote 1 and \(\mu _B=e \hbar / (2 m c)\) is the particle magneton, with m its mass. The g and d dimensionless factors are also referred to as the gyromagnetic and gyroelectric ratios. The interaction of magnetic and electric dipole moments with external electromagnetic fields causes the change of the particle spin direction. The experimental setup to measure this effect relies on three main elements: (i) a source of polarized particles whose direction and polarization degree are known; (ii) an intense electromagnetic field able to induce a sizable spin precession angle during the lifetime of the particle; (iii) the detector to measure the final polarization vector by analysing the angular distribution of the particle decays.

2.1 \(\varLambda \) and \(\overline{\varLambda }\) case

A large amount of \(\varLambda \) baryons is produced directly from the LHC \(p\) \(p\) collisions via strong interactions. The initial polarization direction is perpendicular to the production plane, defined by the proton beam and \(\varLambda \) momentum directions, due to parity conservation. The level of polarization increases with the transverse momentum with respect to the beam direction. Thus a significant initial polarization could be achieved by selecting events within specific kinematic regions [36].

In contrast, weak decays of heavy baryons (charm and beauty), mostly produced in the forward/backward directions at LHC, can induce large longitudinal polarization due to parity violation. For example, the decay of unpolarized \({\varLambda } ^+_{c} \) baryons to the \({\varLambda } {{\pi } ^+} \) final state [37], produces \(\varLambda \) baryons with longitudinal polarization \({\approx }-\) \(90\%\), since the decay asymmetry parameter is \(\alpha _{{\varLambda } {{\pi } ^+}} = -0.91 \pm 0.15\) [38]. Another example is the \({{\varLambda } ^0_{b}} \rightarrow {\varLambda } {{J /\psi }} \) decay where \(\varLambda \) baryons are produced almost 100% longitudinally polarized [39, 40].

The spin-polarization vector \(\mathbf s\) of an ensemble of \(\varLambda \) baryons can be analysed through the angular distribution of the \({\varLambda } \rightarrow {p} {{\pi } ^-} \) decay [41, 42],

$$\begin{aligned} \frac{dN}{d\varOmega '} \propto 1 + \alpha \mathbf s \cdot \hat{\mathbf k} ~, \end{aligned}$$
(1)

where \(\alpha = 0.642\, \pm \,0.013\) [38] is the decay asymmetry parameter. The \(C\!P\) invariance in the \(\varLambda \) decay implies \(\alpha = -\overline{\alpha }\), where \(\overline{\alpha }\) is the decay parameter of the charge-conjugate decay. The unit vector \(\hat{\mathbf k} = (\sin \theta ' \cos \phi ' , \sin \theta ' \sin \phi ' , \cos \theta ' )\) indicates the momentum direction of the proton in the \(\varLambda \) helicity frame, with \(\varOmega ' = (\theta ' ,\phi ' )\) the corresponding solid angle, as illustrated in (left) Fig. 1. We can consider the \(\varLambda \) momentum either in the heavy hadron (H) helicity frame, \(\mathrm S_{H}\), shown in (center) Fig. 1, or in the laboratory frame, \(\mathrm S_{\mathrm L}\), defined in (right) Fig. 1. This offers two possible options for the \(\varLambda \) helicity frame, as seen from the \(\mathrm S_{H}\) or the \(\mathrm S_{\mathrm L}\) frames and referred to as \(\mathrm S_{{\varLambda }}\) or \(\mathrm S_{{\varLambda } \mathrm L}\), respectively, the latter sketched in (left) Fig. 1.

Fig. 1
figure 1

Left \(\varLambda \) helicity frame (\(\mathrm S_{{\varLambda } \mathrm L}\)), center heavy baryon (\(\mathrm S_{H}\)), and right laboratory frame (\(\mathrm S_{\mathrm L}\)). The proton and \(\varLambda \) angles, \((\theta ' ,\phi ' )\) and \((\theta ,\phi )\) are defined in the \(\mathrm S_{{\varLambda } \mathrm L}\) and the \(\mathrm S_{H}\) frames, respectively. The z axis in \(\mathrm S_{{\varLambda } \mathrm L}\) is defined by the \(\varLambda \) momentum in \(\mathrm S_{\mathrm L}\), and the x axis is along the normal to the \(\varLambda \) production plane, defined by the \(\varLambda \) and \(H\) momenta in \(\mathrm S_{\mathrm L}\) frame. The z axis in \(\mathrm S_{H}\) is given by the heavy hadron momentum in \(\mathrm S_{\mathrm L}\), and the x axis is parallel to the normal to its production plane. The proton beam momentum is taken along the z axis and the vertical direction by the y axis in the \(\mathrm S_{\mathrm L}\) frame

The dynamics of the spin vector in presence of external electromagnetic fields is given by the T-BMT equation [43,44,45] (see Appendix A). For a neutral particle in a magnetic field \(\mathbf B\) in the laboratory with negligible field gradient effects, the general solution as a function of the \(\varLambda \) flight length l is described in Sect. A.1. For the particular case of the \(\varLambda \) and \(H\) baryons flying along the z axis in \(\mathrm S_{\mathrm L}\) frame, an initial longitudinal polarization \(s_0\), \(\mathrm{i.e.\,} \)  \(\mathbf s_0=(0,0,s_0)\), and \(\mathbf B = (0,B_y,0)\), the solution is

$$\begin{aligned}&\mathbf s = \left\{ \begin{array}{l} s_x = - s_{0} \sin \varPhi \\ s_y = - s_{0} \dfrac{d \beta }{g} \sin \varPhi \\ s_z = s_{0} \cos \varPhi \\ \end{array} \right. \nonumber \\&\text {where~} {\varPhi = \frac{D_y\mu _B}{\beta \hbar c} \sqrt{d^2 \beta ^2 + g^2} \approx \frac{g D_y \mu _B}{\beta \hbar c} }, \end{aligned}$$
(2)

with \(D_y\equiv D_y(l) = \int _0^l B_y dl'\) the integrated magnetic field along the \(\varLambda \) flight path. The polarization vector precesses in the xz plane, normal to the magnetic field, with the precession angle \(\varPhi \) proportional to the gyromagnetic factor of the particle. The presence of an EDM introduces a non-zero \(s_y\) component perpendicular to the precession plane of the MDM, otherwise not present. At LHCb, with a tracking dipole magnet providing an integrated field \(D_y \approx \pm 4~\mathrm {T \mathrm { \,m}}\) [46], the maximum precession angle for particles traversing the entire magnetic field region yields \(\varPhi _\mathrm{max} \approx \pm \pi /4\), and allows to achieve about 70% of the maximum \(s_y\) component. Moreover, a test of \(C\!PT\) symmetry can be performed by comparing the g and \(-\bar{g}\) factors for \(\varLambda \) and \(\overline{\varLambda }\) baryons, respectively, which precess in opposite directions as g and d change sign from particle to antiparticle.

Contrarily to the past fixed-target EDM experiments where the momentum direction in the laboratory frame was fixed and perpendicular to the magnetic field [9, 47], in this case the \(\varLambda \) momentum varies being the particle produced from heavy baryon decays. As a consequence, the polarization vector is not fixed to be perpendicular to the magnetic field and the signature of the EDM becomes the variation of the \(s_y\) component of the polarization vector before and after the magnetic field. To avoid the dilution introduced by the rotation of the \(\varLambda \) production plane, the change of the polarization has to be determined separately for ensembles of \(\varLambda \) baryons with similar initial polarization, selected according to the kinematics of the decay. In particular, the projection of the \(\varLambda \) trajectory in the xy plane in \(\mathrm S_{\mathrm L}\) at the z position of the \(H\) production vertex can be used to select events with similar polarization, as discussed in Sect. A.1.2.

2.2 \({\varLambda } ^+_{c} \) and \({\varXi } ^+_{c} \) case

The \({\varLambda } ^+_{c} \) and the \({\varXi } ^+_{c} \) baryon EDM can be extracted by measuring the precession of the polarization vector of channeled particles in a bent crystal. There, a positively-charged particle channeled between atomic planes moves along a curved path under the action of the intense electric field between crystal planes. In the instantaneous rest frame of the particle the electromagnetic field causes the spin rotation. The signature of the EDM is a polarization component perpendicular to the initial baryon momentum and polarization vector, otherwise not present, similarly to the case of the \(\varLambda \) baryon.

The phenomenon of spin precession of positively-charged particles channeled in a bent crystal was firstly observed by the E761 collaboration, which measured the MDM of the strange \(\varSigma ^+\) baryon [48]. The possibility to measure the MDM of short-lived charm baryons using channeling in bent crystals, in the momentum range of hundreds of \({\mathrm {\,GeV\!/}c}\), is discussed in Refs. [49, 50]. The feasibility of the measurement at LHC energies is studied in Ref. [34] and offers clear advantages with respect to lower beam energies since the estimated number of produced charm baryons that are channeled into the crystal is proportional to \(\gamma ^{3/2}\) where \(\gamma \) is the Lorentz factor of the particles.

Fig. 2
figure 2

Left production plane of the \({\varLambda } ^+_{c} \) baryon defined by the proton and the \({\varLambda } ^+_{c} \) momenta. The initial polarization vector \(\mathbf s_0\) is perpendicular to the production plane, along the y axis, due to parity conservation in strong interactions. Right deflection of the baryon trajectory and spin precession in the yz and xy plane induced by the MDM and the EDM, respectively. The red (dashed) arrows indicate the (magnified) \(s_x\) spin component proportional to the particle EDM. \(\varPhi \) is the MDM precession angle and \(\theta _C\) is the crystal bending angle

Charm baryons produced by interaction of protons on a fixed target, \(\mathrm{e.g.\,} \)tungsten target, are polarized perpendicularly to the production plane due to parity conservation in strong interactions [51]. The production plane xz, shown in (left) Fig. 2, is determined by the proton and the charm baryon momenta; the latter defines the z axis. The initial polarization vector \(\mathbf s_0 =(0,s_0,0)\) is perpendicular to the production plane, along the y axis. To induce spin rotation the crystal is bent in the yz plane.

The intense electric field \(\mathbf E\) between the crystal planes which deflects positively-charged particles, transforms into a strong electromagnetic field \(\mathbf {E^*}\approx \gamma \mathbf {E}\), \(\mathbf {B^*}\approx -\gamma \varvec{\beta }\times \mathbf {E} /c\) in the particle rest frame and induces the spin precession, as it is described in detail in Refs. [52, 53] and illustrated in (right) Fig. 2. The crystal bending angle is defined as \(\theta _C=L/\rho _0\), where L is the length of the crystal arc and \(\rho _0\) the curvature radius. The precession angle \(\varPhi \) is defined as the angle between the polarization vector and the y axis, as shown in (right) Fig. 2. In the limit of large boost with Lorentz factor \(\gamma \gg 1\), the precession angle in the yz plane induced by the MDM is [54]

$$\begin{aligned} \varPhi \approx \frac{g-2}{2}\gamma \theta _C, \end{aligned}$$
(3)

where g is the gyromagnetic factor.

In presence of a non-zero EDM, the spin precession is no longer confined to the yz plane, originating a \(s_x\) component proportional to the particle EDM represented by the red (dashed) arrows in (right) Fig. 2. The integration of the equation of motion in presence of EDM is described in Appendix A, as well as the approximations used to solve the equations analytically. The polarization vector, after channeling through the crystal is

$$\begin{aligned} \mathbf s = \left\{ \begin{array}{l} s_{x} \approx s_0 \dfrac{d}{g-2} (\cos {\varPhi }-1) \\ s_{y} \approx s_{0} \cos \varPhi \\ s_{z} \approx s_{0} \sin \varPhi \end{array} \right. , \end{aligned}$$
(4)

where \(\varPhi \) is given by Eq. (3). The polarization can be determined, as in the case of the \(\varLambda \) EDM described in Sect. 2.1, by studying the angular distribution of the final state particles. The angular distribution for non-channeled particles allows to determine the initial polarization along the y axis, which compared to the final polarization allows to extract the gyromagnetic and gyroelectric factors. The same method applies to both \({\varLambda } ^+_{c} \) and \({\varXi } ^+_{c} \) baryons.

For \({\varLambda } ^+_{c} \) decaying to two-body final states such as \(p{{K} ^{*0}} \), \({\varDelta } ^{++}{{\pi } ^-} \), \({\varLambda } (1520){{\pi } ^+} \) and \({\varLambda } {{\pi } ^+} \), the angular distribution is described by Eq. (1), where \(\alpha \) is a parity violating coefficient depending on the final state, \(\hat{\mathbf k}\) the direction of the final state baryon in the \({\varLambda } ^+_{c} \) helicity frame, and \(\mathbf{s}\) the \({\varLambda } ^+_{c} \) polarization vector. In the case of the \({{\varLambda } ^+_{c}} \rightarrow {\varLambda } {{\pi } ^+} \) decay, the \(\alpha \) parameter is measured to be \(\alpha _{{\varLambda } {{\pi } ^+}} = -0.91\pm 0.15\) [38]. For other \({\varLambda } ^+_{c} \) decays no measurements are available but an effective \(\alpha \) parameter can be calculated from a Dalitz plot analysis of \({{\varLambda } ^+_{c}} \rightarrow {p} {{K} ^-} {{\pi } ^+} \) decays [55], as discussed in Appendix B and summarized in Table 4. Eventually, a Dalitz plot analysis would provide the ultimate sensitivity to the EDM measurement. The initial polarization \(s_0\) of \({\varLambda } ^+_{c} \) particles produced from the interaction of 7\(\mathrm {\,TeV}\) protons on a fixed target has not been measured. However, a measurement of \({\varLambda } ^+_{c} \) polarization from 40–70 \({\mathrm {\,MeV\!/}c}\) neutron on carbon target gives \(s_0=0.5\pm 0.2\) [56], and a measurement from interaction of 230 \({\mathrm {\,MeV\!/}c}\)  \({\pi } ^-\) on copper target yields \(s_0=-0.65^{+0.22}_{-0.18}\) [57].

3 Sensitivity studies

3.1 \(\varLambda \) and \(\overline{\varLambda }\) case

To identify the most copious \(\varLambda \) production channels from heavy baryons, we consider decays containing only charged particles in the final state, with at least one originated from the heavy baryon decay vertex. No other long-living particles besides the \(\varLambda \) baryon, except an intermediate \(\varXi ^-\) baryon decaying into the \({\varLambda } {{\pi } ^-} \) final state, are considered. These conditions are required to reconstruct the production and the decay vertex of the \(\varLambda \) particle and eventually exploit this information in the event reconstruction. The number of \(\varLambda \) particles produced can be estimated as

$$\begin{aligned} N_{\varLambda }= & {} 2 \mathcal {L} \sigma _{{q} {\overline{{q}}}} f({q} \rightarrow H)\mathcal {B} (H \rightarrow {\varLambda } X') \nonumber \\&\times \mathcal {B} ({\varLambda } \rightarrow {p} {{\pi } ^-})\mathcal {B} (X'\rightarrow \mathrm {charged}) , \end{aligned}$$
(5)

where \(\mathcal {L}\) is the total integrated luminosity, \(\sigma _{{q} {\overline{{q}}}} \) (\({q} ={c},{b} \)) are the heavy quark production cross sections from \(p\) \(p\) collisions at \(\sqrt{s}\approx 14\) TeV [58,59,60,61], and f is the fragmentation fraction into the heavy baryon \(H\)  [62,63,64,65]. All branching fractions \(\mathcal {B}\) are taken from Ref. [38], and where they are given relative to other decays all the known decay modes are assumed to sum the total width. In Table 1 the dominant production channels and the estimated yields are summarised. Overall, there are about \(1.5\times 10^{11}\) \(\varLambda \) baryons per \(\text{ fb }^{-1}\) produced directly from heavy baryon decays (referred hereafter as short-lived, or SL events), and \(3.8\times 10^{11}\) from charm baryons decaying through an intermediate \(\varXi ^-\) particle (long-lived, or LL events). The yield of \(\varLambda \) baryons experimentally available can then be evaluated as \(N_{\varLambda } ^\mathrm{reco} = \epsilon _\mathrm{geo} \epsilon _\mathrm{trigger} \epsilon _\mathrm{reco} N_{\varLambda } \), where \(\epsilon _\mathrm{geo}\), \(\epsilon _\mathrm{trigger}\) and \(\epsilon _\mathrm{reco}\) are the geometric, trigger and reconstruction efficiencies of the detector system.

Table 1 Dominant \(\varLambda \) production mechanisms from heavy baryon decays and estimated yields produced per \(\text{ fb }^{-1}\) at \(\sqrt{s}=13\) \(\mathrm {\,TeV}\), shown separately for SL and LL topologies. The \(\varLambda \) baryons from \(\varXi ^-\) decays, produced promptly in the \(p\) \(p\) collisions, are given in terms of the unmeasured production cross section
Fig. 3
figure 3

Left sketch of the simplified geometry of the LHCb tracking system in the yz plane. The crosswise lines represent the angular acceptance. The tracking layers and the limits of the R\(_1\) and R\(_2\) regions are shown as solid and dotted thick lines, respectively. The magnet is divided in three regions by thin dotted lines. A simulated \({{\varLambda } ^+_{c}} \rightarrow {\varLambda } (p {{\pi } ^-}){{\pi } ^+} \) decay with corresponding \({{\pi } ^+} \) (green), \({{\pi } ^-} \) (blue) and p (red) tracks is overlaid. Right decay products from \(\varLambda \) baryons decaying in the last region of the magnet, M\(_3\)

Table 2 Geometric efficiencies (in %) for \(\varLambda \) baryons decaying in different regions of the LHCb detector, for several charm baryon decays produced at \(\sqrt{s}=13\) \(\mathrm {\,TeV}\)

The geometric efficiency for SL topology has been estimated using a Monte Carlo simulation of \({p} {p} \) collisions at \(\sqrt{s}=13\) \(\mathrm {\,TeV}\) and the decay of heavy hadrons, using Pythia [66] and EvtGen [67] standalone toolkits, together with a simplified geometrical model of the LHCb detector [46]. Tracking devices upstream of the dipole magnet (VErtex LOcator and Tracker Turicensis) and downstream the magnet (T stations) are modelled to have rectangular shape. The height and width of the tracking layers along the beam axis are determined by the detector angular acceptance, between 10 and 250 \(\mathrm { \,mrad}\) (300 \(\mathrm { \,mrad}\)) in the vertical (horizontal) direction, as illustrated in (left) Fig. 3. Particle trajectories are approximated by straight lines defined by the momentum directions.

Table 2 summarizes the geometric efficiencies for \(\varLambda \) baryons decaying in different regions of the detector volume, for three different SL topologies. Region R\(_1\) is defined such that the z position of the \(\varLambda \) decay vertex is in the range [0–40] \(\mathrm { \,cm}\) from the collision point and the decay products are within the detector acceptance. Events in the R\(_2\) region have a \(\varLambda \) decay z position in the range [40–800]\(\mathrm { \,cm}\). Charged particles produced together with the \(\varLambda \) baryon are required to be within the VELO and T1–T3, or the VELO and TT acceptances, to insure a precise reconstruction of the \(\varLambda \) origin vertex. Events in the R\(_1\) region provide the measurement of the initial \(\varLambda \) polarization vector; events in the R\(_2\) region allow to determine the polarization as a function of the \(\varLambda \) decay length in the magnetic field region. Among the latter, \(\varLambda \) baryons decaying towards the end of the magnet (M\(_3\) region in Table 2) provide most of the sensitivity to the EDM and MDM. These events are sketched in (right) Fig. 3. The total geometric efficiency for R\(_1\) and R\(_2\) regions is about 16%, with small differences among SL topologies, and about \(2.4\times 10^{10}\) \(\varLambda \) baryons per \(\text{ fb }^{-1}\) can be reconstructed.

Fig. 4
figure 4

Left dependence of the \(\varLambda \) gyroelectric factor uncertainty with the initial polarization for \(N_{\varLambda } ^\mathrm{reco}=10^6\) events, and right as a function of the integrated luminosity assuming reconstruction efficiency of 0.2% and 1%

To assess the EDM sensitivity, pseudo-experiments have been generated using a simplified detector geometry that includes an approximate LHCb magnetic field mapping [46, 68]. The angular distribution and spin dynamics have been simulated using Eq. (1) and the general solution as a function of the \(\varLambda \) flight length described in Sect. A.1, respectively. For this study the initial polarization vector \(\mathbf s_0 = (0,0,s_0)\), with \(s_0\) varying between 20 and 100%, and factors \(g=-1.458\) [38] and \(d=0\), were used. Each generated sample was fitted using an unbinned maximum likelihood method with d, g and \(\mathbf s_0\) (or \(\alpha \mathbf s_0\)) as free parameters. The d-factor uncertainty scales with the number of events \(N_{\varLambda } ^\mathrm{reco}\) and the initial longitudinal polarization \(s_0\) as \(\sigma _d \propto 1/(s_0 \sqrt{N_{\varLambda } ^\mathrm{reco}} )\). The sensitivity saturates at large values of \(s_0\), as shown in (left) Fig. 4, and it partially relaxes the requirements on the initial polarizations. Similarly, (right) Fig. 4 shows the expected sensitivity on the EDM as a function of the integrated luminosity, summing together SL and LL events, assuming global trigger and reconstruction efficiency \(\epsilon _\mathrm{trigger} \epsilon _\mathrm{reco}\) of 1% (improved LHCb software-based trigger and tracking for the upgrade detector [69, 70]) and 0.2% (current detector [46]), where the efficiency estimates are based on a educated guess. An equivalent sensitivity is obtained for the gyromagnetic factor. Therefore, with 8 \(\text{ fb }^{-1}\) a sensitivity \(\sigma _d \approx 1.5\times 10^{-3}\) could be achieved (current detector), to be compared to the present limit, \(1.7\times 10^{-2}\) [9]. With 50 \(\text{ fb }^{-1}\) (upgraded detector) the sensitivity on the gyroelectric factor can reach \({\approx }\, 3\times 10^{-4}\).

The reconstruction of long-lived \(\varLambda \) baryons decaying inside and after the magnet represents a challenge for the LHCb experiment, introducing significant backgrounds and a limited resolution on the measurement of the \(\varLambda \) momentum and decay point. Events can be reconstructed by exploiting the kinematics of exclusive decays and the determination of the production and the decay vertex of the \(\varLambda \). According to simulation studies even with relatively poor resolutions, the EDM and MDM measurements do not degrade significantly.

Table 3 Definitions and estimated values of the relevant quantities for charm baryon EDM and MDM sensitivity studies, for a tungsten (W) target

3.2 \({\varLambda } ^+_{c} \) and \({\varXi } ^+_{c} \) case

We propose to search for charm baryon EDMs in a dedicated fixed-target experiment at the LHC to be installed in front of the LHCb detector, as close as possible to the VELO detector. The target should be attached to the crystal to maximize the yield of short-lived charm baryons to be channeled. The rate of \({\varLambda } ^+_{c} \) baryons produced with 7\(\mathrm {\,TeV}\) protons on a fixed target can be estimated as

$$\begin{aligned} \frac{dN_{{\varLambda } ^+_{c}}}{dt} = \frac{F}{A}\sigma ({p} {p} \rightarrow {{\varLambda } ^+_{c}} X) N_T , \end{aligned}$$
(6)

where F is the proton rate, A the beam transverse area, \(N_T\) the number of target nucleons, and \(\sigma ({p} {p} \rightarrow {{\varLambda } ^+_{c}} X)\) is the cross-section for \({\varLambda } ^+_{c} \) production in \(p \) \(p \) interactions at \(\sqrt{s}=114.6\mathrm {\,GeV} \) center-of-mass energy. The number of target nucleons is \(N_T=N_A\rho A T A_N/A_T \), where \(N_A\) is the Avogadro number, \(\rho \) (T) is the target density (thickness), and \(A_T\) (\(A_N\)) is the atomic mass (atomic mass number). The rate of \({\varLambda } ^+_{c} \) particles channeled in the bent crystal and reconstructed in the LHCb detector is estimated as

$$\begin{aligned} \frac{dN_{{\varLambda } ^+_{c}} ^\mathrm{reco}}{dt} = \frac{dN_{{\varLambda } ^+_{c}}}{dt} \mathcal {B} ({{\varLambda } ^+_{c}} \rightarrow f)\varepsilon _\mathrm{CH} \varepsilon _\mathrm{DF}({{\varLambda } ^+_{c}})\varepsilon _\mathrm{det}, \end{aligned}$$
(7)

where each quantity and the corresponding estimated value is defined in Table 3. A 6.5\(\mathrm {\,TeV}\) proton beam was extracted from the LHC beam halo by channeling protons in bent crystals [71]. A beam with intensity of \(5\times 10^8~\text {proton/s}\), to be directed on a fixed target, is attainable with this technique [72]. An alternative experimental setup to be considered is a target-crystal system positioned in the vacuum pipe of the LHC where collisions with protons of the beam halo can be reached at comparable rates. Both solutions should be studied very accurately to be compliant with machine protection and safety requirements. Recent results from the UA9 collaboration [71], relative to crystal collimation tests, demonstrated that a similar setup is technically viable and can be installed successfully in the LHC. Fixed-target collision events can be recorded in short dedicated runs or in parallel to the \(p \) \(p \) data taking, if the background caused by the insertion of a fixed target in the beam halo is negligible with respect to \(p \) \(p \) collisions. Both solutions have to be studied in detail using ad-hoc simulations.

The \({\varLambda } ^+_{c} \) cross section can be estimated from the total charm production cross section measured by the PHENIX experiment in proton-proton collisions at \(\sqrt{s} = 200\mathrm {\,GeV} \) [73], \(\sigma _{c\overline{c}} = (567 \pm 57_\mathrm{stat.} \pm 193_\mathrm{syst.}) {\mathrm { \,\upmu b}} \), rescaled to \(\sqrt{s} = 114.6 \mathrm {\,GeV} \) assuming a linear dependence on \(\sqrt{s}\). By applying the \({\varLambda } ^+_{c} \) fragmentation function used in Ref. [73], \(\sigma _{{{\varLambda } ^+_{c}}}/\sigma _{c\overline{c}}\approx 5.6\%\), compatible with theoretical predictions [74], the \({\varLambda } ^+_{c} \) cross section is \(\sigma _{{{\varLambda } ^+_{c}}} \approx 18.2 {\mathrm { \,\upmu b}} \).

The channeling efficiency in silicon crystals, including both channeling angular acceptance and dechanneling effects, is estimated to be \(\varepsilon _\mathrm{CH} \approx 10^{-3}\) [75], while the fraction of \({\varLambda } ^+_{c} \) baryons decaying after the crystal is \(\varepsilon _\mathrm{DF}({{\varLambda } ^+_{c}})\approx 19\%\), for \(\gamma = 1000\) and 10\(\mathrm { \,cm}\) crystal length. The geometrical acceptance for \({{\varLambda } ^+_{c}} \rightarrow {p} {{K} ^-} {{\pi } ^+} \) decaying into the LHCb detector is \(\varepsilon _\mathrm{geo} \approx 25\%\) according to simulation studies. For \({\varLambda } ^+_{c} \) to \(\varLambda \) decays, \(\mathrm{e.g.\,} \) \({{\varLambda } ^+_{c}} \rightarrow {\varLambda } (p{{\pi } ^-}){{\pi } ^+} \), the geometrical efficiency is reduced by about a factor 50 since most \(\varLambda \) baryons decay after the detector tracking volume. The LHCb software-based trigger for the upgrade detector [69] is expected to have efficiency for charm hadrons comparable to the current high level trigger [46], \(\mathrm{i.e.\,} \) \(\varepsilon _\mathrm{trigger} \approx 80\%\). A specific trigger scheme for the fixed-target experiment can be adopted to enhance the trigger efficiency for \({\varLambda } ^+_{c} \) decays close to \(100\%\). For example, a trigger based on the energy loss in a instrumented silicon crystal was used in the E761 experiment to enhance the rate of reconstructed channeled \(\varSigma ^+\) baryons [48]. The tracking efficiency is estimated to be \(70\%\) per track, leading to an efficiency \(\varepsilon _\mathrm{track} \approx 34\%\) for a \({\varLambda } ^+_{c} \) decay with three charged particles. The detector reconstruction efficiency, \(\varepsilon _\mathrm{det}= \varepsilon _\mathrm{geo} \varepsilon _\mathrm{trigger} \varepsilon _\mathrm{track} \), is estimated to be

$$\begin{aligned}&\varepsilon _\mathrm{det}({p} {{K} ^-} {{\pi } ^+}) \approx 5.4\times 10^{-2}, \text { for } {{\varLambda } ^+_{c}} \rightarrow pK^-\pi ^+, \\&\varepsilon _\mathrm{det}({\varLambda } {{\pi } ^+}) \approx 1.0\times 10^{-3}, \text { for } {{\varLambda } ^+_{c}} \rightarrow {\varLambda } {{\pi } ^+}. \end{aligned}$$
Fig. 5
figure 5

Dependence of the (left) d and (right) g uncertainties for the \({\varLambda } ^+_{c} \) baryon, reconstructed in the \({\varDelta } ^{++}{{K} ^-} \) final state, with the number of protons on target. One month of data taking corresponds to \(1.3\times 10^{15}\) incident protons (dashed line), according to the estimated quantities listed in Table 3

The initial \({\varLambda } ^+_{c} \) polarization will be eventually measured using non-channeled \({\varLambda } ^+_{c} \) particles. Few \({\varLambda } ^+_{c} \) decay asymmetry parameters are known, the only one relevant for our experiment is that associated to \({{\varLambda } ^+_{c}} \rightarrow \varLambda (p\pi ^-)\pi ^+\), \(\alpha _{\varLambda {{\pi } ^+}} = -0.91\pm 0.15\) [38]. Asymmetry parameters for different \({\varLambda } ^+_{c} \) decays can be measured precisely at LHCb in the future. At present, they can be computed from existing \({{\varLambda } ^+_{c}} \rightarrow {p} {{K} ^-} {{\pi } ^+} \) amplitude analysis results [55] (see Appendix B), yielding \(\alpha _{{\varDelta } ^{++}{{K} ^-}} = -0.67 \pm 0.30\) for the \({{\varLambda } ^+_{c}} \rightarrow {\varDelta } ^{++}{{K} ^-} \) decay.

For the sensitivity studies we assume \(s_0=0.6\) and \((g-2)/2 = 0.3\), according to experimental results and available theoretical predictions, respectively, quoted in Ref. [50]. The d and \(g-2\) values can be derived from Eq. (4) as

$$\begin{aligned} d\approx & {} \frac{(g-2)A_x}{\alpha s_0 \left( \cos \varPhi -1 \right) },\nonumber \\ g-2\approx & {} \frac{2}{\gamma \theta _C}\arccos {\left( \frac{A_y}{\alpha s_0}\right) } \nonumber \\\approx & {} \frac{2}{\gamma \theta _C}\arcsin {\left( \frac{A_z}{\alpha s_0}\right) } , \end{aligned}$$
(8)

where the quantity \(A_{x,y,z}=\alpha s_{x,y,z}\) is measured from a fit to the angular distribution of the decay products. The main contribution to the statistical uncertainty on d and g, in the limit \(\gamma \gg 1\), can be estimated as

$$\begin{aligned}&\sigma _d&\approx \frac{g-2}{\alpha s_0\left( \cos \varPhi -1 \right) }\frac{1}{\sqrt{N_{{\varLambda } ^+_{c}} ^\mathrm{reco}}},\nonumber \\&\sigma _{g}&\approx \frac{2}{\alpha s_0 \gamma \theta _C }\frac{1}{\sqrt{N_{{\varLambda } ^+_{c}} ^\mathrm{reco}}}, \end{aligned}$$
(9)

where \(N_{{\varLambda } ^+_{c}} ^\mathrm{reco}\) is the number of channeled and reconstructed \({\varLambda } ^+_{c} \), as given in Eq. (7), and \(\varPhi \approx 3 \mathrm { \,rad} \) is the precession angle defined in Eq. (3) estimated using the quantities reported in Table 3. The estimate assumes negligibly small uncertainties on \(\theta _C\), \(\gamma \) and the initial \({\varLambda } ^+_{c} \) polarization, \(s_0\), the latter to be measured with large samples of non-channeled \({\varLambda } ^+_{c} \) decays.

Given the estimated quantities reported in Table 3, we obtain

$$\begin{aligned} \frac{dN_{{\varLambda } ^+_{c}} ^\mathrm{reco}}{dt} \approx 5.9 \times 10^{-3}~{{\mathrm{s}^{-1}}} = 21.2~{{\mathrm{h}^{-1}}} \end{aligned}$$

for \({{\varLambda } ^+_{c}} \rightarrow {\varDelta } ^{++}{{K} ^-} \), and

$$\begin{aligned} \frac{dN_{{\varLambda } ^+_{c}} ^\mathrm{reco}}{dt} \approx 8.3 \times 10^{-5}~{{\mathrm{s}^{-1}}} = 0.3~{{\mathrm{h}^{-1}}} \end{aligned}$$

for \({{\varLambda } ^+_{c}} \rightarrow {\varLambda } (p{{\pi } ^-}){{\pi } ^+} \). For reaching a sensitivity of \(\sigma _d=0.01\), corresponding to \(\sigma _g = 2.1\times 10^{-17} e\) cm, we need, inverting Eq. (9), \(5.6\times 10^3\) \({{\varLambda } ^+_{c}} \rightarrow {\varDelta } ^{++}{{K} ^-} \) or \(3.0\times 10^3\) \({{\varLambda } ^+_{c}} \rightarrow {\varLambda } (p{{\pi } ^-}){{\pi } ^+} \) events, recorded during a data taking time t of

$$\begin{aligned} t&\approx 265 {~\mathrm h} \approx 11 {~\mathrm days}, \text { for }{{\varLambda } ^+_{c}} \rightarrow {\varDelta } ^{++}{{K} ^-} , \\ t&\approx 1.0\times 10^{4}{~\mathrm h} \approx 420 {~\mathrm days}, \text { for }{{\varLambda } ^+_{c}} \rightarrow {\varLambda } (p{{\pi } ^-}){{\pi } ^+} . \end{aligned}$$

Therefore, a measurement of \({\varLambda } ^+_{c} \) EDM is feasible in \({\varLambda } ^+_{c} \) quasi two-body decays at LHCb, while it is difficult in \({\varLambda } ^+_{c} \) to \(\varLambda \) final states.

Considering only \({{\varLambda } ^+_{c}} \rightarrow {\varDelta } ^{++}{{K} ^-} \) events, the uncertainties scale as

$$\begin{aligned} \sigma _d\approx & {} 6.1\times 10^{-3}\frac{1}{\sqrt{t(\mathrm{month})}}, \nonumber \\ \sigma _g\approx & {} 4.0\times 10^{-3}\frac{1}{\sqrt{t(\mathrm{month})}}, \end{aligned}$$
(10)

corresponding to

$$\begin{aligned} \sigma _{\delta }\approx & {} 1.3 \times 10^{-17} e \mathrm { \,cm} \frac{1}{\sqrt{t(\mathrm{month})}},\nonumber \\ \sigma _{\mu }\approx & {} 4.2 \times 10^{-27} \text {erg/G} \frac{1}{\sqrt{t(\mathrm{month})}}, \end{aligned}$$
(11)

where the time t of the data taking period is expressed in months. The dependence of the sensitivity to \({\varLambda } ^+_{c} \) EDM and MDM as a function of the number of incident protons on the target is shown in Fig. 5.

Estimating the \({\varXi } ^+_{c} \) baryon production and the absolute \({{\varXi } ^+_{c}} \rightarrow {p} {{K} ^-} {{\pi } ^+} \) branching fraction as described in Sect. 2.1, we obtain the ratio

$$\begin{aligned} \frac{\sigma _{{\varXi } ^+_{c}} \mathcal {B} ({{\varXi } ^+_{c}} \rightarrow pK^-\pi ^+)}{\sigma _{{\varLambda } ^+_{c}} \mathcal {B} ({{\varLambda } ^+_{c}} \rightarrow pK^-\pi ^+)} \approx 18 \%, \end{aligned}$$
(12)

while the fraction of \({\varXi } ^+_{c} \) baryons decaying after the crystal is \(\varepsilon _\mathrm{DF}({{\varXi } ^+_{c}}) \approx 47\%\). Assuming decay asymmetry parameters and initial polarization similar to the \({\varLambda } ^+_{c} \) baryon, the expected statistical uncertainty on the \({\varXi } ^+_{c} \) MDM and EDM is

$$\begin{aligned} \sigma _{\delta }\approx & {} 2.0 \times 10^{-17} e \mathrm { \,cm} \frac{1}{\sqrt{t(\mathrm{month})}},\nonumber \\ \sigma _{\mu }\approx & {} 6.3 \times 10^{-27} \text {erg/G} \frac{1}{\sqrt{t(\mathrm{month})}}. \end{aligned}$$
(13)

The experimental reconstruction of \({\varLambda } ^+_{c} \) and \({\varXi } ^+_{c} \) baryons with energies of about 1 \(\mathrm {\,TeV}\) features reduced invariant mass and vertex resolution with respect to hadrons originated from nominal \({p} {p} \) collisions in LHCb. This is mainly due to the fact that decay products have relatively large momentum of few hundreds \({\mathrm {\,GeV\!/}c}\) and small opening angles. According to preliminary simulations, signal events can be reconstructed with no showstopper for the analysis.

4 Conclusions

The unique possibility to search for the EDM of strange and charm baryons at LHC is discussed, based on the exploitation of large statistics of baryons with large Lorentz boost and polarization. The \(\varLambda \) strange baryons are selected from weak charm baryon decays produced in \(p \) \(p \) collisions at \(\approx \) 14 TeV center-of-mass energy, while \({\varLambda } ^+_{c} \) and \({\varXi } ^+_{c} \) charm baryons are produced in a fixed-target experiment to be installed in the LHC, in front of the LHCb detector. Signal events can be reconstructed using the LHCb detector in both cases. The sensitivity to the EDM and the MDM of the strange and charm baryons arises from the study of the spin precession in intense electromagnetic fields. The long-lived \(\varLambda \) precesses in the magnetic field of the detector tracking system. Short-lived charm baryons are channeled in a bent crystal attached to the target and the intense electric field between atomic planes induces the spin precession. Sensitivities for the \(\varLambda \) EDM at the level of \( 1.3 \times 10^{-18}~e\mathrm { \,cm} \) can be achieved using a data sample corresponding to an integrated luminosity of 50 \(\text{ fb }^{-1}\) to be collected during the LHC Run 3. A test of \(C\!PT\) symmetry can be performed by measuring the MDM of \(\varLambda \) and \(\overline{\varLambda }\) baryons with a precision of about \(4\times 10^{-4}\) on the g factor. The EDM of the \({\varLambda } ^+_{c} \) (\({\varXi } ^+_{c} \)) can be searched for with a sensitivity of \(1.3~(2.0)\times 10^{-17}/\sqrt{t(\mathrm{month})}~e\mathrm { \,cm} \) with dedicated runs or running in synergetic mode with the LHCb experiment, in parallel to \(p \) \(p \) collisions. Both solutions have to be studied in details using ad-hoc simulations. The proposed experiment would allow about two orders of magnitude improvement in the sensitivity for the \(\varLambda \) EDM and the first search for the charm baryon EDM, expanding the search for new physics through the EDM of fundamental particles.