Abstract
In the framework of mean-field based transport approaches, we discuss recent results concerning collective motion and low-energy heavy ion reactions involving neutron-rich systems. We focus on aspects which are particularly sensitive to the isovector terms of the nuclear effective interaction and the corresponding symmetry energy. As far as collective excitations are concerned, we discuss the mixed nature of dipole oscillations in neutron-rich systems. On the other hand, for reactions close to the Coulomb barrier, we investigate the structure of pre-equilibrium collective dipole oscillations, focusing on their sensitivity to the symmetry energy behavior below normal density. Nucleon emission is also considered within the same context. The possible impact of other relevant terms of the nuclear effective interaction on these mechanisms is also examined. From this analysis we expect to put further constraints on the nuclear Equation of State, of crucial importance also in the astrophysical context.
1. Introduction
Collective patterns exhibited by complex systems can bear important information on relevant properties of the particle interaction. In nuclei, the investigation of the giant resonances, whose collective nature is well established, is therefore of primary importance []. A prominent example in this context is the giant dipole resonance (GDR), which can be described in terms of protons and neutrons oscillating as a whole against each other.
Stimulated by the advent of new radioactive beam facilities, a large amount of research has been devoted in recent years to the features of unstable nuclei and their collective multipole response. In the case of nuclei with some neutron excess, a strong fragmentation of strength has been observed in the isovector dipole response, mainly located at lower energy with respect to the GDR [–]. These low-lying excitations, which are referred in the literature as Pygmy Dipole Resonance (PDR), have been the object of intense discussion [–], proving to be, like the GDR, an important probe of crucial information of the nuclear effective interaction, especially concerning its isovector component and the corresponding contribution to the Equation of State (EoS) [–], namely the symmetry energy.
Collective oscillations of neutrons against protons might occur also in low-energy reactions involving charge-asymmetric systems, at least during the pre-equilibrium stage. If the N/Z ratios of the reaction partners are appreciably different, then neutron and proton centers of mass of the involved composite system do not coincide in the early phase of the fusion path and charge equilibration mechanisms take place. As a result, together with the incoherent exchange of nucleons between the reacting ions, a dynamical dipole (DD) mode, also known as pre-equilibrium GDR [–], is observed along the symmetry axis of the dinuclear system.
Since the transient composite system might experience large prolate deformation with respect to the equilibrium configuration of the final compound nucleus, the corresponding pre-equilibrium radiation carries out fundamental information about the density distribution and the shape of the di-nuclear complex. It is worth noting that this mechanism may also provide a cooling effect, which could favor superheavy element formation [, ].
Apart from the strong influence of different parameters, such as mass and charge asymmetry, collision centrality and energy [, , ], collective oscillations which characterize the DD turn out to be mainly ruled by the isovector channel of the nuclear effective interaction, which yields once again the restoring force. However, within the selected beam energy (around 10 MeV/A), where the DD mechanism is better evidenced, other pre-equilibrium effects, such as nucleon and light particle emission, can occur, leading to a reduction of the initial charge asymmetry of the colliding nuclei and contributing to cooling down the system. Likewise the DD mechanism, also the N/Z ratio of the pre-equilibrium nucleon emission, has been proposed as a probe of the symmetry energy behavior below normal density [, , ].
In this article we review recent studies devoted to the investigation, within a semi-classical transport approach, of collective excitations in isolated nuclei and of pre-equilibrium effects, such as dipole radiation and nucleon emission, occurring in nuclear reactions at low beam energy [, ]. The nuclear effective interaction is described by Skyrme-like parameterizations, which are mainly tuned on the features of selected nuclei, especially in spin-isospin channels [, ]. We will explore the sensitivity of the mechanisms considered to specific properties of the effective interaction and, in the case of nuclear reactions, also to the strength of two-body (n-n) collision cross section. In particular, from our combined analysis, we aim to get a consistent picture of the impact of the density dependence of the symmetry energy on dipole excitations in neutron-rich systems and on the features of pre-equilibrium DD oscillations, together with nucleon emission. We stress the general interest of this study, considering the leading role played by the symmetry energy in nuclear structure problems (the neutron skin thickness, for instance) [–] and its impact in the astrophysical context [, ].
2. Theoretical Framework
We adopt here the same theoretical and numerical treatments illustrated in Zheng et al. [] and Zheng et al. [], namely calculations are based on the semi-classical Boltzmann-Nordheim-Vlasov (BNV) model [, ].
Within such a framework, the evolution of the system is investigated by solving the two dynamical coupled equations []:
where fq and ϵq, with q = n, p, are the distribution functions and the single particle energies of neutrons and protons, respectively. In the spirit of the density functional theory, the single particle energy, which includes the mean-field potential, can be derived from an energy density functional, []. The latter quantity, in the case of Skyrme-like interactions, is written as Raduta et al. []:
where (ρ = ρn + ρp, ρ3 = ρn − ρp) and (τ = τn + τp, τ3 = τn − τp) denote isoscalar and isovector density and kinetic energy densities, respectively, and the standard Skyrme parameters have been properly combined into the coefficients C.., D... In the calculations, the Coulomb contribution is also included []. The effect of the residual two-body correlations is taken into account in the collision integral, Icoll[fn, fp], employing the isospin, energy- and angular-dependent free nucleon-nucleon cross section. The test-particle (t.p.) method [] is adopted to integrate Equation (1). However, the finite number of t.p. considered requires setting a maximum cutoff of 50 mb for the n-n cross section [, ], to quench spurious collisions that may originate from an inaccurate evaluation of Pauli blocking effects.
The model illustrated here is able to describe nuclear dynamics at low beam energies, from fusion to quasi-fission and deep-inelastic processes [, , ]. Moreover, the features of zero-sound excitations are well reproduced, both in nuclear matter and finite nuclei [, , ], though quantum effects, such as shell effects, cannot be accounted for.
Among the different channels of the effective interaction, we are mainly interested in the isovector terms. Thus, we introduce the definition of the symmetry energy per nucleon, , where I = ρ3/ρ is the asymmetry parameter. The coefficient C(ρ) can be expressed in terms of the Skyrme coefficients:
where εF denotes the Fermi energy and m is the nucleon mass.
In our calculations we will employ the recently introduced SAMi-J Skyrme effective interactions. The details of the SAMi fitting protocol and the derivation of the corresponding parameters can be found in Roca-Maza et al. []. As a key feature, the SAMi-J family has been produced to allow for different values of the symmetry energy at normal density, J = C(ρ0), from 27 to 35 MeV, but keeping the same optimal values of the main isoscalar nuclear matter properties and of the main features of selected finite nuclei. In this way, these interactions mainly differ in the isovector channel and are thus well suited to explore the impact of isovector terms on a given observable.
It is worth noting here that, by construction, the Skyrme mean-field potential Uq is associated with a quadratic dependence on the momentum. This behavior is a good approximation for low momenta, such as in the situation explored in our study []. Actually, for the SAMi-J interactions, a rather flat momentum dependence is observed for the symmetry potential, according to the small splitting, , between neutron and proton effective masses. Moreover, an effective isoscalar mass m*(I = 0) = 0.67m MeV is predicted by these interactions.
In the following, we will consider three SAMi-J parametrizations: SAMi-J27, SAMi-J31 and SAMi-J35 []. Since the fitting procedure involves the properties of finite nuclei, the coefficient C(ρ) gets the same value, i.e., C(ρc) ≈ 22 MeV at the density ρc ≈ 0.6ρ0, that approximately represents the average density of nuclei of intermediate mass. Consequently, each parameterization is characterized by a different symmetry energy value, J, at normal density, as indicated in the corresponding interaction name.
The values of the slope parameter are reported in Table 1. The corresponding density dependence of C(ρ) is displayed in Figure 1A.
Table 1
| Effective interaction | J [MeV] | L [MeV] | Effective interaction | J [MeV] | L [MeV] |
|---|---|---|---|---|---|
| asy-soft | 30 | 14.8 | SAMi-J27 | 27 | 29.9 |
| asy-stiff | 30.5 | 79 | SAMi-J31 | 31 | 74.5 |
| asy-superstiff | 30.5 | 106 | SAMi-J35 | 35 | 115.2 |
The values of the symmetry energy J and its slope L at normal density are reported for the Skyrme interactions adopted in our study.
Figure 1
We will also adopt momentum-independent Skyrme interactions (Ceff = Deff = 0, m* = m), characterized by an incompressibility modulus K = 200 MeV [] and widely employed in the literature [, , ].
In order to distinguish these interactions from the SAMi-J family introduced above, which is momentum-dependent (MD), in the following we will indicate them as momentum-independent (MI) interactions. Concerning the symmetry energy, several trends are considered, as shown in Figure 1, lower, leading to different values of the slope L, but close values of the symmetry energy at normal density (J ≈ 30 MeV) (see also Table 1) []. As we will show in the following, the possibility to probe several interactions in the transport dynamics allows one to define the density regime explored in low-energy nuclear reactions and to test the impact of density dependent terms, such as the symmetry energy, on reaction observables.
The ground state of the considered nuclei is determined by solving Equation (1) in the stationary limit. Since we work with test particles which are usually associated with wave packets of finite width, some surface contributions are already implicitly taken into account, both in the initialization and in the dynamics, in addition to the surface terms of the SAMI-J interactions. In our case, in particular, we adopt triangular functions []. We find that the optimal reproduction of the experimental features (binding energy and charge radius) of selected nuclei is attained when suppressing the explicit surface terms, i.e., Csurf = Dsurf = 0. Therefore, we will adopt this choice in the following.
3. Dipole Excitations in Neutron-Rich Systems
For the study of collective motion in nuclei, we neglect the collision integral in Equation (1). Thus, we are lead to consider the Vlasov equation, which represents the semi-classical limit of Time-Dependent Hartree-Fock (TDHF) and, for small oscillations, of the Random Phase Approximation (RPA) equations. In our calculations, a number of 1, 500 t.p. per nucleon is considered, ensuring a good spanning of the phase space. We will consider the following neutron-rich nuclei, spanning three mass regions: 68Ni (N/Z = 1.43), 132Sn (N/Z = 1.64), 208Pb (N/Z = 1.54).
3.1. Ground State Properties
The numerical procedure that we adopt to define the ground state gives charge radius and binding energy values which agree rather well with the predictions of Hartree-Fock calculations [] and allow getting a reasonable reproduction of experimental values [], as one observes in Table 2 in the case of the MD parameterizations. MI calculations overestimate the binding energy if one imposes having similar neutron and proton density profiles as obtained in the MD case (see Table 2). The values reported for the neutron skin thickness are also in good agreement with previous results obtained with the Sly4 Skyrme interaction [], though the latter predicts smaller neutron and proton radii with respect to our results. On the other hand, a more diffuse neutron skin is observed in the case of the Relativistic Mean Field (RMF) calculations reported in Sarriguren et al. [] (see also [, ] for further details). Isoscalar and isovector density profiles are shown in Figure 2 for the system 132Sn and the three SAMi-J parameterizations adopted in our analysis. As evidenced in the left panel, a more diffuse density profile is obtained when increasing the slope parameter L. From the inspection of the isovector density (right panel), it appears that this effect can be ascribed to the development of a neutron skin. Indeed a larger slope L (see for instance the SAMi-J35 parametrization) is associated with a steeper variation of the symmetry energy around normal density, thus favoring the migration of the neutron excess toward the low-density nuclear surface. Similar behavior is seen for the 68Ni and 208Pb ground state configuration and also in the case of the MI interactions []. The trend observed for the dependence of the neutron skin thickness on the symmetry energy features is in agreement with previous investigations with other models [, –].
Table 2
| [fm] | [fm] | [fm] | [MeV] | |
|---|---|---|---|---|
| 68Ni | ||||
| asy-stiff | 4.104 | 3.907 | 0.197 | 10.905 |
| SAMi-J31 | 4.102 | 3.898 | 0.204 | 9.050 |
| Exp | — | 3.857 | — | 8.682 |
| 132Sn | ||||
| asy-stiff | 5.062 | 4.781 | 0.281 | 10.365 |
| SAMi-J31 | 5.035 | 4.741 | 0.294 | 8.552 |
| Exp | — | 4.709 | — | 8.354 |
| 208Pb | ||||
| asy-stiff | 5.793 | 5.592 | 0.201 | 9.826 |
| SAMi-J31 | 5.735 | 5.536 | 0.199 | 8.042 |
| Exp | — | 5.501 | — | 7.867 |
Neutron and proton root mean square radii, their difference, and binding energy for the three systems considered in our study, as obtained with asy-stiff (MI) and SAMi-J31 (MD) interaction.
The experimental values, for charge radius and binding energy, are also indicated [–].
Figure 2
3.2. Collective Dipole Response: Isoscalar-Isovector Mixing
We concentrate our analysis on the E1 (isoscalar and isovector) response of nuclear systems. Thus, we inject at the initial time the instantaneous excitation , at t = t0, along the z direction [, ], following the time evolution of the system until t = tmax. Here indicates the operator inducing dipole excitations of isoscalar or isovector type (k = S or V, respectively):
where τi = 0(1) for neutrons (protons) and 〈r2 〉 refers to the mean square radius of the system under study. It should be noticed that, in the general case of asymmetric systems (with different N and Z numbers), the operator also contains an isoscalar component.
The strength function Sk(E) is evaluated considering the Fourier transform of Dk(t), which is the expectation value of the time-dependent dipole moment:
where , with E = ℏω.
Introducing a gentle perturbation on the ground state of the considered nucleus, we follow the time oscillations of the dipole moment, solving Equation (1), until the final time tmax = 1, 800 fm/c. A filtering procedure, as described in Reinhard et al. [], was applied in order to cure the problems connected to the finite calculation time. To this purpose, a smooth cut-off function was introduced such that .
As discussed in Zheng et al. [], whereas in symmetric matter one can isolate pure isoscalar and isovector excitations, in asymmetric systems a mixing is generally observed, owing to the different amplitude of neutron and proton oscillations. It is quite interesting to try to get deeper insight into this effect and its dependence on the features of the effective interaction employed.
3.3. Sensitivity to System Size and Effective Interaction
Let us first discuss how the response of the system evolves in the three mass regions considered in this work. For the results shown in the following, we only consider the IS(IV) response generated by a corresponding IS(IV) perturbation.
Figure 3 represents dipole oscillations and corresponding strength, as a function of the excitation energy E, as obtained for the system 132Sn and the SAMi-J31 interaction. The (Figures 3A–D) correspond to an initial IS perturbation with ηS = 0.5 MeV fm−2, whereas an initial IV perturbation with ηV = 25 MeV has been considered in Figures 3E–H.
Figure 3
The mixing between isoscalar and isovector excitations is rather evident. Indeed, the IS perturbation (Figures 3A,B) also excites oscillations of the IV dipole moment (Figures 3C,D). In a similar way, when an IV perturbation is applied (Figures 3E,F), one also gets an isoscalar response (Figures 3G,H).
Let us start our discussion by looking at the features of the isovector response (Figure 3F). Here we easily identify the IV GDR peak, with EGDR ≈ 14 MeV. On the left, the low-energy region (the so called PDR region) is moderately populated, with some strength located between E1 = 9 MeV and E2 = 11 MeV. Interestingly, the contribution of the latter region is enhanced when looking at the IS projection (Figure 3H), where the corresponding strength now acquires a similar amplitude as compared to the GDR. This observation already suggests that these low-energy modes are mostly of isoscalar nature and is confirmed by the results obtained considering an initial IS perturbation (Figures 3A,B). Indeed, in Figure 3B one observes two important peaks with energies close to E1 and E2, together with a moderate strength contribution in the IV GDR region (EGDR ≈ 14 MeV). Considerable strength is also located in the high-energy region of the spectrum (E ≈ 29 MeV) and can be attributed to the IS GDR mode. One can notice that this mode also has some mixed character. In fact, a sizable (negative) contribution appears, at this energy, also in the IS projection corresponding to an initial IV perturbation (see Figure 3F).
From the results discussed above, one can conclude that in asymmetric systems the normal modes are of quite mixed nature, so that they can be excited, though with different strength, by both IS and IV perturbations. Thus it is appropriate to discuss essentially in terms of isoscalar-like (i.e., mostly isoscalar) and isovector-like (i.e., mostly isovector) modes. In particular, our analysis suggests that the modes located in the PDR region are isoscalar-like; they contribute to the IV response because of their mixed character []. The dependence of these effects on the features of the nuclear effective interaction is discussed in the next section.
In Figure 4 we show, for the SAMi-J31 effective interaction, the strength function corresponding to the IS (Figure 4, Left) and IV (Figure 4, Right) dipole response as a function of the excitation energy E. As a general remark, we observe that the response is shifted to lower energy regions when increasing the system mass. Let us start discussing the IS response, whose spectrum is generally characterized by two main region of contributions: a large peak in the strength, which is associated with the compressional IS GDR mode and is located at high energy, above 25 MeV for all the nuclei under study, and a quite fragmented response, which is observed in the low-energy domain, in all cases below 15 MeV. The isoscalar-like nature of the isolated high-energy mode is considered well established, while its spreading width is still under investigation, although a significant dependence on the size of the nucleus is already evidenced in Figure 4. Concerning the low-lying energy modes, despite of the fragmentation, one can generally observe the emergence of two main peaks of comparable height with respect to the strength of the IS GDR, in agreement with previous results deduced within other semi-classical studies [], where these excitations have been predominantly interpreted as surface modes. However, it is worth noting here that, owing to the coupling induced in neutron-rich systems discussed above, these oscillations are then responsible also for the strength observed in the PDR region of the isovector response. This correspondence holds for the three nuclei considered, including the largest system, 208Pb, where the low-lying IS peaks tend to merge together. The features regarding the low-energy part of the dipole spectrum can therefore be addressed by looking also at the IV response. In this case, however, we observe—for all nuclei—that the IV projection of the PDR is quite smaller than the IV GDR (about one order of magnitude), in agreement with previous RPA calculations []. We conclude that the PDR region is mainly populated by isoscalar low-energy modes, which generally involve mostly nucleons belonging to the nuclear surface []. Thus the position and the relative importance of the different low-lying energy modes may reflect the shape (i.e., the volume/surface relative contributions) of the density profile of the nucleus considered.
Figure 4
The reliability of our results is demonstrated by the good reproduction of the experimental data related to the IV GDR. Also the PDR region is reasonably reproduced, though a systematic overestimation is present in our calculations. This discrepancy might be attributable to the semi-classical treatment of surface effects. Indeed, these low-lying energy modes are mostly related to the oscillations of the most external nucleons. An improvement within the semi-classical framework can probably be achieved through a fine tuning of the coefficients Csur f and Dsur f in the Skyrme parameterizations.
To better explore this issue, in Figure 5 we compare the IV dipole response extracted within our semi-classical Vlasov model to the results of standard TDHF calculations []. Despite the general good agreement, especially for the heavier systems, of the main IV GDR peak energy resulting from semi-classical and quantal approaches, significant differences between the two calculations are observed for the low-lying dipole modes. This comparison supports the conclusion that Vlasov results, in the PDR region, are affected by our numerical treatment of surface effects and by the lack of gradient terms of intrinsic quantal nature. A deeper investigation of the detailed structure of these excitations, both in semi-classical and quantal approaches, must thus be envisaged. However, though the exact energy location of the PDR region is not well reproduced, it is still worth examining the dependence of the response in this region on the effective interaction employed.
Figure 5
We will focus on the description of 132Sn. Let us look in particular at the change introduced in the spectrum when employing different SAMi-J parameterizations. We offer a reminder that this allows us to appreciate the sensitivity to the isovector channel of the interaction. Qualitatively, looking at the central panel of Figure 6, it appears that the different peaks arising in the low-energy region of the IV dipole response become higher for larger L values. Moreover, the left panel indicates also that the IS response of the lowest energy mode increases with L. This is expected on the basis that a larger symmetry energy slope L leads to a larger coupling between isoscalar and isovector modes, as pointed out by calculations in asymmetric nuclear matter []. Moreover, as seen in Figure 2, stiffer symmetry energy leads to thicker neutron skin. Thus, surface modes become more important, with a sizable isovector component as well, owing to the neutron enrichment of the surface region. We conclude that the strength of the dipole response located in the PDR region is quite sensitive to the symmetry energy parameterization and, in particular, to its slope L. On the other hand, almost no sensitivity to the isovector channel is seen for the energy position of the PDR strength, as it is expected for IS-like excitations.
Figure 6
Other interesting features of the IV response can be discussed by taking into account also the results related to the MI Skyrme interactions, represented in the right panel of Figure 6, still for 132Sn. As far as the energy of the IV GDR is concerned, one can see that it does not evolve much in the SAMi-J case, whereas it shows a pronounced sensitivity to the interaction in the MI case, being smaller for the asy-superstiff parameterization. This suggests that the GDR energy reflects the value of the symmetry energy close to ρc = 0.6 ρ0, which indeed can be taken as the average density of medium-heavy nuclei. In fact, the three SAMi-J interactions have equal symmetry energy at ρc (see Figure 1A), whereas in the MI case (Figure 1B) the symmetry energy is smaller for the stiffer interaction. It is of particular interest to compare the results of the asy-stiff and SAMi-J31 parameterizations, which show a close density behavior of the symmetry energy (see Figure 1). In spite of this, one observes a higher frequency in the MD case. This can be ascribed to momentum-dependent effects, thus evidencing an interesting interplay between symmetry energy and other terms of the effective interaction in shaping the features of the nuclear response.
Finally, a quite pronounced IV peak is observed in the energy region above the GDR, whose strength looks sensitive to the stiffness of the interaction. As confirmed by the analysis of the transition densities, discussed below, this peak is associated with volume IV excitations.
3.4. Transition Densities
Additional information on the nature of the nuclear excitations, namely on the mixing of IS and IV components, is gained through the study of the associated transition densities. The latter describe how neutrons and protons move in response to the external perturbation, thus helping to identify the volume/surface character of the different modes []. The transition densities, δρq, essentially correspond to the density oscillations, around the ground state configuration, induced by the initial perturbation. They can be calculated separately for neutrons and protons. Exploiting the cylindrical symmetry of the system and making the same assumptions (linear response regime) as in Urban [], one can write: δρq(r, cos θ, t) = δρq(r, t)cos θ. Thus, at each time step, the transition densities can be finally extracted by performing an angular average, just as a function of the radial coordinate r.
As discussed above, different modes are excited by the delta function perturbation, Vext, associated with the operator . Thus, the observed transition densities will reflect the superposition of the different oscillations. As explained in Zheng et al. [], in order to extract the contribution of a given mode, of energy E, to the transition densities, one can consider the Fourier transform of δρq(r, t):
Owing to the finite calculation time, the sine function is multiplied by the same damping factor, as considered for the strength function Sk(E).
Since we are dealing with asymmetric systems, it is convenient to evaluate, for each mode, isoscalar and isovector transition densities (or, equivalently, neutron and proton transition densities). Then we expect the isoscalar component to be dominant in the case of isoscalar-like excitation, where neutrons and protons oscillate in phase, though with different amplitude. On the other hand, isovector-like oscillations should be associated with a dominant isovector component of the transition density.
In Figure 7, we display the isoscalar and isovector transition densities associated with the lowest energy peak observed in the IS response (the PDR peak), for 132Sn and the three SAMi-J parameterizations adopted. As shown in Figure 4, this excitation contributes also to the IV response. The curves in Figure 7 clearly manifest the isoscalar nature of the pygmy mode. Indeed, the amplitude of the isoscalar density fluctuation is predominant overall, in contrast with the isovector one, which is generally rather small. However, a significant isovector density oscillation seems to involve the external part, where its contribution equals the isoscalar one. In fact, as a consequence of the neutron skin development, in this radial region only the neutron practically oscillate, and this fluctuation is responsible for the observed IV projection in the PDR region. Moreover, the interactions characterized by a large slope L lead to an increase of both IS and IV transition densities (which practically coincide) in the surface region. Indeed, as discussed above, for increasing L, the system develops a thicker neutron skin, thus the surface region becomes rather neutron rich and isovector effects are correspondingly enhanced in this region. Hence, the analysis of the transition densities confirms the interpretation of the strength observed in the PDR region discussed above.
Figure 7
Figure 8 shows the results relative to other modes, which give a relevant contribution to the dipole strength, for the system 132Sn and the SAMi-J31 interaction. IS (IV) initial perturbations are considered in the left (right) panel. Let us comment first the features of the low-energy modes. As one can see in Figure 8, Left A, the transition densities of the second peak identified in the IS response (around E2 = 11 MeV, see Figure 6) indicate that neutron and proton densities are in phase, though some mismatch is present. The mode has an isoscalar-like character, but with some mixing, leading to a contribution also to the IV response.
Figure 8
Indeed, when the same energy region is explored through IV excitations (Figure 8, Right A), the difference between protons and neutrons becomes stronger. We note that the feature described above, namely the splitting of the PDR strength into an isoscalar contribution (at lower energy) and an isovector (more energetic) component has been reported also in recent theoretical and experimental analyses [–].
Turning now to discuss the highest energy isoscalar mode (Figure 8D) we observe transition densities of considerable amplitude also in the internal part of the system. This confirms that this mode can actually be associated with the IS dipole compression mode, of robust isoscalar nature. By comparing left and right panels, one can also notice that the features of this mode are practically not affected by the initial perturbation.
Concerning the modes that are essentially isovector-like (Figures 8B,C), one can see that also in this case the transition densities associated with IS or IV excitations are quite similar (compare left and right panels). The transition densities nicely indicate that, whereas the standard GDR (Figure 8B) is essentially a surface mode, the higher energy mode (Figure 8C) deeply involves the nucleons belonging to the internal part, exhibiting a double oscillation, typical of Steinwedel-Jensen volume modes.
4. Results for Reaction Dynamics
Let us move now to discuss reaction dynamics between charge asymmetric systems, where charge equilibration takes place []. Full BNV calculations are performed for the reaction 132Sn+58Ni at 10 MeV/A, considering several impact parameters leading to incomplete fusion and employing both MD and MI interactions. A proper number of t.p. (600 t.p./nucleon) is adopted to ensure a reasonable spanning of fq in phase space, as well as an acceptable computing time. Considering 132Sn as a projectile induces a sizable charge asymmetry in the entrance channel, also allowing to explore possible reaction effects related to the neutron skin thickness, whose dependence on L has been already addressed in the previous section [] and shown in Figure 2. Moreover, considering projectile (P) and target (T) with different N/Z ratios [ and in the case considered], neutron and proton centers of mass do not coincide, thus creating an initial dipole moment which may trigger DD oscillations along the rotating reaction symmetry axis. Another interesting aspect of nuclear reactions at energies just above the Coulomb barrier is that pre-equilibrium nucleon emission starts to take place. The two pre-equilibrium effects, namely nucleon and γ-ray emission, will be therefore addressed in the following [].
4.1. The DD Emission
We first discuss dipole oscillations following a collective bremsstrahlung analysis [, , ]. Similarly to Equation 5, the dipole moment is defined in coordinate space as:
where A = AT + AP is the total mass of the dinuclear system and Z = ZT + ZP (N = NT + NP) is the proton (neutron) number. Rn and Rp denote the center of mass of neutrons and protons, respectively. For the system considered, when the two nuclei touch each other, the dipole moment is equal to Di = 45.1 fm. The DD emission probability of photons with energy Eγ is given by (Eγ = ℏω):
where D″(ω) is the Fourier transform of the dipole acceleration D″(t) []. We will show the results for the average DD evolution, as obtained by considering 10 events for each specific calculation. This allows one to avoid spurious oscillations caused by the numerical noise associated with the finite number of test particles. We expect that, according to both mean-field and two-body collisional effects, the DD oscillations will be damped. Since microscopic calculations suggest that in-medium effects quench the two-body nucleon cross section [], we also performed simulations multiplying the latter quantity by a global factor fcs. Two cases (i.e., fcs = 0, corresponding to Vlasov calculations, and fcs = 0.5) are considered. Hence in the following we will explore, in addition to mean-field effects, also the influence of collisional damping and nucleon emission on the dipole oscillations.
In Figure 9, Left, we plot the time evolution of the DD, at b = 2 fm, as obtained for the SAMi-J31 EoS and with the different choices of fcs. The initial dipole moment is quite large because at the initial time considered, the distance between the centers of mass of the two nuclei is 14 fm. First of all, it can be interesting to compare the results related to the Vlasov case (fcs = 0) with the oscillations displayed in Figure 3. A considerable damping of the dipole oscillations is observed in Figure 9, owing to possible non-linear effects and to nucleon emission, which cools down the system and reduces the initial charge asymmetry. One can also notice that the DD oscillations with the free n-n cross section (fcs = 1) are damped faster than in the calculations associated with smaller fcs values. As a result, when neglecting the in-medium suppression of the n-n cross section, dipole oscillations are fully damped within about 600 fm/c. The corresponding power spectrum, |D″(ω)|2, which enters the expression of the photon emission probability [see Equation (9)] is also represented in Figure 9, Right. One can notice that the peak centroid is located at smaller energy with respect to the IV GDR in isolated nuclei. This is due to the elongated shape of the system at the initial reaction stage. Moreover, one observes that the centroids of the power spectra do not depend too much on the cross section choice. However, a slight shift to lower energies is observed especially in the fcs = 0 case, indicating that, in absence of two-body collisions, the systems maintain the elongated shape for a longer time. Dissipation effects are clearly larger for the larger cross section and the DD strength correspondingly decreases. The calculations for the other impact parameters indicate that the DD signal is quenched in more peripheral events, though similar features are observed with respect to the results for central collisions discussed above.
Figure 9
4.2. Sensitivity to the Effective Interaction
Once the effect of the two-body collisional damping is clarified, we can look at the role of different ingredients of the nuclear effective interaction, in determining the energy spectrum, analogous with the study of collective dipole modes carried out in section Sensitivity to System Size and Effective Interaction. The value of fcs = 1 will be fixed therefore in the following, while all interactions introduced in Figure 1 will be taken into account. The results are displayed in Figure 10, for MI (Figure 10, Left) and MD (Figure 10, Right) interactions, respectively.
Figure 10
Several differences emerge between the power spectra obtained in the two cases which are actually consistent with what we have already discussed in section Sensitivity to System Size and Effective Interaction. Indeed the peak centroid is insensitive to the effective interaction in the MD case, where a sizable dependence is observed in the MI case. As already stressed above, the restoring force of IV dipole oscillations is essentially ruled by the symmetry energy. Thus, we conclude that the DD features are also determined by the symmetry energy value in the density region around ρc ≈ 0.6 ρ0, the crossing point of the SAMi-J interactions. The lower frequency of the oscillations, with respect to the standard GDR, is due to the elongated shape of the system. In the MI case, the frequency of the power spectrum is higher for the asy-soft case, in connection to the larger value taken by the symmetry energy below normal density (see Figure 1B). One can also observe that a higher peak energy is also associated with a higher magnitude of the power spectrum, consistent with previous studies []. Moreover, the MD calculations are characterized by higher magnitude and frequency of the power spectrum, with respect to the MI results. As already stressed in the case of the small amplitude excitations (see Figure 6), this result is related to the influence of the effective mass on the features of collective dipole modes, for which it is well known that MD interactions yield a larger EWSR, which can better reproduce the experimental data []. As discussed in Zheng et al. [], the results illustrated above can be grasped, in a rather schematic manner, in terms of the damped oscillator model. Indicated by τ and ω0 damping time and oscillation frequency, respectively, the Fourier transform of the dipole acceleration can be written as Baran et al. []:
From the above equation, it is clear that the DD emission is governed by the size of the initial dipole Di, as expected, but it also reflects the amplitude of the oscillation frequency ω0. This is fully in line with the results of Figure 10, Left: a larger energy centroid also corresponds to a larger strength. One can also realize that a short damping rate, 1/τ, leading to strong two-body collisional effects, quenches the strength [see the denominator of Equation (10)], as is shown in Figure 9.
4.3. Nucleon Emission
Let us discuss here the nucleon emission characterizing the pre-equilibrium stage. As we will show in the following, nucleon and light particle emission may carry out important information on selected properties of the nuclear effective interaction. Indeed, for reactions at Fermi energies, the isotopic content of such emission has already been shown to reflect the behavior of the symmetry energy at subsaturation density. This is in fact the density region explored during the expansion phase of the nuclear composite system, when this emission mainly occurs [, , ]. The nucleons emitted escape from the dense compact system, so that they can be traced back just looking at the particles belonging to low-density regions (ρ < 0.01 fm−3) at the final calculation time, tmax = 600 fm/c. One observes that a larger nucleon emission is associated with MD interactions. This can be attributed to the fact that the most energetic particles feel a less attractive mean-field potential when momentum-dependent effects are turned on. Larger fcs values also lead to an enhanced pre-equilibrium emission, owing to the increased n-n collision number. On the contrary, the N/Z ratio decreases in the MD case and also in the calculations associated with a larger n-n cross section. Thus, it appears that, on top of the expected sensitivity to the symmetry energy, the N/Z ratio also reflects some isoscalar features, such as effective mass and n-n collisions. When the pre-equilibrium emission becomes more abundant, the relative importance of isospin effects may become smaller and hence the N/Z ratio approaches 1. This dependence on several aspects of the effective interation is seen also for the DD emission, as discussed above. However, one can pin down the sensitivity to the symmetry energy by inspecting in deeper detail the results corresponding to the three parameterizations considered in our study, either in the MI or MD case, with a given fcs choice. For instance, in the case of fcs = 1, the N/Z ratio obtained in central collisions evolves from 2.049 to 1.774 when going from the asy-soft to the asy-superstiff symmetry energy parameterizations, whereas it changes from 1.433 to 1.687 in correspondence to the SAMi-J27 and SAMi-J35 interactions, respectively. Bearing in mind that the N/Z ratio increases with the symmetry energy value (owing to the increased neutron repulsion), one can conclude that, in low-energy nuclear reactions, the pre-equilibrium emission mainly tests the density region in between ρc and ρ0. This statement is corroborated by the opposite trend of N/Z with respect to the slope, as obtained in the MI and MD cases. Indeed, the crossing point of the symmetry energy (see Figure 1) is different for the two types of interactions, so that the ordering of the symmetry energy strength with L, in the density region considered, is opposite in the two cases. Therefore, combining the study of DD and pre-equilibrium nucleon emission, one has the possibility to probe the effective interaction in the low-density range [ρc, ρ0].
4.4. Global Analysis
The results discussed above indicate that the selected pre-equilibrium mechanisms, i.e., γ radiation and nucleon emission, are influenced by several ingredients of the nuclear effective interaction. A deeper insight into this issue is obtained by performing a more global analysis in order to explore mutual correlations between the features of the interaction and some proper observables []. The correlation coefficient CXY between the variable X and observable Y is usually defined as:
where cov(X, Y) denotes the covariance, Ā and s(A) indicate average value and variance of A (=X or Y), respectively. These quantities are calculated from the considered set of MI and MD calculations, with different symmetry energy and n-n cross section parameterizations. It should be noticed that a linear correlation between X and Y leads to CXY = ±1, whereas, in the absence of correlations, one gets CXY = 0. Inspired by the results that have been presented so far, we select three observables, which can be investigated experimentally [, ]: the centroid (Ecentr) and the integral of the DD emission strength, and the N/Z ratio of the nucleons emitted. Three model parameters, which have been shown to impact significantly these observables, are considered: the symmetry energy slope L, the effective mass m* and the fcs value. The correlations between the model parameters and the observables are presented in Figure 11 (see the solid bars) []. In this analysis, we intend to use the L parameter to characterize the low-density behavior of the symmetry energy. Thus, the correlation functions have been evaluated excluding the SAMi-J27 and SAMi-J35 interactions. Indeed, within the SAMi-J family, the symmetry energy takes the same value below normal density (at ρ = ρc, see Figure 1), and this feature could blur the sensitivity to L. One can see that an appreciable sensitivity to the slope L can be identified for the DD centroid energy and for the N/Z of the pre-equilibrium emission. A negative correlation is observed (denoted by the blue color); indeed, both the centroid energy and the pre-equilibrium N/Z decrease for stiffer (i.e., with larger L) interactions, as discussed above. On the other hand, for the integral of the DD power spectrum, the sensitivity to L is overwhelmed by the huge dependence of the results on effective mass and n-n cross section. For the sake of completeness, calculations have been performed also considering the full set of MD and MI interaction (dashed bars in Figure 11). As one can see, this does not change the conclusions drawn above. This analysis underlines and better quantifies the concurrent impact of several aspects of the effective interaction on observables which, by construction, should be particularly sensitive to the isovector terms, namely to the symmetry energy. In particular, it appears that it would be rather difficult to pin down the sensitivity to the symmetry energy without fixing, at the same time, other ingredients, such as the effective mass and the strength of the residual interaction.
Figure 11
5. Conclusions
To summarize, in this work we have performed, within a semi-classical transport approach, a combined study of collective modes in neutron-rich nuclei and pre-equilibrium dipole radiation in low-energy nuclear reactions. Within the latter framework, pre-equilibrium nucleon emission has also been discussed. The aim of our investigation was indeed to assess the sensitivity of several observables involved in the reaction dynamics to the main ingredients of the nuclear effective interaction as well as to the in-medium n-n cross section.
Concerning the analysis of collective excitations in neutron rich-systems, we aimed at gaining better insight into the features of the low-lying IV dipole response which is experimentally observed in several neutron-rich nuclei []. Interesting features of the E1 nuclear response have been evidenced by exploring three mass regions and considering a variety of effective interactions, mainly differing in the isovector channel. Inspecting both IS and IV response of the systems investigated, our analysis indicates the emergence, in neutron rich systems, of a strong mixing between isoscalar and isovector excitations, analogous to the one discussed for infinite nuclear matter [], in agreement with previous semi-classical [] or RPA [] investigations. PDR excitations are mainly of isoscalar nature; however, because of mixing effects, some strength is observed also in the IV response, whose amplitude is rather sensitive to the slope L of the symmetry energy around saturation. This observation is associated with the appearance of a thicker neutron skin, in neutron-rich systems, for increasing L values. Hence our analysis confirms the important contribution of the study of low-lying dipole excitations to the symmetry energy debate. At last, it is also worth noting that our results give a centroid energy of the IV GDR that is quite close to the experimental value as well as to RPA calculations [], also reflecting the value of the symmetry energy below normal density.
Moreover, the results discussed here, in particular the link between the PDR strength, the neutron enrichment of the nuclear surface and the IS/IV mixing of the collective excitations, can be useful for the experimental search of the PDR and for a more accurate estimate of the corresponding strength. The latter, in turn, can provide information about the neutron skin thickness of the nucleus considered, complementary to more direct measurements.
As far as the study of pre-equilibrium dipole and nucleon emission is concerned, considering a variety of effective interactions in the calculations, we have asserted that these observables are also quite appropriate to explore the symmetry energy behavior below normal density, in the density range [0.6ρ0, ρ0]. However, though sensitive to the isovector channel of the interaction, we have brought out that these mechanisms strongly depend also on isoscalar terms, such as n-n correlations and momentum-dependent terms. A significant dependence on the latter terms, i.e., on the effective mass, is observed also for the collective excitations investigated in section 3 []. Moreover, the sensitivity to several ingredients of the effective interaction has been recently pointed out also for the competition between fusion and quasifission processes in heavy ion reactions close to the Coulomb barrier [].
Our analysis suggests that a consistent study of collective excitations and low-energy reaction dynamics may open the possibility, by a parallel investigation of different observables, to constrain at once the details of effective interaction and n-n cross section, together with the low-density behavior of the symmetry energy.
Statements
Author contributions
SB and HZ have been mostly involved in the semi-classical calculations leading to the results discussed in the manuscript. MC has supervised this work. SB, MC, and HZ have contributed to the writing of the manuscript and to the general interpretation of the results.
Acknowledgments
This project has received funding from the European Union's Horizon 2020 research and innovation program under grant agreement N. 654002.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1.
HarakehMNvan der WoudeA. Giant Resonances, Oxford Studies in Nuclear Physics, Vol. 24. New York, NY: Oxford University (2001).
2.
HartmannTEndersJMohrPVogtKVolzSZilgesA. Measurement of the dipole and electric quadrupole strength distributions up to 10 MeV in the doubly magic nuclei 40Ca and 48Ca. Phys Rev Lett. (2000) 85:274–7. 10.1103/PhysRevLett.85.274
3.
HartmannTEndersJMohrPVogtKVolzSZilgesA. Dipole and electric quadrupole excitations in 40, 48Ca. Phys Rev C. (2002) 65:034301. 10.1103/PhysRevC.65.034301
4.
KondoTUtsunomiyaHGorielySDaoutidisIIwamotoCAkimuneHet al. Total and partial photoneutron cross sections for Pb isotopes. Phys Rev C. (2012) 86:014316. 10.1103/PhysRevC.86.014316
5.
AdrichPKlimkiewiczAFallotMBoretzkyKAumannTCortina-GilDet al. Evidence for pygmy and giant dipole resonances in 130Sn and 132Sn. Phys Rev Lett. (2005) 95:132501. 10.1103/PhysRevLett.95.132501
6.
KlimkiewiczAPaarNAdrichPFallotMBoretzkyKAumannTet al. Nuclear symmetry energy and neutron skins derived from pygmy dipole resonances. Phys Rev C. (2007) 76:051603(R). 10.1103/PhysRevC.76.051603
7.
CarboneAColòGBraccoACaoL-GBortignonPFCameraFet al. Constraints on the symmetry energy and neutron skins from pygmy resonances in 68Ni and 132Sn. Phys Rev C. (2010) 81:041301(R). 10.1103/PhysRevC.81.041301
8.
WielandOBraccoA. The pygmy dipole resonance in 68Ni and the neutron skin. Prog Part Nucl Phys. (2011) 66:374–8. 10.1016/j.ppnp.2011.01.037
9.
TamiiAPoltoratskaIvonNeumann-Cosel PFujitaYAdachiTBertulaniCAet al. Complete electric dipole response and the neutron skin in 208Pb. Phys Rev Lett. (2011) 107:062502. 10.1103/PhysRevLett.107.062502
10.
RossiDMAdrichPAksouhFAlvarez-PolHAumannTBenlliureJet al. Measurement of the dipole polarizability of the unstable neutron-rich nucleus 68Ni. Phys Rev Lett. (2013) 111:242503. 10.1103/PhysRevLett.111.242503
11.
ToftHKLarsenACBürgerAGuttormsenMGörgenANyhusHTet al. Evolution of the pygmy dipole resonance in Sn isotopes. Phys Rev C. (2011) 83:044320. 10.1103/PhysRevC.83.044320
12.
AumannTNakamuraT. The electric dipole response of exotic nuclei. Phys Scr. (2013) T 152:014012. 10.1088/0031-8949/2013/T152/014012
13.
SavranDAumannTZilgesA. Experimental studies of the Pygmy Dipole Resonance. Prog Part Nucl Phys. (2013) 70:210–45. 10.1016/j.ppnp.2013.02.003
14.
PaarNVretenarDKhanEColòG. Exotic modes of excitation in atomic nuclei far from stability. Rep Prog Phys. (2007) 70:691–793. 10.1088/0034-4885/70/5/R02
15.
LanzaEGVitturiAAndrésMVCataraFGambacurtaD. Excitations of pygmy dipole resonances in exotic and stable nuclei via Coulomb and nuclear fields. Phys Rev C. (2011) 84:064602. 10.1103/PhysRevC.84.064602
16.
CrespiFCLBraccoANicoliniRMengoniDPellegriLLanzaEGet al. Isospin character of low-lying pygmy dipole states in 208Pb via inelastic scattering of 17O Ions. Phys Rev Lett. (2014) 113:012501. 10.1103/PhysRevLett.113.012501
17.
LanzaEGVitturiALitvinovaESavranD. Dipole excitations via isoscalar probes: The splitting of the pygmy dipole resonance in 124Sn. Phys Rev C. (2014) 89:041601. 10.1103/PhysRevC.89.041601
18.
EndresJSavranDButlerPAHarakehPAHarissopulosSHerzbergR-Det al. Structure of the pygmy dipole resonance in 124Sn. Phys Rev C. (2012) 85:064331. 10.1103/PhysRevC.85.064331
19.
RepkoAReinhardP-GNesterenkoVOKvasilJ. Toroidal nature of the low-energy E1 mode. Phys Rev C. (2013) 87:024305. 10.1103/PhysRevC.87.024305
20.
PapakonstantinouPHergertHRothR. Isoscalar and neutron modes in the E1 spectra of Ni isotopes and the relevance of shell effects and the continuum. Phys Rev C. (2015) 92:034311. 10.1103/PhysRevC.92.034311
21.
KnappFLo IudiceNVeselỳPAndreozziFDe GregorioGPorrinoA. Dipole response in 208Pb within a self-consistent multiphonon approach. Phys Rev C. (2015) 92:054315. 10.1103/PhysRevC.92.054315
22.
KnappFLo IudiceNVeselỳPAndreozziFDe GregorioGPorrinoA. Dipole response in 132Sn within a self-consistent multiphonon approach. Phys Rev C. (2014) 90:014310. 10.1103/PhysRevC.90.014310
23.
AuerbachNStoyanovCAndersMRShlomoS. Isoscalar and isovector dipole strength distributions in nuclei and the Schiff moment. Phys Rev C. (2014) 89:014335. 10.1103/PhysRevC.89.014335
24.
BaranVColonnaMGrecoVDi ToroM. Reaction dynamics with exotic nuclei. Phys Rep. (2005) 410:335–466. 10.1016/j.physrep.2004.12.004
25.
LiB-AChenL-WKoCM. Recent progress and new challenges in isospin physics with heavy-ion reactions. Phys Rep. (2008) 464:113–281. 10.1016/j.physrep.2008.04.005
26.
GiulianiGZhengHBonaseraA. The many facets of the (non-relativistic) Nuclear Equation of State. Prog Part Nucl Phys. (2014) 76:116–64. 10.1016/j.ppnp.2014.01.003
27.
BaranVBrinkDMColonnaMDi ToroM. Collective Dipole Bremsstrahlung in Fusion Reactions. Phys Rev Lett. (2001) 87:182501. 10.1103/PhysRevLett.87.182501
28.
SimenelCChomazPhde FranceG. Quantum Calculation of the Dipole Excitation in Fusion Reactions. Phys Rev Lett. (2001) 86:2971–4. 10.1103/PhysRevLett.86.2971
29.
UmarASOberackerVE. 64Ni+132Sn fusion within the density-constrained time-dependent Hartree-Fock formalism. Phys Rev C. (2007) 76:014614. 10.1103/PhysRevC.76.014614
30.
WuHLTianWDMaYGCaiXZChenJGFangDQet al. Dynamical dipole γ radiation in heavy-ion collisions on the basis of a quantum molecular dynamics model. Phys Rev C. (2010) 81:047602. 10.1103/PhysRevC.81.047602
31.
FlibotteSChomazPhColonnaMCromazMDeGraafJDrakeTEet al. Pre-equilibrium effects in the population of giant dipole resonances. Phys Rev Lett. (1996) 77:1448–51. 10.1103/PhysRevLett.77.1448
32.
PapaMTianWGiulianiGAmoriniFCardellaGDi PietroAet al. Pre-equilibrium γ-ray emission induced in the 40Ca+48Ca system at 10 MeV/nucleon and isospin equilibration processes. Phys Rev C. (2005) 72:064608. 10.1103/PhysRevC.72.064608
33.
PierroutsakouDMartinBAgodiCAlbaRBaranVBoianoAet al. Dynamical dipole mode in fusion reactions at 16 MeV/nucleon and beam energy dependence. Phys Rev C. (2009) 80:024612. 10.1103/PhysRevC.80.024612
34.
GiazACorsiABarliniSKravchukVLWielandOColonnaMet al. Measurement of dynamical dipole γ-ray emission in the N/Z-asymmetric fusion reaction 16O+116Sn at 12 MeV/nucleon. Phys Rev C. (2014) 90:014609. 10.1103/PhysRevC.90.014609
35.
SimenelCChomazPhde FranceG. Fusion process studied with a preequilibrium giant dipole resonance in time-dependent Hartree-Fock theory. Phys Rev C. (2007) 76:024609. 10.1103/PhysRevC.76.024609
36.
ParascandoloCPierroutsakouDAlbaRDel ZoppoAMaiolinoCSantonocitoDet al. Evidence of dynamical dipole excitation in the fusion-evaporation of the 40Ca+152Sm heavy system. Phys Rev C. (2016) 93: 044619. 10.1103/PhysRevC.93.044619
37.
BaranVCabibboMColonnaMDi ToroMTsonevaN. The dynamical dipole mode in dissipative heavy-ion collisions. Nucl Phys A. (2001) 679:373–92. 10.1016/S0375-9474(00)00365-1
38.
BaranVRizzoCColonnaMDi ToroMPierroutsakouD. Dynamical dipole mode in fusion reactions with exotic nuclear beams. Phys Rev C. (2009) 79:021603(R). 10.1103/PhysRevC.79.021603
39.
ZhengHBurrelloSColonnaMBaranV. Dipole response in neutron-rich nuclei with new Skyrme interactions. Phys Rev C. (2016) 94:014313. 10.1103/PhysRevC.94.014313
40.
ZhengHBurrelloSColonnaMBaranV. Pre-equilibrium effects in charge-asymmetric low-energy reactions. Phys Lett B. (2017) 769:424–9. 10.1016/j.physletb.2017.04.002
41.
Roca-MazaXColòGSagawaH. New Skyrme interaction with improved spin-isospin properties. Phys Rev C. (2012) 86:031306(R).
42.
Roca-MazaXBrennaMAgrawalBKBortignonPFColòGLi-Gang CaoNPet al. Giant quadrupole resonances in 208Pb, the nuclear symmetry energy, and the neutron skin thickness. Phys Rev C. (2013) 87:034301. 10.1103/PhysRevC.87.034301
43.
PiekarewiczJAgrawalBKColòGNazarewiczWPaarNReinhardP-Get al. Electric dipole polarizability and the neutron skin. Phys Rev C. (2012) 85:041302(R). 10.1103/PhysRevC.85.041302
44.
GaidarovMKAntonovANSarrigurenPMoya de GuerraE. Surface properties of neutron-rich exotic nuclei: a source for studying the nuclear symmetry energy. Phys Rev C. (2011) 84:034316. 10.1103/PhysRevC.84.034316
45.
GaidarovMKAntonovANSarrigurenPMoya de GuerraE. Symmetry energy of deformed neutron-rich nuclei. Phys Rev C. (2012) 85:064319. 10.1103/PhysRevC.85.064319
46.
SteinerAWPrakashMLattimerJMEllisPJ. Isospin asymmetry in nuclei and neutron stars. Phys Rep. (2005) 411:325–75. 10.1016/j.physrep.2005.02.004
47.
BurrelloSGulminelliFAymardFColonnaMRadutaAdR. Heat capacity of the neutron star inner crust within an extended nuclear statistical equilibrium model. Phys Rev C. (2015) 92:055804. 10.1103/PhysRevC.92.055804
48.
BonaseraAGulminelliFMolitorisJ. The Boltzmann equation at the borderline. A decade of Monte Carlo simulations of a quantum kinetic equation. Phys Rep. (1994) 243:1–124. 10.1016/0370-1573(94)90108-2
49.
GuarneraAColonnaMChomazPh. 3D stochastic mean-field simulations of the spinodal fragmentation of dilute nuclei. Phys Lett B. (1996) 373:267–74. 10.1016/0370-2693(96)00152-9
50.
ColonnaMDi ToroMLarionovAB. Collective modes in asymmetric nuclear matter. Phys Lett B. (1998) 428:1–7. 10.1016/S0370-2693(98)00383-9
51.
RadutaAdRAymardFGulminelliF. Clusterized nuclear matter in the (proto-)neutron star crust andthe symmetry energy. Eur Phys J A. (2014) 50:24. 10.1140/epja/i2014-14024-y
52.
WongCY. Dynamics of nuclear fluid. VIII. Time-dependent Hartree-Fock approximation from a classical point of view. Phys Rev C. (1982) 25:1460–75. 10.1103/PhysRevC.25.1460
53.
BaranVColonnaMDi ToroMGrecoVZielinska-PfabécMWolterHH. Isospin effects in nuclear fragmentation. Nucl Phys A. (2002) 703:603–32. 10.1016/S0375-9474(01)01671-2
54.
RizzoJColonnaMOnoA. Comparison of multifragmentation dynamical models. Phys Rev C. (2007) 76:024611. 10.1103/PhysRevC.76.024611
55.
RizzoCBaranVColonnaMCorsiADi ToroM. Symmetry energy effects on fusion cross sections. Phys Rev C. (2011) 83:014604. 10.1103/PhysRevC.83.014604
56.
RizzoCColonnaMBaranVDi ToroM. Exotic break-up modes in heavy ion reactions at low energies. Phys Rev C. (2014) 90:054618. 10.1103/PhysRevC.90.054618
57.
UrbanM. Pygmy resonance and torus mode within Vlasov dynamics. Phys Rev C. (2012) 85:034322. 10.1103/PhysRevC.85.034322
58.
LiX-HGuoW-JLiB-AChenL-WFattoyevFJNewtonWG. Neutron–proton effective mass splitting in neutron-rich matter at normal density from analyzing nucleon–nucleus scattering data within an isospin dependent optical model. Phys Lett B. (2015) 743:408–14. 10.1016/j.physletb.2015.03.005
59.
BaranVColonnaMDi ToroMCroitoruADumitruD. Connecting the pygmy dipole resonance to the neutron skin. Phys Rev C. (2013) 88:044610. 10.1103/PhysRevC.88.044610
60.
De VriesHDe JagerCWDe VriesC. Nuclear charge-density-distribution scattering. Atom Data Nucl Data Tables. (1987) 36:495–536. 10.1016/0092-640X(87)90013-1
61.
AudiGWapstraAH. The 1995 update to the atomic mass evaluation. Nucl Phys A. (1995) 595:409–80.
62.
AngeliIMarinovaKP. Table of experimental nuclear ground state charge radii: an update. Atom Data Nucl Data Tables (2013) 99:69–95. 10.1016/j.adt.2011.12.006
63.
SarrigurenPGaidarovMKde GuerraEMAntonovAN. Nuclear skin emergence in Skyrme deformed Hartree-Fock calculations. Phys Rev C. (2007) 76:044322. 10.1103/PhysRevC.76.044322
64.
CalvayracFReinhardPGSuraudE. Spectral signals from electronic dynamics in sodium clusters. Ann Phys. (1997) 225:125–62. 10.1006/aphy.1996.5654
65.
BaranVFrecusBColonnaMDi ToroM. Pygmy dipole resonance: Collective features and symmetry energy effects. Phys Rev C. (2012) 85:051601. 10.1103/PhysRevC.85.051601
66.
ReinhardP-GStevensonPDAlmehedDMaruhnJAStrayerMR. Role of boundary conditions in dynamic studies of nuclear giant resonances and collisions. Phys Rev E. (2006) 73:036709. 10.1103/PhysRevE.73.036709
67.
Roca-MazaXPozziGBrennaMMizuyamaKColòG. Low-lying dipole response: isospin character and collectivity in 68Ni, 132Sn, and 208Pb. Phys Rev C. (2012) 85:024601. 10.1103/PhysRevC.85.024601
68.
LacroixDAyikSChomazP. Nuclear collective vibrations in extended mean-field theory. Prog Part Nucl Phys. (2004) 52:497–563. 10.1016/j.ppnp.2004.02.002
69.
LiGQMachleidtR. Microscopic calculation of in-medium nucleon-nucleon cross sections. Phys Rev C. (1993) 48:1702–12. 10.1103/PhysRevC.48.1702
70.
ZhangYXTsangMBLiZX. Covariance analysis of symmetry energy observables from heavy ion collision. Phys Lett B. (2015) 749:262–6. 10.1016/j.physletb.2015.07.064
71.
KongH-YXuJChenL-WLiB-AMaY-G. Constraining simultaneously nuclear symmetry energy and neutron-proton effective mass splitting with nucleus giant resonances using a dynamical approach. Phys Rev C. (2017) 95:034324. 10.1103/PhysRevC.95.034324
72.
ZhengHBurrelloSColonnaMLacroixDScampsG. Connecting the nuclear equation of state to the interplay between fusion and quasifission processes in low-energy nuclear reactions. Phys Rev C. (2018) 98:024622. 10.1103/PhysRevC.98.024622
Summary
Keywords
symmetry energy, EoS, collective motion, low-energy reaction dynamics, charge equilibration
Citation
Burrello S, Colonna M and Zheng H (2019) The Symmetry Energy of the Nuclear EoS: A Study of Collective Motion and Low-Energy Reaction Dynamics in Semiclassical Approaches. Front. Phys. 7:53. doi: 10.3389/fphy.2019.00053
Received
21 December 2018
Accepted
20 March 2019
Published
17 April 2019
Volume
7 - 2019
Edited by
Angela Bonaccorso, National Institute of Nuclear Physics, Section of Pisa, Italy
Reviewed by
Theodoros Gaitanos, Aristotle University of Thessaloniki, Greece; Mitko Gaidarov, Institute for Nuclear Research and Nuclear Energy (BAS), Bulgaria
Updates
Copyright
© 2019 Burrello, Colonna and Zheng.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Maria Colonna colonna@lns.infn.it
This article was submitted to Nuclear Physics, a section of the journal Frontiers in Physics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.