Using $448.0 \times 10^6$~$\psi(3686)$ events collected with the BESIII
detector, an amplitude analysis is performed for
$\psi(3686)\to\gamma\chi_{c1}$, $\chi_{c1}\to\eta\pi^+\pi^-$ decays. The most
dominant two-body structure observed is $a_0(980)^{\pm}\pi^{\mp}$;
$a_0(980)^{\pm}\to\eta\pi^{\pm}$. The $a_0(980)$ line shape is modeled using a
dispersion relation, and a significant non-zero $a_0(980)$ coupling to the
$\eta^{\prime}\pi$ channel is measured. We observe $\chi_{c1}\to a_2(1700)\pi$
production for the first time, with a significance larger than 17$\sigma$. The
production of mesons with exotic quantum numbers, $J^{PC}=1^{-+}$, is
investigated, and upper limits for the branching fractions $\chi_{c1}\to
\pi_1(1400)^{\pm}\pi^{\mp}$, $\chi_{c1}\to \pi_1(1600)^{\pm}\pi^{\mp}$, and
$\chi_{c1}\to\pi_1(2015)^{\pm}\pi^{\mp}$, with subsequent $\pi_1(X)^{\pm} \to
\eta\pi^{\pm}$ decay, are determined.

We present a comparative analysis of two models for the $A( {\gamma,\,\pi
N})B$ reaction, which take into account the isobar configurations in the ground
state of the nuclei: the $\Delta N$ correlation model and the quasifree pion
photoproduction model. The considered models differ in their descriptions of
the nucleus states. The $\Delta N$ correlation model takes into account the
dynamic correlations of the nucleon and isobar formed in the virtual transition
$NN \to \Delta N$, and in the quasifree pion photoproduction model, isobars and
nucleons in the nucleus are considered as independent constituents. The
predictions of the models are considered for two reactions ${}^{16}$O$( {\gamma
,\,\pi ^{ +} p} ){}^{15}$C and ${}^{16}$O$( {\gamma ,\,\pi ^{ -} p} ){}^{15}$O.
It is shown, that the two models predict the differential cross section
significantly differing both in absolute values, and in the shape of the
angular dependence. We compare the results of the $\Delta N$ correlation model
for the $( {\gamma ,\,\pi N} )$ and $( {\gamma ,\,\pi NN})$ reactions with the
${}^{16}$O$( {\gamma ,\, \pi ^{ -} p})$ reaction data measured at BNL. Our
results give support to the $\Delta N$ correlation model.