model [@gibuu; @leitner_2009], these are not currently implemented in neutrino interaction generators. In these models, a mean field potential due to the presence of other nucleons within the nucleus is calculated, which will generally depend on the position and momentum of the struck nucleon. These models are not discussed further as they cannot be easily implemented in the [<span style="font-variant:small-caps;">NEUT</span>]{}neutrino interaction generator [@hayato:neut].

Although alternative nuclear models modify the cross section as a function of the outgoing lepton kinematics significantly, they do not change the total CCQE cross section significantly as a function of neutrino energy [@garvey_review_2014]. Additional nuclear effects are also likely to be required to explain the current global dataset. Multinucleon interaction (2p2h) models such as those from Nieves [*et al.*]{} [@nieves; @nievesExtension] and Martini [*et al.*]{} [@martini] go beyond the impulse approximation and include diagrams where two nucleons are involved in the interaction. This adds significant strength to the CCQE-like cross section and explains the difference in normalization observed in the [MiniBooNE]{}data, which was previously modelled with a large axial mass [@nieves_MB_nu_2011; @martini_MB_nu_2011; @nieves_MB_anu_2013; @martini_MB_anu_2013]. Because these 2p2h models are not two-body processes, they are expected to lead to significant biases in the neutrino energy reconstruction from 