Observation of Transverse Polarization and Determination of Electromagnetic Form Factor of $\Lambda$ Hyperon at $\sqrt{s}=3.773$ GeV

Abstract

Using a 20.3 fb^-1 of $e^{+}e^{-}$ collision data sample collected by the BESIII detector at the BEPCII collider, we present an observation of transverse polarization and a complete determination of the electromagnetic form factor of the $\Lambda$ hyperon in $e^{+}e^{-}\to{\Lambda}\bar{\Lambda}$ decay with the entangled ${\Lambda}-\bar{\Lambda}$ pair at $\sqrt{s}=3.773$ GeV. The relative phase between the electric and magnetic form factors is determined to be $\Delta\Phi=(1.53\pm 0.36\pm 0.03)$ rad with a significance of 5.5 $\sigma$ taking into account systematic uncertainty. This result indicates a non-zero phase between the transition amplitudes of the ${\Lambda}\bar{\Lambda}$ helicity states. Additionally, we measure the angular distribution parameter and the modulus of the ratio between the electric and the magnetic form factor is found to be $\eta=0.86\pm 0.05\pm 0.03$ and $R(s)=|G_{E}(s)/G_{M}(s)|=0.47\pm 0.08\pm 0.05$, where the first uncertainty is statistical and the second systematic.

Understanding the structure of baryons is a major goal of contemporary particle physics Brodsky:1974vy ; Geng:2008mf ; Green:2014xba ; Wang:2022zyc . In the context of Quantum Chromodynamics (QCD), electromagnetic form factors (EMFFs) serve as important observables that connect measurable quantities to theoretical predictions. In the 1960s, Cabibbo et al. Cabibbo:1961sz first proposed that timelike EMFFs could be studied at $e^{+}e^{-}$ experiments by measuring the production cross sections of baryon-antibaryon pairs. A large amount of research has been carried out regarding nucleon and strange hyperon EMFFs in the timelike momentum transfer regions ($s>0$) ParticleDataGroup:2024cfk , where $s$ is the square of the center-of-mass (c.m.) energy. Among them, the proton, being a stable particle, can serve as a suitable target for studying its spacelike EMFFs through scattering experiments. A recent study also revealed that for large $s$, the $d$-quark contributions to the proton EMFFs are reduced relative to the $u$-quark contributions Cates:2011pz . This sets the proton apart from unstable hyperons with finite lifetimes, which are not suitable for such scattering experiments. Instead, the interaction in $e^{+}e^{-}$ collisions allows access to timelike EMFFs of hyperons due to virtual photon production in the $e^{+}e^{-}$ interaction which facilitates the quantitative assessment of the electromagnetic structure. Experimentally accessible timelike EMFFs are connected with the spacelike quantities, such as charge and magnetization densities, through the dispersion relation Belushkin:2006qa . The pair production of spin-$1/2$ baryons can then be parametrized by the electric form factor $G_{E}(s)$ and the magnetic form factor $G_{M}(s)$ Wang:2022zyc ; Qian:2021neg , which are analytic functions of the momentum transfer squared. In the timelike region, the EMFFs are complex and have a relative phase $\Delta\Phi=\Delta\Phi_{E}-\Delta\Phi_{M}$, i.e. $G_{E}(s)/G_{M}(s)=Re^{i\Delta\Phi}$ with the modulus of EMFFs ratio $R=|G_{E}(s)/G_{M}(s)|$. This relative phase $\Delta\Phi$ reflects interfering production amplitudes and has a polarizing effect on the final state even if the initial state is unpolarized Dubnickova:1992ii . This provides a handle to study the asymptotic properties of the timelike EMFFs related to the intrinsic structure of hyperons at large $s$, where the spacelike and timelike EMFFs should converge to the same value. For protons, the onset of this scale can be studied by measuring spacelike and timelike EMFFs. For ground-state hyperons, on the other hand, the weak parity-violating decays provide straightforward access to their polarization.

