# A Systematic Study of Magnetic Field in Relativistic Heavy-Ion Collisions in the RHIC and LHC Energy Regions

Authors: Yang Zhong; Chun-Bin Yang; Xu Cai; Sheng-Qin Feng
Journal: Advances in High Energy Physics
 (2014)

---

## Abstract

The features of magnetic field in relativistic heavy-ion collisions are systematically studied by using a modified magnetic field model in this paper. The features of magnetic field distributions in the central point are studied in the RHIC and LHC energy regions. We also predict the feature of magnetic fields at LHCs
N
N
=
900, 2760, and 7000 GeV based on the detailed study at RHIC s
N
N
=
62.4, 130, and 200 GeV. The dependencies of the features of magnetic fields on the collision energies, centralities, and collision time are systematically investigated, respectively.


---

## Body

## 1. Introduction

Collisions of two heavy nuclei at high energy serve as a means for creating and exploring strongly interacting matter at high possible energy densities where a new extreme state of matter, the deconfined quark-gluon plasma (QGP), is expected to be formed [1–3]. Besides the formation of QGP, relativistic heavy-ion collisions create also extremely strong (electro)magnetic field due to the relativistic motion of the colliding heavy ions carrying large positive electric charge [3–6].Now we turn to a key question: can the Chiral magnetic effect (CME) occur in heavy-ion collisions? The answer to the question seems to be positive. Two elements are needed for the CME to occur: an external magnetic field and a locally nonzero axial charge density. The relativistically moving heavy ions, typically with large positive charges (e.g., +79e, for Au), carry strong magnetic (and electric) fields with them. In the short moments before/during/after the impact of two ions in noncentral collisions, there is a very strong magnetic field in the reaction zone [5–7]. In fact, such a magnetic field is estimated to be of the order of m
π
2
≈
1
0
18 Gauss [4, 8, 9], probably the strongest, albeit transient, magnetic field in the present universe. The other required element, a locally nonvanishing axial charge density, can also be created in the reaction zone during the collision process through sphaleron transitions (see, e.g., [10] for discussions and references therein). As such, it appears at least during the very early stage of a heavy-ion collision and there can be both strong magnetic field and nonzero axial charge density in the created hot matter.It is suggested by [5, 6, 11–13] that off-central heavy-ion collisions can create strong transient magnetic fields due to the fast, oppositely directed motion of two colliding ions. Thus, heavy-ion collisions provide a unique terrestrial environment to study QCD in strong magnetic fields [14–17]. It has been shown that a strong magnetic field can convert topological charge fluctuations in the QCD vacuum into global electric charge separation with respect to the reaction plane. This so-called Chiral magnetic effect may serve as a sign of the local P and CP violation of QCD.Experimentally, the STAR [18–21], PHENIX [22], and ALICE [23] collaborations have presented the measurements of CME by the two-particle or three-particle correlations of charged particles with respect to the reaction plane, which are qualitatively consistent with the CME. A clear signal compatible with a charge dependent separation relative to the reaction plane is observed, which shows little or no collision energy dependence when compared to measurements at RHIC energies. This provides a new insight into understanding the nature of the charge-dependent azimuthal correlations observed at RHIC and LHC energies.Kharzeev, Mclerran, and Warringa (KMW) [5] suggested that, in the vicinity of the deconfinement phase transition and under the influence of the strong magnetic field generated by the colliding nuclei, the quark spin alignment along the direction of the angular momentum and the imbalance of the left- and right-handed quarks generate an electromagnetic current. The experimental search of those effects has intensified recently, following the realization that the consequent quark fragmentation into charged hadrons results in a charge separation along the direction of the magnetic field and perpendicular to the reaction plane.The magnitude of the effect either should not change or should decrease with increasing energy as long as a deconfined state of matter is formed in a heavy-ion collision. In addition, in [24, 25] it is also suggested that there should be no energy dependence between the top RHIC and the LHC energies, based on arguments related to the universality of the underlying physical process, without, however, explicitly quantifying what the contribution of the different values and time evolution of the magnetic field for different energies will be. On the other hand, it is argued [26] that the CME should strongly decrease at higher energies because the magnetic field decays more rapidly. Such a spread in the theoretical expectations makes it important to measure the charge-dependent azimuthal correlations at the LHC, where the collision energy is an order of magnitude higher compared to the RHIC.Charge separation needs a symmetry axis along which the separation can take place. The only symmetry axis in a heavy-ion collision is along the angular momentum which is perpendicular to the reaction plane. In central collisions there is no symmetry axis, so in that case charge separation should vanish. The strong magnetic field and the QCD vacuum can both completely be produced in the noncentral nucleus-nucleus collision.In [27], we used the Wood-Saxon nucleon distribution instead of a uniform one to improve the magnetic field calculation of the magnetic field for noncentral collision. It was argued that the magnitude of the magnetic field decreases at higher energies. The detailed research of the magnetic field dependencies on collision energies, impact parameter, and collision times is presented in this paper.The paper is organized as follows. The improved calculation of magnetic field and the comparison of our new results with that given by KMW are described in Section2, along with the predicted results of LHC energy region. A summary is given in Section 3.


