Abstract
Permanent electric dipole moments (EDMs) of fundamental particles provide powerful probes for physics beyond the Standard Model. We propose to search for the EDM of strange and charm baryons at LHC, extending the ongoing experimental program on the neutron, muon, atoms, molecules and light nuclei. The EDM of strange \(\varLambda \) baryons, selected from weak decays of charm baryons produced in \(p \) \(p \) collisions at LHC, can be determined by studying the spin precession in the magnetic field of the detector tracking system. A test of \(C\!PT\) symmetry can be performed by measuring the magnetic dipole moment of \(\varLambda \) and \(\overline{\varLambda }\) baryons. For short-lived \({\varLambda } ^+_{c} \) and \({\varXi } ^+_{c} \) baryons, to be produced in a fixed-target experiment using the 7 \(\mathrm {\,TeV}\) LHC beam and channeled in a bent crystal, the spin precession is induced by the intense electromagnetic field between crystal atomic planes. The experimental layout based on the LHCb detector and the expected sensitivities in the coming years are discussed.
Similar content being viewed by others
Explore related subjects
Discover the latest articles, news and stories from top researchers in related subjects.Avoid common mistakes on your manuscript.
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],
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.
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
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.
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]
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
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
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.
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.
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.
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
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
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
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
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
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
for \({{\varLambda } ^+_{c}} \rightarrow {\varDelta } ^{++}{{K} ^-} \), and
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
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
corresponding to
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
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
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.
Notes
The spin-polarization vector is defined such as \(\mathbf s = 2 \langle \mathbf S \rangle / \hbar \), where \(\mathbf S\) is the spin operator.
References
E.M. Purcell, N.F. Ramsey, Phys. Rev. 78, 807 (1950). doi:10.1103/PhysRev.78.807
J.H. Smith, E.M. Purcell, N.F. Ramsey, Phys. Rev. 108, 120 (1957). doi:10.1103/PhysRev.108.120
J. Baron et al., Science 343, 269 (2014). doi:10.1126/science.1248213
G.W. Bennett et al., Phys. Rev. D 80, 052008 (2009). doi:10.1103/PhysRevD.80.052008
K. Inami et al., Phys. Lett. B 551, 16 (2003). doi:10.1016/S0370-2693(02)02984-2
J.M. Pendlebury et al., Phys. Rev. D 92(9), 092003 (2015). doi:10.1103/PhysRevD.92.092003
W.C. Griffith, M.D. Swallows, T.H. Loftus, M.V. Romalis, B.R. Heckel, E.N. Fortson, Phys. Rev. Lett. 102, 101601 (2009). doi:10.1103/PhysRevLett.102.101601
V.F. Dmitriev, R.A. Sen’kov, Phys. Rev. Lett. 91, 212303 (2003). doi:10.1103/PhysRevLett.91.212303
L. Pondrom, R. Handler, M. Sheaff, P.T. Cox, J. Dworkin, O.E. Overseth, T. Devlin, L. Schachinger, K.J. Heller, Phys. Rev. D 23, 814 (1981). doi:10.1103/PhysRevD.23.814
J. Grange, et al., Muon (g-2) Technical Design Report. (2015). arXiv:1501.06858
N. Saito, AIP Conf. Proc. 1467, 45 (2012). doi:10.1063/1.4742078
V. Anastassopoulos et al., A storage ring experiment to detect a proton electric dipole moment, Rev. Sci. Instrum. 87(11), 115116 (2016). doi:10.1063/1.4967465
J. Pretz, Phys. Scr. 2015(T166), 014035 (2015). doi:10.1088/0031-8949/2015/T166/014035. http://stacks.iop.org/1402-4896/2015/i=T166/a=014035
I.B. Khriplovich, Phys. Lett. B 444, 98 (1998). doi:10.1016/S0370-2693(98)01353-7
J. Engel, M.J. Ramsey-Musolf, U. van Kolck, Prog. Part. Nucl. Phys. 71, 21 (2013). doi:10.1016/j.ppnp.2013.03.003
T. Fukuyama, Int. J. Mod. Phys. A 27, 1230015 (2012). doi:10.1142/S0217751X12300153
K. Jungmann, Ann. Phys. 525(8–9), 550 (2013). doi:10.1002/andp.201300071
M. Pospelov, A. Ritz, Ann. Phys. 318, 119 (2005). doi:10.1016/j.aop.2005.04.002
Y. Semertzidis, J Phys. Conf. Ser. 335, 012012 (2011). doi:10.1088/1742-6596/335/1/012012
Y.K. Semertzidis, EPJ Web Conf. 118, 01032 (2016). doi:10.1051/epjconf/201611801032
CJG Onderwater, J. Phys. Conf. Ser 295(1), 012008 (2011). doi:10.1088/1742-6596/295/1/012008. http://stacks.iop.org/1742-6596/295/i=1/a=012008
R.D. Peccei, H.R. Quinn, Phys. Rev. Lett. 38, 1440 (1977). doi:10.1103/PhysRevLett.38.1440
S. Weinberg, Phys. Rev. Lett. 40, 223 (1978). doi:10.1103/PhysRevLett.40.223
F. Wilczek, Phys. Rev. Lett. 40, 279 (1978). doi:10.1103/PhysRevLett.40.279
F.K. Guo, U.G. Meissner, JHEP 12, 097 (2012). doi:10.1007/JHEP12(2012)097
D. Atwood, A. Soni, Phys. Lett. B 291, 293 (1992). doi:10.1016/0370-2693(92)91048-E
A. Pich, E. de Rafael, Nucl. Phys. B 367, 313 (1991). doi:10.1016/0550-3213(91)90019-T
B. Borasoy, Phys. Rev. D 61, 114017 (2000). doi:10.1103/PhysRevD.61.114017
F. Sala, JHEP 03, 061 (2014). doi:10.1007/JHEP03(2014)061
J. DiSciacca et al., Phys. Rev. Lett. 110(13), 130801 (2013). doi:10.1103/PhysRevLett.110.130801
R.S. Van Dyck, P.B. Schwinberg, H.G. Dehmelt, Phys. Rev. Lett. 59, 26 (1987). doi:10.1103/PhysRevLett.59.26
G.W. Bennett et al., Phys. Rev. Lett. 92, 161802 (2004). doi:10.1103/PhysRevLett.92.161802
S. Ulmer, Y. Yamazaki, C. Smorra, K. Blaum, K. Franke, Y. Matsuda, H. Nagahama, W. Quint, J. Walz, A. Mooser, G. Schneider, Technical Design Report BASE. Tech. Rep. CERN-SPSC-2013-002. SPSC-TDR-002 (2013). http://cds.cern.ch/record/1503514
V.G. Baryshevsky, Phys. Lett. B 757, 426 (2016). doi:10.1016/j.physletb.2016.04.025
L. Burmistrov, G. Calderini, Y. Ivanov, L. Massacrier, P. Robbe, W. Scandale, A. Stocchi, Measurement of short living Baryon magnetic moment using bent crystals at SPS and LHC. Tech. Rep. CERN-SPSC-2016-030. SPSC-EOI-012, CERN, Geneva (2016). https://cds.cern.ch/record/2194564
K.J. Heller et al., Phys. Rev. Lett. 41, 607 (1978). doi:10.1103/PhysRevLett. 41.607. [Erratum: Phys. Rev. Lett.45(1980)]
J.M. Link et al., Phys. Lett. B 634, 165 (2006). doi:10.1016/j.physletb.2006.01.017
C. Patrignani, Chin. Phys. C 40(10), 100001 (2016). doi:10.1088/1674-1137/40/10/100001
R. Aaij et al., Phys. Lett. B 724, 27 (2013). doi:10.1016/j.physletb.2013.05.041
G. Aad et al., Phys. Rev. D 89(9), 092009 (2014). doi:10.1103/PhysRevD.89.092009
T.D. Lee, C.N. Yang, Phys. Rev. 108, 1645 (1957). doi:10.1103/PhysRev.108.1645
J.D. Richman, An experimenter’s guide to the helicity formalism. Tech. Rep. CALT-68-1148, Calif. Inst. Technol., Pasadena, CA (1984). http://cds.cern.ch/record/153636
L.H. Thomas, Nature 117, 514 (1926). doi:10.1038/117514a0
L.H. Thomas, Philos. Mag. 3, 1 (1927)
V. Bargmann, L. Michel, V.L. Telegdi, Phys. Rev. Lett. 2, 435 (1959). doi:10.1103/PhysRevLett.2.435
R. Aaij et al., Int. J. Mod. Phys. A 30, 1530022 (2015). doi:10.1142/S0217751X15300227
L. Schachinger et al., Phys. Rev. Lett. 41, 1348 (1978). doi:10.1103/PhysRevLett.41.1348
D. Chen et al., Phys. Rev. Lett. 69, 3286 (1992). doi:10.1103/PhysRevLett.69.3286
V.V. Baublis et al., Nucl. Instrum. Methods B 90, 112 (1994). doi:10.1016/0168-583X(94)95524-7
V.M. Samsonov, Nucl. Instrum. Methods B 119, 271 (1996). doi:10.1016/0168-583X(96)00348-5
M. Jacob, G.C. Wick, Ann. Phys. 7, 404 (1959). doi:10.1016/0003-4916(59)90051-X
V.G. Baryshevsky, Spin rotation and depolarization of high-energy particles in crystals at Hadron Collider (LHC) and Future Circular Collider (FCC) energies and the possibility to measure the anomalous magnetic moments of short-lived particles (2015). arXiv:1504.06702
I.J. Kim, Nucl. Phys. B 229, 251 (1983). doi:10.1016/0550-3213(83)90363-2
V.L. Lyuboshits, Sov. J. Nucl. Phys. 31, 509 (1980)
E.M. Aitala et al., Phys. Lett. B 471, 449 (2000). doi:10.1016/S0370-2693(99)01397-0
J. Szwed, Phys. Lett. B 105, 403 (1981). doi:10.1016/0370-2693(81)90788-7
M. Jezabek, K. Rybicki, R. Rylko, Phys. Lett. B 286, 175 (1992). doi:10.1016/0370-2693(92)90177-6
R. Aaij et al., JHEP 03, 159 (2016). doi:10.1007/JHEP03(2016)(159). [Erratum: JHEP09,013(2016)]
M. Cacciari, FONLL heavy quark production. http://www.lpthe.jussieu.fr/~cacciari/fonll/fonllform.html. Accessed 17 May 2016
R. Aaij et al., Phys. Lett. B 694, 209 (2010). doi:10.1016/j.physletb.2010.10.010
R. Aaij et al., JHEP 10, 172 (2015). doi:10.1007/JHEP10(2015)172
M. Lisovyi, A. Verbytskyi, O. Zenaiev, Eur. Phys. J. C 76(7), 397 (2016). doi:10.1140/epjc/s10052-016-4246-y
L. Gladilin, Eur. Phys. J. C 75(1), 19 (2015). doi:10.1140/epjc/s10052-014-3250-3
Y. Amhis et al., Averages of \(b\)-hadron, \(c\)-hadron, and \(\tau \)-lepton properties as of summer 2014 (2014). arXiv:1412.7515
M. Galanti, A. Giammanco, Y. Grossman, Y. Kats, E. Stamou, J. Zupan, JHEP 11, 067 (2015). doi:10.1007/JHEP11(2015)067
T. Sjostrand, S. Mrenna, P.Z. Skands, JHEP 05, 026 (2006). doi:10.1088/1126-6708/2006/05/026
D.J. Lange, Nucl. Instrum. Methods A 462, 152 (2001). doi:10.1016/S0168-9002(01)00089-4
A. Hicheur, G. Conti, in Proceedings, 2007 IEEE Nuclear Science Symposium and Medical Imaging Conference (NSS/MIC 2007): Honolulu, Hawaii, October 28–November 3, 2007, pp. 2439–2443 (2007). doi:10.1109/NSSMIC.2007.4436650
LHCb collaboration, LHCb Trigger and Online Technical Design Report (2014). LHCb-TDR-016. https://cds.cern.ch/record/1701361
LHCb collaboration, LHCb Tracker Upgrade Technical Design Report (2014). LHCb-TDR-015. https://cds.cern.ch/record/1647400
W. Scandale et al., Phys. Lett. B 758, 129 (2016). doi:10.1016/j.physletb.2016.05.004
J.P. Lansberg et al., PoS QNP2012, 049 (2012)
A. Adare et al., Phys. Rev. Lett. 97, 252002 (2006). doi:10.1103/PhysRevLett.97.252002
B.A. Kniehl, G. Kramer, Phys. Rev. D 71, 094013 (2005). doi:10.1103/PhysRevD.71.094013
V.M. Biryukov et al., Crystal Channeling and Its Application at High-Energy Accelerators (Springer, Berlin, 1997)
J.D. Jackson, Classical Electrodynamics (Wiley, New York, 1998)
E. Leader, Spin in Particle Physics, Cambridge Monogr. Part. Phys., Nucl. Phys. Cosmol. 15, 1 (2001)
T. Fukuyama, A.J. Silenko, Int. J. Mod. Phys. A 28, 1350147 (2013). doi:10.1142/S0217751X13501479
A.J. Silenko, Phys. Scr. 90(6), 065303 (2015). doi:10.1088/0031-8949/90/6/065303
R.H. Good, Phys. Rev. 125, 2112 (1962). doi:10.1103/PhysRev.125.2112
E.M. Metodiev, Thomas-BMT equation generalized to electric dipole moments and field gradients (2015). arXiv:1507.04440
D.P. Grosnick et al., Nucl. Instrum. Methods A 290, 269 (1990). doi:10.1016/0168-9002(90)90541-D
Acknowledgements
We express our gratitude to our colleagues of the LHCb collaboration. The authors would like to thank G. Cavoto, M. Ferro-Luzzi, G. Graziani, M. Nebot, M. Schiller, A. Pich, V. Vagnoni and G. Wilkinson for interesting discussions. We acknowledge support from INFN (Italy), MinECo and GVA (Spain).
Author information
Authors and Affiliations
Corresponding authors
Appendices
Appendix A: Spin precession and time evolution equations
The time evolution of the spin-polarization vector for a particle with charge q in an electromagnetic field, as a function of the proper time \(\tau \), is given by the Thomas–Bargmann–Michel–Telegdi (T-BMT) equation [43,44,45],
where \(F^{\mu \nu }\) is the electromagnetic tensor, \(a^\mu = (a^0,\mathbf a)\) is the spin 4-pseudovector, and \(p^\mu = m u ^\mu = \left( E/c,\mathbf p \right) \) is the momentum 4-vector. For homogeneous fields, the velocity derivative is given by the Lorentz force,
In the rest frame of the particle, \(a^\mu =(0, \mathbf s)\), \(p^\mu =(mc,\varvec{0})\), where \(\mathbf s\) is the non-relativistic spin-polarization vector. Therefore, in any frame \(a^\mu p_\mu =0\) and \(a_\mu a^\mu =-{\mathbf s}^2\).
In a frame comoving with respect to the particle rest frame where the particle has velocity \(\varvec{\beta }= \mathbf p / m \gamma \), \(\mathrm{e.g.\,} \)the laboratory frame, \(a^\mu \) is given by [76, 77]
where the components of the momentum 4-vector are \(p^0 = \gamma m c^2\) and \(\mathbf p =\gamma m \varvec{\beta }c\). Substituting in the covariant Eq.(14), the spin precession equation is [76,77,78,79],
where t is the time in the laboratory frame, and the precession angular velocity vector \(\varvec{\Omega }\) has been split into three contributions,
corresponding to the MDM, EDM and Thomas precession. The electric and magnetic fields, \(\mathbf E\) and \(\mathbf B\), respectively, are expressed in the laboratory frame.
For a neutral particle (\(q=0\)) the Thomas precession term, arising from Lorentz forces, does not contribute and we obtain the classical equation, \(d \mathbf s/ d \tau = \varvec{\mu }\times \mathbf B^* + \varvec{\delta }\times \mathbf E^*\), where \(\mathbf E^*\) and \(\mathbf B^*\) are the external fields in the rest frame of the particle [76]. Equations (17) and (18) can be generalized to account for field gradient effects as described in Ref. [80, 81].
1.1 A.1: Spin time evolution for the \(\varLambda \) case
For \(\mathbf E = 0\) and \(q=0\), Eqs. (17) and (18) simplify to
where \(\varvec{\beta }\) is the particle velocity in the laboratory frame. This system of homogeneous first order linear differential equations can be solved analytically with the approximation that the precession of the particle depends only on the integrated magnetic field along its flight path. Given the initial condition \(\mathbf s(0)=\mathbf s_0\), the time evolution of the polarization is
where \(\varOmega = |\varvec{\Omega }|\) and \(\varvec{\omega } = \varvec{\Omega } / \varOmega \), with the precession angular velocity given by Eq. (20).
The polarization in terms of the experimentally measured \(\varLambda \) flight length \(l=\beta c t\), \(\mathbf s (l)\), has similar form,
where \(\varPhi = |\varvec{\Phi }|\) and \(\varvec{\omega }' =\varvec{\Phi }/\varPhi \). The precession angle vector is
with \(\mathbf D \approx \mathbf {\overline{B}} l = \int _0^l \mathbf B(\mathbf r_0 + \varvec{\beta }l'/\beta )dl'\) the integrated magnetic field along the \(\varLambda \) flight path.
1.1.1 A.1.1: Magnetic field gradients
The inhomogeneities of the magnetic field are not expected to introduce significant effects in the spin precession. The spin equation of motion including first-order field gradients is derived in Ref. [81] to be
In LHCb the ratio of the field gradient terms to the homogeneous field ones can be estimated as
with \(\beta \simeq 1\) and \(\gamma \gg 1\), and where \(| \varvec{\nabla } B | = 1.14~\text {T\,m}^{-1}\) and \(B \equiv | \mathbf B | =1~\text {T}\) are the maximum values within the detector acceptance as extracted from the LHCb field mapping [46, 68]. Therefore, this effect is negligibly small at LHCb.
1.1.2 A.1.2: Spin rotations
The variation of the \(\varLambda \) momentum direction in the laboratory frame results in an initial polarization vector which is not fixed to be perpendicular to the magnetic field. The relative orientation of the spin and magnetic field vectors is determined by two rotations. On one hand, the polarization vector from the equation of motion is given in the comoving rest frame reached from the laboratory frame, \(\mathrm S_{\mathrm L}\), by a pure boost. This is usually referred to as canonical frame [77]. However, the analyser, given by Eq. (1), is defined in the particle helicity frame. The two rest frames, canonical and helicity, are related by the rotation between the \(\mathrm S_{\mathrm L}\) and \(\mathrm S_{{\varLambda } \mathrm L}\) frames, defined by the \(\varLambda \) and \(H\) momentum directions in \(\mathrm S_{\mathrm L}\) (see Fig. 1). One the other hand, the choice of the \(\mathrm S_{{\varLambda } \mathrm L}\) frame induces a second rotation of the polarization components with respect to the \(\mathrm S_{{\varLambda }}\) frame, where the \(\varLambda \) longitudinal polarization is maximal, as illustrated in Fig. 6. This is known in the literature as Wick rotation. To avoid dilution effects, the change of the polarization has to be analysed as a function of the kinematics of the decay. For example, a longitudinally polarized \(\varLambda \) with polarization \(s_0\) along z in \(\mathrm S_{{\varLambda }}\) would have a transverse component in \(\mathrm S_{{\varLambda } \mathrm L}\) of magnitude \(s_0 \sin \alpha \), with \(\sin \alpha = (m_{\varLambda }/ m_{H}) (p_{H}^{(\mathrm L)} / p_{\varLambda } ^{\mathrm (L)}) \sin \theta \) [77]. As shown in Fig. 6, the \(\varLambda \) helicity angle \(\theta \) and the spin direction are related to the \(\varLambda \) impact parameter in the laboratory [82]. The relation can be exploited to define ensembles of \(\varLambda \) particles having similar initial polarization, therefore improving the sensitivity to detect the spin change.
For the sensitivity studies, the rotation of the magnetic field into the \(\mathrm S_{{\varLambda } \mathrm L}\) frame and the Wick rotation are neglected. The first is expected to have a negligible impact on our study since \(\varLambda \) baryons have momenta largely along the z axis, and the main component of the magnetic field is along the vertical direction (\(B_y\)), thus mostly perpendicular to the \(\varLambda \) motion. Instead, the effect of the Wick rotation is not relevant when measuring the spin change of ensembles of \(\varLambda \) particles having similar initial polarization.
1.2 A.2: Spin time evolution for the \({\varLambda } ^+_{c} \) and \({\varXi } ^+_{c} \) case
For \(\mathbf B = 0\) and \(q = +1\), Eq. (18) simplifies to
where
is the anomalous magnetic moment for a spin-1/2 particle. Since we are dealing with ultra relativistic \({\varLambda } ^+_{c} \) with \(\gamma \approx 437\) at \(1 \mathrm {\,TeV} \) energy, in first approximation the terms \(\propto 1/\gamma \) are neglected.
We describe the particle trajectory in a bent crystal using radial coordinates [52], as shown in Fig. 7,
where \(\varOmega \) is the revolution frequency for the particle traversing the bent crystal. In our ultra-relativistic case it is well approximated by \(\varOmega \approx c/\rho _0\), where \(\rho _0\) is the crystal curvature radius. The radius of the trajectory as a function of time is
where a, \(\varOmega _k\) and \(\delta \) are the oscillation amplitude, frequency and phase, respectively; a and \(\delta \) depend on the particle energy and incident angle, while \(\varOmega _k\) depends on the crystal potential and particle energy. The radial equilibrium position \(\rho '_0\) differs from the electric potential minimum position \(\rho _0\), due to the centrifugal potential, avoiding periodical cancellations and therefore inducing spin precession [53]. The electric potential in the crystal around the minimum can be approximated as an harmonic potential,
and the corresponding electric field is
where the oscillation frequency of the particle around its equilibrium position \(\rho '_0\) is \(\varOmega _k=\sqrt{kc^2/eW}\) with W being the particle energy. Typical values for the relevant quantities are \(\rho _0 \sim 30 \mathrm { \,m} \), \(\varOmega \approx c/\rho _0 \sim 10^7 {\mathrm { \,Hz}} \), \(a \sim 10^{-10} \mathrm { \,m} \), \(k = 4\times 10^{17} \mathrm {\,eV}/{\mathrm { \,cm}}^2 \) for a Si crystal, yielding \(\varOmega _k \sim 10^{13} {\mathrm { \,Hz}} \) for \(1 \mathrm {\,TeV} \) particles.
Substituting the radial coordinates and applying the ultra-relativistic approximation to Eq. (26) we obtain:
In absence of EDM, \(\mathrm{i.e.\,} \) \(d=0\), the spin precession inside the bent crystal occurs in the yz plane with the following spin time evolution [52],
for the initial condition \({\mathbf s_0} = \left( 0, s_0, 0\right) \) and where \(\omega \approx 2\mu 'E(\rho '_0)/\hbar \) is the precession frequency. The spin precession angle defined in Eq. (3) is \(\varPhi = \omega \overline{t}\), where \(\overline{t}\) is the time needed to traverse the crystal. In presence of a non-zero EDM the spin precession is no longer confined to the yz plane, generating a \(s_x\) spin component otherwise not present,
To derive Eq. (34), EDM effects are assumed to be small compared to the MDM effects, \(\mathrm{i.e.\,} \) \(d \ll (g-2)\), and therefore \(\varOmega _y,\varOmega _z\ll \varOmega _x\). We neglect terms of order \(\dot{\rho }/c\) where
since the second term of Eq. (34) is about \( \dot{\rho }\rho \varOmega /c^2 \sim \dot{\rho }/c \sim 3 \times 10^{-4}\) and the third term is about \(\dot{\rho }^2/c^2 \sim 9 \times 10^{-8}\). We demonstrate that \(\varOmega \ll \omega \) by requiring the electric force to be identical to the centripetal force,
and obtain \(\omega \approx \frac{2\mu '}{\hbar }E(\rho '_0) \sim 10^{10} {\mathrm { \,Hz}} \gg \varOmega \sim 10^7 {\mathrm { \,Hz}} \).
Then, Eq. (35) simplifies as
and the time evolution is
The second integral is negligibly small since \(\varOmega _k \gg \omega \) and its fast oscillation averages the integral to zero. The calculation can be decomposed into two analytically integrable terms proportional to \(\sin (\varOmega _k t')\sin (\omega t')\) and \(\cos (\varOmega _k t')\sin (\omega t')\). Assuming \(\varOmega _k \gg \omega \), the maximum value of this integral is
where \(\xi = \mu 'ka/\hbar e\varOmega _k\lesssim 10^{-2}\) and terms proportional to \(\xi \) were neglected to derive Eq. (33) [52]. Finally we obtain the time evolution of the polarization vector in presence of a non-negligible EDM,
1.2.1 A.2.1: Electric field gradients
The equations describing the particle trajectory and its spin precession in an electromagnetic field, including first-order electromagnetic field gradients, as well as a particle EDM contributions, are derived in [81]. In absence of magnetic fields the spin precession vector \(\mathbf \Omega = \mathbf \Omega _\mathrm{MDM} + \mathbf \Omega _\mathrm{EDM} + \mathbf \Omega _\mathrm{TH}\) is
with unchanged Thomas precession component. Using the harmonic potential approximation we obtain
and employing the values used in this appendix, the ratio of the field gradient terms to the homogeneous field ones is estimated to be
which is negligibly small in the ultra-relativistic regime.
When including electric field gradient effects, in absence of magnetic fields, the particle trajectory equation becomes
where the first term is the Lorentz force and the following two terms are the MDM and EDM contributions. In our experimental setup the initial spin vector is orthogonal to \(\mathbf E\times \varvec{\beta }\), hence the MDM component is negligible. The typical magnitude of the ratio between the EDM electric field gradient term and the Lorentz force contribution is \(\sim d\gamma \times 10^{-3}\) which can be close to 1 for \(\gamma \sim 1000\) only if \(d\sim 1\), \(\mathrm{i.e.\,} \)similar EDM and MDM magnitudes. However, we assume the EDM magnitude to be tiny with respect to the MDM one, as already assumed in the derivation of the spin equation of motion. In case of a large EDM, this term would make the spin precession frequency dependent on the spin direction.
Appendix B: Asymmetry parameter \(\alpha \) for quasi two-body final states in \({{\varLambda } ^+_{c}} \!\rightarrow {p} {{K} ^-} {{\pi } ^+} \) decays
The angular distribution for a spin \(1/2\!\rightarrow 1/2\; 0\) baryon decay is given by Eq. (1). The parameter \(\alpha \) characterizes the parity violation in the decay and determines the sensitivity to the initial polarization. The effective \(\alpha \) parameter for \({{\varLambda } ^+_{c}} \!\rightarrow {{\overline{K}{}} {}^{*0}} ({{K} ^-} {{\pi } ^+}) {p} \), \({{\varLambda } ^+_{c}} \!\rightarrow \varDelta ^{++} \left( {p} {{\pi } ^+} \right) {{K} ^-} \) and \({{\varLambda } ^+_{c}} \!\rightarrow \varLambda (1520)\left( {p} {{K} ^-} \right) {{\pi } ^+} \) quasi two-body decays can be calculated using the results of an amplitude analysis for \({{\varLambda } ^+_{c}} \!\rightarrow {p} {{K} ^-} {{\pi } ^+} \) decays reported in Ref. [55]. The angular distribution for those decays is determined by the helicity amplitudes. A similar angular distribution to Eq. (1) is obtained for the above quasi two-body decays when integrating over all the decay angles, except for the helicity angle of the baryon daughter of the \({\varLambda } ^+_{c} \). The computed \(\alpha \) parameters are listed in Table 4.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
Funded by SCOAP3
About this article
Cite this article
Botella, F.J., Garcia Martin, L.M., Marangotto, D. et al. On the search for the electric dipole moment of strange and charm baryons at LHC. Eur. Phys. J. C 77, 181 (2017). https://doi.org/10.1140/epjc/s10052-017-4679-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1140/epjc/s10052-017-4679-y