Charge Symmetry Breaking in dd->4He{\pi}0 with WASA-at-COSY

Charge symmetry breaking (CSB) observables are a suitable experimental tool to examine effects induced by quark masses on the nuclear level. Previous high precision data from TRIUMF and IUCF are currently used to develop a consistent description of CSB within the framework of chiral perturbation theory. In this work the experimental studies on the reaction dd->4He{\pi}0 have been extended towards higher excess energies in order to provide information on the contribution of p-waves in the final state. For this, an exclusive measurement has been carried out at a beam momentum of p=1.2 GeV/c using the WASA-at-COSY facility. The total cross section amounts to sigma(tot) = (118 +- 18(stat) +- 13(sys) +- 8(ext)) pb and first data on the differential cross section are consistent with s-wave pion production.

the contribution of p-waves in the final state. For this, an exclusive measurement has been carried out at a beam momentum of p d = 1.2 GeV/c using the WASAat-COSY facility. The total cross section amounts to σ tot = (118 ± 18 stat ± 13 sys ± 8 ext ) pb and first data on the differential cross section are consistent with s-wave pion production.

Introduction
Within the Standard Model there are two sources of isospin violation 8 , namely the electro-magnetic interaction and the differences in the masses of the lightest quarks [1,3]. Especially in situations where one is able to disentangle these two sources, the observation of isospin violation in hadronic reactions is a direct window to quark mass ratios [3,4,5].
The effective field theory for the Standard Model in the MeV range is chiral perturbation theory (ChPT). It maps all symmetries of the Standard Model onto hadronic operators -their strength then needs to be fixed either from experiment or from lattice QCD calculations. At leading order the only parameters are the pion mass and the pion decay constant which are the basis for a series of famous low energy theorems in hadron-hadron scattering (see, for example, Ref. [6]). Although at subleading orders the number of a priori unknown parameters increases, the theory still provides non-trivial links between different operators. A very interesting example is the close link between the quark mass induced proton-neutron mass difference, ∆M qm pn , and, at leading order, isospin violating πN scattering, the Weinberg term. In general, it is difficult to get access to quark mass effects in low energy hadron physics: by far the largest isospin violating effect is the pion mass difference, which also drives the spectacular energy dependence of the π 0 -photoproduction amplitude near threshold (see Ref. [7] and references therein). Thus, it is important to use observables where the pion mass difference does not contribute. An example are charge symmetry breaking (CSB) observables -charge symmetry is an isospin rotation by 180 degrees that exchanges up and down quarks -as the pion mass term is invariant under this rotation. For this case, the impact of soft photons has been studied systematically [8,9,10,11,12] and can be controlled. Already in 1977 Weinberg predicted a huge effect (up to 30% difference in the scattering lengths for pπ 0 and nπ 0 ) of CSB in π 0 N scattering [1] (see also Ref [2] for the recent extraction of these quantities from pionic atoms data).
While the π 0 p scattering length might be measurable in polarized neutral pion photoproduction very near threshold [14], it is not possible to measure the nπ 0 channel. As an alternative access to CSB pion-nucleon scattering in Ref. [13] it was suggested to use N N induced pion production instead. There have been two successful measurements of corresponding CSB observables, namely a measurement of A f b (pn → dπ 0 ) [15] -the forward-backward asymmetry in pn → dπ 0 -as well as of the total cross section of dd → 4 Heπ 0 close to the reaction threshold [16]. The first experiment was analyzed using ChPT in Ref. [17] (see also Ref. [18]), where it was demonstrated that A f b (pn → dπ 0 ) is directly proportional to ∆M qm pn , while the effect of π − η mixing, previously believed to completely dominate this CSB observable [19], was shown to be subleading. The value for ∆M qm pn extracted turned out to be consistent with other, direct calculations of this part based on dispersive analyses [3,20,21] and from lattice. See Ref. [22] for the latest review. In order to cross check the systematics and to eventually reduce the uncertainties, additional experimental information needs to be analyzed.
The first theoretical results for dd → 4 Heπ 0 are presented in [23,24]. The studies show that the relative importance of the various charge symmetry breaking mechanisms is very different compared to pn → dπ 0 . For example, soft photon exchange may significantly enhance the cross sections for dd → 4 Heπ 0 [25]. Furthermore, a significant sensitivity of the results to the nuclear potential model was reported in Ref. [26], which called for a simultaneous analysis of CSB in the N N scattering length and in dd → 4 Heπ 0 [26]. Thus, as part of a consistent investigation of CSB in the two nucleon sector, pn → dπ 0 and dd → 4 Heπ 0 should help to further constrain the relevant CSB mechanisms.
The main challenge in the calculation of dd → 4 Heπ 0 is to get theoretical control over the initial state interactions: high accuracy wave functions are needed for dd → 4N in low partial waves at relatively high energies. One prerequisite to control this is the earlier WASA-at-COSY measurement of dd → 3 Henπ 0 [27], which is allowed by charge symmetry and partially shares the same initial state as dd → 4 Heπ 0 . In addition, higher partial waves are predicted to be very sensitive to the CSB N N → N ∆ transition potential that is difficult to access in other reactions. In leading order in chiral perturbation theory this potential is known. Thus, a measurement of, for example, p-waves provides an additional, non-trivial test of our current understanding of isospin violation in hadronic systems. Future theoretical CSB studies for dd → 4 Heπ 0 can be based on recent developments in effective field theories for few-nucleon systems [28] as well as for the reaction N N → N N π [29,30,31], thus promising a modelindependent analysis of the data.
While the previous measurements of dd → 4 Heπ 0 close to reaction threshold were limited to the total cross section [16], in order to extract constraints on higher partial waves any new measurement at higher excess energies in addition has to provide information on the differential cross section. For this, an exclusive measurement detecting the 4 He ejectile as well as the two decay photons of the π 0 has been carried out utilizing the same setup used for dd → 3 Henπ 0 [27]. The latter reaction was also used for normalization.

