Abstract
Across the physics disciplines, the 186Pb nucleus is the only known system, where the two first excited states, together with the ground state, form a triplet of zero-spin states assigned with prolate, oblate and spherical shapes. Here we report on a precision measurement where the properties of collective transitions in 186Pb were determined in a simultaneous in-beam γ-ray and electron spectroscopy experiment employing the recoil-decay tagging technique. The feeding of the \({0}_{2}^{+}\) state and the interband \({2}_{2}^{+}\to {2}_{1}^{+}\) transition have been observed. We also present direct measurement of the energies of the electric monopole transitions from the excited 0+ states to the 0+ ground state. In contrast to the earlier understanding, the obtained reduced transition probability \(B(E2;{2}_{1}^{+}\to {0}_{2}^{+})\) value of 190(80) W.u., the transitional quadrupole moment \(| {Q}_{t}({2}_{1}^{+}\to {0}_{2}^{+})| =7.7\)(33) eb and intensity balance arguments provide evidence to reassign the \({0}_{2}^{+}\) and \({0}_{3}^{+}\) states with predominantly prolate and oblate shape, respectively. Our work demonstrates a step-up in experimental sensitivity and paves the way for systematic studies of electric monopole transitions in this region. These electric monopole transitions probe the nuclear volume in a unique manner and provide unexploited input for development of the next-generation energy density functional models.
Similar content being viewed by others
Introduction
The governing interactions in atomic nuclei give rise to diverse quantum phenomena such as single-particle motion, nucleon pairing and collectivity, and can result in a configuration that prefers a specific shape. At one extreme, competing shape-driving configurations can appear within a small energy range, a phenomenon commonly known as shape coexistence.
One of the richest regions of shape coexistence is formed by very neutron-deficient nuclei with proton number Z close to the magic number 82 and neutron number N close to 104 (mid-shell)1,2. Shape coexistence in this region was first established in the early 1970s in laser spectroscopy experiments that discovered a sudden change in the nuclear charge distribution of neutron-deficient Hg isotopes3. Over the years, a whole arsenal of spectroscopic techniques has been employed to study the competing structures in this region. For example, rotational bands have been investigated via in-beam γ-ray spectroscopy4, β-decay and α-decay fine structure measurements have probed the level energies of the band-head 0+ states5,6,7, whereas lifetime and Coulomb excitation experiments have shed light on the collectivity of transitions connecting the low-spin states8,9,10,11. The lifetimes of excited 0+ states have been measured for 190,192,194Pb nuclei12.
The discovery of two low-lying excited 0+ states in the N = 104 neutron mid-shell nucleus 186Pb in an α-decay fine-structure experiment established a feature in the atomic nucleus that remains unique across different physics domains5. While in its ground state the 186Pb nucleus is spherical13, its first and second excited states were associated with predominantly π(2p − 2h) and π(4p − 4h) excitations across the Z = 82 shell gap and assigned with deformed oblate and prolate shape, respectively5. This triplet of 0+ states, which lie within an excitation energy range of 700 keV, represents a system where the perturbation of the ground state by less than half a per mille of its binding energy results in configuration changes that are associated with different shapes. However, despite all efforts, the linking transitions from rotational bands to the deformed band-head 0+ states have remained unobserved. Results from previous in-beam spectroscopy experiments on 186Pb have been reported in refs. 14,15,16,17,18. In contrast to ref. 5, many theoretical works18,19,20,21,22,23,24,25,26,27,28,29 have proposed opposite shape assignments for the 0+ states. The band-head 0+ states are the lowest excited states, which due to angular momentum conservation, can only de-excite via electric monopole (E0) transitions. Therefore, in-beam observation of conversion electrons is required.
In the present work, results from an in-beam spectroscopic study of 186Pb employing the sage spectrometer30 are presented. sage is a state-of-the-art instrument that builds on years of endeavour put into research and development of γ-ray31 and solenoidal electron32,33 spectrometers. It allowed for the needed gain in sensitivity to probe the transitions between the low-spin states in 186Pb and provided the first rigorous experimental evidence that supports the assignment of the \({0}_{2}^{+}\) state with predominantly prolate shape. The present work responds to an increasing demand, also recalled in the recent review article34, to advance our understanding on electric monopole transitions.
Results
Recoil-gated, α-tagged (RDT from now on) background subtracted prompt γ-ray and conversion-electron energy spectra are shown in Figs. 1 and 2, respectively. The most prominent γ-ray peaks have been labelled according to their transition energies. The region covering the non-yrast band, which is associated with an oblate shape, is expanded in the inset of Fig. 1. In the electron-energy spectrum in Fig. 2, the most prominent conversion-electron lines with corresponding K-, L- and M-components have been marked and labelled accordingly. The low-energy region below 150 keV is dominated by the δ-electron background.
Deconvolution of the electron energy spectrum has been performed and overlaid in Fig. 2. The blue line illustrates the total fit made for the electron energy spectrum. The dotted lines correspond to the calculated K-, L, and M-electron components of the most prominent E2 and M1 transitions, also including contaminant events. After subtracting the calculated components from the total fit, there remains an excess of electron counts, indicated by the red line in Fig. 2. This excess will be discussed below.
The high-statistics RDT γ − γ and γ − e− coincidence data allowed new transitions in 186Pb to be established. The properties of the new transitions and transitions relevant to the present work are listed in Table 1. Based on coincidence relations and intensity balances, the new transitions have been placed in the level scheme as shown in Fig. 3. It is interesting to note that, despite the same target and beam combination as in refs. 17,18, the slightly higher beam energy of the present work resulted in much higher relative feeding of the non-yrast states in 186Pb.
De-excitation path passing through the \({0}_{2}^{+}\) state
An electron peak at 447 keV stands out in the RDT electron-energy spectrum, see Fig. 2. The peak cannot be associated with any previously known γ-ray transition in 186Pb, suggesting the transition has a strong E0 component. If this is considered to be a K-conversion electron component, the corresponding transition energy would be 535 keV and the L-conversion electron component would have an energy of 520 keV. Indeed, an electron peak at that energy is observed in the RDT electron energy spectrum in Fig. 2. These are in a perfect agreement with the \({0}_{2}^{+}\to {0}_{1}^{+}\) transition energy extracted in the α-decay of 190Po5, confirming the earlier discovery of the \({0}_{2}^{+}\) state and the assignment of electron peaks in Fig. 2 with a transition energy of 535 keV.
The de-excitation path passing through the \({0}_{2}^{+}\) state at 535 keV is investigated in Fig. 4. An RDT γ-ray energy spectrum in coincidence with the 447 keV prompt electrons is presented in Fig. 4a. In addition to the known yrast-band transitions, a peak at 125.0(24) keV is observed. The energy matches with the level energy difference of 127(2) keV between the \({2}_{1}^{+}\) and \({0}_{2}^{+}\) states, thus the transition has been placed in the level scheme accordingly. The groups of counts appearing at around 160 and 350 keV could not be assigned with any transitions in 186Pb. The total intensity value of \({I}_{{{{{{{{\rm{total}}}}}}}}}({2}_{1}^{+}\to {0}_{2}^{+})=31(10)\) was extracted from the weighted average of the \({0}_{2}^{+}\to {0}_{1}^{+}\) K-electron events normalised with the K-electron components of the \({2}_{1}^{+}\to {0}_{1}^{+}\), \({6}_{1}^{+}\to {4}_{1}^{+}\) and \({8}_{1}^{+}\to {6}_{1}^{+}\) transitions observed in coincidence with the \({4}_{1}^{+}\to {2}_{1}^{+}\;\gamma\)-ray transition. The corresponding spectrum is shown in Fig. 4b, where the presence of the K-conversion electrons of the 535 keV transition and peaks used for normalisation are evident. The spectrum also includes events at ~113 keV that could correspond to the L-conversion electrons originating from the \({2}_{1}^{+}\to {0}_{2}^{+}\) transition as expected from the high L-conversion coefficient of αL ≈ 1.4. However, they overlap with the δ-electron background preventing precise extraction of the energy and intensity. The γ-ray intensity value of \({I}_{\gamma }({2}_{1}^{+}\to {0}_{2}^{+})=9(3)\) was obtained using the value αtotal(E2) = 2.35(4)35. Such a small intensity is beyond the observational limit in the RDT γ-ray singles events.
Observation of the E0 transition de-exciting the \({0}_{3}^{+}\) state
A group of events is found in the RDT electron-energy spectrum around 574 keV (Fig. 2) where the K-conversion electron component of the \(E2({2}_{1}^{+}\to {0}_{1}^{+})\) transition is expected. Based on the observed γ-ray intensity and the calculated K-conversion coefficient35, the number of corresponding K-conversion events in the peak has been determined. As a result, the intensity of the peak at ~574 keV cannot be solely explained by the \({2}_{1}^{+}\to {0}_{1}^{+}\) transition or by other observed γ-ray transitions. It is also noteworthy, that the full width at half maximum resolution of the peak is 17 keV, which is larger than the expected resolution of ~14 keV. Thus, it is considered the peak comprises two components, one that is assigned with the \({2}_{1}^{+}\to {0}_{1}^{+}\) transition and the other to the excess events. The composition of the peak with the corresponding fits is visualised in Fig. 5. The centroid of the excess events is at 571 keV. If this was an L-conversion electron component of a transition, the corresponding more intense K-component should be observed at 498 keV, which is clearly not the case. The non-observation of a corresponding γ-ray transition suggests a strong E0 component. The suggested energy of 659(4) keV is in excellent agreement with the \({0}_{3}^{+}\to {0}_{1}^{+}\) transition [650(20) keV] first observed in the α-decay of 190Po5. Consequently, the 571 keV electron peak is assigned as the K-conversion electron component of the \({0}_{3}^{+}\to {0}_{1}^{+}\) transition with \({I}_{{{{{{{{\rm{K}}}}}}}}}({0}_{3}^{+}\to {0}_{1}^{+})\) = 11(7). Unfortunately, transitions feeding the \({0}_{3}^{+}\) state are beyond the observational limit.
Discovery of the interband \({2}_{2}^{+}\to {2}_{1}^{+}\) transition
The interband \({2}_{2}^{+}\to {2}_{1}^{+}\) transition has not been observed previously. The missing γ-ray intensity analysis based on γ-γ coincidences, similar to that performed for 188Pb which revealed considerable E0 components in the interband J → J transitions36, did not indicate any strong E0 components in 186Pb17,18. However, this could be explained by insufficient γ-γ statistics as the intensity balances listed in Table 1 of ref. 18 do not rule out an additional transition feeding the \({2}_{1}^{+}\) state.
The present RDT electron data show evidence for the K-conversion electron component of a \({2}_{2}^{+}\to {2}_{1}^{+}\) 283 keV transition as demonstrated in the RDT electron energy spectrum in Fig. 2. The electron peak overlaps partly with the K-conversion electron component of the contaminant (11−) → (10−) 276 keV transition in 186Tl. Based on the recoil-gated γ-γ and γ-e− coincidence data, a value of \({\alpha }_{{{{{{{{\rm{K}}}}}}}},\exp }(276)=0.33(7)\) was obtained. This is in-line with the calculated value of αK,M1(276) = 0.404(6) for an M1 transition35. Accordingly, the K-conversion electron component of the contaminant 276 keV M1 transition is plotted (green dashed line) in Fig. 2. The total fit (blue line) shows the sum of the excess electrons in 186Pb (red line) and the M1 transition in 186Tl and matches well with the observed electrons.
Further evidence for the highly-converted \({2}_{2}^{+}\to {2}_{1}^{+}\) transition is presented in Fig. 6, where RDT γ-ray (a) and electron (b) energy spectra in coincidence with 662 keV γ rays stemming for the \({2}_{1}^{+}\to {0}_{1}^{+}\) transition are shown. The γ-ray energy spectrum in Fig. 6a presents peaks associated with yrast-band transitions, whereas the γ-ray transition at 283 keV (see the inset) is absent. The upper limit for the γ-ray intensity of the \({2}_{2}^{+}\to {2}_{1}^{+}\) transition was estimated to be Iγ < 20. It is noteworthy, that apart from the \({8}_{2}^{+}\to {6}_{2}^{+}\) transition at 424 keV with Iγ = 61(18), other non-yrast transitions are below the observational limit in Fig. 6a. In Fig. 6b, in addition to the two preceding yrast-band transitions, an electron peak emerges at 196(2) keV corresponding to the energy of the K-conversion electron component of the 283 keV \({2}_{2}^{+}\to {2}_{1}^{+}\) transition. The intensity of the K-conversion electron component \({I}_{{{{{{{{\rm{K}}}}}}}}}({2}_{2}^{+}\to {2}_{1}^{+})\) = 10(5) given in Table 1 is the weighted average of intensities extracted from RDT singles electrons (Fig. 2) and from RDT electrons with a gate on the \({2}_{1}^{+}\to {0}_{1}^{+}\;\gamma\)-ray transition (Fig. 6). Based on obtained intensities, the lower limit of αK > 0.5 for the conversion coefficient can be extracted.
Electron-conversion component of the \({4}_{2}^{+}\to {4}_{1}^{+}\) transition
The γ-e− coincidence statistics were not sufficient to obtain information about the conversion-electron component of the interband \({4}_{2}^{+}\to {4}_{1}^{+}\) transition by gating on the feeding transition. Since the \({4}_{2}^{+}\to {4}_{1}^{+}\) transition energy overlaps with the \({8}_{1}^{+}\to {6}_{1}^{+}\) transition energy, gating with a transition below the \({4}_{1}^{+}\) state included both transitions and resulted in too few statistics for deconvolution analysis. However, RDT γ-ray and electron singles data were sufficient for the deconvolution of the electron energy spectrum based on the extracted γ-ray intensities. This is illustrated in Fig. 7, where a peak at the corresponding K-conversion electron energy has been divided into its components. The properties of these five components, as identified in Fig. 7, are based on the observed E2 γ-ray transitions. A comparison of the sum of the calculated intensities of the electron components to the measured electron energy spectrum reveals a clear excess of electrons. In order to match with the experimental observations, an additional electron excess component (red line in Fig. 7) with free fit parameters has been included, leading to a measured intensity of I = 12(8). If the electron excess was of L-electron origin, the corresponding K-electron line at 252 keV should be observed, which is not the case. An αK = 0.32(30) conversion coefficient was extracted for the \({4}_{2}^{+}\to {4}_{1}^{+}\) transition.
Discussion
The αK > 0.5 value obtained for the \({2}_{2}^{+}\to {2}_{1}^{+}\) transition is higher than expected for an M1, E2 or mixed E2/M1 transition (αK,M1 = 0.409, αK,E2 = 0.075), suggesting the presence of an E0 component. This is only possible for a ΔK = 0 transition between states having the same spin and parity. Consequently, the level at 945 keV can be firmly assigned as the \({J}_{i}^{\pi }={2}_{2}^{+}\) state. Using similar arguments for the \({4}_{2}^{+}\to {4}_{1}^{+}\) transition, the tentative assignment of the \(({4}_{2}^{+})\) state at 1337 keV in refs. 17,18 can be now fixed. Based on these assignments and the angular distribution information in refs. 17,18, levels at 1738 and 2162 keV can be also firmly assigned as \({J}_{i}^{\pi }={6}_{2}^{+}\) and \({J}_{i}^{\pi }={8}_{2}^{+}\) states, respectively. Accordingly, spins and parities of the non-yrast band up to the \({8}_{2}^{+}\) state are fixed. It also follows from the ΔK = 0 criteria that the M1 components for the interband transitions are ruled out.
The E0 transition probability is defined as
where ρ2 is the transition monopole strength and Ω(E0) is the electronic factor37,38. For the monopole strength, the following relation can be extracted:
where IK(E0) is the K-conversion electron intensity of the E0 transition, Iγ(E2) and Wγ(E2) are the γ-ray intensity and transition rate of the competing E2 transition, respectively, and ΩK(E0) is the electronic factor for the K-electron conversion of the E0 transition.
Although the information on the transition rates is lacking, it is intriguing to use estimates to evaluate the monopole strengths of the interband transitions. A \(B(E2;{2}_{1}^{+}\to {0}_{1}^{+})=6\) W.u. has been extracted from the lifetime of the \({2}_{1}^{+}\) state8,9, while \(B(E2;{2}_{2}^{+}\to {0}_{1}^{+})=3.9\) W.u. has been measured for 184Hg in a Coulomb excitation experiment10. In both cases, transitions are between states of different intrinsic configurations. Assuming \(B(E2;{2}_{2}^{+}\to {0}_{1}^{+})=5\) W.u., a ρ2 monopole strength value of 100(60) × 10−3 can be extracted for the E0 component of the \({2}_{2}^{+}\to {2}_{1}^{+}\) transition. Similarly, ρ2 = 40(30) × 10−3 can be obtained from the present data for the \({4}_{2}^{+}\to {4}_{1}^{+}\) transition assuming \(B(E2;{4}_{2}^{+}\to {2}_{2}^{+})=100\) W.u., a typical value for an oblate band8,9.
A large monopole strength is a fingerprint of mixing between two coexisting states with different deformation39. Since the spherical Jπ states are expected to lie well above the energy of the prolate and oblate Jπ states, they can be omitted in the following mixing calculations. According to the two-level mixing model, the monopole strength can be expressed as:
where Z is the atomic number, a and b are the mixing amplitudes with the relation a2 + b2 = 1 and β2,i is the quadrupole deformation parameter of the \({J}_{i}^{+}\) state. Quadrupole deformation values of ∣β2∣ = 0.29(5) and ∣β2∣ = 0.17(3) have been extracted for the prolate and oblate bands, respectively, from lifetime measurements performed in this region8,9. Adopting these values, the amount of mixing needed to produce the monopole strength estimated for the \({2}_{2}^{+}\to {2}_{1}^{+}\) transition is of the order of 10% (a2 = 0.9), while it is only ~4% (a2 = 0.96) for the \({4}_{2}^{+}\to {4}_{1}^{+}\) transition. It is noteworthy, that the presence of a possible M1 component in the interband transitions would reduce the value obtained for the monopole strength, further decreasing the amount of mixing. Accordingly, small prolate-oblate mixing can be expected for the yrast \({2}_{1}^{+}\) and \({4}_{1}^{+}\) states, which is in line with the results obtained in the lifetime measurements8,9. Calculations performed in the IBM framework predict more pronounced mixing40.
The intra-band transition intensities of the non-yrast band (oblate) were extracted to be around 10% of the corresponding yrast-band (prolate) transitions and with decreasing spin the majority of the de-excitation of the non-yrast states ends up to the yrast band. The feeding of the excited 0+ states is weak as the feeding transitions have to compete against higher energy transitions to the 0+ ground state. In the present work, the first and second excited 0+ states were fed at the level of ~2.2% and ~1.1% of the total intensity feeding the ground state, respectively.
The only observed feeding of the \({0}_{2}^{+}\) state is via the \({2}_{1}^{+}\to {0}_{2}^{+}\) transition. Based on the branching ratios and the lifetime of the \({2}_{1}^{+}\) state [18(5) ps]8,9, a reduced transition probability \(B(E2;{2}_{1}^{+}\to {0}_{2}^{+})=190(80)\) W.u. can be obtained. For comparison, a weighted average of B(E2) values measured for higher-spin prolate band transitions in 186,188Pb nuclei is 390(40) W.u. The transition quadrupole moments ∣Qt∣ extracted from these B(E2) values under the assumption of a rotating quadrupole-deformed nucleus are 7.7(33) eb and 9.0(5) eb, respectively, suggesting that the \({0}_{2}^{+}\) state has a predominant prolate component.
The de-excitation of the \({0}_{2}^{+}\) state with \({I}_{{{{{{{{\rm{total}}}}}}}}}({0}_{2}^{+}\to {0}_{1}^{+})=26(6)\) extracted from the corresponding K-electron peak in the singles RDT electron energy spectrum of Fig. 2 is fully covered by the feeding transition with the intensity of \({I}_{{{{{{{{\rm{total}}}}}}}}}({2}_{1}^{+}\to {0}_{2}^{+})=31(10)\) extracted from the electron energy spectrum of Fig. 4 gated with 261 keV γ rays. Nevertheless, it is interesting to estimate the potential side-feeding within the given margins. As discussed earlier, the \({2}_{2}^{+}\) state can be considered almost pure oblate state. Provided the \({0}_{2}^{+}\) state were predominantly oblate, the \(B(E2;{2}_{2}^{+}\to {0}_{2}^{+})\) value could be estimated to be ~50 W.u. Assuming \(B(E2;{2}_{2}^{+}\to {0}_{1}^{+})=5\) W.u., an intensity value of \({I}_{{{{{{{{\rm{total}}}}}}}}}({2}_{2}^{+}\to {0}_{2}^{+})=19(2)\) can be obtained. Such a large side-feeding does not fit within the given margins, suggesting the \({0}_{2}^{+}\) state does not have a notable oblate admixture.
The \({0}_{3}^{+}\)/\({0}_{2}^{+}\) mixing can be further assessed using monopole strengths obtained for the transitions from the \({0}_{3}^{+}\) state. Applying Eq. (3) and typical prolate and oblate quadrupole deformation values of ∣β2∣ = 0.29(5) and ∣β2∣ = 0.17(3)8,9, respectively, the value of \({\rho }^{2}({0}_{3}^{+}\to {0}_{2}^{+})=290\times 1{0}^{-3}\) is obtained for maximal prolate-oblate mixing (i.e. a2 = b2 = 0.5). Based on the lifetimes of the predominantly oblate \({0}_{2}^{+}\) states in 190,192,194Pb nuclei12, a ρ2 = 10 × 10−3 value is estimated for the \({0}_{3}^{+}\to {0}_{1}^{+}\) transition. Consequently, by using \({I}_{{{{{{{{\rm{total}}}}}}}}}({0}_{3}^{+}\to {0}_{1}^{+})=13(8)\) obtained in the present work, \({I}_{{{{{{{{\rm{total}}}}}}}}}({0}_{3}^{+}\to {0}_{2}^{+})=170(100)\) can be extracted. In the case of smaller prolate-oblate mixing (e.g. a2 = 0.05), \({I}_{{{{{{{{\rm{total}}}}}}}}}({0}_{3}^{+}\to {0}_{2}^{+})=34(20)\) is obtained. While there is not much room for side-feeding for the \({0}_{2}^{+}\) state, the small amount of mixing between the excited 0+ states is better in line with the extracted intensity balances.
The feeding of the \({0}_{3}^{+}\) state was beyond the observational limit. Taking \(B(E2;{2}_{2}^{+}\to {0}_{1}^{+})=5\) W.u. as above and assuming the feeding of the \({0}_{3}^{+}\) state is solely from a predominantly oblate \({2}_{2}^{+}\) state, a reduced transition probability of \(B(E2;{2}_{2}^{+}\to {0}_{3}^{+})=200(120)\) W.u. can be obtained. This shows that the \({2}_{2}^{+}\to {0}_{3}^{+}\) transition can be a collective transition within a band.
To conclude, the B(E2) value shows that the \({2}_{1}^{+}\to {0}_{2}^{+}\) transition is a collective transition, whereas the transition quadrupole moment, albeit with large uncertainty, suggests a transition within the prolate band rather than a transition between different configurations. The intensity balances between the low-spin states imply small mixing between the excited 0+ states. In particular, the oblate admixture in the \({0}_{2}^{+}\) state has been found to be small. Consequently, the \({0}_{2}^{+}\) state is assigned as the band-head of the yrast band and associated with predominantly prolate shape. Accordingly, the \({0}_{3}^{+}\) state is considered to be the band-head of the non-yrast band associated with predominantly oblate shape. In Fig. 8, level energy systematics in the neutron-deficient even-mass Pb nuclei has been plotted including the new assignments of the 0+ states. The parabolic behaviour of the intruder states is now more evident at the minima that take place at N = 104 and N = 106 for prolate and oblate configurations, respectively. This is also well in-line with the very recent study on the odd-mass 187Pb nucleus with N = 105 which shows that the prolate minimum is below the oblate one in energy41.
Many theoretical calculations predict the prolate minimum in 186Pb to be lower in energy than the oblate one (see refs. 18,19,20,21,22,23,24,25,26,27,28,29), although opposite results have also been obtained5,42,43. Moreover, the lowest oblate minimum in energy is predicted at N = 106 by several calculations20,22,24,26,27,28,29 and experimental findings support the oblate assignment of the \({0}_{2}^{+}\) state in 188Pb44. Based on the IBM calculations for 186Pb, the \({0}_{2}^{+}\) [\({0}_{3}^{+}\)] state is a mixture of configurations assigned with spherical, prolate and oblate shapes with admixtures of 3% [3%], 63% [38%] and 34% [58%] weights in their wave function, respectively40. Experimental results obtained on mixing of the 0+ states are only available for the oblate and spherical states in the heavier 190,192,194Pb isotopes45.
However, the fact that the assignments of the excited 0+ states presented here are in contrast to what was proposed by Andreyev et al. warrants further discussion. Their arguments were based on (i) the level energy systematics of the oblate 0+ states, (ii) the extrapolation from the high-spin members of the prolate band and (iii) reduced 190Po α-decay widths5. The competition between the intruding oblate and prolate minima in the neutron-deficient Pb isotopes is most striking at the N = 104 mid-shell. In the light of shape staggering between different isotopes46 and configuration mixing45 observed in the region, it is questionable whether the level energy systematics can be used as a valid argument for association of 0+ states with specific shape, particularly when direct measurements of shapes, such as the spectroscopic quadrupole moments, are lacking. Secondly, extrapolating the band-head energy from higher spin members of the yrast band is not straightforward in the presence of configuration mixing, particularly when all states involved in mixing are not known. Finally, there is no firm experimental evidence regarding the shape of the ground state in 190Po and the existing theoretical predictions are conflicting. Beyond mean-field calculations present the ground state as predominantly oblate29, whereas particle-core model calculations suggest predominantly prolate shape47,48. While the ground state is a mixture of spherical, oblate and prolate configurations with unknown admixtures, reduced α-decay widths cannot be unambiguously used as a tool to assign the configuration of the \({0}_{2}^{+}\) state in 186Pb. It is noteworthy, that many experimental findings providing supporting arguments for the proposed new assignments, such as lifetime measurements of prolate intruder states in 186Pb8,9 and in-beam γ-ray spectroscopy of 186Pb and 190Po17,18,48 nuclei, have been published after ref. 5.
Conclusions
The N = 104 neutron mid-shell nucleus 186Pb has been studied in a simultaneous in-beam γ-ray and electron spectroscopy experiment employing the sage spectrometer and the recoil-decay tagging technique at the Accelerator Laboratory of the University of Jyväskylä. The nucleus of interest is unique as it possesses a triplet of 0+ states within 700 keV of the ground state.
The existence of the excited 0+ states has been confirmed and their level energies have been determined with higher precision compared to earlier work through direct observation of the \({0}_{2}^{+}\to {0}_{1}^{+}\) and \({0}_{3}^{+}\to {0}_{1}^{+}\;E0\) transitions. The observation of the \({2}_{1}^{+}\to {0}_{2}^{+}\) transition allowed the collectivity of the transition feeding the band-head \({0}_{2}^{+}\) state to be assessed. The \(B(E2;{2}_{1}^{+}\to {0}_{2}^{+})\) value obtained suggests the \({0}_{2}^{+}\) state has a large prolate admixture. To the best of our knowledge, this marks the first observation of a transition feeding an excited predominantly prolate 0+ state in Pb nuclei.
Simultaneous measurement of conversion electrons and γ rays and intensity analysis revealed E0 components in the interband \({2}_{2}^{+}\to {2}_{1}^{+}\) and \({4}_{2}^{+}\to {4}_{1}^{+}\) transitions. This allowed for firm assignments of the \({2}_{2}^{+},{4}_{2}^{+},{6}_{2}^{+}\) and \({8}_{2}^{+}\) states that form part of the predominantly oblate band.
The two-level mixing calculations performed in the present work indicate that the 2+ and 4+ states are close to pure states with predominant components of ~90% and ~96%, respectively.
These findings call for a Coulomb excitation experiment to be performed. The spede spectrometer49 in conjunction with the Miniball spectrometer50 with radioactive ion beams from HIE-ISOLDE51 could provide information on diagonal and transitional matrix elements in 186Pb. Together with the data obtained in the present work, the monopole strength for interband transitions could be determined. Moreover, the intrinsic configuration of intruder states in this region could be probed in transfer-reaction experiments, e.g., at the ISOLDE Solenoidal Spectrometer52. The recoil-shadow method53 has been proposed54 to measure the lifetimes of the 0+ states and to extract the monopole strength of the E0(0+ → 0+) transitions. In addition, simultaneous in-beam γ-ray and electron spectroscopy could shed light on the configuration assignments of the 0+ band head states in the 188Pb nucleus.
Results obtained in the present work provide new insight into the nucleus at the heart of triple shape coexistence. It forms the basis for a systematic study of electric monopole transitions in the region and beyond with in-beam spectroscopic methods which have only now become possible.
Methods
The experiment was performed at the Accelerator Laboratory of the University of Jyväskylä. The use of sage+ritu+great30,55,56 instrumentation, described below, allowed for collection of simultaneous in-beam conversion-electron and γ-ray data employing the recoil-decay tagging method57,58. Nuclei of interest were produced via the 106Pd(83Kr,3n)186Pb reaction with a beam energy of 365 MeV. The 83Kr beam impinged on a 1-mg cm−2 thick target with intensities varying between 4 and 5 particle nA. During 108 h of beam on target, 6.34 × 105 α decays were recorded and correlated with 186Pb nuclei.
Prompt γ rays and conversion electrons were observed with the sage spectrometer. sage consisted of 10 eurogam Phase I and 24 eurogam Clover-type Compton-suppressed germanium detectors59 around the target for the detection of γ rays, and an annular segmented silicon detector placed upstream of the target for detection of conversion electrons. Transportation of conversion electrons was made by means of a solenoid coil that was operated with a current of 800 A. The δ-electron flux arising from atomic collisions in the target was suppressed by employing a high-voltage barrier at −35 kV. The full width at half maximum resolution values of 2.9 keV and 4.5 keV for 261 keV and 662 keV γ rays, respectively, and 11 keV and 14 keV for 173 keV and 447 keV conversion electrons, respectively, were determined from RDT γ-ray and electron singles data.
The ritu gas-filled separator was employed to separate fusion evaporation residues (referred to as recoils in the present work) from the primary beam particles. The great spectrometer, deployed at the focal plane of ritu, was used to identify recoils implanted in the double-sided silicon strip detector (DSSD) based on the deposited energy and characteristic α-decay energy. The maximum searching time between recoil implantation and the 186Pb α decay within the same DSSD pixel was 15 s which is equal to approximately three times the half-life of 186Pb (t1/2 = 4.83(3)s)60. The shallow implantation depth of recoils allowed for ~45% of α particles to escape from the DSSD without full energy deposition. Approximately 19% of the escaped α particles were detected by the PIN diode box of great installed upstream of the DSSD. The multiwire-proportional counter (MWPC) of great, immediately upstream of the PIN diode box, was used to measure the energy-loss of the recoils and scattered beam particles in a gas medium and to extract the time-of-flight information between the MWPC and the DSSD.
Data were collected using the fully digital Total Data Readout (TDR) acquisition system61. TDR operates without a common hardware trigger and events are time-stamped with a global 100 MHz clock. The GRAIN software package62 was used for temporal and spatial correlation of events. The final analysis was completed using the ROOT framework63.
Due to the long α-decay searching time, random correlations with 186Tl, 187Tl, 184Hg and 186Hg nuclei gave rise to contaminant events. The false correlation level was ~6% of all recoils. The contaminant nuclei are well known and their contribution was taken into account in the analysis.
Data availability
The data obtained in the present work and the corresponding metadata are available from https://doi.org/10.23729/a6444894-5fe7-4683-a9d5-a8ff1e28208b.
References
Heyde, K. & Wood, J. L. Shape coexistence in atomic nuclei. Rev. Mod. Phys. 83, 1467–1521 (2011).
Julin, R., Helariutta, K. & Muikku, M. Intruder states in very neutron-deficient Hg, Pb and Po nuclei. J. Phys. G: Nucl. Part. Phys. 27, R109–R139 (2001).
Bonn, J., Huber, G., Kluge, H.-J., Kugler, L. & Otten, E. W. Sudden change in the nuclear charge distribution of very light mercury isotopes. Phys. Lett. B 38, 308–311 (1972).
Julin, R., Grahn, T., Pakarinen, J. & Rahkila, P. In-beam spectroscopic studies of shape coexistence and collectivity in the neutron-deficient Z ≈ 82 nuclei. J. Phys. G: Nucl. Part. Phys. 43, 024004 (2016).
Andreyev, A. N. et al. A triplet of differently shaped spin-zero states in the atomic nucleus 186Pb. Nature 405, 430–433 (2000).
Van Duppen, P. et al. Observation of low-lying Jπ = 0+ states in the single-closed-shell nuclei 192−198Pb. Phys. Rev. Lett. 52, 1974–1977 (1984).
Van Duppen, P., Coenen, E., Deneffe, K., Huyse, M. & Wood, J. L. Low-lying Jπ=0+ states in 190,192Pb populated in the α-decay of 194,196Po. Phys. Lett. B 154, 354–357 (1985).
Grahn, T. et al. Collectivity and configuration mixing in 186,188Pb and 194Po. Phys. Rev. Lett. 97, 062501 (2006).
Grahn, T. et al. Lifetimes of intruder states in 186Pb, 188Pb and 194Po. Nucl. Phys. A 801, 83–100 (2008).
Wrzosek-Lipska, K. et al. Electromagnetic properties of low-lying states in neutron-deficient Hg isotopes: Coulomb excitation of 182Hg, 184Hg, 186Hg and 188Hg. Eur. Phys. J. A 55, 130 (2019).
Gaffney, L. P. et al. Shape coexistence in neutron-deficient Hg isotopes studied via lifetime measurements in 184,186Hg and two-state mixing calculations. Phys. Rev. C 89, 24307 (2014).
Dendooven, P. et al. Life time measurements of 0+ intruder states in 190,192,194Pb. Phys. Lett. B 226, 27–30 (1989).
De Witte, H. et al. Nuclear charge radii of neutron-deficient lead isotopes beyond N = 104 midshell investigated by in-source laser spectroscopy. Phys. Rev. Lett. 98, 112502 (2007).
Baxter, A. M. et al. Spectroscopy of 186Pb with mass identification. Phys. Rev. C 48, R2140–R2143 (1993).
Heese, J. et al. Evidence for low-lying prolate bands in 188Pb and 186Pb. Phys. Lett. B 302, 390–395 (1993).
Reviol, W. et al. Spectroscopy of 186Pb and 186Tl via evaporation residue detection. Phys. Rev. C 68, 054317 (2003).
Pakarinen, J. et al. Evidence for oblate structure in 186Pb. Phys. Rev. C 72, 011304 (2005).
Pakarinen, J. et al. Investigation of nuclear collectivity in the neutron mid-shell nucleus 186Pb. Phys. Rev. C 75, 014302 (2007).
May, F. R., Pashkevich, V. V. & Frauendorf, S. A prediction on the shape transitions in very neutron-deficient even-mass isotopes in the lead region. Phys. Lett. B 68, 113–116 (1977).
Nazarewicz, W. Variety of shapes in the mercury and lead isotopes. Phys. Lett. B 305, 195–201 (1993).
Tajima, N., Flocard, H., Bonche, P., Dobaczewski, J. & Heenen, P.-H. Diabatic effects in 186Pb: a generator-coordinate analysis. Nucl. Phys. A 551, 409–433 (1993).
Chasman, R. R., Egido, J. L. & Robledo, L. M. Persistence of deformed shapes in the neutron-deficient Pb region. Phys. Lett. B 513, 325–329 (2001).
Duguet, T., Bender, M., Bonche, P. & Heenen, P.-H. Shape coexistence in 186Pb: beyond-mean-field description by configuration mixing of symmetry restored wave functions. Phys. Lett. B 559, 201–206 (2003).
Fossion, R., Heyde, K., Thiamova, G. & Van Isacker, P. Intruder bands and configuration mixing in lead isotopes. Phys. Rev. C 67, 024306 (2003).
Bender, M., Bonche, P., Duguet, T. & Heenen, P.-H. Configuration mixing of angular momentum projected self-consistent mean-field states for neutron-deficient Pb isotopes. Phys. Rev. C 69, 064303 (2004).
Rodríguez-Guzmán, R. R., Egido, J. L. & Robledo, L. M. Beyond mean field description of shape coexistence in neutron-deficient Pb isotopes. Phys. Rev. C 69, 054319 (2004).
Egido, J. L., Robledo, L. M. & Rodríguez-Guzmán, R. R. Unveiling the origin of shape coexistence in lead isotopes. Phys. Rev. Lett. 93, 082502 (2004).
Nomura, K., Rodríguez-Guzmán, R., Robledo, L. M. & Shimizu, N. Shape coexistence in lead isotopes in the interacting boson model with a gogny energy density functional. Phys. Rev. C 86, 034322 (2012).
Yao, J. M., Bender, M. & Heenen, P.-H. Systematics of low-lying states of even-even nuclei in the neutron-deficient lead region from a beyond-mean-field calculation. Phys. Rev. C 87, 034322 (2013).
Pakarinen, J. et al. The SAGE spectrometer. Eur. Phys. J. A 50, 53 (2014).
Eberth, J. & Simpson, J. From Ge(Li) detectors to gamma-ray tracking arrays–50 years of gamma spectroscopy with germanium detectors. Prog. Part. Nucl. Phys. 60, 283–337 (2008).
Butler, P. A. et al. Electron spectroscopy using a multi-detector array. Nucl. Instrum. Methods A 381, 433–442 (1996).
Kankaanpää, H. et al. In-beam electron spectrometer used in conjunction with a gas-filled recoil separator. Nucl. Instrum. Methods A 534, 503–510 (2004).
Kibédi, T., Garnsworthy, A. B. & Wood, J. L. Electric monopole transitions in nuclei. Prog. Part. Nucl. Phys. 123, 103930 (2022).
Kibédi, T., Burrows, T. W., Trzhaskovskaya, M. B., Davidson, P. M. & Nestor, C. W. Evaluation of theoretical conversion coefficients using BrIcc. Nucl. Instrum. Methods A 589, 202–229 (2008).
Dracoulis, G. D. et al. Isomer bands, E0 transitions, and mixing due to shape coexistence in \({}_{82}^{188}{{{{{{{{\rm{Pb}}}}}}}}}_{106}\). Phys. Rev. C 67, 051301 (2003).
Church, E. L., Rose, M. E. & Weneser, J. Electric-monopole directional-correlation experiments. Phys. Rev. 109, 1299–1306 (1958).
Dowie, J. T. H., Kibédi, T., Eriksen, T. K. & Stuchbery, A. E. Table of electronic factors for E0 electron and electron–positron pair conversion transitions. Atomic Data Nuclear Data Tables 131, 101283 (2020).
Wood, J. L., Zganjar, E. F., De Coster, C. & Heyde, K. Electric monopole transitions from low energy excitations in nuclei. Nucl. Phys. A 651, 323–368 (1999).
Hellemans, V., De Baerdemacker, S. & Heyde, K. Configuration mixing in the neutron-deficient 186−196Pb isotopes. Phys. Rev. C 77, 064324 (2008).
Zhang, W. Q. et al. First observation of a shape isomer and a low-lying strongly-coupled prolate band in neutron-deficient semi-magic 187Pb. Phys. Lett. B 829, 137129 (2022).
Smirnova, N., Heenen, P.-H. & Neyens, G. Self-consistent approach to deformation of intruder states in neutron-deficient Pb and Po. Phys. Lett. B 569, 151–158 (2003).
Möller, P., Sierk, A. J., Bengtsson, R., Sagawa, H. & Ichikawa, T. Global calculation of nuclear shape isomers. Phys. Rev. Lett. 103, 212501 (2009).
Van de Vel, K. et al. Fine structure in the α decay of 188,192Po. Phys. Rev. C 68, 054311 (2003).
Van Duppen, P., Huyse, M. & Wood, J. L. Mixing of intruder and normal states in Pb nuclei. J. Phys. G: Nucl. Part. Phys. 16, 441–450 (1990).
Marsh, B. A. et al. Characterization of the shape-staggering effect in mercury nuclei. Nat. Phys. 14, 1163–1167 (2018).
Oros, A. M. et al. Shape coexistence in the light Po isotopes. Nucl. Phys. A 645, 107–142 (1999).
Van de Vel, K. et al. In-beam γ-ray spectroscopy of 190Po: First observation of a low-lying prolate band in Po isotopes. Eur. Phys. J. A 17, 167–171 (2003).
Papadakis, P. et al. The SPEDE spectrometer. Eur. Phys. J. A 54, 42 (2018).
Warr, N. et al. The Miniball spectrometer. Eur. Phys. J. A 49, 40 (2013).
Kadi, Y. et al. Post-accelerated beams at ISOLDE. J. Phys. G: Nucl. Part. Phys. 44, 084003 (2017).
Tang, T. L. et al. First exploration of neutron shell structure below lead and beyond n = 126. Phys. Rev. Lett. 124, 062502 (2020).
Backe, H. et al. In-beam spectroscopy of low energy conversion electrons with a recoil shadow method — a new possibility for subnanosecond lifetime measurements. Z. Phys. A 285, 159–169 (1978).
Pakarinen, J. et al. Recoil-shadow electron spectroscopy of low-lying 0+ states in 186Pb and 194Po. Proposal to the JYFL-PAC (2006).
Leino, M. et al. Gas-filled recoil separator for studies of heavy elements. Nucl. Instrum. Methods B 99, 653–656 (1995).
Page, R. D. et al. The GREAT spectrometer. Nucl. Instrum. Methods B 204, 634–637 (2003).
Paul, E. S. et al. In-beam γ-ray spectroscopy above \({}^{100}{{{{{{{\rm{Sn}}}}}}}}\) using the new technique of recoil decay tagging. Phys. Rev. C 51, 78–87 (1995).
Simon, R. S. et al. Evidence for nuclear shape coexistence in 180Hg. Z. Phys. A 325, 197–202 (1986).
Beck, F. A. EUROBALL: Large gamma ray spectrometers through european collaborations. Prog. Part. Nucl. Phys. 28, 443–461 (1992).
Wauters, J. et al. Alpha decay of 186Pb and 184Hg: The influence of mixing of 0+ states on α-decay transition probabilities. Phys. Rev. C 50, 2768–2773 (1994).
Lazarus, I. et al. The GREAT triggerless total data readout method. IEEE Trans. Nucl. Sci. 48, 567–569 (2001).
Rahkila, P. Grain—a java data analysis system for total data readout. Nucl. Instrum. Methods A 595, 637–642 (2008).
Brun, R. & Rademakers, F. ROOT — an object oriented data analysis framework. Nucl. Instrum. Methods A 389, 81–86 (1997).
Acknowledgements
This work has been supported through the Academy of Finland under the Finnish Centre of Excellence Programme 2012–2017 (Nuclear and Accelerator Based Physics Programme at JYFL), through the Academy Research Fellow funding (contract number 257562) and the European Gamma-Ray Spectroscopy pool. Support from the Science and Technology Facilities Council (UK) and EU 7th framework programme Integrating Activities - Transnational Access, Project No. 262010 (ENSAR) are acknowledged. AH would like to thank the Slovak Research and Development Agency under contract No. APVV-20-0532, and Slovak grant agency VEGA (contract No. 2/0067/21).
Author information
Authors and Affiliations
Contributions
J. Pakarinen, P. Papadakis and D.G.J. prepared the proposal for the experiment. J. Pakarinen, P. Papadakis, J. Sorri, M.S. and D.M.C. set up the instrumentation. K.A., H.B., D.M.C., P.J.D., T.G., P.T.G., J. Henderson, A.H., R-.D.H., J. Hilton, U.J., D.G.J., D.T.J., R.J., S.J., T.K., J.K., G.J.L., M.L., J.L., C.G.M., O.N., R.D.P., J. Pakarinen, P. Papadakis, E.P., J. Partanen., P. Peura, P. Rahkila, J.R., P. Ruotsalainen, M.S., J. Sarén, C.S., J. Sorri, S.S., J.U., A.W. and R.W. monitored the detector, data acquisition and beam systems. J.O., J. Pakarinen, P. Papadakis, R.J. and R.D.P. carried out the data analysis and interpretation of the data, and J. Pakarinen, J.O., P. Papadakis, R.J., R.D.P., K.A., M.L., P.T.G., R-.D.H., E.P., P. Peura, J. Sorri, J.U. and R.W. prepared the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Communications Physics thanks the anonymous reviewers for their contribution to the peer review of this work.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Ojala, J., Pakarinen, J., Papadakis, P. et al. Reassigning the shapes of the 0+ states in the 186Pb nucleus. Commun Phys 5, 213 (2022). https://doi.org/10.1038/s42005-022-00990-4
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s42005-022-00990-4
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.