Observation of $J/\psi$ decays to $e^{+}e^{-}e^{+}e^{-}$ and $e^{+}e^{-}\mu^{+}\mu^{-}$

Using a data sample of $4.481\times 10^8$ $\psi(3686)$ events collected with the BESIII detector, we report the first observation of the four-lepton-decays $J/\psi\to e^+e^-e^+e^-$ and $J/\psi\to e^+e^-\mu^+\mu^-$ utilizing the process $\psi(3686)\to \pi^+\pi^- J/\psi$. The branching fractions are determined to be $[4.32\pm0.26~(\rm stat)\pm0.19~(\rm syst)]\times 10^{-5}$ and $[2.45~\pm0.21~(\rm stat)\pm0.10~(\rm syst)]\times 10^{-5}$, respectively. The results deviate from theoretical predictions calculated assuming lepton flavor universality, by 2.8 and 5.2 $\sigma$, respectively. No significant signal is observed for $J/\psi\to \mu^+\mu^-\mu^+\mu^-$, and an upper limit on the branching fraction is set at $1.6\times 10^{-6}$ at the 90$\%$ confidence level. A CP asymmetry observable is constructed for the first two channels, which is measured to be $(-0.019\pm0.067\pm0.025)$ and $(-0.016\pm0.081\pm0.003)$, respectively. No evidence for CP violation is observed in this process.