Experimentally, the first determination of the effective form factor of the $\Lambda$ hyperon was reported by the BABAR experiment using the initial state radiation (ISR) method Aubert:2007uf . Subsequently, the CLEO-c experiment reported the measurements of the timelike EMFFs of several baryons Dobbs:2014ifa ; Dobbs:2017 . Results yielded timelike EMFFs and emphasized the importance of diquark correlations Jaffe:2003sg under the assumption that one-photon exchange dominates the production process and the charmonia decay contributions are negligible. Here, the spin formalism introduced in Ref. Faldt:2017kgy is also appropriate for the vicinity of vector charmonia. In this case, the form of the hadron current matrix element for the charmonia process is the same as for the virtual photon process. Here, the timelike EMFFs referring to the amplitudes for $e^{+}e^{-}\to\Lambda\bar{\Lambda}$ process represent hadronic form factors. Recently, the BESIII collaboration performed a pioneering measurement of the relative phase, the modulus of the timelike EMFFs ratio and spin polarizations of the hyperons near threshold BESIII:2019nep ; BESIII:2020uqk ; BESIII:2022kzc ; BESIII:2023ynq , around the resonances of vector charmonia BESIII:2018cnd ; BESIII:2020fqg ; BESIII:2021ypr ; BESIII:2022qax ; BESIII:2022lsz ; BESIII:2023lkg ; BESIII:2023jhj ; BESIII:2024nif and above the open charm threshold BESIII:2019cuv ; BESIII:2021ccp ; BESIII:2021cvv ; BESIII:2023rse ; BESIII:2021cvv ; BESIII:2024umc by considering quantum entanglement of hyperon and antihyperon. This resulted in increased activity within the theoretical community, encompassing a variety of approaches including hyperon-antihyperon final state interactions Haidenbauer:2020wyp ; Qian:2022whn , vector meson dominance Yang:2019mzq ; Li:2021lvs ; Yan:2023yff ; Chen:2023oqs , the covariant spectator model Ramalho:2019koj ; Ramalho:2024wxp , and dispersive calculations Lin:2022baj ; Mangoni:2021qmd . In particular, Ref. Mangoni:2021qmd proposed a method for examining zero crossings and extracting the charge radius of the $\Lambda$ hyperon by combining comprehensive single-energy timelike EMFFs measured by the BESIII experiment BESIII:2019nep with the partial timelike EMFF measured at two energy points by the BABAR experiment Aubert:2007uf . However, due to the limited data and large uncertainties in the BABAR measurement, a definitive solution was not possible. With a larger data sample, it is anticipated that the energy dependence of the timelike EMFFs phase will provide constraints for calculations and potentially allow for the extraction of the charge radius of hyperons. In this Letter, we report an observation of spin polarization and determination of timelike EMFFs of the $\Lambda$ hyperon using a multi-dimensional angular distribution analysis with a complete decomposition of the spin structure of the process of $e^{+}e^{-}\to\Lambda\bar{\Lambda}$. The data set corresponding to an integrated luminosity of 20.3 fb^-1 collected at $\sqrt{s}=3.773$ GeV Ablikim:2014gna ; BESIII:2024lbn by the BESIII detector Wang:2007tv at the BEPCII collider BESIII is analyzed, which is about seven times larger than the one used in the previous study BESIII:2021cvv .

To describe the process of $e^{+}e^{-}\to{\Lambda}\bar{\Lambda}\to{p\bar{p}\pi^{+}\pi^{-}}$, it is essential to acquire information about each particle in the coordinate system of the parent particle. Here, a right-handed coordinate system is used according to that was provided in Ref. zhangzhe to describe hyperon decay and orientation of $p/\bar{p}$, see Fig. 1.

Figure 1: The coordinate system utilized to describe the $e^{+}e^{-}\to\Lambda\bar{\Lambda}\to p\bar{p}\pi^{+}\pi^{-}$ process. The $\hat{z}_{e}$ axis is defined as the direction of $e^{+}$. The $\hat{z}_{\Lambda/\bar{\Lambda}}$ axis is determined as the direction of the $\Lambda/\bar{\Lambda}$ particle emission. The $\hat{y}_{\Lambda/\bar{\Lambda}}=\hat{z}_{e}\times\hat{z}_{\Lambda/\bar{\Lambda}}$, and $\hat{x}_{\Lambda/\bar{\Lambda}}=\hat{y}_{\Lambda/\bar{\Lambda}}\times\hat{z}_{\Lambda/\bar{\Lambda}}$. The $\theta_{\Lambda}$ is the angle between the $\Lambda$ hyperon and $e^{+}$ in $e^{+}e^{-}$ c.m. The angles $\theta_{p/\bar{p}}$ and $\phi_{p/\bar{p}}$ are the polar and azimuthal angles of the $p/\bar{p}$ momentum direction in the $\Lambda/\bar{\Lambda}$ rest frame, respectively. The $v_{\Lambda/\bar{\Lambda}}^{LT}$ is the Lorentz transformation with the velocity of $v_{\Lambda/\bar{\Lambda}}$.

