Measurement of the cross section for $e^{+}e^{-}\rightarrow\Lambda\bar\Lambda$ and evidence of the decay $\psi(3770)\rightarrow\Lambda\bar\Lambda$

The Born cross section of the process $e^{+}e^{-}\rightarrow\Lambda\bar\Lambda$ is measured at 33 center-of-mass energies between 3.51 and 4.60 GeV using data corresponding to the total integrated luminosity of 20.0 fb$^{-1}$ collected with the BESIII detector at the BEPCII collider. Describing the energy dependence of the cross section requires a contribution from the $\psi(3770)\rightarrow\Lambda\bar\Lambda$ decay, which is fitted with a significance of 4.6-4.9$\sigma$ including the systematic uncertainty. The lower bound on its branching fraction is $2.4\times10^{-6}$ at the 90% confidence level (C.L.), at least an order of magnitude larger than expected from predictions using a scaling based on observed electronic widths. This result indicates the importance of effects from vector charmonium(-like) states when interpreting data in terms of e.g., electromagnetic structure observables. The data do not allow for definite conclusions on the interplay with other vector charmonium(-like) states, and we set 90% C.L.upper limits for the products of their electronic widths and the branching fractions.

(Dated: December 1, 2021) The Born cross section of the process e + e − → ΛΛ is measured at 33 center-of-mass energies between 3.51 and 4.60 GeV using data corresponding to the total integrated luminosity of 20.0 fb −1 collected with the BESIII detector at the BEPCII collider. Describing the energy dependence of the cross section requires a contribution from the ψ(3770) → ΛΛ decay, which is fitted with a significance of 4.6−4.9σ including the systematic uncertainty. The lower bound on its branching fraction is 2.4 × 10 −6 at the 90% confidence level (C.L.), at least an order of magnitude larger than expected from predictions using a scaling based on observed electronic widths. This result indicates the importance of effects from vector charmonium(-like) states when interpreting data in terms of e.g., electromagnetic structure observables. The data do not allow for definite conclusions on the interplay with other vector charmonium(-like) states, and we set 90% C.L. upper limits for the products of their electronic widths and the branching fractions.
Two-body baryonic decays of vector (J P C = 1 −− ) charmonium(-like) resonances provide a testing ground of predictions from quantum chromodynamics [1,2]. The ψ(3770) vector meson is believed to be a conventional cc state located above the open-charm threshold and is expected to decay into a DD meson pair with a branching fraction of at least 99% [3]. However, the decay modes to light quark systems would be considerably enhanced if the ψ(3770) would include gluonic or light quark and antiquark constituents [4]. In 2003, the BES collaboration observed the first non-DD decay of ψ(3770) into J/ψπ + π − [5,6]. Subsequently, the CLEO collaboration confirmed the observation and found more non-DD decays of ψ(3770) [7][8][9] and the first decay into lightquark hadrons ψ(3770) → φη [10].
The production of light quark baryon-antibaryon (BB) final states leads to relatively simple topologies. In an early study, the BESIII collaboration found evidence for the interference effect in e + e − → pp in the vicinity of ψ(3770) [11]. However, the data did not allow to uncover the mechanism of ψ(3770) charmless decays. Thus, the experimental study of e + e − → BB will be a good search-ground for clarifying the nature of the charmless decays and even non-DD decays of ψ(3770) [12][13][14].
In this paper, we present a measurement of the Born cross section for the e + e − → ΛΛ process using data corresponding to a total integrated luminosity of 20.0 fb −1 [36][37][38] collected at center-of-mass (c.m.) energies √ s between 3.51 and 4.60 GeV with the BESIII detector [39,40] at the BEPCII collider [41]. We extract the Λ effective form factor and report an evidence of the ψ(3770) → ΛΛ process by fitting the e + e − → ΛΛ dressed cross section.
Candidates for e + e − → ΛΛ events are reconstructed using the Λ → pπ − andΛ →pπ + decay modes. The detection efficiency is determined by Monte Carlo (MC) simulations. A sample of 100,000 events is simulated for each of the 33 c.m. energy points. The production process is simulated by the KKMC generator [42,43] that includes corrections for ISR effects. The Λ andΛ decays are handled by the EVT-GEN [44,45] program. The response of the BESIII detector is modeled with MC simulations using a framework based on GEANT4 [46,47].
Tracks of charged particles are reconstructed in the multilayer drift chamber with a helical fit requiring a good quality [48]. These tracks should be within | cos θ| < 0.93, where θ is the polar angle with respect to the e + beam direction. Events with two successfully reconstructed negatively charged and two positively charged particles are kept for further analysis.
To reconstruct Λ(Λ) candidates, we apply a secondary vertex fit [49] to all pairs of positive and negative charged particles.
The corresponding χ 2 value is required to be less than 500.
The track combination with mini- is the invariant mass of the pπ − (pπ + ) pair, and m Λ is the world-average Λ mass value from the Particle Data Group (PDG) [15]. To further suppress background from non-Λ processes, the Λ decay length is required to be larger than zero, where the observed negative decay lengths are caused by the limited detector resolution. Here the misreconstruction ratio for Λ particle is found to be less than 1% based on the study of MC simulation.
To further suppress the background and to improve the mass resolution, a four-constraint (4C) kinematic fit imposing energy-momentum conservation is applied for the ΛΛ hypothesis. The χ 2 4C of the fit is required to be less than 200, which can improve the resolutions of signals significantly in addition to suppressing backgrounds for a soft π 0 and radiated photons event. Figure 1 shows the distribution of M pπ − versus Mp π + of the accepted candidates from all data samples. Clear peaks around the Λ known mass can be discerned. The invariant mass M pπ − is required to be within 5 MeV/c 2 of the known Λ(Λ) mass (signal region marked by S). After applying the above selection criteria, the survived background events are mainly from non-Λ(Λ) events, such as e + e − → π + π − pp. The background yield in the signal region is estimated using four sideband regions B i , where i = 1, 2, 3, 4, each with the same area as the signal region. The regions are shown in Fig. 1, and the exact ranges are given in the Supplemental Material [36]. The signal yield N obs for e + e − → ΛΛ events at each energy point can then be extracted by subtracting the number of events in the sideband regions from the number of events in the signal region, N S : and they are listed in Table I. The ISR corrected ("dressed") cross section σ dr for the process e + e − → ΛΛ is defined as where L is the integrated luminosity at given c.m. energy √ s, (1 + δ) is the ISR correction factor [43,50], ǫ is the detection efficiency, and the branching fraction B(Λ → pπ − ) = (63.9 ± 0.5)% is taken from PDG. The ISR correction factor is obtained using the calculation described in Ref. [51], where the dressed cross sections are adopted as initial input and are iterated to obtain stable result. The dressed cross section is related to the Born cross section via the vacuum polarization factor 1 |1−Π| 2 [52] as σ dr = σ B /|1 − Π| 2 (further details are provided in the Supplemental Material [36]).
Systematic uncertainties on the cross section measurement mainly come from the luminosity measurement, the Λ reconstruction, the 4C kinematic fit, the branching fraction for the decay Λ → pπ − , the line-shape description, and the physical model dependence. The uncertainty due to the vacuum polarization is negligible. The integrated luminosity is measured by e + e − → (γ)e + e − events with a similar method to Ref. [37] with an uncertainty of 1.0%. The systematic uncertainty of the Λ(Λ) reconstruction incorporating the tracking, the mass window of Λ(Λ), and the decay length of Λ(Λ) is studied using a control sample of ψ(3686) → ΛΛ decay (∼ 20000 events) with the same method as introduced in Refs. [53][54][55][56][57][58][59]. The signal MC sample is simulated using a DIY model [44] implementing the joint angular distribution from Refs. [60,61]. The efficiency difference between data and MC simulation is found to be 0.5% for the Λ reconstruction and 1.5% for theΛ reconstruction. The uncertainty from the 4C kinematic fit is studied using the control sample of ψ(3686) → ΛΛ decays with and without performing a 4C kinematic fit. The relative change of 1.0% is assigned as the systematic uncertainty. The uncertainty of the branching fraction for Λ → pπ − from the PDG [15] is 0.8%, and is propagated to the final result. The uncertainty from the line-shape description is estimated with an alternative input cross section line shape based on a simple power-law function. The change of the efficiency, 2.6%, is taken as the systematic uncertainty. The uncertainty due to the physical model dependence is estimated to be 2.5% by comparing the efficiencies between phase space and the DIY model incorporating the Λ transverse polarization and spin correlation based on the control sample of ψ(3770) → ΛΛ decays. Assuming all sources are independent, the total systematic uncertainty on the cross section measurement is determined to be 4.3% by adding these sources in quadrature. The correlations for the different points are negligible due to the limited statistics. The extracted Born cross sections at each energy point are listed in Table I and shown in Fig. 2(top) together with the CLEO-c results at 3.770 and 4.160 GeV [34,35]. Our results are significantly lower than those of CLEO-c at both energy points. Figure 2 (bottom) shows the extracted energy dependence of the Λ effective form factor G eff (s) defined as [63] G eff (s) = 3sτ σ B 2πα 2 β(2τ + 1) , where α is the fine-structure constant, β = (τ − 1)/τ is the Λ velocity and τ = s/(4m 2 Λ ). The dressed cross section for the continuum e + e − → ΛΛ process is expected to have an asymptotic power-law behavior ∝ s −n with the exponent n ≈ 10 [63][64][65]. A least-χ 2 fit including statistical and systematic uncertainties to the powerlaw distribution describes the data points reasonably well, as shown with the dashed line in Fig. 3. The fitted value of the exponent n is not close to 10 within the uncertainty of 1σ, as shown in the column "Fit I" in Table II. A fit with the coherent sum of the power-law function and a Breit-Wigner (BW) function is applied, where M is the ψ(3770) mass, σ 0 is the value of the continuum cross section at ψ(3770) and φ is the relative phase between the continuum and the resonance. The BW function is where Γ and Γ ee are the total and the electronic width of the ψ(3770) resonance, respectively, and B denotes the branching fraction to ΛΛ. The solid line in Fig. 3 and the column "Fit II" in Table II shows the result of the fit with two solutions, where the mass and width of ψ(3770) are fixed according to the PDG values [15], and σ 0 , n, φ and σ ψ parameters are free. The improvement of the χ 2 value gives a significance of 4.6−4.9σ for the hypothesis with the ψ(3770) resonance including the systematic uncertainty. The correlation coefficient between the resonance cross section σ ψ and the phase φ is almost equal to one. In the Fit II, two solutions are expected according to mathematical calculation [66,67], but fits give the consistent results within the uncertainty of 1σ due to statistics limitation. Figure 4 shows the contour of B and φ on the distribution of χ 2 values for each set of parameters. Our results can be summarized by giving 90% C.L. intervals 24 < σ ψ < 1800 fb and 2.4 × 10 −6 < B < 1.8 × 10 −4 . This represents the first evidence of the decay ψ(3770) → ΛΛ. This result is larger by at least an order of magnitude than the prediction based on a scaling from the electronic branching fraction value. This implies, that the ψ(3770) resonance needs to be considered when interpreting the CLEO-c data. Note that the systematic uncertainties due to beam energy, mass and width of the ψ(3770) resonance have been considered by varying the known value within one standard deviation, and they turn out to be negligible.
Finally, we have included an additional charmonium(-like) state (i.e. ψ(4040), ψ(4160), Y (4260), or ψ(4415)) in the fit, one at a time. It turns out that the exponent and significance for the ψ(3770) state are consistent with a single reso-  TABLE II. Results of the fit to the dressed cross section for the e + e − → ΛΛ process, where two solutions in "Fit II" are provided. The fitting procedure includes both statistical and systematic uncertainties except for the c.m. energy calibration. B is the branching fraction of the decay ψ(3770) → ΛΛ measured assuming Bee = 9.7 × 10 −6 , the central value of the world average [15], and Γee = Γ ψ(3770) Bee = (261.2 ± 21.3) eV. nance assumption. Since the significance of each mentioned state is smaller than 3σ, we quote upper limits at the 90% C.L. for the Γ ee B products: < 5.5 × 10 −3 eV for ψ(4040), < 0.7 × 10 −3 eV for ψ(4160), < 0.8 × 10 −3 eV for Y (4260) and < 1.8 × 10 −3 eV for ψ(4415) including the systematic uncertainty. These results provide important information to understand the nature of charmonium(-like) states above the open charm threshold. In particular this concerns their coupling to the BB final states and insight into the puzzle of large non-DD component of ψ(3770).