Detailed study of the decay of 32Ar
Abstract
In an experiment performed at the SPIRAL1 facility of GANIL, the decay of 32Ar has been studied by means of the ”Silicon Cube” device associated with germanium clover detectors from the EXOGAM array. Beta-delayed protons and rays have been observed and allowed the determination of all relevant decay branches. The Gamow-Teller strength distribution is compared to shell-model calculations and excellent agreement is found. The Fermi strength is inline with expectations. A quasi-complete decay scheme of 32Ar is established.
1 Introduction
Nuclear decay is a powerful tool for the study of nuclear structure for nuclei far from stability. It is a high-precision tool that allows the extraction of effects linked e.g. to the pairing interaction, nuclear deformation, and questions related to the exact formulation of the weak interaction. Beyond its importance for nuclear structure and fundamental interactions, nuclear decay plays also an important role for basically all types of stars at any moment of their evolution. The access to nuclear decay information is facilitated by the fact that -decay experiments are among the most simple-to-implement experiments in nuclear physics. Therefore, beyond the observation of a new nuclide, the first investigations are usually carried out by studying its decay properties blank08review .
Detectors | 1 | 2 | 3 | 4 | 5 | 6 |
---|---|---|---|---|---|---|
DSSSD thickness (m) | 300 | 287 | 270 | 64 | 1000 | 288 |
Large-area detector thickness (m) | 300 | 300 | 500 | 1473 | 150 | 1498 |
The present paper describes a study of the -decay properties of 32Ar. This nucleus is perhaps the most proton-rich isotope of argon bound by the strong interaction, as 31Ar, its neighbour, is predicted by some mass models to be (close to) proton unbound pape88 ; comay88a ; duflo95 ; moeller16 . This nucleus is also a candidate for precision studies of beyond-standard-model contributions to the weak interaction via the search for scaler and tensor currents schardt93 ; adelberger99 ; araujo20 . Its -decay properties have been studied four times in experiments, three times at ISOLDE schardt93 ; hagberg77 ; bjornstad85 and a fourth time in two experiments at NSCL and at ISOLDE bhattacharya08 . While the first experiment of Hagberg et al. hagberg77 only observed protons from the isobaric analogue state (IAS) and verified the isobaric multiplet mass equation (IMME), the second experiment of Björnstad et al. bjornstad85 was mainly concerned with the giant Gamow-Teller (GT) resonance. The third work of Schardt and Riisager schardt93 extracted limits of exotic currents of the weak interaction. The last work executed by Bhattacharya et al. bhattacharya08 dealt with the Fermi strength to determine as precisely as possible the decay strength, the value, for the super-allowed decay of 32Ar in order to deduce the isospin-impurity correction experimentally and compare it to a theoretical prediction. Unfortunately, the authors of this last work were not interested in extracting the distribution of the GT strength, which allows a relatively precise comparison with shell-model calculations. The main drawback for the rather detailed data of Björnstad et al. bjornstad85 is probably the very low -ray detection efficiency, which did not allow the authors to attribute all proton groups observed correctly to decays to the first excited state in 31S (see below). We will compare our results only to the data of Björnstad et al. and Bhattacharya et al., because Schardt et al. gave information only for a few proton lines and Hagberg et al. observed only the decay of the IAS.
The experiment described here was performed at the identification station of SPIRAL1 at GANIL. The data were taken at the same time as data for 31,33Ar already published matea09 ; adimi10 . For the purpose of the present experiment, part of the identification station was dismounted and replaced by the ”Silicon Cube” array matea09 . Although we will also investigate the Fermi strength in the decay of 32Ar, the focus of the paper is on the GT strength distribution, which we will compare to predictions of the nuclear shell model.
2 Experimental procedure and set-up
The 32Ar3+ beam was produced by projectile fragmentation of a 36Ar primary beam at 95 MeV/A accelerated by the CSS cyclotrons of GANIL, which impinged on the SPIRAL1 graphite target. The 32Ar atoms diffusing out of the target and reaching the NANOGAN-III ECR ion source of SPIRAL1 were subsequently ionised and sent to the SPIRAL1 identification station with an energy of 30 keV, where they were intercepted by a 0.9m thick mylar foil. The beam line optics between the target-ion source and the detection set-up was regularly optimised by means of a stable 40Ar beam also produced by the ion source.
The experimental set-up consisted of a cube of six double-sided silicon strip detectors (DSSSDs, 16 X and 16 Y strips with a width of 3 mm) backed by large-area silicon detectors (LASDs, 50 mm x 50 mm). The different thicknesses of the detectors are given in table 1 and their position with respect to the incoming beam is shown in figure 1.

All DSSSD channels were read-out by means of 16-channel pre-amplifier cards mounted directly on the vacuum chamber of the detection set-up and connected to the detectors via a printed circuit board. The LASDs were connected to their pre-amplifiers by means of LEMO cables. The total detection efficiency of this set-up can reach up to 54% matea09 , if all detectors are fully operational. In the present paper, we used for the main analysis only half of the detectors (see in the following paragraph). The Silicon Cube was surrounded by three EXOGAM clover detectors.
Detectors DSSSD 1, 3, 6 worked without major problems (one Y strip missing for DSSSD 3, DSSSD 5, and DSSSD 6 ) and were used for the main data analysis, except for the highest proton energies (see below). Technical problems, which we could not solve prior to the data taking, prevented DSSSD 2 from functioning correctly. It was therefore excluded from the data analysis. The thickness of detector 4 was initially chosen to be thin for the study of the -delayed 2p decay of 31Ar. For the mainly high-energy protons in the present study, it was too thin and also excluded from the data analysis. DSSSD 5 has, due to its thickness, a large contribution to the spectrum from particles. This is not a problem for high-energy protons, but ”pollutes” the spectrum at low proton energies. However, protons with energies above 6000 keV punch through the 300m silicon detectors. Therefore, we used DSSSD 5 for protons above 5700 keV (proton groups 29-31 in table 5), instead of detectors 1, 3, and 6.
The large-area detectors, which were meant to be used in anti-coincidence with the DSSSDs in front of them, were calibrated with particles prior to the data taking. However, we did not use them in the present analysis, because the probability of having the particle and the proton in the same pixel of a DSSSD (see below) is so small that this effect was negligible. We did not use them either to detect high-energy protons punching through the 300m DSSSD, because of their unknown dead-layer and the fact that only an -particle calibration was available.
The energy calibration of the DSSSDs was performed with well-known energies from the literature nndcA=32 . Precise laboratory energies of 14 proton groups from 1210 keV to 6483 keV were calculated from known excited states in either 32Cl or 33Cl and used to calibrate the proton spectra taken for the decays of 32Ar and 33Ar. The precision of the energy of these calibration proton lines ranges from 0.4 keV for the lowest proton energies to 8 keV for the highest energies. As it is not possible to determine an energy uncertainty from the error propagation based on 128 different strips for the four DSSSDs used in the present analysis, we determined the difference between the expected energy and the one determined in the fit of the experimental spectra for four proton groups used in the calibration. We found an RMS deviation of 4 keV, which we have added quadratically to all proton energy uncertainties given in the present work.
Absolute proton detection efficiencies are not needed, because all proton emission branching ratios will be normalised with respect to the most intense proton group from the IAS. The only assumption is that the detection efficiency is the same for all proton energies, which is a good assumption, because the proton detection efficiency should be the solid angle covered by the detectors used with respect to 4.
The energy resolution obtained for the DSSSD sum spectrum from detectors 1, 3, and 6 was of the order of 50 keV (FWHM). To accept an event the signals in the front and back side had to be within 125 keV of each other. Once this condition verified, the energy was taken from the better of the two sides of each DSSSD (Y side for DSSSDs 1 and 6, X side for the others).
The EXOGAM clover detectors were calibrated in efficiency with the standard calibration sources 60Co, 133Ba, 137Cs, and 207Bi. Each of the four elements of the clover detectors was calibrated independently and the spectra were summed after calibration. One segment of clover 1 had a drifting gain and was removed from the data analysis. No add-back was used. The efficiency curve was fit with a straight line on a double-logarithmic scale. We found an efficiency of 3.29(25)% at 1 MeV. This efficiency curve is correct down to energies of 150 - 200 keV, where the entrance window thickness of the germanium detectors and other effects start playing a role and where the straight line on a double-logarithmic scale is no longer valid. The 89.9 keV line from the decay of 32Ar was cut by the ADC threshold in some of the germanium crystals (see below).
The energy calibration was performed crystal by crystal with 11 on-line -ray lines from 511 keV to 4772 keV. The precision of the calibration was checked with the sum spectrum from the 11 segments out of 12 kept for the analysis. A RMS deviation of 0.4 keV was found for the 11 calibration lines and added quadratically to the energy uncertainty for all rays from the decay of 32Ar and of its daughter nuclei. Due to the fact that the 89.9 keV line was partially cut, we refrain from giving an energy value for the present experiment.
For the present analysis, the only data acquisition triggers used were triggers from the DSSSDs. Therefore, the trigger efficiency is different for -delayed proton events (higher efficiency) and for -delayed rays (lower trigger efficiency).
Runs on 31,32,33,34Ar were performed during an on-line data taking of 7 days. A total of 26h 26 min was devoted to the study of 32Ar. The average detection rate of 32Ar was about 100 pps.
3 Experimental results
In the following, we will first discuss the -ray spectra with rays observed in the -decay daughter nucleus 32Cl, in the p daughter 31S, and in the decay of 32Cl with rays from 32S. In a second step, we will turn our attention to the emission of protons with and without -ray coincidences in order to determine proton emission from excited states in 32Cl to the ground and first excited states of 31S. In contrast to Bhattacharya et al., proton emission to higher excited states could not be observed in the present work. We neglect the very weak p (0.026%) and (0.054%) decay channels of 32Cl firestone , which might very slightly alter the total proton emission branching ratio.
For the determination of average values, we use in the calculations one digit more than given in the tables. Therefore, a difference of 1 in the last digit given may occur for calculations with the numbers in the table. If the different values to be averaged are more apart from each other than allowed by errors, we increase the error of the average value by the square root of the normalised as prescribed by the Particle Data Group pdg18 .

