CMS-PAS-TOP-23-001 | ||
Probing entanglement in top quark production with the CMS detector | ||
CMS Collaboration | ||
27 March 2024 | ||
Abstract: Entanglement is an intrinsic property of quantum mechanics and its measurement probes the current understanding of the underlying quantum nature of elementary particles at a fundamental level. A measurement of the extent of entanglement in top quark-antiquark ($ \mathrm{t\bar{t}} $) events produced in proton-proton collisions at a center-of-mass energy of 13 TeV is performed on the data recorded by the CMS experiment at the CERN LHC in 2016, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. The events are selected based on the presence of two oppositely charged high transverse momentum leptons. An entanglement-sensitive observable $ D $ is derived from the top quark spin-dependent parts of the $ \mathrm{t\bar{t}} $ production density matrix. Values of $ D $ smaller than $ -$1/3 are evidence of entanglement and, within a particular phase space, $ D $ is measured to be $-$0.478 $ ^{+0.025}_{-0.027} $. With an expected (observed) significance of 5.1 (5.7) standard deviations, this provides observation for quantum mechanical entanglement within $ \mathrm{t\bar{t}} $ pairs in this phase space. This measurement provides a new quantum probe of the inner workings of the standard model at the highest energy ever tested. | ||
Links:
CDS record (PDF) ;
Physics Briefing ;
CADI line (restricted) ;
These preliminary results are superseded in this paper, Submitted to ROPP. The superseded preliminary plots can be found here. |
Figures | |
png pdf |
Figure 1:
Subfigure (a) shows the Feynman diagram for the leading order QCD gluon-gluon fusion process and (b) shows the quark-antiquark fusion process for top quark pair production. |
png pdf |
Figure 2:
Predicted values of $ -(1 + \Delta)/ 3 $ (see Eq. (4)) obtained from calculations of the $ \mathrm{t} \overline{\mathrm{t}} $ production cross section [7], without accounting for detector effects, are shown as a function of $ m({\mathrm{t}\overline{\mathrm{t}}} ) $ and the top quark scattering angle $ \Theta $ for the gluon-gluon fusion production. The black solid lines represent the $ D= -1/ 3$ boundary for entanglement while the black dashed line corresponds to $ m({\mathrm{t}\overline{\mathrm{t}}} )= $ 400 GeV. The minimum value on the z axis corresponds to the boundary below which $ \Delta = tr[\mathbf{C}] - 1$. Top quarks with no spin correlations correspond to a value of 0. |
png pdf |
Figure 3:
Reconstruction-level $ m({\mathrm{t}\overline{\mathrm{t}}} ) $ (left) and $ p_{\mathrm{T}}(\mathrm{t}/\overline{\mathrm{t}}) $ (right) distributions of the combined signal model in the full phase space comparing the modeling of the data by MC when not including $ \eta_{\mathrm{t}} $ (red line in upper ratio pad), no $ p_{\mathrm{T}}(\mathrm{t}/\overline{\mathrm{t}}) $ reweighting is applied (blue line in upper ratio pad), or neither of those (pink line in upper ratio pad). The lower ratio pad compares the data to POWHEGV2+HERWIG (pink line in lower ratio pad), to MG5_AMC@NLO (FXFX)+PYTHIA8 (blue line in lower ratio pad), and finally to the nominal MC including $ \eta_{\mathrm{t}} $ and $ p_{\mathrm{T}}(\mathrm{t}/\overline{\mathrm{t}}) $ reweighting (green line in lower ratio pad). |
png pdf |
Figure 3-a:
Reconstruction-level $ m({\mathrm{t}\overline{\mathrm{t}}} ) $ (left) and $ p_{\mathrm{T}}(\mathrm{t}/\overline{\mathrm{t}}) $ (right) distributions of the combined signal model in the full phase space comparing the modeling of the data by MC when not including $ \eta_{\mathrm{t}} $ (red line in upper ratio pad), no $ p_{\mathrm{T}}(\mathrm{t}/\overline{\mathrm{t}}) $ reweighting is applied (blue line in upper ratio pad), or neither of those (pink line in upper ratio pad). The lower ratio pad compares the data to POWHEGV2+HERWIG (pink line in lower ratio pad), to MG5_AMC@NLO (FXFX)+PYTHIA8 (blue line in lower ratio pad), and finally to the nominal MC including $ \eta_{\mathrm{t}} $ and $ p_{\mathrm{T}}(\mathrm{t}/\overline{\mathrm{t}}) $ reweighting (green line in lower ratio pad). |
png pdf |
Figure 3-b:
Reconstruction-level $ m({\mathrm{t}\overline{\mathrm{t}}} ) $ (left) and $ p_{\mathrm{T}}(\mathrm{t}/\overline{\mathrm{t}}) $ (right) distributions of the combined signal model in the full phase space comparing the modeling of the data by MC when not including $ \eta_{\mathrm{t}} $ (red line in upper ratio pad), no $ p_{\mathrm{T}}(\mathrm{t}/\overline{\mathrm{t}}) $ reweighting is applied (blue line in upper ratio pad), or neither of those (pink line in upper ratio pad). The lower ratio pad compares the data to POWHEGV2+HERWIG (pink line in lower ratio pad), to MG5_AMC@NLO (FXFX)+PYTHIA8 (blue line in lower ratio pad), and finally to the nominal MC including $ \eta_{\mathrm{t}} $ and $ p_{\mathrm{T}}(\mathrm{t}/\overline{\mathrm{t}}) $ reweighting (green line in lower ratio pad). |
png pdf |
Figure 4:
Reconstruction-level distribution (left) of $ \cos \varphi $ with 345 $ < m({\mathrm{t}\overline{\mathrm{t}}} ) < $ 400 GeV and 0.0 $ < \beta < $ 0.9. On the right, the $ p_{\mathrm{T}}(\mathrm{t}/\overline{\mathrm{t}}) $ distributions is shown with the same model comparison as was done in Fig. 3. |
png pdf |
Figure 4-a:
Reconstruction-level distribution (left) of $ \cos \varphi $ with 345 $ < m({\mathrm{t}\overline{\mathrm{t}}} ) < $ 400 GeV and 0.0 $ < \beta < $ 0.9. On the right, the $ p_{\mathrm{T}}(\mathrm{t}/\overline{\mathrm{t}}) $ distributions is shown with the same model comparison as was done in Fig. 3. |
png pdf |
Figure 4-b:
Reconstruction-level distribution (left) of $ \cos \varphi $ with 345 $ < m({\mathrm{t}\overline{\mathrm{t}}} ) < $ 400 GeV and 0.0 $ < \beta < $ 0.9. On the right, the $ p_{\mathrm{T}}(\mathrm{t}/\overline{\mathrm{t}}) $ distributions is shown with the same model comparison as was done in Fig. 3. |
png pdf |
Figure 5:
Reconstruction-level distribution (top) of the combined $ \mathrm{t} \overline{\mathrm{t}} $ + $ \eta_{\mathrm{t}} $ signal model in mixtures of the noSC combined signal sample, more details in the text. Template variations are shown as a function of $ \cos \varphi $ in the phase space of 345 $ < m({\mathrm{t}\overline{\mathrm{t}}} ) < $ 400 GeV and 0.0 $ < \beta < $ 0.9 for a $+$50% mixture of SC and noSC (bottom left) and a $-$50% mixture (bottom right). |
png pdf |
Figure 5-a:
Reconstruction-level distribution (top) of the combined $ \mathrm{t} \overline{\mathrm{t}} $ + $ \eta_{\mathrm{t}} $ signal model in mixtures of the noSC combined signal sample, more details in the text. Template variations are shown as a function of $ \cos \varphi $ in the phase space of 345 $ < m({\mathrm{t}\overline{\mathrm{t}}} ) < $ 400 GeV and 0.0 $ < \beta < $ 0.9 for a $+$50% mixture of SC and noSC (bottom left) and a $-$50% mixture (bottom right). |
png pdf |
Figure 5-b:
Reconstruction-level distribution (top) of the combined $ \mathrm{t} \overline{\mathrm{t}} $ + $ \eta_{\mathrm{t}} $ signal model in mixtures of the noSC combined signal sample, more details in the text. Template variations are shown as a function of $ \cos \varphi $ in the phase space of 345 $ < m({\mathrm{t}\overline{\mathrm{t}}} ) < $ 400 GeV and 0.0 $ < \beta < $ 0.9 for a $+$50% mixture of SC and noSC (bottom left) and a $-$50% mixture (bottom right). |
png pdf |
Figure 5-c:
Reconstruction-level distribution (top) of the combined $ \mathrm{t} \overline{\mathrm{t}} $ + $ \eta_{\mathrm{t}} $ signal model in mixtures of the noSC combined signal sample, more details in the text. Template variations are shown as a function of $ \cos \varphi $ in the phase space of 345 $ < m({\mathrm{t}\overline{\mathrm{t}}} ) < $ 400 GeV and 0.0 $ < \beta < $ 0.9 for a $+$50% mixture of SC and noSC (bottom left) and a $-$50% mixture (bottom right). |
png pdf |
Figure 6:
Postfit detector level distribution of $ \cos \varphi \otimes m({\mathrm{t}\overline{\mathrm{t}}} ) $ for 345 $ < m({\mathrm{t}\overline{\mathrm{t}}} ) < $ 400 GeV and 0.0 $ < \beta < $ 0.9. |
png pdf |
Figure 7:
Result of the scan of the quantity $ -2\ln L $ resulting from a maximum likelihood fit as a function of the parameter of interest, $ D $. The region where entanglement is not implied ($ D > -1/ $ 3) is indicated by the pink shading. |
png pdf |
Figure 8:
Summary of the measurement of the entanglement proxy $ D $ compared with the SM expectations including (filled) or not including (open) contributions from the hypothetical toponium state. |
Tables | |
png pdf |
Table 1:
The number of observed and predicted events categorized into signal and background contributions. All selection requirements are applied, i.e., the final event yield includes all requirements. Finally, we also apply an electroweak correction to the MC signal samples of $ \mathrm{t} \overline{\mathrm{t}} $ and $ \eta_{\mathrm{t}} $. The uncertainties reflect those originating from the limited MC sample size. The ``Only $ \eta_{\mathrm{t}} $" contribution is not added to the total MC prediction since it is included in the combined signal contribution. |
png pdf |
Table 2:
The number of predicted and observed events in the finally selected phase space, before the fit to the data (pre-fit) and with their best fit normalizations (post-fit). The uncertainties in the pre-fit yields originate from limited MC statistics, while for the uncertainties of post-fit yields total uncertainties including systematic ones are provided (but without correlations). Please note that the ``Only $ \eta_{\mathrm{t}} $" contribution is not added to the total MC prediction since it is included in the combined signal contribution. |
png pdf |
Table 3:
Breakdown into leading ten systematic uncertainties in the entanglement proxy $ D $ at the postfit level. |
Summary |
A measurement of the entanglement of top quark pairs $ \mathrm{t} \overline{\mathrm{t}} $ utilizing the spin correlation variable $ D $ is presented. This entanglement proxy is measured using events containing two oppositely charged leptons (including semileptonic decays of $ \tau $ leptons) produced in proton-proton collisions at a center-of-mass energy of 13 TeV. The extent to which top quarks are entangled is measured by means of a binned profile likelihood fit of the parameter of interest $ D $ from the distribution of $ \cos\varphi $ in the most sensitive kinematic phase space of 345 $ < m({\mathrm{t}\overline{\mathrm{t}}} ) < $ 400 GeV and 0.0 $ < \beta < $ 0.9. The value of the entanglement proxy $ D $ itself is measured by a negative log-likelihood scan of the parameter of interest $ D $ and yields an expected value of $ D = -$0.465 $ ^{+0.025}_{-0.028} $ and an observed value of $ D = -$0.478 $ ^{+0.025}_{-0.027} $. This result has an expected (observed) significance of 5.1 (5.7) standard deviations, corresponding to the observation of top quark entanglement. This measurement with CMS data of the quantum mechanics phenomena of entanglement present in top quark events provides a new and fundamental probe to beyond the standard model contributions and represents a first step towards quantum tomography and coherence in top quark events. |
References | ||||
1 | J. S. Bell | On the Einstein Podolsky Rosen paradox | Physics Physique Fizika 1 (1964) 195 | |
2 | F. M. Rafael Aoude, Eric Madge and L. Mantani | Quantum SMEFT tomography: Top quark pair production at the LHC | PRD 106 (2022) | |
3 | Y. Afik and J. R. M. n. de Nova | Entanglement and quantum tomography with top quarks at the LHC | Eur. Phys. J. Plus 136 (2021) 907 | 2003.02280 |
4 | B. Fuks, K. Hagiwara, K. Ma, and Y.-J. Zheng | Signatures of toponium formation in LHC run 2 data | PRD 104 (2021) 034023 | 2102.11281 |
5 | W.-L. Ju et al. | Top quark pair production near threshold: single/double distributions and mass determination | JHEP 06 (2020) 158 | 2004.03088 |
6 | CMS Collaboration | Measurement of the top quark polarization and $ \mathrm{t\bar{t}} $ spin correlations using dilepton final states in proton-proton collisions at $ \sqrt{s} = $ 13 TeV | PRD 100 (2019) 072002 | CMS-TOP-18-006 1907.03729 |
7 | W. Bernreuther, D. Heisler, and Z.-G. Si | A set of top quark spin correlation and polarization observables for the LHC: Standard Model predictions and new physics contributions | JHEP 12 (2015) 026 | 1508.05271 |
8 | A. Peres | Separability criterion for density matrices | PRL 77 (1996) 1413 | |
9 | M. Horodecki, P. Horodecki, and R. Horodecki | Separability of mixed states: necessary and sufficient conditions | Physics Letters A 223 (1996) 1 | |
10 | CMS Collaboration | Projection of the top quark spin correlation measurement and search for top squark pair production at the HL-LHC | CMS Physics Analysis Summary, 2022 CMS-PAS-FTR-18-034 |
CMS-PAS-FTR-18-034 |
11 | CMS Collaboration | Measurement of the differential cross section for top quark pair production in pp collisions at $ \sqrt{s} = 8\,\text {TeV} $ | EPJC 75 (2015) 542 | CMS-TOP-12-028 1505.04480 |
12 | D0 Collaboration | Measurement of the top quark mass using dilepton events. | PRL 80 (1998) 2063 | hep-ex/9706014 |
13 | CMS Collaboration | The CMS experiment at the CERN LHC | JINST 3 (2008) S08004 | |
14 | CMS Collaboration | Development of the CMS detector for the CERN LHC Run 3 | Accepted by JINST, 2023 | CMS-PRF-21-001 2309.05466 |
15 | CMS Collaboration | Particle-flow reconstruction and global event description with the CMS detector | JINST 12 (2017) P10003 | CMS-PRF-14-001 1706.04965 |
16 | CMS Collaboration | Pileup mitigation at CMS in 13 TeV data | JINST 15 (2020) P09018 | CMS-JME-18-001 2003.00503 |
17 | CMS Collaboration | Jet energy scale and resolution in the CMS experiment in pp collisions at 8 TeV | JINST 12 (2017) P02014 | CMS-JME-13-004 1607.03663 |
18 | CMS Collaboration | Performance of missing transverse momentum reconstruction in proton-proton collisions at $ \sqrt{s} = $ 13 TeV using the CMS detector | JINST 14 (2019) P07004 | CMS-JME-17-001 1903.06078 |
19 | CMS Collaboration | Electron and photon reconstruction and identification with the CMS experiment at the CERN LHC | JINST 16 (2021) P05014 | CMS-EGM-17-001 2012.06888 |
20 | CMS Collaboration | ECAL 2016 refined calibration and Run2 summary plots | CMS Detector Performance Summary CMS-DP-2020-021, 2020 CDS |
|
21 | CMS Collaboration | Performance of the CMS muon detector and muon reconstruction with proton-proton collisions at $ \sqrt{s}= $ 13 TeV | JINST 13 (2018) P06015 | CMS-MUO-16-001 1804.04528 |
22 | S. Frixione, P. Nason, and G. Ridolfi | A Positive-weight next-to-leading-order Monte Carlo for heavy flavour hadroproduction | JHEP 09 (2007) 126 | 0707.3088 |
23 | P. Nason | A New method for combining NLO QCD with shower Monte Carlo algorithms | JHEP 11 (2004) 040 | hep-ph/0409146 |
24 | S. Frixione, P. Nason, and C. Oleari | Matching NLO QCD computations with Parton Shower simulations: the POWHEG method | JHEP 11 (2007) 070 | 0709.2092 |
25 | S. Alioli, P. Nason, C. Oleari, and E. Re | A general framework for implementing NLO calculations in shower Monte Carlo programs: the POWHEG BOX | JHEP 06 (2010) 043 | 1002.2581 |
26 | M. Czakon and A. Mitov | Top++: A Program for the Calculation of the Top-Pair Cross-Section at Hadron Colliders | Comput. Phys. Commun. 185 (2014) 2930 | 1112.5675 |
27 | J. Alwall et al. | The automated computation of tree-level and next-to-leading order differential cross sections, and their matching to parton shower simulations | JHEP 07 (2014) 079 | 1405.0301 |
28 | P. Artoisenet, R. Frederix, O. Mattelaer, and R. Rietkerk | Automatic spin-entangled decays of heavy resonances in Monte Carlo simulations | JHEP 03 (2013) 015 | 1212.3460 |
29 | R. Frederix and S. Frixione | Merging meets matching in MC@NLO | JHEP 12 (2012) 061 | 1209.6215 |
30 | NNPDF Collaboration | Parton distributions for the LHC Run II | JHEP 04 (2015) 040 | 1410.8849 |
31 | J. Alwall et al. | Comparative study of various algorithms for the merging of parton showers and matrix elements in hadronic collisions | EPJC 53 (2008) 473 | 0706.2569 |
32 | T. Sj\"ostrand et al. | An introduction to PYTHIA 8.2 | Comput. Phys. Commun. 191 (2015) 159 | 1410.3012 |
33 | CMS Collaboration | Investigations of the impact of the parton shower tuning in Pythia 8 in the modelling of $ \mathrm{t\overline{t}} $ at $ \sqrt{s}= $ 8 and 13 TeV | CMS Physics Analysis Summary, 2016 CMS-PAS-TOP-16-021 |
CMS-PAS-TOP-16-021 |
34 | CMS Collaboration | Event generator tunes obtained from underlying event and multiparton scattering measurements | EPJC 76 (2016) 155 | CMS-GEN-14-001 1512.00815 |
35 | M. Bahr et al. | Herwig++ Physics and Manual | EPJC 58 (2008) 639 | 0803.0883 |
36 | S. Gieseke, C. Rohr, and A. Siodmok | Colour reconnections in Herwig++ | EPJC 72 (2012) 2225 | 1206.0041 |
37 | GEANT4 Collaboration | Geant4-a simulation toolkit | Nuclear Instruments and Methods in Physics Research Section A 506 (2003) 250 | |
38 | CMS Collaboration | Measurements of $ \mathrm{t\overline{t}} $ differential cross sections in proton-proton collisions at $ \sqrt{s}= $ 13 TeV using events containing two leptons | JHEP 02 (2019) 149 | CMS-TOP-17-014 1811.06625 |
39 | J. A. Aguilar-Saavedra and J. A. Casas | Improved tests of entanglement and Bell inequalities with LHC tops | EPJC 82 (2022) 666 | 2205.00542 |
40 | CMS Collaboration | Measurement of the $ t\bar t $ production cross section in the dilepton channel in pp collisions at $ \sqrt{s} = $ 7 TeV | JHEP 11 (2012) 067 | CMS-TOP-11-005 1208.2671 |
41 | CMS Collaboration | Precision luminosity measurement in proton-proton collisions at $ \sqrt{s} = $ 13 TeV in 2015 and 2016 at CMS | EPJC 81 (2021) 800 | CMS-LUM-17-003 2104.01927 |
42 | CMS Collaboration | Combination of measurements of the top quark mass from data collected by the ATLAS and CMS experiments at $ \sqrt{s}= $ 7 and 8 TeV | CMS-TOP-22-001 2402.08713 |
|
43 | CMS Collaboration | Measurement of the shape of the b quark fragmentation function using charmed mesons produced inside b jets from $ \mathrm{t}\bar{\mathrm{t}} $ pair decays | CMS Physics Analysis Summary, 2021 CMS-PAS-TOP-18-012 |
CMS-PAS-TOP-18-012 |
44 | Particle Data Group Collaboration | Review of Particle Physics (2017 update) | Chin. Phys. C 40 (2016) 100001 | |
45 | CMS Collaboration | Measurement of the top quark Yukawa coupling from $ \mathrm{t\bar{t}} $ kinematic distributions in the dilepton final state in proton-proton collisions at $ \sqrt{s}= $ 13 TeV | PRD 102 (2020) 092013 | CMS-TOP-19-008 2009.07123 |
Compact Muon Solenoid LHC, CERN |