Experiment
The experiment was carried out at the Institute for Nuclear Physics of the Forschungszentrum Jülich in Germany using the Cooler Synchrotron COSY [32] together with the WASA detection system [33]. For the measurement of dd → 4 Heπ 0 at an excess energy of Q ≈ 60 MeV a deuteron beam with a momentum of 1.2 GeV/c was scattered on frozen deuterium pellets provided by an internal pellet target. The 4 He ejectile and the two photons from the π 0 decay were detected by the Forward Detector and the Central Detector of the WASA facility, respectively. The experimental setup and trigger conditions were the same as described in Ref. [27].

Data Analysis
The basic analysis leading to event samples with one helium and two photons in final state follows the strategy used for dd → 3 Henπ 0 outlined in Ref. [27]. Compared to this reaction, however, the charge symmetry breaking reaction dd → 4 Heπ 0 has a more than four orders of magnitude smaller cross section. The only other channel with 4 He and two photons in final state is the double radiative capture reaction dd → 4 Heγγ. The cross sections for both reactions are not enough large to provide a visual signature for 4 He in the previously used ∆E − ∆E plots from the Forward Detector. Thus, all 3 He and 4 He candidates together with the two photons have been tested against the hypotheses dd → 4 Heγγ (" 4 He hypothesis") and dd → 3 Henγγ (" 3 He hypothesis") by means of a kinematic fit. Besides the overall energy and momentum conservation no other constraints have been included. Especially, there is no constraint on the invariant mass of the two photons in order to leave a decisive missing-mass plot and not to introduce a fake 4 Heπ 0 signal.
For final event classification the cumulative probabilities P (χ 2 , n.d.f.) for the two hypotheses have been plotted as probability for the 4 He hypothesis versus the probability for the 3 He hypothesis (see Fig. 1). The data (right plot) have been compared to Monte-Carlo generated samples of dd → 4 Heπ 0 events (left plot) and dd → 3 Henπ 0 events (middle plot). from dd → 4 Heπ 0 populate the low probability region for the 3 He hypothesis and form a uniform distribution for the 4 He hypothesis. As there is no pion constraint in the fit, events from the double radiative capture reaction show the same signature. For dd → 3 Henπ 0 the situation is opposite. The indicated cut is based on the Monte-Carlo simulations, but has been optimized by maximizing the statistical significance of the π 0 signal in final missing mass plot. In addition, it has been checked that the result is stable within the statistical errors against variations of the probability cut. For the simulations the standard Geant3 [34] based WASA Monte-Carlo package has been used, which includes the full detector setup and which has already been benchmarked against a wide range of reactions from the WASA-at-COSY physics program. After this analysis step the contribution from misidentified 3 He was reduced by about four orders of magnitude.
In a next step, the resulting four momenta based on the fit hypothesis dd → 4 Heγγ have been used to calculate the missing mass m X in dd → 4 HeX as a function of the center-of-mass scattering angle θ * of the particle X. Figure 2 shows a peak at the pion mass on top of a broad background. In order to extract the number of signal events the background in the peak region has to be described and subtracted. Instead of a (rather arbitrary) fit using a polynomial,  the shape of signal and background has been reproduced using a composition of physics reactions with a double charged nucleus and two photons in the final state. Any further sources of background -physics as well as instrumentalhave already been eliminated by the analysis steps described in Ref. [27] and the subsequent kinematic fit. The signal has then been extracted by fitting a linear combination of the corresponding Monte-Carlo generated high-statistics template distributions for the three reactions • dd → 4 Heγγ (double radiative capture) using 3-body phase space (green dashed), plus • dd → 3 Henπ 0 using the model described in Ref. [27] (blue dotted) for which the 3 He is falsely identified as 4 He, plus • dd → 4 Heπ 0 using 2-body phase space (i.e. plain s-wave, red solid).
Please note that in Fig. 2 as well as in Fig. 3 the comulated distributions are shown, e.g. the red solid curve represents the sum of all contributions. For the differential cross section the data have been divided into four angular bins within the detector acceptance (−0.85 ≤ cos θ * ≤ 0.75). Independent fits of the different contributions listed above have been performed for each bin to address possible anisotropies. In the course of the fit two systematic effects have been observed, which are discussed in the following.
First, the background originating from misidentified 3 He is slightly shifted compared to the Monte-Carlo simulations. The effect is angular dependent and is largest at forward angles. Possible reasons are a mismatch in the actual beam momentum, a different amount of insensitive material in Monte-Carlo simulations compared to the real experiment or systematic differences in the simulated detector response for 3 He and 4 He -the limited statistics did not allow for a detailed study of the origin of that effect. The background stemming from dd → 3 Henπ 0 is sensitive to these effects as the energy losses from a (true) 3 He ejectile are used for energy reconstruction of a (falsely identified) 4 He. The mismatch can be compensated by introducing an angular dependent scaling factor on the missing mass axis for the 3 Henπ 0 background, which has been included in the fit as additional free parameter. For the angular bins from backward to forward these factors are 1.0, 0.99, 0.97 and 0.94, respectively. As the resulting fits describe the shape of the data especially in the region of the pion peak, no additional systematic error has been assigned to this effect.
The second systematic effect concerns a mismatch in the low mass range m ≤ 0.11 GeV/c 2 in the most backward angular bin. According to the fit only events from the reaction dd → 4 Heγγ contribute in this mass region. The model used for this channel was 3-body phase space, which was not expected to provide a perfect description. However, with the dominating background from dd → 3 Henπ 0 in a wide mass range, it is currently not possible to disentagle the two contributions precisely enough in order to verify any more advanced theoretical model -this issue will be addressed in a follow-up experiment, see below. Consequently, the final fit excludes the corresponding missing mass range (consistently in all angular bins). Based on the difference to the fit with the low mass region included a corresponding systematic uncertainty for this effect has been assigned in the result. Figure 3 shows the fitted missing mass spectra for the different bins in cos θ * together with the fit result. The chosen ansatz provides a good overall description of the full data set. Any tests for further systematic effects (according to the definition in Ref. [35]), for example concerning rate effects and selection cuts in the basic analysis (see Ref. [27]), did not reveal any additional systematic uncertainties.

