Ablikim et al

We study the process $e^{+}e^{-} \rightarrow \pi^{+}\pi^{-}\psi(3686)$ using 5.1 $fb^{-1}$ of data collected at 16 center-of-mass energy ($\sqrt{s}$) points from 4.008 to 4.600$\sim$GeV by the BESIII detector operating at the BEPCII collider. The measured Born cross sections for $e^{+}e^{-}\rightarrow \pi^{+}\pi^{-}\psi(3686)$ are consistent with previous results, but with much improved precision. A fit to the cross section shows contributions from two structures: the first has $M=4209.5\pm7.4\pm1.4$$\sim$MeV/$c^{2}$ and $\Gamma=80.1\pm24.6\pm2.9$$\sim$MeV, and the second has $M=4383.8\pm4.2\pm0.8$$\sim$MeV/c$^{2}$ and $\Gamma=84.2\pm12.5\pm2.1$$\sim$MeV, where the first errors are statistical and the second systematic. The lower-mass resonance is observed in the process $e^{+}e^{-}\rightarrow \pi^{+}\pi^{-}\psi(3686)$ for the first time with a statistical significance of $5.8\sigma$. A charged charmonium-like structure is observed in the $\pi^{\pm}\psi(3686)$ invariant mass spectrum for data at $\sqrt{s} = 4.416~$GeV. A fit with an $S$-wave Breit-Wigner function yields a mass $M=4032.1\pm2.4$~$MeV$^{2}$, where the errors are statistical only. However, there are still unresolved discrepancies between the fit model and data. The width of the intermediate state varies in a wide range for different kinematic regions within the data set. Therefore no simple interpretation of the data has been found, and a future data sample with larger statistics and more theoretical input will be required to better understand this issue.