document: the analysis of modern accelerator - based neutrino oscillation experiments requires good control over the intermediate - energy neutrino - nucleus scattering cross section @xcite . in particular 
   the importance of multi - nucleon events has been suggested in many calculations of charge - changing quasielastic cross sections @xmath2 , at typical neutrino energies of @xmath3 gev  @xcite . 
   the contribution of two - particle - two - hole ( 2p-2h ) excitations is now thought to be essential for a proper description of data  @xcite . 
   thus a growing interest has arisen in including 2p-2h models into the monte carlo event generators used by the neutrino collaborations @xcite . 
   the only 2p-2h model implemented up to date in some of the monte carlo neutrino event generators corresponds to the so - called ific valencia model  @xcite , which has been incorporated in genie  @xcite . 
   there are also plans to incorporate the lyon model  @xcite in genie , while phenomenological approaches like the effective transverse enhancement model of  @xcite are implemented , for instance , in nuwro generator  @xcite . 
   one of the main problems to implementing the 2p-2h models is the high computational time . 
   this is due to the large number of nested integrals involved in the evaluation of the inclusive hadronic tensor with sums over the final 2p-2h states . to speed up the calculations , 
   several approximations can be made , such as choosing an average momentum for the nucleons in the local fermi gas  @xcite , neglecting the exchange matrix elements , or reducing the number of integrations to two nested integrals by performing a non - relativistic expansion of the current operators  @xcite . 
   the latter approach is only useful for some pieces of the elementary 2p-2h response .    in this work we present a fast and very efficient method to calculate the inclusive 2p-2h responses in the relativistic fermi gas model ( rfg ) 
   this approach , denoted as the frozen nucleon approximation , was first explored in @xcite but restricted to the analysis of the 2p-2h phase - space . here 
   it is extended to the evaluation of the full hadronic tensor assuming that the initial momenta of the two struck nucleons can be neglected for high enough energy and momentum transfer , @xmath4 . 
   the frozen nucleon approximation was found to work properly in computing the phase space function for two - particle emission in the range of momentum transfers of interest for neutrino experiments with accelerators . 
   here we investigate the validity of the frozen approximation beyond the phase - space study by including the electroweak meson - exchange current ( mec ) model of @xcite . 
   we find that the presence of virtual delta excitations requires one to introduce a `` frozen '' @xmath1-propagator , designed by a convenient average over the fermi sea . 
   the main advantage of the frozen approximation consists in reducing the number of nested integrals needed to evaluate the inclusive 2p-2h electroweak responses from 7 ( full calculation ) to 1 . 
   thus it is well - suited to computing the 2p-2h neutrino cross sections folded with the neutrino flux , and it can be of great help in order to implement the 2p-2h models in the monte carlo codes currently available . 
   the plan of this work is as follows : in section [ sec_form ] we review the formalism of neutrino scattering and describe mathematically the frozen approximation approach . in section [ sec_results ] 
   we validate the nucleon frozen approximation by computing the 2p-2h response functions and by comparing with the exact calculation . 
   finally , in section [ sec_conclusions ] we summarize our conclusions . 
   the double - differential inclusive @xmath5 or @xmath6 cross section is given by @xmath7 \ 
   , , \end{aligned}\ ] ] where the sign @xmath8 is positive for neutrinos and negative for antineutrinos . 
   the term @xmath9 in eq . 
   ( [ cross ] ) represents the elementary neutrino scattering cross section with a point nucleon , while the @xmath10 are kinematic factors that depend on lepton kinematic variables . 
   their explicit expressions can be found in @xcite . 
   the relevant nuclear physics is contained in the five nuclear response functions @xmath11 , where @xmath12 is the momentum transfer , defining the @xmath13 direction , and @xmath14 is the energy transfer . 
   they are defined as suitable combinations of the hadronic tensor @xmath15    in this work we compute the inclusive hadronic tensor for two - nucleon emission in the relativistic fermi gas , given by @xmath16 where @xmath17 by momentum conservation , @xmath18 is the nucleon mass , @xmath19 is the volume of the system and we have defined the product of step functions @xmath20 with @xmath21 the fermi momentum .    finally the function @xmath22 is the elementary hadron tensor for the 2p-2h transition of a nucleon pair with given initial and final momenta , summed up over spin and isospin , @xmath23 which is written in terms of the antisymmetrized two - body current matrix elements @xmath24 the factor @xmath25 in eq . 
   ( [ elementary ] ) accounts for the antisymmetry of the two - body wave function .    for the inclusive responses considered in this work 
   there is a global axial symmetry , so we can fix the azimuthal angle of one of the particles . 
   we choose @xmath26 , and consequently the integral over @xmath27 gives a factor @xmath28 . 
   furthermore , the energy delta function enables analytical integration over @xmath29 , and so the integral in eq . 
   ( [ hadronic ] ) can be reduced to 7 dimensions ( 7d ) . in the `` exact '' results shown in the next section , this 
   7d integral has been computed numerically using the method described in @xcite . 
   the frozen nucleon approximation consists in assuming that the momenta of the initial nucleons can be neglected for high enough values of the momentum transfer . 
   thus , in the integrand of eq . 
   ( [ hadronic ] ) , we set @xmath30 , and @xmath31 . 
   we roughly expect this approximation to become more accurate as the momentum transfer increases . 
   the integration over @xmath32 is trivially performed and the response function @xmath33 , with @xmath34 , is hence approximated by @xmath35 where @xmath36 and @xmath37 are the elementary response functions for a nucleon pair excitation , which are defined similarly to eqs . 
   ( [ rcc][rtprima ] ) . 
   the integral over @xmath38 can be done analytically by using the delta function for energy conservation , and the integral over @xmath39 gives again a factor of @xmath28 . 
   thus only an integral over the polar angle @xmath40 remains : @xmath41 where the sum runs over the , in general two , possible values of the momentum of the first particle for given emission angle @xmath42 . 
   these are obtained as the positive solutions @xmath43 of the energy conservation equation @xmath44 the explicit values of the solutions of the above equation can be found in the appendix of  @xcite . 
   care is needed in performing the integral over @xmath42 because the denominator inside the integral can be zero for some kinematics . 
   the quadrature in these cases can be done with the methods explained in @xcite . 
   to investigate the validity of the frozen nucleon approximation , we have to choose a specific model for the two - body current matrix elements @xmath45 entering in the elementary 2p-2h response functions , eqs . 
   ( [ elementary],[anti ] ) . 
   here we use the relativistic model of electroweak mec operators developed in @xcite . 
   the mec model can be summarized by the feynman diagrams depicted in fig . 
   [ fig_feynman ] . 
   it comprises several contributions coming from the pion production amplitudes of @xcite . 
   the seagull current , corresponding to diagrams ( a , b ) , is given by the sum of vector and axial - vector pieces @xmath46_{1'2',12 } \frac{f^2_{\pi nn}}{m^2_\pi } \frac{\bar{u}_{s^\prime_1}({{\bf p}}^\prime_1)\,\gamma_5   { \not{\!k}}_{1 } \ , u_{s_1}({{\bf h}}_1)}{k^2_{1}-m^2_\pi}\ ,   \nonumber\\   & \times &    \bar{u}_{s^\prime_2}({{\bf p}}^\prime_2 )   \left [ f^v_1(q^2)\gamma_5 \gamma^\mu   + \frac{f_\rho\left(k_{2}^2\right)}{g_a}\,\gamma^\mu     \right ] u_{s_2}({{\bf h}}_2 )     \nonumber\\     & + &     ( 1\leftrightarrow2 ) \ , , \label{seacur }   \end{aligned}\ ] 
   ] where @xmath47 corresponds to the @xmath8-components of the two - body isovector operator @xmath48 $ ] . 
   the @xmath49 sign refers to neutrino ( antineutrino ) scattering . 
   the four - vector @xmath50 is the momentum carried by the exchanged pion and @xmath51 . 
   the @xmath52 ( @xmath53 ) and axial ( @xmath54 ) couplings , and the form factors ( @xmath55 , @xmath56 ) have been taken from the pion production amplitudes of @xcite . 
   the pion - in - flight current corresponding to diagram ( c ) is purely vector and is given by @xmath57_{1'2',12 }   \frac{f^2_{\pi nn}}{m^2_\pi}\ ,   \frac { f^v_1(q^2 ) \ , \left(k^\mu_{1}-k^\mu_{2}\right)}{\left(k^2_{1}-m^2_\pi\right )   \left(k^2_{2}-m^2_\pi\right)}\ , \nonumber   \\ & \times &   \bar{u}_{s^\prime_1}({{\bf p}}^\prime_1)\,\gamma_5   { \not{\!k}}_{1 } \ 
   , u_{s_1}({{\bf h}}_1 ) \bar{u}_{s^\prime_2}({{\bf p}}^\prime_2)\,\gamma_5   { \not{\!k}}_{2 } \ , u_{s_2}({{\bf h}}_2 ) \ 
   , , \end{aligned}\ ] ] where @xmath58 is the momentum of the pion absorbed by the second nucleon . 
   the pion - pole current corresponds to diagrams ( d , e ) and is purely axial , given by @xmath59_{1'2',12 } \frac{f^2_{\pi nn}}{m^2_\pi}\ , \frac{f_\rho\left(k_{1}^2\right)}{g_a}\;q^\mu \bar{u}_{s^\prime_1}({{\bf p}}^\prime_1){\not{\!q}}u_{s_1}({{\bf h}}_1 ) \nonumber\\ & \times & \frac { \bar{u}_{s^\prime_2}({{\bf p}}^\prime_2)\,\gamma_5 { \not{\!k}}_{2 } \ 
   , u_{s_2}({{\bf h}}_2 )   } { \left(k^2_{2}-m^2_\pi\right)\left(q^2-m^2_\pi\right ) } + ( 1\leftrightarrow2 ) \,.\end{aligned}\ ] ]    finally the @xmath1 current corresponds in fig . 
   [ fig_feynman ] to diagrams ( f , g ) for the forward and ( h , i ) for the backward @xmath1 propagations , respectively . 
   the current matrix elements are given by @xmath60_{1'2',12 } \frac{\bar{u}_{s^\prime_2}({{\bf p}}^\prime_2)\,\gamma_5   { \not{\!k}}_{2 } \ , u_{s_2}({{\bf h}}_2)}{k^2_{2}-m^2_\pi } \nonumber\\ & \times & k^\alpha_{2}\,\bar{u}_{s^\prime_1}({{\bf p}}^\prime_1)g_{\alpha\beta}(h_1+q ) \gamma^{\beta\mu}(h_1,q)u_{s_1}({{\bf h}}_1 ) \nonumber\\ & + & ( 1\leftrightarrow2 ) \label{delta_forward } \\ 
   j^\mu_{\delta,\rm b}&= & \frac{f^ * f_{\pi nn}}{m^2_\pi}\ , \left[u_{\rm b}^{\pm}\right]_{1'2',12}\ ;   \frac{\bar{u}_{s^\prime_2}({{\bf p}}^\prime_2)\,\gamma_5   { \not{\!k}}_{2 } \ , u_{s_2}({{\bf h}}_2)}{k^2_{2}-m^2_\pi } \nonumber\\ & \times & k^\beta_{2}\,\bar{u}_{s^\prime_1}({{\bf p}}^\prime_1 ) 
   \hat{\gamma}^{\mu\alpha}(p^\prime_1,q ) g_{\alpha\beta}(p^\prime_1-q)u_{s_1}({{\bf h}}_1 ) \nonumber\\ & + & ( 1\leftrightarrow2)\label{delta_backward } \ , .\end{aligned}\ ] ] the @xmath61 coupling is @xmath62 . 
   the forward , @xmath63 , and backward , @xmath64 , isospin transition operators have the following cartesian components @xmath65 where @xmath66 and @xmath67 are the isovector transition operators from isospin @xmath68 to @xmath69 or vice - versa , respectively . 
   the @xmath70 operator is for neutrino ( antineutrino ) scattering . 
   the @xmath1-propagator , @xmath71 , is given by @xmath72 where @xmath73 is the projector over spin-@xmath68 on - shell particles , @xmath74\end{aligned}\ ] ] and whose denominator has been obtained from the free propagator for stable particles , @xmath75 , with the replacement @xmath76 to take into account the finite decay width of the @xmath77 . 
   the tensor @xmath78 in the forward current is the weak @xmath79 transition vertex a combination of gamma matrices with vector and axial - vector contributions : @xmath80\gamma_5 \\ \gamma^{\beta\mu}_a(p , q)&= & \frac{c^a_3}{m_n}\left(g^{\beta\mu}{\not{\!q}}- q^\beta\gamma^\mu\right)\nonumber\\ & + & \frac{c^a_4}{m^2_n}\left(g^{\beta\mu}q\cdot \left(p+q\right)- q^\beta \left(p+q\right)^\mu\right ) \nonumber\\ & + & c^a_5 g^{\beta\mu}+\frac{c^a_6}{m^2_n }   q^\beta q^\mu \ , .\end{aligned}\ ] ] for the backward current , we take @xmath81^{\dagger } \gamma^0 \ , .\ ] ] finally , it is worth noting that the form factors @xmath82 are taken from @xcite . 
   we refer to that work for further details of the model . 
   the evaluation of the relevant elementary responses requires one to contract the electroweak two - body mec with themselves by spin - isospin summation . 
   this leads to the squares of each of the diagrams depicted in fig . 
   [ fig_feynman ] plus all their interferences . 
   the validity of the frozen nucleon approximation relies on the fact that the integrand inside the 2p-2h response is a function that depends slowly on the momenta of the two initial nucleons inside the fermi sea . in that case 
   the mean - value theorem applied to the resolution of the integrals provides very precise results . 
   this is so for all of the diagrams of the mec except for the forward @xmath1 diagram , which shows a sharp maximum for kinematics around the @xmath1 peak for pion emission , located at @xmath83 . 
   this is due to the denominator in the @xmath1 propagator , @xmath84 where @xmath85 is the momentum of the hole that gets excited to a @xmath1 .    in these cases 
   the integrand changes very significantly with a small variation of the momentum of the holes and consequently , the frozen approximation can not properly describe the integrand . on the contrary , it only provides a general estimation of the order of magnitude . 
   to get rid of these difficulties we have developed a prescription to deal with the forward @xmath1-propagator appearing in eq . 
   ( [ delta_forward ] ) . 
   this procedure is based on the use of an effective propagator ( `` frozen '' ) for the @xmath1 , conveniently averaged over the fermi gas . 
   this average is an analytical complex function , which is used instead of the `` bare '' propagator inside the frozen approximation , recovering the precision of the rest of diagrams . 
   the `` frozen '' prescription amounts to the replacement : @xmath86 where the frozen denominator is defined by @xmath87    taking the non - relativistic limit for the energies of the holes ( @xmath88 ) , which is justified because hole momenta are below the fermi momentum , itself a value far below the nucleon rest mass , we can write : @xmath89 where @xmath90    assuming the @xmath1 width ( @xmath91 ) to be constant , we can integrate eq . 
   ( [ averaged_prop ] ) over the angles , getting @xmath92 \,.\ ] ] note the complex logarithm inside the integral , which provides the needed kinematical dependence of the averaged propagator , differing from the bare lorentzian shape . finally the integral over the momentum @xmath93 can also be performed , resulting in @xmath94 \right\ } \ , . 
   \nonumber\end{aligned}\ ] ]    by comparing the response functions evaluated in the frozen approximation , _ 
   i.e. , _ substituting the denominator of the @xmath1 propagator in eq . 
   ( [ delta_prop ] ) for the frozen expression in eq . 
   ( [ averaged_denom ] ) , with the exact results , we find that the shapes around the @xmath1 peak are similar , but with slightly different width and position of the center of the peak . 
   we have checked that the differences can be minimized by changing the parameters @xmath95 with respect to the `` bare '' ones , given by eqs . 
   ( [ a],[b ] ) . 
   this is because we have computed the averaged denominator without taking into account the current matrix elements appearing in the exact responses , although the functional form and kinematical dependence is the appropriate one .    in practice , we adjust @xmath91 and apply a shift in the expression for @xmath96 in eq . 
   ( [ a ] ) in order to obtain the best approximation to the exact results . 
   the effective `` frozen '' parameters we actually introduce in eq . 
   ( [ averaged_denom ] ) , are given by @xmath97 we consider @xmath98 and the frozen shift , @xmath99 , to be tunable parameters depending on the momentum transfer @xmath100 . 
   we have adjusted these parameters for different @xmath100-values and we provide them in table  [ table1 ] . 
   .values of the free parameters of the fermi - averaged @xmath1-propagator for different kinematic situations corresponding to different values of the momentum transfer @xmath100 . [ cols="^,^,^",options="header " , ]     -propagator in frozen approximation compared to the average propagator . in this evaluation 
   we have taken @xmath101 mev / c and @xmath102 mev.,width=302 ] 
   in this section we validate the frozen approximation by computing the approximate 2p-2h response functions and comparing with the exact results in the rfg . 
   we consider the case of the nucleus @xmath103c with fermi momentum @xmath104 mev / c , and show the different response functions for low to high values of the momentum transfer . 
   for other nuclei with different @xmath21 the frozen parameters of table 1 should be determined again , and we expect their values change slightly .    in fig . 
   [ fig_denom ] we show the modulus squared of the @xmath1 propagator , given by the @xmath105 function defined in eq . 
   ( [ denominator ] ) , computed for @xmath106 , as a function of @xmath14 for @xmath107 gev / c . 
   it presents the typical lorentzian shape corresponding to width @xmath102 mev . 
   we observe a narrow peak around @xmath108 mev . 
   this corresponds to the @xmath1-peak position for @xmath107 gev / c . in the same figure 
   we also show the square of the frozen average @xmath109 ( solid line ) . 
   the resulting peak is quenched and broadened as compared to the lorentzian shape , reducing its strength and enlarging its width . 
   this behavior of the averaged @xmath1-propagator drives the actual shape of the exact 2p-2h nuclear responses , being more realistic than the simple lorentzian shape of the frozen approximation without the average , as we will see below .     of @xmath103c for different momentum transfers @xmath100 . 
   the exact results are compared to the frozen approximation with and without the averaged @xmath1 propagator . 
   , width=283 ]    in fig . 
   [ fig_rt ] we show the weak transverse 2p-2h response function of @xmath103c for four different values of the momentum transfer ranging from 300 to 1500 mev / c . 
   the curves correspond to different calculations or approximations made in the evaluation of the responses , as labeled in the legend . 
   the solid line corresponds to the seven - dimensional calculation with no approximations . 
   the other two curves refer to the different frozen nucleon approximations developed in this work : the dashed line is obtained with eq . 
   ( [ frozen_eq ] ) but performing the replacement expressed in ( [ replacement_rule ] ) for the forward @xmath1-excitation terms in the evaluation of the current matrix elements ; on the contrary , the dotted line corresponds to the same frozen nucleon approximation , eq . 
   ( [ frozen_eq ] ) , but without the fermi - average of the @xmath1-propagator in the forward terms .    as it can be seen from fig . 
   [ fig_rt ] , for those values of the momentum transfer for which the @xmath1-peak is not reached ( the panel with @xmath110 mev / c ) , there is really little difference between averaging or not the @xmath1 propagator . 
   this is certainly not the case when the @xmath1-peak is fully reached , as shown in the other panels . 
   in this situation there is a dramatic difference between performing the fermi - average of the @xmath1-propagator or not . 
   this difference is in consonance with the results shown in the previous fig . 
   [ fig_denom ] , and it shows how crucial is the treatment of the @xmath1-propagator to obtain accurate results for the 2p-2h responses in the frozen nucleon approximation , i.e , with only one integration .    the results in fig . 
   3 have been obtained after fitting the parameters @xmath111 for the fermi - averaged @xmath1-propagator at the different values of the momentum transfer quoted in table  [ table1 ] . 
   it is also worth noting that there is no way of converting the dotted line into the dashed one by only a suitable fitting of these parameters , _ 
   i.e. , _ without averaging the @xmath1-propagator . 
   c within different models for two values of the momentum transfer . 
   the exact rfg results and the frozen approximation are compared with the shell model ( sm ) results of @xcite . 
   the total shell model results ( 1p-1h ) + ( 2p-2h ) are also shown for comparison . 
   , width=302 ]    c. several values of the momentum transfer are displayed : @xmath112 and 2000 mev / c . 
   , width=302 ]    in fig . 
   [ fig_shellmodel ] we show results for the transverse electromagnetic 2p-2h response function . 
   the frozen and exact ( 7d ) @xmath113 response of the rfg are compared with the results obtained in the shell model 2p-2h calculation of @xcite . 
   this was one of the first computations of the 2p-2h response within the nuclear shell model . 
   the total nuclear response in the shell model , obtained by adding the 1p-1h to the 2p-2h channel , is also shown to appreciate the relative size of the 2p-2h contribution to the total result .    as shown in fig . 
   [ fig_shellmodel ] , the fermi gas results ( either in frozen approximation or not ) are similar to the shell model ones . the small discrepancy between them can not be attributed to relativistic effects because of the low momentum transfer values considered , but to the different coupling constants and form factors used in the model of the @xmath1 meson - exchange current considered in @xcite and the present approach . 
   we can remark the slightly different threshold effects between both calculations . 
   these effects are , as expected , very sensitive to the treatment of the nuclear ground state . 
   note also that the frozen approximation describes reasonably well this low momentum @xmath110 mev / c , considering the simplifications involved .    finally in fig . 
   [ fig5 ] we show that the frozen approximation works notably well in a range of momentum transfer from low to high values of @xmath100 . 
   we compare the @xmath113 , @xmath114 , and @xmath115 2p-2h responses in frozen approximation with the exact results obtained computing numerically the 7d integral of the hadronic tensor . the accord is particularly good for the two transverse responses which dominate the cross section . 
   a slight disagreement occurs for very low energy transfer at threshold where the response functions are anyway small . in the case of the cc response function 
   some tiny differences are observed . 
   however , note that this response is small because the dominant @xmath1 current is predominantly transverse 
   . moreover , its global contribution to the cross section is not very significant because it is partially canceled with the contribution of the cl and ll responses . 
   [ [ physical - interpretation - of - the - frozen - approximation ] ] physical interpretation of the frozen approximation + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    the validity of the frozen approximation led us to conclude that , in the inclusive responses for two - particle emission , the detailed information about the momenta carried out by the two nucleons is lost . 
   this is because the energy and momentum transfer @xmath116 are shared by the two nucleons in multiple ways . 
   this is reminiscent from the phase - space kinematical dependence ( which can be obtained setting the elementary response @xmath117 to unity ) already seen in @xcite . 
   the soft dependence of the elementary response on the initial momenta makes the same argument applicable to the full responses with the exception of the @xmath1 forward current that requires one to soften and average the rapid variations of the @xmath1 propagator . 
   only the low - energy region where the sharing is highly restricted and the cross section is therefore very small , is found to be sensitive to the details of the initial state . 
   this is also supported by the comparison between the shell model and the rfg .    on the other hand , in the 2p-2h model of @xcite , an average momentum @xmath118 
   was determined by imposing quasi - deuteron kinematics . 
   note that this condition is similar to the present frozen approach , but this only guarantees that the total momentum of the two holes is zero , corresponding to selecting back - to - back pair configurations in the ground state only . 
   in this work we have introduced and validated the frozen nucleon approximation for a fast and precise calculation of the inclusive 2p-2h response functions in a relativistic fermi gas model . 
   this approximation neglects the momentum dependence of the two holes in the ground state and requires the use of an effective propagator for the @xmath1 resonance conveniently averaged over the fermi sphere , for which we have provided a simple analytical expression . for momentum transfers above the fermi momentum 
   this approximation makes it possible to compute the responses with only a one - dimensional integral . taking into account all the uncertainties in modeling the two - nucleon emission reactions 
   , this approach can be used instead of the full 7d integral , obtaining very satisfactory results . 
   although we have used a specific model of mec to prove the validity of the approximation , it is reasonable to expect that the frozen approach is also valid for other 2p-2h models . 
   this can be of great interest when implementing 2p-2h models in monte carlo event generators , which up to now have relied on parameterizations from external calculations . in summary , 
   the frozen approximation enables one to make 2p-2h calculations very efficiently and rapidly , instead of interpolating pre - calculated tables , including allowing the parameters of the models to be modified inside the codes , if desired . 
   finally , in the near future this study will be extended to an exploration of how the 2p-2h mec responses depend on nuclear species @xcite . 
   this work has been partially supported by the spanish ministerio de economia y competitividad and erdf ( european regional development fund ) under contracts fis2014 - 59386-p , fis2014 - 53448-c2 - 1 , by the junta de andalucia ( grants no . 
   fqm-225 , fqm160 ) , by the infn under project manybody , and part ( twd ) by the u.s . 
   department of energy under cooperative agreement de - fc02 - 94er40818 . 
   irs acknowledges support from a juan de la cierva fellowship from mineco ( spain ) . 
   gdm acknowledges support from a junta de andalucia fellowship ( fqm7632 , proyectos de excelencia 2011 ) . 
   i.  ruiz simo , j.  e.  amaro , m.  b.  barbaro , a.  de pace , j.  a.  caballero and t.  w.  donnelly , arxiv:1604.08423 [ nucl - th ] . 
   j.  e.  amaro , m.  b.  barbaro , j.  a.  caballero , t.  w.  donnelly and c.  maieron , phys . 
   c * 71 * ( 2005 ) 065501 .
summary: we present a fast and efficient method to compute the inclusive two - particle two - hole ( 2p-2h ) electroweak responses in the neutrino and electron quasielastic inclusive cross sections . 
   the method is based on two approximations . 
   the first neglects the motion of the two initial nucleons below the fermi momentum , which are considered to be at rest . 
   this approximation , which is reasonable for high values of the momentum transfer , turns out also to be quite good for moderate values of the momentum transfer @xmath0 . 
   the second approximation involves using in the `` frozen '' meson - exchange currents ( mec ) an effective @xmath1-propagator averaged over the fermi sea . within the resulting `` frozen nucleon approximation '' , 
   the inclusive 2p-2h responses are accurately calculated with only a one - dimensional integral over the emission angle of one of the final nucleons , thus drastically simplifying the calculation and reducing the computational time . 
   the latter makes this method especially well - suited for implementation in monte carlo neutrino event generators . 
   neutrino scattering , meson - exchange currents , 2p-2h . 
   25.30.pt , 25.40.kv , 24.10.jv