3.1 Gamma-ray spectra
The first step in the data analysis is the determination of the -ray energies and their branching ratios. Figure 2 shows the spectrum obtained from the eleven germanium segments used in the analysis. The peaks labelled with letters belong to the decay of 32Ar. Labels , are from levels in 32Cl, in the p daughter 31S (peak ) or from the decay of the daughter nucleus 32Cl (). Peak is not visible in the -ray singles spectrum. It becomes visible only in the coincidence spectrum. In addition, it is partially cut by the ADC thresholds for some of the crystals. Therefore, it is difficult to determine its exact energy and branching ratio. We refrain from doing this, as its branching ratio is the same as the one of the ray at 1078.6 keV (see decay scheme in figure 6). Peak is visible in the singles spectrum, but much better evidenced in a proton--ray coincidence spectrum, as shown in the inset of the upper part of figure 2. The other insets show expanded views of the singles spectrum to better visualise the rays of interest.
Table 2 summarizes the -ray data from this work and from Björnstad et al. bjornstad85 for states in the -decay daughter 32Cl populated by a Gamow-Teller decay and for the ray in the p daughter 31S. Reasonable agreement is obtained for all energies. As the lines from the decay of 32Ar were not included in the energy calibration, we can build the averages from the present work and from Björnstad et al. For the branching ratios, good agreement is also obtained, except for the ray at 1249 keV. The value from our work is overestimated, because this ray is always in coincidence with protons, for which the trigger probability of the data acquisition is significantly higher than in the case of -delayed rays triggered by the particles in the DSSSDs. We refrain therefore from averaging the two values.
The absolute branching ratios can be obtained from the complement of the total proton emission branching ratio from Bhattacharya et al. of BR = 35.58(22)%. We assume here that the ground-state feeding is negligible (see below). The results of this procedure are given in the last row of table 2.
Three rays were observed in the work of Bhattacharya et al. to de-excite the IAS of the ground state of 32Ar in 32Cl. These are peaks , , and , where we admit that peak would have been difficult to identify without knowledge from the work of Bhattacharya et al. Nevertheless, the energies and branching ratios determined show good agreement with the results of Bhattacharya et al. The energies as well as the absolute branching ratios and their averages are presented in table 3 and will be discussed below.
Gamma rays were not only observed in the p daughter nucleus 31S, but also in the decay of the daughter nucleus 32Cl. These are given in table 4. To obtain their absolute branching ratios, we normalised them with the total branching ratio calculated from the total proton branching ratio as given by Bhattacharya et al. bhattacharya08 . In the table, we also present a comparison to data from the literature nndcA=32 . Good agreement is obtained indicative of the fact that the extrapolated -ray detection efficiencies work well up to energies as high as 5 MeV.
peak | -ray energies (keV) | relative branching ratios (%) | average absolute | ||||
---|---|---|---|---|---|---|---|
this work | Björnstad bjornstad85 | average | this work | Björnstad bjornstad85 | average | branching ratio (%) | |
a | - | 89.9(1) | 89.9(1) | 38.3(32) | 36.6(18) | 37.0(16) | 13.8(7) |
b | 460.7(4) | 461.1(1) | 461.1(1) | 100.0(69) | 100.0(49) | 100.0(40) | 37.2(18) |
c | 707.1(4) | 707.4(2) | 707.3(2) | 100.0(69) | 100.0(49) | 100.0(40) | 37.2(18) |
d | 1078.1(4) | 1078.7(2) | 1078.6(2) | 38.3(32) | 36.6(18) | 37.0(16) | 13.8(7) |
e | 1168.3(5) | 1168.5(2) | 1168.5(2) | 31.8(38) | 36.9(18) | 35.9(20) | 13.4(8) |
f | 1248.5(5) | 1248.4(3) | 1248.4(3)11footnotemark: 1 | 17.1(17)2221 | 5.8(6) | 5.8(6) | 2.2(2) |
1 In nndcA=31 , this energy is given as 1248.87(9) keV. We use in the following text 1249 keV
22footnotetext: 2 2 p triggered ray, not used in averaging procedure due to detection efficiency difference, see text
peak | -ray energies (keV) | absolute branching ratios (%) | ||||
---|---|---|---|---|---|---|
this work | Bhattacharya bhattacharya08 | average | this work | Bhattacharya bhattacharya08 | average | |
j | 2838.7(34) | 2836.0(10) | 2836.2(10) | 0.50(37) | 0.24(3) | 0.24( 3) |
l | 3877.7(42) | 3877.5( 3) | 3877.5( 3) | 1.03(22) | 1.58(8) | 1.52(18) |
p | 5047.5(50) | 5046.3( 4) | 5046.3( 4) | 0.22(14) | 0.10(2) | 0.10( 2) |
peak | -ray energies (keV) | absolute branching ratios (%) | ||
---|---|---|---|---|
this work | nndcA=32 | average | ||
g | 1548.0(20) | 2.4(4) | 2.7(5) | 2.5(3) |
h | 2230.5( 2) | 81.5(44) | 70.0(30) | 73.6(54) |
i | 2463.8(10) | 3.0(7) | 3.0(3) | 3.0(3) |
k | 3317.5(15) | 1.3(4) | 1.9(3) | 1.7(3) |
m | 4281.5(15) | 1.3(5) | 2.0(1) | 2.0(1) |
n | 4433.0(20) | 0.9(6) | 0.6(2) | 0.6(2) |
o | 4770.0(15) | 14.3(17) | 15.5(15) | 15.0(11) |
3.2 Proton spectra
Figure 3 shows the singles proton spectrum from DSSSDs 1, 3 and 6. A spectrum was also created with the same proton spectrum, however, in coincidence with a ray at 1249 keV, the ray de-exciting the first excited state of 31S to its ground state. The result is also shown in figure 3. A total of 30 proton lines were observed in the present work. A very weak peak at 912 keV observed by Bhattacharya et al. bhattacharya08 was not seen in the present work nor in the work of Björnstad et al. bjornstad85 .