## 2. The Improved Calculation of Magnetic Field

The situation with the experimental search for the local strong parity violation drastically changed once it was noticed [28–30] that in noncentral nuclear collisions it would lead to the asymmetry in the emission of positively and negatively charged particle perpendicular to the reaction plane. Such a charge separation is a consequence of the difference in the number of quarks with positive and negative helicities positioned in the strong magnetic field of a noncentral nuclear collision, the so-called Chiral magnetic effect.We begin with a charged particle moving along the direction of thez axis. The magnetic field around it can be given by
(1)
B
→
=
1
c
2
v
→
×
E
→If the movement is relativistic, at the timet
=
0, the charge is the origin of the coordinate. The magnitude of the magnetic field B
→ is given by
(2)
B
=
1
4
π
ε
0
c
2
q
v
(
1
-
β
2
)
sin
θ
r
2
(
1
-
β
2
sin
θ
)
3
/
2
.Now we consider a particle with chargeZ and rapidity Y traveling along the z axis. At t
=
0 the particle can be found at position x
→
⊥
′; the magnetic field at the position x
→
=
(
x
→
⊥
,
z
) caused by the particle is given by
(3)
e
B
→
(
x
→
)
=
Z
α
EM
sinh
⁡
Y
×
(
x
→
⊥
′
-
x
→
⊥
)
×
e
→
z
[
(
x
→
⊥
′
-
x
→
⊥
)
2
+
(
t
sinh
⁡
Y
-
z
cosh
⁡
Y
)
2
]
3
/
2
.Now we suppose that two similar nuclei with chargeZ and radius R are traveling in the positive and negative z direction with rapidity Y
0. At t
=
0 they have a noncentral collision with impact parameter b at the origin point. We take the center of the two nuclei at x
=
±
b
/
2 at time t
=
0 so that the direction of b lies along the x axis.As the nuclei are nearly traveling with the speed of light in typical heavy-ion collision experiments, the Lorentz contraction factorγ is so large that we can consider the two included nuclei as pancake shaped. As a result, the nucleon’s number density of each nucleus at x
→
′
=
(
x
→
⊥
′
,
z
) can be given by
(4)
ρ
s
±
(
x
→
⊥
′
)
=
2
4
/
3
π
R
3
R
2
-
(
x
→
⊥
′
±
b
→
2
)
2
.As a result, it seems that the nucleon distribution on average in a nucleus is an approximate result before considering the Lorentz contraction. In [1], KMW model used the uniform nuclear distribution as the nuclear distribution. But for a real situation, the nucleon distribution is not strictly uniform. It seems more reasonable to use the Wood-Saxon distribution in place of the uniform distribution. We use the Wood-Saxon distribution in this paper:
(5)
n
A
(
r
)
=
n
0
1
+
exp
⁡
(
(
r
-
R
)
/
d
)
,
where  n
0
=
0.17
fm
-
3, d
=
0.54 fm, and the radius R
=
1.12
A
1
/
3
fm. Considering the Lorentz contraction, the density in the two-dimensional plane can be given by
(6)
ρ
±
(
x
→
⊥
′
)
=
N
·
∫
-
R
R
d
z
′
×
n
0
1
+
exp
⁡
(
(
(
x
′
∓
b
/
2
)
2
+
y
′
2
+
z
′
2
-
R
)
/
d
)
,
where N is the normalization constant. The number densities should be normalized as
(7)
∫
‍
d
x
→
⊥
′
ρ
±
(
x
→
⊥
′
)
=
1
.We now estimate the strength of the magnetic field at positionx
→
=
(
x
→
⊥
,
z
) caused by the two traveling nuclei. We are only interested in the time t
>
0, that is, just after the collision. Then we can split the contribution of particles to the magnetic field in the following way:
(8)
B
→
=
B
→
s
+
+
B
→
s
-
+
B
→
p
+
+
B
→
p
-
,
where B
→
s
± and B
→
p
± are the contributions of the spectators and the participants moving in the positive or negative z direction, respectively. For spectators, we assume that they do not scatter at all and that they keep travelling with the beam rapidity Y
0. According to (3), we use the density above and find
(9)
e
B
→
s
±
(
τ
,
η
,
x
→
⊥
)
=
±
Z
α
EM
sinh
⁡
⁡
(
Y
0
∓
η
)
∫
‍
d
2
x
→
⊥
′
ρ
±
(
x
→
⊥
′
)
×
[
1
-
θ
∓
(
x
→
⊥
′
)
]
(
x
→
⊥
′
-
x
→
⊥
)
×
e
→
z
[
(
x
→
⊥
′
-
x
→
⊥
)
2
+
τ
2
sinh
⁡
(
Y
0
∓
η
)
2
]
3
/
2
,
where τ
=
(
t
2
-
z
2
)
1
/
2 is the proper time, η
=
(
1
/
2
)
ln
⁡
[
(
t
+
z
)
/
(
t
-
z
)
] is the space-time rapidity, and
(10)
θ
∓
(
x
→
⊥
′
)
=
θ
[
R
2
-
(
x
→
⊥
′
±
b
→
2
)
2
]
.Here, we would like to neglect the contribution from the particles created by the interactions and so we just need to take into account the contribution of the participants that were originally there. The distribution of participants that remain travelling along the beam axis is given by(11)
f
(
Y
)
=
a
2
sinh
⁡
⁡
(
a
Y
0
)
e
a
Y
,
-
Y
0
≤
Y
≤
Y
0
.
Experimental data shows that a
≈
1
/
2, consistent with the baryon junction stopping mechanism. The contribution of the participants to the magnetic field can be also given by
(12)
e
B
→
p
±
(
τ
,
η
,
x
→
⊥
)
=
±
Z
α
EM
∫
‍
d
2
x
→
⊥
′
∫
‍
d
Y
f
(
Y
)
sinh
⁡
⁡
(
Y
∓
η
)
×
ρ
±
(
x
→
⊥
′
)
θ
∓
(
x
→
⊥
′
)
(
x
→
⊥
′
-
x
→
⊥
)
×
e
→
z
[
(
x
→
⊥
′
-
x
→
⊥
)
2
+
τ
2
sinh
⁡
⁡
(
Y
∓
η
)
2
]
3
/
2
.We calculate the magnetite of the magnetic field at the origin (η
=
0
,
x
→
⊥
=
0) in which case it is pointing in the y direction. We took a Au-Au collision with different beam rapidity and different impact parameters.Figure1 shows that the dependencies of magnetic fields of the central point on impact parameters at different proper times τ at different center-of-mass energies for RHIC and LHC, respectively. Figures 1(a), 1(b), and 1(c) are for the results in RHIC energy region. It is shown from Figures 1(a), 1(b), and 1(c) that the magnitudes of the magnetic fields increase with the increasing impact parameter ay proper time τ
≤
1.0 fm. Figures 1(d), 1(e), and 1(f) are for the results in LHC energy region. Figures 1(d) and 1(e) show that the magnitudes of the magnetic field increase with the increasing impact parameter only at proper time τ
≤
0.05 fm but start to decrease with increasing impact parameter when proper time τ
>
0.05 fm. Figure 1(f) shows that the magnitudes of the magnetic field decrease with the increasing impact parameter at proper time τ
≥
0.02 fm at s
N
N
=
7000 GeV. In comparison with that of RHIC energy region, we find that the magnitudes of the magnetic fields of τ
≥
2 fm fall to zero more rapidly at LHC energy region. The variation characteristics of magnetic field with impact parameter at RHIC energy region are different from that of LHC energy region.The dependencies of magnetic fields on impact parameters at different proper timesτ at different energies of central mass system of s
N
N
=
62.4 GeV (a), 130 GeV (b), 200 GeV (c), 900 GeV (d), 2760 GeV (e), and 7000 GeV (f), respectively.
(a)
(b)
(c)
(d)
(e)
(f)The dependencies of the magnetic field on proper time and impact parameters for Au-Au collisions at RHIC energy region withs = 200 GeV, s = 130 GeV, and s = 64 GeV are shown in Figures 2(a), 2(b), and 2(c), respectively. The same tendency of the dependencies of the magnetic field on proper time and impact parameters at RHIC is observed from Figure 2. The maximum position of magnetic field is located at much smaller proper time τ
~
0.02
fm and much larger impact parameter b
≃
12
fm shown in Figure 2. It is observed that the magnitude of magnetic field decreases sharply with increasing proper time τ and decreases with decreasing impact parameter b.The dependencies of the magnetic field on proper time and impact parameters for Au-Au collisions at RHIC collision energies withs
=
200 GeV (a), s
=
130 GeV (b), and s
=
64 GeV (c), respectively.
(a)
(b)
(c)Figures3(a), 3(b), and 3(c) show that the maximum of magnetic field decrease with the increase of collision energies when s from 900 GeV to 7000 GeV at LHC energy region. It is found that the maximum of magnetic field are about 3.0
×
1
0
3
MeV
2, 4.0
×
1
0
2
MeV
2, and 80.0
MeV
2 at s
=
900, s
=
2760, and s
=
7000 GeV, respectively. These are much smaller than that in RHIC energy region. The same variation characteristics of magnetic field with impact parameter and proper time are observed at LHC energy region and that of RHIC.The dependencies of the magnetic field on proper time and impact parameters for Pb-Pb collisions at LHC collision energies withs
=
900 GeV (a), s
=
2760 GeV (b), and s
=
7000 GeV (c), respectively.
(a)
(b)
(c)The dependencies of the magnetic field on collision energiess at RHIC energy region at different impact parameters for Au-Au collisions at different proper time τ are shown in Figure 4. It is found that the dependencies of magnetic field on collision energies at different impact parameter are all on the rising trend at τ
=
0.02 fm and τ
=
0.05 fm as shown in Figures 4(a) and 4(b). Figure 4(c) (τ
=
0.2 fm) shows that the magnetic fields increase with the increase of collision energies when s
≤
30 GeV but then decrease with the increase of collision energies when s
>
30 GeV. As the proper time τ increases to τ
≥
0.5 as shown in Figures 4(d), 4(e), 4(f), and 4(g), the magnetic fields decrease with the increase of collision energies and decrease more rapidly of more off-central collisions of b
=
12 than the more central collisions of b
=
2 and 4 fm.The dependencies of the magnetic field on collision energies of central mass system less than 200 GeV at RHIC energy region at different impact parameters for Au-Au collisions at different proper times forτ
=
0.02 fm (a), τ
=
0.05 fm (b), τ
=
0.2 fm (c), τ
=
0.5 fm (d), τ
=
1.0 fm (e), τ
=
2.0 fm (f), and τ
=
3.0 fm (g), respectively.
(a)
(b)
(c)
(d)
(e)
(f)
(g)Figure5 shows that the dependencies of the magnetic field on collision energies of central mass system of LHC energy region of Pb-Pb collisions at different proper times at different impact parameters for b
=
2 fm, b
=
4 fm, b
=
8 fm, and b
=
12 fm with different proper time τ. It is found that the magnitudes of magnetic fields decrease with increasing energy at LHC energy region, and the speed of decreasing at large impact parameter such as b
=
12 fm (shown in Figure 5(d)) is more quick than that of more central collision such as b
=
2 fm (shown in Figure 5(a)). The magnetic field approaches zero as τ
≥
1.0 fm when b
=
12 fm.The dependencies of the magnetic field on collision energies of central mass system of LHC energy region of Pb-Pb collisions at different proper times at different impact parameters forb
=
2 fm (a), b
=
4 fm (b), b
=
8 fm (c), and b
=
12 fm (d), respectively.
(a)
(b)
(c)
(d)The dependencies of the magnetic field on collision energies of central mass system and impact parameters at different proper times forτ
=
0.02 fm, 0.05 fm, 1.0 fm, 2.0 fm, and 3.0 fm are shown in Figures 6 and 7, respectively. The three-dimensional schemes of different proper time show that the large magnetic fields are produced during small τ, large impact parameter, and (
s
N
N
)
≤
500 GeV energy regions. The magnetic fields decrease with increasing collision energy during (
s
N
N
)
≥
500 GeV. The maximum of magnetic field is e
B
≃
50
Me
V
2 when τ
=
2.0 fm, but that of τ
=
0.02 fm is e
B
≃
1
0
5
Me
V
2.The dependencies of the magnetic field on collision energies of central mass system and impact parameters at different proper times forτ
=
0.02 fm (a) and τ
=
0.05 fm (b).
(a)
(b)The dependencies of the magnetic field on collision energies of central mass system and impact parameters at different proper times forτ
=
1.0 fm (a), τ
=
2.0 fm (b), and τ
=
3.0 fm (c), respectively.
(a)
(b)
(c)


