Isospin splitting of nucleon effective mass and isospin effect in heavy-ion collisions

Isospin splitting of nucleon effective mass and isospin effect in heavy-ion collisions 1 corresp first-author Corresponding author. E-mail address:fengzhq@impcas.ac.cn Corresponding author. E-mail address:fengzhq@impcas.ac.cn fengzhq@impcas.ac.cn Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou 730000, China Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou 730000, China Two different isospin splittings of nucleon effective mass in nuclear medium as the form of <math id="M1"><mrow><msubsup><mi>m</mi><mi>n</mi><mo>*</mo></msubsup><mo>></mo><msubsup><mi>m</mi><mi>p</mi><mo>*</mo></msubsup></mrow></math> and <math id="M2"><mrow><msubsup><mi>m</mi><mi>n</mi><mo>*</mo></msubsup><mo><</mo><msubsup><mi>m</mi><mi>p</mi><mo>*</mo></msubsup></mrow></math> have been implemented in an isospin and momentum dependent transport model. Their impacts on the isospin emission in heavy-ion collisions is investigated thoroughly. It is found that the yield ratios of energetic neutrons to protons squeezed out during the compression stage of two colliding nuclides are sensitive to the isospin splitting. The elliptic flows of free nucleons are also to be promising observables for extracting the nucleon effective mass splitting. Further experimental measurements are being expected, in particular at the CSR-CEE platform in Lanzhou. Several observables are proposed for constraining the density dependence of symmetry energy, such as the transverse flow difference of neutrons and protons, double ratios of n / p and π − / π + , excitation functions of π − / π + and K 0 / K + . Isospin splitting Nucleon effective mass Collective flows Symmetry energy 1 Introduction 1 It has been well known that the masses of neutron and proton are equal each other in vacuum (about 1 GeV), and the effective mass in nuclear matter or finite nuclei deviates from its vacuum value [ 1 1 , 2 2 ] . Moreover, a splitting of neutron and proton effective mass exists in neutron-rich nuclear matter, which increases with the isospin asymmetry and nucleon density. Predictions of the mass splitting based on nuclear many-body theories differ widely. Based on the realistic nucleon-nucleon interactions, the Brueckner-Hartree- Fock (BHF) and Dirac-Brueckner-Hartree-Fock (DBHF) calculations predict a neutron-proton mass splitting of <math id="M3"><mrow><msubsup><mi>m</mi><mi>n</mi><mo>*</mo></msubsup></mrow></math> > <math id="M4"><mrow><msubsup><mi>m</mi><mi>p</mi><mo>*</mo></msubsup></mrow></math> [ 3 3 - 6 6 3-6 ] . The models based on relativistic mean-field (RMF) theory give a contrary result after the inclusion of the isovector mesons ρ and δ [ 7 7 , 8 8 , 9 9 ] . The Skyrme-Hartree-Fock (SHF) predicts both splittings of <math id="M5"><mrow><msubsup><mi>m</mi><mi>n</mi><mo>*</mo></msubsup></mrow></math> > <math id="M6"><mrow><msubsup><mi>m</mi><mi>p</mi><mo>*</mo></msubsup></mrow></math> and <math id="M7"><mrow><msubsup><mi>m</mi><mi>n</mi><mo>*</mo></msubsup></mrow></math> < <math id="M8"><mrow><msubsup><mi>m</mi><mi>p</mi><mo>*</mo></msubsup></mrow></math> exist with different Skyrme parameters [ 10 10 ] . More realistically, the mass splitting of neutron and proton has been constrained from the energy dependence of the Lane (symmetry) potential by the nucleon-nucleus scattering experimental data, which moderately supports the mass splitting of <math id="M9"><mrow><msubsup><mi>m</mi><mi>n</mi><mo>*</mo></msubsup></mrow></math> > <math id="M10"><mrow><msubsup><mi>m</mi><mi>p</mi><mo>*</mo></msubsup></mrow></math> [ 11 11 ] . Further constraints on the effective mass splitting from heavy-ion collisions are still necessary and its influence on reaction dynamics would be very interesting. The knowledge of the density dependence of the nuclear symmetry energy is still poorly known, in particular at high baryon densities, but which has important application in nuclear physics itself and also in astrophysics. Heavy-ion reactions with neutron-rich beams provide a unique opportunity to explore the density-dependent symmetry energy in a broad domain of density. In this work, we present systematic investigations of the effective mass splitting of neutron and proton and its influence on reaction dynamics. An extraction of the high-density symmetry energy from the preequilibrium nucleon emission and meson production ( π and K ) is performed with an isospin and momentum dependent transport model [ 12 12 ] . 2 Model description 1 The time evolutions of the baryons and mesons in the process of two colliding partners under the self-consistently generated mean-field are governed by Hamilton's equations of motion. The Hamiltonian of baryons consists of the relativistic energy, the effective interaction potential and the momentum dependent part [ 12 12 ] . A Skyrme-type momentum-dependent nucleon -nucleon force distinguishing isospin effect is parameterized and implemented in the Lanzhou Quantum Molecular Dynamics (LQMD) model, which leads to a splitting of nucleon effective mass in neutron-rich matter. A Skyrme-type form for the momentum-dependent potential in the Hamiltonian by distinguishing isospin effect is taken in the LQMD model, which is expressed as <math display="block" id="M11"><mtable><mtr><mtd><msub><mi>U</mi><mrow><mtext>mom</mtext></mrow></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mn>2</mn><msub><mi>ρ</mi><mn>0</mn></msub></mrow></mfrac><mstyle displaystyle="true"><munder><mo>∑</mo><mrow><mi>τ</mi><msup><mi>τ</mi><mo>′</mo></msup></mrow></munder><mrow><msub><mi>C</mi><mrow><mi>τ</mi><msup><mi>τ</mi><mo>′</mo></msup></mrow></msub><mstyle displaystyle="true"><mrow><mo>∭</mo><mrow><mi>d</mi><mi>p</mi><mi>d</mi><msup><mi>p</mi><mo>′</mo></msup><mi>d</mi><mi>r</mi><msub><mi>f</mi><mi>τ</mi></msub><mo stretchy="false">(</mo><mi>r</mi><mo>,</mo><mi>p</mi><mo stretchy="false">)</mo></mrow></mrow></mstyle></mrow></mstyle></mtd></mtr><mtr><mtd><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mo>×</mo><msup><mrow><mo>[</mo><mrow><mi>ln</mi><mo stretchy="false">(</mo><mi>ε</mi><msup><mrow><mo stretchy="false">(</mo><mi>p</mi><mo>−</mo><msup><mi>p</mi><mo>′</mo></msup><mo stretchy="false">)</mo></mrow><mn>2</mn></msup><mo>+</mo><mn>1</mn><mo stretchy="false">)</mo></mrow><mo>]</mo></mrow><mn>2</mn></msup><msub><mi>f</mi><msup><mi>τ</mi><mo>′</mo></msup></msub><mo stretchy="false">(</mo><mi>r</mi><mo>,</mo><msup><mi>p</mi><mo>′</mo></msup><mo stretchy="false">)</mo><mo>.</mo></mtd></mtr></mtable></math> Here <math id="M12"><mrow><msub><mi>C</mi><mrow><mi>τ</mi><mi>τ</mi></mrow></msub><mo>=</mo><msub><mi>C</mi><mrow><mi mathvariant="normal">mom</mi></mrow></msub><mo stretchy="false">(</mo><mn>1</mn><mo>+</mo><mi>x</mi><mo stretchy="false">)</mo><msup><mo>,</mo><mrow/></msup><msub><mi>C</mi><mrow><mi>τ</mi><msup><mi>τ</mi><mo>′</mo></msup></mrow></msub><mo>=</mo><msub><mi>C</mi><mrow><mi mathvariant="normal">mom</mi></mrow></msub><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><mi>x</mi><mo stretchy="false">)</mo></mrow></math> ( τ ≠ τ ') and the isospin symbols τ ( τ ') represent proton or neutron. The sign of x determines different mass splitting of proton and neutron in nuclear medium, e.g. positive signs corresponding to the case of <math id="M13"><mrow><msubsup><mi>m</mi><mi>n</mi><mo>*</mo></msubsup></mrow></math> < <math id="M14"><mrow><msubsup><mi>m</mi><mi>p</mi><mo>*</mo></msubsup></mrow></math> . The parameters C mom and ε were determined by fitting the real part of optical potential as a function of incident energy from the proton-nucleus elastic scattering data. In the calculation, we take the values of 1.76 MeV, 500 c 2 /GeV 2 for C mom and ε , respectively, which result in the effective mass m * / m =0.75 in nuclear medium at saturation density for symmetric nuclear matter. The parameter x can be adjusted as the strength of the mass splitting, and the values of −0.65 and 0.65 are respective to the cases of <math id="M15"><mrow><msubsup><mi>m</mi><mi>n</mi><mo>*</mo></msubsup></mrow></math> > <math id="M16"><mrow><msubsup><mi>m</mi><mi>p</mi><mo>*</mo></msubsup></mrow></math> and <math id="M17"><mrow><msubsup><mi>m</mi><mi>n</mi><mo>*</mo></msubsup></mrow></math> < <math id="M18"><mrow><msubsup><mi>m</mi><mi>p</mi><mo>*</mo></msubsup></mrow></math> . One can get the contribution of momentum-dependent interaction to symmetry energy for a cold nuclear with Fermi distribution for a phase-space density <math id="M19"><mrow><msub><mi>f</mi><mi>τ</mi></msub><mo stretchy="false">(</mo><mi>r</mi><mo>,</mo><mi>p</mi><mo stretchy="false">)</mo><mo>=</mo><msub><mi>ρ</mi><mi>τ</mi></msub><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo><mi>Θ</mi><mo stretchy="false">(</mo><msub><mi>p</mi><mi>F</mi></msub><mo stretchy="false">(</mo><mi>τ</mi><mo stretchy="false">)</mo><mo>−</mo><mrow><mo>|</mo><mi>p</mi><mo>|</mo></mrow><mo stretchy="false">)</mo><mo>/</mo><mo stretchy="false">(</mo><mn>4</mn><mi>π</mi><msubsup><mi>p</mi><mi>F</mi><mn>3</mn></msubsup><mo stretchy="false">(</mo><mi>τ</mi><mo stretchy="false">)</mo><mo>/</mo><mn>3</mn><mo stretchy="false">)</mo></mrow></math> with a Fermi momentum <math id="M20"><mrow><msub><mi>p</mi><mi>F</mi></msub><mo stretchy="false">(</mo><mi>τ</mi><mo stretchy="false">)</mo><mo>=</mo><mi>ℏ</mi><msup><mrow><mo stretchy="false">(</mo><mn>3</mn><msup><mi>π</mi><mn>2</mn></msup><msub><mi>ρ</mi><mi>τ</mi></msub><mo stretchy="false">)</mo></mrow><mrow><mn>1</mn><mo>/</mo><mn>3</mn></mrow></msup><mo>.</mo></mrow></math> The symmetry energy per nucleon in the LQMD model is composed of three parts, namely the kinetic energy, the local part and the momentum dependence of the potential energy as <math display="block" id="M21"><mrow><msub><mi>E</mi><mrow><mtext>sym</mtext></mrow></msub><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mn>1</mn><mn>3</mn></mfrac><mfrac><mrow><msup><mi>ℏ</mi><mn>2</mn></msup></mrow><mrow><mn>2</mn><mi>m</mi></mrow></mfrac><msup><mrow><mrow><mo>(</mo><mrow><mfrac><mn>3</mn><mn>2</mn></mfrac><msup><mi>π</mi><mn>2</mn></msup><mi>ρ</mi></mrow><mo>)</mo></mrow></mrow><mrow><mn>2</mn><mo>/</mo><mn>3</mn></mrow></msup><mo>+</mo><msubsup><mi>E</mi><mrow><mtext>sym</mtext></mrow><mrow><mtext>loc</mtext></mrow></msubsup><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo><mo>+</mo><msubsup><mi>E</mi><mrow><mtext>sym</mtext></mrow><mrow><mtext>mom</mtext></mrow></msubsup><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo><mo>.</mo></mrow></math> After the inclusion of momentum dependent contribution, a density, isospin and momentum dependent single-nucleon potential is obtained as follows: <math display="block" id="M22"><mtable><mtr><mtd><msub><mi>U</mi><mi>τ</mi></msub><mo stretchy="false">(</mo><mi>ρ</mi><mo>,</mo><mi>δ</mi><mo>,</mo><mi>p</mi><mo stretchy="false">)</mo><mo>=</mo><mi>α</mi><mfrac><mi>ρ</mi><mrow><msub><mi>ρ</mi><mn>0</mn></msub></mrow></mfrac><mo>+</mo><mi>β</mi><msup><mrow><mo>(</mo><mrow><mfrac><mi>ρ</mi><mrow><msub><mi>ρ</mi><mn>0</mn></msub></mrow></mfrac></mrow><mo>)</mo></mrow><mi>γ</mi></msup><mo>+</mo><mfrac><mn>8</mn><mn>3</mn></mfrac><msub><mi>g</mi><mi>τ</mi></msub><msup><mrow><mo>(</mo><mrow><mfrac><mi>ρ</mi><mrow><msub><mi>ρ</mi><mn>0</mn></msub></mrow></mfrac></mrow><mo>)</mo></mrow><mrow><mn>5</mn><mo>/</mo><mn>3</mn></mrow></msup><mo>+</mo><msubsup><mi>E</mi><mrow><mi mathvariant="normal">sym</mi></mrow><mrow><mi mathvariant="normal">loc</mi></mrow></msubsup><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo><msup><mi>δ</mi><mn>2</mn></msup></mtd></mtr><mtr><mtd><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mo>+</mo><mfrac><mrow><mo>∂</mo><msubsup><mi>E</mi><mrow><mi mathvariant="normal">sym</mi></mrow><mrow><mi mathvariant="normal">loc</mi></mrow></msubsup><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo></mrow><mrow><mo>∂</mo><mi>ρ</mi></mrow></mfrac><mi>ρ</mi><msup><mi>δ</mi><mn>2</mn></msup><mo>+</mo><msubsup><mi>E</mi><mrow><mi mathvariant="normal">sym</mi></mrow><mrow><mi mathvariant="normal">loc</mi></mrow></msubsup><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo><mi>ρ</mi><mfrac><mrow><mo>∂</mo><msup><mi>δ</mi><mn>2</mn></msup></mrow><mrow><mo>∂</mo><msub><mi>ρ</mi><mi>τ</mi></msub></mrow></mfrac></mtd></mtr><mtr><mtd><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mo>+</mo><mfrac><mn>1</mn><mrow><msub><mi>ρ</mi><mn>0</mn></msub></mrow></mfrac><mi>C</mi><mrow><mo>∫</mo><mrow><mi>d</mi><msup><mi>p</mi><mo>′</mo></msup><msub><mi>f</mi><mi>τ</mi></msub><mo stretchy="false">(</mo><mi>r</mi><mo>,</mo><mi>p</mi><mo stretchy="false">)</mo><msup><mrow><mrow><mo>[</mo><mrow><mi>ln</mi><mo stretchy="false">(</mo><mi>ε</mi><msup><mrow><mo stretchy="false">(</mo><mi>p</mi><mo>−</mo><msup><mi>p</mi><mo>′</mo></msup><mo stretchy="false">)</mo></mrow><mn>2</mn></msup><mo>+</mo><mn>1</mn><mo stretchy="false">)</mo></mrow><mo>]</mo></mrow></mrow><mn>2</mn></msup></mrow></mrow></mtd></mtr><mtr><mtd><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mo>+</mo><mfrac><mn>1</mn><mrow><msub><mi>ρ</mi><mn>0</mn></msub></mrow></mfrac><mi>C</mi><mrow><mo>∫</mo><mrow><mi>d</mi><msup><mi>p</mi><mo>′</mo></msup><msub><mi>f</mi><msup><mi>τ</mi><mo>′</mo></msup></msub><mo stretchy="false">(</mo><mi>r</mi><mo>,</mo><mi>p</mi><mo stretchy="false">)</mo><msup><mrow><mrow><mo>[</mo><mrow><mi>ln</mi><mo stretchy="false">(</mo><mi>ε</mi><msup><mrow><mo stretchy="false">(</mo><mi>p</mi><mo>−</mo><msup><mi>p</mi><mo>′</mo></msup><mo stretchy="false">)</mo></mrow><mn>2</mn></msup><mo>+</mo><mn>1</mn><mo stretchy="false">)</mo></mrow><mo>]</mo></mrow></mrow><mn>2</mn></msup></mrow></mrow></mtd></mtr></mtable></math> Here τ ≠ τ ', ∂δ 2 / ∂ρ n =4 δρ p / ρ 2 , and ∂δ 2 / ∂ρ p = −4 δρ n / ρ 2 . The nucleon effective (Landau) mass in nuclear matter of isospin asymmetry δ =( ρ n − ρ p )/ ρ with ρ n and ρ p being the neutron and proton density, respectively, is calculated through the potential as <math id="M23"><mrow><msubsup><mi>m</mi><mi>τ</mi><mo>∗</mo></msubsup><mo>=</mo><msub><mi>m</mi><mi>τ</mi></msub><mo>/</mo><mrow><mo>(</mo><mrow><mn>1</mn><mo>+</mo><mfrac><mrow><msub><mi>m</mi><mi>τ</mi></msub></mrow><mrow><mrow><mo>|</mo><mi>p</mi><mo>|</mo></mrow></mrow></mfrac><mrow><mo>|</mo><mrow><mfrac><mrow><mi>d</mi><msub><mi>U</mi><mi>τ</mi></msub></mrow><mrow><mi>d</mi><mi>p</mi></mrow></mfrac></mrow><mo>|</mo></mrow></mrow><mo>)</mo></mrow></mrow></math> with the free mass m τ at Fermi momentum p = p F . Therefore, the nucleon effective mass only depends on the momentum- dependent term of the nucleon optical potential. Figure 1 Figure 1 is a comparison of the momentum dependence of single-nucleon optical potential with the mass splittings of <math id="M24"><mrow><msubsup><mi>m</mi><mi>n</mi><mo>*</mo></msubsup></mrow></math> > <math id="M25"><mrow><msubsup><mi>m</mi><mi>p</mi><mo>*</mo></msubsup></mrow></math> in the left window and <math id="M26"><mrow><msubsup><mi>m</mi><mi>n</mi><mo>*</mo></msubsup></mrow></math> < <math id="M27"><mrow><msubsup><mi>m</mi><mi>p</mi><mo>*</mo></msubsup></mrow></math> in the right window for the hard and supersoft symmetry energies, respectively. One should note that a cross of neutron and proton single-nucleon potentials appears in the case of <math id="M28"><mrow><msubsup><mi>m</mi><mi>n</mi><mo>*</mo></msubsup></mrow></math> > <math id="M29"><mrow><msubsup><mi>m</mi><mi>p</mi><mo>*</mo></msubsup></mrow></math> at high momentum, but which does not take place in the mass splitting of <math id="M30"><mrow><msubsup><mi>m</mi><mi>n</mi><mo>*</mo></msubsup></mrow></math> < <math id="M31"><mrow><msubsup><mi>m</mi><mi>p</mi><mo>*</mo></msubsup></mrow></math> and a broader separation exists with increasing the nucleon momentum. The difference of single-particle potentials for both mass splitting affects the isospin emission in heavy-ion collisions, in particular, in the domain of high-baryon densities. The scattering in two-particle collisions is performed by using a Monte Carlo procedure, in which the probability to be a channel in a collision is calculated by its contribution of the channel cross section to the total cross section. The primary products in nucleon-nucleon ( NN ) collisions in the region of 1 A GeV energies are the resonances of Δ(1232), N * (1440), N * (1535) and the pions. We have included the reaction channels as follows: <math display="block" id="M34"><mtable columnalign="left"><mtr><mtd><mi>N</mi><mi>N</mi><mo>↔</mo><mi>N</mi><mi>Δ</mi><mo>,</mo><mtext> </mtext><mi>N</mi><mi>N</mi><mo>↔</mo><mi>N</mi><msup><mi>N</mi><mo>*</mo></msup><mo>,</mo><mtext> </mtext><mi>N</mi><mi>N</mi><mo>↔</mo><mi>Δ</mi><mi>Δ</mi><mo>,</mo></mtd></mtr><mtr><mtd><mi>Δ</mi><mo>↔</mo><mi>N</mi><mi>π</mi><mo>,</mo><mtext> </mtext><msup><mi>N</mi><mo>*</mo></msup><mo>↔</mo><mi>N</mi><mi>π</mi><mo>,</mo><mtext> </mtext><mi>N</mi><mi>N</mi><mo>↔</mo><mi>N</mi><mi>N</mi><mi>π</mi><mo stretchy="false">(</mo><mi>s</mi><mo>−</mo><mi mathvariant="normal">state</mi><mo stretchy="false">)</mo><mo>,</mo></mtd></mtr><mtr><mtd><msup><mi>N</mi><mo>*</mo></msup><mo stretchy="false">(</mo><mn>1535</mn><mo stretchy="false">)</mo><mo>→</mo><mi>N</mi><mi>η</mi><mo>.</mo></mtd></mtr></mtable></math> At the considered energies, there are mostly Δ resonances which disintegrate into a π and a nucleon in the evolutions. However, the N * yet gives considerable contribution to the energetic pion yields. The strangeness is created as the secondary products in inelastic hadron-hadron collisions [ 13 13 , 14 14 ] . We included the channels as follows: <math display="block" id="M35"><mtable columnalign="left"><mtr><mtd><mi>B</mi><mi>B</mi><mo>→</mo><mi>B</mi><mi>Y</mi><mi>K</mi><mo>,</mo><mtext> </mtext><mi>B</mi><mi>B</mi><mo>→</mo><mi>B</mi><mi>B</mi><mi>K</mi><mover accent="true"><mi>K</mi><mo>¯</mo></mover><mo>,</mo><mtext> </mtext><mi>B</mi><mi>π</mi><mo>→</mo><mi>Y</mi><mi>K</mi><mo>,</mo></mtd></mtr><mtr><mtd><mi>B</mi><mi>π</mi><mo>→</mo><mi>N</mi><mi>K</mi><mover accent="true"><mi>K</mi><mo>¯</mo></mover><mo>,</mo><mtext> </mtext><mi>Y</mi><mi>π</mi><mo>→</mo><mi>B</mi><mover accent="true"><mi>K</mi><mo>¯</mo></mover><mo>,</mo><mtext> </mtext><mi>B</mi><mover accent="true"><mi>K</mi><mo>¯</mo></mover><mo>→</mo><mi>Y</mi><mi>π</mi><mo>,</mo></mtd></mtr><mtr><mtd><mi>Y</mi><mi>N</mi><mo>→</mo><mover accent="true"><mi>K</mi><mo>¯</mo></mover><mi>N</mi><mi>N</mi><mo>.</mo></mtd></mtr></mtable></math> Here the B strands for ( N , Δ, N * ) and Y(Λ, ∑), K ( K 0 , K + ) and <math id="M101"><mover accent="true"><mi>K</mi><mo>¯</mo></mover></math> ( K 0 , K − ). The elastic scattering between strangeness and baryons is considered through the channels of K B → K B , Y B → Y B and <math id="M102"><mover accent="true"><mi>K</mi><mo>¯</mo></mover></math> B → <math id="M103"><mover accent="true"><mi>K</mi><mo>¯</mo></mover></math> B . The charge-exchange reactions between the K N → K N and Y N → Y N channels are included by using the same cross sections with the elastic scattering, such as K 0 p → K + n , K + n → K 0 p etc. Correction of effective mass of kaons in nuclear medium on the elementary cross section is considered through the threshold energy, which results in the reduction of kaon and the enhancement of anti-kaon yields in heavy-ion collisions. The evolution of mesons (here mainly pions and kaons) is also determined by a Hamiltonian. Here, the Coulomb interaction is considered for pion and kaon propagation. The kaon and anti-kaon energy in the nuclear medium is computed from the chiral Lagrangian. A repulsive kaon-nucleon ( KN ) potential and attractive antikaon-nucleon potential with the values of 25.5 MeV and −96.8 MeV are used in the model, which result in the effective mass with <math id="M36"><mrow><msubsup><mi>m</mi><mi>K</mi><mo>*</mo></msubsup></mrow></math> / <math id="M37"><mrow><msub><mi>m</mi><mi>K</mi></msub></mrow></math> =1.05 and <math id="M38"><mrow><msubsup><mi>m</mi><mover accent="true"><mi>K</mi><mo>¯</mo></mover><mo>*</mo></msubsup></mrow></math> / <math id="M39"><mrow><msub><mi>m</mi><mover accent="true"><mi>K</mi><mo>¯</mo></mover></msub></mrow></math> =0.8 at saturation baryon density. 3 Results and discussion 1 The preequilibrium nucleons in high-energy heavy-ion collisions are mostly produced during a compression stage of two colliding partners within a very short time. Therefore, the high-density information of nuclear phase diagram is expected to be extracted from the nucleons or light complex particles, which can be constrained from the longitudinal rapidity distributions and the azimuthal emissions. Shown in Fig.2 Fig.2 is a comparison of transverse emission ratios of neutron/proton within the rapidity selection of | y / y proj | < 0.25 in the 197 Au+ 197 Au reaction at the incident energy of 400 A MeV for the near central ( b =1 fm) and semi-central ( b =6 fm) collisions, but with different mass splittings. One can see that the influence of the symmetry energy appears at low transverse momentum and the situation does not be changed with the mass splitting. A hard symmetry energy enforces a strongly repulsive force on neutrons in the high-density domain, furthermore, squeezes out more neutrons in the preequilibrium stage of dynamical evolution. From the negative contribution of the momentum-dependent interaction of <math id="M40"><mrow><msubsup><mi>m</mi><mi>n</mi><mo>*</mo></msubsup></mrow></math> > <math id="M41"><mrow><msubsup><mi>m</mi><mi>p</mi><mo>*</mo></msubsup></mrow></math> to the potential part of the symmetry energy in Ref.[ 12 12 ], one notices that the momentum-dependent potential leads to an attractive force on neutrons, in particular at high densities, which reduces the n / p yields at high transverse momentum. Opposite contributions take place and a flat distribution appears in the case of <math id="M42"><mrow><msubsup><mi>m</mi><mi>n</mi><mo>*</mo></msubsup></mrow></math> < <math id="M43"><mrow><msubsup><mi>m</mi><mi>p</mi><mo>*</mo></msubsup></mrow></math> . The spectrum of n / p has a slightly different structure with the results in Ref.[ 15 15 ] that is caused from the momentum dependent contribution of symmetry potential. Collective flow has been verified as a nice approach to reconstruct the reaction plane and to study azimuthal correlation of the fireball formed in heavy-ion collisions. More sensitive observable can be seen from the flow difference between neutron and proton as shown in Fig.3 Fig.3 . The <math id="M44"><mrow><msubsup><mi>m</mi><mi>n</mi><mo>*</mo></msubsup></mrow></math> > <math id="M45"><mrow><msubsup><mi>m</mi><mi>p</mi><mo>*</mo></msubsup></mrow></math> case is always larger, in particular in the domain of high transverse momentum. Basically, the symmetry energy does not change the spectrum structure. We calculated the difference of neutron and proton directed flows as a function of rapidity distribution in the 197 Au+ 197 Au collisions ( b =6 fm) at the energy of 400 A MeV with the hard and supersoft symmetry energies, but for the different mass splitting as shown in Fig.4 Fig.4 . One can see that an anti-flow spectrum appears in both mass splitting in comparison with nucleon flows [ 16 16 , 17 17 ] . Moreover, the supersoft symmetry energy gives a larger flow difference, in particular around the region of projectile (target) rapidity. The phenomena can not be influenced by changing different mass splitting. The flows of fast nucleon emissions would be a nice probe for extracting the high-density behavior of nuclear symmetry energy. Kaon meson is produced at the early stage in heavy-ion collision and promptly emitted after production, which can get directly the information of high-density phase diagram. A pronounced effect of the stiffness of symmetry energy on kaon production can be observed from the spectrum of isospin ratio. The isospin effects appear at deep subthreshold energies as shown in Fig.5 Fig.5 . The in-medium potential slightly changes the K 0 / K + value because of its influence on the kaon propagation and also on the charge-exchange reactions. At the considered energies, the channel of N Δ→ NYK contributes the main part for the kaon yields due to the larger production cross sections and the higher invariant energy, and the NN → NYK as well as πN → YK have about one third contributions. One notices that a hard symmetry energy always has the larger values of the isospin ratios than the supersoft case in the domain of subthreshold energies ( E th ( K )=1.58 GeV). 4 Conclusion 1 In summary, within the transport model (LQMD) we have investigated the impact of the momentum- dependent interaction on fast nucleon emissions in heavy-ion collisions. Two different mass splittings, i.e., <math id="M48"><mrow><msubsup><mi>m</mi><mi>n</mi><mo>*</mo></msubsup><mo>></mo><msubsup><mi>m</mi><mi>p</mi><mo>*</mo></msubsup></mrow></math> and <math id="M49"><mrow><msubsup><mi>m</mi><mi>n</mi><mo>*</mo></msubsup><mo><</mo><msubsup><mi>m</mi><mi>p</mi><mo>*</mo></msubsup></mrow></math> , are chosen in the calculations. The momentum-dependent potential plays a significant role on the fast nucleon emissions. Specifically, the neutron/proton ratio at high transverse momenta is influenced, in which a flat spectrum appears in the case of <math id="M50"><mrow><msubsup><mi>m</mi><mi>n</mi><mo>*</mo></msubsup><mo><</mo><msubsup><mi>m</mi><mi>p</mi><mo>*</mo></msubsup></mrow></math> . The elliptic flow difference between neutron and proton sensitively depends on the mass splitting, in particular at high transverse momenta. The neutron/proton ratio and the transverse flow difference of neutron and proton via transverse momentum distribution in the domain of mid-rapidity are sensitive to the stiffness of nuclear symmetry energy, which can be promising observables for extracting the high-density behavior of symmetry energy. The K 0 / K + ratio of neutron-rich heavy system in the domain of subthreshold energies is also sensitive to the stiffness of nuclear symmetry energy. References Jeukenne J P , Lejeune A , Mahaux C . Phys Rep , 1976 , 25 : 83 - 174 . Phys Rep Mahaux C , Bortignon P F , Broglia R A , et al . Phys Rep , 1985 , 120 : 1 - 274 . Phys Rep van Dalen E N E , Fuchs C , Faessler A . Phys Rev Lett , 2005 , 95 : 022302 . Phys Rev Lett van Dalen E N E , Fuchs C , Faessler A . Phys Rev C , 2005 , 72 : 065803 . Phys Rev C Sammarruca F , Barredo W , Krastev P . Phys Rev C , 2005 , 71 : 064306 . Phys Rev C Zuo W , Bombaci I , Lombardo U . Phys Rev C , 1999 , 60 : 024605 . Phys Rev C Frick T , Gad K , Müther H , et al . Phys Rev C , 2002 , 65 : 034321 . Phys Rev C Hassaneen Kh S A and Müther H . Phys Rev C , 2004 , 70 : 054308 . Phys Rev C Liu B , Greco V , Baran V , et al . Phys Rev C , 2002 , 65 : 045201 . Phys Rev C Baran V , Colonna M , Greco V , et al . Phys Rep , 2005 , 410 : 335 - 466 . Phys Rep Li B A . Phys Rev C , 2004 , 69 : 064602 . Phys Rev C Feng Z Q . Phys Rev C , 2011 , 84 : 024610 . Phys Rev C Feng Z Q and Jin G M . Phys Rev C , 2010 , 82 : 044615 . Phys Rev C Feng Z Q . Phys Rev C , 2011 , 83 : 067604 . Phys Rev C Li Q , Li Z , Soff S , et al . Phys Rev C , 2005 , 72 : 034613 . Phys Rev C Feng Z Q . Phys Lett B , 2012 , 707 : 83 - 87 . Phys Lett B Feng Z Q . Phys Rev C , 2012 , 85 : 014604 . Phys Rev C 10.13538/j.1001-8042/nst.2013.05.004 Supported by National Natural Science Foundation of China (NSFC) projects (No. 11175218) and the Advancement Society of Young Innovation of Chinese Academy of Sciences. 2013-06-27 2013-10-01 2021-04-11 2013 Nuclear Science and Techniques 5 24