  It has been observed in multiple lattice determinations of isovector axial
and pseudoscalar nucleon form factors, that, despite the fact that the partial
conservation of the axialvector current is fulfilled on the level of
correlation functions, the corresponding relation for form factors (sometimes
called the generalized Goldberger-Treiman relation in the literature) is broken
rather badly. In this work we trace this difference back to excited state
contributions and propose a new projection method that resolves this problem.
We demonstrate the efficacy of this method by computing the axial and
pseudoscalar form factors as well as related quantities on ensembles with two
flavors of improved Wilson fermions using pion masses down to 150 MeV. To this
end, we perform the $z$-expansion with analytically enforced asymptotic
behaviour and extrapolate to the physical point.