In Ref. zhangzhe , the timelike EMFFs ratio $(R)$, the relative phase $\Delta\Phi$ and the angular distribution parameter $\eta$ are used to describe the process $e^{+}e^{-}\to{\gamma}^{*}/\psi\to{\Lambda}\bar{\Lambda}$. For $e^{+}e^{-}\to{\gamma^{*}/\psi}\to{\Lambda}\bar{\Lambda}\to{p\bar{p}\pi^{+}\pi^{-}}$, the joint decay angular distribution of this process is expressed in terms of the parameters of $\Delta\Phi$ and $\eta$ as

$$\mathcal{W}(\bm{\xi};\bm{\Omega})=\mathcal{F}_{0}(\xi)+\eta\mathcal{F}_{5}(\xi)$$ (1)

$$\hskip 10.00002pt+\alpha_{\Lambda}\alpha_{\bar{\Lambda}}[\mathcal{F}_{1}(\xi)+\sqrt{1-\eta^{2}}\cos(\Delta\Phi)\mathcal{F}_{2}(\xi)+\eta\mathcal{F}_{6}(\xi)]$$

$$\hskip 10.00002pt+\sqrt{1-\eta^{2}}\sin(\Delta\Phi)[\alpha_{\Lambda}\mathcal{F}_{3}(\xi)+\alpha_{\bar{\Lambda}}\mathcal{F}_{4}(\xi)],$$

where the angular functions $\mathcal{F}_{j}(\xi)$ $(j=0,1,...6)$ are defined as:

$${\cal{F}}_{0}= 1,$$

$${\cal{F}}_{1}= \sin^{2}\!\theta_{\Lambda}\sin\!\theta_{p}\sin\!\theta_{\bar{p}}\cos\!\phi_{p}\cos\!\phi_{\bar{p}}-\cos^{2}\!\theta_{\Lambda}\cos\!\theta_{p}\cos\!\theta_{\bar{p}},$$

$${\cal{F}}_{2}= \sin\!\theta_{\Lambda}\cos\!\theta_{\Lambda}(\sin\!\theta_{p}\cos\!\theta_{\bar{p}}\cos\!\phi_{p}-\cos\!\theta_{p}\sin\!\theta_{\bar{p}}\cos\!\phi_{\bar{p}}),$$

$${\cal{F}}_{3}= -\sin\!\theta_{\Lambda}\cos\!\theta_{\Lambda}\sin\!\theta_{p}\sin\!\phi_{p},$$ (2)

$${\cal{F}}_{4}= \sin\!\theta_{\Lambda}\cos\!\theta_{\Lambda}\sin\!\theta_{\bar{p}}\sin\!\phi_{\bar{p}},$$

$${\cal{F}}_{5}= \cos^{2}\!\theta_{\Lambda},$$

$${\cal{F}}_{6}= \sin^{2}\!\theta_{\Lambda}\sin\!\theta_{p}\sin\!\theta_{\bar{p}}\sin\!\phi_{p}\sin\!\phi_{\bar{p}}-\cos\!\theta_{p}\cos\!\theta_{\bar{p}}.$$

Furthermore, $\alpha_{\Lambda(\bar{\Lambda})}$ represents the decay parameters of $\Lambda(\bar{\Lambda})\to p\pi^{-}(\bar{p}\pi^{+})$ and $\eta$ denotes the scattering angle distribution parameter related to $R$ by

$$\eta=\frac{\tau-R^{2}}{\tau+R^{2}},$$ (3)

