Measurement of matter-antimatter differences in beauty baryon decays

Differences in the behaviour of matter and antimatter have been observed in $K$ and $B$ meson decays, but not yet in any baryon decay. Such differences are associated with the non-invariance of fundamental interactions under the combined charge-conjugation and parity transformations, known as $C\!P$ violation. Using data from the LHCb experiment at the Large Hadron Collider, a search is made for $C\!P$-violating asymmetries in the decay angle distributions of $\Lambda^0_b$ baryons decaying to $p\pi^-\pi^+\pi^-$ and $p\pi^-K^+K^-$ final states. These four-body hadronic decays are a promising place to search for sources of $C\!P$ violation both within and beyond the Standard Model of particle physics. We find evidence for $C\!P$ violation in $\Lambda^0_b$ to $p\pi^-\pi^+\pi^-$ decays with a statistical significance corresponding to 3.3 standard deviations including systematic uncertainties. This represents the first evidence for $C\!P$ violation in the baryon sector.

The asymmetry between matter and antimatter is related to the violation of the CP symmetry (CPV), where C and P are the charge-conjugation and parity operators. CP violation is accommodated in the Standard Model (SM) of particle physics by the Cabibbo-Kobayashi-Maskawa (CKM) mechanism that describes the transitions between up-and down-type quarks [1,2], in which quark decays proceed by the emission of a virtual W boson and where the phases of the couplings change sign between quarks and antiquarks. However, the amount of CPV predicted by the CKM mechanism is not sufficient to explain our matter-dominated Universe [3,4] and other sources of CPV are expected to exist. The initial discovery of CPV was in neutral K meson decays [5], and more recently it has been observed in B 0 [6,7], B + [8][9][10][11], and B 0 s [12] meson decays, but it has never been observed in the decays of any baryon. Decays of the Λ 0 b (bud) baryon to final states consisting of hadrons with no charm quarks are predicted to have non-negligible CP asymmetries in the SM, as large as 20% for certain three-body decay modes [13]. It is important to measure the size and nature of these CP asymmetries in as many decay modes as possible, to determine whether they are consistent with the CKM mechanism or, if not, what extensions to the SM would be required to explain them [14][15][16].
The decay processes studied in this article, Λ 0 b → pπ − π + π − and Λ 0 b → pπ − K + K − , are mediated by the weak interaction and governed mainly by two amplitudes, expected to be of similar magnitude, from different diagrams describing quark-level b → uud transitions, as shown in Fig. 1. Throughout this paper the inclusion of charge-conjugate reactions is implied, unless otherwise indicated. CPV could arise from the interference of two amplitudes with relative phases that differ between particle and antiparticle decays, leading to differences in the Λ 0 b and Λ 0 b decay rates. The main source of this effect in the SM would be the large relative phase (referred to as α in the literature) between the product of the CKM matrix elements V ub V * ud and V tb V * td , which are present in the different diagrams depicted in Fig. 1. Parity violation (PV) is also expected in weak interactions, but has never been observed in Λ 0 b decays.
The two diagrams show the transitions that contribute most strongly to Λ 0 b → pπ − π + π − and Λ 0 b → pπ − K + K − decays. In both cases, a pair of π + π − (K + K − ) is produced by gluon emission from the light quarks (u,d). The difference is in the b quark decay that happens on the left through a virtual W − boson emission ("tree diagram") and on the right as a virtual W − boson emission and absorption together with a gluon emission ("loop diagram"). The magnitudes of the two amplitudes are expected to be comparable, and each is proportional to the product of the CKM matrix elements involved, which are shown in the figure.
To search for CP -violating effects one needs to measure CP -odd observables, which can be done by studying asymmetries in the T operator. This is a unitary operator that reverses both the momentum and spin three-vectors [17,18], and is different from the antiunitary time-reversal operator T [19,20] that also exchanges initial and final states. A non-zero CP -odd observable implies CP violation and similar considerations apply to P -odd observables and parity violation [21]. Furthermore, different values of P -odd observables for a decay and its charge conjugate would imply CPV. In this paper, scalar triple products of final-state particle momenta in the Λ 0 b centre-of-mass frame are studied to search for P -and CP -violating effects in four-body decays. These are defined In the latter case there is an inherent ambiguity in the choice of the pion for h 1 that is resolved by taking that with the larger momentum in the Λ 0 b rest frame, referred to as π fast . The following asymmetries may then be defined [22,23]: where N and N are the numbers of Λ 0 b and Λ 0 b decays. These asymmetries are P -odd and T -odd and so change sign under P or T transformations, i.e.
The P -and CP -violating observables are defined as and a significant deviation from zero would signal PV or CPV, respectively. Searches for CPV with triple-product asymmetries are particularly suited to Λ 0 b fourbody decays to hadrons with no charm quark [24] thanks to the rich resonant substructure, dominated by ∆(1232) ++ → pπ + and ρ(770) 0 → π + π − resonances in the Λ 0 b → pπ − π + π − final state. The observable a T -odd CP is sensitive to the interference of T -even and T -odd amplitudes with different CP -odd ("weak") phases. Unlike the overall asymmetry in the decay rate that is sensitive to the interference of T -even amplitudes, a T -odd CP does not require a non-vanishing difference in the CP -invariant ("strong") phase between the contributing amplitudes [19,25]. The observables A T , A T , a T -odd P and a T -odd CP are, by construction, largely insensitive to particle-antiparticle production asymmetries and detector-induced charge asymmetries [26].
This article describes measurements of the CP -and P -violating asymmetries introduced in equation (3) The asymmetries are measured first for the entire phase space of the decay, integrating over all possible final state configurations, and then in different regions of phase space so as to enhance sensitivity to localised CPV. The analysis is performed using proton-proton collision data collected by the LHCb detector, corresponding to 3.0 fb −1 of integrated luminosity at centre-of-mass energies of 7 and 8 TeV, and exploits the copious production of Λ 0 b baryons at the LHC, which constitutes around 20% of all b hadrons produced [27]. Control samples of Λ 0 b → pK − π + π − and Λ 0 b → Λ + c π − decays, with Λ + c decaying to pK − π + , pπ − π + , and pK − K + final states, are used to optimise the event selection and study systematic effects; the most abundant control sample consists of Λ 0 b → Λ + c (→ pK − π + )π − decays mediated by b → c quark transitions in which no CPV is expected [28]. To avoid introducing biases in the results, all aspects of the analysis, including the selection, phase space regions, and procedure used to determine the statistical significance of the results, were fixed before the data were examined.
The LHCb detector [29,30] is designed to collect data of b-hadron decays produced from proton-proton collisions at the Large Hadron Collider. It instruments a region around the proton beam axis, covering the polar angles between 10 and 250 mrad, where approximately 24% of the b-hadron decays occur [31]. The detector includes a highprecision tracking system with a dipole magnet, providing measurements of the momentum and decay vertex position of particle decays. Different types of charged particles are distinguished using information from two ring-imaging Cherenkov detectors, a calorimeter and a muon system. Simulated samples of Λ 0 b signal modes and control samples are used in this analysis to verify the experimental method and to study certain systematic effects. These simulated events model the experimental conditions in detail, including the proton-proton collision, the decays of the particles, and the response of the detector. The software used is described in Refs. [32,[34][35][36]38]. The online event selection is performed by a trigger system that takes fast decisions about which events to record. It consists of a hardware stage, based on information from the calorimeter and muon systems, followed by a software stage, which applies a full event reconstruction. The software trigger requires Λ 0 b candidates to be consistent with a b-hadron decay topology, with tracks originating from a secondary vertex detached from the primary pp collision point. The mean Λ 0 b lifetime is 1.5 ps [39], which corresponds to a typical flight distance of a few millimetres in LHCb.
The Λ 0 b → pπ − h + h − candidates are formed by combining tracks identified as protons, pions, or kaons that originate from a common vertex. The proton or antiproton identifies the candidate as a Λ 0 b or Λ 0 b . There are backgrounds from b-hadron decays to charm hadrons that are suppressed by reconstructing the appropriate two-or three-body invariant masses, and requiring them to differ from the known charm hadron masses by at least three times the experimental resolution. For the Λ 0 b → Λ + c π − control mode, only the Λ 0 b → ph + h − π − events with reconstructed ph + h − invariant mass between 2.23 and 2.31 GeV/c 2 are retained.
A boosted decision tree (BDT) classifier [40] is constructed from a set of kinematic variables that discriminate between signal and background. The signal and background training samples used for the BDT are derived from the Λ 0 b → pK − π + π − control sample, since its kinematics and topology are similar to the decays under study; background in this sample is subtracted with the sP lot technique [41], a statistical technique to disentangle signal and background contributions. The background training sample consists of candidates that lie far from the signal mass peak, between 5.85 and 6.40 are used to optimise the particle identification criteria for the signal mode with the same final state. For events in which multiple candidates pass all selection criteria for a given mode, one candidate is retained at random and the rest discarded.
Unbinned extended maximum likelihood fits to the pπ − π + π − and the pπ − K + K − invariant mass distributions are shown in Fig. 2. The invariant mass distribution of the Λ 0 b signal is modelled by a Gaussian core with power law tails [42] with the mean and the width of the Gaussian determined from the fit to data. The combinatorial background is modelled by an exponential distribution with the rate parameter extracted from data. All other parameters of the fit model are taken from simulations except the yields. Partially reconstructed Λ 0 b decays are described by an empirical function [43] convolved with a Gaussian function to account for resolution effects. The shapes of backgrounds from other b-hadron decays due to incorrectly identified particles, e.g. kaons identified as pions or protons identified as kaons, are modelled using simulated events. These consist mainly of Λ 0 b → pK − π + π − and B 0 → K + π − π − π + decays for the Λ 0 b → pπ − π + π − sample and of similar final states for the Λ 0 b → pπ − K + K − sample, as shown in Fig. 2. The yields of these contributions are obtained from fits to data reconstructed under the appropriate mass hypotheses for the final-state particles. The signal yields of Λ 0 b → pπ − π + π − and Λ 0 b → pπ − K + K − are 6646 ± 105 and 1030 ± 56, respectively. This is the first observation of these decay modes.
Signal candidates are split into four categories according to Λ 0 b or Λ 0 b flavour and the sign of C T or C T in order to calculate the asymmetries defined in equations (1) and (2). The reconstruction efficiency for signal candidates with C T > 0 is identical to that with C T < 0 within the statistical uncertainties of the control sample, and likewise for C T , which indicates that the detector and the reconstruction program do not bias this measurement. This check is performed both on the Λ 0 b → Λ + c (pK − π + )π − data control sample and on large samples of simulated events, using yields about 30 times those found in data, which are generated with no CP asymmetry. The CP asymmetry measured in the control sample is a T -odd categories, and are found to be uncorrelated. Corresponding asymmetries for each of the background components are also measured in the fit; they are found to be consistent with zero, and do not lead to significant systematic uncertainties in the signal asymmetries. The values of a T -odd CP and a T -odd P are then calculated from A T and A T . In four-body particle decays, the CP asymmetries may vary over the phase space due to resonant contributions or their interference effects, possibly cancelling when integrated over the whole phase space. Therefore, the asymmetries are measured in different regions of phase space for the Λ 0 b → pπ − π + π − decay using two binning schemes, defined before examing the data. Scheme A, defined in Table 1, is designed to isolate regions of phase space according to their dominant resonant contributions. Scheme B exploits in more detail the interference of contributions which could be visible as a function of the angle Φ between the decay planes formed by the pπ − fast and the π − slow π + systems, as illustrated in Fig. 3. Scheme B has 10 non-overlapping bins of width π/10 in |Φ|. For every bin in each of the schemes, the Λ 0 b efficiencies for C T > 0 and C T < 0 are compared and found to be equal within uncertainties, and likewise the Λ 0 b efficiencies for C T > 0 and C T < 0. The analysis technique is validated on the Λ 0 b → Λ + c (pK − π + )π − control sample, for which the angle Φ is defined by the decay planes of the pK − and π + π − pairs, and on simulated signal events.
The asymmetries measured in Λ 0 b → pπ − π + π − decays with these two binning schemes are shown in Fig. 4 and reported in Table 2, together with the integrated measurements. For each scheme individually, the compatibility with the CP -symmetry hypothesis is evaluated by means of a χ 2 test, with χ 2 = R T V −1 R, where R is the array of a T -odd CP measurements and V is the covariance matrix, which is the sum of the statistical and systematic covariance matrices. An average systematic uncertainty, whose evaluation is discussed below, is assigned for all bins. The systematic uncertainties are assumed to be fully correlated; their contribution is small compared to the statistical uncertainties. respectively. The overall significance for CPV in Λ 0 b → pπ − π + π − decays from the results Table 1: Definition of binning scheme A for the decay mode Λ 0 b → pπ − π + π − . Binning scheme A is defined to exploit interference patterns arising from the resonant structure of the decay. Bins 1-4 focus on the region dominated by the ∆(1232) ++ → pπ + resonance. The other eight bins are defined to study regions where pπ − resonances are present (5-8) on either side of the ρ(770) 0 → π + π − resonances (5-12). Further splitting for |Φ| lower or greater than π/2 is done to reduce potential dilution of asymmetries, as suggested in Ref. [19]. Masses are in units of GeV/c 2 .
For the Λ 0 b → pπ − K + K − decays, the smaller purity and signal yield of the sample do not permit PV and CPV to be probed with the same precision as for Λ 0 b → pπ − π + π − , and therefore only two regions of phase space are considered. One spans 1.43 < m(pK − ) < 2.00 GeV/c 2 (bin 1) and is dominated by excited Λ resonances decaying to pK and the other covers the remaining phase space, 2.00 < m(pK − ) < 4.99 GeV/c 2 (bin 2). The observables measured in these regions are given in Table 2 and are consistent with CP and P symmetry.
The main sources of systematic uncertainties for both pπ − π + π − and pπ − K + K − decays are experimental effects that could introduce biases in the measured asymmetries. This is tested by measuring the asymmetry a T -odd CP , integrated over phase space and in various phase space regions, using the control sample Λ 0 b → Λ + c (→ pK − π + )π − , which is expected to exhibit negligible CPV. The results are in agreement with the CP -symmetry hypothesis; an uncertainty of 0.31% is assigned as a systematic uncertainty for the a T -odd CP and a T -odd P integrated measurements; an uncertainty of 0.60%, the largest asymmetry from a fit to scheme B measurements using a range of efficiency and fit models, is assigned for the corresponding phase space measurements. The systematic uncertainty arising from the experimental resolution in the measurement of the triple products C T and C T , which could introduce a migration of events between the bins, is estimated from simulated samples of Λ 0 b → pπ − π + π − and Λ 0 b → pπ − K + K − decays where neither P -nor CP -violating effects are present. The difference between the reconstructed and generated asymmetry is taken as a systematic uncertainty due to this effect and is less than 0.06% in all cases. To assess the uncertainty associated with the fit models, alternative functions are used; these tests lead only to small changes in the asymmetries, the largest being 0.05%. For Λ 0 b → pπ − K + K − decays, this contribution is larger, about 0.28% for the a T -odd CP and a T -odd P asymmetries.
Further cross-checks are made to investigate the stability of the results with respect to different periods of recording data, different polarities of the spectrometer magnet, the choice made in the selection of multiple candidates, and the effect of the trigger and selection criteria. Alternative binning schemes are studied as a cross-check, such as using 8 or 12 bins in |Φ| for Λ 0 b → pπ − π + π − decays. For these alternative binning schemes, the significance of the CPV measurement of the modified scheme B is reduced to below 3σ. Nonetheless, the overall significance of the combination of these two additional binnings with schemes A and B remains above three standard deviations, with a p-value of 1.8 × 10 −3 (3.1σ), consistent with the 3.3σ result seen in the baseline analysis. An independent analysis of the data based on alternative selection criteria confirmed the results. It used a similar number of events, of which 73.4% are in common with the baseline analysis, and gave p-values for CP symmetry of 3.4 × 10 −3 (2.9σ) for scheme A and 1.4 × 10 −4 (3.8σ) for scheme B.
In conclusion, a search for P and CP violation in Λ 0 b → pπ − π + π − and Λ 0 b → pπ − K + K − decays is performed on signal yields of 6646 ± 105 and 1030 ± 56 events. This is the first observation of these decay modes. Measurements of asymmetries in the entire phase space do not show any evidence of P or CP violation. Searches for localised P or CP violation are performed by measuring asymmetries in different regions of the phase space. The results are consistent with CP symmetry for Λ 0 b → pπ − K + K − decays but evidence for CP violation at the 3.3σ level is found in Λ 0 b → pπ − π + π − decays. No significant P violation is found. This represents the first evidence of CP violation in the baryon sector and indicates an asymmetry between baryonic matter and antimatter.
LHCb collaboration c P. N