Using data samples collected at center of mass energies of $\sqrt{s}$ = 4.009, 4.226, 4.257, 4.358, 4.416 and 4.599 GeV with the BESIII detector operating at the BEPCII storage ring, we search for the isospin violating decay $Y(4260)\rightarrow J/\psi \eta \pi^{0}$. No signal is observed, and upper limits on the cross section $\sigma(e^{+}e^{-}\rightarrow J/\psi \eta \pi^{0})$ at the 90\% confidence level are determined to be 3.6, 1.7, 2.4, 1.4, 0.9 and 1.9 pb, respectively.
Results on $e^{+}e^{-}\rightarrow J/\psi\eta\pi^{0}$. Listed in the table are the integrated luminosity $\cal{L}$, radiative correction factor (1+$\delta^{r}$) taken from QED calculation assuming the $Y(4260)$ cross section follows a Breit$-$Wigner line shape, vacuum polarization factor (1+$\delta^{v}$), average efficiency ($\epsilon^{ee}{\cal B}^{ee}$ + $\epsilon^{\mu\mu}{\cal B}^{\mu\mu}$), number of observed events $N^\text{obs}$, number of estimated background events $N^\text{bkg}$, the efficiency corrected upper limits on the number of signal events $N^\text{up}$, and upper limits on the Born cross section $\sigma^\text{Born}_\text{UL}$ (at the 90 $\%$ C.L.) at each energy point.