with $\tau=s/4m^{2}_{\Lambda}$. Since the production process is either strong or electromagnetic and thus parity conserving, if the initial state is unpolarized, non-zero transverse polarization can only occur in the transverse, or $y$ direction. The spin polarization is defined as

$$P_{y}=\frac{\sqrt{1-\eta^{2}}\sin\theta_{\Lambda}\cos\theta_{\Lambda}}{1+\eta\cos^{2}\theta_{\Lambda}}\sin(\Delta\Phi).$$ (4)

For further analysis, fully reconstructed $e^{+}e^{-}\to{\Lambda}\bar{\Lambda}$ events with $\Lambda\to{p\pi^{-}}$ and $\bar{\Lambda}\to{\bar{p}\pi^{+}}$ are selected. To determine the detection efficiency and perform the unbinned maximum likelihood fit, $10^{7}$ Monte Carlo (MC) events are generated using kkmc kkmc ; kkmc-01 , which includes the ISR effect. The $e^{+}e^{-}\to\Lambda\bar{\Lambda}$ and ${\Lambda}(\bar{\Lambda})$ to $p\pi^{-}(\bar{p}\pi^{+})$ decays are simulated according to a phase space (PHSP) model using evtgen evt2 ; evt2-01 . The response of the BESIII detector is modeled with simulation using a framework based on geant4 geant4 ; geant4-01 .

Charged tracks are reconstructed in the multi-layer drift chamber within its angular coverage, $|\cos\theta|<0.93$, where $\theta$ is the polar angle with respect to the $e^{+}$ beam direction in the laboratory system. The numbers of negatively and positively charged tracks of events are both larger than one. Charged tracks with momenta greater than 0.6 GeV/$c$ are identified as $p(\bar{p})$, while others are assigned as $\pi^{+}(\pi^{-})$.

To reconstruct $\Lambda(\bar{\Lambda})$ candidates, a vertex fit and a secondary vertex fit XUM are applied to all combinations of one $p(\bar{p})$ track and one $\pi^{-}(\pi^{+})$ track. From all possible combinations, the one with the minimum value of $\sqrt{|M_{p\pi^{-}}-m_{\Lambda}|^{2}+|M_{\bar{p}\pi^{+}}-m_{\bar{\Lambda}}|^{2}}$ is selected. Here, $M_{p\pi^{-}(\bar{p}\pi^{+})}$ denotes the invariant mass of the $p\pi^{-}(\bar{p}\pi^{+})$ pair, and $m_{\Lambda(\bar{\Lambda})}$ represents the nominal mass of $\Lambda(\bar{\Lambda})$ ParticleDataGroup:2024cfk . To further suppress background contributions from non-$\Lambda(\bar{\Lambda})$ events, the decay lengths of $\Lambda$ and $\bar{\Lambda}$ are both required to be greater than zero, where the negative decay lengths are due to detector resolution.

After $\Lambda(\bar{\Lambda})$ reconstruction, a four-constraint (4C) kinematic fit is applied to all $\Lambda(\bar{\Lambda})$ hypotheses, enforcing energy-momentum conservation from the initial $e^{+}e^{-}$ to the final $\Lambda(\bar{\Lambda})$ state and combined with the requirement of $\chi^{2}_{\rm{4C}}<100$. Figure 2 shows the distribution of $M_{\bar{p}\pi^{+}}$ versus $M_{p\pi^{-}}$ after applying all above selection.

Figure 2: Two-dimensional distribution of $M_{\bar{p}\pi^{+}}$ versus $M_{p\pi^{-}}$ for data, where the red box indicates the signal region, and the green boxes show the selected sideband regions.

