Observation of a neutral charmoniumlike state $Z_c(4025)^0$ in $e^{+} e^{-} to (D^{*} bar{D}^{*})^{0} pi^0$

We report a study of the process $e^{+} e^{-} to (D^{*} bar{D}^{*})^{0} pi^0$ using $e^+e^-$ collision data samples with integrated luminosities of $1092 rm{pb}^{-1}$ at $sqrt{s}=4.23 rm{GeV}$ and $826 rm{pb}^{-1}$ at $sqrt{s}=4.26 rm{GeV}$ collected with the BESIII detector at the BEPCII storage ring. We observe a new neutral structure near the $(D^{*} bar{D}^{*})^{0}$ mass threshold in the $pi^0$ recoil mass spectrum, which we denote as $Z_{c}(4025)^{0}$. Assuming a Breit-Wigner line shape, its pole mass and pole width are determined to be $(4025.5^{+2.0}_{-4.7}pm3.1) rm{MeV}/c^2$ and $(23.0pm 6.0pm 1.0) rm{MeV}$, respectively. The Born cross sections of $e^{+}e^{-}to Z_{c}(4025)^{0} pi^0to (D^{*} bar{D}^{*})^{0}pi^0$ are measured to be $(61.6pm8.2pm9.0) rm{pb}$ at $sqrt{s}=4.23 rm{GeV}$ and $(43.4pm8.0pm5.4) rm{pb}$ at $sqrt{s}=4.26 rm{GeV}$. The first uncertainties are statistical and the second are systematic.