Using the quantum molecular dynamics model, we study the nuclear dynamics at the balance energy of mass asymmetric colliding nuclei by keeping the total mass of the system fixed and by varying the mass asymmetry (η=(AT-AP)/(AT+AP), where AT and AP are the masses of the target and projectile, resp.) of the reaction. In particular, we study the various quantities like average and maximum density and temperature, collision rate, participant-spectator matter, anisotropic ratio, relative momentum, and their mass asymmetry and mass dependence. Our results indicate sizeable effect of mass asymmetry on these quantities.
1. Introduction
It is now well established that collective transverse flow of the nucleons is a signature of the interaction and can provide information about the equation of state (EoS) as well as nucleon-nucleon (nn) cross-section of the nuclear matter. Extensive studies have been done over the past three decades on the sensitivity of collective transverse flow towards the nuclear EoS, nn cross-section, and entrance channel parameters such as incident energy of projectile, size of the system (ATOT=AT+AP), and colliding geometry (i.e., impact parameter) [1–9]. The disappearance of flow is predicted to appear at some incident energy, which is termed as balance energy (Ebal) [10]. This balance energy has been subjected to extensive theoretical and experimental calculations to know its accurate value as well as its mass and impact parameter dependence [11–22].
Recently, the sensitivity of collective transverse flow and Ebal towards the mass asymmetry of the reaction has been carried out at different colliding geometries by keeping the total mass of the system fixed [20–22]. It has been found that almost independent of total system mass as well as colliding geometry, mass asymmetry has a uniform effect on the collective transverse flow and its disappearance [20]. This study was motivated from the recent observations by FOPI collaboration on the mass asymmetric reactions of Ca40 + Au197 and Au197 + Ca40 [23, 24]. They observe that flow in asymmetric reactions is a key observation for investigating the reaction dynamics. The difference between the reaction dynamics for symmetric and asymmetric reactions is attributed to the different role played by excitation energy in these reactions. In symmetric reactions most of the excitation energy is deposited in the form of compressional energy, whereas an asymmetric reaction deposits it in the form of thermal energy [25]. Therefore, in the present paper, this study is further extended for the central collisions to see the effect of mass asymmetry of the reaction on the nuclear dynamics at the balance energy. Similar work has been carried out by Sood and Puri [26] for the nearly symmetric and symmetric reactions at the balance energy, but the role of mass asymmetry in the participant-spectator matter, average and maximum density and temperature, net collisions, anisotropic ratio, relative momentum, and the mass dependence of these quantities is not taken care of. This has been taken care of in the present study. The study is made within the framework of quantum molecular dynamics (QMD) model [12–18, 20–22, 25–44], which is explained in Section 2. Results and discussion are explained in Section 3 and finally we summarize the results in Section 4.
2. The Model
The quantum molecular dynamics model [12–18, 20–22, 25–44] simulates the reaction on an event by event basis. Here each nucleon i is represented by a Gaussian wave packet with a width of L centered around the mean position r→i(t) and mean momentum p→i(t). Here each nucleon is represented by a coherent state of the form
(1)ϕi(r→,p→,t)=1(2πL)3/4e[-{r→-r→i(t)}2/4L]e[ip→i(t)·r→/ħ].
The Wigner distribution of a system with AT+AP nucleons is given by
(2)f(r→,p→,t)=∑i=1AT+AP1(πħ)3e[-{r→-r→i(t)}2/2L]e[-{p→-p→i(t)}22L/ħ2],
with L=1.08fm2.
The center of each Gaussian (in the coordinate and momentum space) is chosen by the Monte Carlo procedure. The momentum of nucleons (in each nucleus) is chosen between zero and local Fermi momentum [=2miVi(r→); Vi(r→) is the potential energy of nucleon i]. Naturally, one has to take care that the nuclei, thus generated, have right binding energy and proper root mean square radii.
The centroid of each wave packet is propagated using the classical equations of motion:
(3)dr→idt=dHdp→i,dp→idt=-dHdr→i,
where the Hamiltonian is given by
(4)H=∑ip→i22mi+Vtot.
Our total interaction potential Vtot reads as
(5)Vtot=VLoc+VYuk+VCoul+VMDI
with
(6)VLoc=t1δ(r→i-r→j)+t2δ(r→i-r→j)δ(r→i-r→k),VYuk=t3e-|r→i-r→j|/μ(|r→i-r→j|μ),
with μ=1.5 fm and t3=-6.66 MeV.
The static (local) Skyrme interaction can further be parameterized as
(7)VLoc=α(ρρo)+β(ρρo)γ.
Here α, β, and γ are the parameters that define equation of state. The momentum-dependent interaction is obtained by parameterizing the momentum dependence of the real part of the optical potential. The final form of the potential reads as
(8)VMDI≈t4ln2[t5(p→i-p→j)2+1]δ(r→i-r→j).
Here t4 = 1.57 MeV and t5=5×10-4MeV-2. A parameterized form of the local plus momentum-dependent interaction (MDI) potential (at zero temperature) is given by
(9)V=α(ρρ0)+β(ρρ0)γ+δln2[ϵ(ρ/ρ0)2/3+1]ρρ0.
The parameters α, β, and γ in above equation must be readjusted in the presence of momentum-dependent interactions so as to reproduce the ground state properties of the nuclear matter. The set of parameters corresponding to different equations of state can be found in [37].
3. Results and Discussion
For the present analysis, we simulated the reactions of O817+Na1123 (η=0.15), N714+Mg1226 (η=0.3), B510+Si1430 (η=0.5), and Li36+S1634 (η=0.7) for ATOT = 40, Ar1836+Ca2044 (η=0.1), Si1428+Cr2452 (η=0.3), Ne1020+Ni2860 (η=0.5), and B510+Ge3270 (η=0.7) for ATOT = 80, Ge3270+Zr4090 (η=0.125), Fe2654+Cd48106 (η=0.325), Ca2040+Te52120 (η=0.5), and Mg1224+Ce58136 (η=0.7) for ATOT = 160, and Cd48108+Ba56132 (η=0.1), Sr3884+Dy66156 (η=0.3), Ni2860+W74180 (η=0.5), and Ar1836+Pb82204 (η=0.7) for ATOT = 240, at their corresponding theoretical balance energies (taken from [20]). The balance energies at which these reactions are simulated were calculated using a momentum-dependent soft equation of state with standard energy-dependent cugnon cross-section and reduced impact parameter (b^ = b/bmax, where bmax=R1+R2; Ri is the radius of projectile or target) of 0.25. The values of parameters used in (9) for momentum-dependent soft equation of state are compressibility (K(MeV))=200, α(MeV)=-390, β(MeV)=320, γ=1.14, δ(MeV)=1.57, and ϵ=21.54. The reactions are followed uniformly up to 500 fm/c.
In Figure 1, we display the 〈ρavg〉/ρ0 (left column) and 〈ρmax〉/ρ0 (right column) as a function of the reaction time. The values of 〈ρavg〉 and 〈ρmax〉 are calculated within a sphere of radius 2 fm around the center of mass. The density is then computed at each time step during the reaction using
(10)ρ(r→,t)=∑i=1AT+AP1(2πL)3/2e(-(r→-r→i(t))2/2L).
The results are displayed by varying η from 0.1 to 0.7 for different mass ranges; that is, ATOT = 40, 80, 160, and 240. It is clear from the figure that the maximal value of 〈ρavg〉/ρ0 decreases with increase in η. A similar trend can be seen for the evolution of 〈ρmax〉/ρ0 at all mass ranges. This is due to the decrease in the interaction region with increase in mass asymmetry of the reaction. Further, one should note that for nearly symmetric reactions with η = 0.1, the values of 〈ρavg〉/ρ0 and 〈ρmax〉/ρ0 are comparable for the heavy systems, indicating the wide and uniform formation of dense matter in the central region of 2 fm radius. This is similar to what predicted in [26]. However, a nonhomogeneous nature of the dense matter is seen for lighter systems at all asymmetries. For each mass range, as η increases and for each η, as ATOT decreases, the incident energy increases. Therefore, for each ATOT, the reactions with larger asymmetries and for each η, the systems with smaller ATOT, finishes much earlier. Similarly, the peak values of the densities are also delayed for reactions of heavy colliding nuclei and smaller asymmetries. Also, one should note that the time of saturation increases with increase in η.
The time evolution of the average density 〈ρavg〉/ρ0 (a) and maximum density 〈ρmax〉/ρ0 (b) reached in a central sphere of radius 2 fm for different system masses. The results for different asymmetries η = 0.1, 0.3, 0.5, and 0.7 are represented, respectively, by the solid, dashes, dotted, and dashed-dotted lines.
Similar to Figure 1, we display the net collision rate as a function of the reaction time in Figure 2. It is directly connected with the density. Due to the decrease in interaction volume with increase in mass asymmetry of colliding nuclei, the net collision rate decreases with increase in η and the interactions among nucleons also cease earlier. This is observed for all system mass ranges. As expected, with increase in system mass, opposite trend is observed for each η. This behavior is also evident from the density profile (see Figure 1).
Same as Figure 1, but for the rate of allowed collisions versus reaction time.
It is well known that balance energy represents a counterbalancing between the attractive and repulsive forces; therefore, this fact should also be reflected in quantities like spectator and participant matter. In Figure 3, we display the normalized spectator (a, b) and participant (c, d) matter as a function of reaction time. The results are displayed for different asymmetries by keeping the total mass of the system fixed as ATOT = 80 and 240. All nucleons having experienced at least one collision are termed as participant and the remaining matter is termed as spectator matter. One can also define spectator and participant matter in terms of rapidity distribution; however, the results are similar for both definitions, as shown in [26]. From the figure, one clearly sees that at the start of the reaction, all nucleons constitute spectator matter; therefore, no participant matter exists at initial time (i.e., t=0 fm/c). Due to decrease in the number of nn collisions with increase in η, the spectator (participant) matter increases (decreases) with increase in η. Similar behavior is seen for different system masses.
The time evolution of the normalized spectator matter (a, b) and participant matter (c, d). The results for total system mass ATOT = 80 and 240 are displayed in left and right columns. Lines have the same meaning as in Figure 1.
In Figure 4, we display the time evolution of relative momentum 〈KR〉 (a) and anisotropy ratio 〈Ra〉 (b) for different η and ATOT. The anisotropy ratio 〈Ra〉 is defined as [38–44]
(11)〈Ra〉=〈px2〉+〈py2〉2〈pz2〉.
This anisotropy ratio 〈Ra〉 is an indicator of the global equilibrium of the system as it does not depend on the local density. The full global equilibrium averaged over large number of events will correspond to 〈Ra〉 = 1. The second quantity, the relative momentum 〈KR〉 of two colliding spheres, is defined as [38–44]
(12)〈KR〉=〈|P→P(r→,t)-P→T(r→,t)|ħ〉,
where
(13)P→i(r→,t)=∑j=1AP→j(t)ρj(r→,t)ρj(r→,t),i=1,2.
Here P→j and ρj are the momentum and density experienced by the jth particle and i stands either for target or projectile. As noted, this quantity measures deviation from a single Fermi sphere and hence represents local equilibrium. Such a concept of local equilibrium is commonly used in the hydrodynamical model. Obviously, with the passage of the time, density in a central sphere will decrease due to lesser and lesser nucleons and, as a result, 〈KR〉 will also decrease. On the other hand, no such density dependence exists for 〈Ra〉. The anisotropy ratio 〈Ra〉 will saturate after the finishing of the reaction. From the figure, we see that 〈KR〉 decreases as the reaction proceeds, while 〈Ra〉 ratio increases and saturates after the high dense phase is over. Due to increase in incident energy with η, thermalization is little bit better achieved for larger η. Similar behavior is observed for each fixed system mass. However, for each η, due to high density obtained in heavier systems, the thermalization is better achieved compared to lighter colliding nuclei.
Same as Figure 1, but for the relative momentum 〈KR〉 (a) and anisotropy ratio 〈Ra〉 (b).
Temperature is one of the associated quantities linked with a dense matter. In principle, one can define true temperature only for a thermalized and equilibrated matter. Since in heavy-ion collisions the matter is nonequilibrated, one cannot talk of temperature. One can only look in terms of the local environment. In our case, we follow the description of the temperature as given in [45–47]. The extraction of the temperature T is based on the local density approximation; that is, one deduces the temperature in a volume element surrounding the position of each particle at a given time step [45–47]. In present case each local volume element of nuclear matter in coordinate space and time has some temperature defined by the diffused edge of the deformed Fermi distribution consisting of two colliding Fermi spheres, which is typical for a nonequilibrium momentum distribution in heavy-ion collisions. In this formalism which is dubbed the hot Thomas-Fermi approach [45–47], one determines extensive quantities such as the density and kinetic energy as well as entropy with the help of momentum distributions at a given temperature. Using this formalism, we also extracted the average and maximum temperature within a central sphere of 2 fm radius as described in the case of density.
In Figure 5, we display the mass dependence of the maximal values of 〈Tavg〉, 〈Tmax〉, 〈ρavg〉, and 〈ρmax〉 as well as the final stage value of (allowed) nn collisions and spectator/participant matter. For the nearly symmetric reaction having η = 0.1, the dependence is similar to that shown in [26]. All the quantities follow a power law behavior ∝ATOTτ. The values of power factor τ are displayed in the figure. It shows that the dependence of temperature reached in the central region is in sharp contradiction to the evolution of density reached in the central region because the dependence of maximal value of 〈ρavg〉 and 〈ρmax〉 is very weak on system mass for all η while weak dependence on system mass is seen in case of temperature. This is due to the fact that the temperature depends also on the kinetic energy (i.e., excitation energy) of the system [45–47]. For a given colliding geometry, the maximal value of 〈Tavg〉 and 〈Tmax〉 also depends on the bombarding energy rather than on the size of the interacting source [48]. Since, in our case, we have carried out the study at balance energy, which varies with the system mass, this in turn leads to dependence of temperature on the system mass due to the variation in the bombarding energy. At a fixed bombarding energy, dependence could be different. For nn collisions, one sees a linear enhancement with the system mass for each η. Naturally, at a fixed energy, the nn collisions should scale as ATOT [26] and it is seen for each η. The spectator and participant matter are nearly independent of the mass of the system for each η.
The maximal value of the average temperature 〈Tavg〉max (a), maximum temperature 〈Tmax〉max (b), average density 〈ρavg〉max (c), maximum density 〈ρmax〉max (d), total number of the allowed collisions obtained at the final stage (e), and final saturated spectator/participant matter (f) as a function of total mass of the system. The results for different asymmetries η = 0.1, 0.3, 0.5, and 0.7 are represented, respectively, by the open squares, circles, triangles, and inverted triangles. The lines are power law (∝ATOTτ) fits to the calculated results. The values of the power factor τ are displayed in the figure for various quantities.
In Figure 6, we display the mass asymmetry dependence of the maximal values of 〈Tavg〉, 〈Tmax〉, 〈ρavg〉, and 〈ρmax〉 as well as the final stage value of (allowed) nn collisions and spectator/participant matter. The results are displayed for different mass ranges. Lines are the linear fits (∝mη). We find a clear dependence of quantities on η. The dependence of maximal value of 〈ρavg〉 and 〈ρmax〉 on η decreases with increase in system mass while in case of 〈Tavg〉max and 〈Tmax〉max; different trend is seen which is due to the same reason as explained in earlier paragraph, while for nn collisions an opposite trend is observed. The nn collisions also decrease with increase in η, but due to larger interaction volume in heavier systems, heavier colliding nuclei show large η dependence. Obviously, the spectator matter and the participant matter will behave oppositely. For each system mass, the spectator (participant) matter increases (decreases) with increase in η and similar to density; the η dependence decreases with increase in total mass of the system. From this figure, it is clear that mass asymmetry of the reaction has a significant role in the nuclear dynamics of the reaction. Therefore, while studying various phenomenon in the intermediate energy heavy-ion collisions, one should take care of the mass asymmetry of the reaction.
Same as Figure 5, but as a function of mass asymmetry of the reaction. The results for different system masses ATOT = 40, 80, 160, and 240 are represented, respectively, by the solid squares, circles, triangles, and inverted triangles. Lines are the linear fits (∝mη).
4. Summary
We studied the nuclear dynamics (particularly, average and maximum temperature and density, collision rate, participant-spectator matter, anisotropic ratio, relative momentum, and their mass asymmetry and mass dependence) at the balance energy of mass asymmetric reactions by keeping the total mass of the system fixed as 40, 80, 160, and 240 and varying the mass asymmetry of the reaction from 0.1 to 0.7. A sizeable effect of mass asymmetry on these quantities is observed. Therefore one cannot ignore the presence of mass asymmetry of reaction while studying the various phenomena of intermediate energy heavy-ion reactions.
Conflict of Interests
The author declares that there is no conflict of interests regarding the publication of this paper.
ScheidW.MüllerH.GreinerW.Nuclear shock waves in heavy-ion collisions197432741745GustafssonH. A.GutbrodH. H.KolbB.Collective flow observed in relativistic nuclear collisions198452159010.1103/PhysRevLett.52.1590OgilvieC. A.CebraD. A.ClaytonJ.Transverse collective motion in intermediate-energy heavy-ion collisions198940, article 259210.1103/PhysRevC.40.2592BlB.KochV.LangA.WeberK.CassingW.MoselU.Origin of transverse momentum in relativistic heavy-ion collisions: microscopic study1991432728AndronicA.ReisdorfW.HerrmannN.Directed flow in Au+Au, Xe+CsI, and Ni+Ni collisions and the nuclear equation of state20036703490710.1103/PhysRevC.67.034907 J. ŁukasikaG. AugerbM.L. Begemann-BlaichDirected and elliptic flow in 197Au + 197Au at intermediate energies20056083-422323010.1016/j.physletb.2004.12.076ZhangY.LiZ.Elliptic flow and system size dependence of transition energies at intermediate energies20067401460210.1103/PhysRevC.74.014602HongB.KimY. J.KangD. H.Proton and deuteron rapidity distributions and nuclear stopping in 96Ru(96Zr) + 96Ru(96Zr) collisions at 400AMeV20026603490110.1103/PhysRevC.66.034901LukasikJ.TrautmannW.Collective flow in heavy ion collisions at intermediate energies2Proceedings of the 23rd International Nuclear Physics Conference (INPC '07)2007513KrofcheckD.BauerW.CrawleyG. M.Disappearance of flow in heavy-ion collisions198963202810.1103/PhysRevLett.63.2028MagestroD. J.BauerW.BjarkiO.Disappearance of transverse flow in Au+Au collisions2000612021602(R)10.1103/PhysRevC.61.021602KumarS.SharmaM. K.PuriR. K.SinghK. P.GovilI. M.Impact parameter dependence of the disappearance of flow and in-medium nucleon-nucleon cross section1998586349434992-s2.0-0032267002SoodA. D.PuriR. K.Study of balance energy in central collisions for heavier nuclei20045943-426026410.1016/j.physletb.2004.05.053SoodA. D.PuriR. K.Mass dependence of disappearance of transverse in-plane flow200469505461210.1103/PhysRevC.69.0546122-s2.0-10144245113SoodA. D.PuriR. K.Influence of momentum-dependent interactions on balance energy and mass dependence200630357157710.1140/epja/i2006-10143-42-s2.0-33845985313SoodA. D.PuriR. K.Systematic study of the energy of vanishing flow: role of equations of state and cross sections200673067602SoodA. D.PuriR. K.Participant-spectator matter at the energy of vanishing flow20097906461810.1103/PhysRevC.79.064618ChughR.PuriR. K.Importance of momentum dependent interactions at the energy of vanishing flow201082014603MagestroD. J.BauerW.WestfallG. D.Isolation of the nuclear compressibility with the balance energy20006244160314160342-s2.0-0034288543GoyalS.PuriR. K.On the sensitivity of the energy of vanishing flow towards mass asymmetry of colliding nuclei2011853116417210.1016/j.nuclphysa.2011.01.022GoyalS.Role of the mass asymmetry of reaction on the geometry of vanishing flow201185615416110.1016/j.nuclphysa.2011.03.004GoyalS.Role of colliding geometry on the balance energy of mass-asymmetric systems20118304760410.1103/PhysRevC.83.047604HartmannO.2003Technische Universität DarmstadtEskolaK. J.Quark and gluon production in high energy nucleus-nucleus collisions19893231375210.1016/0550-3213(89)90586-5SinghJ.2001Chandigarh, IndiaPanjab UniversitySoodA. D.PuriR. K.Nuclear dynamics at the balance energy20047003461110.1103/PhysRevC.70.034611LehmannE.PuriR. K.FaesslerA.BatkoG.HuangS. W.Consequences of a covariant description of heavy-ion reactions at intermediate energies1995512113PuriR. K.HartnackC.AichelinJ.Early fragment formation in heavy-ion collisions1996541R28R3110.1103/PhysRevC.54.R28HartnackC.PuriR. K.AichelinJ.Modelling the many-body dynamics of heavy ion collisions: present status and future perspective1998115116910.1007/s100500050045PuriR. K.AichelinJ.Simulated annealing clusterization algorithm for studying the multifragmentation2000162124526610.1006/jcph.2000.6534KumarS.KumarS.PuriR. K.Medium mass fragment production due to momentum dependent interactions20087806460210.1103/PhysRevC.78.064602VermaniY. K.PuriR. K.Microscopic approach to the spectator matter fragmentation from 400 to 1000AMeV200985562001VermaniY. K.R. K. PuriMass dependence of the onset of multifragmentation in low energy heavy-ion collisions2009361010510310.1088/0954-3899/36/10/105103VermaniY. K.DhawanJ. K.GoyalS.PuriR. K.AichelinJ.Study of fragmentation using clusterization algorithm with realistic binding energies201037101510510.1088/0954-3899/37/1/015105KumarS.KumarS.PuriR. K.Elliptical flow and isospin effects in heavy-ion collisions at intermediate energies20108101461110.1103/PhysRevC.81.014611KumarS.KumarS.PuriR. K.Effect of the symmetry energy on nuclear stopping and its relation to the production of light charged fragments20108101460110.1103/PhysRevC.81.014601AichelinJ.‘Quantum’ molecular dynamics−a dynamical microscopic n-body approach to investigate fragment formation and the nuclear equation of state in heavy ion collisions199120223336010.1016/0370-1573(91)90094-3DhawanJ. K.2007Chandigarh, IndiaPanjab UniversityGautamS.ChughR.SoodA. D.PuriR. K.HartnackC. H.AichelinJ.Isospin effects on the energy of vanishing flow in heavy-ion collisions201037808510210.1088/0954-3899/37/8/085102GautamS.Sood A. D.PuriR. K.AichelinJ.Isospin effects in the disappearance of flow as a function of colliding geometry20118301460310.1103/PhysRevC.83.014603GautamS.SoodA. D.PuriR. K.AichelinJ.Sensitivity of the transverse flow to the symmetry energy20118303460610.1103/PhysRevC.83.034606GautamS.Density and temperature of neutron-rich systems at the energy of vanishing flow in heavy-ion collisions20118306460410.1103/PhysRevC.83.064604GautamS.PuriR. K.Participant-spectator matter and thermalization of neutron-rich systems at the energy of vanishing flow20128506760110.1103/PhysRevC.85.067601KaurS.PuriR. K.Isospin effects on the energy of peak mass production20138701462010.1103/PhysRevC.87.014620KhoaD. T.Ohtsuka N.FaesslerA.MatinM. A.HuangS. W.LehmannE.LoftyY.Microscopic study of thermal properties of the nuclear matter formed in heavy-ion collisions1992542467169810.1016/0375-9474(92)90263-JKhoaD. T.OhtsukaN.MatinM. A.In-medium effects in the description of heavy-ion collisions with realistic NN interactions199254810213010.1016/0375-9474(92)90079-YPuriR. K.OhtsukaN.LehmannE.Temperature-dependent mean field and its effect on heavy-ion reactions19945754733765PuriR. K.LehmannE.FaesslerA.HuangS. W.Study of non-equilibrium effects and thermal properties of heavy-ion collisions using a covariant approach199420111817182810.1088/0954-3899/20/11/0102-s2.0-0040095957