The invariant mass of $p\pi^{-}(\bar{p}\pi^{+})$ is required to be within 5 MeV/$c^{2}$ of the $\Lambda(\bar{\Lambda})$ mass taken from ParticleDataGroup:2024cfk ($|M_{p\pi^{-}(\bar{p}\pi^{+})}-m_{\Lambda(\bar{\Lambda})}|<$ 5 MeV/$c^{2}$). The signal region, denoted by $S$ in Fig. 2, is determined and optimized using the figure of merit ${S}/{\sqrt{S+B}}$ derived from MC sample. Here, $S$ represents the number of signal MC events, and $B$ corresponds to the expected number of background events from the inclusive MC simulation of $e^{+}e^{-}\to$ hadron events. The background comes from non-$\Lambda(\bar{\Lambda})$ events, such as $e^{+}e^{-}\to\pi^{+}\pi^{-}p\bar{p}$, that can be estimated using the corner method, i.e., $\sum^{4}_{i=1}B_{i}/4$ for $M_{p\pi^{-}}$ and $M_{\bar{p}\pi^{+}}$ windows. Here, the definition for regions ($B_{i}$) displayed in Fig. 2 is the same as the one used in Ref. BESIII:2021cvv . The number of background events estimated from the aforementioned corner method is $4\pm 2$, which is negligible. To estimate peaking background contributions such as $e^{+}e^{-}\to\gamma\psi(3686)\to\gamma\Lambda\bar{\Lambda}$, an inclusive MC sample of $\psi(3770)$ is employed. The number of background events of ISR $\psi(3686)$ is estimated to be $39\pm 6$ events with a background level of approximately $1.8\%$ of the signal yield. The number of observed events in data is determined to be $2194\pm 48$ with a signal MC efficiency of (37.00$\pm$0.03)%.

To determine the set of $\Lambda$ spin polarization parameters $\bm{\Omega}=\{\Delta\Phi,\eta\}$, an unbinned maximum likelihood fit is performed. In the fit, $\alpha_{\Lambda/\bar{\Lambda}}$ is fixed to $\pm 0.7542$ by referring to BESIII:2022qax assuming charge-parity conservation. The likelihood function $\mathcal{L}$ is constructed from the probability density function (PDF), ${\cal{P}}({\bm{\xi}}_{i})$, for the event $i$ characterized by the measured angles $\bm{\xi}_{i}$ as

$$\mathcal{L}=\prod_{i=1}^{N}\mathcal{P}(\bm{\xi}_{i},\mathbf{\Omega})=\prod_{i=1}^{N}\mathcal{C}\mathcal{W}(\bm{\xi}_{i},\mathbf{\Omega})\epsilon(\bm{\xi}_{i}),$$ (5)

where $N$ is the number of events in the signal region. The joint angular distribution ${\cal{W}}({\bm{\xi}}_{i},{\bm{\Omega}})$ is given in Eq. (LABEL:w), and $\epsilon(\bm{\xi}_{i})$ is the detection efficiency. The normalization factor $\mathcal{C}^{-1}=\frac{1}{N_{\mathrm{MC}}}\sum_{j=1}^{N_{\mathrm{MC}}}{\cal{W}}({\bm{\xi}}^{j},{\bm{\Omega}})$ is calculated as a sum of the corresponding amplitudes $\cal{W}$ from the accepted PHSP MC events $N_{\mathrm{MC}}$, applying the same event selection criteria as to the data. The minimization of the objective function defined as

$$\mathit{S}=-\mathrm{ln}\mathcal{L}_{\rm S}+\mathrm{ln}\mathcal{L}_{\rm B},$$ (6)

is conducted using the MINUIT package from the ROOT library James:1975dr . In Eq. (6), $\mathcal{L}_{\rm S}$ and $\mathcal{L}_{\rm B}$ represent the likelihood function for events chosen in the signal region and sideband regions. Figure 3 shows the distributions of the five moments ${{F}}_{k}$ for $k=\{1,2,3,4,6\}$ defined in Eq. (Observation of Transverse Polarization and Determination of Electromagnetic Form Factor of $\Lambda$ Hyperon at $\sqrt{s}=3.773$ GeV) and $\Lambda$ angular distribution (${{F}}_{0}+\eta{{F}}_{5}$) with respect to $\cos\theta_{\Lambda}$ in 10 intervals.

Figure 3: The $F_{k}(k=1,2,...6)$ moment and $F_{0}+\eta F_{5}$ distribution (last panel) with respect to $\cos\theta_{\Lambda}$. The dots with error bars are the data, and the red line is the weighted PHSP MC corrected by the fitting results. The blue dashed line is the distributions for unweighted simulated PHSP events.

