The MORA project

The MORA (Matter's Origin from the RadioActivity of trapped and oriented ions) project aims at measuring with unprecedented precision the D correlation in the nuclear beta decay of trapped and oriented ions. The D correlation offers the possibility to search for new CP-violating interactions, complementary to searches done at the LHC and with Electric Dipole Moments. Technically, MORA uses an innovative in-trap orientation method which combines the high trapping efficiency of a transparent Paul trap with laser orientation techniques. The trapping, detection, and laser setups are under development, for first tests at the Accelerator laboratory, JYFL, in the coming years.


Introduction
Why are we living in a world of matter? What is the reason for the strong matter antimatter asymmetry we observe in the Universe?
The MORA (Matter's Origin from the RadioActivity of trapped and oriented ions) project aims at searching for possible hints. In 1967, A. Sakharov expressed the 3 conditions which should be fulfilled for the baryogenesis process [1]. These conditions are: (i) a large C and CP violation; (ii) a violation of the baryonic number, (iii) a process out of thermal equilibrium. A large CP violation has therefore to be discovered to account for this large matterantimatter asymmetry, at a level beyond the CP violation predicted to occur in the Standard Model via the quark-mixing mechanism.
MORA aims at measuring the D correlation [2] in the nuclear beta decay of trapped and oriented ions with unprecedented precision. There is a large window in which D correlation and Electric Dipole Moment (EDM) measurements can contribute to the search for other sources of CP violation at a much higher level than predicted by the Standard Model [3]. The D correlation offers the possibility to search for new CP-violating interactions in a region that is less accessible by EDM searches, in particular via the Leptoquark model [4,5,6]. Leptoquarks are heavy bosons coupling leptons and quarks which appear in many extensions of the Standard Model. Their decay out of thermal equilibrium in the framework of Grand Unified Theories came rather evidently as the first hypothesis for the baryon and lepton number non conservation [7,8,9]. As such, Leptoquarks still play a peculiar role among the many models now formulated for baryogenesis. Remarkably, they were lately invoked to explain recent lepton flavor non conservation in the decay of B mesons [3]. With sensitivity on D close to 10 −5 , the MORA apparatus should additionally permit a probe of the Final State Interactions (FSI) effects for the first time. The latter can mimic a non-zero D correlation of the order of 10 −5 to 10 −4 , depending on the beta decay transition which is observed. Technically, MORA uses an innovative in-trap orientation method which combines the high trapping efficiency of a transparent Paul trap with laser orientation techniques. The MORA setup is schematically shown in Fig. 1. The project will first focus on the proof-ofprinciple of the in-trap laser orientation technique, before the measurement of the D correlation in the decay of 23 Mg + ions is undertaken firstly at JYFL and then later, at GANIL.
The measurement of correlations appearing in the nuclear beta decay spectrum provides stringent constraints on non standard model interactions (see Figure 1: In trap optical polarization and detection setup. The polarization is obtained by superimposing the laser beam with the injected ion beam. Helmoltz coils maintain a preferential axis for the B-field in the trap. The Helmoltz coils used to cancel transverse stray B fields in the azimuthal plane are not shown for more visibility. The polarization is monitored by a difference of count rates in the two annular Si detectors on the axis of the trap. The D correlation is measured in the azimuthal plane of the trap by an emiT like [26] arrangement of beta and recoil ion detectors, consisting respectively of phoswhich and Micro-Channel Plate (MCP) detectors. e.g. [11,12] and ref therein). The competitiveness of such measurements with respect to LHC searches for new boson interactions has recently been investigated (see e.g. [11,12,1]) using an Effective Field Theory (EFT) framework, which encompasses a wide class of New Physics (NP) models. Bounds for exotic vector, axial-vector, tensor and scalar interactions were compared. For some of the nonstandard interactions LHC searches provide limits that are competitive with (or sometimes much stronger than) beta decays. In contrast, the LHC cannot offer competitive bounds for the specific CP-violating interaction affecting the D correlation, which thus presents an excellent NP potential. For future reference, the correlations in the decay of oriented nuclei are listed in Tab. 1 together with their dependence on the kinematics parameters of the decay and their intrinsic symmetries: m e , p e , E e are the electron mass, momentum, and total energy, p ν , E ν the momentum and energy of the neutrino, and J the spin of the nucleus. P stands for the parity transformation, T for time reversal. The correlations with capital letters refer to nuclei of oriented spin, while the Fierz interference term b and the β − ν angular correlation parameter a βν are also encountered in decays of non-oriented nuclei. Other types of correlations appear when measuring the polarization of electrons or positrons, like the triple R Table 1: Correlations in β decay for oriented nuclei [2]. See text for details.

Name of the correlation Correlation term Symmetry
Beta neutrino angular correlation a βν · pe Ee · pν Eν P even T even Fierz interference term b · me 3 Laser polarization of the 23 Mg + ions So far, trap-based experiments aiming at measuring observables with polarized nuclei have been achieved with atoms using Magneto-Optical-Traps (MOTs), sometimes in combination with other types of atom traps, see e.g. [14] and [15]. In [14] the degree of polarization is determined from the photo-ionization of atoms in a non-polarized state, while in [15] a sophisticated method of retrapping atoms in a MOT was used. In both setups a degree of polarization above 99% was achieved and determined with a precision in the per-mille range.
In contrast to radioactive atom trapping in MOTs, radioactive ion trapping can be an extremely efficient process (more than 10%), especially since the advent of ion cooling techniques in RFQ cooler bunchers [16]. The degree of polarization can be expected to be equally high (close to 100%) because of the extended exposure of the confined ions to the laser light. The proof-of-principle of the laser polarization will be performed at JYFL at the IGISOL-4 facility [17], with a pulsed laser system. The outlined combination of an ion trap and laser polarization is a new concept, taking advantage of the high trapping efficiency of the transparent Paul trap and the high degree of polarization from optical pumping. This method has never been used for such correlation measurements in nuclear β decay, although laser pumping has recently been applied in a RFQ cooler-buncher for collinear laser spectroscopy [18], and in linear Paul traps, where ion laser cooling allows the study of quantum phase transitions in crystals of a few trapped ions or the development of new generation atomic clocks [19]. We will first focus on 23 Mg + , which is a suitable candidate for both laser pumping and β decay correlation measurements. 39 Ca + will be considered at a later stage. The trap will be inspired from the LPCTrap design [20], but installed into a dedicated chamber with ports for the laser beams, and allowing a different detection configuration for probing the polarization degree of the trapped ions, as shown in Fig. 1. A RFQ cooler buncher is readily available in the IGISOL-4 beamlines [17] and will provide bunched and cooled beams of 23 Mg for efficient trapping. The optical pumping schemes to polarize Mg ions are well established [21,22,23]. Recent simulations indicate that the optimal laser system will consist of a single broadband pulsed Ti:Sa laser, whose frequency will be tripled to cover the transitions leading to the maximum magnetic substates M F = ±2 shown in the left inset of Fig. 2. Such a laser system is readily available at IGISOL-4. The laser beam will be circularly polarized using a linear polarizer and quarter wave plate, and guided via suitable transport optics to the trap. The simulations take into account the Doppler effect due to the thermal motion of trapped ions. Assuming a typical repetition rate of 10 kHz and 20 µJ per pulse (0.2 W average power) the achievable nuclear polarization degree after 1 ms (10 pulses) of light exposure would essentially be limited by the degree to which the laser light could be circularly polarized, saturating at a level above 99%. Similar results are obtained after 10 ms for a lower power of 5 µJ per pulse (see right inset of Fig. 2). Ions are cooled via collisions with the He buffer gas, at a typical pressure of 10 −5 mbar. Owing to the expected low collision rate with the most abundant 4 He atoms (kHz), for which no depolarization in collisions with the ground state (S=I=J=0) is expected, possible depolarization due to collisions with the residual gas between the pulses will be marginal. Special care will therefore be taken to monitor the degree of polarization of the laser light by standard optical methods. The polarization will be easily inverted by changing the light circular polarization from σ+ to σ−. Helmoltz coils placed outside the vacuum chamber will maintain a B field of the order of a fraction of a mT directed along the axis of the trap to prevent depolarization. The polarization degree will be determined by measuring an asymmetry in the counts recorded by two annular β detectors placed along and against the axis of the trap as shown in the middle inset of Figure 1, taking the form: where N ↑ β + and N ↓ β + are the numbers of positrons detected parallel and antiparallel to the polarization direction, respectively, and P the polarization degree. The Standard Model value of A β is tabulated for some of the mirror isotopes in [24] and is known to 0.3% for 23 Mg: A β = −0.5584 ± 0.0017. The use of coincidences with recoils measured in the azimuthal plane of the trap, with the recoil ion detectors serving also for the D correlation measurement (right inset of Fig. 1 and section 4) will a priori be preferred to single β events, as coincidences permit a high rejection of background [25]. The coincidence with the recoil ions will induce constraints on the neutrino momentum of the detected events, because of the law of momentum conservation. In the case of β's conditioned by the coincident detection of a recoil ion, the asymmetry of Eq. 1 can be reformulated as a linear combination of the β and neutrino asymmetry coefficients: Here N ↑ β + coinc and N ↓ β + coinc are the numbers of coincidences between recoil ions detected in the azimuthal plane of the trap, and positrons detected parallel and anti-parallel to the polarization axis in the two annular detectors, and α and β are coefficients depending on the detection solid angle and decay parameters. These coefficients can be precisely determined by a Monte Carlo simulation, and B ν is known to 0.54% precision: B ν = −0.7404 ± 0.0040 [24]. With this method the precision aimed at in the determination of the degree of polarization will be of the order of 2%, which is more than sufficient for looking for a non-zero D correlation. The systematic uncertainties arising in the measurement of P from non-uniform detection efficiencies will be efficiently reduced by regularly inverting the direction of polarization. Eq. (2) shows that in first approximation the relative accuracy on P will scale as the square root of the number of β recoil coincidences recorded in the different detectors: Other optical methods (e.g. monitoring the fluorescence signal induced by a laser probe) for the precise determination of the degree of polarization will additionally be investigated. An independent measurement of the degree of polarization would for instance enable the precise measurement of the A β and B ν correlations (see section 2), which could then be used for other tests of the Standard Model, for example the determination of V ud from mirror transitions [12].

D correlation measurement
The D correlation will be maximal in the detection plane of MORA shown in Fig. 1, the product p e × p ν being equivalent to p e × p r , where p r is the ion recoil momentum, because of momentum conservation. The magnitude of the correlation as a function of the angle between the electron and the recoil θ er in the plane perpendicular to the polarization axis is shown in the upper inset of Fig. 3. It is maximal for an angle approaching 165 • , similar to the neutron decay. The detector arrangement will therefore be close to the one of the emiT experiment [26], alternating every 45 • a detector for the recoil ion and a detector for the β particle, as shown in Fig. 1. Such arrangement is magnified in Fig.  4. The lower inset of Fig. 3 shows the D correlation as would be observed by such a detection setup, assuming an octagonal arrangement of 6 cm diameter recoil ion and positron detectors placed at 10 cm from the trap center. Given the polarization direction, the D correlation can be inferred from an asymmetry in the number of coincidences recorded at average θ er angles of +45 • , +135 • on one side and -45 • , -135 • on the other side, the sign of θ er being defined clockwise with respect to the spin direction (see Fig. 4).
N +45 • coinc , N +135 • coinc , N −45 • coinc , N −135 • coinc refer to the number of recorded coincidences at a given averageθ er , and δ is a constant depending on the detection solid angle and decay parameters. δ can be precisely determined by a Monte Carlo simulation of the β decay spectrum and detection setup. In the conditions assumed in this first simulation (point-like ion cloud, no RF disturbance of the recoil, no detection threshold) it is equal to 0.775(1) and has to be compared to the factor |K| ∼ 0.3 for the emiT experiment quoted in [26]. A more realistic simulation will not give a large difference in the δ value (not more than a few percent), as the listed approximations correspond to second order effects: the size of the ion cloud is of the order of 1 mm of radius in comparison to the 4 cm radius detectors. The RF disturbance can be corrected by recording the RF phase as presently done for the measurement of a βν . As is presently the case for the a βν measurement, low detection thresholds will be achievable for both the electron detectors and the MCPs detecting the recoiling ions, which will be equipped with grids to post-accelerate the ions from a few tens of eV to a several keV Figure 3: Contribution of the D correlation to the β decay spectrum of 23 Mg + ion. θ er is the angle between the electron and recoil in the plane perpendicular to the polarization axis. The D coefficient is arbitrarily magnified to 1 for more visibility. Black line: integral β decay spectrum as a function of θ er , which contains all additional relevant correlations: a βν , A β and B ν from [24]. Blue dotted line: same, without the contribution of the D correlation. Red dotted line: contribution of the D correlation. Upper inset: no phase space cut. Bottom inset: same but filtered by the detection setup as shown in Fig. 1 and 4. to reach the efficiency plateau [27]. The sizeable gain of sensitivity compared to the emiT experiment is attributed to the selection of the plane of the decay, perpendicular to the polarization axis, which becomes possible because of the confinement of the emitting source (ions in our case) in a small volume in the center of the trap. Systematic uncertainties arising from non-uniform detection efficiencies and polarization will be reduced by regularly inverting the direction of polarization, and using the symmetries of the octagonal arrangement. A more elaborate method, such as the ones described in [26,28], for inferring limits on D while cancelling the possible contributions of other correlations caused by these non-uniformities will be developped for the purpose of MORA. In experiments aiming at setting limits on a non-zero D correlation, the precise knowledge of P is obviously not required, as long as P can be maintained over 70-80%. P will be monitored during the D correlation measurement using coincidences of the recoiling ions with βs detected in the annular Si detectors, as described in section 3. The sensitivity on a non -zero D correlation (eq. 4) will scale as: Despite its apparent simplicity Eq. 5 is in agreement with the statistical considerations given in e.g. [26]. The expected rates are of the order of 5 · 10 8 pps for 23 Mg from the upgraded SPIRAL facility. An overall ion beam cooling, injection and trapping efficiency of the order of few tens of % can reasonably be achieved, provided that the trap is not filled to its maximum space charge capacity, which is presently of the order of 5 · 10 5 ions [29]. One believes that this capacity could be enlarged by one order of magnitude with dedicated R&D. DESIR, which includes a High Resolution Separator and a number of purification traps, will be the ideal place for providing a clean beam for MORA. It will host a laser setup (LUMIERE) presently missing in the low energy beam lines of GANIL for the polarization of 23 Mg + ions. In the case of JYFL, the statistics  23 Mg could be produced with 2 µA of p, while for safe operation, the primary proton beam would be essentially limited to 10 µA. With a 23 Na target, the second reaction resulted in a level of contamination from 23 Na at least one order of magnitude higher than for the first one: a few 10 8 pps compared to 10 7 pps, respectively. Even by choosing the appropriate target, the level of contamination would be very difficult to handle in the IGISOL RFQ beam cooler used to prepare the 23 Mg bunches prior to their injection into the transparent Paul trap of MORA. With a typical capacity of 5 · 10 5 ions/bunch, the number of 23 Mg injected in the trap would be limited to about 10 3 ions per cycle. We plan to reduce the contamination by an order of magnitude by using a new dedicated RF sextupole ion guide for the extraction of beams from IGISOL, which would allow one to trap 10 4 23 Mg + ions per cycle. The use of the future MR-ToF-MS of IGISOL will additionally be considered to further purify the beam. Using equations 3 and 5, and the results of the Monte Carlo simulation of the β spectrum, the relative accuracy on P and sensitivity to a non -zero D can be estimated for the different phases of MORA. Table 2 summarizes the results of such estimates. A typical efficiency for the recoil-ion detector (Micro-Channel Plates) of 50% was assumed, as well as realistic detection solid angles obtained from latest designs of the trap and detection chamber (see Fig. 5. Three weeks of beam time, which could be organized in successive campaigns at DESIR, should be sufficient to reach a sensitivity well below 10 −4 .

Sensitivity to New Physics
In beta decay, the triple correlations can only be motion-reversal odd. As such, the triple D and R correlations are, under CPT invariance, sensitive to new sources of CP violation [30]. At higher energy, CP violation has been observed in the decay of the kaon and B mesons. CP violation is incorporated in the Standard Model via the quark mixing mechanism, but at a level which cannot account for the large matter antimatter asymmetry observed in the Universe. D and R are sensitive probes which offer a large window to search for a much larger CP violation. This window will eventually be limited by the level at which the effects of electromagnetic FSI, which mimic the existence of non-zero D and R correlations, can be calculated [1]. The FSI effects have never been observed so far at the present level of precision on D and R. Such effects are typically of the order of 10 −4 to 10 −5 depending on the system, and have been recently calculated to an absolute accuracy of 10 −7 for the neutron [31]. D F SI values of the order of −1.3 · 10 −4 and 0.47 · 10 −4 can be estimated using for instance [32] for 23 Mg and 39 Ca, respectively. The bounds on imaginary Tviolating scalar and tensor interactions are obtained from measurements of the R correlation in the decay of the neutron and 8 Li respectively [11,12]. Stronger bounds are obtained from the LHC using the same analysis as for the real parts of these exotic interactions, in such a way that an improvement in future R measurements will be necessary to keep them competitive. The situation is quite different for the D correlation, which sets bounds on a CP violating V, A interaction which are virgin from competition with the LHC using the analysis of [20]. The best constraints so far on such T-violating interaction arise from the neutron decay, with a measurement yielding D n = (−0.9 ± 2.1 · 10 −4 ) [28]. Weaker constraints were obtained for nuclear decays using the decay of 19 Ne, with a combination of measurements yielding a limit D19 N e < 6 · 10 −4 [33]. In terms of the coupling constants from the effective Hamiltonian of beta decay, the present limits on D n translate into a limit on the relative phase between the vector and axial couplings: Im(C V /C A ) = (1.6 ± 6.3) · 10 −4 (90% C.L.) [11]. In practice the neutron and nuclear β decays have a different sensitivity to the relative phase between the vector and axial couplings, which would bring evidence for NP. This sensitivity can be expressed as where X is the neutron or β decaying nucleus, and F (X) is a factor which depends on the nuclear spin and Gamow Teller to Fermi ratio [2]. F (X) is displayed in Tab. 3 for the neutron and different isotopes. The sensitivity of D to NP is found to be superior for some nuclear β decays, by 50% in the case of 23 Mg and 39 Ca, which can be in-trap polarized in the MORA setup. The measurement of the D correlation using the MORA setup presents additional advantages thanks to the confinement of ions: • A 2.5-fold increase in sensitivity to the D correlation by the selection of the most sensitive plane of detection (see δ factor in Eq. 4 in section 4 and the subsequent discussion).
• An arguable suppression of the systematic effects listed in Tab. 4 of ref. [26] for the neutron, thanks to the well controlled characteristics of the ion cloud (small cloud size, homogeneous polarization in the plane of detection).
EDMs in nuclei or in the neutron are other sensitive probes to CP violation. In full generality, these searches are complementary and they can only be compared within specific NP scenarios. In a general EFT framework it is easily shown that they access different combinations of parameters. For some leptoquark models that are not covered by the EFT approach current D measurements of neutron beta decay represent the strongest constraints, even taking into account nEDM bounds [6]. Moreover, the D correlation is theoretically a much cleaner probe to NP than EDM, which relies on hadronic and nuclear computations [5,34]. In summary, the development in the coming years of the in-trap polarization of radioactive ion beams with the MORA setup will permit to attain new frontiers of sensitivity to NP, searching for possible sources of CP violation. With a precision down to a few 10 −5 , the sensitivity to the imaginary phase Im(C V /C A ) will be improved by one order of magnitude compared to the present bounds. The FSI effects mimic a non-zero D correlation at a level which is below the precision of the current experiments. Such effects will eventually be probed for the first time, if no signature of NP is found at a higher level, resulting in a stringent test of these calculations. We will first focus on 23 Mg + which is easier to polarize. This isotope is well produced at JYFL for performing the proof-ofprinciple of the in-trap polarization, and the first D correlation measurement. This first phase at JYFL will permit a study of all the parameters that affect the sensitivity of the MORA setup, before the measurements of the D correlation with highest sensitivity to NP are performed at DESIR, where 23 Mg and 39 Ca will be available with copious intensities from the SPIRAL upgrade and S3-LEB facilities.