Search for $\eta_{1}(1855)$ in $\chi_{cJ}\to\eta\eta\eta^{\prime}$ decays

Based on a sample of $2.7\times10^{9}$ $\psi(3686)$ events collected by the BESIII detector operating at the BEPCII collider, an analysis of the decay $\psi(3686)\to\gamma\chi_{cJ}, \chi_{cJ}\to\eta\eta\eta^{\prime}$ is performed. The decay modes $\chi_{c1}$ and $\chi_{c2}\to\eta\eta\eta^{\prime}$ are observed for the first time, and their corresponding branching fractions are determined to be $\mathcal{B}(\chi_{c1}\to\eta\eta\eta^{\prime}) = (1.39 \pm 0.13(\text{stat.}) \pm 0.09(\text{sys.})) \times 10^{-4}$ and $\mathcal{B}(\chi_{c2}\to\eta\eta\eta^{\prime}) = (4.42 \pm 0.86(\text{stat.}) \pm 0.37(\text{sys.})) \times 10^{-5}$. An upper limit on the branching fraction of $\chi_{c0}\to\eta\eta\eta^{\prime}$ is set as $2.64 \times 10^{-5}$ at 90\% confidence level (CL). A partial wave analysis (PWA) of the decay $\chi_{c1}\to\eta\eta\eta^{\prime}$ is performed to search for the $1^{-+}$ exotic state $\eta_1(1855)$. The PWA result indicates that the structure in the $\eta\eta^{\prime}$ mass spectrum is mainly attributed to the $f_0(1500)$, while in the $\eta\eta$ mass spectrum, it is primarily the $0^{++}$ phase space. The upper limit of $\mathcal{B}(\chi_{c1}\to\eta_{1}(1855)\eta) \cdot \mathcal{B}(\eta_{1}(1855)\to\eta\eta^{\prime})< 9.79 \times 10^{-5}$ is set based on the PWA at 90\% CL.