Results
For the acceptance correction an isotropic angular distribution has been assumed. For absolute normalization the reaction dd → 3 Henπ 0 has been used. The resulting differential cross sections extracted from In general, only statistical errors are given, except for the first bin where the uncertainty caused by the systematic effect in the low mass region has been included. A systematic error of 10% for luminosity determination and 7% for the normalization to external data is common to all numbers. Integrating the individual results, the (partial) total cross section within the detector acceptance amounts to σ acc tot = (94 ± 14 stat ± 10 sys ± 6 ext ) pb (5) with the systematic error originating from luminosity determination and the uncertainty from the different fit methods. The external normalization error has been propagated from the luminosity determination for dd → 3 Henπ 0 (see Ref. [27]). Extrapolation to full phase space by assuming an isotropic distribution yields σ tot = (118 ± 18 stat ± 13 sys ± 8 ext ) pb.
This result can be compared with the values measured close to threshold by dividing out phase space (see Fig. 4). A constant value could be interpreted as a dominating s-wave, but one has to keep in mind that the energy dependence of the formation of a 4 He in the 4N final state might have some influence here, too. Figure 5 shows the differential cross section. Due to the identical particles in the initial state, odd and even partial waves do not interfere and the angular distribution is symmetric with respect to cos θ * = 0. As the p-wave and s-d interference terms contribute to the quadratic term and the p-wave also adds to the constant term, the different partial waves cannot be directly disentangled. However, a fit including the Legendre polynomials P 0 (cos θ * ) and P 2 (cos θ * ) -although not excluding -does not show any evidence for contributions of higher partial waves: dσ dΩ = (9.8 ± 2.6) pb/sr · P 0 (cos θ * ) + (9.5 ± 7.4) pb/sr · P 2 (cos θ * ).
Here, the two coefficients are strongly correlated with the correlation parameter 0.85, i.e. the reader should not interpret the two contributions as independent results. Based on the fit results a first estimate of the total cross section of dd → 4 Heγγ has been extracted assuming a homogeneous 3-body phase space. It amounts to σ tot = (0.92 ± 0.07 stat ± 0.10 sys ± 0.07 norm ) nb.
It should be noted that this result depends on the underlying models for the reactions dd → 3 Henπ 0 and dd → 4 Heγγ. This model dependence is not included In the first bin the additional systematic uncertainty from the fit has been added (see text). The blue dashed line represents the total cross section given in the text assuming an isotropic distribution, the solid red curve shows the fit with the Legendre polynomials P 0 and P 2 .
in the given systematic error.

Summary and Conclusions
In this letter results were presented for a measurement of the charge symmetry breaking reaction dd → 4 Heπ 0 at an excess energy of 60 MeV. The energy dependence of the square of the production amplitude might indicate the on-set of higher partial waves or some unusual energy dependence of the s-wave amplitude -given the current statistical error, no conclusion on the strength of the higher partial waves is possible from the differential cross section.
However, since within chiral perturbation theory the leading and next-toleading p-wave contribution does not introduce any new free parameter (it is expected to be dominated by the Delta-isobar), the data on the strength of higher partial waves presented in this work will still provide a non-trivial constraint for future theoretical analyses.
The results presented here are based on a two-week run using the standard WASA-at-COSY setup. Based on the experiences gained during this experiment another 8 week measurement with a modified detector setup optimized for a time-of-flight measurement of the forward going ejectiles has been performed recently. In total, an increase of statistics by nearly a factor of 10 and significantly reduced systematic uncertainties can be expected. In particular, the experiment has been designed to provide a better discrimination of background events from dd → 3 Henπ 0 .