scieee AI-readable full text Open interactive document viewer

Two-particle transfer processes as a signature of shape phase transition in Zirconium isotopes

Lay Valera, José Antonio; Vitturi, Andrea; Fortunato, Lorenzo; Tsunoda, Yusuke; Togashi, Tomoaki; Otsuka, T.

Abstract

We explore two-particle transfer reactions as a crucial probe of the occurrence of shape coexistence in shape phase transitions. The (t,p) reactions to the ground state and to excited 0+ states are calculated for the isotope chain of even-even Zirconium isotopes starting from stable nuclei up to beyond current experimental limits. Two-particle spectroscopic factors derived from Monte Carlo Shell Model calculations are used, together with the sequential description of the two-particle transfer reaction mechanism. The calculation shows a clear signature for a shape phase transition between 98Zr and 100Zr, which displays coexistence of a deformed ground state with an excited spherical 0+ state. Furthermore, we show that there is a qualitative difference with respect to the case of a normal shape phase transition that can be discriminated with two-neutron transfer reactions.

Full text

Physics Letters B 838 (2023) 137719 Contents lists available at ScienceDirect Physics Letters B journal homepage: www.elsevier.com/locate/physletb Two-particle transfer processes as a signature of shape phase transition in Zirconium isotopes J.A. Lay a,b,∗, A. Vitturi c,d, L. Fortunato c,d, Y. Tsunoda e,f, T. Togashi f, T. Otsuka f,g,h,i aDepartamento de Física Atómica, Molecular y Nuclear, Facultad de Física, Universidad de Sevilla, Apartado 1065, E-41080 Sevilla, Spain bInstituto Interuniversitario Carlos I de Física Teórica y Computacional (iC1), Apdo. 1065, E-41080 Sevilla, Spain cDipartimento di Fisica e Astronomia “Galileo Galilei”, Univ. Padova, via Marzolo, 8, I-35131 Padova, Italy dINFN, Sezione di Padova, via Marzolo, 8, I-35131 Padova, Italy eCenter for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki 305-8577, Japan fCenter for Nuclear Study, University of Tokyo, 7-3-1 Hongo, Bunkyo, Tokyo 113-0033, Japan gDepartment of Physics, University of Tokyo, Hongo, Bunkyo-ku, Tokyo 113-0033, Japan hRIKEN Nishina Center, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan iInstituut voor Kern- en Stralingsfysica, KU Leuven, B-3001 Leuven, Belgium a r t i c l e i n f o a b s t r a c t Article history: Received 27 October 2022 Received in revised form 6 December 2022 Accepted 20 January 2023 Available online 24 January 2023 Editor: B. Balantekin We explore two-particle transfer reactions as a crucial probe of the occurrence of shape coexistence in shape phase transitions. The (t,p) reactions to the ground state and to excited 0+states are calculated for the isotope chain of even-even Zirconium isotopes starting from stable nuclei up to beyond current experimental limits. Two-particle spectroscopic factors derived from Monte Carlo Shell Model calculations are used, together with the sequential description of the two-particle transfer reaction mechanism. The calculation shows a clear signature for a shape phase transition between 98Zr and 100Zr, which displays coexistence of a deformed ground state with an excited spherical 0+state. Furthermore, we show that there is a qualitative difference with respect to the case of a normal shape phase transition that can be discriminated with two-neutron transfer reactions. ©2023 The Author(s). Published by Elsevier B.V. This is an open access article under the CC BY license (http://creativecommons .org /licenses /by /4 .0/). Funded by SCOAP3. 1. Introduction The phases and transitions between them are prominent features of many-body systems. The atomic nucleus, in some cases, clearly exhibits these features. This is a unique and precious situation because the nucleus is an isolated system, and the transition occurs as a consequence of certain changes of its ingredients rather than due to a change of external environment. We note that because the nucleus is a finite quantal system, quantum phase transitions (QPT) may occur [1–4]. The phase transition can take place in different ways. Typically the change occurs as a function of control parameters as the excitation energy (i.e. the temperature in a thermodynamical framework) or the angular momentum. But equally important are the transitions taking place in the shape of the ground state along a chain of isotopes (or isotones), where the discrete control parameter is the number of neutrons (or protons). Meaningful order *Corresponding author at: Departamento de Física Atómica, Molecular y Nuclear, Facultad de Física, Universidad de Sevilla, Apartado 1065, E-41080 Sevilla, Spain. E-mail address: [email protected] (J.A. Lay). parameters systematically used in such shape evolution are, in the case of even-even nuclei, the energy of the first 2+state, the ratio E4/E2and the strength of the electromagnetic E2 transition connecting ground state and the first excited 2+state. However, both the excitation energies and the E2 transition depend also on the structure of other states. It is of great and broad interest whether and how one can directly see the structure change between the ground states. Two-particle transfer processes, e.g. reactions populating the 0+ states, can however provide a complementary but crucial clear-cut signature of the occurrence of the phase transition. In particular, in the presence of a sharp or abrupt transition, one expects a sudden weakening of the usual dominance of ground-to-ground A → A+2 transitions and a corresponding abnormally strong population of one (or more) excited 0+states. The behavior of the full pair response (as defined below) should indicate not only the occurrence of a shape phase transition, but also the nature of this transition. A correct description of the process implies, besides a proper reaction model, a proper microscopic description of the nuclear wave functions and in this Letter we present novel microscopic calculations of the pair transfer process based on nuclear structure inputs obtained within the Monte Carlo Shell Model for the https://doi.org/10.1016/j.physletb.2023.137719 0370-2693/©2023 The Author(s). Published by Elsevier B.V. This is an open access article under the CC BY license (http://creativecommons .org /licenses /by /4 .0/). Funded by SCOAP3. J.A. Lay, A. Vitturi, L. Fortunato et al. Physics Letters B 838 (2023) 137719 chain of Zirconium isotopes [5]. Before that, for the sake of clarity, we discuss schematic behaviors that can be obtained within simplified algebraically-based approaches [6,7] and are meant as a guideline. After a two-year hiatus, we have resumed this work, that only made it to the arXiv without being published [8], because we think that it is still very timely and interesting and we believe that other research groups, in parallel [9]or right after [10,11]our pioneering efforts in this field [12], have come to the same or very similar conclusions and that these conclusions are relevant for disentangling the onset of the usual QPT scenario from the more involved interplay of intruder configurations. In particular, the comparison of our results (see pictures and table) with Fig. 1 in Ref. [9]reveals differences, namely that there is a clear case for opening up new experimental investigations: we find that the fourth 0+states takes up all the transfer strengths, but the shell model calculations suggest that the second 0+state in 100Zr is the analogous to the ground state of 98Zr. It would be highly desirable to perform twoneutron transfer reactions with sufficient energy resolution and analyze the outcome along the lines proposed in the present paper. The intensity for the two-neutron transfer according to the structure model of [9]is presented in [10] where one of the present authors is also involved. However, this is, to our knowledge, the first time that second-order Distorted Wave Born Approximation (DWBA) calculations of the actual cross section for two-neutron transfer are discussed. The problem of how the intruder configurations come down in energy and alter the onset of a quantum shape phase transition is clearly an open one and an interesting one, because the scheme can be repeated in other shells. Several improvements to Figures, comments and new calculations have been added to the present final version of the manuscript. 2. Schematic description We first consider the case of a series of isotopes described within the Interacting Boson Model (IBM), with a model Hamiltonian that abruptly performs a transition from U(5) to SU(3), i.e. from a spherical behavior to the situation of deformed axial symmetry, so as to simulate the structure of the ground state of the MCSM calculation [5]. As an example, considering a core with A =90, we can assume the U(5) Hamiltonian for systems with a number of bosons, N, ranging from 1 (92Zr) to 4 (98Zr) and the SU(3) Hamiltonian starting from N=5 (100Zr). The system will therefore be “spherical” up to N=4 and “axially-symmetric deformed” from N=5 on. This choice of core and valence bosons is not unique (see for example Ref. [13]), but a different core would not qualitatively change the physical picture within the Interacting Boson Model IBM-1 model. Within the IBM the pair creation operator, in leading order, is given just by the s†operator and one can evaluate the corresponding pair addition intensities obtained by taking the square of the matrix element connecting the ground state in system N with the ground and excited state in N+1. See for example Eq. (2.1) and Eq. (2.2) of Ref. [6]. These intensities are shown in the upper frame of Fig. 1and display a clear “anomaly” for the transition across the change of phase, i.e. the one connecting a “spherical” system with a “deformed” one. The full pair response is shown in the three lower frames panels of Fig. 1for a pair addition within the spherical phase (N=3→4), across the phase transition (N=4→5) and within the deformed phase (N=5→6). The pair strength, normally concentrated in the ground to ground transition, appears completely fragmented in correspondence of the critical point as seen in Fig. 1. A different physical scenario is that of shape coexistence, where different shape phases occur within the same nucleus at similar excitation energies [14]. We may face the situation of a (slow or Fig. 1. Upper panel: two-particle transfer intensities to the ground state (solid) and to the summed excited 0+states (dashed) in the case of a sequence of isotopes making the sharp transition from U(5) to SU(3) (cf. text). In the lower panels, the full pair response spectrum is shown for the indicated cases, showing fragmentation of the strength at the critical point. The excitation energy E∗is in arbitrary units since the actual energies of the different states will vary on each particular case. Fig. 2. Two-particle transfer intensities to the ground state and to the excited 0+ state in the case of a sequence of isotopes in the presence of shape coexistence. From N=4 to N=5 there is the exchange of the spherical ground-state configuration with the “intruder” deformed configuration. In the lower panels, the full pair response spectrum is shown for the indicated cases. As in Fig. 1, the excitation energy E∗is in arbitrary units. rapid) progressive mixing of the spherical and deformed phases, eventually leading to the interchange of the dominant phase in the ground state. Again we can have a first guess of the consequences of this situation on the pair-transfer processes within a simplified IBM-like framework. Following the idea of ref. [15,16]we can assume for each system a possible mixing of a “spherical” state, characterized by N valence bosons, obtained within an IBM U(5) Hamiltonian with another “deformed” 0+state obtained within a SU(3) Hamiltonian with N+2 bosons, microscopically originated by a 2p-2h core excitation. In this case the pair creation operator will be given, in leading order, by (s†+s), since we can either add a valence-like boson or destroy the “hole-like” boson. Assuming a sharp transition with increasing number of particles from a fully spherical ground state to a fully deformed ground state we obtain, still with a transition taking place passing from N=4 to N=5 (cf. inset in Fig. 2), the pair transfer intensities shown in Fig. 2. As in the previous case a clear discontinuity appears at the tran- 2 J.A. Lay, A. Vitturi, L. Fortunato et al. Physics Letters B 838 (2023) 137719 Table 1 Two-particle transfer amplitudes for the different reactions connecting even-even Zr isotopes, for the most relevant single-particle orbits. For each case the largest component is evidenced. Notice that, depending on the phase convention, an extra (−1)factor should be added to these amplitudes. 90-92gs 92-94gs 94-96gs 96-98gs 98-100gs 98-100(0+ 4) 100-102gs 1d5/20.74 0.86 0.86 0.13 ∼0.0 0.16 0.08 2s1/20.10 0.08 0.10 0.90 ∼0.0 0.16 0.05 1d3/20.13 0.18 0.16 0.07 ∼0.0 0.90 0.04 0h11/20.22 0.20 0.19 0.08 ∼0.0 0.14 0.55 sition point. However, at variance with the previous case, the pair strength is always practically concentrated in a single state, without the fragmentation illustrated in Fig. 1(lower central panel). Therefore, while the discontinuity just signals the occurrence of a phase transition, it is precisely the presence or the absence of this fragmentation that characterizes the physical scenario: normal phase transition vs. phase transition driven by coexistence. 3. Numerical results We move now from schematic models to a fully microscopic calculation, for both reaction mechanism and structure. We take the case of (t,p) reactions on even-A Zirconium isotopes, where experimental data at E=20 MeV are available [19]at least for the lighter systems, i.e. up to 96Zr(t,p). Novel interest on Zirconium isotopes has arisen from the recent shell model calculations [5] that indicate a possible case of shape coexistence in these nuclei with a sharp transition occurring between 98Zr and 100Zr. The situation seems to resemble the schematic case of shape coexistence displayed in Fig. 2, although the amount of the particle-hole excitations exceeds the, proton in this case, 2p-2h picture as emphasized in [5]. On the other hand, the amount of the particle-hole excitations may not matter as the normal (intruder) states are connected by the present transfer reactions to the normal (intruder) states. In other words, the cross section of transferring two neutrons will be large if the rest of nucleons remain in the same configuration but it will not depend if this configuration is 0p-0h, neutron or proton 4p-4h, or even a mixture of all possible ones. The deformation of Zr isotopes were already studied within the Shell Model in Ref. [20] and within the IBM in Ref. [10,11,21,22]. Recent experimental studies can be found searching for evidences of shape coexistence in 98Zr [23,24] and 96Zr [25]. We have therefore calculated the two-particle transfer probabilities across the phase transition up to 100Zr(t,p) to the ground and excited 0+states. In parallel with the detailed microscopic structure description described above, also the reaction process has been described in microscopic terms. In particular, we have performed second order DWBA calculations with the code FRESCO [26, 27]. Therefore the reaction mechanism includes the “correlated” sequential single-particle transfer through all intermediate states in the A+1 odd system, the simultaneous transfer of the two neutrons and non-orthogonality terms. Preliminary calculations with only the sequential transfer can be found in [12]. Single-particle wave functions for the construction of the single-particle form factors have been generated within a Saxon-Woods potential adjusted to yield the proper single-particle energy. Optical potential for protons is taken from Chapel-Hill 89 [28], for deuterons is taken from Lohr-Haeberli as retrieved from [17], and for tritons from Bechetti- Greenless [17]. Such choice, together with the structure input from MCSM calculations allows for a good reproduction of the available data from [19]. As an example we show the case of 96Zr(t,p)98Zr in Fig. 3(a) together with other reasonable options for the different optical potentials. Fig. 3(a) also shows several calculations with different choices of optical potentials as explain in the caption of the figure. It can be seen that the overall description is preserved and it has been tested that present conclusions are not affected by the choice of optical potentials. Fig. 3. Upper panel: Angular distribution for the reaction 96Zr(t,p)98Zr at 20 MeV when calculated with the present choice of optical potential (black thick solid line, see text). In red dashed, the result of changing the proton optical potential by the global parametrization of Menet et al. as retrieved from [17]; in blue dot-dashed, the result of changing the deuteron potential by that of Perey-Perey [17]; and in green dotted the result of changing the triton potential by that of X. Li et al. [18]. Lower panel: Angular distributions when the two particles are transferred in a pure configuration. The thick solid line gives the result in the case of correlated wave function according to the value of two-particle amplitudes, the largest being given in the table I. Experimental values are taken from [19]. This reaction mechanism generates a dynamical dependence on each specific orbit on which the pair is transferred. The origin of this behavior in the case of (t, p)or (p, t)reactions comes from the different content of (0s)relative n −nmotions, that is associated to each two-particle configuration (cf. how the Talmi- Moshinsky brackets enter into the calculation of reaction probabilities [29]). The transfer probabilities become therefore sensitive not only to the value of the “global” pair strength, but also to the details of microscopic wave functions [30]. This is better evidenced in Fig. 3(b) where the cross sections associated with single-particle orbits are reported in the case of the 96Zr(t,p)98Zr reaction. The collective effects in the pair transfer process comes from the correlations present in both initial and final states that induce a coherent 3 J.A. Lay, A. Vitturi, L. Fortunato et al. Physics Letters B 838 (2023) 137719 Fig. 4. Differential cross sections at the maximum around 40◦for (t,p) reactions on the different isotopes and specific final states. Experimental values, when available, are also given [19]. Total cross section would be preferable although not available. Explicit values for error bars are not given in [19]. and constructive interference of all the sequential paths (cf. refs. [29,31] and references therein). In our case this coherence is obtained by using the two-particle spectroscopic amplitudes provided by the Monte Carlo Shell Model calculation [5]. The most important two-particle addition amplitudes or two-neutron amplitudes (TNA) are reported in Table 1, but also the smaller contributions from the other orbits included in the model space have been used in the reaction calculation. The constructive effect of the residual pairing-like interaction is evidenced by the enhancement of the correlated cross section (also shown in Fig. 3(b), full black line) with respect the single particle estimates. The comparison with the experimental data is not straightforward since explicit values are not given in the only available reference [19]. In Fig. 4that summarizes the results for the full sequence of transfer reactions, the black bars correspond to experimental data that we obtained by combining relative strengths from Fig. 10 of Ref. [19]with the differential cross section at the first maximum (excluding θ=0◦) given in Table 4 of the same reference. The overall behavior of our calculations (blue bars in Fig. 4) reproduces the experimental trend, when this is available. As expected from the amplitudes given in the table, in the case of 98Zr(t,p)100Zr the calculation predicts a large population of the fourth 0+state in 100Zr, which displays a “spherical” behavior as the target 98Zr(gs). Such fourth 0+state appears in the NCSM calculation at an excitation energy of 1.47 MeV [5]. Continuing beyond the critical point, we predict again a relatively weak population of the ground state in 100Zr(t,p)102Zr, although the reaction connects now two deformed systems with practically the same deformation. In this case, this is not due to small two-particle spectroscopic amplitudes (cf. the rather large spectroscopic amplitude associated with the h11/2orbit in Table 1), but to the reaction mechanism that does not favor the transfer of a pair to the h11/2 single-particle level characterized by a large single-particle orbital angular momentum, having a small overlap with a 0swavefunc- tion in the relative n −nmotion. In this sense, reaction dynamics modulates the pattern of Figs. 1and 2, but such modulation is well under control in second-order DWBA. The fact that the small cross section comes from the small overlap of the particular orbit to the 0srelative n −nmotion can be shown just multiplying the intensity to the ground state in Fig. 1by the square of the Talmi-Moshinsky bracket corresponding to the largest amplitude of each case, renormalizing all values to one experimental cross section, 94Zr(t,p)96Zr in this case (red dashed line in Fig. 4). It can be seen that this approximation follows the experimental trend and explains reasonably well the small cross section in 100Zr(t,p)102Zr. The same information is encoded in the TNA: if one multiplies the largest two-neutron transfer amplitude by the square of the corresponding Talmi-Moshinsky bracket, (green dotted line) the same behavior is found as can be seen in Fig. 4(green dotted line) also renormalizing to the case of 94Zr(t,p)96Zr. 4. Summary and conclusions In summary, reaction calculations for (t, p)processes between even-even Zr isotopes have been performed using two-particle transition amplitudes provided by state-of-the-art MCSM calculations. The outcome reproduces the trends of available experimental data and indicates a sharp change in the pattern of twoparticle cross-sections between 98Zr and 100Zr: the ground-to- ground cross-section drops dramatically, in correspondence with a strong population of an excited 0+state. This confirms that twoparticle transfer reactions can be nicely used as an additional, but crucial probe to pinpoint the occurrence and the nature of quantum shape phase transitions. The higher resemblance to the IBM with configuration mixing appears to be consistent with the actual situation in the MCSM calculation. Finally, we stress that the pair transfer reaction is a very crucial tool to look into the structure of the wave functions of unstable nuclei, such us 98Zr and 100Zr, and it will play more central roles in the near future as the Rare- Isotope beams are becoming more intense, being more suitable for this kind of reactions. Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. Data availability The authors are unable or have chosen not to specify which data has been used. Acknowledgements The MCSM calculations were performed on the K computer at RIKEN AICS (hp150224, hp160211). This work was supported in part by the HPCI Strategic Program (The origin of matter and the universe) and “Priority Issue on post-K computer” (Elucidation of the Fundamental Laws and Evolution of the Universe) from MEXT and JICFuS. J.A.L. acknowledges funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement N. 654002, from grants PID2020-114687GB-I00 and FIS2017-88410-P by MCIN/AEI/10.13039/501100011033, the project PAIDI 2020 with Ref. P20_01247 by the Consejería de Economía, Conocimiento, Empresas y Universidad, Junta de Andalucía (Spain), and by ERDF A way of making Europe. References [1] S. Sachdev, Quantum Phase Transitions, Cambridge Univ. Press, Cambridge, 2011. [2] R. Casten, Shape phase transitions and critical-point phenomena in atomic nuclei, Nat. Phys. 2 (2006) 811, https://doi .org /10 .1038 /nphys451. [3] P. Cejnar, J. Jolie, R.F. Casten, Quantum phase transitions in the shapes of atomic nuclei, Rev. Mod. Phys. 82 (2010) 2155, https://doi .org /10 .1103 /RevModPhys .82 . 2155. [4] L. Fortunato, Quantum phase transitions in algebraic and collective models of nuclear structure, Prog. Part. Nucl. Phys. 121 (2021) 103891, https://doi .org / 10 .1016 /j .ppnp .2021.103891, https://www.sciencedirect .com /science /article /pii / S0146641021000508. 4 J.A. Lay, A. Vitturi, L. Fortunato et al. Physics Letters B 838 (2023) 137719 [5] T. Togashi, Y. Tsunoda, T. Otsuka, N. Shimizu, Phys. Rev. Lett. 117 (2016) 172502, https://doi .org /10 .1103 /PhysRevLett .117.172502, https://link.aps .org / doi /10 .1103 /PhysRevLett .117.172502. [6] R. Fossion, C.E. Alonso, J.M. Arias, L. Fortunato, A. Vitturi, Phys. Rev. C 76 (2007) 014316, https://doi .org /10 .1103 /PhysRevC .76 .014316, https://link.aps . org /doi /10 .1103 /PhysRevC .76 .014316. [7] Y. Zhang, F. Iachello, Phys. Rev. C 95 (2017) 034306, https://doi .org /10 .1103 / PhysRevC .95 .034306, https://link.aps .org /doi /10 .1103 /PhysRevC .95 .034306. [8] J.A. Lay, A. Vitturi, L. Fortunato, Y. Tsunoda, T. Togashi, T. Otsuka, arXiv:1905 . 12976, https://arxiv.org /abs /1905 .12976. [9] J.E. García-Ramos, K. Heyde, Quest of shape coexistence in zr isotopes, Phys. Rev. C 100 (2019) 044315, https://doi .org /10 .1103 /PhysRevC .100 .044315, https://link.aps .org /doi /10 .1103 /PhysRevC .100 .044315. [10] J. García-Ramos, J. Arias, A. Vitturi, Two-neutron transfer reactions as a tool to study the interplay between shape coexistence and quantum phase transitions, Chin. Phys. C 44 (2020) 124101, https://doi .org /10 .1088 /1674 -1137 /abb4ca. [11] N. Gavrielov, A. Leviatan, F. Iachello, Zr isotopes as a region of intertwined quantum phase transitions, Phys. Rev. C 105 (2022) 014305, https://doi .org /10 .1103 /PhysRevC .105 .014305, https://link.aps .org /doi /10 .1103 / PhysRevC .105 .014305. [12] A. Vitturi, L. Fortunato, I. Inci, J.A. Lay, JPS Conf. Proc. 23 (2018) 01201, https:// doi .org /10 .7566 /JPSCP.23 .012013. [13] S. Lalkovski, P.V. Isacker, Phys. Rev. C 79 (2009) 044307, https://doi .org /10 .1103 / PhysRevC .79 .044307, https://link.aps .org /doi /10 .1103 /PhysRevC .79 .044307. [14] A. Frank, P. Van Isacker, C.E. Vargas, Phys. Rev. C 69 (2004) 034323, https://doi . org /10 .1103 /PhysRevC .69 .034323, https://link.aps .org /doi /10 .1103 /PhysRevC .69 . 034323. [15] P.D. Duval, B.R. Barrett, Nucl. Phys. A 376 (1982) 213, https://doi .org /10 . 1016 /0375 -9474(82 )90061 -6, http://www.sciencedirect .com /science /article /pii / 0375947482900616. [16] J.E. García-Ramos, K. Heyde, Phys. Rev. C 92 (2015) 034309, https://doi . org /10 .1103 /PhysRevC .92 .034309, https://link.aps .org /doi /10 .1103 /PhysRevC .92 . 034309. [17] C. Perey, F. Perey, Compilation of phenomenological optical-model parameters 1954–1975, At. Data Nucl. Data Tables 17 (1976) 1–101, https:// doi .org /10 .1016 /0092 -640X(76 )90007 -3, https://www.sciencedirect .com / science /article /pii /0092640X76900073. [18] X. Li, C. Liang, C. Cai, Global triton optical model potential, Nucl. Phys. A 789 (2007) 103–113, https://doi .org /10 .1016 /j .nuclphysa .2007.03 .004, https://www. sciencedirect .com /science /article /pii /S0375947407002291. [19] E. Flynn, J. Beery, A. Blair, Nucl. Phys. A 218 (1974) 285, https://doi .org /10 . 1016 /0375 -9474(74 )90005 -0, http://www.sciencedirect .com /science /article /pii / 0375947474900050. [20] P. Federman, S. Pittel, Phys. Rev. C 20 (1979) 820, https://doi .org /10 .1103 / PhysRevC .20 .820, https://link.aps .org /doi /10 .1103 /PhysRevC .20 .820. [21] N. Gavrielov, A. Leviatan, F. Iachello, Intertwined quantum phase transitions in the zr isotopes, Phys. Rev. C 99 (2019) 064324, https://doi .org /10 .1103 / PhysRevC .99 .064324, https://link.aps .org /doi /10 .1103 /PhysRevC .99 .064324. [22] José-Enrique García-Ramos, Kris Heyde, EPJ Web Conf. 178 (2018) 05005, https://doi .org /10 .1051 /epjconf /201817805005. [23] P. Singh, et al., Phys. Rev. Lett. 121 (2018) 192501, https://doi .org /10 . 1103 /PhysRevLett .121.192501, https://link.aps .org /doi /10 .1103 /PhysRevLett .121. 192501. [24] V. Karayonchev, J. Jolie, A. Blazhev, A. Dewald, A. Esmaylzadeh, C. Fransen, G. Häfner, L. Knafla, J. Litzinger, C. Müller-Gatermann, J.-M. Régis, K. Schomacker, A. Vogt, N. Warr, A. Leviatan, N. Gavrielov, Tests of collectivity in 98Zr by absolute transition rates, Phys. Rev. C 102 (2020) 064314, https://doi .org /10 .1103 / PhysRevC .102 .064314, https://link.aps .org /doi /10 .1103 /PhysRevC .102 .064314. [25] C. Kremer, S. Aslanidou, S. Bassauer, M. Hilcker, A. Krugmann, P. von Neumann-Cosel, T. Otsuka, N. Pietralla, V.Y. Ponomarev, et al., Phys. Rev. Lett. 117 (2016) 172503, https://doi .org /10 .1103 /PhysRevLett .117.172503, https:// link.aps .org /doi /10 .1103 /PhysRevLett .117.172503. [26] I.J. Thompson, Comput. Phys. Rep. 7 (1988) 167, https://doi .org /10 .1016 / 0167 -7977(88 )90005 -6, http://www.sciencedirect .com /science /article /pii / 0167797788900056. [27] I.J. Thompson, Reaction Mechanisms of Pair Transfer, World Scientific, 2013, p. 455, Ch. 34, https://www.worldscientific .com /doi /abs /10 .1142 / 9789814412490 _0034. [28] R. Varner, W. Thompson, T. McAbee, E. Ludwig, T. Clegg, A global nucleon optical model potential, Phys. Rep. 201 (1991) 57–119, https:// doi .org /10 .1016 /0370 -1573(91 )90039 -O, https://www.sciencedirect .com / science /article /pii /037015739190039O. [29] R. Broglia, O. Hansen, C. Riedel, Adv. Nucl. Phys. 6 (1973) 287. [30] J.A. Lay, L. Fortunato, A. Vitturi, Phys. Rev. C 89 (2014) 034618, https://doi . org /10 .1103 /PhysRevC .89 .034618, https://link.aps .org /doi /10 .1103 /PhysRevC .89 . 034618. [31] W. von Oertzen, A. Vitturi, Rep. Prog. Phys. 64 (2001) 1247, https://doi .org /10 . 1088 /0034 -4885 /64 /10 /202. 5