A nuclear density probe: isobaric yield ratio difference

A nuclear density probe: isobaric yield ratio difference 1 2 first-author 3 2 4 2 2 corresp Corresponding author, machunwang@126.com Corresponding author, machunwang@126.com machunwang@126.com Xinxiang Medical University, Xinxiang 453003, China Xinxiang Medical University, Xinxiang 453003, China Institute of Particle and Nuclear Physics, Henan Normal University, Xinxiang 453007, China Institute of Particle and Nuclear Physics, Henan Normal University, Xinxiang 453007, China Henan University of Urban Construction, Pingdingshan 467036, China Henan University of Urban Construction, Pingdingshan 467036, China Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800, China Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800, China We report our recent progress on the nuclear symmetry energy probe, which is called the isobaric yield ratio difference (IBD), and its application in neutron density determination in experiments. The results obtained by the IBD, from which the isobaric yields in the measured 140 A MeV 40, 48 Ca + 9 Be and 58,64 Ni + 9 Be reactions, and the calculated 80 A MeV 38-52 Ca + 12 C reactions by using a modified statistical abrasion-ablation model, show the sensitivity of the IBD to the density differences between reactions. Isobaric yield ratio difference Symmetry energy Neutron density Neutron-rich nucleus I. INTRODUCTION 1 The symmetry energy of asymmetric nuclear matter, especially neutron-rich matter, has been an important topic both in theory and experiment with the development of radioactive nuclear beam (RNB) facilities. The new or updated RNB facilities can provide nuclear beams with larger neutron-richness, which are supposed to promote the study of neutron-rich nuclear matter and extend the chart of nuclides to the more neutron-rich side. In heavy-ion collisions above the intermediate energy, nuclear matter from sub-saturation to supra-saturation densities can be produced. The symmetry energy of nuclear matter depends on density and temperature. In the past decade, though many probes have been proposed, there are still difficulties in measuring the symmetry energy of nuclear matter at different nuclear densities. The isobaric yield ratio (IYR) method was first proposed by Huang et al . to extract the symmetry-energy coefficient for fragments with small isospins, which also have some temperature since they are measured in heavy-ion collisions [ 1 1 , 2 2 ]. In the IYR, some terms deciding the fragment yield cancel out, which makes it useful to study the symmetry energy term of the fragment with a low temperature. The IYR method is extended to study the symmetry energy of fragment with larger neutron richness [ 3 3 - 6 6 3-6 ]. The results obtained by the IYR method are compared to a similar method, i.e., isoscaling, which extracts the symmetry energy of the colliding source, showing a large difference between them [ 7 7 - 10 10 7-10 ]. We studied the results of the IYR and isoscaling, and found that though the IYR and isoscaling methods both aim at the extraction of nuclear symmetry energy, the obtained results are for different nuclear matter. A new probe was constructed via IYR between two reactions, which is called the isobaric yield ratio difference (IBD), and the results of the IBD and isoscaling were found to be similar because they are both the probes for the chemical potential of neutrons and protons in the reactions [ 11 11 , 12 12 ]. Moreover, the IBD probe is also sensitive to the density differences between the reactions [ 13 13 , 14 14 ]. In this article, we will summarize the recent progress on the IBD method. II. THEORETICAL DESCRIPTIONS 1 In thermodynamics theory, an equilibrium state is assumed for the reaction, which introduces the concept of temperature. The temperature affects the yield of fragment. With high temperature, light particles like neutrons and protons can escape from the binding of hot fragment. The formation time of the final fragments, which correspond to the measured ones in experiments, is called the "chemical freeze-out" (when the fragment ceases to emit light particles, like neutrons, protons or α, and the neutron and proton numbers do not change anymore). In general, the yield of fragment in heavy-ion collisions above the Fermi energy obeys an exponential law. For example, in the grand-canonical ensembles theory within the grand canonical limit, the yield of a fragment with mass numbers A and neutron-excess I ( I = N – Z is the difference between the numbers of neutrons and protons) is written as [ 15 15 , 16 16 ], <math display="block" id="M1"><mrow><mi>σ</mi><mo stretchy="false">(</mo><mi>I</mi><mn>,</mn><mi>A</mi><mo stretchy="false">)</mo><mo>=</mo><mi>c</mi><msup><mi>A</mi><mi>τ</mi></msup><mi>exp</mi><mo>{</mo><mo stretchy="false">[</mo><mi>F</mi><mo stretchy="false">(</mo><mi>I</mi><mn>,</mn><mi>A</mi><mo stretchy="false">)</mo><mo>+</mo><mi>N</mi><msub><mi>μ</mi><mi>n</mi></msub><mo>+</mo><mi>Z</mi><msub><mi>μ</mi><mi>p</mi></msub><mo stretchy="false">]</mo><mo>/</mo><mi>T</mi><mo>}</mo><mn/></mrow></math> where T is the temperature of the system. F ( I , A ) is the free energy of the fragment which depends on T . μ n ( μ p ) is the chemical potential of the neutrons (protons) in the equilibrium system, which depends on the density and T of the systems. In the isoscaling method [ 17 17 , 18 18 ], the difference between the μ n ( μ p ) of two reactions can be extracted from the isotopic (isotonic) ratio, <math display="block" id="M2"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mi>R</mi><msub><mrow><mo stretchy="false">(</mo><mi>I</mi><mn>,</mn><mi>A</mi><mo stretchy="false">)</mo></mrow><mrow><mn>21</mn></mrow></msub><mo>=</mo><mfrac><mrow><mi>σ</mi><msub><mrow><mo stretchy="false">(</mo><mi>I</mi><mn>,</mn><mi>A</mi><mo stretchy="false">)</mo></mrow><mn>2</mn></msub></mrow><mrow><mi>σ</mi><msub><mrow><mo stretchy="false">(</mo><mi>I</mi><mn>,</mn><mi>A</mi><mo stretchy="false">)</mo></mrow><mn>1</mn></msub></mrow></mfrac></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mo>∝</mo><mi>exp</mi><mo>{</mo><mo stretchy="false">[</mo><mi>N</mi><mo stretchy="false">(</mo><msub><mi>μ</mi><mrow><mi>n</mi><mn>2</mn></mrow></msub><mo>−</mo><msub><mi>μ</mi><mrow><mi>n</mi><mn>1</mn></mrow></msub><mo stretchy="false">)</mo><mo>+</mo><mi>Z</mi><mo stretchy="false">(</mo><msub><mi>μ</mi><mrow><mi>p</mi><mn>2</mn></mrow></msub><mo>−</mo><msub><mi>μ</mi><mrow><mi>p</mi><mn>1</mn></mrow></msub><mo stretchy="false">)</mo><mo stretchy="false">]</mo><mo>/</mo><mi>T</mi><mo>}</mo><mn>,</mn></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow/></mtd></mtr></mtable></mrow></math> with indices 1 and 2 denoting the reactions. In Eq. (2), the T of the two reactions is assumed to be the same. For the chemical freeze-out fragments this can be true since they both have a very long evolution time. For the isotopic ratio, <math display="block" id="M3"><mrow><mi>ln</mi><mi>R</mi><msub><mrow><mo stretchy="false">(</mo><mi>I</mi><mn>,</mn><mi>A</mi><mo stretchy="false">)</mo></mrow><mrow><mn>21</mn></mrow></msub><mo>=</mo><mo stretchy="false">[</mo><mi>N</mi><mo stretchy="false">(</mo><msub><mi>μ</mi><mrow><mi>n</mi><mn>2</mn></mrow></msub><mo>−</mo><msub><mi>μ</mi><mrow><mi>n</mi><mn>1</mn></mrow></msub><mo stretchy="false">)</mo><mo>+</mo><mi>Z</mi><mo stretchy="false">(</mo><msub><mi>μ</mi><mrow><mi>p</mi><mn>2</mn></mrow></msub><mo>−</mo><msub><mi>μ</mi><mrow><mi>p</mi><mn>1</mn></mrow></msub><mo stretchy="false">)</mo><mo stretchy="false">]</mo><mo>/</mo><mi>T</mi><mn>,</mn></mrow></math> Z does not change. The isotopic ratio shows scaling phenomena obeying the linear function. The fitted slope is defined as α ( μ n 2 - μ n 1 )/ T . Similarly, for the isotonic ratio, N is a constant. The isotonic ratio also shows scaling phenomena obeying the linear function. The fitted slope is defined as β ( μ p 2 - μ p 1 )/ T . α and β are used to extract the symmetry energy of the system. But the symmetry energy in the isoscaling method is an indirect result obtained from the yield of fragments. Starting from Eq. (1), assuming the isobars differing 2 units in I have the same temperature, the isobaric yield ratio can be defined as, <math display="block" id="M4"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mi>R</mi><msub><mrow><mo stretchy="false">(</mo><mi>I</mi><mo>+</mo><mn>2,</mn><mi>I</mi><mn>,</mn><mi>A</mi><mo stretchy="false">)</mo></mrow><mrow><mn>21</mn></mrow></msub><mo>=</mo><mfrac><mrow><mi>σ</mi><mo stretchy="false">(</mo><mi>I</mi><mo>+</mo><mn>2,</mn><mi>A</mi><mo stretchy="false">)</mo></mrow><mrow><mi>σ</mi><mo stretchy="false">(</mo><mi>I</mi><mn>,</mn><mi>A</mi><mo stretchy="false">)</mo></mrow></mfrac></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mo>=</mo><mi>exp</mi><mo>{</mo><mo stretchy="false">[</mo><mi>F</mi><mo stretchy="false">(</mo><mi>I</mi><mo>+</mo><mn>2,</mn><mi>A</mi><mo stretchy="false">)</mo><mo>−</mo><mi>F</mi><mo stretchy="false">(</mo><mi>I</mi><mn>,</mn><mi>A</mi><mo stretchy="false">)</mo><mo>+</mo><msub><mi>μ</mi><mi>n</mi></msub><mo>−</mo><msub><mi>μ</mi><mi>p</mi></msub><mo stretchy="false">]</mo><mo>/</mo><mi>T</mi><mo>}</mo><mn>,</mn></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow/></mtd></mtr></mtable></mrow></math> To make cancellation of F , we adopt a reference reaction system with similar measurements as the isoscaling method does, and assume that the free energies of one fragment in the two reactions are the same. Defined as the difference between the IYRs of two reactions, the IBD is [ 11 11 , 12 12 ], <math display="block" id="M5"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mi>Δ</mi><mi>ln</mi><mi>R</mi><msub><mrow><mo stretchy="false">(</mo><mi>I</mi><mo>+</mo><mn>2,</mn><mi>I</mi><mn>,</mn><mi>A</mi><mo stretchy="false">)</mo></mrow><mrow><mn>21</mn></mrow></msub></mrow></mtd><mtd columnalign="left"><mrow><mo>=</mo><mi>ln</mi><msub><mi>R</mi><mn>2</mn></msub><mo stretchy="false">(</mo><mi>I</mi><mo>+</mo><mn>2,</mn><mi>I</mi><mn>,</mn><mi>A</mi><mo stretchy="false">)</mo><mo>−</mo><mi>ln</mi><msub><mi>R</mi><mn>1</mn></msub><mo stretchy="false">(</mo><mi>I</mi><mo>+</mo><mn>2,</mn><mi>I</mi><mn>,</mn><mi>A</mi><mo stretchy="false">)</mo></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mo>=</mo><mo stretchy="false">[</mo><mo stretchy="false">(</mo><msub><mi>μ</mi><mrow><mi>n</mi><mn>2</mn></mrow></msub><mo>−</mo><msub><mi>μ</mi><mrow><mi>n</mi><mn>1</mn></mrow></msub><mo stretchy="false">)</mo><mo>−</mo><mo stretchy="false">(</mo><msub><mi>μ</mi><mrow><mi>p</mi><mn>2</mn></mrow></msub><mo>−</mo><msub><mi>μ</mi><mrow><mi>p</mi><mn>1</mn></mrow></msub><mo stretchy="false">)</mo><mo stretchy="false">]</mo><mo>/</mo><mi>T</mi></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mo>=</mo><mo stretchy="false">(</mo><mi>Δ</mi><msub><mi>μ</mi><mrow><mi>n</mi><mn>21</mn></mrow></msub><mo>−</mo><mi>Δ</mi><msub><mi>μ</mi><mrow><mi>p</mi><mn>21</mn></mrow></msub><mo stretchy="false">)</mo><mo>/</mo><mi>T</mi><mo>=</mo><mi>Δ</mi><msub><mi>μ</mi><mrow><mn>21</mn></mrow></msub><mo>/</mo><mi>T</mi><mn>.</mn></mrow></mtd></mtr></mtable></mrow></math> The indices also denote the reaction systems. From the definition of the isoscaling parameters, the following correlation can be obtained, <math display="block" id="M6"><mrow><mi>Δ</mi><mi>ln</mi><mi>R</mi><msub><mrow><mo stretchy="false">(</mo><mi>I</mi><mo>+</mo><mn>2,</mn><mi>I</mi><mn>,</mn><mi>A</mi><mo stretchy="false">)</mo></mrow><mrow><mn>21</mn></mrow></msub><mo>=</mo><mi>α</mi><mo>−</mo><mi>β</mi><mn>.</mn></mrow></math> Eq. (5) can also be written as, <math display="block" id="M7"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mi>Δ</mi><mi>ln</mi><mi>R</mi><msub><mrow><mo stretchy="false">(</mo><mi>I</mi><mo>+</mo><mn>2,</mn><mi>I</mi><mn>,</mn><mi>A</mi><mo stretchy="false">)</mo></mrow><mrow><mn>21</mn></mrow></msub></mrow></mtd><mtd columnalign="left"><mrow><mo>=</mo><mo stretchy="false">[</mo><mo stretchy="false">(</mo><msub><mi>μ</mi><mrow><mi>n</mi><mn>2</mn></mrow></msub><mo>−</mo><msub><mi>μ</mi><mrow><mi>p</mi><mn>1</mn></mrow></msub><mo stretchy="false">)</mo><mo>−</mo><mo stretchy="false">(</mo><msub><mi>μ</mi><mrow><mi>n</mi><mn>1</mn></mrow></msub><mo>−</mo><msub><mi>μ</mi><mrow><mi>p</mi><mn>1</mn></mrow></msub><mo stretchy="false">)</mo><mo stretchy="false">]</mo><mo>/</mo><mi>T</mi></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mo>=</mo><mo stretchy="false">(</mo><mi>Δ</mi><msub><mi>μ</mi><mrow><mi>n</mi><mi>p</mi><mn>2</mn></mrow></msub><mo>−</mo><mi>Δ</mi><msub><mi>μ</mi><mrow><mi>n</mi><mi>p</mi><mn>1</mn></mrow></msub><mo stretchy="false">)</mo><mo>/</mo><mi>T</mi><mn>.</mn></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow/></mtd></mtr></mtable></mrow></math> The chemical potentials of neutrons and protons are related to the densities of neutrons ( ρ n ) and protons ( ρ p ) respectively, [ 14 14 , 19 19 ], <math display="block" id="M8"><mrow><msub><mi>μ</mi><mi>n</mi></msub><mo>/</mo><mi>T</mi><mo>=</mo><mi>ln</mi><msub><mi>ρ</mi><mi>n</mi></msub><mn>,</mn><mtext> </mtext><mtext>and</mtext><mtext> </mtext><msub><mi>μ</mi><mi>p</mi></msub><mo>/</mo><mi>T</mi><mo>=</mo><mi>ln</mi><msub><mi>ρ</mi><mi>p</mi></msub><mn>,</mn></mrow></math> which result in α = ln ρ n 2 - ln ρ n 1 = ln( ρ n 2 / ρ n 1 ), and β = ln ρ p 2 - ln ρ p 1 = ln( ρ p 2 / ρ p 1 ). We can also write the relationship between the IBD results and the densities of the reaction systems as [ 14 14 ], <math display="block" id="M9"><mrow><mi>Δ</mi><mi>ln</mi><mi>R</mi><msub><mrow><mo stretchy="false">(</mo><mi>I</mi><mo>+</mo><mn>2,</mn><mi>I</mi><mn>,</mn><mi>A</mi><mo stretchy="false">)</mo></mrow><mrow><mn>21</mn></mrow></msub><mo>=</mo><mi>ln</mi><mo stretchy="false">(</mo><msub><mi>ρ</mi><mrow><mi>n</mi><mn>2</mn></mrow></msub><mo>/</mo><msub><mi>ρ</mi><mrow><mi>n</mi><mn>1</mn></mrow></msub><mo stretchy="false">)</mo><mo>−</mo><mi>ln</mi><mo stretchy="false">(</mo><msub><mi>ρ</mi><mrow><mi>p</mi><mn>2</mn></mrow></msub><mo>/</mo><msub><mi>ρ</mi><mrow><mi>p</mi><mn>1</mn></mrow></msub><mo stretchy="false">)</mo><mn>,</mn></mrow></math> <math display="block" id="M10"><mrow><mi>Δ</mi><mi>ln</mi><mi>R</mi><msub><mrow><mo stretchy="false">(</mo><mi>I</mi><mo>+</mo><mn>2,</mn><mi>I</mi><mn>,</mn><mi>A</mi><mo stretchy="false">)</mo></mrow><mrow><mn>21</mn></mrow></msub><mo>=</mo><mi>ln</mi><mo stretchy="false">(</mo><msub><mi>ρ</mi><mrow><mi>n</mi><mn>2</mn></mrow></msub><mo>/</mo><msub><mi>ρ</mi><mrow><mi>p</mi><mn>2</mn></mrow></msub><mo stretchy="false">)</mo><mo>−</mo><mi>ln</mi><mo stretchy="false">(</mo><msub><mi>ρ</mi><mrow><mi>n</mi><mn>1</mn></mrow></msub><mo>/</mo><msub><mi>ρ</mi><mrow><mi>p</mi><mn>1</mn></mrow></msub><mo stretchy="false">)</mo><mn>.</mn></mrow></math> Eqs. (9) and (10) are written according to Eqs. (5) and (7), respectively. The densities of protons can be determined by experimental measurement. Thus the relationships in Eqs. (9) and (10) can be used to determine the neutron density of the neutron-rich system. For example, if the reference reaction system is a neutron-proton symmetric one, in which ρ n 1 and ρ p 1 can be assumed to be the same and ρ p 2 is measurable, ρ n 2 can be determined from the IBD results. In Eq. (9), if ρ p 2 and ρ p 1 can be assumed to be the same (for isotopic projectiles) and ρ p 1 ρ n 1 , Δ R ( I + 2, I , A ) 21 = ln ( ρ n 2 / ρ p 1 ). Similarly, in Eq. (10), the reference reaction is neutron-proton symmetric and ρ p 1 ρ n 1 , by knowing ρ p 2 , Δ R ( I + 2, I , A ) 21 = ln ( ρ n 2 / ρ p 2 ) and ρ n 2 can be determined from the IBD. III. RESULTS AND DISCUSSION 1 Now we will show the IBD results obtained from the measured 140 A MeV 40, 48 Ca + 9 Be and 58,64 Ni + 9 Be projectile fragmentation reactions [ 20 20 ]. These reactions were measured by M. Mocko et al . at the National Superconducting Cyclotron Laboratory (NSCL) at Michigan State University. The measured fragments range from Z = 5 to Z proj , which provides high quality data for the test of the IBD methods. The IBD results has been previously In Fig. 1 Fig. 1 , the results of the IBD and isoscaling methods are re-plotted for comparison [ 11 11 ]. Four groups of reactions are analyzed, the 48 Ca/ 40 Ca + 9 Be and 64 Ni/ 58 Ni + 9 Be reactions of the isotopic projectiles, the 58 Ni/ 40 Ca + 9 Be reactions with the projectiles both being neutron-proton symmetric, and the 48 Ca/ 64 Ni + 9 Be reactions with the projectiles both being neutron-rich. The reactions using the relative neutron-proton symmetric projectile was set as the reference reaction denoted by index 1. Theoretically, the relationship between the IBD and isoscaling is described in Eq. (6). The fragment belongs to one isotopic chain and one isotonic chain, thus it can be related to α and β according to the respective isotopic and isotonic ratio fitting. For the isobaric chains with I from -1 to 2, the IBD and isoscaling results were found to be similar, except for some fragments with a large A . The IBD results were the same in the I = 0 and 1 isobaric chains for the 48 Ca/ 40 Ca reactions. The same phenomena can be found in the I = 1 and 2 isobaric chains for the 64 Ni/ 58 Ni reactions, and the I = -1 isobaric chains for the 58 Ni/ 40 Ca reactions. It was also found that the height of the IBD-plateaus in the four groups were different, about 2, 1.4, 0.5 and 0.4, respectively [ 11 11 ]. And the lengths of the plateaus are also different. In the IBD of neutron-rich isobars, the plateau can disappear [ 11 11 - 14 14 11-14 ]. The IBD results denote the difference between the neutrons and protons for the two reactions according to Eqs. (9) and (10). If the IBD results change very small, it reflects that the difference between the densities of neutrons and protons for the reactions are small. The increasing IBD with A of the fragment indicates an enlarged difference between the densities of the neutrons and protons [ 21 21 , 22 22 ]. Taking the IBD results for the 48 Ca/ 40 Ca reactions as an example, ρ p 1 and ρ n 1 for 40 Ca can be assumed to be the same, and assuming that ρ p 1 ρ p 2 , Δ R ( I + 2, I , A ) 21 = ln( ρ n 2 / ρ p 1 ). ρ n 2 of the 48 Ca reaction can be determined from the IBD results. For the 48 Ca/ 40 Ca reaction, <math id="M11"><mrow><mi>Δ</mi><mi>R</mi><msub><mrow><mo stretchy="false">(</mo><mi>I</mi><mo>+</mo><mn>2,</mn><mi>I</mi><mn>,</mn><mi>A</mi><mo stretchy="false">)</mo></mrow><mrow><mn>21</mn></mrow></msub><mo>=</mo><mi>ln</mi><mo stretchy="false">(</mo><mfrac><mrow><msub><mi>ρ</mi><mrow><mi>n</mi><mn>2</mn></mrow></msub></mrow><mrow><msub><mi>ρ</mi><mrow><mi>p</mi><mn>1</mn></mrow></msub></mrow></mfrac><mo stretchy="false">)</mo><mo>=</mo><mn>2</mn></mrow></math> which results in <math id="M12"><mrow><mfrac><mrow><msub><mi>ρ</mi><mrow><mi>n</mi><mn>2</mn></mrow></msub></mrow><mrow><msub><mi>ρ</mi><mrow><mi>p</mi><mn>1</mn></mrow></msub></mrow></mfrac><mo>=</mo><mn>7.4</mn></mrow></math> . It is very large compared to the isoscaling method. It should be noted that 48 Ca has a much larger neutron excess than 40 Ca, and α∼β, of which β cannot be omitted in Eq. (6). This makes Δln R ( I +2, I , A ) 21 = 2 α , which results in the overestimation of <math id="M13"><mrow><mfrac><mrow><msub><mi>ρ</mi><mrow><mi>n</mi><mn>2</mn></mrow></msub></mrow><mrow><msub><mi>ρ</mi><mrow><mi>p</mi><mn>1</mn></mrow></msub></mrow></mfrac></mrow></math> . Similar explanations can be found for the other three reaction groups. It should be noted that the densities determined by the IBD results are for the time of freeze-out, assuming thermal equilibrium. The relationship between the density of the projectile and the system at the freeze-out should be further studied in transport models. The IBD-plateau was found to be sensitive to differences between the densities of the reaction systems. The asymmetry of these systems (N/Z) 21 are calculated and the results are plotted in Fig. 2 Fig. 2 . The asymmetry of the reaction systems (N/Z) 21 is defined as [(N p +N t )/(Z p +Z t )] 2 /[(N p +N t )/(Z p +Z t )] 1 . From high to low, (N/Z) 21 is in the order of 48 Ca/ 40 Ca, 48 Ca/ 58 Ni, 64 Ni/ 40 Ca, 64 Ni/ 58 Ni, 48 Ca/ 64 Ni, and 58 Ni/ 40 Ca. The (N/Z) 21 of the 48 Ca/ 64 Ni and 58 Ni/ 40 Ca reactions are very similar. In Fig. 3(a) Fig. 3(a) , the IBD-results for the I = 0 isobaric chain are re-plotted. From high to low, the IBD-plateau is in the order of 48 Ca/ 40 Ca, 48 Ca/ 58 Ni, 64 Ni/ 40 Ca, 64 Ni/ 58 Ni, 58 Ni/ 40 Ca, and 48 Ca/ 64 Ni. The IBD-plateaus for the 58 Ni/ 40 Ca and 48 Ca/ 64 Ni reactions are very similar. The order of the (N/Z) 21 and the IBD-plateau are consistent, except for the 58 Ni/ 40 Ca and 48 Ca/ 64 Ni reactions. The same phenomena can also found in the IBD results for the I = 1 isobaric chain in the reactions as plotted in Fig. 3(b) Fig. 3(b) . To test the sensitivity of the IBD-plateau to the density difference between the reactions, we calculate the fragment yield in the 80 A MeV 38-52 Ca + 12 C reactions (only the even-A calcium isotopes are used) by using a modified statistical abrasion-ablation (SAA) model [ 6 6 , 13 13 , 14 14 , 21 21 - 27 27 21-27 ]. The statistical abrasion-ablation model will not be described in this work. Readers can refer to Refs. [ 21 21 - 27 27 21-27 ]. The fragments produced in the 40 Ca + 12 C reaction is set as reaction 1, which is the reference reaction. The IBD results for the isobaric chains from I =-1 to 2 are plotted in Fig. 3 Fig. 3 . Since the analyzed fragments are the prefragments which are determined by the neutron and proton densities, as well as the nucleus-nucleus reaction cross sections, the yields of prefragments are sensitive the density changes of the projectile [ 13 13 ]. It is found that, as the projectile ( X Ca, with X denotes the mass number of projectile) becomes more neutron-rich, the IBD results increases regularly. The IBD obtained from the prefragments are sensitive to the difference between the neutron density of the X Ca and 40 Ca reactions if we assume that the proton densities in X Ca are the same or vary very little [ 13 13 , 14 14 ]. We calculated the average values ( < Δ μ / T > ) of the IBD-plateau for the isobars plotted in Fig. 4 Fig. 4 . The < Δ μ / T > are plotted as a function of I in Fig. 5 Fig. 5 . It is shown that for the relative symmetry fragments, < Δ μ / T > changes very small with I. In the inserted figure of Fig. 4 Fig. 4 , the difference between the neutron densities (Δ ln ρ n ) of X Ca and 40 Ca is plotted. When r < 2 fm, Δ ln ρ n are flat and change very small. It has been illustrated that the yield of fragments with small mass numbers are mostly produced in the central collisions [ 21 21 , 22 22 , 26 26 , 27 27 ]. The flat distribution of Δ ln ρ n between the core part of the calcium isotopes can account for the phenomena of the IBD-plateau. The enlarged Δ ln ρ n with r ( r > 2 fm for different calcium isotopes) can also account for the increasing of Δ μ /T with A , as well as the disappearance of the IBD-plateau in the neutron-rich isobars [ 11 11 - 14 14 11-14 ]. Both from the experiment results and the calculated results by the statistical abrasion-ablation model [ 11 11 - 14 14 11-14 ], it is shown that the IBD probe is sensitive to the density difference between reactions. The IBD probe provides a new method to determine the neutron density in heavy-ion collisions. But as we have noted, the density is for the time of the chemical freeze-out, assuming there is a thermal equilibrium in the reaction. Actually, the process of the reaction is dynamic and changes with reaction time. The IBD results should be verified by a dynamic evolution process. Previously we have reported the IBD result for the 58 Ni + 12 C reaction calculated by the anti-symmetric molecular dynamical (AMD) model to verify the sensitivity of the IBD to density [ 11 11 ]. In dynamical model calculations, the evolution of the IBD results with time can be studied. We have finished the calculations and the IBD analysis of the 58, 64 Ni + 12 C reactions and hope to report the AMD results in near future to see the IBD for primary fragments and fragments after the decay process. Moreover, the Shannon information entropy uncertainty is introduced to explain the physical meaning of the IBD probe [ 28 28 , 29 29 ]. In the dynamical model, the construct of fragments are important [ 30 30 , 31 31 ], and at the same time the temperature for heavy fragments should be carefully treated [ 32 32 - 39 39 32-39 ]. It is also suggested that the IBD can be used to study the evolution of chemical potential and the density of the reactions. We will also study the evolution of the IBD results by using the AMD model. IV. SUMMARY 1 In summary, the method of isobaric yield ratio difference is explained in this work. The IBD method is a more direct method than the isoscaling method to extract the symmetry energy or density of nuclear matter in heavy-ion collisions. It is proven that the IBD probe is sensitive to the difference between the densities of the reactions, and can be used to study the neutron density of the reactions. References M Huang , Z Chen , S Kowalski , et al . Isobaric yield ratios and the symmetry energy in heavy-ion reactions near the Fermi energy . Phys Rev C , 2010 , 81 : 044620 . DOI: 10.1103/PhysRevC.81.044620 http://doi.org/10.1103/PhysRevC.81.044620 Isobaric yield ratios and the symmetry energy in heavy-ion reactions near the Fermi energy M Huang A Bonasera , Z Chen et al . Isospin dependence of the nuclear equation of state near the critical point . Phys Rev C , 2010 , 81 : 044618 . DOI: 10.1103/PhysRevC.81.044618 http://doi.org/10.1103/PhysRevC.81.044618 Isospin dependence of the nuclear equation of state near the critical point C W Ma , F Wang , Y G Ma , et al . Isobaric yield ratios in heavy-ion reactions, and symmetry energy of neutron-rich nuclei at intermediate energies . Phys Rev C , 2011 , 83 : 064620 . DOI: 10.1103/PhysRevC.83.064620 http://doi.org/10.1103/PhysRevC.83.064620 Isobaric yield ratios in heavy-ion reactions, and symmetry energy of neutron-rich nuclei at intermediate energies C W Ma , J Pu , H L Wei , et al . Symmetry energy extracted from fragments in relativistic energy heavy-ion collisions induced by ^124,136Xe . Eur Phys J A , 2012 , 48 : 78 . DOI: 10.1140/epja/i2012-12078-5 http://doi.org/10.1140/epja/i2012-12078-5 Symmetry energy extracted from fragments in relativistic energy heavy-ion collisions induced by ^124,136Xe C W Ma , J Pu , S S Wang , et al . The Symmetry Energy from the Neutron-Rich Nucleus Produced in the Intermediate-Energy ^40,48Ca and ^58,64Ni Projectile Fragmentation . Chin Phys Lett , 2012 , 29 : 062101 . DOI: 10.1088/0256-307X/29/6/062101 http://doi.org/10.1088/0256-307X/29/6/062101 The Symmetry Energy from the Neutron-Rich Nucleus Produced in the Intermediate-Energy ^40,48Ca and ^58,64Ni Projectile Fragmentation C W Ma , H L Song , J Pu , et al . Symmetry energy from neutron-rich fragments in heavy-ion collisions, and its dependence on incident energy, and impact parameters . Chin Phys C , 2013 , 37 : 024102 . DOI: 10.1088/1674-1137/37/2/024102 http://doi.org/10.1088/1674-1137/37/2/024102 Symmetry energy from neutron-rich fragments in heavy-ion collisions, and its dependence on incident energy, and impact parameters P Marini , A Bonasera , A Mclntosh , et al . Constraining the symmetry term in the nuclear equation of state at subsaturation densities and finite temperatures . Phys Rev C , 2012 , 85 : 034617 . DOI: 10.1103/PhysRevC.85.034617 http://doi.org/10.1103/PhysRevC.85.034617 Constraining the symmetry term in the nuclear equation of state at subsaturation densities and finite temperatures P Marini , A Bonasera , G A Souliotis , et al . Systematic study of the symmetry energy within the approach of the statistical multifragmentation model . Phys Rev C , 2013 , 87 : 024603 . DOI: 10.1103/PhysRevC.87.024603 http://doi.org/10.1103/PhysRevC.87.024603 Systematic study of the symmetry energy within the approach of the statistical multifragmentation model S Mallik and G Chaudhuri . Symmetry energy from nuclear multifragmentation . Phys Rev C , 2013 , 87 : 011602 . DOI: 10.1103/PhysRevC.87.011602 http://doi.org/10.1103/PhysRevC.87.011602 Symmetry energy from nuclear multifragmentation S Mallik and G Chaudhuri . Effect of particle fluctuation on isoscaling and isobaric yield ratio of nuclear multifragmentation . Phys Lett B , 2013 , 727 : 282 - 286 . DOI: 10.1016/j.physletb.2013.10.011 http://doi.org/10.1016/j.physletb.2013.10.011 Effect of particle fluctuation on isoscaling and isobaric yield ratio of nuclear multifragmentation C W Ma , S S Wang , Y L Zhang , et al . Isobaric yield ratio difference in heavy-ion collisions, and comparison to isoscaling . Phys Rev C , 2013 , 87 : 034618 . DOI: 10.1103/PhysRevC.87.034618 http://doi.org/10.1103/PhysRevC.87.034618 Isobaric yield ratio difference in heavy-ion collisions, and comparison to isoscaling C W Ma , S S Wang , Y L Zhang , et al . Chemical properties of colliding sources in ^124,136Xe and 112,124Sn induced collisions in isobaric yield ratio difference and isoscaling methods . J Phys G Nucl Partic , 2013 , 40 : 125106 . DOI: 10.1088/0954-3899/40/12/125106 http://doi.org/10.1088/0954-3899/40/12/125106 Chemical properties of colliding sources in ^124,136Xe and 112,124Sn induced collisions in isobaric yield ratio difference and isoscaling methods C W Ma , J Yu , X M Bai , et al . Isobaric yield ratio difference and neutron density difference in calcium isotopes . Phys Rev C , 2014 , 89 : 057602 . DOI: 10.1103/PhysRevC.89.057602 http://doi.org/10.1103/PhysRevC.89.057602 Isobaric yield ratio difference and neutron density difference in calcium isotopes C W Ma , X M Bao , J Yu , et al . Neutron density distributions of neutron-rich nuclei studied with the isobaric yield ratio difference . Eur Phys J A , 2014 , 50 : 139 . DOI: 10.1140/epja/i2014-14139-1 http://doi.org/10.1140/epja/i2014-14139-1 Neutron density distributions of neutron-rich nuclei studied with the isobaric yield ratio difference C B Das , S D Gupta , X D Liu , et al . Comparison of canonical and grand canonical models for selected multifragmentation data . Phys Rev C , 2001 , 64 : 044608 . DOI: 10.1103/PhysRevC.64.044608 http://doi.org/10.1103/PhysRevC.64.044608 Comparison of canonical and grand canonical models for selected multifragmentation data M B Tsang , W G Lynch , W A Friedman , et al . Fragmentation cross sections and binding energies of neutron-rich nuclei . Phys Rev C , 2007 , 76 : 041302(R) . DOI: 10.1103/PhysRevC.76.041302 http://doi.org/10.1103/PhysRevC.76.041302 Fragmentation cross sections and binding energies of neutron-rich nuclei H S Xu , M B Tsang , T X Liu , et al . Isospin fractionation in nuclear multifragmentation . Phys Rev Lett , 2000 , 85 : 716 . DOI: 10.1103/PhysRevLett.85.716 http://doi.org/10.1103/PhysRevLett.85.716 Isospin fractionation in nuclear multifragmentation M B Tsang , Friedman , C K Gelbke , et al . Isotopic Scaling in Nuclear Reactions . Phys Rev Lett , 2001 , 86 : 5023 . DOI: 10.1103/PhysRevLett.86.5023 http://doi.org/10.1103/PhysRevLett.86.5023 Isotopic Scaling in Nuclear Reactions E Geraci , M Bruno , M D’Agostino , et al . Isoscaling in central 124Sn+64Ni, 112Sn+58Ni collisions at 35 A MeV . Nucl Phys A , 2004 , 732 : 173 - 201 . DOI: 10.1016/j.nuclphysa.2003.11.055 http://doi.org/10.1016/j.nuclphysa.2003.11.055 Isoscaling in central 124Sn+64Ni, 112Sn+58Ni collisions at 35 A MeV M Mocko , M B Tsang , L Andronenko , et al . Projectile fragmentation of 40Ca, 48Ca, 58Ni, and 64Ni at 140 MeV/nucleon . Phys Rev C , 2006 , 74 : 054612 . DOI: 10.1103/PhysRevC.74.054612 http://doi.org/10.1103/PhysRevC.74.054612 Projectile fragmentation of 40Ca, 48Ca, 58Ni, and 64Ni at 140 MeV/nucleon C W Ma , S S Wang , H L Wei , et al . Re-examination of finite-size effects in isobaric yield ratios using a statistical Abrasion-Ablation model . Chin Phys Lett , 2013 , 30 : 052101 . DOI: 10.1088/0256-307X/30/5/052501 http://doi.org/10.1088/0256-307X/30/5/052501 Re-examination of finite-size effects in isobaric yield ratios using a statistical Abrasion-Ablation model C W Ma , H L Wei and Y G Ma . Neutron-skin effects in isobaric yield ratios for mirror nuclei in a statistical abrasion-ablation model . Phys Rev C , 2013 , 88 , 044612 . DOI: 10.1103/PhysRevC.88.044612 http://doi.org/10.1103/PhysRevC.88.044612 Neutron-skin effects in isobaric yield ratios for mirror nuclei in a statistical abrasion-ablation model J J Gaimard and S H Schmidt . A reexamination of the abrasion-ablation model for the description of the nuclear fragmentation reaction . Nucl Phys A , 1991 , 531 : 709 - 745 . DOI: 10.1016/0375-9474(91)90748-U http://doi.org/10.1016/0375-9474(91)90748-U A reexamination of the abrasion-ablation model for the description of the nuclear fragmentation reaction T Brohm and K H Schmidt . Statistical abrasion of nucleons from realistic nuclear-matter distributions . Nucl Phys A , 1994 , 569 : 821 - 832 . DOI: 10.1016/0375-9474(94)90386-7 http://doi.org/10.1016/0375-9474(94)90386-7 Statistical abrasion of nucleons from realistic nuclear-matter distributions D Q Fang , W Q Shen , J Feng , et al . Isospin effect of fragmentation reactions induced by intermediate energy heavy ions and its disappearance . Phys Rev C , 2000 , 61 : 044610 . DOI: 10.1103/PhysRevC.61.044610 http://doi.org/10.1103/PhysRevC.61.044610 Isospin effect of fragmentation reactions induced by intermediate energy heavy ions and its disappearance C W Ma , H L Wei , J Y Wang , et al . Isospin dependence of projectile-like fragment production at intermediate energies . Phys Rev C , 2009 , 79 : 034606 . DOI: 10.1103/PhysRevC.79.034606 http://doi.org/10.1103/PhysRevC.79.034606 Isospin dependence of projectile-like fragment production at intermediate energies C W Ma , H L Wei , G J Liu , et al . Systematic behavior in the isospin dependence of projectile fragmentation of mirror nuclei with A = 20-60 . J Phys G Nucl Partic , 2010 , 37 : 015104 . DOI: 10.1088/0954-3899/37/1/015104 http://doi.org/10.1088/0954-3899/37/1/015104 Systematic behavior in the isospin dependence of projectile fragmentation of mirror nuclei with A = 20-60 C W Ma and H L Wei . Isotopic ratio, isotonic ratio, isobaric ratio and Shannon information uncertainty . Commun Theor Phys , 2014 , 62 : 717 - 720 . Isotopic ratio, isotonic ratio, isobaric ratio and Shannon information uncertainty C W Ma , H L Wei , S S Wang , et al . Isobaric yield ratio difference and Shannon information entropy . Phys Lett B , 2015 , 742 : 19 - 22 . DOI: 10.1016/j.physletb.2015.01.015 http://doi.org/10.1016/j.physletb.2015.01.015 Isobaric yield ratio difference and Shannon information entropy W Lin , R Wada , M Huang , et al . Experimental reconstruction of primary fragments with kinematical focusing method . Nucl Sci Tech , 2013 , 24 : 050511 . Experimental reconstruction of primary fragments with kinematical focusing method R Wada , M Hang , W Lin , et al . IMF production and symmetry energy in heavy ion collisions near Fermi energy . Nucl Sci Tech , 2013 , 24 : 050501 . IMF production and symmetry energy in heavy ion collisions near Fermi energy C Ma , C Qiao , S Wang , et al . Temperature and symmetry energy of neutron-rich fragments in the 1A GeV ^124,136Xe+Pb reactions . Nucl Sci Tech , 2013 , 24 : 050510 . Temperature and symmetry energy of neutron-rich fragments in the 1A GeV ^124,136Xe+Pb reactions J Su and F S Zhang . Isotopic dependence of nuclear temperatures . Phys Rev C , 2011 , 84 : 037601 . DOI: 10.1103/PhysRevC.84.037601 http://doi.org/10.1103/PhysRevC.84.037601 Isotopic dependence of nuclear temperatures J Su , L Zhu , W J Xie , et al . Nuclear temperatures from kinetic characteristics . Phys Rev C , 2012 , 85 : 017604 . DOI: 10.1103/PhysRevC.85.017604 http://doi.org/10.1103/PhysRevC.85.017604 Nuclear temperatures from kinetic characteristics C W Ma , J Pu , Y G Ma , et al . Temperature determined by isobaric yield ratios in heavy-ion collisions . Phys Rev C , 2012 , 86 : 054611 . DOI: 10.1103/PhysRevC.86.054611 http://doi.org/10.1103/PhysRevC.86.054611 Temperature determined by isobaric yield ratios in heavy-ion collisions C W Ma , X L Zhao , J Pu , et al . Temperature determined by isobaric yield ratios in a grand-canonical ensemble theory . Phys Rev C , 2013 , 88 : 014609 . DOI: 10.1103/PhysRevC.88.014609 http://doi.org/10.1103/PhysRevC.88.014609 Temperature determined by isobaric yield ratios in a grand-canonical ensemble theory C W Ma , S S Wang , J Pu , et al . Temperature of heavy fragments in heavy-ion collisions . Commun Theor Phys , 2013 , 59 : 95 - 98 . Temperature of heavy fragments in heavy-ion collisions H Y Wu , Z G Xiao , G M Jin , et al . Evidence of “slow” relaxation of isospin degree of freedom . Phys Lett B , 2002 , 538 : 39 - 44 . DOI: 10.1016/S0370-2693(02)01971-8 http://doi.org/10.1016/S0370-2693(02)01971-8 Evidence of “slow” relaxation of isospin degree of freedom C W Ma , Y L Zhang , C Y Qiao , et al . Target effects in isobaric yield ratio differences between projectile fragmentation reactions . Phys Rev C , 2015 , 91 : 014615 . DOI: 10.1103/PhysRevC.91.014615 http://doi.org/10.1103/PhysRevC.91.014615 Target effects in isobaric yield ratio differences between projectile fragmentation reactions 10.13538/j.1001-8042/nst.26.S20503 Supported by the Program for Science & Technology Innovation Talents in Universities of Henan Province (No. 13HASTIT046) and the Young Teacher Project in Henan Normal University 2015-01-02 2015-01-26 2015-04-20 2015-04-20 2021-04-11 2015 Nuclear Science and Techniques 2 26