Figure 4 shows the fit result of the $M(\cos\theta_{\Lambda})$ distribution. It is consistent with the behavior described by Eq. (4) when compared to the data. The significance of the spin polarization signal, including systematic uncertainties, is determined to be 5.5$\sigma$ by comparing the likelihoods with and without spin polarization. The moment is expressed as

$$M(\cos\theta_{\Lambda})=-\frac{m}{N}\sum_{i}^{N(m)}(\sin\!\theta_{p}\sin\!\phi_{p}+\sin\!\theta_{\bar{p}}\sin\!\phi_{\bar{p}}),$$ (7)

which is related to the spin polarization and calculated for $m=10$ intervals in $\cos\theta_{\Lambda}$. Here, $N(m)$ denotes the number of events within each $\cos\theta_{\Lambda}$ interval, and $m$ represents the number of bins.

Figure 4: The moments $M(\cos\theta_{\Lambda})$ as a function of $\cos\theta_{\Lambda}$. The dots with error bars are data, and the red line is the weighted PHSP MC corrected by the results of the global fit. The blue dashed line is the distributions from simulated PHSP events.

Systematic uncertainties in the measurement of $\Lambda$ hyperon polarization arise from various sources, including background contributions, $\Lambda$ reconstruction, kinematic fit, beam transverse polarization, decay parameters of $\Lambda\to p\pi$, and the fit method. The background contributions include the sideband region which describes the non-$\Lambda(\bar{\Lambda})$ background and ISR $\psi(3686)$ background. The uncertainty due to background candidates is estimated by comparing the fits with and without the background contributions. To estimate the uncertainty related to the $\Lambda$ reconstruction including the tracking, the requirement on the mass window and decay length of $\Lambda$, it is studied from a control sample of $\psi(3686)\to\Lambda\bar{\Lambda}$ events. The uncertainty arising from this source is evaluated using the same method as in Ref. BESIII:2022lsz . The discrepancy between the nominal and average values obtained from variations is taken as the systematic uncertainty. The uncertainty associated with the kinematic fit is assigned as the results with and without track helix parameter corrections ref5 . The systematic uncertainty originating from the transverse beam polarization is estimated by changing the joint decay angular distribution $\mathcal{W}$ according to Ref. shuliu . The difference of results between the nominal and released parameters of transverse beam polarization is taken as the systematic uncertainty. The uncertainty caused by the fixed decay parameters of $\alpha_{\Lambda/\bar{\Lambda}}$ is estimated by varying mean values obtained from averaging results in BESIII:2022qax within $\pm 1\sigma$. The change of the result is negligible and thereby the related uncertainty is neglected. The reliability of the fit results is validated by performing an input and output check based on 300 pseudoexperiments using the helicity amplitude formula from Ref. BESIII:2022qax . The mean value of polarization parameters measured in the analysis ($\eta$=0.86, $\Delta\Phi$=1.53) are used as input in the formula, and the number of events in each generated MC sample is ten times of the data sample. The difference between the input and output results is taken as the systematic uncertainty. Assuming all sources to be independent, the total systematic uncertainty is calculated as the square root of their quadratic sum. All systematic uncertainties are listed in Table 1.

Table 1: The absolute systematic uncertainties in the measurement of the $\Lambda$ hyperon polarization parameters.

Source	$\eta$	$\Delta\Phi$ (rad)	$R$
Backgrounds	0.02	0.01	0.04
$\Lambda$ reconstruction	0.02	0.01	0.04
Kinematic fit	0.00	0.01	0.00
Beam transverse polarization	0.00	0.02	0.00
Fit method	0.00	0.02	0.00
Total	0.03	0.03	0.05

In summary, we report the observation of spin polarization and complete determination of timelike EMFFs of the $\Lambda$ hyperon in the $e^{+}e^{-}\to{\Lambda}\bar{\Lambda}$ process at $\sqrt{s}=3.773$ GeV, using a data sample corresponding to an integrated luminosity of 20.3 fb^-1 collected by the BESIII detector, offering a higher precision compared with the previous measurement BESIII:2021cvv . The relative phase, the angular distribution parameter and the modulus of the timelike EMFFs ratio are measured to be $\Delta\Phi=(1.53\pm 0.36\pm 0.03)$ rad, $\eta=0.86\pm 0.05\pm 0.03$ and $R=0.47\pm 0.08\pm 0.05$, respectively, where the first uncertainty is statistical and the second systematic. The measured relative phase differs significantly from zero with a significance of 5.5$\sigma$ for the first time at the high momentum transfer region ($s>14$ GeV²), taking into account the systematic uncertainty. A comparison of spin polarization and the modulus of the timelike EMFFs ratio between this work and previous measurements at different c.m. energies BESIII:2019nep ; BESIII:2021cvv ; BESIII:2022qax ; BESIII:2023euh are illustrated in Fig. 5.