## 3. Summary and Conclusion

It is found that the relativistically moving heavy ions, typically with large positive charges, carry strong magnetic (and electric) fields with them. In the short moments before/during/after the impact of two ions in noncentral collisions, there is a very strong magnetic field in the reaction zone. In fact, such a magnetic field is estimated to be of the order ofm
π
2
≈
1
0
18 Gauss, probably the strongest magnetic field in the present universe. It appears at least during the very early stage of a heavy-ion collision and there can be both strong magnetic field and nonzero axial charge density in the created hot matter.It is suggested that off-central heavy-ion collisions can create strong transient magnetic fields due to the fast, oppositely directed motion of two colliding ions. Thus, a unique terrestrial environment to study QCD in strong magnetic fields is provided in relativistic heavy-ion collisions. A strong magnetic field can convert topological charge fluctuations in the QCD vacuum into global electric charge separation with respect to the reaction plane. This so-called Chiral magnetic effect may serve as a sign of the local P and CP violation of QCD.The features of magnetic fields at LHCs
N
N
=
900, 2760, and 7000 GeV and at RHIC s
N
N
=
62.4, 130, and 200 GeV are systematically discussed. The dependencies of the features of magnetic fields on the collision energies, centralities, and collision time are systematically investigated, respectively.We show that an enormous magnetic field can indeed be created in off-central heavy-ion collisions. The magnitude of the field is quite large, especially just after the collision, and decreases rapidly with time. The drop velocity increases with the collision energy increase. It is shown that the magnitudes of magnetic fields decrease with increasing energy at LHC energy region, and the speed of decreasing at large impact parameter such asb
=
12 fm is more quick than that of more central collision such as b
=
2 fm.The dependencies of the magnetic field on proper time and impact parameters for at RHIC and LHC energy regions, respectively. In comparison with that of RHIC energy region, we find that the magnitudes of the magnetic fields with proper time fall more rapidly at LHC energy region. The variation characteristics of magnetic field with impact parameter at RHIC energy region are different from that of LHC energy region. The maximum position is located in the small proper time (τ
~
0.02
fm), more off-central collisions of (b
≃
12
fm), and s
N
N
~
200 GeV. The maximum of magnetic field in our calculation is about e
B
≃
1
0
5
Me
V
2 when τ
=
0.02, b
≃
12
fm, and s
N
N
~
200 GeV.When surveying the dependencies of magnetic field on the proper timeτ, we find that the magnitude of the magnetic field decreases with the increasing energy when s
N
N
≥
200 GeV. The systematic research of magnetic field is consistent with the suggestion given by [24, 25] that the Chiral magnetic effect either should not change or should decrease with increasing energy as long as a deconfined state of matter is formed in a heavy-ion collision.



---

*Source: 193039-2014-08-26.xml*