We investigate an improved method to extract nucleon matrix elements from lattice 3-point functions using a generalized eigenvalue problem (GEVP) with nucleon and pion-nucleon interpolating fields. Our method avoids the computation of the costly three-point functions that have pion-nucleon interpolators at both source and sink. We demonstrate that excited state contamination from $N\pi$ is minimized in nucleon matrix elements of the scalar, vector, pseudoscalar, axial, and tensor currents and discuss our results based on a physical-point ensemble with a pion mass value of 131 MeV.
We find that the GEVP is most significant for the isovector pseudoscalar and axial currents.