Figure 5: Two dimensional distribution of $\Delta\Phi$ and $R$ between this work and previous BESIII measurements at different c.m. energies BESIII:2019nep ; BESIII:2021cvv ; BESIII:2022qax ; BESIII:2023euh . The uncertainty combines from both statistical and systematic uncertainties. The inner, intermediate and outer contours in orange represent 68.2%, 95.4%, and 99.7% confidence level, respectively.

The measured polarization in this work is consistent with and more precise than the previous measurements at $\sqrt{s}=3.773$ GeV BESIII:2022qax , and $\Delta\Phi$ is also roughly consistent with the results at other c.m. energies with an uncertainty of $(1-2)\sigma$. However, noticeable differences exist in the modulus of the timelike EMFFs ratio between this work and other c.m. energies, indicating potential variations in production mechanisms at different c.m. energies.

Spin-1/2 hyperons produced in a hypreon-antihyperon pair can have either the same or opposite helicity. The nonvanishing relative phase $\Delta\Phi$ between the transition amplitudes of these helicity states implies the contributions not only from the $S$-wave but also $D$-wave amplitudes to ${\Lambda}\bar{\Lambda}$ production. Since the modulus of the timelike EMFFs ratio suggests an energy dependence of this value, more data samples at various c.m. energies are needed for a detailed study of the phase dependence on the momentum transfer squared, $s$. The clear and prominent signal enhances our understanding of the ${\Lambda}\bar{\Lambda}$ production mechanism within the $e^{+}e^{-}\to{\gamma/\Psi}\to{\Lambda}\bar{\Lambda}$ process, providing valuable insights into the structure of baryons.

Acknowledgement— The BESIII Collaboration thanks the staff of BEPCII (https://cstr.cn/31109.02.BEPC) and the IHEP computing center for their strong support. This work is supported in part by National Key R&D Program of China under Contracts Nos. 2023YFA1606000, 2023YFA1606704; National Natural Science Foundation of China (NSFC) under Contracts Nos. 12075107, 12247101, 11635010, 11935015, 11935016, 11935018, 12025502, 12035009, 12035013, 12061131003, 12192260, 12192261, 12192262, 12192263, 12192264, 12192265, 12221005, 12225509, 12235017, 12361141819; the Fundamental Research Funds for the Central Universities No. lzujbky-2024-jdzx06; the Natural Science Foundation of Gansu Province (No. 22JR5RA389, No.25JRRA799; by the ‘111 Center’ under Grant No. B20063; the Chinese Academy of Sciences (CAS) Large-Scale Scientific Facility Program; CAS under Contract No. YSBR-101; 100 Talents Program of CAS; The Institute of Nuclear and Particle Physics (INPAC) and Shanghai Key Laboratory for Particle Physics and Cosmology; Agencia Nacional de Investigación y Desarrollo de Chile (ANID), Chile under Contract No. ANID PIA/APOYO AFB230003; German Research Foundation DFG under Contract No. FOR5327; Istituto Nazionale di Fisica Nucleare, Italy; Knut and Alice Wallenberg Foundation under Contracts Nos. 2021.0174, 2021.0299; Ministry of Development of Turkey under Contract No. DPT2006K-120470; National Research Foundation of Korea under Contract No. NRF-2022R1A2C1092335; National Science and Technology fund of Mongolia; National Science Research and Innovation Fund (NSRF) via the Program Management Unit for Human Resources & Institutional Development, Research and Innovation of Thailand under Contract No. B50G670107; Polish National Science Centre under Contract No. 2024/53/B/ST2/00975; Swedish Research Council under Contract No. 2019.04595; U. S. Department of Energy under Contract No. DE-FG02-05ER41374.