Measurement of cross sections of from = 4$\textit{.}$682 to 4$\textit{.}$951 GeV

The process is studied by analyzing $\textit{e}$ $\textit{e}$ collision data samples collected at eight center-of-mass energies ranging from 4.682 to 4.951 GeV with the BESIII detector operating at the BEPCII collider, corresponding to an integrated luminosity of 4$\textit{.}$1 fb . Observation of the process is found for the first time with a statistical significance of 6$\textit{.}$3$\textit{σ}$, and the cross sections at each center-of-mass energy are measured. The ratio of cross sections of relative to $\textit{e}$ $\textit{e}$ $\textit{→ K}$ $\textit{K}$ $\textit{ψ}$(3686) is determined to be , which is consistent with the prediction based on isospin symmetry. The uncertainty includes both statistical and systematic contributions. Additionally, the invariant mass distribution is found to be consistent with three-body phase space. The significance of a contribution beyond three-body phase space is only 0$\textit{.}$8$\textit{σ}$.