Nuclear dynamics induced by antiprotons

Nuclear dynamics induced by antiprotons 1 2 corresp first-author Corresponding author, fengzhq@impcas.ac.cn Corresponding author, 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 State Key Laboratory of Theoretical Physics and Kavli Institute for Theoretical Physics China, Chinese Academy of Sciences, Beijing 100190, China State Key Laboratory of Theoretical Physics and Kavli Institute for Theoretical Physics China, Chinese Academy of Sciences, Beijing 100190, China Reaction dynamics in collisions of antiprotons on nuclei is investigated within the Lanzhou quantum molecular dynamics model. The reaction channels of elastic scattering, annihilation, charge exchange and inelastic collisions of antiprotons on nucleons have been included in the model. Dynamics on particle production, in particular pions, kaons, antikaons and hyperons, is investigated in collisions of <math id="M1"><mrow><mover accent="true"><mi>p</mi><mo stretchy="true">¯</mo></mover></mrow></math> on 12 C, 20 Ne, 40 Ca and 181 Ta from low to high incident momenta. It is found that the annihilations of <math id="M2"><mrow><mover accent="true"><mi>p</mi><mo stretchy="true">¯</mo></mover></mrow></math> on nucleons are of importance to the dynamics of particle production in phase space. Hyperons are mainly produced via meson induced reactions on nucleons and strangeness exchange collisions, which lead to delayed emission in antiproton-nucleus collisions. LQMD transport model Antiproton-nucleus collisions Strangeness production I. INTRODUCTION 1 The dynamics of antiproton-nucleus collisions is a complex process, which is associated with the mean-field potentials of antinucleons and produced particles in a nuclear medium and with a number of reaction channels, i.e., the annihilation channels, charge-exchange reaction, elastic and inelastic collisions. A more localized energy deposition is able to be produced in antiproton-nucleus collisions in comparison to heavy-ion collisions due to the annihilations. Searching for cold quark-gluon plasma (QGP) with antiproton beams has been a hot topic both in experiments and in theory calculations over the past several decades. The large yields of strange particles may be produced in antiproton induced reactions, which have an advantage in comparison to proton-nucleus and heavy-ion collisions. The particle production in collisions of antiprotons on nuclei has been investigated by using the intranuclear cascade (INC) model [ 1 1 , 3 3 , 4 4 ] and the kinetic approach [ 5 5 ]. A number of experimental data was nicely explained within these approaches. A self-consistent description of dynamical evolutions and collisions of antiprotons on a nucleus with transport models is still very necessary, in particular the fragmentation in collisions of antiprotons on a nucleus to form hypernuclei. Strangeness production in antiproton induced nuclear reactions has been investigated thoroughly with the Giessen Boltzmann-Uehling-Uhlenbeck (GiBUU) transport model [ 6 6 , 7 7 ] and the Lanzhou quantum molecular dynamics (LQMD) approach [ 8 8 ]. The production of hypernuclei is associated with the reaction channels of hyperons and also hyperon-nucleon (HN) potential. From comparison of kinetic energy or momentum spectra of hyperons to experimental data, the HN potential can be extracted. Also the antinucleon-nucleon potential is able to be constrained from particle production. The dynamical mechanism on strange particle production can be explored from the analysis of reaction channels and comparison to experimental spectra. II. MODEL DESCRIPTION 1 In the Lanzhou quantum molecular dynamics (LQMD) model, the dynamics of resonances ( ∆ (1232), N *(1440), N *(1535) etc), hyperons ( Λ , ∑ , Ξ ) and mesons ( π , K , η , <math id="M3"><mrow><mover accent="true"><mi>K</mi><mo stretchy="true">¯</mo></mover></mrow></math> , ρ , ω ) is described via hadron-hadron collisions, decays of resonances and mean-field potentials in a nuclear medium [ 9 9 - 18 18 9-18 ]. The evolutions of baryons (nucleons, resonances and hyperons), anti-baryons and mesons in the collisions are governed by Hamilton’s equations of motion which read as <math display="block" id="M4"><mrow><msub><mstyle mathvariant="bold" mathsize="normal"><mover accent="true"><mi>p</mi><mo>˙</mo></mover></mstyle><mi>i</mi></msub><mo>=</mo><mo>−</mo><mfrac><mrow><mo>∂</mo><mi>H</mi></mrow><mrow><mo>∂</mo><msub><mstyle mathvariant="bold" mathsize="normal"><mi>r</mi></mstyle><mi>i</mi></msub></mrow></mfrac><mn>,</mn><mtext> </mtext><msub><mstyle mathvariant="bold" mathsize="normal"><mover accent="true"><mi>r</mi><mo>˙</mo></mover></mstyle><mi>i</mi></msub><mo>=</mo><mfrac><mrow><mo>∂</mo><mi>H</mi></mrow><mrow><mo>∂</mo><msub><mstyle mathvariant="bold" mathsize="normal"><mi>p</mi></mstyle><mi>i</mi></msub></mrow></mfrac><mn>.</mn></mrow></math> The Hamiltonian of baryons consists of the relativistic energy, the effective interaction potential and the momentum dependent part as follows <math display="block" id="M5"><mrow><msub><mi>H</mi><mi>B</mi></msub><mo>=</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder></mstyle><msqrt><mrow><msubsup><mi>p</mi><mi>i</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>m</mi><mi>i</mi><mn>2</mn></msubsup></mrow></msqrt><mo>+</mo><msub><mi>U</mi><mrow><mtext>int</mtext></mrow></msub><mo>+</mo><msub><mi>U</mi><mrow><mtext>mom</mtext></mrow></msub><mn>.</mn></mrow></math> Here the p i and m i represent the momentum and the mass of the baryons. The effective interaction potential U int is composed of the Coulomb interaction and the local potential [ 13 13 - 18 18 13-18 ]. A Skyrme-type momentum dependent interaction has been used in the evaluation of the local potential energy for nucleons and resonances. The effect of the momentum dependence of the symmetry potential in heavy-ion collisions was also investigated with the isospin-dependent Boltzmann Uehling Uhlenbeck transport model [ 19 19 - 21 21 19-21 ]. The hyperon mean-field potential is constructed on the basis of the light-quark counting rule. The self-energies of Λ and ∑ are assumed to be two thirds of that experienced by nucleons. And the Ξ self-energy is one third of a nucleon’s. Thus, the in-medium dispersion relation reads <math display="block" id="M6"><mrow><msub><mi>ω</mi><mi>B</mi></msub><mo stretchy="false">(</mo><msub><mi>p</mi><mi>i</mi></msub><mn>,</mn><msub><mi>ρ</mi><mi>i</mi></msub><mo stretchy="false">)</mo><mo>=</mo><msqrt><mrow><msup><mrow><mo stretchy="false">(</mo><msub><mi>m</mi><mi>B</mi></msub><mo>+</mo><msubsup><mi>Σ</mi><mi>S</mi><mi>B</mi></msubsup><mo stretchy="false">)</mo></mrow><mn>2</mn></msup><mo>+</mo><msubsup><mi>p</mi><mi>i</mi><mn>2</mn></msubsup></mrow></msqrt><mo>+</mo><msubsup><mi>Σ</mi><mi>V</mi><mi>B</mi></msubsup><mn>,</mn></mrow></math> e.g., for hyperons <math id="M7"><mrow><msubsup><mi>Σ</mi><mi>S</mi><mrow><mi>Λ</mi><mn>,</mn><mi>Σ</mi></mrow></msubsup><mo>=</mo><mn>2</mn><msubsup><mi>Σ</mi><mi>S</mi><mi>N</mi></msubsup><mo>/</mo><mn>3</mn></mrow></math> , <math id="M8"><mrow><msubsup><mi>Σ</mi><mi>V</mi><mrow><mi>Λ</mi><mn>,</mn><mi>Σ</mi></mrow></msubsup><mo>=</mo><mn>2</mn><msubsup><mi>Σ</mi><mi>V</mi><mi>N</mi></msubsup><mo>/</mo><mn>3</mn></mrow></math> , <math id="M9"><mrow><msubsup><mi>Σ</mi><mi>S</mi><mi>Ξ</mi></msubsup><mo>=</mo><msubsup><mi>Σ</mi><mi>S</mi><mi>N</mi></msubsup><mo>/</mo><mn>3</mn></mrow></math> and <math id="M10"><mrow><msubsup><mi>Σ</mi><mi>V</mi><mi>Ξ</mi></msubsup><mo>=</mo><msubsup><mi>Σ</mi><mi>V</mi><mi>N</mi></msubsup><mo>/</mo><mn>3</mn></mrow></math> . The antibaryon energy is computed from the G-parity transformation of baryon potential as <math display="block" id="M11"><mrow><msub><mi>ω</mi><mrow><mover accent="true"><mi>B</mi><mo stretchy="true">¯</mo></mover></mrow></msub><mo stretchy="false">(</mo><msub><mi>p</mi><mi>i</mi></msub><mn>,</mn><msub><mi>ρ</mi><mi>i</mi></msub><mo stretchy="false">)</mo><mo>=</mo><msqrt><mrow><msup><mrow><mo stretchy="false">(</mo><msub><mi>m</mi><mrow><mover accent="true"><mi>B</mi><mo stretchy="true">¯</mo></mover></mrow></msub><mo>+</mo><msubsup><mi>Σ</mi><mi>S</mi><mrow><mover accent="true"><mi>B</mi><mo stretchy="true">¯</mo></mover></mrow></msubsup><mo stretchy="false">)</mo></mrow><mn>2</mn></msup><mo>+</mo><msubsup><mi>p</mi><mi>i</mi><mn>2</mn></msubsup></mrow></msqrt><mo>+</mo><msubsup><mi>Σ</mi><mi>V</mi><mrow><mover accent="true"><mi>B</mi><mo stretchy="true">¯</mo></mover></mrow></msubsup><mn>,</mn></mrow></math> with <math id="M12"><mrow><msubsup><mi>Σ</mi><mi>S</mi><mrow><mover accent="true"><mi>B</mi><mo stretchy="true">¯</mo></mover></mrow></msubsup><mo>=</mo><msubsup><mi>Σ</mi><mi>S</mi><mi>B</mi></msubsup></mrow></math> and <math id="M13"><mrow><msubsup><mi>Σ</mi><mi>V</mi><mrow><mover accent="true"><mi>B</mi><mo stretchy="true">¯</mo></mover></mrow></msubsup><mo>=</mo><mo>−</mo><msubsup><mi>Σ</mi><mi>V</mi><mi>B</mi></msubsup></mrow></math> . The nuclear scalar <math id="M14"><mrow><msubsup><mi>Σ</mi><mi>S</mi><mi>N</mi></msubsup></mrow></math> and vector <math id="M15"><mrow><msubsup><mi>Σ</mi><mi>V</mi><mi>N</mi></msubsup></mrow></math> self-energies are computed from the well-known relativistic mean-field model with the NL3 parameter ( <math id="M16"><mrow><msubsup><mi>g</mi><mrow><mi>σ</mi><mi>N</mi></mrow><mn>2</mn></msubsup></mrow></math> =80.82, <math id="M17"><mrow><msubsup><mi>g</mi><mrow><mi>ω</mi><mi>N</mi></mrow><mn>2</mn></msubsup></mrow></math> =155). The optical potential of baryon or antibaryon is derived from the in-medium energy as <math id="M18"><mrow><msub><mi>V</mi><mrow><mi>o</mi><mi>p</mi><mi>t</mi></mrow></msub><mo stretchy="false">(</mo><mi>p</mi><mn>,</mn><mi>ρ</mi><mo stretchy="false">)</mo><mo>=</mo><mi>ω</mi><mo stretchy="false">(</mo><mi>p</mi><mn>,</mn><mi>ρ</mi><mo stretchy="false">)</mo><mo>−</mo><msqrt><mrow><msup><mi>p</mi><mn>2</mn></msup><mo>+</mo><msubsup><mi>m</mi><mi>B</mi><mn>2</mn></msubsup></mrow></msqrt></mrow></math> . A factor ξ is introduced in the evaluation of antinucleon optical potential to mimic the antiproton-nucleus scattering [ 22 22 ] and the real part of phenomenological antinucleon-nucleon optical potential [ 23 23 ] as <math id="M19"><mrow><msubsup><mi>Σ</mi><mi>S</mi><mrow><mover accent="true"><mi>N</mi><mo stretchy="true">¯</mo></mover></mrow></msubsup><mo>=</mo><mi>ξ</mi><msubsup><mi>Σ</mi><mi>S</mi><mi>N</mi></msubsup></mrow></math> and <math id="M20"><mrow><msubsup><mi>Σ</mi><mi>V</mi><mrow><mover accent="true"><mi>N</mi><mo stretchy="true">¯</mo></mover></mrow></msubsup><mo>=</mo><mo>−</mo><mi>ξ</mi><msubsup><mi>Σ</mi><mi>V</mi><mi>N</mi></msubsup></mrow></math> with ξ =0.25, which leads to the optical potential <math id="M21"><mrow><msub><mi>V</mi><mrow><mover accent="true"><mi>N</mi><mo stretchy="true">¯</mo></mover></mrow></msub></mrow></math> =–164 MeV for an antinucleon at the zero momentum and normal nuclear density ρ 0 =0.16 fm -3 . Fig. 1 Fig. 1 shows a comparison of the baryon and antibaryon energies in a nuclear medium and the optical potentials as a function of baryon density. The empirical Λ potential extracted from hypernuclei experiments [ 24 24 , 25 25 ] is well reproduced with the approach. The antihyperons exhibit strongly attractive potentials in a nuclear medium. The optical potentials of hyperons will affect the production of hypernuclei. The values of optical potentials at saturation density are -32 MeV, -16 MeV, -164 MeV, -436 MeV and -218 MeV for Λ ( ), Ξ , <math id="M22"><mrow><mover accent="true"><mi>N</mi><mo stretchy="true">¯</mo></mover></mrow></math> , <math id="M23"><mrow><mover accent="true"><mi>Λ</mi><mo stretchy="true">¯</mo></mover></mrow></math> and <math id="M24"><mrow><mover accent="true"><mi>Ξ</mi><mo stretchy="true">¯</mo></mover></mrow></math> , respectively. The evolution of mesons is also determined by the Hamiltonian, which is given by <math display="block" id="M25"><mrow><msub><mi>H</mi><mi>M</mi></msub><mo>=</mo><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><msub><mi>N</mi><mi>M</mi></msub></mrow></munderover></mstyle><mrow><mo>(</mo><mrow><msubsup><mi>V</mi><mi>i</mi><mrow><mi>C</mi><mi>o</mi><mi>u</mi><mi>l</mi></mrow></msubsup><mo>+</mo><mi>ω</mi><mo stretchy="false">(</mo><msub><mi>p</mi><mi>i</mi></msub><mn>,</mn><msub><mi>ρ</mi><mi>i</mi></msub><mo stretchy="false">)</mo></mrow><mo>)</mo></mrow><mn>.</mn></mrow></math> Here the Coulomb interaction is given by <math display="block" id="M26"><mrow><msubsup><mi>V</mi><mi>i</mi><mrow><mi>C</mi><mi>o</mi><mi>u</mi><mi>l</mi></mrow></msubsup><mo>=</mo><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mrow><msub><mi>N</mi><mi>B</mi></msub></mrow></munderover></mstyle><mfrac><mrow><msub><mi>e</mi><mi>i</mi></msub><msub><mi>e</mi><mi>j</mi></msub></mrow><mrow><msub><mi>r</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub></mrow></mfrac><mn>,</mn></mrow></math> where the N M and N B are the total number of mesons and baryons, including charged resonances. The kaon and anti-kaon energies in the nuclear medium distinguish isospin effects based on the chiral Lagrangian approach in [ 26 26 - 30 30 26-30 ] <math display="block" id="M27"><mrow><msub><mi>ω</mi><mi>K</mi></msub><mo stretchy="false">(</mo><msub><mi>p</mi><mi>i</mi></msub><mn>,</mn><msub><mi>ρ</mi><mi>i</mi></msub><mo stretchy="false">)</mo><mo>=</mo><msup><mrow><mrow><mo>[</mo><mrow><msubsup><mi>m</mi><mi>K</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>p</mi><mi>i</mi><mn>2</mn></msubsup><mo>−</mo><msub><mi>a</mi><mi>K</mi></msub><msubsup><mi>ρ</mi><mi>i</mi><mi>S</mi></msubsup><mo>−</mo><msub><mi>τ</mi><mn>3</mn></msub><msub><mi>c</mi><mi>K</mi></msub><msubsup><mi>ρ</mi><mrow><mi>i</mi><mn>3</mn></mrow><mi>S</mi></msubsup><mo>+</mo><msup><mrow><mo stretchy="false">(</mo><msub><mi>b</mi><mi>K</mi></msub><msub><mi>ρ</mi><mi>i</mi></msub><mo>+</mo><msub><mi>τ</mi><mn>3</mn></msub><msub><mi>d</mi><mi>K</mi></msub><msub><mi>ρ</mi><mrow><mi>i</mi><mn>3</mn></mrow></msub><mo stretchy="false">)</mo></mrow><mn>2</mn></msup></mrow><mo>]</mo></mrow></mrow><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup><mo>+</mo><msub><mi>b</mi><mi>K</mi></msub><msub><mi>ρ</mi><mi>i</mi></msub><mo>+</mo><msub><mi>τ</mi><mn>3</mn></msub><msub><mi>d</mi><mi>K</mi></msub><msub><mi>ρ</mi><mrow><mi>i</mi><mn>3</mn></mrow></msub><mn>,</mn></mrow></math> and <math display="block" id="M28"><mrow><msub><mi>ω</mi><mrow><mover accent="true"><mi>K</mi><mo stretchy="true">¯</mo></mover></mrow></msub><mo stretchy="false">(</mo><msub><mi>p</mi><mi>i</mi></msub><mn>,</mn><msub><mi>ρ</mi><mi>i</mi></msub><mo stretchy="false">)</mo><mo>=</mo><msup><mrow><mrow><mo>[</mo><mrow><msubsup><mi>m</mi><mrow><mover accent="true"><mi>K</mi><mo stretchy="true">¯</mo></mover></mrow><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>p</mi><mi>i</mi><mn>2</mn></msubsup><mo>−</mo><msub><mi>a</mi><mrow><mover accent="true"><mi>K</mi><mo stretchy="true">¯</mo></mover></mrow></msub><msubsup><mi>ρ</mi><mi>i</mi><mi>S</mi></msubsup><mo>−</mo><msub><mi>τ</mi><mn>3</mn></msub><msub><mi>c</mi><mi>K</mi></msub><msubsup><mi>ρ</mi><mrow><mi>i</mi><mn>3</mn></mrow><mi>S</mi></msubsup><mo>+</mo><msup><mrow><mo stretchy="false">(</mo><msub><mi>b</mi><mi>K</mi></msub><msub><mi>ρ</mi><mi>i</mi></msub><mo>+</mo><msub><mi>τ</mi><mn>3</mn></msub><msub><mi>d</mi><mi>K</mi></msub><msub><mi>ρ</mi><mrow><mi>i</mi><mn>3</mn></mrow></msub><mo stretchy="false">)</mo></mrow><mn>2</mn></msup></mrow><mo>]</mo></mrow></mrow><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup><mo>−</mo><msub><mi>b</mi><mi>K</mi></msub><msub><mi>ρ</mi><mi>i</mi></msub><mo>−</mo><msub><mi>τ</mi><mn>3</mn></msub><msub><mi>d</mi><mi>K</mi></msub><msub><mi>ρ</mi><mrow><mi>i</mi><mn>3</mn></mrow></msub><mn>,</mn></mrow></math> respectively. Here the <math id="M29"><mrow><msub><mi>b</mi><mi>K</mi></msub><mo>=</mo><mn>3</mn><mo>/</mo><mo stretchy="false">(</mo><mn>8</mn><msubsup><mi>f</mi><mi>π</mi><mrow><mo>∗</mo><mn>2</mn></mrow></msubsup><mo stretchy="false">)</mo><mo>≈</mo></mrow></math> 0.333 GeV fm 3 , the a K and <math id="M30"><mrow><msub><mi>a</mi><mrow><mover accent="true"><mi>K</mi><mo stretchy="true">¯</mo></mover></mrow></msub></mrow></math> , are 0.18 GeV 2 fm 3 and 0.31 GeV 2 fm 3 , respectively, which results in the strength of repulsive kaon-nucleon ( K N) potential and of attractive antikaon-nucleon <math id="M31"><mrow><mover accent="true"><mi>K</mi><mo stretchy="true">¯</mo></mover></mrow></math> N potential with values of 27.8 MeV and –100.3 MeV at saturation baryon density for isospin symmetric matter, respectively. The τ 3 =1 and -1 for the isospin pair K + ( <math id="M32"><mrow><msup><mrow><mover accent="true"><mi>K</mi><mo stretchy="true">¯</mo></mover></mrow><mn>0</mn></msup></mrow></math> ) and K 0 (K - ), respectively. The parameters c K =0.0298 GeV 2 fm 3 and d K = 0.111 GeV fm 3 determine the isospin splitting of kaons in neutron-rich nuclear matter. The optical potential of the kaon is derived from the in-medium energy as <math id="M33"><mrow><msub><mi>V</mi><mrow><mi>o</mi><mi>p</mi><mi>t</mi></mrow></msub><mo stretchy="false">(</mo><mi>p</mi><mn>,</mn><mi>ρ</mi><mo stretchy="false">)</mo><mo>=</mo><mi>ω</mi><mo stretchy="false">(</mo><mi>p</mi><mn>,</mn><mi>ρ</mi><mo stretchy="false">)</mo><mo>−</mo><msqrt><mrow><msup><mi>p</mi><mn>2</mn></msup><mo>+</mo><msubsup><mi>m</mi><mi>K</mi><mn>2</mn></msubsup></mrow></msqrt></mrow></math> . The values of <math id="M34"><mrow><msubsup><mi>m</mi><mi>K</mi><mo>*</mo></msubsup><mo>/</mo><msub><mi>m</mi><mi>K</mi></msub></mrow></math> =1.056 and <math id="M35"><mrow><msubsup><mi>m</mi><mrow><mover accent="true"><mi>K</mi><mo stretchy="true">¯</mo></mover></mrow><mo>*</mo></msubsup><mo>/</mo><msub><mi>m</mi><mrow><mover accent="true"><mi>K</mi><mo stretchy="true">¯</mo></mover></mrow></msub></mrow></math> =0.797 at normal baryon density concluded with the parameters in isospin symmetric nuclear matter. The effective mass is used to evaluate the threshold energy for kaon and antikaon production, e.g., the threshold energy in the pion-baryon collisions <math id="M36"><mrow><msqrt><mrow><msub><mi>s</mi><mrow><mi>t</mi><mi>h</mi></mrow></msub></mrow></msqrt><mo>=</mo><msubsup><mi>m</mi><mi>Y</mi><mo>∗</mo></msubsup><mo>+</mo><msubsup><mi>m</mi><mi>K</mi><mo>∗</mo></msubsup></mrow></math> . Based on hadron-hadron collisions in heavy-ion reactions in the region of 1–2 A GeV energies [ 13 13 - 18 18 13-18 ], we have further included the annihilation channels, charge-exchange reaction, elastic and inelastic scattering in antinucleon-nucleon collisions: <math id="M37"><mrow><mover accent="true"><mi>N</mi><mo stretchy="true">¯</mo></mover><mi>N</mi><mo>→</mo><mover accent="true"><mi>N</mi><mo stretchy="true">¯</mo></mover><mi>N</mi></mrow></math> , <math id="M38"><mrow><mover accent="true"><mi>N</mi><mo stretchy="true">¯</mo></mover><mi>N</mi><mo>→</mo><mover accent="true"><mi>B</mi><mo stretchy="true">¯</mo></mover><mi>B</mi></mrow></math> , <math id="M39"><mrow><mover accent="true"><mi>N</mi><mo stretchy="true">¯</mo></mover><mi>N</mi><mo>→</mo><mover accent="true"><mi>Y</mi><mo stretchy="true">¯</mo></mover><mi>Y</mi></mrow></math> and <math id="M40"><mrow><mover accent="true"><mi>N</mi><mo stretchy="true">¯</mo></mover><mi>N</mi><mo>→</mo><mi>a</mi><mi>n</mi><mi>n</mi><mi>i</mi><mi>h</mi><mi>i</mi><mi>l</mi><mi>a</mi><mi>t</mi><mi>i</mi><mi>o</mi><mi>n</mi><mo stretchy="false">(</mo><mi>π</mi><mn>,</mn><mi>η</mi><mn>,</mn><mi>ρ</mi><mn>,</mn><mi>ω</mi><mn>,</mn><mi>K</mi><mn>,</mn><mover accent="true"><mi>K</mi><mo stretchy="true">¯</mo></mover><mn>,</mn><msup><mi>K</mi><mo>∗</mo></msup><mn>,</mn><msup><mrow><mover accent="true"><mi>K</mi><mo stretchy="true">¯</mo></mover></mrow><mo>∗</mo></msup><mn>,</mn><mi>ϕ</mi><mo stretchy="false">)</mo></mrow></math> . Here the B strands for ( N ,, N ∗ ) and Y ( , , ), K ( K 0 , K + ) and <math id="M41"><mrow><mover accent="true"><mi>K</mi><mo stretchy="true">¯</mo></mover></mrow></math> ( <math id="M42"><mrow><mover accent="true"><mrow><msup><mi>K</mi><mn>0</mn></msup></mrow><mo stretchy="true">¯</mo></mover></mrow></math> , K - ). The overline of B ( Y ) is its antiparticle. The cross sections of these channels are based on the parametrization of experimental data [ 31 31 ]. The <math id="M43"><mrow><mover accent="true"><mi>N</mi><mo stretchy="true">¯</mo></mover><mi>N</mi></mrow></math> annihilation is described by a statistical model with SU(3) symmetry [ 32 32 ], which includes various combinations of possible emitted mesons with the final state up to six particles [ 7 7 ]. A hard core scattering is assumed in two-particle collisions by the Monte Carlo procedures, in which the scattering of two particles is determined by a geometrical minimum distance criterion <math id="M44"><mrow><mi>d</mi><mo>≤</mo><msqrt><mrow><msub><mi>σ</mi><mrow><mi>t</mi><mi>o</mi><mi>t</mi></mrow></msub><mo>/</mo><mi>π</mi></mrow></msqrt></mrow></math> fm weighted by the Pauli blocking of the final states. Here, the total cross section σ tot in mb is the sum of all possible channels. The probability of reaching a channel in a collision is calculated by its contribution of the channel cross section to the total cross section as P ch = ch / tot . The channel is randomly chosen based on the weight of the probability. III. RESULT AND DISCUSSION 1 The emission mechanism of particles produced in antiproton induced reactions is significant in understanding contributions of different reaction channels associated with antiprotons on nucleons and secondary collisions. Fig. 2 Fig. 2 shows the temporal evolutions of pions, kaons, anti-kaons and hyperons in the reaction <math id="M45"><mrow><mover accent="true"><mi>p</mi><mo stretchy="true">¯</mo></mover></mrow></math> + 12 C at an incident momentum of 1 GeV/c. It is shown that the kaons are emitted immediately after the annihilation of antinucleons and nucleons. The secondary collisions of pions and anti-kaons on nucleons retard the emissions in the reaction dynamics, i.e., πN , <math id="M46"><mrow><mover accent="true"><mi>K</mi><mo stretchy="true">¯</mo></mover><mi>N</mi><mo>→</mo><mi>π</mi><mi>Y</mi></mrow></math> etc, which lead to the production of hyperons. At the considered momentum below its threshold energy, e.g., the reaction <math id="M47"><mrow><mover accent="true"><mi>N</mi><mo stretchy="true">¯</mo></mover><mi>N</mi><mo>→</mo><mover accent="true"><mi>Λ</mi><mo stretchy="true">¯</mo></mover><mi>Λ</mi></mrow></math> (p threshold =1.439 GeV/c), hyperons are mainly contributed by secondary collisions after annihilations, i.e., πN KY and <math id="M48"><mrow><mover accent="true"><mi>K</mi><mo stretchy="true">¯</mo></mover><mi>N</mi><mo>→</mo><mi>π</mi><mi>Y</mi></mrow></math> . Phase-space distributions of particles produced in heavy-ion collisions were used to investigate the in-medium properties of hadrons, in particular for strange particles [ 26 26 - 28 28 26-28 ]. It has been shown that the mean-field potentials influence the spectrum structures, i.e., the kinetic energy spectra of inclusive cross sections, rapidity and transverse momentum distributions of particles etc. The phase-space structure of particle emission in antiproton induced reactions is also investigated in this work. Fig. 3 Fig. 3 shows the transverse momentum distributions of neutral particles in collisions of antiprotons on different targets. The hyperon emission is coupled to the strangeness exchange reactions. The yields increase with the mass number of the target. Shown in Fig. 4 Fig. 4 are the rapidity distributions. The target effect is pronounced in the domain of the antiproton-nucleon frame ( y =0.46) for π 0 , K 0 and <math id="M50"><mrow><msup><mrow><mover accent="true"><mi>K</mi><mo stretchy="true">¯</mo></mover></mrow><mn>0</mn></msup></mrow></math> emissions, where a heavier target increases the production of pions and kaons. It leads to an opposite contribution for antikaon production because of the strangeness exchange <math id="M51"><mrow><mover accent="true"><mi>K</mi><mo stretchy="true">¯</mo></mover><mi>N</mi><mo>→</mo><mi>π</mi><mi>Y</mi></mrow></math> . The peak of hyperon emission is close to the target rapidity and more sensitive to the mass. IV. CONCLUSION 1 Dynamics on particle production in antiproton induced reactions, in particular pions, kaons, antikaons and hyperons, has been investigated within the LQMD model. The production of pions, kaons and antikaons are mainly contributed by the annihilation of antiprotons on nucleons. Hyperons are dominated via the meson-nucleon collisions and strangeness exchange reactions. Kaons are emitted promptly after production in nuclear dynamics. Secondary collisions of particles and nucleons retard the dynamical emission of pions, antikaons and hyperons, which further reduces the yields of antikaons and enhances hyperons with increasing the mass number of the target. References M R Clover , R M DeVries , N J DiGiacomo , et al . Low energy antiproton-nucleus interactions . Phys Rev C , 1982 , 26 : 2138 - 2151 . DOI: 10.1103/PhysRevC.26.2138 http://doi.org/10.1103/PhysRevC.26.2138 Low energy antiproton-nucleus interactions A S Ilynov , V I Nazaruk , S E Chigrinov . Nuclear absorption of stopped antiprotons: Multipion-nucleus interactions . Nucl Phys A , 1982 , 382 : 378 - 400 . DOI: 10.1016/0375-9474(82)90352-9 http://doi.org/10.1016/0375-9474(82)90352-9 Nuclear absorption of stopped antiprotons: Multipion-nucleus interactions J Cugnon , P Deneye and J Vandermeulen . Multipion dynamics following antiproton annihilation on nuclei . Nucl Phys A , 1989 , 500 : 701 - 730 . DOI: 10.1016/0375-9474(89)90236-4 http://doi.org/10.1016/0375-9474(89)90236-4 Multipion dynamics following antiproton annihilation on nuclei J Cugnon , P Deneye , J Vandermeulen . Strangeness production in antiproton annihilation on nuclei . Phys Rev C , 1990 , 41 : 1701 - 1718 . DOI: 10.1103/PhysRevC.41.170 http://doi.org/10.1103/PhysRevC.41.170 Strangeness production in antiproton annihilation on nuclei C M Ko and R Yuan . Lambda production from anti-proton annihilation in nuclei . Phys Lett B , 1987 , 192 : 31 - 34 . DOI: 10.1016/0370-2693(87)91136-1 http://doi.org/10.1016/0370-2693(87)91136-1 Lambda production from anti-proton annihilation in nuclei T Gaitanos , A B Larionov , H Lenske , et al . Formation of double-Λ hypernuclei at PANDA . Nucl Phys A , 2012 , 881 : 240 - 254 . DOI: 10.1016/j.nuclphysa.2011.12.010 http://doi.org/10.1016/j.nuclphysa.2011.12.010 Formation of double-Λ hypernuclei at PANDA A B Larionov , T Gaitanos and U Mosel . Kaon and hyperon production in antiproton-induced reactions on nuclei . Phys Rev C , 2012 , 85 : 024614 . DOI: 10.1103/PhysRevC.85.024614 http://doi.org/10.1103/PhysRevC.85.024614 Kaon and hyperon production in antiproton-induced reactions on nuclei Z Q Feng and H Lenske . Particle production in antiproton-induced nuclear reactions . Phys Rev C , 2014 , 89 : 044617 . DOI: 10.1103/PhysRevC.89.044617 http://doi.org/10.1103/PhysRevC.89.044617 Particle production in antiproton-induced nuclear reactions Z Q Feng and G M Jin . Pion production in heavy-ion collisions in the 1 A GeV region . Chin Phys Lett , 2009 , 26 : 062501 . DOI: 10.1088/0256-307X/26/6/062501 http://doi.org/10.1088/0256-307X/26/6/062501 Pion production in heavy-ion collisions in the 1 A GeV region Z Q Feng and G M Jin . Dynamics of pion production in heavy-ion collisions around 1 A GeV energies . Phys Rev C , 2010 , 82 : 044615 . DOI: 10.1103/PhysRevC.82.044615 http://doi.org/10.1103/PhysRevC.82.044615 Dynamics of pion production in heavy-ion collisions around 1 A GeV energies Z Q Feng and G M Jin . Dynamics of strangeness production in heavy-ion collisions near threshold energies . Phys Rev C , 2010 , 82 : 057901 . DOI: 10.1103/PhysRevC.82.057901 http://doi.org/10.1103/PhysRevC.82.057901 Dynamics of strangeness production in heavy-ion collisions near threshold energies Z Q Feng and G M Jin . Probing high-density behavior of symmetry energy from pion emission in heavy-ion collisions . Phys Lett B , 2010 , 683 : 140 - 144 . DOI: 10.1016/j.physletb.2009.12.006 http://doi.org/10.1016/j.physletb.2009.12.006 Probing high-density behavior of symmetry energy from pion emission in heavy-ion collisions Z Q Feng . Constraining the high-density behavior of the nuclear equation of state from strangeness production in heavy-ion collisions . Phys Rev C , 2011 , 83 : 067604 . DOI: 10.1103/PhysRevC.83.067604 http://doi.org/10.1103/PhysRevC.83.067604 Constraining the high-density behavior of the nuclear equation of state from strangeness production in heavy-ion collisions Z Q Feng . Momentum dependence of the symmetry potential and its influence on nuclear reactions . Phys Rev C , 2011 , 84 : 024610 . DOI: 10.1103/PhysRevC.84.024610 http://doi.org/10.1103/PhysRevC.84.024610 Momentum dependence of the symmetry potential and its influence on nuclear reactions Z Q Feng . Nuclear in-medium effects and collective flows in heavy-ion collisions at intermediate energies . Phys Rev C , 2012 , 85 : 014604 . DOI: 10.1103/PhysRevC.85.014604 http://doi.org/10.1103/PhysRevC.85.014604 Nuclear in-medium effects and collective flows in heavy-ion collisions at intermediate energies Z Q Feng . Effective mass splitting of neutron and proton and isospin emission in heavy-ion collisions . Nucl Phys A , 2012 , 878 : 3 - 13 . DOI: 10.1016/j.nuclphysa.2012.01.014 http://doi.org/10.1016/j.nuclphysa.2012.01.014 Effective mass splitting of neutron and proton and isospin emission in heavy-ion collisions Z Q Feng . Transverse emission of isospin ratios as a probe of high-density symmetry energy in isotopic nuclear reactions . Phys Lett B , 2012 , 707 : 83 - 87 . DOI: 10.1016/j.physletb.2011.12.001 http://doi.org/10.1016/j.physletb.2011.12.001 Transverse emission of isospin ratios as a probe of high-density symmetry energy in isotopic nuclear reactions Z Q Feng . Isospin splitting of nucleon effective mass and isospin effect in heavy-ion collisions . Nucl Sci Tech , 2013 , 24 : 050504 . Isospin splitting of nucleon effective mass and isospin effect in heavy-ion collisions Y Gao , L Zhang , H F Zhang , et al . Effect of the momentum dependence of the nuclear symmetry potential on the π-/π+ ratio in heavy-ion collisions . Phys Rev C , 2011 , 83 : 047602 . DOI: 10.1103/PhysRevC.83.047602 http://doi.org/10.1103/PhysRevC.83.047602 Effect of the momentum dependence of the nuclear symmetry potential on the π-/π+ ratio in heavy-ion collisions Y Gao , L Zhang , W Zuo , et al . Probing the nuclear symmetry energy using single and double π-/π+ ratios from mirror reaction systems . Phys Rev C , 2012 , 86 , 034611 . DOI: 10.1103/PhysRevC.86.034611 http://doi.org/10.1103/PhysRevC.86.034611 Probing the nuclear symmetry energy using single and double π-/π+ ratios from mirror reaction systems L Zhang , Y Gao , Y Du , et al . Effect of the momentum dependence of nuclear symmetry potential on the transverse and elliptic flows . Eur Phys J A , 2012 , 48 : 30 . DOI: 10.1140/epja/i2012-12030-9 http://doi.org/10.1140/epja/i2012-12030-9 Effect of the momentum dependence of nuclear symmetry potential on the transverse and elliptic flows A B Larionov , I A Pshenichnov , I N Mishustin , et al . Antiproton-nucleus collisions simulation within a kinetic approach with relativistic mean fields . Phys Rev C , 2009 , 80 : 021601(R) . DOI: 10.1103/PhysRevC.80.021601 http://doi.org/10.1103/PhysRevC.80.021601 Antiproton-nucleus collisions simulation within a kinetic approach with relativistic mean fields J Côté , M Lacombe , B Loiseau , et al . Nucleon-antinucleon optical potential . Phys Rev Lett , 1982 , 48 : 1319 - 1322 . DOI: 10.1103/PhysRevLett.48.1319 http://doi.org/10.1103/PhysRevLett.48.1319 Nucleon-antinucleon optical potential D J Millener , C B Dover and A Gal . Λ-nucleus single-particle potentials . Phys Rev C , 1988 , 38 : 2700 - 2708 . DOI: 10.1103/PhysRevC.38.2700 http://doi.org/10.1103/PhysRevC.38.2700 Λ-nucleus single-particle potentials B E Gibson and III E V Hungerford . A survey of hypernuclear physics . Phys Rep , 1995 , 257 : 349 - 388 . DOI: 10.1016/0370-1573(94)00114-I http://doi.org/10.1016/0370-1573(94)00114-I A survey of hypernuclear physics Z Q Feng . Dynamics of strangeness and collective flows in heavy-ion collisions near threshold energies . Nucl Phys A , 2013 , 919 : 32 - 45 . DOI: 10.1016/j.nuclphysa.2013.10.005 http://doi.org/10.1016/j.nuclphysa.2013.10.005 Dynamics of strangeness and collective flows in heavy-ion collisions near threshold energies Z Q Feng . Nuclear in-medium and isospin effects on subthreshold kaon production in heavy-ion collisions . Phys Rev C , 2013 , 87 : 064605 . DOI: 10.1103/PhysRevC.87.064605 http://doi.org/10.1103/PhysRevC.87.064605 Nuclear in-medium and isospin effects on subthreshold kaon production in heavy-ion collisions Z Q Feng , W J Xie and G M Jin . Nuclear in-medium effects of strange particles in proton-nucleus collisions . Phys Rev C , 2014 , 90 : 064604 . DOI: 10.1103/PhysRevC.90.064604 http://doi.org/10.1103/PhysRevC.90.064604 Nuclear in-medium effects of strange particles in proton-nucleus collisions D B Kaplan and A E Nelson . Strange goings on in dense nucleonic matter . Phys Lett B , 1986 , 175 : 57 - 63 . DOI: 10.1016/0370-2693(86)90331-X http://doi.org/10.1016/0370-2693(86)90331-X Strange goings on in dense nucleonic matter J Schaffner-Bielich , I N Mishustin and J Bondorf . In-medium kaon production at the mean-field level . Nucl Phys A , 1997 , 625 : 325 - 346 . DOI: 10.1016/S0375-9474(97)81464-9 http://doi.org/10.1016/S0375-9474(97)81464-9 In-medium kaon production at the mean-field level O Buss , T Gaitanos , K Gallmeister , et al . Transport-theoretical description of nuclear reactions . Phys Rep , 2012 , 512 : 1 - 124 . DOI: 10.1016/j.physrep.2011.12.001 http://doi.org/10.1016/j.physrep.2011.12.001 Transport-theoretical description of nuclear reactions E S Golubeva , A S Iljinov , B V Krippa , et al . Effects of mesonic resonance production in annihilation of stopped antiprotons on nuclei . Nucl. Phys. A , 1992 , 537 : 393 – 417 . DOI: 10.1016/0375-9474(92)90362-N http://doi.org/10.1016/0375-9474(92)90362-N Effects of mesonic resonance production in annihilation of stopped antiprotons on nuclei 10.13538/j.1001-8042/nst.26.S20512 Supported by the Major State Basic Research Development Program in China (Nos. 2014CB845405 and 2015CB856903), the National Natural Science Foundation of China Projects (Nos. 11175218 and U1332207) and the Youth Innovation Promotion Association of Chinese Academy of Sciences 2015-01-05 2015-01-28 2015-04-20 2015-04-20 2021-04-11 2015 Nuclear Science and Techniques 2 26