Probing quantum entanglement with generalized parton distributions at the Electron-Ion Collider

Within the collinear factorization framework based on generalized parton distributions, we calculate the spin density matrix of exclusively produced quark and antiquark pairs u u ¯ , d d ¯ , s s ¯ , c c ¯ , b b ¯ in electron-proton scattering. The presence of both real and imaginary parts in the scattering amplitudes leads to a rich pattern of entanglement between the quark and the antiquark. We map out kinematical regions where the pairs exhibit entanglement, Bell nonlocality, and nonstabilizerness. We also predict that massive quarks and antiquarks are transversely polarized, similar to the well-known transverse hyperon polarization in unpolarized collisions. In strangeness, charm, and bottom productions, the polarization can reach 50%–80% in certain kinematic regions in the low-energy runs of the Electron-Ion Collider.