All proton lines were fit by a Gaussian and a low-energy tail to take into account the possibility that the full charge was not always collected in the silicon detectors or that interference effects create a low-energy tail. One single Gaussian plus a tail was sufficient for all proton groups to get a good description of the experimental peak, except for the line at 3353 keV, which is due to the emission of a proton from the IAS, the strongest peak in the spectrum of figure 3. In this case (peak 14), a second Gaussian was added. The integral of the Gaussians yielded the number of counts in the different peaks and allowed the determination of relative branching ratios. In table 5, we present the energy values of these proton lines and their relative branching ratios normalised to 1000 for the strongest proton line from the IAS. The data from the present work are compared to the data available from the work of Björnstad et al. bjornstad85 and Bhattacharya et al. bhattacharya08 , this latter work giving only proton lines with energies below about 4 MeV if they feed the first excited state in 31S (except for peak 1), and all proton lines above 4 MeV.
proton laboratory energies | relative branching ratios | |||||||||
---|---|---|---|---|---|---|---|---|---|---|
peak | this | Bhattacharya | Björnstad | average | this | Bhattacharya | Björnstad | average | CM | absolute |
work | et al. bhattacharya08 | et al. bjornstad85 | work | et al. bhattacharya08 | et al. bjornstad85 | energies | BR | |||
1 | 606.8( 60) | 610.0(100) | 607.0(100) | 607.5( 46) | 32.3(165) | 18.8( 4) | 17.0(31) | 18.76( 39) | 627.3( 46) | 0.385( 8) |
2∗ | - | 912.0( 50) | - | 912.0( 50) | - | 0.70(40) | - | 0.70( 40) | 941.6( 50) | 0.014( 8) |
3∗ | 1205.5( 41) | 1218.0( 50) | 1214.0(100) | 1210.8( 42) | 14.68(78) | 19.0(22) | 17.9(18) | 15.57(112) | 1250.1( 42) | 0.319(23) |
4∗ | 1676.6(123) | - | - | 1676.6(123) | 1.43(65) | - | - | 1.43( 65) | 1731.1(123) | 0.029(13) |
5 | 2121.3( 40) | - | 2124.0(100) | 2121.7( 37) | 176.6(30) | - | 174.8(204)3332 | 176.6( 30) | 2190.6( 37) | 3.62(7) |
6∗ | 2146.0( 52) | 2145.0( 50) | - | 2145.5( 36) | 18.7(40 | 12.80(40) | 12.9(18)3332 | 12.86( 41) | 2215.2( 36) | 0.264( 9) |
7∗ | 2395.9( 41) | 2394.0( 50) | - | 2395.1( 32) | 9.7(47) | 5.6(11) | - | 5.8(11) | 2472.9( 32) | 0.119(22) |
8 | 2422.9( 40) | - | 2423.0(100) | 2422.9( 37) | 352.8(47) | - | 358.1(393) | 352.8(47) | 2501.6( 37) | 7.24(11) |
9∗ | 2510.0( 40) | 2515.0( 50) | 2511.0(100)g | 2511.9( 30) | 36.2(12) | 29.3(11) | 39.3(87) | 32.6(25) | 2593.5( 30) | 0.668(51) |
10∗ | 2615.7( 86) | - | - | 2615.7( 86) | 2.47(63) | - | - | 2.47( 63) | 2700.6( 86) | 0.051(13) |
11∗ | 2778.5( 65) | 2870.0( 50)44footnotemark: 4 | 2768.0(100)g | 2775.4(54) | 3.83(50) | 30.0(100)44footnotemark: 4 | 5.68( 87) | 4.29( 80) | 2865.6(54) | 0.088(16) |
12 | 3117.4( 47) | - | 3113.0(100) | 3116.6( 43) | 1.32(13) | - | 1.67( 39)3332 | 1.36( 12) | 3217.8( 43) | 0.028( 2) |
13∗ | 3117.0(105) | - | - | 3117.0(105) | 0.76(41) | - | 0.95( 27)3332 | 0.90( 22) | 3218.3(105) | 0.018( 5) |
14 | 3352.7( 40) | 3353.0( 50) | 3353.5( 30) | 3353.2( 22) | 1000.0(41) | 1000.0(59) | 1000.0(699) | 1000.0( 34) | 3462.1( 22) | 20.51(17) |
15∗ | 3584.4(239) | 3581.0( 50) | 3592.0(100) | 3583.2( 44) | 3.5(19) | 2.40(40) | 3.48( 68)3332 | 2.70( 34) | 3699.6( 44) | 0.055( 7) |
16 | 3605.2( 82) | - | - | 3605.2( 82) | 3.56(45) | - | 4.82( 69)3332 | 3.94( 58) | 3722.3( 82) | 0.081(12) |
17∗ | 3651.8( 68) | 3649.0( 50) | 3643.0(100) | 3649.0( 37) | 3.21(37) | 3.20(30) | 3.93( 87) | 3.25( 23) | 3767.5( 37) | 0.067( 5) |
18 | 3725.9( 48) | - | 3732.0(100) | 3727.1( 44) | 4.03(17) | - | 6.29( 83)3332 | 4.13( 46) | 3848.1( 44) | 0.085( 9) |
19∗ | 3778.2( 63) | 3785.0( 50) | - | 3782.4( 39) | 2.01(33) | 5.20(50) | 3.32( 80)3332 | 3.00( 98) | 3905.2( 39) | 0.062(20) |
20 | 3997.8( 50) | 3984.0( 50) | 3994.0(100) | 3991.2( 46) | 9.04(21) | 11.0(10) | 11.79( 87) | 9.26( 50) | 4120.8( 46) | 0.190(10) |
21∗ | 4347.5( 96) | 4340.0( 50) | 4341.0(100)g | 4341.5( 41) | 1.28(43) | 2.3(16)3332 | 1.44( 53)3332 | 1.38( 33) | 4482.5( 41) | 0.028( 7) |
22 | 4344.9( 41) | 4340.0( 50) | 4341.0(100)g | 4342.7( 30) | 4.94(18) | 4.8(14)3332 | 5.55( 91)3332 | 4.96( 18) | 4483.8( 30) | 0.102( 4) |
23 | 4524.6( 41) | 4529.0( 50) | 4521.0(100)e | 4525.9( 30) | 4.45(13) | 5.40(50) | 5.24( 44) | 4.57( 21) | 4672.8( 30) | 0.094( 4) |
24∗ | 4622.0( 46) | 4630.0( 50) | 4621.0(100) | 4625.2( 32) | 1.5( 7) | 1.60(50) | 1.75( 44) | 1.65( 30) | 4775.4( 32) | 0.034( 6) |
25∗ | 4862.5(180) | 4869.0( 50) | 4858.0(100) | 4866.5( 43) | 2.0(14) | 2.60(30) | 2.62( 44) | 2.59( 24) | 5024.6( 43) | 0.053( 5) |
26 | 4966.7( 79) | 4997.0(100) | 4975.0(100) | 4977.5( 89) | 0.38(17) | 1.00(20) | 0.61( 17) | 0.63( 17) | 5139.1( 89) | 0.013( 4) |
27 | 5556.1( 41) | 5567.0( 50) | 5552.0(100) | 5559.8( 40) | 5.39(13) | 7.60(80) | 5.68( 44) | 5.47( 25) | 5740.3( 40) | 0.112( 5) |
28 | 5673.4(124) | 5699.0(100) | 5675.0(100) | 5683.7( 85) | 0.20( 5) | 1.80(80) | 0.48( 31) | 0.22( 8) | 5868.3( 85) | 0.005( 2) |
29 | 5815.3( 44) | 5833.0( 50) | 5817.0(100) | 5822.5( 60) | 4.00(21) | 5.40(50) | 4.37( 44) | 4.23( 32) | 6011.6( 60) | 0.087( 7) |
30 | 6066.2( 52) | 6097.0(100) | 6060.0(100) | 6070.5( 88) | 1.15(11) | 1.10(20) | 0.92( 17) | 1.09( 8) | 6267.7( 88) | 0.022( 2) |
31 | 6358.6( 73) | 6396.0(100) | 6347.0(100) | 6365.2(133) | 0.56( 8) | 0.60(20) | 0.61( 17) | 0.57( 7) | 6572.0(133) | 0.012( 1) |
3 only one peak observed, sharing between ground and first excited states recalculated according to branching ratios from other experiments, see text
33footnotetext: 24 not used in the averaging procedure, see text
In the following, we quickly discuss all proton lines from table 5 individually. For this discussion, we will use the average proton energies. The numbers correspond to the peak numbers used in figure 3 and in table 5.
-
•
1: In the present work, we can see the 608 keV peak only under the condition that we require a signal in two different strips (one being from the proton, one from the associated particle) and keep only the higher energy. This condition drastically reduces the background at low energies and makes this peak visible. For a correct normalisation for this peak, the 3353 keV line was also analysed under the same condition. It is a proton to the ground state of 31S.
-
•
2: The 912 keV line is not identified in our data. Only Bhattacharya et al. have a tiny proton branch at this energy to the first excited state.
-
•
3: The 1211 keV line is seen in the three experiments as a proton emission to the first excited state in 31S with comparable intensities.
-
•
4: The authors from previous work did not identify a proton line at 1677 keV, but we have clear evidence for such a line in the singles and in the coincidence spectra. Björnstad et al. have also a small peak at this energy, however, they do not mention it. We attribute it to the decay to the first excited state of 31S.
-
•
5-6: Björnstad et al. identify a single peak at 2124 keV. However, in our data we have clear evidence for a ground-state decay at 2122 keV and a decay to the first excited state at 2146 keV. This is in agreement with the work of Bhattacharya et al. who observe also a peak to the first excited state at 2146 keV. The branching ratio to the ground state and to the first excited state as determined in the present work add up to the branching ratio given by Björnstad et al. So we determine from the average of our work and the one of Bhattacharya et al. the relative branching ratios to the ground and first excited states and use these value to distribute the branching ratio found by Björnstad et al. to the ground and first excited states for these authors.
-
•
7: Björnstad et al. have no peak at 2395 keV. In the present work and in Bhattacharya et al., this peak is visible as a decay to the first excited state, although with almost a factor of 2 difference in branching ratio.
-
•
8: The peak at 2423 keV is observed in the present work and by Björnstad et al. as a decay to the ground state in 31S.
-
•
9: The 2512 keV line is observed in the three experiments with comparable intensity. As Bhattacharya et al., we see it as a decay to the first excited state. Björnstad et al. identify it as a decay to the ground state. Due to the clear evidence from the present work and from Bhattacharya et al., we attribute it to the decay to the first excited state.
-
•
10: No 2616 keV peak is observed in the work of Bhattacharya et al. and of Björnstad et al., but there is a visible shoulder in the present data and in those of Björnstad et al., which we identify as a proton emission to the first excited state. The observation of this peak may be questionned due to the fact that there is no corresponding decay to the ground state of 31S. This needs to be investigated in a new experiment.
-
•
11-13: At 2775 and 3117 keV, Björnstad et al. have two peaks as in the present work, but these authors do not give them in coincidence with the ray at 1249 keV, which could be due to the small branching ratio and their small -ray detection efficiency. We clearly see this coincidence and attribute these two peaks to decays to the first excited state. However, in our data, the 3117 keV peak is also in the ground-state spectrum with a branching ratio of 1.32% (peak 12). Bhattacharya et al. see in this region a broad structure at 2870 keV with a much higher relative intensity (30(10)%) as compared to 5.9(7)% for the present work and 8.3(10)% for Björnstad et al. The data of Bhattacharya et al. for the 2870 keV peak were not used for the averaging of the proton energies nor for the branching ratios. Like for peaks 5 and 6, we distribute the branching ratio of Björnstad et al. for the 3117 keV line to the ground-state and excited-state decay.
-
•
14: The 3353 keV peak is the most prominent peak for all three data sets. It is therefore taken for the normalisation. As will be discussed below, a small contribution in coincidence with the ray de-exciting the first excited state in 31S cannot be completely excluded.
-
•
15-16: Around 3600 keV, we see a somewhat broader peak with one part decaying to the ground state and one to the first excited state of 31S. We attribute the 3583 keV line to the decay to the first excited state and the 3605 keV line to the decay to the ground state. Bhattacharya et al. see the decay to the first excited state, whereas Björnstad et al. see one peak ”predominantly” decaying to the first excited state, indicative that in their data there might be also decay strength to the ground state. We therefore distribute their decay branching ratio to both decays according to the relative branching ratios from our data and from Bhattacharya et al.
-
•
17: Nice agreement is obtained for all three data sets for the peak at 3649 keV as a decay to the first excited state.
-
•
18-19: The peak at 3732 keV in Björnstad et al. has at least two components. We find a ground-state peak at 3727 keV and an excited-state peak at 3782 keV, in agreement with Bhattacharya et al. If we sum the two peaks observed here for decays to the ground and first excited states, we come close to the integral of Björnstad et al.. We distribute their decay branching ratio to both decays according to the relative branching ratios from our data and from Bhattacharya et al.
-
•
20: The 3991 keV peak is observed in all three experiments with comparable intensities as a decay to the ground state.
-
•
21-22: For the peak at about 4340 keV, we clearly see it for the ground state decay (22) and in the coincidence spectrum (21). So we share its activity between the ground state and the first excited state also for Bhattacharya et al. and Björnstad et al.
-
•
23: The 4526 keV peak is seen in the three experiments, but Björnstad et al. claim to observe it predominantly in coincidence with the 1249 keV ray, whereas Bhattacharya et al. see it mainly as a ground-state decay. We can only observe a ground-state branch at this energy. We attribute the strong ground-state decay observed by Bhattacharya et al. (0.54(5)% with respect to the IAS decay) to this peak 23. We do the same for the peak of Björnstad et al. (0.12(1)% absolute branching ratio). We neglect the small possible excited-state contribution (0.03% with respect to the IAS peak 14) from Bhattacharya et al.
-
•
24-25: A proton group at 4625 keV is seen in the three experiments with comparable intensities as a decay to the first excited state. The 4867 keV peak decays also to the first excited state.
-
•
26-31: The proton groups at 4978, 5560, 5684, 5823, 6071, and 6365 keV are identified as decays to the ground state in 31S in all three experiments.
The last three peaks were evaluated with DSSSD 5. We found that the branching ratio determined with detectors 1, 3, and 6 was the same as with detector 5 within the error bars for proton group 29, but factors of 2 and 6 lower for the high-energy proton groups 30 and 31, respectively.
We note that the proton peak energies of Bhattacharya et al. are systematically higher than those given for our data and those from Björnstad et al. above about 4500 keV. Schardt and Riisager schardt93 mention that the energies of Björnstad et al. should be corrected, because these authors used an -particle energy calibration at high energies. However, the correction proposed moves the proton energies of Björnstad et al. even above the values of Bhattacharya et al. and far away from our proton energies. For the present paper, we keep the energy values from Björnstad et al.
The doublets that we observe at around 3130 keV, 3595 keV, and 4340 keV have been given particular attention. We first fit a single proton peak in the -coincidence spectrum and in the ground-state decay spectrum, this latter being obtained by subtracting the spectrum for the decay to the first excited state with the correct normalisation from the proton singles spectrum. This procedure was verified by fitting a doublet to the proton singles spectrum by keeping one of the two energies and the width of the peaks fixed. It allowed us also to determine the relative intensities of the ground-state and excited-state decays and to observe that none of the two alone gives the intensity of the singles-spectrum peak, when fitted with only one proton line.
The relative branching ratios were averaged (with the restrictions mentioned above) as given in table 5. With the technique used in the present work, we cannot determine the number of implanted 32Ar in any precise way to transform the relative branching ratios into absolute proton emission branching ratios. In order to determine absolute branching ratios, we use again the total proton emission branching ratio of Bhattacharya et al. However, if we just add the integrals of the identified peaks to determine the total number of protons detected in the present work, the procedure is not correct. Some activity is not identified in peaks (e.g. activity between peaks 26 and 27). We therefore determined the total number of protons by fitting an exponential to the low-energy -particle contribution in the upper spectrum of figure 3 and subtracted this contribution before determining the total number of counts in the spectrum, which corresponds to the total number of protons detected. If we compare this number to the sum of the integrals of all identified proton groups, we find that a small correction factor of 1.033(14) is needed. With this correction and the total proton emission branching ratio from Bhattacharya et al. of BR = 35.58(22)%, we can determine the absolute branching ratios given in the table.
From the number of protons in the singles spectrum after the subtraction of the -particle contribution at low energies and from the number of protons in the spectrum conditioned by a ray at 1249 keV after subtraction of random coincidences (see below), we can determine the proton branching ratio to the first excited state of 31S to be 2.5(3)%, which can be compared to the value obtained by Bhattacharya et al. of 2.3(4)%. Björnstad et al. found a value of 1.9(2)%, whereas we determine 2.1(2)% with a slightly different absolute normalisation for the latter data as already discussed by Bhattacharya et al. (see table 2).
3.3 Proton- coincidence data
Up to now, we have treated the proton and the -ray data almost independently. Only to better evidence certain proton groups most likely being due to a decay to the first excited state in 31S, we used proton--ray coincidences to produce the lower-statistics spectrum shown in figure 3. However, the fact that a peak shows up in this -ray coincidence spectrum does not necessarily mean that the protons are really in coincidence with the ray at 1249 keV, although, as we will show, this is the case most of the time. We note also that in our experiment, due to the total number of decays observed and due to our efficiencies, the detection limit for proton- coincidences with one count corresponds to a branching ratio of 0.0015%. However, for higher-intensity proton peaks the accidental coincidence probability is much higher.

In order to evidence proton--ray coincidences, we overlay in figure 4 the proton singles spectrum from DSSSD 1, 3, and 6 as well as from DSSSD 5 with a normalisation factor of 3 with the proton spectrum generated from the same detectors requiring a coincidence with a ray between 1240 keV and 1254 keV (we do not subtract for the moment random coincidence, see below). However, as the -ray detection efficiency decreases the number of counts observed in this coincidence spectrum, we multiplied the number of counts in each channel by a factor of 35.1(19) determined as the average of the inverse of the -ray detection efficiency at 1249 keV and the ratios of the number of counts of the 1211 keV and the 2512 keV proton peaks in singles and with the -ray coincidence. As just multiplying this spectrum with a factor of 35.1 would yield a spectrum difficult to ”read”, we ”randomised” the multiplied spectrum by folding each event with a Gaussian with a sigma of 2 keV to smoothen the spectrum. The result is shown in figure 4 together with the proton singles spectra and the original -ray coincidence proton spectrum.
In these spectra, one clearly sees that for some proton groups one recovers almost the full proton singles intensity with the scaled coincidence spectrum (e.g. protons at 1211 keV, 2512 keV, 4867 keV). However, there is also significant proton activity under the peak at 3353 keV, which is supposed to be the decay of the IAS in 32Cl to the ground state of 31S. So this peak should not be in coincidence with rays at 1249 keV.

In order to determine the ”real” coincidence counts and remove random coincidences, we show in figure 5 the same -ray coincidence spectrum with the 1249 keV ray and overlay a spectrum generated with a coincidence left and right of the 1249 keV peak. This latter spectrum should consist basically only of random coincidences. In addition, we determine the ratios between the number of counts in both spectra for 21 proton-peak regions, which include 25 of our 31 proton groups. For the analysis of the ratios in figure 5, we centered the regions, given by the vertical lines, at the energies of the proton groups (mainly for decays to the first excited state) with a width of 40 keV in the low-energy region and 50 keV for higher proton energies (the last 4 regions). Peaks 9-10, 12-13, 15-16, and 21-22 were evaluated within a single region. For some peaks, the intervals are smaller than 40 keV because two neighbouring peaks are too close in energy.
This evidences that we deal with almost the same number of counts e.g. for the region of the 3353 keV peak in both spectra. For regions with real coincidences (e.g. the peaks at 1211 keV, 2146 keV, 2775 keV, etc.), this ratio is small. However, often ground-state decays and excited-state proton groups overlay or are very close. A good ratio limit seems to be a value of 0.2 (see below), where real coincidence peaks have a smaller and random coincidence peaks have a larger value. From these ratios, we determine that peaks with numbers 3, 4, 6, 7, 9, 10, 11, 13, 15, 17, 19, 21, 24, and 25 are proton emissions to the first excited state of 31S. For peak 2, the situation is more difficult to analyse due to the presence of a strong background. According to Bhattacharya et al., it is a rather weak proton emission to the first excited state, which is not observed neither in our work nor by Björnstad et al. We keep it as a proton emission to the first excited state of 31S.
excitation | branching | B(GT)/ | c.m. proton | c.m. proton |
---|---|---|---|---|
energy (keV) | ratio (%) | B(F) | energy g.s. (keV) | energy exc. state (keV) |
1168.5( 2) | 62.57( 29) | 0.4704( 99) | - | - |
2208.4( 46) | 0.145( 31) | 0.0020( 4) | 627.3 [ 1] | - |
3771.7( 30) | 3.635( 68) | 0.1387( 38) | 2190.6 [ 5] | 941.6 [ 2] |
4081.6( 28) | 7.555(112) | 0.3644( 92) | 2501.6 [ 8] | 1250.1 [ 3] |
4561.1(123) | 0.029( 13) | 0.0021( 9) | - | 1731.1 [ 4] |
4799.0( 43) | 0.028( 2) | 0.0024( 2) | 3217.8 [12] | - |
5046.1( 3) | 22.63(25) | 3.991(93) | 3462.1 [14] | 2215.2 [ 6] |
5303.0( 30) | 0.200( 25) | 0.0276( 35) | 3722.3 [16] | 2472.9 [ 7] |
5425.4( 27) | 0.753( 52) | 0.1169( 84) | 3848.1 [18] | 2593.5 [ 9] |
5530.6( 86) | 0.051( 13) | 0.0087( 22) | - | 2700.6 [10] |
5699.3( 35) | 0.278( 19) | 0.0569( 41) | 4120.8 [20] | 2865.6 [11] |
6063.6( 45) | 0.120( 6) | 0.0365( 19) | 4483.8 [22] | 3218.3 [13] |
6254.0( 31) | 0.094( 4) | 0.0355( 18) | 4672.8 [23] | - |
6529.6( 44) | 0.055( 7) | 0.0293( 38) | - | 3699.6 [15] |
6597.5( 38) | 0.067( 5) | 0.0385( 28) | - | 3767.5 [17] |
6732.8( 56) | 0.074( 20) | 0.0514(141) | 5139.1 [26] | 3905.2 [19] |
7317.0( 45) | 0.140( 8) | 0.2269(144) | 5740.3 [27] | 4482.5 [21] |
7449.4( 85) | 0.005( 2) | 0.0090( 34) | 5868.3 [28] | - |
7602.9( 73) | 0.121( 9) | 0.3147(244) | 6011.6 [29] | 4778.1 [24] |
7853.5( 39) | 0.075( 5) | 0.3123(229) | 6267.7 [30] | 5024.6 [25] |
8153.1(133) | 0.012( 1) | 0.0903(114) | 6572.0 [31] | - |
Another finding that might deserve particular attention in a future high-resolution study is the fact that, for the 3353 keV peak number 14, the random coincidences are exclusively located in the high-energy part of this peak. If one subtracts the random coincidences, a peak in the low-energy tail of this strong peak remains, which could be indicative of a small proton peak with a real coincidence with the 1249 keV ray. This might be important because it reduces slightly the Fermi strength discussed below.
The ratio of 0.2 used above might seem somewhat arbitrary, but it is meant to be an indication only for the ratio of counts in the different regions between the two spectra. One can also subtract the spectrum with the random coincidence (condition left and right of the 1249 keV peak) from the coincidence spectrum (condition on the 1249 keV peak). Unfortunately, the statistics becomes rather limited. However, if we fit this difference spectrum, which should contain only real coincidences, with energies fixed to the values of table 5, the numbers of counts for peaks without real coincidences (i.e. peaks 5 and 8) are compatible with zero and confirm thus the above conclusions. As stated in the previous paragraph, counts remain in the low-energy part of peak 14.
4 Results and discussion
With the data analysed and presented in the previous sections, we can determine the probabilities of -decay feeding of different states in 32Cl and their decay to the ground or first excited state of 31S. From this information, we can build the decay scheme of 32Ar. We will then establish the Gamow-Teller strength distribution and compare it to shell-model predictions. Finally, we will analyse the Fermi strength determined from the proton and decay of the IAS.
In the following, we will assume that the decay branch from 32Ar to the ground state of 32Cl is negligibly small. This is justified, if mirror symmetry holds even only roughly, as the value for the mirror decay is as large as 8.21 nndcA=32 corresponding to a branching ratio for 32Ar decay to the ground state of 32Cl of 510-5. The different USD-type Hamiltonians (see below) give branching ratios ranging from 0.007% (USD brown88 ) to 5.6% (USDA brown06 ).
4.1 Branching ratios to states in 32Cl and decay scheme of 32Ar
In the preceding sections, we have been able to attribute all proton groups to decays from states in 32Cl to two states in 31S, the ground state and the first excited state at 1249 keV. No decay to other states of 31S could be evidenced. Proton energies having a center of mass difference in agreement with a -ray energy of 1249 keV (we required a difference of less than 1.5 times the sum of the uncertainties of the proton energies) stem from the same level in 32Cl, but decay to the ground and first excited states of 31S. Table 6 summarizes this information and gives the excitation energies of states populated in the decay of 32Ar (calculated with a proton separation energy of = 1581.1(5) keV ame2016 ), the branching ratio to these levels and the center-of-mass energies of protons involved in the decay of these levels to the ground or first excited state in 31S.
If we compare the levels deduced in the present work with the adopted levels from the evaluation of states in 32Cl nndcA=32 , we propose one additional level at 5530.6 keV, whereas for four states from nndcA=32 we cannot identify the corresponding proton peaks (excitation energies 2611 keV, 4167 keV, 4439 keV, 5794 keV). The first level comes from a ”private communication”. This level is observed in none of the -decay experiments. However, it would correspond to a level clearly present in the mirror nucleus 32P at 2740 keV. The second and third states come from the work of Björnstad et al. from two proton branches attributed to a ground-state decay. We clearly identify the two proton branches as decays to the first excited state of 31S. The last level comes from the very broad 2870 keV proton line identified by Bhattacharya et al. only and discarded in our analysis (see above). So from the four levels, we will keep only the 2611 keV line.
In addition to the unbound levels, which decay by proton emission, we give also the feeding of the bound 1+ level at 1168.5 keV of 32Cl, which -decays to the ground state of 32Cl. This decay takes place via a -ray cascade that can be determined with the rays observed and attributed to the decay of levels of 32Cl (see table 2).
All this information allows us finally to establish a rather complete decay scheme for 32Ar given in figure 6. Information from the present paper is combined with the results from Björnstad et al. and Bhattacharya et al. and data from the literature e.g. for the half-life and the value nndcA=32 . The branching ratios, the -decay Q values and the half-life of 32Ar allowed us also to determine the values for the different decay branches as given in the decay scheme.
4.2 Gamow-Teller strength distribution
The information gathered to establish the decay scheme can also be used to determine the Gamow-Teller strength distribution B(GT). However, this procedure is only correct, if all decay strength can be included. As stated above, about 3.3% of the total decay strength from unbound levels in 32Cl could not be attributed to a proton peak. This induces certainly only a small error in the B(GT) distribution, in particular at high excitation energy. We will determine this strength distribution in the standard way by using the identified peaks and their branching ratios as well as in a more precise way by using directly the proton spectrum.

For this latter purpose, we use the proton singles spectrum (from detectors 1, 3, and 6 up to 5.70 MeV, from detector 5 above), subtract the -coincident spectrum from the same detectors (multiplied with the factor of 35.1) to obtain the spectrum of protons to the ground state of 31S, take the -coincident spectrum and shift it by 1249 keV, and add it back to the ground-state spectrum with the correct multiplication factor. This new spectrum can then be converted into a spectrum of population of excited states in 32Cl. We next exclude the region of the IAS. The total number of counts in the spectrum obtained together with the total proton branching ratio allow us to give to each event and thus to each channel of this spectrum a branching ratio. The excitation energy, to which each channel corresponds, enables us to determine the value for each channel and then the Gamow-Teller strength B(GT) as
with = 6143.6 s (see e.g. hardy20 ) and the axial-vector coupling constant = -1.2756 pdg20 . is the statistical rate function towner15 depending on the -decay value and is the partial half-life defined as
with BR being the branching ratio attributed to each spectrum channel. The B(GT) value for the feeding of the bound state at an excitation energy of 1169 keV is calculated in a similar way. We used T1/2 = 98(2) ms nndcA=32 .
The result of this procedure is shown in figure 7 together with the B(GT) distribution determined with the identified proton peaks. One can see that starting from about 6 MeV excitation energy, the B(GT) distribution obtained from the identified proton peaks stays below the B(GT) distribution from the full proton spectrum.

4.3 Comparison with shell-model calculations
We compare our results to calculations obtained in the sd-shell model space with protons and neutrons in the set of orbitals. We will use the so-called “universal” sd (USD) type Hamiltonians. Universal means that the same set of single-particle energies and two-body matrix elements can be applied to all states whose wavefunctions are thought to be dominated by sd-shell model configurations. The two-body matrix elements are allowed to have a smooth mass dependence. This is in contrast with the ab-initio in-medium similarity renormalization group (IMSRG) stroberg17 ; stroberg19 method where the Hamiltonian is nucleus dependent. At present, the rms deviation between the experimental and theoretical excitation energies for IMSRG is about 500 keV. Making this universal assumption, one can fine-tune the two-body matrix elements to minimize the deviation between experimental and theoretical energies and reduce the rms deviation down to about 150 keV.
We will use five USD-type Hamiltonians: (i) the original isospin-conserving USD Hamiltonian based on data for sd-shell nuclei with available up to 1983 brown88 ; wildenthal84 , (ii) the improved isospin-conserving USDA and USDB Hamiltonians based on data for nuclei with available up to about 2006 brown06 , and (iii) the most recent isospin-nonconserving USDC and USDI Hamiltonians based on data for all sd-shell nuclei magilligan20 . The isospin-nonconserving part contains the Coulomb interaction and an isotensor strong interaction.
4.3.1 Comparison of experiment and theory for 32P
We start with understanding to what extent the USD-type Hamiltonians are able to describe the experimental spectra. The energy levels for the mirror nucleus 32P are better known, so we use this nucleus for a first comparison of experiment and the USD-type hamiltonians. The experimental excitation energies for 32P and those obtained with the USDC Hamiltonian are shown in figure 8 (left part).


The results for other USD-type Hamiltonians are similar. The experimental energies are from nndcA=32 . The lowest known negative parity state, a 2- state at 3.3 MeV, is shown in blue. Other states with an assigned negative parity above that state are not shown. These negative-parity states would be associated with 1 excitations, being dominated by one nucleon moved from the sd shell to the shell, and are not within our model space.
There is a good match between experiment and theory up to 2.3 MeV. Two states with unassigned spin-parity at 2.3 and 2.6 MeV are indicated by the black dots. These are only seen in one experiment, 31P(n,), and need to be confirmed. Above that there is good agreement between experiment and theory for positive-parity states up to about 5 MeV. The rms deviation between the theoretical and experimental excitation energies up to 5 MeV is about 150 keV. This is true for all the USD-type Hamiltonians and is consistent with the rms deviations observed with these Hamiltonians in other nuclei brown06 ; magilligan20 .
If there are pairs of experiment-theory levels with excitation energies that agree within 200 keV we consider this ”good agreement”. There is a well established low-lying 0+ T=1 state at 512.7 keV in agreement with theory, and the 0+ T=2 at 5.072 MeV is also in agreement with theory. Below 6 MeV there is only one more 0+ T=1 state at 4.81 MeV in comparison with a tentative (0+) state experimentally at 4.62 MeV. The location of these 0+ T=1 states is important for the interpretation of isospin mixing with the 0+ T=2 state signoracci11 .
Below the 0+ T=2 state, there are seven 1+ states in the calculation compared to eight observed experimentally. Starting around 5 MeV the experimental level density becomes larger than theory. For example, there are two 2+ states near 5 MeV that do not have a theoretical counterpart. Thus, starting at about 5 MeV there are levels in experiment that are intruder states relative to the sd-shell model space. These would be associated with 2 excitations, being dominated by two nucleons moved from the sd shell to the shell. This means that a direct comparison of experimental and USD model-space states above about 5 MeV is no longer straight forward.
The USD Hamiltonians provide a good starting point for the wavefunctions. But the true Hamiltonian differs from the USD Hamiltonian by an rms average of 150 keV in the diagonal many-body matrix elements. We might assume that there are off-diagonal matrix elements of a similar size. The reason comes from the assumption made about the USD-type Hamiltonians, i.e. that they apply only to the sd model-space degrees of freedom, and that they are the same (except for some smooth mass scaling) for all states ascribed to the sd model space.
The implication of this is that matrix elements associated with observables such as decay, decay, and spectroscopic factors will have uncertainties due to level mixing. For example, consider the decay of 32Ar to two 1+ states in 32Cl that have a USD energy separation of = 1 MeV with USD B(GT) values of 0.5 and 0.1. The experimental states and will be mixtures of the USD states and with and . If mixing is obtained from a two-dimensional matrix with an off-diagonal matrix element of keV, the resulting error on the B(GT) is approximately 0.07 for both, with the sum of 0.60 being exactly the same. Thus, the error on these quantities is not simply proportional to their values, but is the same for both, 0.07 in this example. The implication is that the theoretical error associated with observables is not proportional to the size of the observable, but is rather closer to a constant value. This is consistent with the comparison between experiment and theory made for a large number of observables in richter08 .
4.3.2 Comparison of experiment and theory for 32Cl
In the decay of 32Ar, seven 1+ states plus the 2611 keV state discussed above are observed below the 0+ T=2 state compared to the 8 expected from levels in 32P (see figure 8, right part). Up to about 3 MeV excitation energy, a one-to-one match of states in 32P and 32Cl exists in experiment and theory. Above that energy too many states observed in 32Cl have still an unknown spin/parity and renders a further comparison difficult. It is reasonable to assume that the first intruder state is the negative parity state at around 3 MeV, as in the case of 32P.
The B(GT) values for our five Hamiltonians are compared to experiment in figure 9 up to 6 MeV on a linear scale. We use a quenching factor of about 0.6 obtained from the ratios of the theoretical over the experimental half-life for the different interactions. All of the USD Hamiltonians give similar results and are in qualitative agreement with experiment, in particular for the strong transitions at 1.168 MeV and 4.082 MeV followed by a gradual rise of the B(GT) sum up to 6 MeV. Above 5 MeV the calculated 1+ states must mix with the intruder 1+ states (see the discussion for 32P). This mixing spreads the strength, but does not change the summed value.

The spread of the five calculations provides a measure of the theoretical uncertainty with the model assumptions of a universal sd-shell Hamiltonian. The experimental results generally lie within this spread. For the weak state near 2.2 MeV, the experimental B(GT) of 0.0020(4) is smaller than most of the calculations with the exception of USDA (0.0010). The weak state predicted near 2.7-2.8 MeV is not observed in experiment. There is a suggested (1+) at 2.611 MeV nndcA=32 , which has a 1+ mirror level in 32P. The experimental B(GT) to a state at this energy has an upper limit of B(GT) 0.0001. Most of the calculations are much larger than this value with the smallest value coming from USD with B(GT) = 0.008. To obtain the small experimental value one has to invoke state mixing beyond the model assumptions as discussed above.
USDC | USDI | exp. | ||||
E (MeV) | (g.s.) / (ex. state) | E (MeV) | (g.s.) / (ex. state) | E (MeV) | BR(g.s.) / BR(ex. state) | |
3.8320 | 13.33 | 3.8640 | 9.41 | 3.7717 | 258(147) | |
4.0990 | 11.73 | 4.1640 | 9.32 | 4.0816 | 22(1) | |
4.8480 | 51.85 | 4.8260 | 58.33 | 4.7990 | 56(56) | |
5.0670 | 47.86 | 5.0560 | 55.00 | 5.0463 | 77(2) |
In figure 7, we show the full B(GT) distribution up to the value limit. Evidently, the match between experiment and theoretical predictions is excellent up the highest states populated in the decay of 32Ar. The total B(GT) strength summed over all final states for all these Hamiltonians is about 7.0. Thus about 60% of the Ikeda Gamow-Teller strength sum rule is observed in the 32Ar decay window.
We have also calculated the proton-decay and -decay widths for the 1+ states with the methods used for reaction rates in richter11 . For the proton-decay width we calculate , where C2S are the spectroscopic factors obtained from the sd-shell wavefunctions, and are the single-particle proton-decay widths obtained from proton scattering from a Woods-Saxon potential with the resonance energy constrained to the experimental value. The states that are unbound to proton decay to the ground state and first excited state in 31S have proton-decay widths that are much larger than the -decay widths. This is consistent with the present observations. The proton-decay branching ratio for the decay to the ground state divided by the branching ratio to the first excited state are compared between experiment and theory in table 7. The largest disagreement is observed for the state near 3.8 MeV that is relatively strong in the 32Ar decay. However, the calculated =0 spectroscopic factors are very small. For USDC they are 0.007 to the 31S 1/2+ ground state and 0.0053 to the 31S 3/2+ excited state. Such small spectroscopic factors are very sensitive to the state mixing beyond the model assumptions.
4.4 The Fermi strength to the IAS
The feeding of the IAS in 32Cl can be determined in part from its decay by proton emission to the ground (proton peak 14) and the first excited (peak 6) states. Although the IAS is largely proton unbound and one would expect that it decays only by proton emission, Bhattacharya et al. have shown that it decays also by -ray emission to three states in the -decay daughter nucleus 32Cl. We can clearly identify the first two rays (labelled and in table 3), whereas we consider our evidence for the highest energy ray only tentative (labelled in the same table). This decay of the unbound IAS is favoured due to the isospin-forbidden character of proton emission from the IAS. Only due to a most likely quite small T=1 impurity of the IAS, this proton emission can happen. The proton-emission branching ratio is determined in our work to be 20.77(17)%, whereas the -ray branches sum up to 1.86(18)% yielding a total branching ratio to the IAS of 22.63(25)%. All these values are in perfect agreement with those already given by Bhattacharya et al. bhattacharya08 . In fact, this is not very astonishing, as they are averages from the present data with those from Bhattacharya et al. and from Björnstad et al. and include in addition the total proton-emission branching ratio from Bhattacharya et al.
This super-allowed branching ratio allows us to determine the value of the super-allowed decay branch to be 1514(35), where we used nndcA=32 a half-life of 98(2) ms, a value of 11134.3(18) keV for the ground-state decay and thus a value to the IAS of 6088.0(18) keV yielding f = 3494.6(58). An electron-capture probability of = 0.067% firestone was also included. To determine the corrected value (see e.g. hardy20 ), we use in addition the theoretical isospin-breaking correction =2.0(4)% and the radiative correction = 1.145(41)% from Bhattacharya et al. and obtain = 1501(35) s. As this is a value for a T=2 transition, it can be compared to the value determined for the well-known 0 0+, T=1 decays compiled regularly by Hardy and Towner hardy20 divided by a factor of 2. The value from the latest compilation is 3072.24(57) s, which shows that the presently determined value is in agreement with expectations.
The isospin non-conservation that is contained in the USDC and USDI Hamiltonians results in isospin mixing between the 0+, T=2 state and 0+, T=1 states. The strongest mixing is with a state 226 (308) keV lower in energy with the USDC (USDI) hamiltonian. This results in a splitting of the total Fermi strength of B(F)=4 into B(F)=3.87 (B(F)=3.92) for the T=2 state and 0.13 (0.07) for the T=1 state below the IAS for the USDC (USDI) interaction. The remaining Fermi strength is distributed over the higher-lying 0+ states in the region of 6 to 10 MeV not observed in the present experiment.
The calculated proton-decay width of this 0+, T=1 state is 36 keV (C2S=0.042) with USDC and 37 keV (C2S= 0.042). A possible candidate for the 0+, T=1 state below the IAS could be the experimentally observed state at 4779 keV. However, the branching ratio to this state is about a factor of 25 lower than predicted. Another candidate could be the state at 4561 keV. However, its branching ratio is about a factor of 8 lower than predicted.
From the branching ratios for the IAS, we determine a ratio of the decay widths of R = / = 11.2(11). The literature gives a total width of the IAS of = + = 20(5) eV nndcA=32 . These values allow us to determine = /(1+R) = 1.8(5). This compares well with the calculated value of = 1.3 eV from both USDC and USDI. The experimental proton width is thus = 18.2(5) eV and has to be compared to the calculated values of 525 eV for USDC and 600 eV for USDI. Since the proton width comes from isospin mixing, the conclusion is that the isospin mixing from the experiment is much smaller than in the calculations, which will have consequences for the interpretation of the isobaric multiplet mass equation for A=32 signoracci11 . This conclusion is confirmed by the determination of the experimental B(F) value calculated from the value of the super-allowed decay branch to be 3.991(93) in agreement with the expected value of 4, but also with the theoretical values of 3.87/3.92.
5 Conclusion
We have presented a complete study of the decay of 32Ar performed at the identification station of SPIRAL1 with the silicon cube detection array and three EXOGAM -ray detectors. The data obtained for -delayed protons and rays were compared to data from Björnstad et al. and Bhattacharya et al. and good agreement was obtained. This allowed us to average the available experimental data, determine feedings of excited states of 32Cl and extract values and the Gamow-Teller strength distribution. This distribution was compared to shell-model calculations using the USD, USDA, USDB, USDC, and USDI interactions and excellent agreement was obtained. Finally, the Fermi strength was determined to be in agreement with the value expected for a super-allowed 0 0+ decay. Although there might be two candidate states for a 0+, T=1 state, which could mix with the 0+, T=2 isobaric analog state, we believe that these states are rather Gamow-Teller fed states.
Acknowledgment
We express our gratitude to the EXOGAM collaboration for providing us with the germanium detectors. We thank the whole GANIL and, in particular, the accelerator staff for their support during the experiment. This work was partly funded by the Conseil régional de la Nouvelle Aquitaine and the EU through the Human Capital and Mobility program. We acknowledge support from CICYT via contract FPA2007-62179. BAB is supported by NSF via grant PHY-1811855.
References
- (1) B. Blank and M. J. G. Borge, Prog. Nucl. Part. Phys. 60, 403 (2008).
- (2) A. Pape and M. S. Antony, At. Data Nucl. Data Tab. 39, 201 (1988).
- (3) E. Comay, I. Kelson, and A. Zidon, At. Data Nucl. Data Tab. 39, 235 (1988).
- (4) J. Duflo and A. Zuker, Phys. Rev. C 52, R23 (1995).
- (5) P. Moeller, A. Sierk, T. Ichikawa, and H. Sagawa, At. Data Nucl. Data Tab. 109, 1 (2016).
- (6) D. Schardt and K. Riisager, Z. Phys. A 345, 265 (1993).
- (7) E. G. Adelberger et al., Phys. Rev. Lett. 83, 1299 (1999).
- (8) V. Araujo-Escalona et al., Phys. Rev. C 101, 055501 (2020).
- (9) E. Hagberg et al., Phys. Rev. Lett. 39, 792 (1977).
- (10) T. Björnstad et al., Nucl. Phys. A443, 283 (1985).
- (11) M. Bhattacharya et al., Phys. Rev. C 77, 065503 (2008).
- (12) I. Matea et al., Nucl. Instr. Meth. A 607, 576 (2009).
- (13) N. Adimi et al., Phys. Rev. C 81, 024311 (2010).
- (14) C. Ouellet and B. Singh, Nucl. Data Sheets 112, 2019 (2011).
- (15) R. Firestone, Table of Isotopes, 8th Edition, John Wiley & Sons, Inc., New York, ISBN (1996).
- (16) M. Tanabashi et al., Phys. Rev. D 98, 1 (2018).
- (17) C. Ouellet and B. Singh, Nucl. Data Sheets 114, 209 (2013).
- (18) B. A. Brown and B. H. Wildenthal, Ann. Rev. Nucl. Part. Sci. 38, 29 (1988).
- (19) B. A. Brown and W. A. Richter, Phys. Rev. C 74, 034315 (2006).
- (20) M. Wang et al., Chinese Physics C 41, 030003 (2017).
- (21) J. C. Hardy and I. S. Towner, Phys. Rev. C 102, 045501 (2020).
- (22) Particle Data Group, P. A. Zyla et al., Prog. Theor. Exp. Phys. 2020, 1 (2020).
- (23) I. S. Towner and J. C. Hardy, Phys. Rev. C 91, 015501 (2015).
- (24) A. Magilligan and B. A. Brown, Phys. Rev. C 101, 064312 (2020).
- (25) S. R. Stroberg et al., Phys. Rev. Lett. 118, 032502 (2017).
- (26) S. R. Stroberg, H. Hergert, S. Bogner, and J. D. Holt, Ann. Rev. Nucl. Part. Sci 69, 307 (2019).
- (27) B. H. Wildenthal, Prog. Part. Nucl. Phys. 11, 5 (1984).
- (28) A. Signoracci and B. A. Brown, Phys. Rev. C 84, 031301(R) (2011).
- (29) W. A. Richter, S. Mkhize, and B. A. Brown, Phys. Rev. C 78, 064302 (2008).
- (30) W. A. Richter, B. A. Brown, A. Signoracci, and M. Wiescher, Phys. Rev. C 83, 065803 (2011).