Center-of-mass effects in electromagnetic two-proton knockout reactions Carlotta Giusti,1 Franco Davide Pacati,1 Michael Schwamb,2,3 and Sigfrido Boffi1 1Dipartimento di Fisica Nucleare e Teorica, Universit`a degli Studi di Pavia Istituto Nazionale di Fisica Nucleare, Sezione di Pavia, I-27100 Pavia, Italy 2 Dipartimento di Fisica, Universit`a degli Studi di Trento and Istituto Nazionale di Fisica Nucleare, Gruppo Collegato di Trento, I-38100 Povo (Trento), Italy 3 European Center for Theoretical Studies in Nuclear Physics and Related Areas (ECT∗), I-38100 Villazzano (Trento), Italy The role of center-of mass (CM) effects in the one-body nuclear current in the description of electromagnetically induced two-nucleon knockout reactions is discussed in connection with the problem of the lack of orthogonality between initial bound states and final scattering states obtained by the use of an energy-dependent optical model potential. Results for the cross sections of the exclusive 16O(e,e′pp)14C and 16O(γ,pp)14C knockout reactions in different kinematics are presented and discussed. In super-parallel kinematics CM effects produce a strong enhancement of 7 the16O(e,e′pp)14Cg.s. crosssection whichstrongly reducesthedestructiveinterferencebetweenthe 0 one-bodyand∆-currentandthesensitivitytothetreatmentofthe∆-currentfoundinpreviouswork. 0 2 PACS numbers: 21.60-n Nuclear structuremodels and methods - 25.20Lj Photoproduction reac- n tions - 25.30Fj Inelastic electron scattering to continuum a J 9 I. INTRODUCTION 1 The investigation of nuclear structure is one of the most important and ambitious aims of hadronic physics. A 1 reasonablestarting point is offeredby the independent particle shell model. However,the incorporationof additional v 4 short-range correlations (SRC) beyond a mean-field description turns out to be inevitably necessary for a proper 5 description of nuclear binding. The most direct reaction to study SRC is naturally electromagnetically induced 0 two-nucleon knockout. Intuitively, the probability that a photon is absorbed by a nucleon pair should be a direct 1 measure for SRC. However, due to competing two-body effects like meson-exchange currents (MEC) or final state 0 interactions(FSI),thissimplepictureneedstobemodifiedinordertoobtainquantitativepredictions. Ideally,therole 7 of MEC and FSI should be small or at least under control in order to extract information on SRC from experiment. 0 / This requires a theoretical approach which should be as comprehensive as possible. An overview over the available h theoretical models till the middle of the 90s can be found in [1]. Presently, different models are available (see [2–4] t - andreferencestherein). Duetotheconceptualcomplexityofthenuclearmany-bodyproblem,variousapproximations l c and simplifying assumptions are needed for practicalcalculations. Thus, usually different treatments of initial bound u and final scattering states are adopted in the models. n In the Pavia model [4, 5] bound and scattering states are, in principle, consistently derived as eigenfunctions of an : v energy-dependentnon-HermitianFeshbach-typeopticalpotential. However,inactualcalculationstheinitialhadronic i stateisobtainedfromarecentcalculationofthetwo-nucleonspectralfunction[6]wheredifferenttypesofcorrelations X areincludedconsistently. Forthe finalhadronicstate,acomplexphenomenologicalopticalpotential,derivedthrough r a fit to nucleon-nucleus scattering data, is used for the description of the FSI between the outgoing nucleons and the a residual nucleus. The mutual nucleon-nucleon interaction (NN-FSI) in the final state can be taken into account at least perturbatively [7, 8]. Independently of the specific prescriptions adopted in the calculations, a conceptual problem arises in the model where the initial and final states, which are eigenfunctions of an energy-dependent optical potential at different energies, are, as such, not orthogonal. Indeed, the process involves transitions between bound and continuum states which must be orthogonal, since they are eigenfunctions of the full nuclear many-body Hamiltionian at different energies. Orthogonality is in general lost in a model when the descripion is restricted to a subspace where other channels are suppressed. The description of direct knockout reactions in terms of the eigenfunctions of a complex energy-dependent optical potential considers only partially the contribution of competing inelastic channels. The remainingeffectsduetooccuringinelasticitiescan,inprinciple,betakenintoaccountbyasuitableeffectivetransition operator, which removes the orthogonality defect of the model wave functions [9]. In practice, however, the usual approach does not make use of an effective operator. The present paper deals with the proper treatment of all the CM effects in the matrix elements of the one-body nuclear current in connection with the problem of the lack of orthogonality between initial and final states in the calculation of the cross section of the electromagnetic two-nucleon knockout reactions. The reaction mechanism and CM effects are discussed in sect. II. Different prescriptions are proposed to cure the spuriosity which may result in the numerical calculations as a consequence of the orthogonality defect. These 2 prescriptionsarediscussedandrelatedtoapropertreatmentofalltheCMeffectsinthetransitionmatrixelements. In sect. III the effects of CM and orthogonality are illustrated, with specific numerical examples in selected kinematics, for the exclusive 16O(e,e′pp)14C and 16O(γ,pp)14C knockout reactions. A summary and some conclusions can be found in sect. IV. II. REACTION MECHANISM AND CENTER-OF-MASS EFFECTS Thebasicingredientsforthecalculationofthecrosssectionofthereactioninducedbyarealorvirtualphoton,with momentum q, where two nucleons, with momenta p′, and p′, are ejected from a nucleus, are given by the transition 1 2 matrix elements of the charge-currentdensity operator between initial and final nuclear states Jµ(q)= <Ψ |Jˆµ(r)|Ψ >eiq·rdr. (1) f i Z Bilinear products of these integrals give the components of the hadron tensor, whose suitable combinations allow the calculation of all the observables available from the reaction process [1, 5]. If the residualnucleus is left in a discrete eigenstate of its Hamiltonian, i.e. for an exclusive process,andunder the assumption of the direct knock-out mechanism, the matrix elements of Eq. (1) can be written as [5, 10]1 Jµ(q)= ψ∗(r ,r )Jµ(r,r ,r )ψ(r ,r )eiq·rdrdr dr . (2) Z f 1 2 1 2 i 1 2 1 2 e e Eq. (2) contains three main ingredients: the two-nucleon scattering wave function ψ , the nuclear current Jµ and f the two-nucleon overlap integral (TOF) ψ between the ground state of the target and the final state of the residual i e nucleus. e The nuclear current Jµ is the sum of a one-body and a two-body contribution, i.e. Jµ(r,r ,r )=J(1)µ(r,r )+J(1)µ(r,r )+J(2)µ(r,r ,r ). (3) 1 2 1 2 1 2 The one-body (OB) part includes the longitudinal charge term and the transverse convective and spin currents, and can be written as J(1)µ(r,r )=j(1)µ(r,σ )δ(r−r ) (4) k k k with k = 1,2. The two-body current is derived from the effective Lagrangian of [11], performing a non relativistic reduction of the lowest-order Feynman diagrams with one-pion exchange. We have thus currents corresponding to the seagull and pion-in-flight diagrams, and to the diagramswith intermediate ∆-isobar configurations [12], i.e. J(2)(r,r ,r )=Jsea(r,r ,r )+ Jπ(r,r ,r )+J∆(r,r ,r ). (5) 1 2 1 2 1 2 1 2 Fortwo-protonemissionthe seagullandpion-in-flightmeson-exchangecurrentsandthecharge-exchangecontribution ofthe ∆-currentarevanishinginthe nonrelativisticlimit. Thesurvivingcomponentsofthe ∆-currentcanbe written as J(2)µ(r,r ,r )=j(2)µ(r ,σ ,τ )δ(r−r )+j(2)µ(r ,σ ,τ )δ(r−r ). (6) 1 2 12 2 2 1 12 1 1 2 with r = r −r . Details of the nuclear current components can be found in [4, 12–14]. More specifically, the 12 1 2 various treatments and parametrizations of the ∆-current used in the calculations are given in [4]. In orderto evaluate the transition amplitude of Eq. (2), forthe three-body system consisting of the two protons, 1 and 2, and of the residual nucleus B, it appears to be natural to work with CM coordinates [5, 15] r =r −r , r =r −r , 1B 1 B 2B 2 B A r = r /(A−2). (7) B i Xi=3 1 Spin/isospinindicesaregenerallysuppressedintheformulasofthispaperforthesakeofsimplicity. 3 The conjugated momenta are given by A−1 1 1 p = p′ − p′ − p , (8) 1B A 1 A 2 A B 1 A−1 1 p = − p′ + p′ − p , (9) 2B A 1 A 2 A B P = p′ +p′ +p , (10) 1 2 B where p =q−p′ −p′ is the momentum of the residual nucleus in the laboratory frame. B 1 2 With the help of these relations, one can cast the transition amplitude (2) into the following form Jµ(q)= ψ∗(r ,r )Vµ(r ,r )ψ(r ,r )dr dr , (11) Z f 1B 2B 1B 2B i 1B 2B 1B 2B with the definition ψ (r ,r ):=ψ (r ,r ) (12) i/f 1B 2B i/f 1 2 and the expression e A−1 1 Vµ(r ,r )=exp iq r exp −iq r (j(1)µ(r ,σ )+j(2)µ(r ,σ ,τ )) +(1↔2) . (13) 1B 2B 1B 2B 1B 1 12 2 2 (cid:18) A (cid:19) (cid:18) A (cid:19) It is generally thought that the contribution of the OB current to two-nucleon knockout is entirely due to the correlations included in the two-nucleon wave function. In fact, a OB operator cannot affect two particles if they are not correlated. It can be seen, however, from Eq. (13) that in the CM frame the transition operator becomes a two-body operator even in the case of a OB nuclear current. Only in the limit A → ∞ CM effects are neglected and the expression in Eq. (11) vanishes for a pure OB current in Eq. (13) sandwiched between orthogonalized single particle (s.p.) wave functions. This means that, due to this CM effect, for finite nuclei the OB current can give a contribution to the cross section of two-particle emission independently of correlations. This effect is similar to the one of the effective charges in electromagnetic reactions [16]. ThematrixelementsofEq.(11)involveboundandscatteringstates,ψ andψ ,whichareconsistentlyderivedfrom i f anenergy-dependentnon-HermitianFeshbach-type Hamiltonianfor the consideredfinal state of the residualnucleus. They are eigenfunctions of this Hamiltonian at negative and positive energy values [1, 5]. However, in practice, it is notpossible to achievethis consistency andthe treatmentofinitialandfinalstatesproceedsseparatelywithdifferent approximations. The two-nucleon overlap function (TOF) ψ contains information on nuclear structure and correlations. Different i approachesare used in [6, 10, 17, 18]. In the presentcalculations the TOF is obtained as in [6], from the most recent calculation of the two-proton spectral function of 16O, where both SRC and long-range correlations are included consistently with a two-step procedure. In the calculation of the final-state wave function ψ only the interaction of each one of the two outgoing nucleons f with the residual nucleus is included. Therefore, the scattering state is written as the product of two uncoupled s.p. distorted wave functions, eigenfunctions of a complex phenomenological optical potential which contains a central, a Coulomb, and a spin-orbit term [19]. The effect of the mutual interaction between the two outgoing nucleons has been studied in [7, 8, 20] and can in principle be included as in [7, 8]. The matrix element of Eq. (11) contains a spurious contribution since it does not vanish when the transition operator V is set equal to 1. This is essentially due to the lack of orthogonality between the initial and the final state wave functions. In the model the use of an effective nuclear current operator removes the orthogonality defect besides taking into account space truncation effects [1, 9]. In the usual approach of Eq. (11), however, the effective operatorisreplacedbythebarenuclearcurrentoperator. Thus,itisthisreplacementwhichmayintroduceaspurious contributionwhichisnotspecificallyduetothedifferentprescriptionsadoptedinpracticalcalculations,butisalready presentinEq.(11),whereψ andψ areeigenfunctionsofanenergy-dependentFeshbach-typeHamiltonianatdifferent i f energies. In the past this spuriousity was cured by subtracting from the transition amplitude the contribution of the OB current without correlations in the nuclear wave functions. In detail, the expression A−1 1 ψ∗(r ,r ) exp iq r exp −iq r j(1)µ(r ,σ )+1↔2 ψ (r ,r )dr dr (14) Z f 1B 2B (cid:18) (cid:18) A 1B(cid:19) (cid:18) A 2B(cid:19) 1B 1 (cid:19) i,noCor 1B 2B 1B 2B wassubtractedfrom(11),whereintheinitialstateψ SRCareignored. Thisprescriptionisdenotedasapproach i,noCor A in the proceeding discussions. In this approach, however, we do not subtract only the spuriosity, but also the CM 4 effectgivenbythetwo-bodyoperatorinEq.(13),whichispresentintheOBcurrentindependentlyofcorrelationsand which is not spurious. The relevance of this effect can be estimated comparing our previous results with the results of a different prescription, that is denoted as approach B in the proceeding discussions, where we subtract from (11) instead of (14) the spurious contribution due to the OB current without correlations and without CM corrections. This can be achieved by putting the limit A→∞ in (14), i.e. by the expression ψ∗(r ,r ) exp(iqr )j(1)µ(r ,σ )+1↔2 ψ (r ,r )dr dr . (15) Z f 1B 2B (cid:16) 1B 1B 1 (cid:17) i,noCor 1B 2B 1B 2B This prescription gives an improved, although still rough, evaluation of the spurious contribution. An alternative and more accurate procedure to get rid of the spuriosity is to enforce orthogonality between the initial and final states by means of a Gram-Schmidt orthogonalization [21]. In this approach each one of the two s.p. distorted wave functions is orthogonalized to all the s.p. shell-model wave functions that are used to calculate the TOF, i.e., for the TOF of [6], to the h.o. states of the basis used in the calculation of the spectral function, which range from the 0s up to the 1p-0f shell. This more accurate procedure, that we denote as approach C in the proceedingdiscussions,allowsustogetridofthespuriouscontributiontotwo-nucleonemissionduetoaOBoperator acting on either nucleon of an uncorrelated pair, which is due to the lack of orthogonality between the s.p. bound and scattering states of the pair. In this approach, in consequence, no OB current contribution without correlations like (14) or (15) needs to be subtracted. Moreover,it allows us to include automatically all CM effects via (13). III. RESULTS The effects ofCMandorthogonalizationhavebeeninvestigatedforthe exclusive16O(e,e′pp)14C and16O(γ,pp)14C reactions. Calculationsperformedindifferentsituationsindicatethattheresultsdependonkinematicsandontheprescriptions adopted to treat the theoretical ingredients of the model. The contribution due to the CM effects in the OB current withoutcorrelations,thatwereneglectedinourpreviouscalculations,areingeneralnonnegligible. Althoughinmany situations this contribution is small and does not change significantly the results, there are also situations where it is largeandproducesimportantquantitativeandqualitativedifferences. Thisisthecaseofthesuper-parallelkinematics, where these effects are maximized. The super-parallel kinematics is therefore of particular interest for our study. In the so-called super-parallel kinematics the two nucleons are ejected parallel and anti-parallel to the momentum transferand,forafixedvalueoftheenergyω andmomentumtransferq,itispossibletoexplore,fordifferentvaluesof the kinetic energiesofthe outgoingnucleons,allpossiblevalues ofthe recoilmomentump . This kinematicalsetting B has been widely investigated in our previous work [4–8, 10, 17] and is of particular interest from the experimental point of view, since it has been realized in the recent 16O(e,e′pp)14C [22] and 16O(e,e′pn)14N [23] experiments at MAMI. The super-parallel kinematics chosen for the present calculations of the 16O(e,e′pp)14C reaction is the same alreadyconsideredinourpreviousworkandrealizedintheexperiment[22]atMAMI,i.e. theincidentelectronenergy is E =855 MeV, ω =215 MeV, and q =316 MeV/c. 0 The cross section of the 16O(e,e′pp)14C reaction to the 0+ ground state of 14C calculated in the super-parallel kinematics is displayed in Fig. 1 for the three different approaches A (dotted), B (dashed) and C (solid). The CM contributionincludedinapproachBproducesalargeenhancementofthecrosssectioncalculatedwiththeOBcurrent. The results are shown in the right panel of the figure, where it can be seen that the enhancement is large for recoil momentum values up to about 300 MeV/c and is a factor of about 5 in the maximum region. A similar result is obtained with orthogonalizedinitial and final states. In this case the OB cross section is a bit larger at low values of the recoil momentum and a bit lower at larger values of p . B The results depicted by the dashed and solidlines in Fig. 1 correspondto the two different proceduresproposedto cure the spuriosity due to the lack orthogonalitybetween initial and final states in the model. In the dashed line the spurious contribution is subtracted, in the solid line orthogonality between the s.p. states is restored. With respect to the previous result, shown by the dotted line, the dashed line includes all the CM effects, the solid line takes into account,inaddition,alsotheeffectduetothelackoforthogonality. Itcanbeclearlyseenfromthecomparisonshown in the figure that the large difference between the old and the new results is mostly due to the CM effects and not to thetreatmentofthespuriosityortotherestorationoforthogonalitybetweentheinitialandfinalstatewavefunctions of the model. The final cross sections given by the sum of the OB and the two-body ∆-currents are compared in the left panel of Fig. 1. Calculations have been performed with the so-called ∆(NN) parametrization [4] for the ∆-current, i.e. the parameters have been fixed considering the NN-scattering in the ∆-region, where a reasonable description of data is achieved with parameters similar to the ones of the full Bonn potential [24]. It was shown and explained in [4] that in the super-parallel kinematics, and for the transition to the ground state of 14C, a regularized prescription for 5 FIG. 1: The differential cross section of the 16O(e,e′pp)14C reaction to the 0+ ground state of 14C as a function of pB in a super-parallelkinematicswithE0=855MeV,electronscatteringangleθe =18◦,ω=215MeV,andq=316MeV/c. Different values of pB are obtained changing the kinetic energies of the outgoing nucleons. Positive (negative) values of pB refer to situations where p is parallel (anti-parallel) to q. The final results given by sum of the one-body and ∆-currents (OB+∆) B are displayed in the left panel, the separate contribution of the one-body (OB) current is shown in the right panel. The TOF from thetwo-proton spectral function of [6]andthe∆(NN)parametrization [4]for the∆-currentare usedin thecalculations. Thedottedlinesgivetheresultsof[4],i.e.ofapproachA.Thedashedandsolid linesrefertoapproach B(improvedtreatment of the CM contribution of the OB current) and approach C (explicit orthogonalization of s.p. bound and scattering states), respectively. FIG.2: Thedifferentialcrosssectionofthe16O(e,e′pp)14Cg.s. reactionasafunctionofpB inthesamesuper-parallelkinematics as in Fig. 1. Calculations are performed with approach C. TOF as in Fig. 1. OB (dotted line), OB+∆(NoReg) (dashed line), OB+∆(NN)(solid line). the ∆, such as, e.g., ∆(NN), produces a destructive interference with the OB current which makes the final cross section lower than the OB one. This reduction is strong in our previous calculation of [4], up to about one order of magnitude. The relevance of the destructive interference depends, however, on the relative weight of the OB and ∆-currentcontributions[4]. ThestrongenhancementoftheOBcurrentcontributionproducedinthenewcalculations by CM effects reduces the destructive interference between the OB and ∆-currents. Thus, only a slight reduction of the OBcurrentcontributionisobtainedinthe new calculationsbythe addtionalincorporationofthe ∆-current. The final cross section is completely dominated by the OB current and at low values of the recoil momentum it is more than one order of magnitude larger than the one obtained in the old calculations. An enhancement factor of about 30 is given in the maximum region. It was shown in [4] that dramatic differences are found in the super-parallel kinematics with different parametri- zations of the ∆-current and with different TOFs. ThecrosssectionscalculatedinapproachCfordifferent∆-parametrizationsareshowninFig.2. Theresultwiththe regularized∆(NN)prescription,alreadyshowninFig.1,iscomparedwiththeonegivenbythesimplerunregularized 6 FIG.3: Thedifferentialcrosssectionofthe16O(e,e′pp)14Cg.s. reactionasafunctionofpB inthesamesuper-parallelkinematics asinFig.1. CalculationsareperformedwiththeTOFfromthesimplerapproachof[10]andwiththe∆(NN)parametrization. The dotted line gives the result of approach A [4], thesolid line is obtained with approach C. approach∆(NoReg)of[4]. In[4]thefinalcrosssectionscalculatedwiththesetwoparametrizationsdifferuptoabout one order of magnitude. Only small differences are obtained in Fig. 2. The cross section with ∆(NN) is a bit lower and the one with ∆(NoReg) a bit higher than the cross section given by the OB current. We note that the orthogonalized wave functions are used also in the calculation of the matrix elements with the ∆-current, where the effect of orthogonalization is anyhow negligible. In practice, in the present calculations the contribution of the ∆-current is the same as in [4]. Thus, the large difference with respect to the results of [4] in Figs. 1 and 2 is due to the CM effects in the OB current and, as a consequence, to the strong reduction of the destructive interference betweenthe OB and the ∆-currentcontributioncalculatedwith the ∆(NN) parametrization. ThecrosssectionsshowninFig.3arecalculatedwiththesimplerTOFof[10],wherethetwo-nucleonwavefunction is given by the product of a coupled and antisymmetrized shell model pair function and of a Jastrow-type central and state independent correlation function, taken from [25]. In this approach only SRC are considered and the final state of the residual nucleus is a pure two-hole state. The ground state of 14C is a (p )−2 hole in 16O. Thus, in the 1/2 orthogonalizedcalculation the s.p. distorted wave functions are orthogonalizedonly to the p state. 1/2 The differences between the results of approaches A and C, which are shown in Fig. 3, are significant, although less dramatic than those with the TOF from the spectral function displayed in Fig. 1, and do not change the main qualitative features of the previous results. It can be noted that in Fig. 3 the differences are largerat largervalues of the recoilmomentum, i.e. in the kinematical regionwhere the differences between the correspondingresults in Fig. 1 are strongly reduced. Thus, the CM effects included in the present calculations drastically reduce the sensitivity to the treatment of the ∆-currentfoundin[4]forthesuper-parallelkinematics. TheseCMeffectsare,however,verysensitivetothetreatment of the TOF. The large differences given in the orthogonalized approachC by the two TOFs in Figs. 1 and 3 confirm that the cross sections are very sensitive to the treatment of correlations in the TOF. This result strongly motivates further research,both from the experimental as well as from the theoretical side, in the field of pp-knockout. Similarcalculationsperformedforthetransitiontothe1+ excitedstateof14Cdonotshowanysignificantdifference with respect to our previous results shown in [4]. The effect of the mutual interaction between the two outgoing protons (NN-FSI) has been neglected in the cal- culations presented till now because it is not relevant to investigate CM effects. NN-FSI has been studied within a perturbative treatment in [7, 8], where it is found that the effect depends on the kinematics and on the type of reaction considered. Since NN-FSI turns out to be particularly strong just for the 16O(e,e′pp)14C reaction and in g.s. the super-parallelkinematics,it canbe interestingto give hereonly one numericalexample forthis case,just to show how our previous results of [8] change in the orthogonalized approach. The effect of the NN-FSI on the cross section of the 16O(e,e′pp)14C reaction in the super-parallel kinematics is g.s. showninFig.4. Theresultsobtainedintheapproachconsideredtillnow(DW),whereonlytheinteractionofeachone of the outgoing nucleons with the residual nucleus is considered, are compared with the results of the more complete treatment(DW-NN)wherealsothemutualinteractionbetweenthetwooutgoingnucleonsisincludedwithinthesame perturbative approach as in [8]. The cross sections given by the separate contributions of the OB and ∆-current, as well as the ones given by the sum OB+∆, are displayed in the figure. These results can be compared with the 7 FIG.4: Thedifferentialcrosssectionofthe16O(e,e′pp)14Cg.s. reactionasafunctionofpB inthesamesuper-parallelkinematics asinFig.1. CalculationsareperformedinapproachCwiththesameTOFand∆-parametrizationasinFig.1. Lineconvention for the left panel: OB+∆ with DW-NN (solid line), OB+∆ with DW (dashed line). Line convention for the right panel: OB with DW-NN(solid line), OBwith DW (dotted line), ∆-current with DW-NN(dashed line), ∆-current with DW (dot-dashed line). corresponding ones presented in Figs. 3 and 4 of [8], which differ not only because the calculations of Fig. 4 are performed with orthogonalizedinitial and final states, but also because a different ∆-parametrizationand a different TOF are used in the two calculations. In fact, the ∆(NN) parametrization is used in Fig. 4 compared to an old prescription of ours in [8]. The TOF of [6] is used in Fig. 4 and the one obtained from the first calculation of the spectral function of [17] in [8]. The different treatment of the various theoretical ingredients produce significant numericaldifferencesinthecalculatedcrosssections. ThecontributionofNN-FSItothefinalcrosssectionis,however, of the same type and of about the same relevance as in [8]. In particular, the considerable enhancement given by NN-FSI for medium and large values of the recoil momentum is confirmed in the present calculations. In contrast to the results of [8], where the enhancement at large p is due the ∆-current contribution, in Fig. 4 it is essentially due B to the OB current, which is always dominant in the cross section for all the values of the recoil momentum. AdifferentkinematicalsituationisconsideredinFig. 5. The16O(e,e′pp)14C and16O(γ,pp)14C crosssections g.s. g.s. are calculated in a coplanar symmetrical kinematics where the two nucleons are ejected at equal energies and equal but opposite angles with respect to the momentum transfer. In this kinematical setting different values of p are B obtained changing the scattering angles of the two outgoing protons. The 16O(e,e′pp)14C cross section displayed in the top panel is calculated with E = 855 MeV, θ = 18◦, and g.s. 0 e ω =215 MeV, i.e. the same values as in the super-parallel kinematics. In this symmetrical kinematics, however, the CM effect included in the orthogonalized approach gives only small differences with respect to the previous result. The cross section is dominated by the OB current and, as in [4], it is not affected by the treatment of the ∆-current. It is, however, sensitive to the treatment of correlations in the TOF [4]. The cross section of the 16O(γ,pp)14C reaction at an incident photon energy E = 400 MeV, displayed in the g.s. γ bottompanelofFig.5,isdominatedbythe∆-current. Thus,theeffectsduetoCMandorthogonalizationincludedin thepresentcalculations,whichmainlyaffecttheOBcurrent,donotaffectthe finalcrosssectionforrecoilmomentum valuesupto∼200MeV/c. Athighervaluesofp ,wherethe∆-currentcontributionisstronglyreduced,theseeffects B produce a large increase of the contribution of the OB current and of the final cross section. In the region where the ∆-current is dominant the sensitivity of the results to the ∆-parametrizationis the same as in [4], i.e. very large. IV. SUMMARY AND CONCLUSIONS Two basic aspects have been discussed within the frame of electromagnetic two-proton knockout reactions, i.e. CM effects and the spuriosity arising from the lacking orthogonality between initial and final state wave functions in connection with the usual treatments of the nuclear current. They have been investigated for the cross sections of the exclusive 16O(e,e′pp)14C and 16O(γ,pp)14C reactions under the traditional conditions of super-parallel and symmetrical kinematics. Different kinematics and transitions to discrete low-lying states of the residual nucleus are known to emphasize either the role of the one-body currents, and thus of correlations,or of the two-body ∆-current. Since in two-nucleon knockout one is primarily interested in studying correlations, it is important to keep all the ingredients of the cross section under control in order to extract the useful information from data. 8 FIG. 5: The differential cross section of the 16O(e,e′pp)14Cg.s. (top panel) and 16O(γ,pp)14Cg.s. (bottom panel) reactions as a function of pB in a coplanar symmetrical kinematics with E0 =855 MeV, θe =18◦, ω = 215 MeV, and q = 316 MeV/c (top panel),Eγ =400MeV(bottompanel). DifferentvaluesofpB areobtainedchangingthescatterimg anglesofthetwooutgoing protons. Positive(negative)valuesofpB refertosituationswherepB isparallel(anti-parallel)toq. Calculationsareperformed intheDWapproachandwiththeTOFand∆-parametrization asinFig.1. ThedottedlinesaretheresultsofapproachA[4], thesolid lines are obtained with approach C. Inourpreviouscalculationsoftwo-nucleonknockoutnotalltheCMeffectswereproperlytakenintoaccount. Inthe CM frame the transitionoperator becomes a two-body operatoreven in the case of a one-body nuclear current. As a consequence,the one-body currentcangivea contributionto the crosssectionoftwo-particleemissionindependently of correlations. This effect is similar to the one of the effective charges in electromagnetic reactions [16]. The effective transition operator entering the transition matrix element is in principle defined consistently with the two-body initial and final state wave functions derived from an energy-dependent non-Hermitian Hamiltonian. In such an approach, no spurious contribution comes from the orthogonality defect of the wave functions [1, 9]. In practice, however, one approximates the transition operator in terms of simple forms of one- and two-body currents, thusintroducingsomespuriosity. Inthepast,thisspuriousitywascuredbysubtractingfromthetransitionamplitude the contribution of the one-body current without correlations in the nuclear wave functions. In this way, however, not only the spuriosity is subtracted, but also the CM effect given by the two-body operator which is present in the one-bodycurrentindependentlyofcorrelationsandwhichisnotspurious. Alternatively,onecanenforceorthogonality between the initial and final states by means of a Gram-Schmidt orthogonalization. The two approaches have been investigated here and shown to give similar results. However, the Gram-Schmidt orthogonalization has been further used in the present investigation because it is preferable in principle and allows us to naturally include all the CM effects. TheCMeffectsduetotheone-bodycurrentwithoutcorrelationsaredifferentindifferentsituationsandkinematics. For the 16O(e,e′pp)14C reaction in the super-parallel kinematics these CM effects produce a strong enhancement g.s. of the contribution of the one-body current. As a consequence, the destructive interference between the one-body and the two-body ∆-current as well as the sensitivity to the treatment of the ∆-current discussed in [4] are strongly reduced. With respecttothe resultsof[4], thecalculatedcrosssectionisenhancedandseemsto betterreproducethe experimentaldataof[22]. Onthe other hand,these CM effects areverysensitive to the treatmentofthe two-nucleon overlap function describing the initial correlated pair of protons. The mutual interaction between the two emerging protons produces a large enhancement of the cross section at medium andlargerecoilmomenta. The effect is ofthe same type andofaboutthe samerelevance asin[8]. However, incontrastto the resultsof[8], wherethe enhancementatlargep is due the ∆-currentcontribution,whenincluding B 9 CM effects it is essentially due to the one-body current, which is always dominant in the super-parallel cross section for all the values of the recoil momentum. In the symmetrical kinematics the 16O(e,e′pp)14C reaction is dominated by the one-body current and is thus g.s. sensitivetothetreatmentofcorrelations,confirmingtheresultfoundin[4]. Incontrast,the16O(γ,pp)14C reaction g.s. is dominated by the ∆-current and is not affected by CM and orthogonalization effects up to recoil momenta of the order of 200 MeV/c. In conclusion, the CM effects investigated in this work depend on kinematics and on the final state of the residual nucleus. Thenumericalexamplesshowninthepresentanalysisindicatethattheseeffectsareparticularlylargeforthe 16O(e,e′pp)14C reaction in the super-parallel kinematics of the MAMI experiment. The extreme sensitivity to the g.s. treatment of the different ingredients of the model and to different effects and contributions makes the super-parallel kinematics very interesting but also not particularly suitable to disentangle and investigate the specific contribution of short-range correlations. More and different situations should be considered to achieve this goal. Two examples have been shown in the symmetricalkinematics where either correlationsor the ∆-current are dominant. In order to disentangle and investigate the different ingredients contributing to the cross sections, experimental data are needed in different kinematics which mutually supplement each other. TheinvestigationofCMeffects andorthogonalitybetweeninitialandfinalstateswillbeextendedinaforthcoming paperto the caseofelectromagneticproton-neutronknockout,asurgentlyneededafterthe recentfirstmeasurements of the 16O(e,e′pn)14C reaction performed at the MAMI microtron in Mainz [23]. g.s. [1] S.Boffi,C.Giusti,F.D.Pacati,M.Radici,Electromagnetic ResponseofAtomicNuclei,OxfordStudiesinNuclearPhysics (Clarendon Press, Oxford, 1996); S. Boffi, C. Giusti, F. D. Pacati, Phys. Rep.226, 1 (1993). [2] M. Anguiano, G. Co’, A.M. Lallena, Nucl. Phys. A 744, 168 (2004). [3] J. Ryckebusch,W. van Nespen, Eur. Phys.J. A 20, 435 (2004). [4] C. Giusti, F. D.Pacati, M. Schwamb,S. Boffi, Eur.Phys. J. A 26, 209 (2005). [5] C. Giusti, F. D.Pacati, Nucl.Phys. A 535, 573 (1991). [6] C. Barbieri, C. Giusti, F. D. Pacati, W. H. Dickhoff,Phys. Rev. C 70, 014606 (2004) [7] M. Schwamb,S. Boffi, C. Giusti, F. D.Pacati, Eur.Phys. J. A 17, 7 (2003). [8] M. Schwamb,S. Boffi, C. Giusti, F. D.Pacati, Eur.Phys. J. A 20, 233 (2004). [9] S.Boffi, F. Cannata, F. Capuzzi, C. Giusti, F. D.Pacati, Nucl. PhysA 379, 509 (1982). [10] C. Giusti, F. D.Pacati, Nucl.Phys. A 615, 373 (1997). [11] R.D. Peccei, Phys.Rev. 176, 1812 (1968); 181, 1902 (1969). [12] C. Giusti, F. D.Pacati, Nucl.Phys. A 641, 297 (1998). [13] P.Wilhelm, H. Arenh¨ovel,C. Giusti, F.D. Pacati, Z. Phys. A 359, 467 (1997). [14] C. Giusti, F. D.Pacati, Phys.Rev.C 67, 044601 (2003). [15] D.F. Jackson, T. Berggren, Nucl.Phys. 62, 353 (1965). [16] S.Boffi, C. Giusti, F. D. Pacati, M. Rosa-Clot, NuovoCimento 76A, 685 (1983). [17] C. Giusti, F. D.Pacati, K.Allaart, W. J. W. Geurts, W. H. Dickhoff,H. M. Mu¨ther, Phys.Rev.C 57, 1691 (1998). [18] D.N. Kadrev,M. V.Ivanov,A.N. Antonov,C. Giusti, F. D. Pacati, Phys. Rev.C 68, 014617 (2003). [19] A.Nadasen et al., Phys. Rev.C23, 1023 (1981). [20] D.Kn¨odler, H. Mu¨ther, P. Czerski, Phys. Rev.C 61, 064603 (2000). [21] S.Boffi, R.Cenni, C. Giusti, F. D.Pacati, Nucl.PhysA 420, 38 (1984). [22] G. Rosner, Prog. Part. Nucl. Phys. 44, 99 (2000). [23] D.G. Middleton et al., Eur. Phys.J. A 29, 261 (2006). [24] R.Machleidt, K.Holinde, Ch. Elster, Phys. Rep. 149, 1 (1987). [25] C. C. Gearhart, Ph.D thesis, Washington University,St.Louis (1994); C. C. Gearhart, W. H. Dickhoff,privatecommunication.