Electromagnetic field effects on nucleon transverse momentum and fragment distribution for heavy ion collisions around 100A MeV

Electromagnetic field effects on nucleon transverse momentum and fragment distribution for heavy ion collisions around 100 A MeV 1 2 first-author 1 3 corresp Corresponding author ygma@sinap.ac.cn Corresponding author ygma@sinap.ac.cn ygma@sinap.ac.cn Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800, China Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800, China University of Chinese Academy of Sciences, Beijing 100049, China University of Chinese Academy of Sciences, Beijing 100049, China ShanghaiTech University, Shanghai 200031, China ShanghaiTech University, Shanghai 200031, China With taking electromagnetic field into account for the transport model of Boltzmann-Uehling-Uhlenbeck, electromagnetic effects are studied for 208 Pb+ 208 Pb collisions around 100 A MeV. Electromagnetic field evolution during the collisions was estimated. It was found that the electric field has an obvious effect on the transverse momentum ( p T ) spectra of nucleons during heavy ion collisions, and leads to different minimum position for the peak of p T spectra of nucleons versus beam energy when the electric field is switched on. For the magnetic field, it affects the z -axis direction distributions of nucleons for central heavy ion collisions at lower energy. Heavy-ion collision Electromagnetic filed Transverse momentum 1 Introduction 1 In heavy-ion collisions, strong electrical fields and magnetic fields should be considered which was firstly pointed out by the Rafelski and Müller in Ref.[ 1 1 ]. With our well knowledge, the Coulomb force causes the repulsive interactions for the charge particles and the Lorentz force can alter the velocities of charge particles. Many efforts have been paid for the Coulomb effects in heavy ion collisions [ 2 2 - 8 8 2-8 ]. The electric fields treated as static electric field (Coulomb fields) has been taken into account in most transport models for the study of nuclear matter of intermediate energy heavy ion collision. In Ref. [ 9 9 ], the electromagnetic fields are considered with the velocities of charge particles in heavy ion collisions at intermediate energies within the Boltzmann-Uehling-Uhlenbeck (BUU) transport model. The magnetic field created in noncentral Au+Au collision at <math id="M1"><mrow><msqrt><mi>s</mi></msqrt><mo>=</mo><mn>200</mn></mrow></math> GeV at the Relativistic Heavy Ion Collider (RHIC) can reach about <math id="M2"><mrow><mi>e</mi><mi>B</mi><mo>∼</mo><msubsup><mi>m</mi><mi>π</mi><mn>2</mn></msubsup><mo>∼</mo><msup><mrow><mn>10</mn></mrow><mrow><mn>18</mn></mrow></msup></mrow></math> Gauss (in Lorentz units, 1 MeV 2 ≈ 5.11 × 10 13 Gauss) [ 9 9 - 13 13 9-13 ]. Experimental and theoretical indications showed that the QCD topological effects such as the Chiral Magnetic Eff ect (CME) in the intense electromagnetic fields may be the evidence of local parity violation in strong interactions [ 15 15 - 21 21 15-21 ]. The properties of strange quark matter in a strong external electric and magnetic field will be also significantly influenced [ 22 22 - 24 24 22-24 ]. Very recently, the strongest votical motion has been deduced from the polarization measurement of Λ hyperon in Au + Au collisions at 200 GeV/ c by the STAR collaboration [ 25 25 ] which is consistent with the prediction of the global polarization of the quark-gluon plasma [ 26 26 ]. But electromagnetic field, especially for the magnetic field, are seldom investigated at low energy for heavy ion collisions where is far from pion production energies. So the aim of this paper is to study the electromagnetic effects at low energies ranging from 40 MeV/nucleon to 120 MeV/nucleon for Pb + Pb heavy ion collisions. The organization of the paper is as follows: In Sec. 2 we give a brief introduction of the simulation model including electromagnetic field. Results of electromagnetic field and electromagnetic field effects are discussed in Sec. 3.1. Finally, a summary is given in the end. 2 Theoretical Framework 1 2.1 Boltzmann-Uehling-Uhlenbeck model 2 The BUU model is a very popular tool for describing intermediate-energy heavy ion collision [ 27 27 , 28 28 ], which is a one-body mean-field theory based upon the Boltzmann equation [ 29 29 ]. The BUU equation reads [ 30 30 ]: <math display="block" id="M3"><mtable columnalign="left"><mtr><mtd><mfrac><mrow><mo>∂</mo><mi>f</mi></mrow><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mo>+</mo><mi>υ</mi><msub><mo>∇</mo><mi>r</mi></msub><mi>f</mi><mo>−</mo><msub><mo>∇</mo><mi>r</mi></msub><mi>U</mi><msub><mo>∇</mo><mi>p</mi></msub><mi>f</mi></mtd></mtr><mtr><mtd><mo>=</mo><mfrac><mn>4</mn><mrow><msup><mrow><mo stretchy="false">(</mo><mn>2</mn><mi>π</mi><mo stretchy="false">)</mo></mrow><mn>3</mn></msup></mrow></mfrac><mstyle displaystyle="true"><mrow><mo>∫</mo></mrow></mstyle><msup><mtext>d</mtext><mn>3</mn></msup><msub><mi>p</mi><mn>2</mn></msub><msup><mtext>d</mtext><mn>3</mn></msup><msub><mi>p</mi><mn>3</mn></msub><mtext>d</mtext><mi>Ω</mi><mfrac><mrow><mtext>d</mtext><msub><mi>σ</mi><mrow><mtext>NN</mtext></mrow></msub></mrow><mrow><mi>d</mi><mi>Ω</mi></mrow></mfrac><msub><mi>υ</mi><mrow><mn>12</mn></mrow></msub><mo>×</mo><mo stretchy="false">[</mo><msub><mi>f</mi><mn>3</mn></msub><msub><mi>f</mi><mn>4</mn></msub><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><mi>f</mi><mo stretchy="false">)</mo></mtd></mtr><mtr><mtd><mo>×</mo><mo stretchy="false">(</mo><mn>1</mn><mtext>​</mtext><mo>−</mo><mtext>​</mtext><msub><mi>f</mi><mn>2</mn></msub><mo stretchy="false">)</mo><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mi>f</mi><msub><mi>f</mi><mn>2</mn></msub><mo stretchy="false">(</mo><mn>1</mn><mtext>​</mtext><mo>−</mo><mtext>​</mtext><msub><mi>f</mi><mn>3</mn></msub><mo stretchy="false">)</mo><mo stretchy="false">(</mo><mn>1</mn><mtext>​</mtext><mo>−</mo><mtext>​</mtext><msub><mi>f</mi><mn>4</mn></msub><mo stretchy="false">)</mo><mo stretchy="false">]</mo><msup><mi>δ</mi><mn>3</mn></msup><mo stretchy="false">(</mo><mi>p</mi><mo>+</mo><msub><mi>p</mi><mn>2</mn></msub><mtext>​</mtext><mo>−</mo><mtext>​</mtext><msub><mi>p</mi><mn>3</mn></msub><mtext>​</mtext><mo>−</mo><mtext>​</mtext><msub><mi>p</mi><mn>4</mn></msub><mo stretchy="false">)</mo><mn>,</mn></mtd></mtr></mtable></math> where f = f (r,p, t ) is the phase-space distribution function. One can solve this equation by the method of Bertsch and Das Gupta [ 31 31 ]. In Eq.(1), <math id="M4"><mrow><mfrac><mrow><mtext>d</mtext><msub><mi>σ</mi><mrow><mtext>NN</mtext></mrow></msub></mrow><mrow><mtext>d</mtext><mi>Ω</mi></mrow></mfrac></mrow></math> and υ 12 are nucleon-nucleon cross section and the relative velocity for the colliding nucleons, respectively. The nucleon-nucleon cross sections include the elastic and inelastic ones, and the well-known parametrization of Cugnon, Mizutani, and Vandermeulun [ 32 32 ] as a function of the available center-of-mass energy has been used. The details of the BUU model were described in Ref. [ 31 31 ]. And the mean-field potential U including the isospin-dependent symmetry energy term can be given: <math display="block" id="M5"><mrow><mi>U</mi><mo stretchy="false">(</mo><mi>ρ</mi><mn>,</mn><msub><mi>τ</mi><mi>z</mi></msub><mo stretchy="false">)</mo><mo>=</mo><mi>a</mi><mo stretchy="false">(</mo><mfrac><mi>ρ</mi><mrow><msub><mi>ρ</mi><mn>0</mn></msub></mrow></mfrac><mo stretchy="false">)</mo><mo>+</mo><mi>b</mi><msup><mrow><mo stretchy="false">(</mo><mfrac><mi>ρ</mi><mrow><msub><mi>ρ</mi><mn>0</mn></msub></mrow></mfrac><mo stretchy="false">)</mo></mrow><mi>κ</mi></msup><mo>+</mo><msub><mi>C</mi><mrow><mtext>sym</mtext></mrow></msub><mo stretchy="false">(</mo><mfrac><mrow><msub><mi>ρ</mi><mtext>n</mtext></msub><mo>−</mo><msub><mi>ρ</mi><mtext>p</mtext></msub></mrow><mrow><msub><mi>ρ</mi><mn>0</mn></msub></mrow></mfrac><mo stretchy="false">)</mo><msub><mi>τ</mi><mi>z</mi></msub><mn>,</mn></mrow></math> where ρ 0 ( ρ 0 = 0.168 fm -3 ), ρ n and ρ p are the normal nuclear matter, neutron and proton densities; and τ z equals 1 or -1 for neutrons and protons, respectively; and C sym is the coefficient of the symmetry energy term (here C sym = 32 MeV is used). The coefficients a , b , and κ are parameters for the nuclear equation of state (EoS). In this work, we use one set of mean-field parameters, the semi-soft EoS with the compressibility K of 235 MeV ( a = -218 MeV, b = 164 MeV, κ = 4/3). One can appropriately implement the electrical and magnetic fields into transport model according to the Liénard-Wiechert potentials at a position <math id="M6"><mover accent="true"><mi>r</mi><mo>→</mo></mover></math> and time t [ 9 9 , 33 33 ]: <math display="block" id="M7"><mrow><mi>e</mi><mover accent="true"><mi>E</mi><mo>→</mo></mover><mo stretchy="false">(</mo><mover accent="true"><mi>r</mi><mo>→</mo></mover><mn>,</mn><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><msup><mi>e</mi><mn>2</mn></msup></mrow><mrow><mn>4</mn><mi>π</mi><msub><mi>ϵ</mi><mn>0</mn></msub></mrow></mfrac><mstyle displaystyle="true"><munder><mo>∑</mo><mi>n</mi></munder></mstyle><msub><mi>Z</mi><mi>n</mi></msub><mfrac><mrow><msup><mi>c</mi><mn>2</mn></msup><mo>−</mo><msubsup><mi>υ</mi><mi>n</mi><mn>2</mn></msubsup></mrow><mrow><msup><mrow><mo stretchy="false">(</mo><mi>c</mi><msub><mi>R</mi><mi>n</mi></msub><mo>−</mo><msub><mover accent="true"><mi>R</mi><mo>→</mo></mover><mi>n</mi></msub><mo>⋅</mo><msub><mover accent="true"><mi>υ</mi><mo>→</mo></mover><mi>n</mi></msub><mo stretchy="false">)</mo></mrow><mn>3</mn></msup></mrow></mfrac><mo stretchy="false">(</mo><mi>c</mi><msub><mover accent="true"><mi>R</mi><mo>→</mo></mover><mi>n</mi></msub><mo>−</mo><msub><mi>R</mi><mi>n</mi></msub><msub><mover accent="true"><mi>υ</mi><mo>→</mo></mover><mi>n</mi></msub><mo stretchy="false">)</mo><mn>,</mn></mrow></math> <math display="block" id="M8"><mrow><mi>e</mi><mover accent="true"><mi>B</mi><mo>→</mo></mover><mo stretchy="false">(</mo><mover accent="true"><mi>r</mi><mo>→</mo></mover><mn>,</mn><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><msup><mi>e</mi><mn>2</mn></msup></mrow><mrow><mn>4</mn><mi>π</mi><msub><mi>ϵ</mi><mn>0</mn></msub><mi>c</mi></mrow></mfrac><mstyle displaystyle="true"><munder><mo>∑</mo><mi>n</mi></munder></mstyle><msub><mi>Z</mi><mi>n</mi></msub><mfrac><mrow><msup><mi>c</mi><mn>2</mn></msup><mo>−</mo><msubsup><mi>υ</mi><mi>n</mi><mn>2</mn></msubsup></mrow><mrow><msup><mrow><mo stretchy="false">(</mo><mi>c</mi><msub><mi>R</mi><mi>n</mi></msub><mo>−</mo><msub><mover accent="true"><mi>R</mi><mo>→</mo></mover><mi>n</mi></msub><mo>⋅</mo><msub><mover accent="true"><mi>υ</mi><mo>→</mo></mover><mi>n</mi></msub><mo stretchy="false">)</mo></mrow><mn>3</mn></msup></mrow></mfrac><msub><mover accent="true"><mi>υ</mi><mo>→</mo></mover><mi>n</mi></msub><mo>×</mo><msub><mover accent="true"><mi>R</mi><mo>→</mo></mover><mi>n</mi></msub><mn>,</mn></mrow></math> where in the left-hand side, an additional charge e is in order to get the electromagnetic fine structure constant α = e 2 /4π = 1/137 (setting ϵ 0 = ħ = c = 1) in the right side of Eqs.(3) and (4). And Z n is the charge number of the n th particle; <math id="M9"><mrow><msub><mover accent="true"><mi>R</mi><mo>→</mo></mover><mi>n</mi></msub><mo>=</mo><mover accent="true"><mi>r</mi><mo>→</mo></mover><mo>−</mo><mover accent="true"><mi>r</mi><mo>→</mo></mover><msub><mo>'</mo><mi>n</mi></msub></mrow></math> , where <math id="M10"><mrow><mover accent="true"><mi>r</mi><mo>→</mo></mover><msub><mo>'</mo><mi>n</mi></msub></mrow></math> is position of charge particle moving with the velocity <math id="M11"><mrow><msub><mover accent="true"><mi>υ</mi><mo>→</mo></mover><mi>n</mi></msub></mrow></math> at retarded time <math id="M12"><mrow><msub><mi>t</mi><mrow><mtext>r</mtext><mi>n</mi></mrow></msub><mo>=</mo><mi>t</mi><mo>−</mo><mo>|</mo><mover accent="true"><mi>r</mi><mo>→</mo></mover><mo>−</mo><mover accent="true"><mi>r</mi><mo>→</mo></mover><msub><mo>'</mo><mi>n</mi></msub><mo stretchy="false">(</mo><msub><mi>t</mi><mrow><mtext>r</mtext><mi>n</mi></mrow></msub><mo stretchy="false">)</mo><mo>|</mo><mo>/</mo><mi>c</mi></mrow></math> . In non-relativistic approximation, when υ < < c , one can obtain: <math display="block" id="M13"><mrow><mi>e</mi><mover accent="true"><mi>E</mi><mo>→</mo></mover><mo stretchy="false">(</mo><mover accent="true"><mi>r</mi><mo>→</mo></mover><mn>,</mn><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><msup><mi>e</mi><mn>2</mn></msup></mrow><mrow><mn>4</mn><mi>π</mi><msub><mi>ϵ</mi><mn>0</mn></msub></mrow></mfrac><mstyle displaystyle="true"><munder><mo>∑</mo><mi>n</mi></munder></mstyle><msub><mi>Z</mi><mi>n</mi></msub><mfrac><mn>1</mn><mrow><msubsup><mi>R</mi><mi>n</mi><mn>3</mn></msubsup></mrow></mfrac><msub><mover accent="true"><mi>R</mi><mo>→</mo></mover><mi>n</mi></msub><mn>,</mn></mrow></math> <math display="block" id="M14"><mrow><mi>e</mi><mover accent="true"><mi>B</mi><mo>→</mo></mover><mo stretchy="false">(</mo><mover accent="true"><mi>r</mi><mo>→</mo></mover><mn>,</mn><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><msup><mi>e</mi><mn>2</mn></msup></mrow><mrow><mn>4</mn><mi>π</mi><msub><mi>ϵ</mi><mn>0</mn></msub><mi>c</mi></mrow></mfrac><mstyle displaystyle="true"><munder><mo>∑</mo><mi>n</mi></munder></mstyle><msub><mi>Z</mi><mi>n</mi></msub><mfrac><mn>1</mn><mrow><msubsup><mi>R</mi><mi>n</mi><mn>3</mn></msubsup></mrow></mfrac><msub><mover accent="true"><mi>υ</mi><mo>→</mo></mover><mi>n</mi></msub><mo>×</mo><msub><mover accent="true"><mi>R</mi><mo>→</mo></mover><mi>n</mi></msub><mn>,</mn></mrow></math> With adding the electrical and magnetic fields, the Hamilton’s equations of propagations of nucleons become: <math display="block" id="M15"><mrow><mfrac><mrow><mtext>d</mtext><msub><mover accent="true"><mi>p</mi><mo>→</mo></mover><mi>i</mi></msub></mrow><mrow><mtext>d</mtext><mi>t</mi></mrow></mfrac><mo>=</mo><mo>−</mo><msub><mo>∇</mo><mi>r</mi></msub><mi>U</mi><mo stretchy="false">(</mo><msub><mover accent="true"><mi>r</mi><mo>→</mo></mover><mi>i</mi></msub><mo stretchy="false">)</mo><mo>+</mo><msub><mover accent="true"><mi>F</mi><mo>→</mo></mover><mtext>C</mtext></msub><mo stretchy="false">(</mo><msub><mover accent="true"><mi>r</mi><mo>→</mo></mover><mi>i</mi></msub><mo stretchy="false">)</mo><mo>+</mo><msub><mover accent="true"><mi>F</mi><mo>→</mo></mover><mtext>L</mtext></msub><mo stretchy="false">(</mo><msub><mover accent="true"><mi>r</mi><mo>→</mo></mover><mi>i</mi></msub><mo stretchy="false">)</mo><mn>,</mn></mrow></math> <math display="block" id="M16"><mrow><mfrac><mrow><mtext>d</mtext><msub><mover accent="true"><mi>r</mi><mo>→</mo></mover><mi>i</mi></msub></mrow><mrow><mtext>d</mtext><mi>t</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><msub><mover accent="true"><mi>p</mi><mo>→</mo></mover><mi>i</mi></msub></mrow><mrow><msqrt><mrow><msup><mi>m</mi><mn>2</mn></msup><mo>+</mo><msubsup><mover accent="true"><mi>p</mi><mo>→</mo></mover><mi>i</mi><mn>2</mn></msubsup></mrow></msqrt></mrow></mfrac><mn>,</mn></mrow></math> where <math id="M17"><mrow><msub><mover accent="true"><mi>F</mi><mo>→</mo></mover><mtext>C</mtext></msub><mo stretchy="false">(</mo><msub><mover accent="true"><mi>r</mi><mo>→</mo></mover><mi>i</mi></msub><mo stretchy="false">)</mo></mrow></math> and <math id="M18"><mrow><msub><mover accent="true"><mi>F</mi><mo>→</mo></mover><mtext>L</mtext></msub><mo stretchy="false">(</mo><msub><mover accent="true"><mi>r</mi><mo>→</mo></mover><mi>i</mi></msub><mo stretchy="false">)</mo></mrow></math> are Coulomb and magnetic force of i -th particle, respectively. We have to mention that for dealing with the initial distance of target and projectile, we did it as the same way which was mentioned in Ref. [ 9 9 ]. 3 Results and discussion 1 3.1 Space-time evolution of the electromagnetic field 2 First we plot the time evolutions of nucleon density (top panels), electric field in the z direction (middle panels) and magnetic field in the y direction (bottom panels) in the x - z plane at impact parameter of b = 0.5× b max ( b max = <math id="M19"><mrow><mn>1.12</mn><mo>*</mo><mo stretchy="false">(</mo><msubsup><mi>A</mi><mtext>P</mtext><mrow><mn>1</mn><mo>/</mo><mn>3</mn></mrow></msubsup><mo>+</mo><msubsup><mi>A</mi><mtext>T</mtext><mrow><mn>1</mn><mo>/</mo><mn>3</mn></mrow></msubsup><mo stretchy="false">)</mo></mrow></math> ) fm for the 40 MeV/nucleon 208 Pb+ 208 Pb collisions in Fig. 1 Fig. 1 . The characteristics of electric fields (middle panels) in the z direction display the geometric configurations of two collision nuclei as shown in top panels of Fig. 1 Fig. 1 . There are two parts in the zone of each nucleus for the electric field eE z at the early stage of reaction, as shown in Fig. 1 (b1) and (b2) Fig. 1 (b1) and (b2) . At later stage, a compound system occurs between nuclei and the shapes of nuclei influence the distribution of electric field which are formed two parts as shown in Fig. 1 (b3) Fig. 1 (b3) . As time increases, they maintain it except for the strength decreasing. The magnetic field eB y forms differently in comparison with electric field eE z . The magnetic field forms a centre and two spectator zones about 20 fm/ c . When the reactions reach maximum compression, the centre zone of magnetic field becomes smaller and strength becomes larger. Compared the strength of electromagnetic field in relativistic heavy ion collisions [ 10 10 ], the strengths (electric field ≈ 5×10 -3 m π and magnetic field ≈ 5×10 -3 m π ) are too small here. Our calculations for the distributions of nucleon density and electromagnetic field are similar to the calculation by Ou and Li [ 9 9 ]. In addition, we give the evolution of magnetic field in the y direction at centre position R (0,0,0) at different impact parameters and incident energies in Fig. 2 Fig. 2 and Fig. 3 Fig. 3 , respectively. The time evolution of average eB y (0,0,0) around zero for the central collision is shown in Fig. 2 Fig. 2 . As the impact parameters increase, the maximum values increase but for the situation of b =1.5× b max fm. It is due to different definition of initial position of target and projectile. In Fig. 3 Fig. 3 , the maximum values of the average eB y (0,0,0) increase as the beam energies increase. One should notice that the strengths of eB y (0,0,0) are negative in both Fig. 2 Fig. 2 and Fig. 3 Fig. 3 . This is caused by the relative position of target and projectile which is shown in the top of Fig. 1 Fig. 1 . It is opposite compared with the results in Ref. [ 9 9 ]. 3.2 Electric effects on the p T spectra 2 One should be acquainted with Coulomb interaction for the EoS in heavy ion collisions. Coulomb effects on flows have been studied around the balance energy in many literatures, eg. Refs.[ 8 8 , 34 34 , 35 35 ]. In our work, the electric field is included in the simulation model with considering the velocities of charge particles in the reaction system. Here electric effects are explored by the p T spectra of nucleons at beam energy of 40 MeV/nucleon and impact parameter of b =0 fm as shown in Fig. 4 Fig. 4 . It is worth to mention here that distributions of p T (= <math id="M20"><mrow><msqrt><mrow><msubsup><mi>p</mi><mi>x</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>p</mi><mi>y</mi><mn>2</mn></msubsup></mrow></msqrt></mrow></math> ) of nucleons will keep stable when the reactions reach freeze-out stage. And we extracted the distributions from the final stage of reactions in Fig. 4 Fig. 4 . One can see that the p T peak shifts to left when including electric field (with E ). It indicates that the beam energy becomes difficult to convert into the transverse energy with the Coulomb interaction is included at 40 MeV/nucleon, which also makes the system hard thermalized at lower beam energy. Two peak values of the solid line (without E ) and dash line (with E ) are extracted by the fitting and shown in Fig. 4 Fig. 4 . Interestingly, a minimum value emerges from the dependence of peak values on beam energy as shown in Fig. 5(a) Fig. 5(a) . From the square-blue line without E, the minimum value emerges around 80 MeV/nucleon. Also as the dot-red line shows, there is another minimum value when the electric field is considered. In the first case, it could be related to the competition between the attractive part and repulsive part of EoS. For the latter, the additional repulsive interaction can be generated by electric field and makes the minimum value shift to lower energy. In contrast to the lower beam energy case (eg. 40 MeV/nucleon), the peak value of p T becomes larger when the Coulomb interaction is on when beam energy is larger (eg. 120 MeV/nucleon), which indicates the beam energy has higher efficiency to convert into transverse energy. Furthermore, differences of peak point values between the cases w/ and w/o Coulomb interaction increase as beam energy as shown in Fig. 5(b) Fig. 5(b) . And the difference is vanishing around 55 MeV/nucleon. Thus, it should be paid more attentions to the physics mechanism at beam energy of ranging 40 MeV/nucleon to 80 MeV/nucleon. To see the how electric field affects on the nucleons, the evolutions of distributions of p T of proton, neutron and all nucleons are shown in Fig. 6 Fig. 6 . Here we take the case without electric field as a reference. At early reaction stage of 40 fm/ c , the protons are directly affected by the electric field as Fig. 6(a1) Fig. 6(a1) shows. And the distribution of p T is shifted to right. The shift becomes smaller as time increases. It, however, is inverse for the neutron. The difference becomes larger at the final stage. Although electric field affects on the protons directly, the neutrons are also affected during collisions. It is indicated that there is an energy transfer between the protons and neutrons. 3.3 Observation for magnetic effects 2 Pion production as an probe has been employed to investigate the magnetic effects in Ref. [ 9 9 ], however, pion production energy is far from the energy range we discussed in this paper. First, magnetic effects are likewise considered with the distribution of p T as shown in Fig. 7 Fig. 7 . But here, the simulations include the electric field. One can see that the magnetic effects are not obvious from distribution of p T . And Fig. 8 Fig. 8 shows the ratios of neutron to proton of all and free nucleons as a function of rapidity at different impact parameters with and without magnetic field. The magnetic effects on the ratios are very tiny at impact parameters of b = 0 fm and b = 0.5× b max fm for 40 MeV/nucleon collisions. They are well in agreement with the results in Ref. [ 9 9 ]. The observables which considered in this paper basing on the momentum space of nucleons are not sensitive to the magnetic field at low energy heavy ion collisions. For the knowledge, however, the directions of charged particles can be changed by the magnetic field and fluctuations can be arisen. Thus, we consider more observables with the z -axis direction distributions of nucleons.In Fig. 9(a) Fig. 9(a) , the z -axis direction distribution of nucleons are almost symmetric in impact parameter of b =0 fm. When we consider the magnetic field, less nucleons in the central region will be. As beam energy increases in Fig. 9(b) Fig. 9(b) , the difference nearly vanishes. It indicates that nucleons with high kinetic energies are not nearly affected by magnetic field generated in heavy ion collisions. 4 Conclusion 1 In summary, coupling with electromagnetic filed in the transport model BUU, the evolutions of electric field in the z direction and magnetic field in the y direction in the x - z plane are given and the impact parameters and incident energies dependence of the magnetic field at the centre of mass are discussed. Further, we consider the electric effects from the p T spectra of nucleons and two minimum values of peaks of p T spectra were found from different situations of with and without electric field. The Coulomb interaction is very important for the study of heavy ion collisions. The magnetic effects are not obvious from the p T spectra of nucleons and the ratio of neutron to proton of all and free nucleons as a function of rapidity. Nevertheless, magnetic filed have some effects on the z -axis direction distribution of nucleons. References J. Rafelski , B. Müller , Magnetic Splitting of Quasimolecular electronic states in strong fields . Phys. Rev. Lett . 36 , 517 ( 1976 ). doi: 10.1103/PhysRevLett.36.517 http://doi.org/10.1103/PhysRevLett.36.517 Magnetic Splitting of Quasimolecular electronic states in strong fields K. G. Libbrecht , S. E. Koonin , Coulomb distortion of pion spectra from heavy-Ion collisions . Phys. Rev. Lett . 43 , 1581 ( 1979 ). doi: 10.1103/PhysRevLett.43.1581 http://doi.org/10.1103/PhysRevLett.43.1581 Coulomb distortion of pion spectra from heavy-Ion collisions M. Gyulassy , S.K. Kauffmann , Coulomb effects in relativistic nuclear collisions . Nucl. Phys. A 362 , 503 ( 1981 ). doi: 10.1016/0375-9474(81)90507-8 http://doi.org/10.1016/0375-9474(81)90507-8 Coulomb effects in relativistic nuclear collisions N. Auerbach , Coulomb effects in nuclear structure . Phys. Rep . 98 , 273 ( 1983 ). doi: 10.1016/0370-1573(83)90008-X http://doi.org/10.1016/0370-1573(83)90008-X Coulomb effects in nuclear structure J. J. Molitoris , J. B. Hoffer , H. Kruse et al ., Microscopic calculations of collective flow probing the short-range nature of the nuclear force . Phys. Rev. Lett . 53 , 899 ( 1984 ). doi: 10.1103/PhysRevLett.53.899 http://doi.org/10.1103/PhysRevLett.53.899 Microscopic calculations of collective flow probing the short-range nature of the nuclear force S. Pratt , Coherence and Coulomb effects on pion interferometry . Phys. Rev. D 33 , 72 ( 1986 ). doi: 10.1103/PhysRevD.33.72 http://doi.org/10.1103/PhysRevD.33.72 Coherence and Coulomb effects on pion interferometry C. A. Bertulani , G. Baur , Electromagnetic processes in relativistic heavy ion collisions . Phys. Rep . 163 , 299 ( 1988 ). doi: 10.1016/0370-1573(88)90142-1 http://doi.org/10.1016/0370-1573(88)90142-1 Electromagnetic processes in relativistic heavy ion collisions Y.M. Zheng , C.M. Ko , B.A. Li et al ., Elliptic flow in heavy-ion collisions near the balance energy . Phys. Rev. Lett . 83 , 2534 ( 1999 ). doi: 10.1103/PhysRevLett.83.2534 http://doi.org/10.1103/PhysRevLett.83.2534 Elliptic flow in heavy-ion collisions near the balance energy L. Ou , B.A. Li , Magnetic effects in heavy-ion collisions at intermediate energies . Phys. Rev. C 84 , 064605 ( 2011 ). doi: 10.1103/PhysRevC.84.064605 http://doi.org/10.1103/PhysRevC.84.064605 Magnetic effects in heavy-ion collisions at intermediate energies V. Skokov , A.Y. Illarionov , V. Toneev , Estimate of the magnetic field strength in heavy-ion collisions . Int. J. Mod. Phys. A 24 , 5925 ( 2009 ). doi: 10.1142/S0217751X09047570 http://doi.org/10.1142/S0217751X09047570 Estimate of the magnetic field strength in heavy-ion collisions M. Asakawa , A. Majumder , B. Müller , Electric charge separation in strong transient magnetic fields . Phys. Rev. C 81 , 064912 ( 2010 ). doi: 10.1103/PhysRevC.81.064912 http://doi.org/10.1103/PhysRevC.81.064912 Electric charge separation in strong transient magnetic fields A. Bzdak , V. Skokov , Event-by-event fluctuations of magnetic and electric fields in heavy ion collisions . Phys. Lett. B 710 , 171 ( 2012 ). doi: 10.1016/j.physletb.2012.02.065 http://doi.org/10.1016/j.physletb.2012.02.065 Event-by-event fluctuations of magnetic and electric fields in heavy ion collisions W. T. Deng , X. G. Huang , Event-by-event generation of electromagnetic fields in heavy-ion collisions . Phys. Rev. C 85 , 044907 ( 2012 ). doi: 10.1103/PhysRevC.85.044907 http://doi.org/10.1103/PhysRevC.85.044907 Event-by-event generation of electromagnetic fields in heavy-ion collisions K. Hattori , X. G. Huang , Novel quantum phenomena induced by strong magnetic fields in heavy-ion collisions . Nucl. Sci. Tech . 28 , 26 ( 2017 ). doi: 10.1007/s41365-016-0178-3 http://doi.org/10.1007/s41365-016-0178-3 Novel quantum phenomena induced by strong magnetic fields in heavy-ion collisions D.E. Kharzeev , L. D. McLerran , H. J. Warringa , The effects of topological charge change in heavy ion collisions: “Event by event P and CP violation" . Nucl. Phys. A 803 , 227 ( 2008 ). doi: 10.1016/j.nuclphysa.2008.02.298 http://doi.org/10.1016/j.nuclphysa.2008.02.298 The effects of topological charge change in heavy ion collisions: “Event by event P and CP violation" B. I. Abelev , M. M. Aggarwal , Z. Ahammed et al . (STAR Collaboration) , Azimuthal charged-particle correlations and possible local strong parity violation . Phys. Rev. Lett . 103 , 251601 ( 2009 ). doi: 10.1103/PhysRevLett.103.251601 http://doi.org/10.1103/PhysRevLett.103.251601 Azimuthal charged-particle correlations and possible local strong parity violation B. I. Abelev , M. M. Aggarwal , Z. Ahammed et al . (STAR Collaboration) , Observation of charge-dependent azimuthal correlations and possible local strong parity violation in heavy-ion collisions . Phys. Rev. C 81 , 054908 ( 2010 ). doi: 10.1103/PhysRevC.81.054908 http://doi.org/10.1103/PhysRevC.81.054908 Observation of charge-dependent azimuthal correlations and possible local strong parity violation in heavy-ion collisions L. Adamczyk , J. K. Adkins , G. Agakishiev et al . (STAR Collaboration) , Beam-energy dependence of charge separation along the Magnetic Field in Au+Au collisions at RHIC . Phys. Rev. Lett . 113 , 052302 ( 2014 ). doi: 10.1103/PhysRevLett.113.052302 http://doi.org/10.1103/PhysRevLett.113.052302 Beam-energy dependence of charge separation along the Magnetic Field in Au+Au collisions at RHIC L. Adamczyk , J. K. Adkins , G. Agakishiev et al . (STAR Collaboration) , Observation of charge asymmetry dependence of pion elliptic flow and the possible chiral magnetic wave in heavy-ion collisions . Phys. Rev. Lett . 114 , 252302 ( 2015 ). doi: 10.1103/PhysRevLett.114.252302 http://doi.org/10.1103/PhysRevLett.114.252302 Observation of charge asymmetry dependence of pion elliptic flow and the possible chiral magnetic wave in heavy-ion collisions G. L. Ma , X. G. Huang , Possible observables for the chiral electric separation effect in Cu+ Au collisions . Phys. Rev. C 91 , 054901 ( 2015 ). doi: 10.1103/PhysRevC.91.054901 http://doi.org/10.1103/PhysRevC.91.054901 Possible observables for the chiral electric separation effect in Cu+ Au collisions Q. Y. Shou , G. L. Ma , Y. G. Ma , Charge separation with fluctuating domains in relativistic heavy-ion collisions . Phys. Rev. C 90 , 047901 ( 2014 ). doi: 10.1103/PhysRevC.90.047901 http://doi.org/10.1103/PhysRevC.90.047901 Charge separation with fluctuating domains in relativistic heavy-ion collisions C. Peng , G.X. Peng , C.J. Xia et al ., Magnetized strange quark matter in the equivparticle model with both confinement and perturbative interactions . Nucl. Sci. Tech . 27 , 98 ( 2016 ). doi: 10.1007/s41365-016-0095-5 http://doi.org/10.1007/s41365-016-0095-5 Magnetized strange quark matter in the equivparticle model with both confinement and perturbative interactions M. Ruggieri , G.X. Peng , Chiral phase transition of quark matter in the background of parallel electric and magnetic fields . Nucl. Sci. Tech . 27 , 130 ( 2016 ). doi: 10.1007/s41365-016-0139-x http://doi.org/10.1007/s41365-016-0139-x Chiral phase transition of quark matter in the background of parallel electric and magnetic fields S.S. Cui , G.X. Peng , Z.Y. Lu , et al ., Properties of color-flavor locked strange quark matter in an external strong magnetic field . Nuclear Science and Techniques. Nucl. Sci. Tech . 26 , 040503 ( 2015 ). doi: 10.13538/j.1001-8042/nst.26.040503 http://doi.org/10.13538/j.1001-8042/nst.26.040503 Properties of color-flavor locked strange quark matter in an external strong magnetic field L. Adamczyk et al . (STAR Collaboration) , Global Λ hyperon polarization in nuclear collisions . Nature 548 , 62 ( 2017 ). doi: 10.1038/nature23004 http://doi.org/10.1038/nature23004 Global Λ hyperon polarization in nuclear collisions Z. T. Liang , X. N. Wang , Globally polarized quark-gluon plasma in noncentral A + A collisions . Phys. Rev. Lett . 94 , 102301 ( 2005 ). doi: 10.1103/PhysRevLett.94.102301 http://doi.org/10.1103/PhysRevLett.94.102301 Globally polarized quark-gluon plasma in noncentral A + A collisions G. F. Bertsch , H. Kruse , S. Das Gupta , Boltzmann equation for heavy ion collisions . Phys. Rev. C 29 , 673 ( 1984 ). doi: 10.1103/PhysRevC.29.673 http://doi.org/10.1103/PhysRevC.29.673 Boltzmann equation for heavy ion collisions H. Kruse , B. V. Jacak , and H. Stöcker , Microscopic theory of pion production and sidewards flow in heavy-ion collisions . Phys. Rev. Lett . 54 , 289 ( 1985 ). doi: 10.1103/PhysRevLett.54.289 http://doi.org/10.1103/PhysRevLett.54.289 Microscopic theory of pion production and sidewards flow in heavy-ion collisions C. Y. Wong , H. H. K. Tang , Extended time-dependent Hartree-Fock approximation with particle collisions . Phys. Rev. Lett . 40 , 1070 ( 1978 ). doi: 10.1103/PhysRevLett.40.1070 http://doi.org/10.1103/PhysRevLett.40.1070 Extended time-dependent Hartree-Fock approximation with particle collisions W. Bauer , G. F. Bertsch , W. Cassing , et al ., Energetic photons from intermediate energy proton-and heavy-ion-induced reactions . Phys. Rev. C 34 , 2127 ( 1986 ). doi: 10.1103/PhysRevC.34.2127 http://doi.org/10.1103/PhysRevC.34.2127 Energetic photons from intermediate energy proton-and heavy-ion-induced reactions G. F. Bertsch , S. Das Gupta , A guide to microscopic models for intermediate energy heavy ion collisions . Phys. Rep . 160 , 189 ( 1988 ). doi: 10.1016/0370-1573(88)90170-6 http://doi.org/10.1016/0370-1573(88)90170-6 A guide to microscopic models for intermediate energy heavy ion collisions J. Cugnon , T. Mizutani , J. Vandermeulen , Equilibration in relativistic nuclear collisions . A Monte Carlo calculation. Nucl. Phys. A 352 , 505 ( 1981 ). doi: 10.1016/0375-9474(81)90427-9 http://doi.org/10.1016/0375-9474(81)90427-9 Equilibration in relativistic nuclear collisions V. Voronyuk , V. D. Toneev , W. Cassing et al ., Electromagnetic field evolution in relativistic heavy-ion collisions . Phys. Rev. C 83 , 054911 ( 2011 ). doi: 10.1103/PhysRevC.83.054911 http://doi.org/10.1103/PhysRevC.83.054911 Electromagnetic field evolution in relativistic heavy-ion collisions A. D. Sooda , R. K. Puri , J. Aichelin , Study of balance energy in central collisions for heavier nuclei . Phys. Lett. B 594 , 260 ( 2004 ). doi: 10.1016/j.physletb.2004.05.053 http://doi.org/10.1016/j.physletb.2004.05.053 Study of balance energy in central collisions for heavier nuclei Y. X. Zhang and Z. X. Li , Elliptic flow and system size dependence of transition energies at intermediate energies . Phys. Rev. C 74 , 014602 ( 2006 ). doi: 10.1103/PhysRevC.74.014602 http://doi.org/10.1103/PhysRevC.74.014602 Elliptic flow and system size dependence of transition energies at intermediate energies 10.1007/s41365-017-0337-1 This work was supported by the National Natural Science Foundation of China (Nos. 11421505, 11305239 and 11220101005), the Major State Basic Research Development Program in China (No. 2014CB845401), the Key Research Program of Frontier Sciences of CAS (No.QYZDJSSW-SLH002), the Strategic Priority Research Program of the Chinese Academy of Sciences (No. XDB16). 2017-10-12 2017-11-02 2017-11-04 2017-11-22 2017-12-01 2022-08-18 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2017 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-017-0337-1 http://dx.doi.org/10.1007/s41365-017-0337-1 . 2017 Nuclear Science and Techniques 12 28