Novel quantum phenomena induced by strong magnetic fields in heavy-ion collisions

Novel quantum phenomena induced by strong magnetic fields in heavy-ion collisions 1 2 first-author 1 3 corresp Corresponding author, huangxuguang@fudan.edu.cn Corresponding author, huangxuguang@fudan.edu.cn huangxuguang@fudan.edu.cn Physics Department and Center for Particle Physics and Field Theory, Fudan University, Shanghai 200433, China. Physics Department and Center for Particle Physics and Field Theory, Fudan University, Shanghai 200433, China. RIKEN-BNL Research Center, Brookhaven National Laboratory, Upton, New York 11973-5000, USA RIKEN-BNL Research Center, Brookhaven National Laboratory, Upton, New York 11973-5000, USA Department of Physics, Brookhaven National Laboratory, Upton, New York 11973-5000, USA Department of Physics, Brookhaven National Laboratory, Upton, New York 11973-5000, USA The relativistic heavy-ion collisions create both hot quark-gluon matter and strong magnetic fields, and provide an arena to study the interplay between quantum chromodynamics (QCD) and quantum electrodynamics (QED). In recent years, it has been shown that such an interplay can generate a number of interesting quantum phenomena in hadronic and quark-gluon matter. In this short review, we first discuss some properties of the magnetic fields in heavy-ion collisions and then give an overview of the magnetic-field induced novel quantum effects. In particular, we focus on the magnetic effect on the heavy-flavor mesons, the heavy quark transports, and the phenomena closely related to chiral anomaly. Relativistic heavy-ion collisions Strong magnetic fields Anomalous transports Quarkonium Heavy quark diffusion dynamics 1 Introduction 1 Understanding the phase structure and, more generally, the many-body physics of quantum chromodynamics (QCD) is one of the most outstanding challenges in contemporary physics. At low temperature and low quark chemical potential, the QCD matter is formed by hadrons; at high temperature and/or quark chemical potentials, we expect the QCD matter where quarks and gluons are liberated from the confinement. The "critical temperature" of the transition from the hadronic matter to the quark-gluon matter is expected to be at the order of the QCD confinement scale, Λ QCD ∼200 MeV, at zero quark chemical potential. The unique terrestrial experiments that can achieve such an extremely high temperature are the relativistic heavy-ion collisions which have been carried out in the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory (BNL) and in the Large Hadron Collider (LHC) at CERN. The experimental data produced in RHIC and LHC have shown many features that support the formation of the hot quark-gluon plasma (QGP) and have also revealed a number of interesting properties of the QGP. In noncentral heavy-ion collisions, strong magnetic fields can also be created [ 1 1 - 4 4 1-4 ]. This is because that the two colliding nuclei generate two electric currents in opposite directions and thus produce a magnetic field perpendicular to the reaction plane (the plane defined by the impact parameter and the beam direction). The typical strength of the magnetic fields can be roughly estimated by using the Biot-Savart formula, <math id="M1"><mrow><mi>e</mi><mi>B</mi><mo>∼</mo><mi>γ</mi><msub><mi>α</mi><mrow><mtext>E</mtext><mi>M</mi></mrow></msub><mi>Z</mi><mo>/</mo><msubsup><mi>R</mi><mi>A</mi><mn>2</mn></msubsup></mrow></math> , where γ is the Lorentz gamma factor associated with the moving nuclei, α EM is the fine structure constant, Z is the charge number of the nucleus, and R A is the radius of the nucleus. Thus one finds that the magnetic field in RHIC Au + Au collisions at <math id="M2"><mrow><msqrt><mi>s</mi></msqrt><mo>=</mo><mn>200</mn></mrow></math> GeV can be as large as 10 19 Gauss and in LHC Pb + Pb collisions at <math id="M3"><mrow><msqrt><mi>s</mi></msqrt><mo>=</mo><mn>2.76</mn></mrow></math> TeV can reach the order of 10 20 Gauss. More detailed numerical investigations indeed confirmed the generation of strong magnetic fields in heavy-ion collisions and also revealed many extreme aspects of the magnetic fields [ 5 5 - 16 16 5-16 ]. As the strength of the magnetic field is comparable to the QCD scale, <math id="M4"><mrow><mi>e</mi><mi>B</mi><mo>∼</mo><msubsup><mi>Λ</mi><mrow><mtext>Q</mtext><mi>C</mi><mi>D</mi></mrow><mn>2</mn></msubsup><mo>∼</mo><msup><mrow><mn>10</mn></mrow><mrow><mn>18</mn></mrow></msup></mrow></math> Gauss, we expect that the magnetic-field induced effects may compose important ingredients in exploring the QCD physics in heavy-ion collisions. Indeed, recently, a number of novel quantum phenomena induced by strong magnetic fields in heavy-ion collisions have been found and have attracted a lot of attention in the high-energy and nuclear physics community. In this article, we will give a short review of the recent progress on the study of these novel quantum phenomena. We will focus on the following subjects: chiral magnetic effect [ 5 5 , 17 17 ], chiral separation effect [ 18 18 , 19 19 ], chiral magnetic wave [ 20 20 - 22 22 20-22 ], magnetic-field effects on the heavy-flavor mesons [ 23 23 - 39 39 23-39 ], and the heavy-quark transport [ 40 40 - 42 42 40-42 ]. We will also discuss the experimental implications of these phenomena. In addition to the above mentioned subjects, strong magnetic fields may also drive a number of other phenomena which will not be discussed in this article. These phenomena include the magnetic catalysis of chiral symmetry breaking [ 43 43 - 46 46 43-46 ], the inverse magnetic catalysis or magnetic inhibition at finite temperature and density [ 47 47 - 63 63 47-63 ], the rotating system under a strong magnetic field [ 64 64 - 66 66 64-66 ], the possible ρ meson condensation in strong magnetic field [ 67 67 - 71 71 67-71 ], the neutral pion condensation in vaccum [ 72 72 ], the magnetic-field induced anisotropic viscosities in hydrodynamic equations [ 73 73 - 78 78 73-78 ], and the magnetic-field induced particle production [ 79 79 - 89 89 79-89 ]. Some of these intriguing effects are reviewed in Refs. [ 90 90 - 92 92 90-92 ]. In Sec. 2, we briefly review some features of the magnetic fields in heavy-ion collisions. In Sec. 3, we discuss the chiral magnetic effect, chiral separation effect, chiral magnetic wave, and the current status of the experimental search of these anomalous transport phenomena. In Sec. 4, we discuss how the magnetic fields modify the properties of the heavy-flavor mesons and the heavy-quark transport in a hot medium. Phenomenological implications are also discussed. We summarize in Sec. 5. 2 Magnetic fields in heavy-ion collisions 1 In relativistic heavy-ion collisions, due to the fast motion of the ions which carry positive charges, large magnetic fields can be generated in the reaction zones. In this section, we briefly review some general properties of such-generated magnetic fields by both numerical and analytic models. More detailed information of the numerical study can be found in the original references and in Ref. [ 93 93 ]. 2.1 Numerical approach 2 A number of numerical simulations have revealed that the magnitudes of the magnetic fields in heavy-ion collisions can be very large [ 5 5 - 16 16 5-16 ]. In Fig. 1 Fig. 1 , the different spatial components of the magnetic fields at r = 0 (i.e., the center of the overlapping region) and t =0 (i.e., the time when the two colliding nuclei completely overlap) are presented as functions of the impact parameter, b . The fields are calculated by using the Lienard-Wiechert potential on the event-by-event basis where the events are generated by the HIJING event generater [ 10 10 ]. Please find the discussions on the quantum correction to the otherwise classical Leinard-Wiechert potential in Refs. [ 16 16 , 93 93 ]. The most important information from Fig. 1 Fig. 1 are three folds. First, it shows that heavy-ion collision experiments at RHIC and LHC generate very strong magnetic fields. The event-averaged magnetic fields at Au + Au collisions at <math id="M6"><mrow><msqrt><mi>s</mi></msqrt><mo>=</mo><mn>200</mn></mrow></math> GeV can reach about 10 19 Gauss, and at Pb + Pb collisions at <math id="M7"><mrow><msqrt><mi>s</mi></msqrt><mo>=</mo><mn>2.76</mn></mrow></math> TeV can reach about 10 20 Gauss. They are much stronger than the magnetic fields of neutron stars, including the magnetars which may have surface magnetic fields of the order of 10 14 -10 15 Gauss [ 94 94 , 95 95 ]. They are also much larger than the masses squared of electron, <math id="M8"><mrow><msubsup><mi>m</mi><mi>e</mi><mn>2</mn></msubsup></mrow></math> , and light quarks, <math id="M9"><mrow><msubsup><mi>m</mi><mi>u</mi><mn>2</mn></msubsup><mn>,</mn><msubsup><mi>m</mi><mi>d</mi><mn>2</mn></msubsup></mrow></math> , and thus are capable of inducing significant quantum phenomena. Second, although we do see that only the y -component ( y axis is set to be perpendicular to the reaction plane) of the magnetic fields remains after the event average (due to the left-right symmetry of the noncentral collisions), we can find that the event-by-event fluctuation of the nuclear distribution results in nonzero x -component of the magnetic fields as well. These are reflected in the averaged absolute values of the fields in Fig. 1 Fig. 1 and is most evident for central collisions [ 8 8 , 10 10 ]. Third, the fields in LHC Pb + Pb at 2.76 TeV are roughly 13.8=2.76/0.2 times larger than that in RHIC Au + Au collisions at 200 GeV. In fact, as investigated in Refs. [ 8 8 , 10 10 ], to high precision, the magnitudes of the magnetic fields linearly depend on the collision energy, <math id="M10"><mrow><msqrt><mi>s</mi></msqrt></mrow></math> . This can be understood from the fact that the fields are proportional to the Lorentz gamma factor of the beam motion which is equal to <math id="M11"><mrow><msqrt><mi>s</mi></msqrt><mo>/</mo><mo stretchy="false">(</mo><mn>2</mn><msub><mi>m</mi><mtext>N</mtext></msub><mo stretchy="false">)</mo></mrow></math> , where m N is the mass of a nucleon. The simulations depicted in Fig. 1 Fig. 1 are only for fields at t =0. As the hot medium generated in heavy-ion collisions is expanding fast, a natural question is: how do the magnetic fields evolve in time in the expanding medium? This is still an unsolved problem. We summarize the recent progress briefly here and in the next subsection. As we explained, the strong initial magnetic fields are consequences of the smallness of the system and the relativistically fast motion of the colliding nuclei. The main contributors to the initial magnetic fields are the spectator protons (which do not participate into the collision) in the colliding uclei. They leave the collision zone very fast in a time scale characterized by the Lorentz-suppressed radius of the nucleus, i.e., <math id="M12"><mrow><msub><mi>τ</mi><mi>B</mi></msub><mo>=</mo><msub><mi>R</mi><mi>A</mi></msub><mo>/</mo><mi>γ</mi><mo>≈</mo><mn>2</mn><msub><mi>m</mi><mi>N</mi></msub><msub><mi>R</mi><mi>A</mi></msub><mo>/</mo><msqrt><mi>s</mi></msqrt></mrow></math> . Thus in relativistic heavy-ion collisions, one can expect that the magnetic fields generated by the spectators will decay very fast after a short time scale, τ B . This fact is indeed verified by numerical simulations. In Fig. 2 Fig. 2 , one can find the numerical results of the time evolutions of the magnetic fields at r =0 in collisions with b =10 fm for Au + Au collision at <math id="M13"><mrow><msqrt><mi>s</mi></msqrt><mo>=</mo><mn>200</mn></mrow></math> GeV and for Pb + Pb collision at <math id="M14"><mrow><msqrt><mi>s</mi></msqrt><mo>=</mo><mn>2.76</mn></mrow></math> TeV [ 10 10 ]. We can see that the magnetic fields decay very fast after the collision reflecting the fact that the spectators are leaving the collision region very fast. Once the spectators are all far away from the collision region, the remnant charge carriers, which participate in the collisions and thus moves much slower than the spectators along the beam direction, become the dominant sources for the magnetic fields, and they essentially slow down the decay of the magnetic fields. 2.2 Analytical approach 2 To understand the basic picture of the creation of the magnetic fields, we compared the numerical results shown in the previous section with a simple analytical model proposed in Ref. [ 28 28 ]. With the relativistic beam energy, the momenta of the colliding nuclei are so large that the nuclear stopping power is not really strong enough to stop them: the incident nuclei pass each other and recede from the collision point with almost the same rapidity as the incident beam. Therefore, in this simple model, we assume that a magnetic field is created by protons freely streaming with the beam rapidity, and thus do not distinguish the spectators and participants. From the comparison between the numerical and analytic modeling, we will also discuss effects of the nuclear stopping power. As we have specified with the classical motion of charged particles, the magnitude of the magnetic field can be estimated by using the Lienard-Wiechert potential <math display="block" id="M17"><mrow><mi>e</mi><msup><mtext>B</mtext><mo>±</mo></msup><mo>=</mo><mfrac><mrow><msub><mi>α</mi><mrow><msub><mrow/><mrow><mtext>E</mtext><mi>M</mi></mrow></msub></mrow></msub><msubsup><mi>Z</mi><mrow><mtext>e</mtext><mi>f</mi><mi>f</mi></mrow><mo>±</mo></msubsup></mrow><mrow><msup><mrow><mo stretchy="false">(</mo><msup><mi>r</mi><mo>±</mo></msup><mo stretchy="false">)</mo></mrow><mn>3</mn></msup></mrow></mfrac><mi>sinh</mi><mo stretchy="false">(</mo><mo>±</mo><msub><mi>Y</mi><mrow><mtext>b</mtext><mi>e</mi><mi>a</mi><mi>m</mi></mrow></msub><mo stretchy="false">)</mo><mo stretchy="false">(</mo><msup><mover accent="true"><mtext>x</mtext><mo>˜</mo></mover><mo>±</mo></msup><mo>×</mo><msub><mstyle mathvariant="bold" mathsize="normal"><mi>e</mi></mstyle><mi>z</mi></msub><mo stretchy="false">)</mo><mtext> </mtext><mn>,</mn></mrow></math> where Y beam is the absolute value of the beam rapidity. The positive and negative signs correspond to the nuclei moving in the positive and negative z -direction, respectively. The coordinate system is illustrated in Fig. 3 Fig. 3 (Left). We have defined the coordinate vector <math id="M18"><mover accent="true"><mtext>x</mtext><mo>˜</mo></mover></math> from the the center of the nucleus (± b /2,0, v ± t ) to the observer at x , and its length in the rest frame <math id="M19"><mrow><msup><mi>r</mi><mo>±</mo></msup><mo>=</mo><msqrt><mrow><msup><mrow><mo stretchy="false">(</mo><mi>x</mi><mo>±</mo><mi>b</mi><mo>/</mo><mn>2</mn><mo stretchy="false">)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mi>y</mi><mn>2</mn></msup><mo>+</mo><msup><mi>γ</mi><mn>2</mn></msup><msup><mrow><mo stretchy="false">(</mo><mi>z</mi><mo>−</mo><msup><mi>v</mi><mo>±</mo></msup><mi>t</mi><mo stretchy="false">)</mo></mrow><mn>2</mn></msup></mrow></msqrt></mrow></math> with the impact parameter b , the velocity of the nucleus v ± =± tanh Y beam , and the gamma factor <math id="M20"><mrow><mi>γ</mi><mo>=</mo><mn>1</mn><mo>/</mo><msqrt><mrow><mn>1</mn><mo>−</mo><msup><mrow><mo stretchy="false">(</mo><msup><mi>v</mi><mo>±</mo></msup><mo stretchy="false">)</mo></mrow><mn>2</mn></msup></mrow></msqrt></mrow></math> . The Lienard-Wiechert potential in Eq. (1) is obtained by boosting the electro-static potential in the rest frame of the nuclei (see Fig. 3 Fig. 3 ). As is well known, the net electro-static potential can be calculated from the effective point-like charge sitting at the center of the sphere. Therefore, the effective charge Z eff appearing in Eq. (1) should be the total amount of the charges inside the sphere, i.e., <math display="block" id="M21"><mrow><msubsup><mi>Z</mi><mrow><mtext>e</mtext><mi>f</mi><mi>f</mi></mrow><mo>±</mo></msubsup><mo stretchy="false">(</mo><mi>t</mi><mn>,</mn><mtext>x</mtext><mo stretchy="false">)</mo><mo>=</mo><mn>4</mn><mi>π</mi><mstyle displaystyle="true"><mrow><msubsup><mo>∫</mo><mn>0</mn><mrow><msup><mi>r</mi><mo>±</mo></msup></mrow></msubsup></mrow></mstyle><mtext>​</mtext><mtext>​</mtext><mi>d</mi><msup><mi>r</mi><mo>'</mo></msup><mtext> </mtext><msup><mi>r</mi><mrow><mo>'</mo><mtext> </mtext><mn>2</mn></mrow></msup><msub><mi>ρ</mi><mrow/></msub><mo stretchy="false">(</mo><msup><mi>r</mi><mo>'</mo></msup><mo stretchy="false">)</mo><mtext> </mtext><mn>,</mn></mrow></math> where ρ ( r ′) is the charge distribution inside the nuclei. Note that the upper boundary of the integral is given by the distance r ± from the center to the observer. When the observer is in the exterior of the nucleus r ± ≥ R A , one simply has the total charge of the nucleus Z eff ( r ± > R A )= Z . The impact parameter and time dependences of the effective charge is shown in the right panel of Fig. 3 Fig. 3 . We find significant dependence on both the impact parameter and the time collapsed after the collision, which will be important to reproduce the results from the numerical simulation shown in the previous section. Plugging the effective charge (2) to the Lienard-Wiechert potential (1), we obtain the magnitude of the magnetic field shown in Figs. 4 Figs. 4 , 5 5 and 6 6 . Here, we have taken the Woods-Saxon form for the charge distribution ρ ( r )= N Z /{1+exp[( r - R A )/ a ]} with <math id="M22"><mrow><msub><mi>N</mi><mi>Z</mi></msub><mo>=</mo><mi>Z</mi><msup><mrow><mo>{</mo><mn>4</mn><mi>π</mi><mstyle displaystyle="true"><mrow><msubsup><mo>∫</mo><mn>0</mn><mrow><msub><mi>R</mi><mi>A</mi></msub></mrow></msubsup></mrow></mstyle><mtext>​</mtext><mi>d</mi><msup><mi>r</mi><mo>'</mo></msup><msup><mi>r</mi><mrow><mo>'</mo><mtext> </mtext><mn>2</mn></mrow></msup><msub><mi>ρ</mi><mrow/></msub><mo stretchy="false">(</mo><msup><mi>r</mi><mo>'</mo></msup><mo stretchy="false">)</mo><mo>}</mo></mrow><mrow><mo>−</mo><mn>1</mn></mrow></msup></mrow></math> where Z = 79, R A = 6.38 fm, and a = 0.535 fm for Au, and Z = 82, R A = 6.62 fm, and a = 0.546 fm for Pb. We do not consider the event-by-event fluctuation of the charge distribution in our analytical modeling so that we in Figs. 4 Figs. 4 , 5 5 and 6 6 compare the analytical modeling with the numerical simulations of event-averaged magnetic fields. Figure 4 Figure 4 shows the time dependence of the magnetic field obtained from the analytic model (solid lines) and the event-averaged numerical simulations (dots). The latter has been already shown in Fig. 2. Fig. 2. We find a good agreement between the results from these two models in the time range before and after the collision except for the late time. The origin of the deviation in the late time is, as can be guessed from the asymmetric time dependence in the numerical result with respect to the origin of time, attributed to the effect of the nuclear stopping power. Because of the deceleration by the nuclear stopping power, participant nucleons (also called remnants) will stay in the collision region longer than in the free streaming case so that they will maintain the magnetic field for a longer time. Figure 5 Figure 5 shows the impact parameter dependence of the magnetic field. Again, the sold lines and dots show the results from the analytic model and the event-averaged numerical simulations, respectively, and we find a good agreement between them. The slight deviation may be originated from the energy loss due to the nuclear stopping power in the numerical simulation. In addition to these results, we show the magnetic field created by point-like nuclei carrying the total charge Z , i.e., Z eff ( t , x )= Z in Eq. (1). We find that this simplest model completely fails to reproduce the numerical results when there is an overlap between the nuclei ( <math id="M24"><mrow><mi>b</mi><mo>≲</mo><mn>2</mn><msub><mi>R</mi><mi>A</mi></msub></mrow></math> ), indicating the importance of taking into account the effective charge (2). Finally, Fig. 6 Fig. 6 shows the beam energy dependence of the magnetic field from the analytic model (left panel) and the numerical simulation (right panel) [ 7 7 ]. While we find a good agreement, the nuclear stopping power will play an important role in the low energy collisions such as those in the beam energy scan program in RHIC, FAIR, and J-PARC. In the above discussions, we omit one important issue which may substantially modify the time evolution of the magnetic fields, that is, the finite electric conductivity, σ , of the produced quark-gluon matter. If σ is large, the Faraday induction effect may be strong as well and the lifetime of the magnetic fields can be significantly prolonged. As a matter of fact, the equilibrium QGP is a very good conductor according to the theoretical and lattice QCD studies [ 96 96 - 105 105 96-105 ]; see also Ref. [ 106 106 ] for review. Now let us understand how the large σ induces the strong Faraday induction and prolongs the lifetime of the magnetic fields; more discussions can be found in review [ 93 93 ]. Let us write down the Maxwell’s equations in the background of a charged flow characterized by velocity field v , <math display="block" id="M25"><mrow><mo>∇</mo><mo>×</mo><mi>E</mi><mo>=</mo><mo>−</mo><mfrac><mrow><mo>∂</mo><mi>B</mi></mrow><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mn>,</mn></mrow></math> <math display="block" id="M26"><mrow><mo>∇</mo><mo>×</mo><mi>B</mi><mo>=</mo><mfrac><mrow><mo>∂</mo><mi>E</mi></mrow><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mo>+</mo><mi>J</mi><mn>,</mn></mrow></math> where J is the electric current which according to the Ohm’s law is determined by <math display="block" id="M27"><mrow><mi>J</mi><mo>=</mo><mi>σ</mi><mrow><mo>(</mo><mrow><mi>E</mi><mo>+</mo><mi>v</mi><mo>×</mo><mi>B</mi></mrow><mo>)</mo></mrow><mn>.</mn></mrow></math> Substituting J into Eq. (4) and eliminating E by using Eq. (3), one can obtain an equation that governs the time evolution of the magnetic field <math display="block" id="M28"><mrow><mfrac><mrow><mo>∂</mo><mi>B</mi></mrow><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mo>=</mo><mo>∇</mo><mo>×</mo><mo stretchy="false">(</mo><mi>v</mi><mo>×</mo><mi>B</mi><mo stretchy="false">)</mo><mo>+</mo><mfrac><mn>1</mn><mi>σ</mi></mfrac><mrow><mo>(</mo><mrow><msup><mo>∇</mo><mn>2</mn></msup><mi>B</mi><mo>−</mo><mfrac><mrow><msup><mo>∂</mo><mn>2</mn></msup><mi>B</mi></mrow><mrow><mo>∂</mo><msup><mi>t</mi><mn>2</mn></msup></mrow></mfrac></mrow><mo>)</mo></mrow><mn>,</mn></mrow></math> where we have regarded the QGP as locally neutral (this is a good assumption for high-energy collisions) so that ∇⋅ E = ρ =0. The Faraday prolongation of the lifetime of the magnetic fields is most easily seen in two limits. The first limit is the no-flow limit, i.e., the limit when v can be neglected. In this case, if the electric conductivity is so large that <math id="M29"><mrow><mi>σ</mi><mo>≫</mo><mn>1</mn><mo>/</mo><msub><mi>t</mi><mtext>c</mtext></msub></mrow></math> with t c being the characteristic time scale over which the field strongly varies, one can neglect the second-order time derivative and render Eq. (6) a diffusion equation [ 79 79 ]: <math display="block" id="M30"><mrow><mfrac><mrow><mo>∂</mo><mi>B</mi></mrow><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mo>=</mo><mfrac><mn>1</mn><mi>σ</mi></mfrac><msup><mo>∇</mo><mn>2</mn></msup><mi>B</mi><mn>.</mn></mrow></math> Clearly, in this case, the decay of the magnetic field is due to diffusion, and the diffusion time is given by <math display="block" id="M31"><mrow><msub><mi>t</mi><mtext>D</mtext></msub><mo>=</mo><msup><mi>L</mi><mn>2</mn></msup><mi>σ</mi><mn>,</mn></mrow></math> where L characterizes the length scale of the QGP system. For sufficiently large σ , the diffusion time, t D , could be larger than τ B [ 79 79 ]. The second limit is the ideal magnetohydrodynamic limit, i.e., when σ is very large and at the same time the flow velocity, v , is also large so that we can neglect the term divided by σ in Eq. (6) and obtain <math display="block" id="M32"><mrow><mfrac><mrow><mo>∂</mo><mi>B</mi></mrow><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mo>=</mo><mo>∇</mo><mo>×</mo><mo stretchy="false">(</mo><mi>v</mi><mo>×</mo><mi>B</mi><mo stretchy="false">)</mo><mn>.</mn></mrow></math> This equation is known to lead to the so-called frozen-in theorem , i.e., the magnetic lines are frozen in the ideally conducting plasma. Thus in this limit, the decay of the mangnetic fields is totally due to the expansion of the QGP. For high-energy heavy-ion collisions, the longitudinal expansion is well described by the Bjorken flow, v z = z / t , which makes the area of cross section of the overlapping region on the reaction plane grow linearly in time. Thus the frozen-in theorem leads us to the conclusion that in the ideal magnetohydrodynamic limit, before the transverse expansion play significant role, the magnetic fields will decay inversely proportional to time, B y ( t )∝( t 0 / t ) B y ( t 0 ), which is much slower than the 1/ t 3 -type decay in the insulating case. The transverse expansion of the system can happen because of the transverse gradient of the fluid pressure. When this transverse expansion becomes fast, the decay of the fields will become fast as well, see discussions in Refs. [ 10 10 , 93 93 ]. So far, we discussed two special limiting cases from which we can analytically extract the effect of Faraday induction. However, the realistic environment produced by heavy-ion collisions is not the case as described by the two limits. Moreover, the system is essentially non-equilibrium. Thus, the precise time evolution of the magnetic fields is really a hard task to achieve, and some recent progresses can be found in Refs. [ 10 10 , 15 15 , 79 79 , 82 82 , 107 107 - 113 113 107-113 ]. In the above discussions, we focus only on collisions of the same spherical nuclei, i.e., the Au + Au collisions and Pb + Pb collisions. The Heavy-ion program at RHIC also performed collisions of different nuclei, like Cu + Au collisions, as well as collisions of prolate nuclei like the U + U collisions. Of course, these collisions can also generate strong magnetic fields; see for example, the studies in Refs [ 12 12 , 114 114 - 118 118 114-118 ] and recent review [ 93 93 ]. It is worth mentioning that, owing to the charge number asymmetry between Cu and Au nuclei, the Cu + Au collisions can generate a strong electric field pointing from Au nucleus to Cu nucleus [ 114 114 - 117 117 114-117 , 119 119 ]. Such E field may lead to a splitting between the directed flow v 1 of positively and negatively charged hadrons, which has been recently test by STAR Collaboration at RHIC [ 120 120 ]. 3 Magnetic-field induced anomalous transports in heavy-ion collisions 1 The strong magnetic fields can induce a variety of novel quantum effects in the quark-gluon plasma which we will discuss in the following sections. In this section, we will focus on the ones that are deeply related to the chiral anomaly in QCD and QED, namely, the chiral magnetic effect (CME) and its relatives, the chiral separation effect, the chiral magnetic wave, etc. All these effects represent currents or collective modes that transport vector charges or axial charges in QGP which appear only at quantum level so that they do not have classical counterparts. The recent reviews of these anomalous transports can be found in Refs. [ 93 93 , 121 121 - 126 126 121-126 ]. 3.1 The chiral magnetic effect 2 The CME represents the generation of vector current induced by an external magnetic field in chirality-imbalanced (P-odd) medium [ 5 5 , 17 17 ]. The CME composes a very general class of anomalous transport phenomena spreading in a wide contexts in physics, ranging from astrophysics, condensed matter physics, nuclear physics, to particle physics. We apologize for not being able to mention all the relevant works. We refer the readers to other review articles as mentioned above for more information. The CME can be expressed as <math display="block" id="M33"><mrow><msub><mi>J</mi><mi>V</mi></msub><mo>=</mo><msub><mi>σ</mi><mrow><mtext>CME</mtext></mrow></msub><mi>B</mi><mn>,</mn></mrow></math> <math display="block" id="M34"><mrow><msub><mi>σ</mi><mrow><mtext>CME</mtext></mrow></msub><mo>=</mo><mfrac><mrow><msup><mi>e</mi><mn>2</mn></msup></mrow><mrow><mn>2</mn><msup><mi>π</mi><mn>2</mn></msup></mrow></mfrac><msub><mi>μ</mi><mi>A</mi></msub><mn>,</mn></mrow></math> where e is the charge of the fermion, J V is the vector current defined by <math id="M35"><mrow><msubsup><mi>J</mi><mi>V</mi><mi>μ</mi></msubsup><mo>=</mo><mi>e</mi><mo>〈</mo><mover accent="true"><mi>ψ</mi><mo>¯</mo></mover><msup><mi>γ</mi><mi>μ</mi></msup><mi>ψ</mi><mo>〉</mo></mrow></math> , and μ A is the chiral chemical potential which parameterizes how large the chirality imbalance of the medium is. Note that, the above expression for CME current is for each species of massless fermions. For QGP at high temperature where the masses of up and down quarks can be neglected, the total CME conductivity is given by <math id="M36"><mrow><msub><mi>σ</mi><mrow><mtext>C</mtext><mi>M</mi><mi>E</mi></mrow></msub><mo>=</mo><msub><mi>N</mi><mi>C</mi></msub><msub><mi>μ</mi><mi>A</mi></msub><mstyle displaystyle="true"><msub><mo>∑</mo><mi>f</mi></msub></mstyle><msubsup><mi>q</mi><mi>f</mi><mn>2</mn></msubsup><mo>/</mo><mo stretchy="false">(</mo><mn>2</mn><msup><mi>π</mi><mn>2</mn></msup><mo stretchy="false">)</mo></mrow></math> with q f the charge of quark of flavor, f , and N c =3 the number of color. Although the CME stems from the profound chiral anomaly in gauge theory, its physical origin can be understood in an intuitive way by adopting the Landau quantization picture. We know that the electrically charged fermions in a magnetic field occupy a tower of Landau levels (i.e., the energy spectrum of charged fermions in a magnetic field is discrete) in such a way that the lowest Landau level contains fermions with only one spin polarization (i.e., the spin of positively (negatively) charged fermions is polarized to be along (opposite to) the direction of the magnetic field) while higher Landau levels are occupied by equal numbers of up-polarized and down-polarized fermions. Now, let us consider the situation that the number of right-handed (RH) and the number of left-handed (LH) fermions are unequal, say, for example, N R > N L . Thus there are more RH fermions than LH fermions confined in the lowest Landau level which means that a current consisting of positively charged fermions will flow along the magnetic field as RH fermions lock their momenta with their spins, while a current consisting of negatively charged fermions will flow opposite to the magnetic field. The net effect is the appearance of an electric current along the magnetic field. This is the CME current. Note that the higher Landau levels do not contribute to the CME current as higher Landau levels are occupied symmetrically by spin-up and spin-down fermions and their contribution to the CME current cancel. Now let us quantify the above intuitive argument to see how we can arrive at Eq. (10). To simplify our lives, we consider a system of massless Dirac fermions with only positive charge, e , in an constant magnetic field, B , pointing along z direction which is supposed to be so strong that all the fermions are confined in the lowest Landau level. The number densities of RH and LH fermions are given by <math display="block" id="M37"><mrow><msub><mi>n</mi><mrow><mi>R</mi><mo>/</mo><mi>L</mi></mrow></msub><mo>≡</mo><mfrac><mrow><msup><mi>d</mi><mn>3</mn></msup><msub><mi>N</mi><mrow><mi>R</mi><mo>/</mo><mi>L</mi></mrow></msub></mrow><mrow><mi>d</mi><mi>x</mi><mi>d</mi><mi>y</mi><mi>d</mi><mi>z</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><mi>e</mi><mi>B</mi></mrow><mrow><mn>2</mn><mi>π</mi></mrow></mfrac><mfrac><mrow><msubsup><mi>p</mi><mi>F</mi><mrow><mi>R</mi><mo>/</mo><mi>L</mi></mrow></msubsup></mrow><mrow><mn>2</mn><mi>π</mi></mrow></mfrac><mn>,</mn></mrow></math> where the prefactor eB /(2 π ) is the density of states of the lowest Landau level in the transverse directions and p F denotes the Fermi momentum (which is the momentum that separates the empty states and occupied states). Moreover, chiral fermions on the lowest Landau level move at the speed of light in a way correlated to their spins which tell us that the RH fermions generate a current along z direction, J R = en R , while LH fermions generate a current along - z direction, J L =- en L . Therefore the net vector current reads <math display="block" id="M38"><mrow><msub><mi>J</mi><mi>V</mi></msub><mo>=</mo><msub><mi>J</mi><mi>R</mi></msub><mo>+</mo><msub><mi>J</mi><mi>L</mi></msub><mo>=</mo><mfrac><mrow><msup><mi>e</mi><mn>2</mn></msup><mi>B</mi></mrow><mrow><mn>4</mn><msup><mi>π</mi><mn>2</mn></msup></mrow></mfrac><mo stretchy="false">(</mo><msubsup><mi>p</mi><mi>F</mi><mi>R</mi></msubsup><mo>−</mo><msubsup><mi>p</mi><mi>F</mi><mi>L</mi></msubsup><mo stretchy="false">)</mo><mn>,</mn></mrow></math> which, once we equate <math id="M39"><mrow><msubsup><mi>p</mi><mi>F</mi><mrow><mi>R</mi><mo>/</mo><mi>L</mi></mrow></msubsup></mrow></math> with μ R/L and μ A with ( μ R - μ L )/2, leads us to Eq. (10) and Eq. (11). Here we use the fact that for massless fermions, the chemical potential is equal to the Fermi momentum. We emphasize that although the above argument is for a non-interacting system and relies on Landau quantization picture, the CME conductivity σ CME = e 2 μ A /(2 π 2 ) is actually fixed by the chiral anomaly equation and universal in the sense that it does not involve concrete microscopic scattering among fermions. This is in parallel to the fact that the chiral anomaly equation itself is universal so that it does not receive perturbative correction from scattering between fermions. Apparently, the occurrence of CME requires an environmental parity violation as being characterized by μ A . A nonzero μ A may be induced by different mechanisms in different physical systems. For example, the electro-weak plasma can have a nonzero μ A because the weak interaction does not respect parity symmetry [ 127 127 ]. Recently, it is realized in condensed matter experiments the so-called Weyl or Dirac semimetals, whose band structures permit level crossing around the low energy excitations, are Weyl or Dirac fermions. The CME was observed in these semimetals by applying parallel electric and magnetic fields, which induce a difference between the chemical potentials of RH and LH fermions via the chiral anomaly [ 128 128 - 131 131 128-131 ]; see also discussions in Refs. [ 92 92 , 93 93 , 126 126 , 132 132 ]. In the hot quark-gluon matter produced in heavy-ion collisions, μ A may appear due to the topological transition between the degenerate vacua of the gluonic sector of QCD characterized by different winding numbers [ 5 5 , 133 133 , 134 134 ]. The rate of such topological transition is exponentially suppressed at zero temperature as two neighbor vacua are separated by a finite energy barrier of the order of QCD confinement scale Λ QCD ∼200 MeV. However, at high temperature, the topological transition can be induced by another classical and thermal excitation called sphaleron [ 135 135 , 136 136 ] which, instead of tunneling through the barrier, can take the vacuum over the barrier at a large rate. At equilibrium, the probability of having a topological transition that drives the winding number to increase is equal to the probability of having an opposite topological transition that decreases the winding number. Thus, what really matters is the fluctuation of the topological transition above global equilibrium which results in local domains in the quark-gluon matter within which the gluonic configuration has a nonzero winding number. The quarks coupled with such gluons via the chiral anomaly can acquire nonzero chirality and thus a nonzero μ A . Furthermore, μ A in heavy-ion collisions may also be induced by other mechanisms. For example, the parallel chromoelectric and chromomagnetic fields, E a · B a ≠0, in the initial glasma can also induce a nonzero μ A [ 137 137 - 139 139 137-139 ]. 3.2 Experimental signals of the chiral magnetic effect 2 If the CME occurs in the quark-gluon matter produced in the heavy-ion collisions, what will be its experimental signals? We have already seen that the magnetic fields in heavy-ion collisions are, on average over events, perpendicular to the reaction plane; thus the CME is expected to induce an electric current which on average is also perpendicular to the reaction plane. Such an electric current will drive a charge dipole in the fireball which, after the hydrodynamic evolution of the fireball, will be converted to a charge dipole in momentum space. However, as μ A can flip the sign from event to event, the charge dipole could also flip its direction from event to event, and thus the dipolar distribution in different events cancel each other after event average and cannot be measured directly. Instead, one should look for an observable which measures the dipolar fluctuation over events. Such an observable can be chosen as the following charge-dependent two-particle azimuthal correlation [ 133 133 , 140 140 ]: <math display="block" id="M40"><mrow><msub><mi>γ</mi><mrow><mi>α</mi><mi>β</mi></mrow></msub><mo>≡</mo><mo>〈</mo><mi>cos</mi><mo stretchy="false">(</mo><msub><mi>ϕ</mi><mi>α</mi></msub><mo>+</mo><msub><mi>ϕ</mi><mi>β</mi></msub><mo>−</mo><mn>2</mn><msub><mi>Ψ</mi><mrow><mtext>R</mtext><mi>P</mi></mrow></msub><mo stretchy="false">)</mo><mo>〉</mo><mn>,</mn></mrow></math> where the indices α and β denote the charge of the hadrons, φ α and φ β are the azimuthal angles of hadrons of charge α and β , respectively, Ψ RP is the reaction plane angle, and < … > denotes average over events; see Fig. 7 Fig. 7 for illustration. It is easy to understand that a charge separation with respect to the reaction plane, will give a positive γ +- (and γ -+ ) and a negative γ ++ or γ – . We will refer to γ +- = γ -+ as γ OS and γ ++ ≈ γ – as γ SS where "OS" stands for "opposite-sign" and "SS" stands for "same-sign". In real experiments, the measurements are done with three-particle correlations where the third hadron (of arbitrary charge) is used to reconstruct the reaction plane. The correlation γ αβ was measured by STAR Collaboration at RHIC for Au + Au and Cu + Cu collisions at <math id="M41"><mrow><msqrt><mi>s</mi></msqrt><mo>=</mo><mn>200</mn></mrow></math> GeV [ 141 141 - 143 143 141-143 ] (see Fig. 8 Fig. 8 ), by PHENIX Collaboration at RHIC for Au + Au collisions at <math id="M42"><mrow><msqrt><mi>s</mi></msqrt><mo>=</mo><mn>200</mn></mrow></math> GeV [ 144 144 ], and by ALICE Collaboration at LHC for for Pb + Pb collisions at <math id="M43"><mrow><msqrt><mi>s</mi></msqrt><mo>=</mo><mn>2.76</mn></mrow></math> TeV [ 145 145 ]. At mid-central collisions, these measurements show clear positive opposite-sign correlation and negative same-sign correlation, as we expect from the picture of CME. Furthermore, the observed correlation increases from zero at central collisions to peripheral collisions. This is consistent with what we learned from the simulation in last section that the averaged magnetic field increases with the centrality. More recently, the STAR Collaboration performed the measurement of γ αβ at different beam energies [ 146 146 , 147 147 ] and observed that γ αβ persists as long as the beam energy is larger than 19.6 GeV; for further lowered beam energies, the difference between γ OS and γ SS steeply falls down, which may be understood by noticing that at lower energies the system is probably in a hadronic phase where the chiral symmetry is broken and the CME is strongly suppressed. Although the experimental data supports the presence of the CME, the interpretation of the data is not unique owing to possible background effects that are not related to the CME. First, the transverse momentum conservation (TMC) can contribute to γ αβ [ 148 148 - 150 150 148-150 ]. The momentum conservation enforces the sum of the momenta projected to the transverse plane (the plane that is perpendicular to the beam direction) of all the particles to be zero at any time during the evolution of the system. Thus, the detection of a hadron of transverse momentum, p ⊥ , requires the presence of other produced hadrons whose transverse momenta must be summed to - p ⊥ . This enforces an intrinsic back-to-back two-particle correlation: <math display="block" id="M44"><mtable columnalign="left"><mtr><mtd><msubsup><mi>γ</mi><mrow><mtext>S</mtext><mi>S</mi></mrow><mrow><mtext>TMC</mtext></mrow></msubsup><mo>=</mo><msubsup><mi>γ</mi><mrow><mtext>O</mtext><mi>S</mi></mrow><mrow><mtext>TMC</mtext></mrow></msubsup><mo>=</mo><mo>〈</mo><mfrac><mrow><mstyle displaystyle="true"><msub><mo>∑</mo><mrow><mi>i</mi><mo>≠</mo><mi>j</mi></mrow></msub><mrow><mi>cos</mi></mrow></mstyle><mo stretchy="false">(</mo><msub><mi>ϕ</mi><mi>i</mi></msub><mo>+</mo><msub><mi>ϕ</mi><mi>j</mi></msub><mo>−</mo><mn>2</mn><msub><mi>Ψ</mi><mrow><mtext>RP</mtext></mrow></msub><mo stretchy="false">)</mo></mrow><mrow><mstyle displaystyle="true"><msub><mo>∑</mo><mrow><mi>i</mi><mo>≠</mo><mi>j</mi></mrow></msub><mrow/></mstyle></mrow></mfrac><mo>〉</mo></mtd></mtr><mtr><mtd><mo>≈</mo><mo>−</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><mo>〈</mo><mi>cos</mi><mo stretchy="false">(</mo><mn>2</mn><msub><mi>ϕ</mi><mi>i</mi></msub><mo>−</mo><mn>2</mn><msub><mi>Ψ</mi><mrow><mtext>RP</mtext></mrow></msub><mo stretchy="false">)</mo><mo>〉</mo><mo>=</mo><mo>−</mo><mfrac><mrow><msub><mi>v</mi><mn>2</mn></msub></mrow><mi>N</mi></mfrac><mn>,</mn></mtd></mtr></mtable></math> where we used that <math id="M45"><mrow><mstyle displaystyle="true"><msubsup><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></msubsup><mrow><mi>cos</mi></mrow></mstyle><mo stretchy="false">(</mo><msub><mi>ϕ</mi><mi>i</mi></msub><mo>−</mo><msub><mi>Ψ</mi><mrow><mtext>RP</mtext></mrow></msub><mo stretchy="false">)</mo><mo>=</mo><mstyle displaystyle="true"><msubsup><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></msubsup><mrow><mi>sin</mi></mrow></mstyle><mo stretchy="false">(</mo><msub><mi>ϕ</mi><mi>i</mi></msub><mo>−</mo><msub><mi>Ψ</mi><mrow><mtext>RP</mtext></mrow></msub><mo stretchy="false">)</mo><mo>=</mo><mn>0</mn></mrow></math> and N is the total multiplicity. From the above relation, we can estimate how large the TMC contribution is. For example, let us consider RHIC Au + Au at 200 GeV, v 2 ∼0.1 while N ∼1000 which gives that <math id="M46"><mrow><msubsup><mi>γ</mi><mrow><mtext>S</mtext><mi>S</mi></mrow><mrow><mtext>TMC</mtext></mrow></msubsup><mo>=</mo><msubsup><mi>γ</mi><mrow><mtext>O</mtext><mi>S</mi></mrow><mrow><mtext>TMC</mtext></mrow></msubsup><mo>∼</mo><mo>−</mo><msup><mrow><mn>10</mn></mrow><mrow><mo>−</mo><mn>4</mn></mrow></msup></mrow></math> which is roughly at the same order as the experimental result shown in Fig. 8 Fig. 8 . Obviously, the TMC is charge blind. Thus we can subtract the contribution from TMC by making a difference between opposite-sign and same-sign correlations: <math display="block" id="M47"><mrow><mi>Δ</mi><mi>γ</mi><mo>≡</mo><msub><mi>γ</mi><mrow><mtext>OS</mtext></mrow></msub><mo>−</mo><msub><mi>γ</mi><mrow><mtext>SS</mtext></mrow></msub><mn>.</mn></mrow></math> However, there are other background effects contributing to Δ γ . The most notable one may be the local charge conservation (LCC) [ 148 148 , 151 151 , 152 152 ] which assumes that the pairs of opposite charges are created locally (i.e., within a small volume) in space at the late stage of the heavy-ion collisions. Then, due to the elliptic flow, opposite-sign pairs are strongly correlated while the same-sign pairs are not correlated; thus LCC gives a finite contribution to Δ γ . The detailed analysis reveals that the LCC contribute to Δ γ the amount of [ 93 93 ] <math display="block" id="M48"><mrow><mi>Δ</mi><msup><mi>γ</mi><mrow><mtext>L</mtext><mi>C</mi><mi>C</mi></mrow></msup><mo>≈</mo><mfrac><mi>M</mi><mi>N</mi></mfrac><msub><mi>v</mi><mn>2</mn></msub><mn>,</mn></mrow></math> where M is the average number of hadrons in a typical local neutral cell. As M is expected to be not very large, so we expect the LCC contribution is at the order of 10 -4 or even larger in, e.g., RHIC Au + Au collisions at 200 GeV. This is comparable to the experimental results in Fig. 8 Fig. 8 . The background effects mask the possible CME signals. To disentangle the possible CME signal and the flow-related backgrounds, one can setup experiments to either vary the backgrounds with the signal fixed or vary the signal with the backgrounds fixed (It is also important to setup delicate analysis methods to extract the elliptic flow on the event-by-event basis [ 153 153 , 154 154 ]). The former approach can be carried out by using the prolate shape of the uraniu m nuclei [ 155 155 ]. In central U + U collisions, there would be sizable v 2 due to the prolate shape of the colliding nuclei while the magnetic field would be very tiny [ 12 12 ]. The STAR Collaboration collected 0-1% most central events from U + U collisions at <math id="M49"><mrow><msqrt><mi>s</mi></msqrt><mo>=</mo><mn>197</mn></mrow></math> GeV in 2012, and indeed found sizable v 2 while Δ γ is consistent with zero [ 146 146 ]. However, it was found later that the total multiplicity is far less correlated to the number of binary collisions than expected [ 156 156 ], so that the shape selection in experiments is very hard; see Refs. [ 118 118 , 157 157 ] for recent discussions. The latter approach was proposed to be achieved by using the collisions of the isobaric nuclei, such as <math id="M50"><mrow><msubsup><mrow/><mrow><mn>44</mn></mrow><mrow><mn>96</mn></mrow></msubsup></mrow></math> Ru and <math id="M51"><mrow><msubsup><mrow/><mrow><mn>40</mn></mrow><mrow><mn>96</mn></mrow></msubsup></mrow></math> Zr [ 155 155 , 158 158 , 159 159 ]. It is expected that Ru + Ru and Zr + Zr collisions at the same beam energy and same centrality will produce similar elliptic flow but a 10% difference in the magnetic fields. Thus if the observed γ correlation contains contribution from CME, the isobaric collisions could have a chance to see it by comparing the same γ correlation in Ru + Ru and Zr + Zr collisions. The detailed numerical simulation was reported in Ref [ 158 158 ] where the authors demonstrated that the two collision types at ps <math id="M52"><mrow><msqrt><mi>s</mi></msqrt><mo>=</mo><mn>200</mn></mrow></math> GeV have more than 10% difference in the CME signal and less than 2% difference in the elliptic-flow-driven backgrounds for the centrality range of 20-60% by assuming that the CME contribution to the γ correlation is 1/3. Such a difference is feasible in current experimental setups of heavy-ion collisions at RHIC. Therefore, the isobaric collisions at RHIC would be very valuable in isolating the chiral magnetic effect from the background sources. Recently, the CMS Collaboration at LHC reported the measurement of γ αβ in high multiplicity events in p + Pb and Pb + Pb collisions at <math id="M53"><mrow><msqrt><mi>s</mi></msqrt><mo>=</mo><mn>5.02</mn></mrow></math> TeV [ 160 160 ]. They found that the γ αβ ’s in p + Pb and in Pb + Pb as functions of multiplicity behave similarly at the same multiplicities. These results give useful information about the CME search in heavy-ion collisions. First, although the magnetic field in p + Pb can be very strong, the magnetic-field direction should not be firmly correlated to the v 2 plane, and the lifetime of the magnetic field should be shorter at higher collision energy, thus the magnetic field is not expected to drive strong γ αβ in p + Pb. The CMS measurement thus gives at least two very different possibilities. One is that, although it was not expected before, the CME still occurs in p + Pb collisions and the magnetic field, although weakly correlated to the v 2 plane, is still able to drive a finite signal which combined with the background effects gives the observed results. The other one is that there is no CME at <math id="M54"><mrow><msqrt><mi>s</mi></msqrt><mo>=</mo><mn>5.02</mn></mrow></math> TeV (or very weak CME) in both p + Pb and Pb + Pb collisions, the γ αβ in both p + Pb and Pb + Pb collisions are due to background effects. This is actually somehow compatible with the recent measurement of STAR collaboration in which the so-called H correlation in measured at various energies through the Beam Energy Scan (BES) program [ 147 147 ]. In the H correlation, the background effects are subtracted under reasonable assumptions; see details in Ref. [ 147 147 ]. If one extrapolates the H correlation of STAR Collaboration to <math id="M55"><mrow><msqrt><mi>s</mi></msqrt><mo>=</mo><mn>5.02</mn></mrow></math> TeV, one can find that H correlation is actually small there. This supports that at higher energy the γ αβ is more dominated by background contributions. Nevertheless, more measurements are definitely needed. Especially the CMS collaboration should also measure the so-called δ -correlation from which we can deduce the H correlation. In addition, the charge-dependent v 2 measurements (see next section) in p + Pb collisions will be also very important. 3.3 Chiral separation effect and chiral magnetic wave 2 By looking at the CME (10), one may wonder whether there is a quantum-anomaly induced axial current in response to the external magnetic field and under what condition such current can appear. The answer is yes, and such a current can be caused by the chiral separation effect (CSE) which is the dual effect to the CME [ 18 18 , 19 19 , 161 161 ]: <math display="block" id="M56"><mrow><msub><mi>J</mi><mi>A</mi></msub><mo>=</mo><msub><mi>σ</mi><mrow><mtext>CSE</mtext></mrow></msub><mi>B</mi><mn>,</mn></mrow></math> <math display="block" id="M57"><mrow><msub><mi>σ</mi><mrow><mtext>CSE</mtext></mrow></msub><mo>=</mo><mfrac><mrow><msup><mi>e</mi><mn>2</mn></msup></mrow><mrow><mn>2</mn><msup><mi>π</mi><mn>2</mn></msup></mrow></mfrac><msub><mi>μ</mi><mi>V</mi></msub><mn>,</mn></mrow></math> where the axial current is defined by <math id="M58"><mrow><msubsup><mi>J</mi><mi>A</mi><mi>μ</mi></msubsup><mo>=</mo><mi>e</mi><mo>〈</mo><mover accent="true"><mi>ψ</mi><mo>¯</mo></mover><msup><mi>γ</mi><mi>μ</mi></msup><msub><mi>γ</mi><mn>5</mn></msub><mi>ψ</mi><mo>〉</mo></mrow></math> , and μ V is the vector chemical potential. Note that, unlike the CME, the occurrence of the CSE does not require a parity-violating environment but a charge-conjugation odd environment. We can understand the CSE from an intuitive picture. In parallel to the discussion of CME in Sec.3.1, let us consider a massless fermionic system under a strong magnetic field, and let us focus on the lowest Landau level on which the spin of positively charged fermions are polarized to be parallel to the magnetic field (we call its direction the "up" direction and its opposite direction the "down" direction) while the spin of negatively charged fermions are polarized to be anti-parallel to the magnetic field. If a positively charged fermion on the lowest Landau level moves upward then it contributes to the RH current while if it moves downward it contributes to the LH current; thus no matter which direction it moves it contributes to an axial current along the magnetic field. Similarly, it is easy to see that the negatively charged fermions contributes an axial current in the direction opposite to the magnetic field. Now, if the system contains more positively charged fermions than negatively charged fermions, namely, if μ V > 0, the downward-moving axial current cannot compensate the upward-moving axial current and thus we get a net axial current along the direction of the magnetic field. This is the CSE current. The fermions occupying the higher Landau levels do not contribute to the CSE current as their spins are not polarized. Let us formulate the above argument in a more precise way. For simplicity, we consider a system of massless Dirac fermions with only positive charge, e , in an constant magnetic field, B , in z direction which is supposed to be so strong that all the fermions are confined in the lowest Landau level. The number densities of RH and LH fermions are given by Eq. (12). As the massless fermions on the lowest Landau level move at the speed of light in a way correlated to their spins which means that the RH fermions generate a current along z direction, J R = en R , while LH fermions generate a current along - z direction, J L =- en L . Therefore the net axial current reads <math display="block" id="M59"><mrow><msub><mi>J</mi><mi>A</mi></msub><mo>=</mo><msub><mi>J</mi><mi>R</mi></msub><mo>−</mo><msub><mi>J</mi><mi>L</mi></msub><mo>=</mo><mfrac><mrow><msup><mi>e</mi><mn>2</mn></msup><mi>B</mi></mrow><mrow><mn>4</mn><msup><mi>π</mi><mn>2</mn></msup></mrow></mfrac><mo stretchy="false">(</mo><msubsup><mi>p</mi><mi>F</mi><mi>R</mi></msubsup><mo>+</mo><msubsup><mi>p</mi><mi>F</mi><mi>L</mi></msubsup><mo stretchy="false">)</mo><mn>,</mn></mrow></math> which, once we equate <math id="M60"><mrow><msubsup><mi>p</mi><mi>F</mi><mrow><mi>R</mi><mo>/</mo><mi>L</mi></mrow></msubsup></mrow></math> with μ R/L and μ V with ( μ R + μ L )/2, leads us to Eq. (18) and Eq. (19). More rigorous derivation of CSE and various aspects regarding CSE can be found in Refs. [ 18 18 , 19 19 , 161 161 - 170 170 161-170 ]. The CME and CSE couple together the the vector and axial densities and currents in the presence of the magnetic field and lead to collective modes called chiral magnetic waves (CMWs). This is most easily seen when we express CME and CSE in the chiral basis: <math display="block" id="M61"><mrow><msub><mi>J</mi><mi>R</mi></msub><mo>=</mo><mfrac><mrow><msup><mi>e</mi><mn>2</mn></msup></mrow><mrow><mn>4</mn><msup><mi>π</mi><mn>2</mn></msup></mrow></mfrac><msub><mi>μ</mi><mi>R</mi></msub><mi>B</mi><mn>,</mn></mrow></math> <math display="block" id="M62"><mrow><msub><mi>J</mi><mi>L</mi></msub><mo>=</mo><mo>−</mo><mfrac><mrow><msup><mi>e</mi><mn>2</mn></msup></mrow><mrow><mn>4</mn><msup><mi>π</mi><mn>2</mn></msup></mrow></mfrac><msub><mi>μ</mi><mi>L</mi></msub><mi>B</mi><mn>,</mn></mrow></math> where the RH and LH quantities are defined as J R/L = J V ± J A )/2 and μ R/L = μ V ± μ A . Substituting J R and J L in to the continuation equations and considering a small fluctuation in the density <math id="M63"><mrow><msubsup><mi>J</mi><mrow><mi>R</mi><mo>/</mo><mi>L</mi></mrow><mn>0</mn></msubsup></mrow></math> , we obtain <math display="block" id="M64"><mrow><msub><mo>∂</mo><mi>t</mi></msub><mi>δ</mi><msubsup><mi>J</mi><mrow><mi>R</mi><mo>/</mo><mi>L</mi></mrow><mn>0</mn></msubsup><mo>±</mo><mfrac><mrow><msup><mi>e</mi><mn>2</mn></msup></mrow><mrow><mn>4</mn><msup><mi>π</mi><mn>2</mn></msup><mi>χ</mi></mrow></mfrac><mi>B</mi><mo>⋅</mo><mo>∇</mo><mi>δ</mi><msubsup><mi>J</mi><mrow><mi>R</mi><mo>/</mo><mi>L</mi></mrow><mn>0</mn></msubsup><mo>=</mo><mn>0,</mn></mrow></math> where <math id="M65"><mrow><mi>χ</mi><mo>=</mo><mo>∂</mo><msup><mi>J</mi><mn>0</mn></msup><mo>/</mo><mo>∂</mo><mi>μ</mi></mrow></math> is the number susceptibility (we assume that the RH and LH chiralities have the same susceptibility). Therefore, we can identify two collective wave modes, one propagating along B and one propagating opposite to the direction of B with the same velocity v χ = e 2 B /(4 π 2 χ ). They are the CMWs [ 20 20 ] (see also Ref. [ 171 171 ] for a derivation of CMW based on the kinetic theory). 3.4 Experimental signals for chiral magnetic wave 2 The CMWs are able to transport both chirality and electric charge and can induce an electric quadrupole in the QGP [ 21 21 , 22 22 ]; see Fig. 9 Fig. 9 for an illustration. Owing to strong in-plane pressure gradient, such an electric quadrupole can cause a splitting between the v 2 of π + and π - [ 22 22 ]; see also Refs. [ 139 139 , 172 172 - 174 174 172-174 ] for more advanced numerical implementations. The CMW-induced charge dependent elliptic flows can be parameterized as <math display="block" id="M66"><mrow><msub><mi>v</mi><mn>2</mn></msub><mo stretchy="false">(</mo><msup><mi>π</mi><mo>±</mo></msup><mo stretchy="false">)</mo><mo>=</mo><msub><mi>v</mi><mn>2</mn></msub><mo>∓</mo><mfrac><mrow><mi>r</mi><msub><mi>A</mi><mrow><mtext>c</mtext><mi>h</mi></mrow></msub></mrow><mn>2</mn></mfrac><mn>,</mn></mrow></math> where v 2 is the average elliptic flow and A ch =( N + - N - )/( N + + N - ) is the net charge asymmetry parameter which is usually small (the events that can generate large A ch are rare). The slop parameter, r , is predicted to be proportional to < ( eB ) 2 cos[2( ψ B -Ψ RP )] > where ψ B is the azimuthal direction of the magnetic field [ 11 11 , 22 22 , 170 170 ]. Recently, the CMW has been tested by STAR Collaboration [ 146 146 , 175 175 , 176 176 ] at RHIC and by ALICE Collaboration [ 177 177 ] at LHC. The data shows an elliptic-flow difference, v 2 ( π - )- v 2 ( π + ), linear in A ch with a positive slope whose centrality dependence can be fitted by the CMW estimation. Furthermore, the slope parameter, r , displays no obvious trend of the beam energy dependence for 10-60% centrality at <math id="M67"><mrow><msqrt><mi>s</mi></msqrt><mo>=</mo><mn>20</mn><mo>−</mo><mn>200</mn></mrow></math> GeV. This can be understood by noticing that the strength of magnetic field times is roughly beam energy independent ( Fig. 10 Fig. 10 ). However, ambiguities exist in the interpretation of the experimental results as there are background effects which can also lead to v 2 splitting between π + and π - . A variety of possible background effects were proposed, see Refs. [ 10 10 , 178 178 - 183 183 178-183 ] and review articles, Refs. [ 93 93 , 126 126 ], for discussions. It should be noted that none of the above backgrounds can successfully explain all the features of the experimental data. It is quite plausible that various effects (including CMW) are coupled to produce the measured pattern. Until we can successfully subtract the background effects, we cannot make a conclusive claim. On the experimental side, it will be important to perform more detailed measurements, such as the measurement of the v 2 splitting versus transverse momentum and rapidity, the measurement for other hadrons like kaons, the measurement of splitting in other harmonic flows like v 3 , and the measurement of the cross-correlation of different observables. On the theoretical side, an anomalous hydrodynamic simulation or a simulation based on kinetic theory for the CMW that incorporate both the dynamically evolving magnetic fields and the realistic initial charge distribution will be very useful. Before we end this section, we note that, in addition to the magnetic-field induced anomalous transports, the electric field and vorticity can also induce certain anomalous transport phenomena and corresponding collective modes. For example, electric field can induce an axial current in the medium with both μ A and μ V nonzero which is called chiral electric separation effect (CESE) and, once coupled with the normal Ohm’s current (the vector current induced by electric field), CESE leads to chiral electric wave and axial and vecor density waves [ 184 184 - 187 187 184-187 ]. The CESE may lead to observable effect in, Cu + Au collisions where persistant electric fields from Au to Cu nuclei exist [ 188 188 ]. The fluid vorticity can induce the so-called chiral vortical effects (CVEs), which are the vortical analogues of CME [ 189 189 - 192 192 189-192 ], and the corresponding collective mode, chiral vortical wave (CVW) [ 193 193 ]. The CVE and CVW may have important consequences in noncentral heavy-ion collisions where nonzero fluid vorticity can be generated as a result of angular momentum conservation [ 194 194 - 197 197 194-197 ]. Phenomenologically, the CVE can induce baryon number separation with respect to the reaction plane [ 198 198 ], and CVW may cause a v 2 splitting between Λ and <math id="M68"><mover accent="true"><mi>Λ</mi><mo>¯</mo></mover></math> [ 193 193 ]. Recently, the STAR Collaboration has reported the first measurement of CVE in Au + Au collisions at <math id="M69"><mrow><msqrt><mi>s</mi></msqrt><mo>=</mo><mn>200</mn></mrow></math> GeV [ 199 199 ], and the data is consistent with the expectation of CVE. But, similar with the case of CME and CMW detections, the data contains contributions from background effects like the transverse momentum conservation and the local baryon number conservation, and we need more efforts to understand the data. 4 Heavy-quark dynamics in magnetic fields 1 In the heavy ion collisions, the dominant process for the heavy-quark production is the initial hard scatterings among the partons from the colliding nuclei. This is because the thermal excitation is suppressed due to the large value of the heavy-quark mass compared to the temperature of the quark-gluon plasma (QGP). Therefore, the heavy quarks will serve as a probe of the dynamics in the magnetic field persisting in the early time of the collision events. In this section, we discuss effects of the magnetic field on the quarkonium spectra in Sec. 4.1 and the heavy-quark transport in a hot medium and magnetic field in Sec. 4.2. The quarkonium spectra in magnetic fields have been investigated by using the Cornell potential model [ 27 27 , 33 33 , 38 38 , 39 39 ], the QCD sum rules [ 26 26 , 29 29 , 30 30 ] with the help of the lattice QCD simulations for the heavy quark potential [ 32 32 , 34 34 ], and the quark and gluon condensates [ 48 48 - 50 50 48-50 , 200 200 , 201 201 ]. The phenomenological consequences for the quarkonium production in the heavy-ion collisions have been discussed in Ref. [ 23 23 - 25 25 23-25 , 27 27 , 28 28 , 37 37 ], which could be measured by the dilepton channel. Most of the heavy quarks created in the hard scatterings will evolve in a hot medium and be observed as open heavy flavor mesons in the final state (see Ref. [ 202 202 ] for a review). Therefore, it is important to understand the transport property of the heavy quarks in a hot medium and magnetic field. In Sec. 4.2, we discuss the diagrammatic computation of the heavy-quark diffusion constant. The study on the heavy-quark dynamics in the magnetized medium has been just initiated [ 40 40 , 42 42 ] and is awaiting for further developments. 4.1 Quarkonia in magnetic fields 2 4.1.1 Spin mixing and level repulsion 3 We shall begin with the simplest description of effects of magnetic fields on quarkonium spectra. In a weak magnetic field, the dominant effects of the magnetic fields can be described by the hadronic degrees of freedom. This contrasts the case of strong magnetic fields where modifications of nonperturbative vacuum structure and of the internal structure of hadrons are important (see Refs. [ 58 58 , 52 52 , 62 62 , 68 68 - 70 70 68-70 , 203 203 - 207 207 203-207 ] and references therein for the ongoing discussions). As long as the magnitude of the magnetic field is much smaller than the heavy quark mass, these effects will be subdominant contributions. In the hadronic degrees of freedom, effects of magnetic fields appear in the mixing of spin eigenstates [ 25 25 , 27 27 , 29 29 , 30 30 , 33 33 , 38 38 , 39 39 ]. This mixing is caused by the breaking of the spatial rotation symmetry in magnetic fields: there remains only the azimuthal rotational symmetry with respect to the direction of the external magnetic field. This indicates that only the spin state along the magnetic field can persist as a good quantum number of the mesons, and thus that there is a mixing between the spin single and triplet states, S total , S z ) = (0,0) and (1,0), which are, for example, η c and the longitudinal component of J/ψ . We examine general mixing patterns among the pseudoscalar, P , vector, V μ , scalar, S , and axial-vector, A μ , quarkonia. Possible interaction vertices among those fields are informed from the Lorentz invariance and the parity and charge-conjugation symmetries. The vertices relevant for interactions among static charmonia are found to be [ 29 29 , 30 30 ] <math display="block" id="M70"><mrow><msub><mi mathvariant="script">L</mi><mrow><msub><mi>γ</mi><mrow><mtext>P</mtext><mi>V</mi></mrow></msub></mrow></msub><mo>=</mo><mfrac><mrow><msub><mi>g</mi><mrow><msub><mrow/><mrow><mtext>P</mtext><mi>V</mi></mrow></msub></mrow></msub></mrow><mrow><msub><mi>m</mi><mn>0</mn></msub></mrow></mfrac><mi>e</mi><msub><mover accent="true"><mi>F</mi><mo>˜</mo></mover><mrow><mi>μ</mi><mi>ν</mi></mrow></msub><mo stretchy="false">(</mo><msup><mo>∂</mo><mi>μ</mi></msup><mi>P</mi><mo stretchy="false">)</mo><msup><mi>V</mi><mi>ν</mi></msup><mtext> </mtext><mn>,</mn></mrow></math> <math display="block" id="M71"><mrow><msub><mi mathvariant="script">L</mi><mrow><msub><mi>γ</mi><mrow><mtext>V</mtext><mi>A</mi></mrow></msub></mrow></msub><mo>=</mo><mi>i</mi><msub><mi>g</mi><mrow><msub><mrow/><mrow><mtext>V</mtext><mi>A</mi></mrow></msub></mrow></msub><mi>e</mi><msub><mover accent="true"><mi>F</mi><mo>˜</mo></mover><mrow><mi>μ</mi><mi>ν</mi></mrow></msub><msup><mi>V</mi><mi>μ</mi></msup><msup><mi>A</mi><mi>ν</mi></msup><mtext> </mtext><mn>,</mn></mrow></math> <math display="block" id="M72"><mrow><msub><mi mathvariant="script">L</mi><mrow><msub><mi>γ</mi><mrow><mtext>S</mtext><mi>A</mi></mrow></msub></mrow></msub><mo>=</mo><mfrac><mrow><msub><mi>g</mi><mrow><msub><mrow/><mrow><mtext>S</mtext><mi>A</mi></mrow></msub></mrow></msub></mrow><mrow><msub><mi>m</mi><mn>1</mn></msub></mrow></mfrac><mi>e</mi><msub><mover accent="true"><mi>F</mi><mo>˜</mo></mover><mrow><mi>μ</mi><mi>ν</mi></mrow></msub><mo stretchy="false">(</mo><msup><mo>∂</mo><mi>μ</mi></msup><mi>S</mi><mo stretchy="false">)</mo><msup><mi>A</mi><mi>ν</mi></msup><mtext> </mtext><mn>,</mn></mrow></math> with m 0 =( m P + m V )/2, m 1 =( m S + m A )/2, and dimensionless effective coupling constants g PV , g VA , and g SA . The vertices proportional to <math id="M73"><mrow><msup><mi>F</mi><mrow><mi>μ</mi><mi>ν</mi></mrow></msup><msub><mo>∂</mo><mi>μ</mi></msub></mrow></math> do not play any role for static charmonia, because they pick up the vanishing spatial derivatives and the temporal component of the vector or axial-vector fields. The nonvanishing vertices are responsible for, radiative decay modes of quarkonia such as J / ψ → η c + γ [ 208 208 , 209 209 ]. It is also worth mentioning that the three-point vertex (25) induces the conversion of η c to J/ψ in a magnetic field. The contribution to the J/ψ production from this conversion process was examined in Ref. [ 24 24 ]. From the above vertices, one can identify two pairs of mixing partners: the mixing between P and V and between S and A are possible when they are at rest in external magnetic fields. Since the P and V are not mixed with the S and A , we shall focus on the mixing between S and V . As mentioned in the beginning of this section, a mixing effect arises only in the longitudinal component of the vector field, V || , and the other two transverse components persist as the energy eigenstates. This can be easily confirmed by identifying the non-vanishing components in Eq. (25). With a vanishing spatial momentum q μ =( ω , 0, 0, 0) and the vertex (25), the coupled equations of motion are found to be <math display="block" id="M74"><mrow><mrow><mo>(</mo><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mo>−</mo><msup><mi>ω</mi><mn>2</mn></msup><mo>+</mo><msubsup><mi>m</mi><mtext>P</mtext><mn>2</mn></msubsup></mrow></mtd><mtd columnalign="left"><mrow><mo>−</mo><mi>i</mi><mfrac><mrow><msub><mi>g</mi><mrow><msub><mrow/><mrow><mtext>P</mtext><mi>V</mi></mrow></msub></mrow></msub></mrow><mrow><msub><mi>m</mi><mn>0</mn></msub></mrow></mfrac><mi>ω</mi><mi>e</mi><mi>B</mi></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow><mi>i</mi><mfrac><mrow><msub><mi>g</mi><mrow><msub><mrow/><mrow><mtext>P</mtext><mi>V</mi></mrow></msub></mrow></msub></mrow><mrow><msub><mi>m</mi><mn>0</mn></msub></mrow></mfrac><mi>ω</mi><mi>e</mi><mi>B</mi></mrow></mtd><mtd columnalign="left"><mrow><mo>−</mo><msup><mi>ω</mi><mn>2</mn></msup><mo>+</mo><msubsup><mi>m</mi><mtext>V</mtext><mn>2</mn></msubsup></mrow></mtd></mtr></mtable></mrow><mo>)</mo></mrow><mrow><mo>(</mo><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mi>P</mi></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow><msub><mi>V</mi><mo>∥</mo></msub></mrow></mtd></mtr></mtable></mrow><mo>)</mo></mrow><mo>=</mo><mn>0.</mn></mrow></math> Following from the equations of motion (28), we obtain the physical mass eigenvalues as <math display="block" id="M75"><mrow><msubsup><mover accent="true"><mi>m</mi><mo>¯</mo></mover><mrow><mi>V</mi><mn>,</mn><mi>P</mi></mrow><mn>2</mn></msubsup><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo stretchy="false">(</mo><msubsup><mi>M</mi><mo>+</mo><mn>2</mn></msubsup><mo>+</mo><mfrac><mrow><msup><mi>γ</mi><mn>2</mn></msup></mrow><mrow><msubsup><mi>m</mi><mn>0</mn><mn>2</mn></msubsup></mrow></mfrac><mo>±</mo><msqrt><mrow><msubsup><mi>M</mi><mo>−</mo><mn>4</mn></msubsup><mo>+</mo><mfrac><mrow><mn>2</mn><msup><mi>γ</mi><mn>2</mn></msup><msubsup><mi>M</mi><mo>+</mo><mn>2</mn></msubsup></mrow><mrow><msubsup><mi>m</mi><mn>0</mn><mn>2</mn></msubsup></mrow></mfrac><mo>+</mo><mfrac><mrow><msup><mi>γ</mi><mn>4</mn></msup></mrow><mrow><msubsup><mi>m</mi><mn>0</mn><mn>4</mn></msubsup></mrow></mfrac></mrow></msqrt><mo stretchy="false">)</mo><mn>,</mn></mrow></math> where <math id="M76"><mrow><msubsup><mi>M</mi><mo>+</mo><mn>2</mn></msubsup><mo>=</mo><msubsup><mi>m</mi><mtext>P</mtext><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>m</mi><mtext>V</mtext><mn>2</mn></msubsup><mn>,</mn><msubsup><mi>M</mi><mo>−</mo><mn>2</mn></msubsup><mo>=</mo><msubsup><mi>m</mi><mtext>V</mtext><mn>2</mn></msubsup><mo>−</mo><msubsup><mi>m</mi><mtext>P</mtext><mn>2</mn></msubsup></mrow></math> , and γ = g PV eB . The upper and lower signs are for the vector and pseudoscalar states, respectively. Expanding the right-hand side of Eq. (29) up to the second order in γ and the leading order in <math id="M77"><mrow><mfrac><mn>1</mn><mn>2</mn></mfrac><mo stretchy="false">(</mo><msub><mi>m</mi><mtext>V</mtext></msub><mo>−</mo><msub><mi>m</mi><mtext>P</mtext></msub><mo stretchy="false">)</mo><mo>/</mo><msub><mi>m</mi><mn>0</mn></msub></mrow></math> , we find <math display="block" id="M78"><mrow><msubsup><mover accent="true"><mi>m</mi><mo>¯</mo></mover><mrow><mi>P</mi><mn>,</mn><mi>V</mi></mrow><mn>2</mn></msubsup><mo stretchy="false">(</mo><mi>B</mi><mo stretchy="false">)</mo><mo>=</mo><msubsup><mi>m</mi><mrow><mi>V</mi><mn>,</mn><mi>P</mi></mrow><mn>2</mn></msubsup><mo>±</mo><mfrac><mrow><msup><mi>γ</mi><mn>2</mn></msup></mrow><mrow><msubsup><mi>M</mi><mo>−</mo><mn>2</mn></msubsup></mrow></mfrac><mn>,</mn></mrow></math> with eigenvectors given by <math display="block" id="M79"><mtable columnalign="left"><mtr><mtd><mo>|</mo><mi>P</mi><msub><mo>〉</mo><mi>B</mi></msub><mtext> </mtext><mtext> </mtext><mo>=</mo><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mfrac><mrow><msup><mi>γ</mi><mn>2</mn></msup></mrow><mrow><msubsup><mi>M</mi><mo>−</mo><mn>4</mn></msubsup></mrow></mfrac><mo stretchy="false">)</mo><mo>|</mo><mi>P</mi><mo>〉</mo><mo>−</mo><mi>i</mi><mfrac><mi>γ</mi><mrow><msubsup><mi>M</mi><mo>−</mo><mn>2</mn></msubsup></mrow></mfrac><mo>|</mo><mi>V</mi><mo>〉</mo><mn>,</mn></mtd></mtr><mtr><mtd><mo>|</mo><mi>V</mi><msub><mo>〉</mo><mi>B</mi></msub><mtext> </mtext><mo>=</mo><mo>−</mo><mi>i</mi><mfrac><mi>γ</mi><mrow><msubsup><mi>M</mi><mo>−</mo><mn>2</mn></msubsup></mrow></mfrac><mo>|</mo><mi>P</mi><mo>〉</mo><mo>+</mo><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mfrac><mrow><msup><mi>γ</mi><mn>2</mn></msup></mrow><mrow><msubsup><mi>M</mi><mo>−</mo><mn>4</mn></msubsup></mrow></mfrac><mo stretchy="false">)</mo><mo>|</mo><mi>V</mi><mo>〉</mo><mtext> </mtext><mn>.</mn></mtd></mtr></mtable></math> The mass eigenvalue (30) clearly shows the level repulsion in the presence of the mixing effect. In case of mixing between η c and J/ψ , the coupling constant, g PV , can be read off from the measured radiative decay width, J / ψ → η c + γ [ 29 29 , 30 30 ]. From the vertex (25), the invariant amplitude of this decay process is written as <math display="block" id="M80"><mtable columnalign="left"><mtr><mtd><msub><mi mathvariant="script">M</mi><mrow><msub><mrow/><mrow><mtext>P</mtext><mi>V</mi></mrow></msub></mrow></msub><mo>=</mo><mo>〈</mo><mi>γ</mi><mtext>P</mtext><mo>|</mo><msub><mi mathvariant="script">L</mi><mrow><mi>γ</mi><mtext>P</mtext><mi>V</mi></mrow></msub><mo>|</mo><mtext>V</mtext><mo>〉</mo></mtd></mtr><mtr><mtd><mo>=</mo><mo>−</mo><mfrac><mrow><mi>e</mi><msub><mi>g</mi><mrow><msub><mrow/><mrow><mtext>P</mtext><mi>V</mi></mrow></msub></mrow></msub></mrow><mrow><msub><mi>m</mi><mn>0</mn></msub></mrow></mfrac><msub><mi>ϵ</mi><mrow><mi>μ</mi><mi>ν</mi><mi>α</mi><mi>β</mi></mrow></msub><msubsup><mi>k</mi><mi>γ</mi><mi>μ</mi></msubsup><msubsup><mi>ϵ</mi><mi>γ</mi><mi>ν</mi></msubsup><msubsup><mi>p</mi><mtext>V</mtext><mi>α</mi></msubsup><msubsup><mi>ϵ</mi><mtext>V</mtext><mi>β</mi></msubsup><mtext> </mtext><mtext> </mtext><mn>,</mn></mtd></mtr></mtable></math> where <math id="M81"><mrow><msubsup><mi>k</mi><mi>γ</mi><mi>μ</mi></msubsup></mrow></math> and <math id="M82"><mrow><msubsup><mi>ϵ</mi><mi>γ</mi><mi>ν</mi></msubsup></mrow></math> ( <math id="M83"><mrow><msubsup><mi>p</mi><mtext>V</mtext><mi>α</mi></msubsup></mrow></math> and <math id="M84"><mrow><msubsup><mi>ϵ</mi><mtext>V</mtext><mi>β</mi></msubsup></mrow></math> ) denote the polarization and momentum of the photon ( J/ψ ). Summing the polarizations of the photon and averaging those of the vector meson, we find <math display="block" id="M85"><mrow><mi>Γ</mi><mo stretchy="false">[</mo><mtext>V</mtext><mo>→</mo><mi>γ</mi><mtext>P</mtext><mo stretchy="false">]</mo><mo>=</mo><mfrac><mn>1</mn><mrow><mn>12</mn></mrow></mfrac><mfrac><mrow><msup><mi>e</mi><mn>2</mn></msup><msubsup><mi>g</mi><mrow><msub><mrow/><mrow><mtext>P</mtext><mi>V</mi></mrow></msub></mrow><mn>2</mn></msubsup><msup><mover accent="true"><mi>p</mi><mo>˜</mo></mover><mn>3</mn></msup></mrow><mrow><mi>π</mi><msubsup><mi>m</mi><mn>0</mn><mn>2</mn></msubsup></mrow></mfrac><mn>.</mn><mtext> </mtext></mrow></math> We have also integrated over the phase-space volume in the two-body final state. The center-of-mass momentum in the final state is given by <math id="M86"><mrow><mover accent="true"><mi>p</mi><mo>˜</mo></mover><mo>=</mo><mo stretchy="false">(</mo><msubsup><mi>m</mi><mtext>V</mtext><mn>2</mn></msubsup><mo>−</mo><msubsup><mi>m</mi><mtext>P</mtext><mn>2</mn></msubsup><mo stretchy="false">)</mo><mo>/</mo><mo stretchy="false">(</mo><mn>2</mn><msub><mi>m</mi><mtext>V</mtext></msub><mo stretchy="false">)</mo></mrow></math> . Using the measured radiative decay width, Γ exp [ J / ψ → γη c ]=1.579 keV, the coupling strength is estimated to be g PV =2.095. A consistent value of the coupling strength is also obtained from computation in the quark degrees of freedom with help of the heavy quark expansion and diagrammatic methods in an external magnetic field (see Appendix. 6). Figure 11 Figure 11 shows the mass shifts in the presence of the mixing effect. We find the level repulsion between the η c and longitudinal J/ψ mass spectra. While the mass of η c decreases as eB increases, the mass of the longitudinal J/ψ (denoted by J / ψ ∥ ) increases. This result clearly provides the simple picture of the effect of weak magnetic fields. The level repulsion occurs most strongly between a pair of adjacent states. Therefore, the other excited states give rise to only subdominant contributions to the mass shifts of η c and the longitudinal J/ψ in a weak magnetic field. However, as the mass of the longitudinal J/ψ keeps increasing with an increasing magnetic field, it will approaches the other excited states. In such a strong magnetic field, many states will get involved in the mixing effects, which will result in an intriguing mixing patterns of the spectra as discussed in the next section. Moreover, the mixing effects induced by the minimal coupling (25) to a magnetic field will not be valid in the strong-field regime. Since the strong magnetic field probes the inner structure of hadrons and modifies the QCD vacuum, we will need to look into the dynamics of the fundamental degrees of freedom. To investigate these effects, one can invoke on the Cornell potential models and QCD sum rules with the help of lattice QCD simulations. We will proceed to exploring these frameworks on top of the simplest framework for the mixing effect discussed in this section. 4.1.2 Cornell potential model in magnetic fields 3 To study the nonrelativistic bound states of the heavy quark and antiquark on the basis of the fundamental degrees of freedom, one can use the Cornell potential model [ 210 210 , 211 211 ] with the help of the lattice QCD simulations. The basic framework for the 1 S states of η c , J/ψ , η b , and Υ, was first elaborated by Alford and Strickland [ 27 27 ]. The authors included the interaction with the background magnetic field through the magnetic moments of constituents and the kinetic terms, as well as the mutual interactions by the linear and Coulomb potentials in vacuum. In this framework, the interactions through the magnetic moments give rise to the mixing between the spin singlet and triplet states discussed in the previous section. Based on this model, a systematic spectroscopy for excited states was recently performed by Suzuki and Yoshida owing to an improvement of computational method [ 38 38 , 39 39 ]. Effects of strong magnetic fields on nonperturbative QCD vacuum were investigated by Bonati et al. who computed the modification of the quark-antiquark potential by lattice QCD simulation [ 32 32 , 34 34 ] and plugged it into the bound-state problem [ 33 33 ]. These studies suggest an opportunity to explore the nature of QCD confinement in the presence of a strong magnetic field. In the following, we briefly summarize the framework for the potential model in a magnetic field and the numerical results for the mass spectra. The Hamiltonian for the two-body bound state in a magnetic field is given by [ 25 25 , 27 27 , 38 38 , 39 39 , 33 33 ] <math display="block" id="M87"><mrow><mi mathvariant="script">H</mi><mo>=</mo><mn>2</mn><msub><mi>m</mi><mi>Q</mi></msub><mo>+</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1,2</mn></mrow></munder></mstyle><msubsup><mi mathvariant="script">H</mi><mn>0</mn><mrow><mo stretchy="false">(</mo><mi>i</mi><mo stretchy="false">)</mo></mrow></msubsup><mo>+</mo><msub><mi mathvariant="script">H</mi><mrow><mover accent="true"><mi>q</mi><mo>¯</mo></mover><mi>q</mi></mrow></msub><mo>+</mo><msub><mi mathvariant="script">H</mi><mi>μ</mi></msub></mrow></math> where <math display="block" id="M88"><mrow><msubsup><mi mathvariant="script">H</mi><mn>0</mn><mrow><mo stretchy="false">(</mo><mi>i</mi><mo stretchy="false">)</mo></mrow></msubsup><mo>=</mo><mfrac><mn>1</mn><mrow><mn>2</mn><msub><mi>m</mi><mi>Q</mi></msub></mrow></mfrac><msup><mrow><mrow><mo>(</mo><mrow><mtext> </mtext><msup><mi>p</mi><mrow><mo stretchy="false">(</mo><mi>i</mi><mo stretchy="false">)</mo></mrow></msup><mo>−</mo><msup><mi>q</mi><mrow><mo stretchy="false">(</mo><mi>i</mi><mo stretchy="false">)</mo></mrow></msup><mi>A</mi><mo stretchy="false">(</mo><msup><mi>r</mi><mrow><mo stretchy="false">(</mo><mi>i</mi><mo stretchy="false">)</mo></mrow></msup><mo stretchy="false">)</mo><mtext> </mtext></mrow><mo>)</mo></mrow></mrow><mn>2</mn></msup></mrow></math> <math display="block" id="M89"><mrow><msub><mi mathvariant="script">H</mi><mrow><mover accent="true"><mi>q</mi><mo>¯</mo></mover><mi>q</mi></mrow></msub><mo>=</mo><msub><mi>σ</mi><mi>s</mi></msub><mi>r</mi><mo>−</mo><mfrac><mi>a</mi><mi>r</mi></mfrac><mo>+</mo><mi>b</mi><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo><mrow><mo>(</mo><mrow><msup><mi>σ</mi><mrow><mo stretchy="false">(</mo><mn>1</mn><mo stretchy="false">)</mo></mrow></msup><mo>⋅</mo><msup><mi>σ</mi><mrow><mo stretchy="false">(</mo><mn>2</mn><mo stretchy="false">)</mo></mrow></msup></mrow><mo>)</mo></mrow></mrow></math> <math display="block" id="M90"><mrow><msub><mi mathvariant="script">H</mi><mi>μ</mi></msub><mo>=</mo><mo>−</mo><mrow><mo>(</mo><mrow><msup><mi>μ</mi><mrow><mo stretchy="false">(</mo><mn>1</mn><mo stretchy="false">)</mo></mrow></msup><mo>+</mo><msup><mi>μ</mi><mrow><mo stretchy="false">(</mo><mn>2</mn><mo stretchy="false">)</mo></mrow></msup></mrow><mo>)</mo></mrow><mo>⋅</mo><mi>Β</mi></mrow></math> with the gauge potential, A ( r ), for the external magnetic field, B , the Pauli matrix, σ , and the heavy-quark mass, m Q . The indices i =1,2 denote the quark and antiquark, respectively, and thus q (1) =- q (2) (≡ q ) for quarkonia. The mutual interaction is described by the Cornell potential, <math id="M91"><mrow><msub><mi mathvariant="script">H</mi><mrow><mover accent="true"><mi>q</mi><mo>¯</mo></mover><mi>q</mi></mrow></msub></mrow></math> , which includes the linear and Coulomb potentials as well as the spin-spin interaction. For the nonrelativistic quark and antiquark, the magnetic moments read <math display="block" id="M92"><mrow><msub><mi mathvariant="script">H</mi><mi>μ</mi></msub><mo>=</mo><mo>−</mo><msub><mi>μ</mi><mi>Q</mi></msub><mfrac><mi>g</mi><mn>2</mn></mfrac><mo stretchy="false">(</mo><msubsup><mi>σ</mi><mi>z</mi><mrow><mo stretchy="false">(</mo><mn>1</mn><mo stretchy="false">)</mo></mrow></msubsup><mo>−</mo><msubsup><mi>σ</mi><mi>z</mi><mrow><mo stretchy="false">(</mo><mn>2</mn><mo stretchy="false">)</mo></mrow></msubsup><mo stretchy="false">)</mo><mi>B</mi></mrow></math> where μ Q = q /(2 m Q ) and <math id="M93"><mrow><mi>g</mi><mo>≃</mo><mn>2</mn></mrow></math> . We took the magnetic field in + z direction without loss of generality. The eigenstates of the Pauli matrices are spanned by the direct product <math id="M94"><mrow><mo>|</mo><mo>±</mo><mo>±</mo><mo>〉</mo><mo>≡</mo><mo>|</mo><mo>±</mo><mo>〉</mo><mo>⊗</mo><mo>|</mo><mo>±</mo><mo>〉</mo></mrow></math> where the left and right states in the product stand for the eigenstates of <math id="M95"><mrow><msubsup><mi>σ</mi><mi>z</mi><mrow><mo stretchy="false">(</mo><mn>1</mn><mo stretchy="false">)</mo></mrow></msubsup></mrow></math> and <math id="M96"><mrow><msubsup><mi>σ</mi><mi>z</mi><mrow><mo stretchy="false">(</mo><mn>2</mn><mo stretchy="false">)</mo></mrow></msubsup></mrow></math> , respectively. In the total-spin basis, the spin singlet state is given by <math id="M97"><mrow><mo>|</mo><mn>0,0</mn><mo>〉</mo><mo>=</mo><mo stretchy="false">(</mo><mo>|</mo><mo>+</mo><mo>−</mo><mo>〉</mo><mo>−</mo><mo>|</mo><mo>−</mo><mo>+</mo><mo>〉</mo><mo stretchy="false">)</mo><mo>/</mo><msqrt><mn>2</mn></msqrt></mrow></math> , while the triplet states by is |1,±1〉=|±±〉 and <math id="M98"><mrow><mo>|</mo><mn>1,0</mn><mo>〉</mo><mo>=</mo><mo stretchy="false">(</mo><mo>|</mo><mo>+</mo><mo>−</mo><mo>〉</mo><mo>+</mo><mo>|</mo><mo>−</mo><mo>+</mo><mo>〉</mo><mo stretchy="false">)</mo><mo>/</mo><msqrt><mn>2</mn></msqrt></mrow></math> . Therefore, the operation of <math id="M99"><mrow><msub><mi mathvariant="script">H</mi><mi>μ</mi></msub></mrow></math> on these states results in <math display="block" id="M100"><mrow><msub><mi mathvariant="script">H</mi><mi>μ</mi></msub><mo>|</mo><mn>1,</mn><mo>+</mo><mn>1</mn><mo>〉</mo><mo>=</mo><msub><mi mathvariant="script">H</mi><mi>μ</mi></msub><mo>|</mo><mn>1,</mn><mo>−</mo><mn>1</mn><mo>〉</mo><mo>=</mo><mn>0</mn><mtext> </mtext><mn>,</mn></mrow></math> <math display="block" id="M101"><mrow><msub><mi mathvariant="script">H</mi><mi>μ</mi></msub><mrow><mo>(</mo><mrow><mtext>​</mtext><mtext>​</mtext><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mo>|</mo><mn>0,0</mn><mo>〉</mo></mrow></mtd><mtd columnalign="left"><mrow/></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow><mo>|</mo><mn>1,0</mn><mo>〉</mo></mrow></mtd><mtd columnalign="left"><mrow/></mtd></mtr></mtable><mtext>​</mtext><mtext>​</mtext></mrow><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mn>0</mn></mtd><mtd columnalign="left"><mrow><mn>2</mn><msub><mi>μ</mi><mi>Q</mi></msub><mi>B</mi></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow><mn>2</mn><msub><mi>μ</mi><mi>Q</mi></msub><mi>B</mi></mrow></mtd><mtd columnalign="left"><mn>0</mn></mtd></mtr></mtable></mrow><mo>)</mo></mrow><mtext>​</mtext><mtext>​</mtext><mrow><mo>(</mo><mrow><mtext>​</mtext><mtext>​</mtext><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mo>|</mo><mn>0,0</mn><mo>〉</mo></mrow></mtd><mtd columnalign="left"><mrow/></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow><mo>|</mo><mn>1,0</mn><mo>〉</mo></mrow></mtd><mtd columnalign="left"><mrow/></mtd></mtr></mtable><mtext>​</mtext><mtext>​</mtext></mrow><mo>)</mo></mrow><mtext> </mtext><mn>.</mn></mrow></math> Consistent with the discussion in the previous section, we find that, while the transverse modes |1,±1〉 are eigenstates of the two-body Hamiltonian, there is a mixing between |0, 0〉 and |1, 0〉 states which need to be diagonalized to the energy eigenstates. A specific issue of the dynamics in magnetic fields is the non-conservation of the transverse momentum. The gauge potential for an external magnetic field will depend on the transverse coordinate as <math id="M102"><mrow><msubsup><mi>A</mi><mrow><mtext>e</mtext><mi>x</mi><mi>t</mi></mrow><mi>μ</mi></msubsup><mo>=</mo><mo stretchy="false">(</mo><mn>0,</mn><msubsup><mi>A</mi><mo>⊥</mo><mi>i</mi></msubsup><mo stretchy="false">(</mo><msub><mtext>x</mtext><mo>⊥</mo></msub><mo stretchy="false">)</mo><mn>,0</mn><mo stretchy="false">)</mo></mrow></math> , and apparently breaks the translational invariance in the transverse plane. Whereas the momentum along the magnetic field is conserved, the transverse momentum is not a conserved quantity and is moreover gauge dependent. For example, one of the transverse components is not conserved in the Landau gauge, while neither of them are conserved in the symmetric gauge. Nevertheless, it is clear that a constant magnetic field has a translational invariance in the transverse plane so that there must be a conserved momentum which works as the generator of the translation. This momentum is called the pseudomomentum, and its generalization to the two-body system is defined as (see Refs. [ 27 27 , 204 204 ] and references therein) <math display="block" id="M103"><mrow><mtext>K</mtext><mo>=</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1,2</mn></mrow></munder></mstyle><mrow><mo>(</mo><mrow><msup><mi>p</mi><mrow><mo stretchy="false">(</mo><mi>i</mi><mo stretchy="false">)</mo></mrow></msup><mo>+</mo><msup><mi>q</mi><mrow><mo stretchy="false">(</mo><mi>i</mi><mo stretchy="false">)</mo></mrow></msup><msup><mi>A</mi><mrow><mo stretchy="false">(</mo><mi>i</mi><mo stretchy="false">)</mo></mrow></msup></mrow><mo>)</mo></mrow><mtext> </mtext><mn>,</mn></mrow></math> where the labels denote the quark and antiquark as before. This momentum is related to the center-of-mass kinetic momentum as 〈 P kinetic 〉= K - qB ×〈 r 〉 with q ≡ q (1) =- q (2) and the relative coordinate r = r 1 - r 2 . The same issue of the momentum nonconservation is seen in diagrammatic calculations in the presence of a magnetic field. In such calculations, the gauge- and translational-symmetry breaking parts are encoded in the so-called Schwinger phase appearing from the gauge-dependent wave function of fermions in a magnetic field [ 212 212 ]. Effects of the Schwinger phase on the Bethe-Salpeter equation for light mesons were discussed in Ref. [ 62 62 ]. The Schrödinger equation for the two-body system (34) can be solved numerically for a fixed value of the conserved momentum, K . The parameters for the Cornell potential are so tuned to reproduce the vacuum spectrum in Ref. [ 27 27 , 38 38 , 39 39 ], and are imported from the lattice simulation in a finite magnetic field in Ref. [ 33 33 ]. In those papers, the numerical codes were carefully tested to reproduce the analytic solution for a harmonic potential. Figure 12 Figure 12 shows the charmonium spectra from Ref. [ 27 27 ] with the potential in vacuum. At the vanishing momentum 〈 P kinetic 〉=0 GeV, one sees that the masses of the η c and longitudinal J/ψ decreases and increases, respectively, as expected from the mixing effect. On the other hand, the modification in the transverse J/ψ is weak without the mixing effect. The masses of the moving charmonia, which were obtained by subtracting the center-of-mass kinetic energy, tend to increase with an increasing in its momentum. The bottomoium spectra shown in Ref. [ 27 27 ] have qualitatively the same behavior as in charmonium spectra, with suppression in magnitude due to the smaller magnetic moments of the bottom quark and antiquark. This study was extended to the excited states of quarkonia, which was made possible by using the cylindrical Gaussian expansion method proposed in Ref. [ 38 38 ]. Since this method does not resort to the imaginary time evolution τ → ∞ to extract the spectra [ 27 27 ], it would be more suitable for investigating the excited states. The result of the systematic spectroscopy is shown in Fig. 13 Fig. 13 . In a weak magnetic field, the level repulsions between the nearest mixing partners give rise to the dominant modification of the mass spectra. As the magnetic field becomes stronger, there appear the level-crossing points where the wave functions are strongly mixed. The density plots in the small windows show the transition of the wave functions. One can see the squeezing of the wave function along the magnetic field in the strong-field regime as expected from the shrink of cyclotron orbits. Including the contributions of the excited states, the η c state is, as it is the lowest-lying state, pushed down additively by all of the mixing partners, and its mass monotonically decreases. This systematic computation was very recently extended to the spectra of bottomonia and open heavy flavors [ 39 39 ]. So far, we have not considered the modification of the Cornell potential by the magnetic field. Effects of a magnetic field on the potential was investigated in Ref. [ 32 32 ] by lattice QCD simulation. The authors extracted the strengths of the potential in two cases where the separation between the quark and antiquark is in parallel and transverse to the magnetic field. The left panel in Fig. 14 Fig. 14 shows the results from the lattice simulation. Compared to the potential in vacuum, the strength of the potential is enhanced and suppressed in the transverse ( XY ) and longitudinal ( Z ) directions, respectively. This behavior may be understood from the fact that the screening due to the vacuum polarization of light-quark pairs is coupled to only the temporal and parallel ( z ) components of the gluon field. The full azimuthal-angle and temperature dependences were recently studied by the lattice QCD simulation [ 34 34 ], and those effects were also discussed by the strong couping method [ 35 35 , 36 36 , 213 213 ]. The right panel in Fig. 14 Fig. 14 shows the contour plot of the potential which was obtained by interpolating the values in the transverse ( XY ) and longitudinal ( Z ) directions with respect to the azimuthal-angle dependence. In Ref. [ 33 33 ], the quarkonium spectra were investigated by including this anisotropic potential as well as the mixing effects. Figure 15 Figure 15 shows the bottomonium spectra with and without the inclusion of the anisotropic potential. The authors argued that effects of the anisotropic potential can be seen in the bottomonium spectra more strongly than in the charmonium spectra, because the magnetic moment of the bottom quark is suppressed by the larger value of the quark mass and the smaller electric charge. Without the anisotropic potential, the spectra barely change due to the suppression of the mixing effects while the effects of the anisotropic potential act to increase the bottomonium masses in all three of the channels. The effect of the anisotropic potential is also much stronger than the momentum dependence of the spectra. On the basis of the potential model, several phenomenological consequences have also been discussed. Prior to moving to the next section, we briefly summarize these studies. • In the presence of the mixing effect, the wave function of the mass eigenstate is the mixture of the two components so that the mass eigenstate shares the decay modes with both of the mixing partners. This means, for example, that the physical η c and η b are allowed to decay into the dilepton while the dilepton decay channels of the physical longitudinal J/ψ and Υ are suppressed [ 27 27 ]. As mentioned in Ref. [ 27 27 ], a similar modification of the positronium decay mode is known. • The Lorentz forces acting on a moving quark and antiquark pair are oriented in the opposite direction. Therefore, the force will act to break up the bound state implying a possible dissociation mechanism in a strong magnetic field [ 23 23 , 27 27 ]. This effect may be also understood as an effect of an electric field appearing in the rest frame of a moving quarkonium, which is connected to the original magnetic field by the Lorentz boost. The electric field will stretch the separation between the quark and antiquark pair and induces a tunneling through the binding potential [ 23 23 ]. • Mass shifts of D mesons were also examined with use of the potential model [ 25 25 ]. In turn, these mass shifts induce the shift of the <math id="M106"><mrow><mi>D</mi><mover accent="true"><mi>D</mi><mo>¯</mo></mover></mrow></math> threshold above the quarkonium states. Based on the color evaporation model, modification of the charmonium production rate was discussed [ 25 25 ]. • A time dependence of the magnetic field was introduced in the potential model [ 37 37 ]. Based on the time-dependent model, the authors investigated the anisotropic production rate of high- p T charmonia. The feed-down effect from the excited states was also examined. 4.1.3 QCD sum rules in magnetic fields 3 Effects of the magnetic field on the meson spectra were also studied by using QCD sum rule (QCDSR) for open heavy flavor mesons [ 26 26 , 31 31 ] and quarkonia [ 29 29 , 30 30 ]. While we will focus here on the QCD sum rule for charmonia, extension to bottomonia is straightforward and extension to other systems will be discussed in the last of this section. In the QCD sum rules [ 214 214 - 217 217 214-217 ], the spectra of bound states are extracted from the time-ordered current correlator defined by <math display="block" id="M107"><mrow><msup><mi>Π</mi><mi>J</mi></msup><mo stretchy="false">(</mo><mi>q</mi><mo stretchy="false">)</mo><mo>=</mo><mi>i</mi><mstyle displaystyle="true"><mrow><mo>∫</mo></mrow></mstyle><mtext>​</mtext><mtext>​</mtext><msup><mi>d</mi><mn>4</mn></msup><mi>x</mi><mtext> </mtext><msup><mi>e</mi><mrow><mi>i</mi><mi>q</mi><mi>x</mi></mrow></msup><mo>〈</mo><mn>0</mn><mo>|</mo><mi>T</mi><mo stretchy="false">[</mo><mi>J</mi><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo><mi>J</mi><mo stretchy="false">(</mo><mn>0</mn><mo stretchy="false">)</mo><mo stretchy="false">]</mo><mo>|</mo><mn>0</mn><mo>〉</mo><mtext> </mtext><mn>,</mn></mrow></math> where the superscript J specifies a channel. The dispersion relation relates the hadronic spectral function in the time-like region ( q 2 =- Q 2 ≥0) to the correlator in the deep Euclidean region ( Q 2 → ∞) as <math display="block" id="M108"><mrow><msup><mi>Π</mi><mi>J</mi></msup><mo stretchy="false">(</mo><msup><mi>Q</mi><mn>2</mn></msup><mo stretchy="false">)</mo><mo>=</mo><mfrac><mn>1</mn><mi>π</mi></mfrac><mstyle displaystyle="true"><mrow><msubsup><mo>∫</mo><mn>0</mn><mi>∞</mi></msubsup></mrow></mstyle><mfrac><mrow><mtext> </mtext><mtext>I</mtext><mi>m</mi><mtext> </mtext><msup><mi>Π</mi><mi>J</mi></msup><mo stretchy="false">(</mo><mi>s</mi><mo stretchy="false">)</mo><mtext> </mtext></mrow><mrow><mi>s</mi><mo>+</mo><msup><mi>Q</mi><mn>2</mn></msup></mrow></mfrac><mtext> </mtext><mi>d</mi><mi>s</mi><mtext> </mtext><mo>+</mo><mtext> </mtext><mo stretchy="false">(</mo><mtext>subtractions</mtext><mo stretchy="false">)</mo><mn>.</mn></mrow></math> One may evaluate the correlator on the left-hand side on the basis of the Operator Product Expansion (OPE) in the deep Euclidean region. Here, we shall focus on the charmonia created by the currents: <math display="block" id="M109"><mrow><msub><mi>J</mi><mi>P</mi></msub><mo>=</mo><mi>i</mi><mover accent="true"><mi>c</mi><mo>¯</mo></mover><msup><mi>γ</mi><mn>5</mn></msup><mi>c</mi><mn>,</mn><mtext> </mtext><msubsup><mi>J</mi><mi>V</mi><mi>μ</mi></msubsup><mo>=</mo><mover accent="true"><mi>c</mi><mo>¯</mo></mover><msup><mi>γ</mi><mi>μ</mi></msup><mi>c</mi><mn>.</mn></mrow></math> Up to the dimension-4 operators, the OPE can be written down as <math display="block" id="M110"><mrow><mi>Π</mi><mo>=</mo><msub><mi>C</mi><mn>0</mn></msub><mo>+</mo><msub><mi>C</mi><mrow><msup><mi>G</mi><mn>2</mn></msup></mrow></msub><mo>〈</mo><msubsup><mi>G</mi><mrow><mi>μ</mi><mi>ν</mi></mrow><mi>a</mi></msubsup><msup><mi>G</mi><mrow><mi>a</mi><mtext> </mtext><mi>μ</mi><mi>ν</mi></mrow></msup><mo>〉</mo><mo>+</mo><msubsup><mi>C</mi><mrow><msup><mi>B</mi><mn>2</mn></msup></mrow><mrow><mi>μ</mi><mi>ν</mi><mi>α</mi><mi>β</mi></mrow></msubsup><mo>〈</mo><msubsup><mi>F</mi><mrow><mi>μ</mi><mi>ν</mi></mrow><mrow><mtext>e</mtext><mi>x</mi><mi>t</mi></mrow></msubsup><msubsup><mi>F</mi><mrow><mi>α</mi><mi>β</mi></mrow><mrow><mtext>e</mtext><mi>x</mi><mi>t</mi></mrow></msubsup><mo>〉</mo><mtext> </mtext><mn>.</mn></mrow></math> The Wilson coefficients, C 0 , C G 2 , and C B 2 , are responsible for the short-distance processes, so that one may perturbatively calculate the diagrams for the hard quark loops with insertions of the soft external gluon/magnetic fields (see Refs. [ 29 29 , 30 30 ] for explicit expressions). This can be done by using standard methods in the Fock-Schwinger gauge, which is also as known as the fixed-point gauge [ 216 216 , 218 218 ]. On the other side of the dispersion relation (43), we assume an Ansatz for the functional form of the spectral function ρ ( s )= π -1 Im∏ J ( s ) to extract the lowest-lying pole. In conventional analyses, the lowest-lying pole is simply separated from the continuum as <math display="block" id="M111"><mrow><msub><mi>ρ</mi><mrow><mtext>v</mtext><mi>a</mi><mi>c</mi></mrow></msub><mo stretchy="false">(</mo><mi>s</mi><mo stretchy="false">)</mo><mo>=</mo><msub><mi>f</mi><mn>0</mn></msub><mi>δ</mi><mo stretchy="false">(</mo><mi>s</mi><mo>−</mo><msup><mi>m</mi><mn>2</mn></msup><mo stretchy="false">)</mo><mo>+</mo><mfrac><mn>1</mn><mi>π</mi></mfrac><mtext>I</mtext><mi>m</mi><mtext> </mtext><msubsup><mi>Π</mi><mrow><mtext>pert</mtext></mrow><mi>J</mi></msubsup><mo stretchy="false">(</mo><mi>s</mi><mo stretchy="false">)</mo><mi>θ</mi><mo stretchy="false">(</mo><mi>s</mi><mo>−</mo><msub><mi>s</mi><mrow><mtext>th</mtext></mrow></msub><mo stretchy="false">)</mo><mn>.</mn></mrow></math> The parameters s th and f 0 are the threshold of the continuum and the coupling strength between the current and the lowest-lying pole, respectively. We, however, need to carefully reexamine the Ansatz (46) in the presence of a magnetic field, because effects of a magnetic field should be consistently taken into account on both sides of the dispersion relation (43) [ 26 26 , 29 29 - 31 31 29-31 ]. The effect in the hadronic degrees of freedom is nothing but the mixing effect discussed in the previous sections. In the language of the spectral function, the mixing effect can be seen as the emergence of a new pole for the mixing partner, which is the longitudinal J/ψ in the pseudoscalar channel and the η c in the vector channel. A schematic picture in Fig. 16 Fig. 16 shows the emergence of the longitudinal J/ψ in the pseudoscalar channel, which induces the repulsive shifts of the poles and the reduction of the residue of η c . It is important to include these effects in order to properly maintain the normalization of the spectral function. In the pseudoscalar channel, the second-order term in eB can be obtained by first converting the current to the psedoscalar meson with the strength f 0 and then using the results for the mixing effect shown in Eqs. (30) and (31). Therefore, the correlator can be written in the hadronic degrees of freedom as [ 29 29 ] <math display="block" id="M112"><mrow><msubsup><mi>Π</mi><mrow><mtext>2</mtext><mi>n</mi><mi>d</mi></mrow><mi>P</mi></msubsup><mo stretchy="false">(</mo><msup><mi>q</mi><mn>2</mn></msup><mo stretchy="false">)</mo><mo>=</mo><msub><mi>f</mi><mn>0</mn></msub><mrow><mo>[</mo><mrow><mfrac><mrow><mo>|</mo><msub><mrow><mo>〈</mo><mi>P</mi><mo>|</mo><msub><mi>η</mi><mi>c</mi></msub><mo>〉</mo></mrow><mi>B</mi></msub><msup><mo>|</mo><mn>2</mn></msup></mrow><mrow><msup><mi>q</mi><mn>2</mn></msup><mo>−</mo><msubsup><mi>m</mi><mrow><msub><mi>η</mi><mi>c</mi></msub></mrow><mn>2</mn></msubsup></mrow></mfrac><mo>+</mo><mfrac><mrow><mo>|</mo><msub><mrow><mo>〈</mo><mi>P</mi><mo>|</mo><mi>J</mi><mo>/</mo><mi>ψ</mi><mo>〉</mo></mrow><mi>B</mi></msub><msup><mo>|</mo><mn>2</mn></msup></mrow><mrow><msup><mi>q</mi><mn>2</mn></msup><mo>−</mo><msubsup><mi>m</mi><mrow><mi>J</mi><mo>/</mo><mi>ψ</mi></mrow><mn>2</mn></msubsup></mrow></mfrac><mtext> </mtext></mrow><mo>]</mo></mrow><mn>.</mn></mrow></math> Inserting Eqs. (30) and (31) into Eq. (47) and expanding up to the second order in eB , we find <math display="block" id="M113"><mrow><mtext> </mtext><msubsup><mi>Π</mi><mrow><mtext>2</mtext><mi>n</mi><mi>d</mi></mrow><mi>P</mi></msubsup><mo stretchy="false">(</mo><msup><mi>q</mi><mn>2</mn></msup><mo stretchy="false">)</mo><mo>=</mo><msub><mi>f</mi><mn>0</mn></msub><mfrac><mrow><msup><mi>γ</mi><mn>2</mn></msup></mrow><mrow><msubsup><mi>M</mi><mo>−</mo><mn>4</mn></msubsup></mrow></mfrac><mrow><mo>[</mo><mrow><mfrac><mn>1</mn><mrow><msup><mi>q</mi><mn>2</mn></msup><mo>−</mo><msubsup><mi>m</mi><mtext>V</mtext><mn>2</mn></msubsup></mrow></mfrac><mo>−</mo><mfrac><mn>1</mn><mrow><msup><mi>q</mi><mn>2</mn></msup><mo>−</mo><msubsup><mi>m</mi><mtext>P</mtext><mn>2</mn></msubsup></mrow></mfrac><mo>−</mo><mfrac><mrow><msubsup><mi>M</mi><mo>−</mo><mn>2</mn></msubsup></mrow><mrow><msup><mrow><mo stretchy="false">(</mo><msup><mi>q</mi><mn>2</mn></msup><mo>−</mo><msubsup><mi>m</mi><mtext>P</mtext><mn>2</mn></msubsup><mo stretchy="false">)</mo></mrow><mn>2</mn></msup></mrow></mfrac></mrow><mo>]</mo></mrow><mn>.</mn></mrow></math> Interpretation of the terms in Eq. (48) are as follows. The first term corresponds to production of the on-shell longitudinal J/ψ from the pseudoscalar current via an off-shell η c . The second term with a negative sign is needed to preserve the normalization because the coupling of η c to the current must be reduced to balance the emergence of the coupling to J/ψ . This is confirmed in Eq. (47), where these two terms come from overlaps between the properly normalized unperturbated and perturbated states obtained as <math id="M114"><mrow><mo>|</mo><mo stretchy="false">(</mo><mi>P</mi><mo>|</mo><mi>η</mi><mi>c</mi><msub><mo stretchy="false">)</mo><mi>B</mi></msub><msup><mo>|</mo><mn>2</mn></msup><mo>∼</mo><mn>1</mn><mo>−</mo><msup><mrow><mo stretchy="false">(</mo><mi>γ</mi><mo>/</mo><msubsup><mi>M</mi><mo>−</mo><mn>2</mn></msubsup><mo stretchy="false">)</mo></mrow><mn>2</mn></msup></mrow></math> and <math id="M115"><mrow><mo>|</mo><mo stretchy="false">(</mo><mi>P</mi><mo>|</mo><mi>J</mi><mo>/</mo><mi>ψ</mi><msub><mo stretchy="false">)</mo><mi>B</mi></msub><msup><mo>|</mo><mn>2</mn></msup><mo>∼</mo><msup><mrow><mo stretchy="false">(</mo><mi>γ</mi><mo>/</mo><msubsup><mi>M</mi><mo>−</mo><mn>2</mn></msubsup><mo stretchy="false">)</mo></mrow><mn>2</mn></msup></mrow></math> . The third term has a double pole on the η c mass with a factor <math id="M116"><mrow><msubsup><mi>M</mi><mo>−</mo><mn>2</mn></msubsup></mrow></math> corresponding to a virtual transition to the longitudinal J/ψ state between on-shell η c states, which is nothing but the origin of the mass shift due to the mixing effect. In Eq. (47), this term comes from an expansion with respect to the mass correction shown in Eq. (30). Clearly, if one includes this mixing term in the phenomenological spectral function, its effect is subtracted out from the total mass shift obtained from the QCDSR and thus can be separated from the residual effects of B -fields, not described in the hadronic level. A similar calculation can be applied to the longitudinal vector channel. On the phenomenological side, those terms should be included as well as the vacuum contributions (46) in both the pseudoscalar and longitudinal vector channels. To make the dispersion integral (43) dominated by the pole contribution and to improve the convergence of the OPE, we use the Borel transform defined by <math display="block" id="M117"><mrow><mi mathvariant="script">B</mi><mo stretchy="false">[</mo><mi>f</mi><mo stretchy="false">(</mo><msup><mi>Q</mi><mn>2</mn></msup><mo stretchy="false">)</mo><mo stretchy="false">]</mo><mo>≡</mo><mtext>​</mtext><mtext>​</mtext><munder><mstyle mathsize="140%" displaystyle="true"><mrow><mi>lim</mi></mrow></mstyle><mrow><mtable><mtr><mtd><mrow><msup><mi>Q</mi><mn>2</mn></msup><mn>,</mn><mi>n</mi><mo>→</mo><mi>∞</mi></mrow></mtd></mtr><mtr><mtd><mrow><msup><mi>Q</mi><mn>2</mn></msup><mo>/</mo><mi>n</mi><mo>=</mo><msup><mi>M</mi><mn>2</mn></msup></mrow></mtd></mtr></mtable></mrow></munder><mfrac><mrow><msup><mrow><mo stretchy="false">(</mo><msup><mi>Q</mi><mn>2</mn></msup><mo stretchy="false">)</mo></mrow><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msup></mrow><mrow><mi>n</mi><mo>!</mo></mrow></mfrac><msup><mrow><mrow><mo>(</mo><mrow><mo>−</mo><mfrac><mrow><mi>d</mi><mtext> </mtext></mrow><mrow><mi>d</mi><msup><mi>Q</mi><mn>2</mn></msup></mrow></mfrac></mrow><mo>)</mo></mrow></mrow><mi>n</mi></msup><mi>f</mi><mo stretchy="false">(</mo><msup><mi>Q</mi><mn>2</mn></msup><mo stretchy="false">)</mo><mn>.</mn></mrow></math> By applying the Borel transform to both sides of Eq. (43), the contribution of the delta function for the lowest-lying pole is picked up to result in a mass formula: <math display="block" id="M118"><mrow><msup><mi>m</mi><mn>2</mn></msup><mo>=</mo><mo>−</mo><mfrac><mrow><mo>∂</mo><mtext> </mtext></mrow><mrow><mo>∂</mo><mo stretchy="false">(</mo><mn>1</mn><mo>/</mo><msup><mi>M</mi><mn>2</mn></msup><mo stretchy="false">)</mo></mrow></mfrac><mi>ln</mi><mo stretchy="false">[</mo><msub><mi mathvariant="script">M</mi><mrow><mtext>OPE</mtext></mrow></msub><mo>−</mo><msub><mi mathvariant="script">M</mi><mrow><mtext>cont</mtext></mrow></msub><mo>−</mo><msub><mi mathvariant="script">M</mi><mrow><mtext>mix</mtext></mrow></msub><mo stretchy="false">]</mo><mtext> </mtext><mn>,</mn></mrow></math> where <math id="M119"><mrow><msub><mi mathvariant="script">M</mi><mrow><mtext>OPE</mtext></mrow></msub></mrow></math> , <math id="M120"><mrow><msub><mi mathvariant="script">M</mi><mrow><mtext>cont</mtext></mrow></msub></mrow></math> , and <math id="M121"><mrow><msub><mi mathvariant="script">M</mi><mrow><mtext>mix</mtext></mrow></msub></mrow></math> are the Borel transform of the OPE (45), the perturbative continuum, and those induced by the magnetic field (48), respectively. More details and explicit forms are given in Refs. [ 29 29 , 30 30 ]. Figure 17 Figure 17 shows the mass spectra from the QCDSR together with those from the hadronic effective model shown in Fig. 11. Fig. 11. First, the results from these two methods are consistent with each other, indicating the dominance of the mixing effect in the shown regime. They are also consistent with the results from the potential model discussed in the previous section [ 27 27 , 38 38 , 39 39 , 33 33 ]. Also, the modification in the transverse components of J/ψ (green line) is much smaller than the others due to the absence of the mixing partner. On top of these observations, one finds deviations between the results from the two methods. Although they are small deviations in the current regime of the magnetic field, they suggest the precursor of other effects, which are not explained by the mixing effect. It would be natural to expect the emergence of such effects in stronger magnetic fields. Possible origin of these deviations were discussed in Ref. [ 30 30 ]. One caveat is that the modification of the gluon condensate in a magnetic field has not been included in the evaluation of the OPE (45). Since the lattice QCD simulation has provided the values of the gluon condensate in finite magnetic field and temperature [ 201 201 ], it will be interesting to investigate effects of these modifications by using the QCDSR. Also, applications of the QCDSR to light mesons and open heavy flavor mesons are an interesting topic because these mesons contain a light quark component, which reflects the modification of the quark condensates. In Refs. [ 26 26 , 31 31 ], the OPE was performed by taking into account not only the modification of <math id="M122"><mrow><mo>〈</mo><mover accent="true"><mi>q</mi><mo>¯</mo></mover><mi>q</mi><mo>〉</mo></mrow></math> but also the emergence of the tensor-type condensate <math id="M123"><mrow><mo>〈</mo><mover accent="true"><mi>q</mi><mo>¯</mo></mover><msup><mi>σ</mi><mrow><mi>μ</mi><mi>ν</mi></mrow></msup><mi>q</mi><mo>〉</mo></mrow></math> , which have been calculated by the lattice QCD simulations [ 48 48 - 50 50 48-50 , 200 200 ]. Furthermore, the mass spectra of the open heavy flavor mesons would serve as a signature of the magnetically induced QCD Kondo effect [ 219 219 ] which was proposed on the basis of a close analogy to the QCD Kondo effect in the dense QCD [ 220 220 , 221 221 ]. The readers are referred to the literature for the interesting analogy between the (1+1)-dimensional low-energy dynamics at high densisy [ 222 222 - 224 224 222-224 ] and in strong magnetic fields [ 43 43 , 53 53 , 225 225 , 226 226 ]. 4.2 Transport of heavy quarks in a hot medium and magnetic field 2 As mentioned in the beginning of this section, the heavy quarks are dominantly created by the initial hard scatterings among the partons from the colliding nuclei (see Fig. 4.2 Fig. 4.2 ). Therefore, the heavy quarks typically have the power law p T spectrum according to the perturbative QCD. This spectrum is however gradually modified by the multiple interactions with the thermal quarks and gluons in the QGP, and the spectrum approaches a thermal distribution, which is realized at a large enough time. However, the QGP has a finite lifetime, and to what extent the thermalization is achieved depends on the transport property of the QGP. Therefore, the heavy quarks will probe the transport property of the QGP. In the RHIC and LHC experiments, the nuclear modification factor, R AA , and the anisotropic spectrum, v 2 , of the open heavy flavors have been measured (see a review [ 202 202 ] and references therein). Both of them are thought to be sensitive to the heavy-quark thermalization and are closely related to each other [ 227 227 ]. However, a consistent theoretical modeling, which simultaneously reproduces R AA and v 2 , is still under investigation. Effects of the magnetic field could be new ingredients to resolve the discrepancy between the theoretical and experimental results. In this section, we summarize a recent study on the heavy quark dynamics in a hot medium and magnetic field on the basis of Langevin equation and the perturbative computation of the heavy-quark diffusion constant. The dynamics of the low-energy heavy particles is modeled as a Brownian motion caused by random kicks by the thermal particles which is described by the Langevin equations [ 227 227 ]: <math display="block" id="M124"><mrow><mfrac><mrow><mi>d</mi><msub><mi>p</mi><mi>z</mi></msub></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mo>−</mo><msub><mi>η</mi><mo>∥</mo></msub><msub><mi>p</mi><mi>z</mi></msub><mo>+</mo><msub><mi>ξ</mi><mi>z</mi></msub><mtext> </mtext><mn>,</mn><mtext> </mtext><mtext> </mtext><mfrac><mrow><mi>d</mi><msub><mi>p</mi><mo>⊥</mo></msub></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mo>−</mo><msub><mi>η</mi><mo>⊥</mo></msub><mtext> </mtext><msub><mi>p</mi><mo>⊥</mo></msub><mo>+</mo><msub><mi>ξ</mi><mo>⊥</mo></msub><mtext> </mtext><mn>.</mn></mrow></math> Since the external magnetic field provides a preferred spatial direction, we have a set of two equations for the heavy quark motions, one parallel and one perpendicular to the magnetic field that is oriented in the z -direction. The random forces are assumed to be white noises; <math display="block" id="M125"><mrow><mo>〈</mo><msub><mi>ξ</mi><mi>z</mi></msub><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><msub><mi>ξ</mi><mi>z</mi></msub><mo stretchy="false">(</mo><msup><mi>t</mi><mo>′</mo></msup><mo stretchy="false">)</mo><mo>〉</mo><mo>=</mo><msub><mi>κ</mi><mo>∥</mo></msub><mi>δ</mi><mo stretchy="false">(</mo><mi>t</mi><mo>−</mo><msup><mi>t</mi><mo>′</mo></msup><mo stretchy="false">)</mo><mtext> </mtext><mn>,</mn><mtext> </mtext><mtext> </mtext><mo>〈</mo><msubsup><mi>ξ</mi><mo>⊥</mo><mi>i</mi></msubsup><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><msubsup><mi>ξ</mi><mo>⊥</mo><mi>j</mi></msubsup><mo stretchy="false">(</mo><msup><mi>t</mi><mo>′</mo></msup><mo stretchy="false">)</mo><mo>〉</mo><mo>=</mo><msub><mi>κ</mi><mo>⊥</mo></msub><msup><mi>δ</mi><mrow><mi>i</mi><mi>j</mi></mrow></msup><mi>δ</mi><mo stretchy="false">(</mo><mi>t</mi><mo>−</mo><msup><mi>t</mi><mo>′</mo></msup><mo stretchy="false">)</mo><mtext> </mtext><mo stretchy="false">(</mo><mi>i</mi><mn>,</mn><mi>j</mi><mo>=</mo><mi>x</mi><mn>,</mn><mi>y</mi><mo stretchy="false">)</mo></mrow></math> and these coefficients, κ ∥ and κ ⊥ , are related to the drag coefficients, ≤ ∥ and ≤ ⊥ , through the fluctuation-dissipation theorem as <math display="block" id="M126"><mrow><msub><mi>η</mi><mo>∥</mo></msub><mo>=</mo><mn>2</mn><msub><mi>M</mi><mi>Q</mi></msub><mi>T</mi><msub><mi>κ</mi><mo>∥</mo></msub><mtext> </mtext><mn>,</mn><mtext> </mtext><mtext> </mtext><msub><mi>η</mi><mo>⊥</mo></msub><mo>=</mo><mn>2</mn><msub><mi>M</mi><mi>Q</mi></msub><mi>T</mi><msub><mi>κ</mi><mo>⊥</mo></msub><mtext> </mtext><mn>.</mn></mrow></math> Here, we would like to mention a recent numerical study [ 42 42 ] where the authors investigated the Lorentz force exerting on the heavy quark with the inclusion of the Lorentz force in the Langevin equation (51). In addition to effects of the Lorentz force exerting on the heavy quark as a direct effect of the magnetic field, the diffusion and drag coefficients will be changed because the light quarks in the hot medium are strongly affected by the magnetic field [ 40 40 , 41 41 ]. At the leading order in g s , the anisotropic momentum diffusion coefficients, κ ∥ and κ ⊥ , can be defined by <math display="block" id="M127"><mrow><msub><mi>κ</mi><mo>∥</mo></msub><mo>=</mo><mstyle displaystyle="true"><mrow><mo>∫</mo></mrow></mstyle><msup><mi>d</mi><mn>3</mn></msup><mi>q</mi><mtext> </mtext><mfrac><mrow><mi>d</mi><mi>Γ</mi><mo stretchy="false">(</mo><mi>q</mi><mo stretchy="false">)</mo></mrow><mrow><msup><mi>d</mi><mn>3</mn></msup><mi>q</mi></mrow></mfrac><mtext> </mtext><msubsup><mi>q</mi><mi>z</mi><mn>2</mn></msubsup><mtext> </mtext><mn>,</mn><mtext> </mtext><mtext> </mtext><msub><mi>κ</mi><mo>⊥</mo></msub><mo>=</mo><mn>12</mn><mstyle displaystyle="true"><mrow><mo>∫</mo></mrow></mstyle><msup><mi>d</mi><mn>3</mn></msup><mi>q</mi><mtext> </mtext><mfrac><mrow><mi>d</mi><mi>Γ</mi><mo stretchy="false">(</mo><mi>q</mi><mo stretchy="false">)</mo></mrow><mrow><msup><mi>d</mi><mn>3</mn></msup><mi>q</mi></mrow></mfrac><msubsup><mi>q</mi><mo>⊥</mo><mn>2</mn></msubsup><mn>.</mn><mtext> </mtext></mrow></math> where q is the amount of the momentum transfer from the thermal particles to the heavy quark, and the static limit ( q 0 → 0) is assumed in the above definitions. The momentum transfer rate <math id="M128"><mrow><mfrac><mrow><mi>d</mi><mi>Γ</mi><mo stretchy="false">(</mo><mi>q</mi><mo stretchy="false">)</mo></mrow><mrow><msup><mi>d</mi><mn>3</mn></msup><mi>q</mi></mrow></mfrac></mrow></math> is computed from the gluon-exchange diagrams shown in Fig. 18. Fig. 18. Effects of a magnetic field appear in two places (highlighted by red). While the gluons are not directly coupled to the magnetic field, (i) the Debye screening mass and (ii) the dispersion relation of the thermal-quark scatterers will be changed as discussed below. The Debye mass is necessary for cutting off the infrared divergence in the forward scattering. In a strong magnetic field, the fermion wave function is strongly squeezed along the magnetic field, corresponding to the small radius of the cyclotron orbit. Indeed, from the Landau level quantization and the Zeeman effect, the dispersion relation of fermions in the lowest Landau level (LLL) becomes the (1+1) dimensional, i.e., <math id="M129"><mrow><msup><mi>ϵ</mi><mn>2</mn></msup><mo>=</mo><msup><mi>m</mi><mn>2</mn></msup><mo>+</mo><msubsup><mi>p</mi><mi>z</mi><mn>2</mn></msubsup></mrow></math> for a massive fermion and ϵ =± p z for massless fermions with right- and left-handed chiralities. To be specific, we focus on the strong-field regime such that the transition from the LLL to the hLL states are suppressed according to a hierarchy T 2 ≪ eB . (1) First, the Debye screening mass is obtained from the gluon self-energy computed by the standard method in finite temperature field theory in the (1+1) dimension. The gluon self-energy is completely factorized into the transverse and longitudinal parts as <math display="block" id="M130"><mrow><msup><mi>Π</mi><mrow><mi>μ</mi><mi>ν</mi></mrow></msup><mo stretchy="false">(</mo><mi>q</mi><mo stretchy="false">)</mo><mo>=</mo><msub><mi>q</mi><mi>f</mi></msub><mfrac><mrow><mi>e</mi><mi>B</mi></mrow><mrow><mn>2</mn><mi>π</mi></mrow></mfrac><msup><mi>e</mi><mrow><mo>−</mo><mfrac><mrow><mo>|</mo><msub><mtext>q</mtext><mo>⊥</mo></msub><msup><mo>|</mo><mn>2</mn></msup></mrow><mrow><mn>2</mn><mo>|</mo><msub><mi>q</mi><mi>f</mi></msub><mi>e</mi><mi>B</mi><mo>|</mo></mrow></mfrac></mrow></msup><msubsup><mi>Π</mi><mrow><mn>1</mn><mo>+</mo><mn>1</mn></mrow><mrow><mi>μ</mi><mi>ν</mi></mrow></msubsup><mo stretchy="false">(</mo><msub><mi>q</mi><mo>∥</mo></msub><mo stretchy="false">)</mo><mtext> </mtext><mn>,</mn></mrow></math> where the factor of <math id="M131"><mrow><msub><mi>q</mi><mi>f</mi></msub><mfrac><mrow><mi>e</mi><mi>B</mi></mrow><mrow><mn>2</mn><mi>π</mi></mrow></mfrac></mrow></math> comes from the degeneracy factor in the transverse phase space, and the Gaussian is the wave function of the LLL state. The longitudinal part is nothing but the (1+1)-dimensional polarization tensor, which is known as the Schwinger model: <math display="block" id="M132"><mrow><msubsup><mi>Π</mi><mrow><mn>1</mn><mo>+</mo><mn>1</mn></mrow><mrow><mi>μ</mi><mi>ν</mi></mrow></msubsup><mo stretchy="false">(</mo><msub><mi>q</mi><mo>∥</mo></msub><mo stretchy="false">)</mo><mo>=</mo><mtext>t</mtext><mi>r</mi><mo stretchy="false">[</mo><msup><mi>t</mi><mi>a</mi></msup><msup><mi>t</mi><mi>a</mi></msup><mo stretchy="false">]</mo><mfrac><mrow><msubsup><mi>g</mi><mi>s</mi><mn>2</mn></msubsup></mrow><mi>π</mi></mfrac><mi>f</mi><mo stretchy="false">(</mo><msup><mi>q</mi><mn>0</mn></msup><mn>,</mn><msub><mi>q</mi><mi>z</mi></msub><mo stretchy="false">)</mo><mo stretchy="false">(</mo><msubsup><mi>q</mi><mo>∥</mo><mn>2</mn></msubsup><msubsup><mi>g</mi><mo>∥</mo><mrow><mi>μ</mi><mi>ν</mi></mrow></msubsup><mo>−</mo><msubsup><mi>q</mi><mo>∥</mo><mi>μ</mi></msubsup><msubsup><mi>q</mi><mo>∥</mo><mi>ν</mi></msubsup><mo stretchy="false">)</mo><mtext> </mtext><mn>.</mn></mrow></math> Here, the longitudinal momentum is defined by <math id="M133"><mrow><msubsup><mi>q</mi><mo>∥</mo><mi>μ</mi></msubsup><mo>=</mo><mo stretchy="false">(</mo><msup><mi>q</mi><mn>0</mn></msup><mn>,0,0,</mn><msub><mi>q</mi><mi>z</mi></msub><mo stretchy="false">)</mo></mrow></math> . There is only one gauge-invariant tensor structure in the (1+1) dimension, so the temperature correction is, if any, contained in a function f ( q 0 , q z ). However, an important observation is that there is no temperature or density correction in the massless case [ 40 40 , 228 228 , 229 229 ], and f ( q 0 , q z ) = 1. Accordingly, in the massive case, a temperature correction is proportional to the quark mass and is suppressed by m q / T (see Ref. [ 40 40 ] for explicit expressions). Therefore, in both cases, the Debye mass in strong magnetic fields ( eB ≫ T 2 ) is given by the vacuum part of the gluon self-energy [ 230 230 - 232 232 230-232 ]. Namely, we find <math id="M134"><mrow><msubsup><mi>m</mi><mi>D</mi><mn>2</mn></msubsup><mo>∼</mo><msub><mi>q</mi><mi>f</mi></msub><mfrac><mrow><mi>e</mi><mi>B</mi></mrow><mrow><mn>2</mn><mi>π</mi></mrow></mfrac><mo>⋅</mo><mfrac><mrow><msubsup><mi>g</mi><mi>s</mi><mn>2</mn></msubsup></mrow><mi>π</mi></mfrac></mrow></math> which is much larger than the usual thermal mass squared, ∼( g s T ) 2 . (2) The change in the dispersion relation of fermion scatterers is important for the kinematics of the Coulomb scattering process [ 40 40 , 233 233 ]. We shall consider a massless fermion scatterer in the lowest Landau level (see Fig. 19 Fig. 19 ). For the massless quarks to be on-shell both in the initial and final states, the energy transfer has to be <math id="M135"><mrow><msup><mi>q</mi><mn>0</mn></msup><mo>=</mo><mo>±</mo><mo stretchy="false">(</mo><msubsup><mi>k</mi><mi>z</mi><mo>'</mo></msubsup><mo>−</mo><msub><mi>k</mi><mi>z</mi></msub><mo stretchy="false">)</mo></mrow></math> and thus q 0 =± q z . Note that the chirality does not change at any perturbative vertex, so that the signs appear only as the overall ones. When taking the static limit q 0 → 0 in the definition of the diffusion constant, one can immediately conclude that the longitudinal-momentum transfer is kinematically prohibited in the massless case. On the other hand, this constraint does not apply to the transverse momentum transfer, so the transverse momentum transfer is allowed. This finite contribution to the diffusion constant can be obtained either by the direct computation of the matrix elements in Fig. 19 Fig. 19 or by using the cutting rule. At the leading order in α s , the transverse momentum diffusion constant in the massless limit is obtained as [ 40 40 ] <math display="block" id="M136"><mrow><msubsup><mi>κ</mi><mo>⊥</mo><mrow><mtext>L</mtext><mi>O</mi></mrow></msubsup><mo>∼</mo><msubsup><mi>α</mi><mi>s</mi><mn>2</mn></msubsup><mi>T</mi><mo stretchy="false">(</mo><mfrac><mrow><mi>e</mi><mi>B</mi></mrow><mrow><mn>2</mn><mi>π</mi></mrow></mfrac><mo stretchy="false">)</mo><mi>ln</mi><msubsup><mi>α</mi><mi>s</mi><mrow><mo>−</mo><mn>1</mn></mrow></msubsup><mtext> </mtext><mn>.</mn></mrow></math> Therefore, strong magnetic fields give rise to an anisotropy of the momentum diffusion constant. Nonvanishing contributions to the longitudinal momentum transfer come from either the finite quark-mass correction or the contribution of the gluon scatterers. In Ref. [ 40 40 ], the contribution of the mass correction was obtained as <math display="block" id="M137"><mrow><msubsup><mi>κ</mi><mo>∥</mo><mrow><mtext>L</mtext><mi>O</mi><mn>,</mn><mtext> </mtext><mi>m</mi><mi>a</mi><mi>s</mi><mi>s</mi><mi>i</mi><mi>v</mi><mi>e</mi></mrow></msubsup><mo>∼</mo><msub><mi>α</mi><mi>s</mi></msub><msubsup><mi>m</mi><mi>q</mi><mn>2</mn></msubsup><msqrt><mrow><msub><mi>α</mi><mi>s</mi></msub><mi>e</mi><mi>B</mi></mrow></msqrt><mtext> </mtext><mn>,</mn></mrow></math> while the gluon contribution is obtained by substituting the Debye mass <math id="M138"><mrow><msubsup><mi>m</mi><mi>D</mi><mn>2</mn></msubsup><mo>∼</mo><msub><mi>α</mi><mi>s</mi></msub><mi>e</mi><mi>B</mi></mrow></math> for the conventional one ( g s T ) 2 in Refs. [ 227 227 , 234 234 ] as <math display="block" id="M139"><mrow><msubsup><mi>κ</mi><mo>∥</mo><mrow><mtext>L</mtext><mi>O</mi><mn>,</mn><mtext> </mtext><mi>g</mi><mi>l</mi><mi>u</mi><mi>o</mi><mi>n</mi></mrow></msubsup><mo>∼</mo><msubsup><mi>α</mi><mi>s</mi><mn>2</mn></msubsup><msup><mi>T</mi><mn>3</mn></msup><mi>ln</mi><msubsup><mi>α</mi><mi>s</mi><mrow><mo>−</mo><mn>1</mn></mrow></msubsup><mtext> </mtext><mn>.</mn></mrow></math> We have assumed a hierarchy α s eB ≪ T 2 ≪ eB . Complete expressions, including the prefactors, are found in Ref. [ 40 40 ]. It is instructive to compare the mass correction (58) with the gluon contribution (59). The ratio is written as <math display="block" id="M140"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mfrac><mrow><msubsup><mi>κ</mi><mo>∥</mo><mrow><mtext>L</mtext><mi>O</mi><mn>,</mn><mtext> </mtext><mi>m</mi><mi>a</mi><mi>s</mi><mi>s</mi><mi>i</mi><mi>v</mi><mi>e</mi></mrow></msubsup></mrow><mrow><msubsup><mi>κ</mi><mo>∥</mo><mrow><mtext>L</mtext><mi>O</mi><mn>,</mn><mtext> </mtext><mi>g</mi><mi>l</mi><mi>u</mi><mi>o</mi><mi>n</mi></mrow></msubsup></mrow></mfrac><mo>∼</mo><mfrac><mrow><msub><mi>α</mi><mi>s</mi></msub><msup><mrow><mo stretchy="false">(</mo><msub><mi>α</mi><mi>s</mi></msub><mi>e</mi><mi>B</mi><mo stretchy="false">)</mo></mrow><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup><msubsup><mi>m</mi><mi>q</mi><mn>2</mn></msubsup></mrow><mrow><msubsup><mi>α</mi><mi>s</mi><mn>2</mn></msubsup><msup><mi>T</mi><mn>3</mn></msup></mrow></mfrac></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow><mo>=</mo><mo stretchy="false">(</mo><mfrac><mrow><msubsup><mi>m</mi><mi>q</mi><mn>2</mn></msubsup></mrow><mrow><msub><mi>α</mi><mi>s</mi></msub><mi>e</mi><mi>B</mi></mrow></mfrac><mo stretchy="false">)</mo><mo stretchy="false">(</mo><mfrac><mrow><msub><mi>α</mi><mi>s</mi></msub><mi>e</mi><mi>B</mi></mrow><mrow><msup><mi>T</mi><mn>2</mn></msup></mrow></mfrac><msup><mo stretchy="false">)</mo><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup><mo stretchy="false">(</mo><mfrac><mrow><mi>e</mi><mi>B</mi></mrow><mrow><msup><mi>T</mi><mn>2</mn></msup></mrow></mfrac><mo stretchy="false">)</mo><mtext> </mtext><mn>.</mn></mrow></mtd></mtr></mtable></mrow></math> While the first two factors are small in our working regime, the last factor can be large. Therefore, the massive contribution <math id="M141"><mrow><msubsup><mi>κ</mi><mo>∥</mo><mrow><mtext>L</mtext><mi>O</mi><mn>,</mn><mtext> </mtext><mi>m</mi><mi>a</mi><mi>s</mi><mi>s</mi><mi>i</mi><mi>v</mi><mi>e</mi></mrow></msubsup></mrow></math> could be in principle as comparably large as <math id="M142"><mrow><msubsup><mi>κ</mi><mo>∥</mo><mrow><mtext>L</mtext><mi>O</mi><mn>,</mn><mtext> </mtext><mi>g</mi><mi>l</mi><mi>u</mi><mi>o</mi><mi>n</mi></mrow></msubsup></mrow></math> , and this happens when <math id="M143"><mrow><mi>e</mi><mi>B</mi><mo>∼</mo><msub><mi>α</mi><mi>s</mi></msub><mo stretchy="false">(</mo><msup><mi>T</mi><mn>6</mn></msup><mo>/</mo><msubsup><mi>m</mi><mi>q</mi><mn>4</mn></msubsup><mo stretchy="false">)</mo></mrow></math> . However, to be consistent with our assumed regime, α s eB ≪ T 2 , we have a constraint of <math id="M144"><mrow><msub><mi>α</mi><mi>s</mi></msub><mo>≪</mo><msubsup><mi>m</mi><mi>q</mi><mn>2</mn></msubsup><mo>/</mo><msup><mi>T</mi><mn>2</mn></msup></mrow></math> , which is not quite likely true in the heavy ion collisions. Hence, the longitudinal momentum diffusion constant is dominated by the gluon contribution <math id="M145"><mrow><msubsup><mi>κ</mi><mo>∥</mo><mrow><mtext>L</mtext><mi>O</mi><mn>,</mn><mtext> </mtext><mi>g</mi><mi>l</mi><mi>u</mi><mi>o</mi><mi>n</mi></mrow></msubsup></mrow></math> . Now that we obtained the leading contributions in Eqs. (57) and (59), the anisotropy in the momentum diffusion constant and the drag force is estimated to be <math display="block" id="M146"><mrow><mfrac><mrow><msubsup><mi>κ</mi><mo>∥</mo><mrow><mtext>L</mtext><mi>O</mi><mn>,</mn><mtext> </mtext><mi>g</mi><mi>l</mi><mi>u</mi><mi>o</mi><mi>n</mi></mrow></msubsup></mrow><mrow><msubsup><mi>κ</mi><mo>⊥</mo><mrow><mtext>L</mtext><mi>O</mi></mrow></msubsup></mrow></mfrac><mo>=</mo><mfrac><mrow><msubsup><mi>η</mi><mo>∥</mo><mrow><mtext>L</mtext><mi>O</mi><mn>,</mn><mtext> </mtext><mi>g</mi><mi>l</mi><mi>u</mi><mi>o</mi><mi>n</mi></mrow></msubsup></mrow><mrow><msubsup><mi>η</mi><mo>⊥</mo><mrow><mtext>L</mtext><mi>O</mi></mrow></msubsup></mrow></mfrac><mo>∼</mo><mfrac><mrow><msup><mi>T</mi><mn>2</mn></msup></mrow><mrow><mi>e</mi><mi>B</mi></mrow></mfrac><mo>≪</mo><mn>1</mn><mtext> </mtext><mn>.</mn></mrow></math> This large anisotropy is induced by the enhancement of the density of states in the quark contribution <math id="M147"><mrow><msubsup><mi>κ</mi><mo>⊥</mo><mrow><mtext>L</mtext><mi>O</mi></mrow></msubsup></mrow></math> , which is proportional to eB . On the other hand, a magnetic field does not change the phase space volume of the thermal gluons in the gluon contribution <math id="M148"><mrow><msubsup><mi>κ</mi><mo>∥</mo><mrow><mtext>L</mtext><mi>O</mi><mn>,</mn><mtext> </mtext><mi>g</mi><mi>l</mi><mi>u</mi><mi>o</mi><mi>n</mi></mrow></msubsup></mrow></math> , but changes only the Debye screening mass. Based on the above findings, a toy model for the dynamical modeling was discussed in Ref. [ 40 40 ] by using the Fokker-Planck equation [ 227 227 ]. Since the friction exerted on the heavy quarks is stronger in the direction transverse to the magnetic field (see Fig. 20 Fig. 20 ), the heavy-quark momentum in the final state will acquire an anisotropy. In the small momentum region, the heavy quarks will be dragged by the flow more strongly in the transverse direction while they will feel a larger resistance in the transverse direction as the momentum increases. Therefore, there is a turnover in the anisotropy of the heavy quark spectrum (see the right panel in Fig. 20 Fig. 20 ). It will be interesting to implement the anisotropic diffusion and drag coefficient in the numerical simulations. As the magnetic field in the heavy-ion collisions is time dependent, it is also necessary to investigate the weak-field limit as the first step in considering the time dependence. These issues are left as open questions. Finally, the drag force in the presence of the anomalous currents was also investigated in Ref. [ 233 233 , 235 235 , 236 236 ]. 5 Summary 1 In summary, we have discussed the recent progresses in understanding the novel quantum phenomena in heavy-ion collisions. We first briefly reviewed some special features of the magnetic fields generated in heavy-ion collisions from both numerical and analytical approaches in Sec. 2. The magnetic fields in RHIC Au + Au and in LHC Pb + Pb collisions can be comparable to or even larger than, the QCD confinement scale, <math id="M149"><mrow><msubsup><mi>Λ</mi><mrow><mtext>Q</mtext><mi>C</mi><mi>D</mi></mrow><mn>2</mn></msubsup></mrow></math> , and thus are able to significantly influence the physics of strongly interacting matter. The interplay between the magnetic fields and the non-trivial topological configuration of QCD leads to a number of anomalous transports, for example, the chiral magnetic effect (CME), chiral separation effect (CSE), and chiral magnetic wave (CMW). These effects, once experimentally confirmed, provide "soft probes" of the topological sector of QCD in heavy-ion collisions. We discussed the basic mechanisms that lead to the occurrence of CME, CSE, and CMW in Sec. 3.1 and Sec. 3.3. The experimental search of CME and CMW in heavy-ion collisions has been performed for years, and the data already shows very encouraging features that are consistent with the expectations of CME and CMW. However, background effects, which are not related to CME and CMW, exist and make the interpretation of the experimental measurements ambiguous. More efforts from both experimental and theoretical sides are certainly needed before we can make any conclusive claims. We discussed the current experimental status of the search of CME and CMW in Sec. 3.2 and Sec. 3.4. In the latter half of this short review, we discussed effects of the magnetic fields on the heavy quark sector, which are expected to be observed by hard probes in the heavy-ion collisions. In Sec. 4.1, we first discussed the quarkonium spectra in magnetic fields, which has close connection to the recent lattice QCD simulations on the properties of the QCD vacuum, especially the confinement potential. Extension to the finite temperature and/or density is left for the future study. Finally, in Sec. 4.2, we discussed a recent study on the heavy-quark transport in a hot medium and magnetic field on the basis of the Langevin equation. The momentum diffusion constant was computed by using the joint resummation for Hard Thermal Loop and strong magnetic field. While the computation has been performed for the lowest Landau level in the strong-field limit so far, extension to other regimes, e.g., weak field, and inclusion of the higher Landau levels are relevant for the phenomenology in the heavy-ion physics. The topics covered by the present review are still fast developing, a number of important issues should be further investigated in future. Here we briefly discuss a few of them. 1) Although there are already thorough numerical simulations of the magnetic fields in the initial stage of the heavy-ion collisions, the precise time evolution of the magnetic fields is still unknown, see discussions in Sec. 2. The knowledge of the precise time evolution of the magnetic fields is very important for quantitative understanding of a number of magnetic-field induced phenomena, like the CME, CMW, etc. Thus detailed simulations of the time dependence of the magnetic field by using, for example, the parton cascade models or by using magnetohydrodynamics are very desirable. 2) The experimental measurements of the CME and CMW contain significant background contributions, see Sec. 3; to extract the CME or CMW signals we have to understand these background contributions at the quantitative level which still remains an unsolved problem and will be one of the most important directions to pursue in this field. 3) Heavy flavors will serve as an alternative probe to explore the dynamics in the strong magnetic field and to constrain the spacetime profile of the magnetic field. Therefore, it will be important to develop the numerical modelling of the heavy-flavor dynamics with the inclusion of the effects of the magnetic field. It is also interesting to investigate other transport coefficients. For example, the electrical conductivity will be an important input for the magnetohydrodynamics, which in turn will be relevant for the estimate of the magnetic-field profile. The techniques and ingredients developed in the computation of the heavy-quark diffusion coefficients have been applied in a few works for the jet quenching parameter [ 237 237 ] and the electrical conductivity [ 238 238 , 239 239 ] in the strong-field limit. The other regime, regarding the strength of the magnetic field, still needs to be explored for the phenomenology of the heavy-ion collisions. References J. Rafelski , B. Muller , 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 D.N. Voskresensky , N. Yu. Anisimov , Properties of a pion condensate in a magnetic field . Sov. Phys. JETP 51 , 13 - 22 ( 1980 ). Properties of a pion condensate in a magnetic field S. Schramm , B. Muller , A.J. Schramm , Quark - anti-quark condensates in strong magnetic fields . Mod. Phys. Lett. A 7 , 973 - 982 ( 1992 ). doi: 10.1142/S0217732392000860 http://doi.org/10.1142/S0217732392000860 Quark - anti-quark condensates in strong magnetic fields S. Schramm , B. Muller , A.J. Schramm , Exact Dirac propagator in a magnetic ’sheet’ . Phys. Lett. A 164 , 28 ( 1992 ). doi: 10.1016/0375-9601(92)90900-7 http://doi.org/10.1016/0375-9601(92)90900-7 Exact Dirac propagator in a magnetic ’sheet’ 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 - 253 ( 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’ V. Skokov , A. Illarionov , V. Toneev , Estimate of the magnetic field strength in heavy-ion collisions . Int. J. Mod. Phys. A 24 , 5925 - 5932 ( 2009 ). doi: 10.1142/S0217751X09047570 http://doi.org/10.1142/S0217751X09047570 Estimate of the magnetic field strength in heavy-ion collisions V. Voronyuk , V.D. Toneev , W. Cassing , et al ., (Electro-)Magnetic 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 (Electro-)Magnetic field evolution in relativistic heavy-ion collisions A. Bzdak , V. Skokov , Event-by-event fluctuations of magnetic and electric fields in heavy ion collisions . Phys. Lett.B 710 , 171 - 174 ( 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 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 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 J. Bloczynski , X.G. Huang , X.L. Zhang , et al ., Azimuthally fluctuating magnetic field and its impacts on observables in heavy-ion collisions . Phys. Lett.B 718 , 1529 - 1535 ( 2013 ). doi: 10.1016/j.physletb.2012.12.030 http://doi.org/10.1016/j.physletb.2012.12.030 Azimuthally fluctuating magnetic field and its impacts on observables in heavy-ion collisions J. Bloczynski , X.G. Huang , X.L. Zhang , et al ., Charge-dependent azimuthal correlations from AuAu to UU collisions . Nucl. Phys. A 939 , 85 - 100 ( 2015 ). doi: 10.1016/j.nuclphysa.2015.03.012 http://doi.org/10.1016/j.nuclphysa.2015.03.012 Charge-dependent azimuthal correlations from AuAu to UU collisions Y. Zhong , C.-B. Yang , X. Cai , et al ., A systematic study of magnetic field in Relativistic Heavy-ion Collisions in the RHIC and LHC energy regions . Adv. High Energy Phys . 2014 , 193039 ( 2014 ). doi: 10.1155/2014/193039 http://doi.org/10.1155/2014/193039 A systematic study of magnetic field in Relativistic Heavy-ion Collisions in the RHIC and LHC energy regions Y. Zhong , C.-B. Yang , X. Cai , et al ., Spatial distributions of magnetic field in the RHIC and LHC energy regions . Chin. Phys. C 39 , 104105 ( 2015 ). doi: 10.1088/1674-1137/39/10/104105 http://doi.org/10.1088/1674-1137/39/10/104105 Spatial distributions of magnetic field in the RHIC and LHC energy regions H. Li , X.L. Sheng , Q. Wang , Electromagnetic fields with electric and chiral magnetic conductivities in heavy ion collisions . ( 2016 ). arXiv: 1602.02223 Electromagnetic fields with electric and chiral magnetic conductivities in heavy ion collisions R. Holliday , K. Tuchin , Classical electromagnetic fields from quantum sources in Heavy-Ion collisions . ( 2016 ). arXiv:1604.04572 Classical electromagnetic fields from quantum sources in Heavy-Ion collisions K. Fukushima , D.E. Kharzeev , H.J. Warringa , The chiral magnetic effect . Phys. Rev. D 78 , 074033 ( 2008 ). doi: 10.1103/PhysRevD.78.074033 http://doi.org/10.1103/PhysRevD.78.074033 The chiral magnetic effect D.T. Son , A.R. Zhitnitsky , Quantum anomalies in dense matter . Phys. Rev. D 70 , 074018 ( 2004 ). doi: 10.1103/PhysRevD.70.074018 http://doi.org/10.1103/PhysRevD.70.074018 Quantum anomalies in dense matter M.A. Metlitski , A.R. Zhitnitsky , Anomalous axion interactions and topological currents in dense matter . Phys. Rev. D 72 , 045011 ( 2005 ). doi: 10.1103/PhysRevD.72.045011 http://doi.org/10.1103/PhysRevD.72.045011 Anomalous axion interactions and topological currents in dense matter D.E. Kharzeev , H.-U. Yee , Chiral magnetic wave . Phys. Rev. D 83 , 085007 ( 2011 ). doi: 10.1103/PhysRevD.83.085007 http://doi.org/10.1103/PhysRevD.83.085007 Chiral magnetic wave E.V. Gorbar , V.A. Miransky , I.A. Shovkovy , Normal ground state of dense relativistic matter in a magnetic field . Phys. Rev. D 83 , 085003 ( 2011 ). doi: 10.1103/PhysRevD.83.085003 http://doi.org/10.1103/PhysRevD.83.085003 Normal ground state of dense relativistic matter in a magnetic field Y. Burnier , D.E. Kharzeev , J.F. Liao , et al ., Chiral magnetic wave at finite baryon density and the electric quadrupole moment of quark-gluon plasma in heavy ion collisions . Phys. Rev. Lett . 107 , 052303 ( 2011 ). doi: 10.1103/PhysRevLett.107.052303 http://doi.org/10.1103/PhysRevLett.107.052303 Chiral magnetic wave at finite baryon density and the electric quadrupole moment of quark-gluon plasma in heavy ion collisions K. Marasinghe , K. Tuchin , Quarkonium dissociation in quark-gluon plasma via ionization in magnetic field . Phys. Rev. C 84 , 044908 ( 2011 ). doi: 10.1103/PhysRevC.84.044908 http://doi.org/10.1103/PhysRevC.84.044908 Quarkonium dissociation in quark-gluon plasma via ionization in magnetic field D.L. Yang , B. Muller , J/ψ Production by Magnetic Excitation of ηc . J. Phys. G 39 , 015007 ( 2012 ). doi: 10.1088/0954-3899/39/1/015007 http://doi.org/10.1088/0954-3899/39/1/015007 J/ψ Production by Magnetic Excitation of ηc C.S. Machado , F.S. Navarra , E.G. de Oliveira , et al ., Heavy quarkonium production in a strong magnetic field . Phys. Rev. D 88 , 034009 ( 2013 ). doi: 10.1103/PhysRevD.88.034009 http://doi.org/10.1103/PhysRevD.88.034009 Heavy quarkonium production in a strong magnetic field C.S. Machado , S.I. Finazzo , R.D. Matheus , et al ., Modification of the B meson mass in a magnetic field from QCD sum rules . Phys. Rev. D 89 , 074027 ( 2014 ). doi: 10.1103/PhysRevD.89.074027 http://doi.org/10.1103/PhysRevD.89.074027 Modification of the B meson mass in a magnetic field from QCD sum rules J. Alford , M. Strickland , Charmonia and bottomonia in a magnetic field . Phys. Rev. D 88 , 105017 ( 2013 ). doi: 10.1103/PhysRevD.88.105017 http://doi.org/10.1103/PhysRevD.88.105017 Charmonia and bottomonia in a magnetic field Y.P. Liu , C. Greiner , C.M. Ko , Spin Asymmetry of J/psi in Peripheral Pb+Pb Collisions at LHC . ( 2014 ). arXiv: 1403.4317 Spin Asymmetry of J/psi in Peripheral Pb+Pb Collisions at LHC S. Cho , K. Hattori , S.H. Lee , et al . QCD sum rules for magnetically induced mixing between ηc and J/ψ . Phys. Rev. Lett . 113 , 172301 ( 2014 ). doi: 10.1103/PhysRevLett.113.172301 http://doi.org/10.1103/PhysRevLett.113.172301 QCD sum rules for magnetically induced mixing between ηc and J/ψ S. Cho , K. Hattori , S.H. Lee , et al . Charmonium spectroscopy in strong magnetic fields by QCD sum rules: S-wave ground states . Phys. Rev. D 91 , 045025 ( 2015 ). doi: 10.1103/PhysRevD.91.045025 http://doi.org/10.1103/PhysRevD.91.045025 Charmonium spectroscopy in strong magnetic fields by QCD sum rules: S-wave ground states P. Gubler , K. Hattori , S.H. Lee , et al . D mesons in a magnetic field . Phys. Rev. D 93 , 054026 ( 2016 ). doi: 10.1103/PhysRevD.93.054026 http://doi.org/10.1103/PhysRevD.93.054026 D mesons in a magnetic field C. Bonati , M. D’Elia , M. Mariti , et al ., Anisotropy of the quark-antiquark potential in a magnetic field . Phys. Rev. D 89 , 114502 ( 2014 ). doi: 10.1103/PhysRevD.89.114502 http://doi.org/10.1103/PhysRevD.89.114502 Anisotropy of the quark-antiquark potential in a magnetic field C. Bonati , M. D’Elia , Massimo , A. Rucci , Heavy quarkonia in strong magnetic fields . Phys. Rev. D 92 , 054014 ( 2015 ). doi: 10.1103/PhysRevD.92.054014 http://doi.org/10.1103/PhysRevD.92.054014 Heavy quarkonia in strong magnetic fields C. Bonati , M. D’Elia , M. Mariti , et al ., Magnetic field effects on the static quark potential at zero and finite temperature . Phys. Rev. D 94 , 094007 ( 2016 ). doi: 10.1103/PhysRevD.94.094007 http://doi.org/10.1103/PhysRevD.94.094007 Magnetic field effects on the static quark potential at zero and finite temperature R. Rougemont , R. Critelli , J. Noronha , Anisotropic heavy quark potential in strongly-coupled 𝒩=4 SYM in a magnetic field . Phys. Rev. D 91 , 066001 ( 2015 ). doi: 10.1103/PhysRevD.91.066001 http://doi.org/10.1103/PhysRevD.91.066001 Anisotropic heavy quark potential in strongly-coupled 𝒩=4 SYM in a magnetic field D. Dudal , T.G. Mertens , Melting of charmonium in a magnetic field from an effective AdS/QCD model . Phys. Rev. D 91 , 086002 ( 2015 ). doi: 10.1103/PhysRevD.91.086002 http://doi.org/10.1103/PhysRevD.91.086002 Melting of charmonium in a magnetic field from an effective AdS/QCD model X.Y. Guo , S.Z. Shi , N. Xu , et al ., Magnetic field effect on charmonium production in high energy nuclear collisions . Phys. Lett. B 751 , 215 - 219 ( 2015 ). doi: 10.1016/j.physletb.2015.10.038 http://doi.org/10.1016/j.physletb.2015.10.038 Magnetic field effect on charmonium production in high energy nuclear collisions K. Suzuki , T. Yoshida , Cigar-shaped quarkonia under strong magnetic field . Phys. Rev. D 93 , 051502 ( 2016 ). doi: 10.1103/PhysRevD.93.051502 http://doi.org/10.1103/PhysRevD.93.051502 Cigar-shaped quarkonia under strong magnetic field T. Yoshida , K. Suzuki , Heavy meson spectroscopy under strong magnetic field . Phys. Rev. D 94 , 074043 ( 2016 ). doi: 10.1103/PhysRevD.94.074043 http://doi.org/10.1103/PhysRevD.94.074043 Heavy meson spectroscopy under strong magnetic field K. Fukushima , K. Hattori , Koichi , H.-U. Yee , et al ., Heavy quark diffusion in strong magnetic fields at weak coupling and implications for elliptic flow . Phys. Rev. D 93 , 074028 ( 2016 ). doi: 10.1103/PhysRevD.93.074028 http://doi.org/10.1103/PhysRevD.93.074028 Heavy quark diffusion in strong magnetic fields at weak coupling and implications for elliptic flow S. I. Finazzo , R. Critelli , R. Rougemont , et al ., Momentum transport in strongly coupled anisotropic plasmas in the presence of strong magnetic fields . ( 2016 ). arXiv: 1605.06061 Momentum transport in strongly coupled anisotropic plasmas in the presence of strong magnetic fields S.K. Das , S. Plumari , S. Chatterjee , et al ., Directed flow of charm quarks induced by the initial magnetic field at LHC . ( 2016 ) arXiv: 1608.02231 Directed flow of charm quarks induced by the initial magnetic field at LHC V.P. Gusynin , V.A. Miransky , I.A. Shovkovy , Dimensional reduction and dynamical chiral symmetry breaking by a magnetic field in (3+1)-dimensions . Phys. Lett.B 349 , 477 - 483 ( 1995 ). doi: 10.1016/0370-2693(95)00232-A http://doi.org/10.1016/0370-2693(95)00232-A Dimensional reduction and dynamical chiral symmetry breaking by a magnetic field in (3+1)-dimensions V.P. Gusynin , V.A. Miransky , I.A. Shovkovy , Catalysis of dynamical flavor symmetry breaking by a magnetic field in (2+1)-dimensions . Phys. Rev. Lett . 73 , 3499 - 3502 ( 1994 ). doi: 10.1103/PhysRevLett.73.3499 http://doi.org/10.1103/PhysRevLett.73.3499 Catalysis of dynamical flavor symmetry breaking by a magnetic field in (2+1)-dimensions V.P. Gusynin , V.A. Miransky , I.A. Shovkovy , Dimensional reduction and catalysis of dynamical symmetry breaking by a magnetic field . Nucl. Phys. B 462 , 249 - 290 ( 1996 ). doi: 10.1016/0550-3213(96)00021-1 http://doi.org/10.1016/0550-3213(96)00021-1 Dimensional reduction and catalysis of dynamical symmetry breaking by a magnetic field I.A. Shovkovy , Magnetic catalysis: A review . Lect. Notes Phys . 871 , 13 - 49 ( 2013 ). doi: 10.1007/978-3-642-37305-3_2 http://doi.org/10.1007/978-3-642-37305-3_2 Magnetic catalysis: A review F. Preis , A. Rebhan , A. Schmitt , Inverse magnetic catalysis in dense holographic matter . JHEP 03 , 033 ( 2011 ). doi: 10.1007/JHEP03(2011)033 http://doi.org/10.1007/JHEP03(2011)033 Inverse magnetic catalysis in dense holographic matter G.S. Bali , F. Bruckmann , G. Endrodi , et al ., The QCD phase diagram for external magnetic fields . JHEP 02 , 044 ( 2012 ). doi: 10.1007/JHEP02(2012)044 http://doi.org/10.1007/JHEP02(2012)044 The QCD phase diagram for external magnetic fields G.S. Bali , F. Bruckmann , G. Endrodi , et al ., QCD quark condensate in external magnetic fields . Phys. Rev.D 86 , 071502 ( 2012 ). doi: 10.1103/PhysRevD.86.071502 http://doi.org/10.1103/PhysRevD.86.071502 QCD quark condensate in external magnetic fields F. Bruckmann , G. Endrodi , T.G. Kovacs , Inverse magnetic catalysis and the Polyakov loop . JHEP 04 , 112 ( 2013 ). doi: 10.1007/JHEP04(2013)112 http://doi.org/10.1007/JHEP04(2013)112 Inverse magnetic catalysis and the Polyakov loop K. Fukushima , J.M. Pawlowski , Magnetic catalysis in hot and dense quark matter and quantum fluctuations . Phys. Rev. D 86 , 076013 ( 2012 ). doi: 10.1103/PhysRevD.86.076013 http://doi.org/10.1103/PhysRevD.86.076013 Magnetic catalysis in hot and dense quark matter and quantum fluctuations K. Fukushima , Y. Hidaka , Magnetic catalysis versus magnetic inhibition . Phys. Rev. Lett . 110 , 031601 ( 2013 ). doi: 10.1103/PhysRevLett.110.031601 http://doi.org/10.1103/PhysRevLett.110.031601 Magnetic catalysis versus magnetic inhibition T. Kojo , N. Su , The quark mass gap in a magnetic field . Phys. Lett.B 720 , 192 - 197 ( 2013 ). doi: 10.1016/j.physletb.2013.02.024 http://doi.org/10.1016/j.physletb.2013.02.024 The quark mass gap in a magnetic field J.Y. Chao , P.C. Chu , M. Huang , Inverse magnetic catalysis induced by sphalerons . Phys. Rev. D 88 , 054009 ( 2013 ). doi: 10.1103/PhysRevD.88.054009 http://doi.org/10.1103/PhysRevD.88.054009 Inverse magnetic catalysis induced by sphalerons L. Yu , H. Liu , M. Huang , Spontaneous generation of local CP violation and inverse magnetic catalysis . Phys. Rev. D 90 , 074009 ( 2014 ). doi: 10.1103/PhysRevD.90.074009 http://doi.org/10.1103/PhysRevD.90.074009 Spontaneous generation of local CP violation and inverse magnetic catalysis B. Feng , D.F. Hou , H.C. Ren , Magnetic and inverse magnetic catalysis in the Bose-Einstein condensation of neutral bound pairs . Phys. Rev. D 92 , 065011 ( 2015 ). doi: 10.1103/PhysRevD.92.065011 http://doi.org/10.1103/PhysRevD.92.065011 Magnetic and inverse magnetic catalysis in the Bose-Einstein condensation of neutral bound pairs L. Yu , J. Van Doorsselaere , M. Huang , Inverse magnetic catalysis in the three-flavor NJL model with axial-vector interaction . Phys. Rev. D 91 , 074011 ( 2015 ). doi: 10.1103/PhysRevD.91.074011 http://doi.org/10.1103/PhysRevD.91.074011 Inverse magnetic catalysis in the three-flavor NJL model with axial-vector interaction G.Q. Cao , L.Y. He , P.F. Zhuang , Collective modes and Kosterlitz-Thouless transition in a magnetic field in the planar Nambu-Jona-Lasino model . Phys. Rev. D 90 , 056005 ( 2014 ). doi: 10.1103/PhysRevD.90.056005 http://doi.org/10.1103/PhysRevD.90.056005 Collective modes and Kosterlitz-Thouless transition in a magnetic field in the planar Nambu-Jona-Lasino model E.J. Ferrer , V. de la Incera , X.J. Wen , Quark antiscreening at strong magnetic field and inverse magnetic catalysis . Phys. Rev. D 91 , 054006 ( 2015 ). doi: 10.1103/PhysRevD.91.054006 http://doi.org/10.1103/PhysRevD.91.054006 Quark antiscreening at strong magnetic field and inverse magnetic catalysis J. Braun , W.A. Mian , S. Rechenberger , Delayed magnetic catalysis . Phys. Lett.B 755 , 265 - 269 ( 2016 ). doi: 10.1016/j.physletb.2016.02.017 http://doi.org/10.1016/j.physletb.2016.02.017 Delayed magnetic catalysis N. Mueller , J.M. Pawlowski , Magnetic catalysis and inverse magnetic catalysis in QCD . Phys. Rev. D 91 , 116010 ( 2015 ). doi: 10.1103/PhysRevD.91.116010 http://doi.org/10.1103/PhysRevD.91.116010 Magnetic catalysis and inverse magnetic catalysis in QCD K. Hattori , T. Kojo , N. Su , Mesons in strong magnetic fields: (I) General analyses . Nucl. Phys. A 951 , 1 - 30 ( 2016 ). doi: 10.1016/j.nuclphysa.2016.03.016 http://doi.org/10.1016/j.nuclphysa.2016.03.016 Mesons in strong magnetic fields: (I) General analyses M. Ruggieri , G.X. Peng , Quark matter in a parallel electric and magnetic field background: Chiral phase transition and equilibration of chiral density . Phys. Rev. D 93 , 094021 ( 2016 ). doi: 10.1103/PhysRevD.93.094021 http://doi.org/10.1103/PhysRevD.93.094021 Quark matter in a parallel electric and magnetic field background: Chiral phase transition and equilibration of chiral density H.L. Chen , K. Fukushima , X.G. Huang , et al ., Analogy between rotation and density for Dirac fermions in a magnetic field . Phys. Rev. D 93 , 104052 ( 2016 ). doi: 10.1103/PhysRevD.93.104052 http://doi.org/10.1103/PhysRevD.93.104052 Analogy between rotation and density for Dirac fermions in a magnetic field K. Hattori , Y. Yin , Charge redistribution from anomalous magneto-vorticity coupling . ( 2016 ). arXivi: 1607.01513 Charge redistribution from anomalous magneto-vorticity coupling S. Ebihara , K. Fukushima , K. Mameda , Boundary effects and gapped dispersion in rotating fermionic matter . ( 2016 ). arXivi: 1608.00336 Boundary effects and gapped dispersion in rotating fermionic matter M.N. Chernodub , Superconductivity of QCD vacuum in strong magnetic field . Phys. Rev. D 82 , 085011 ( 2010 ). doi: 10.1103/PhysRevD.82.085011 http://doi.org/10.1103/PhysRevD.82.085011 Superconductivity of QCD vacuum in strong magnetic field M.N. Chernodub , Spontaneous electromagnetic superconductivity of vacuum in strong magnetic field: evidence from the Nambu-Jona-Lasinio model . Phys. Rev. Lett . 106 , 142003 ( 2011 ). doi: 10.1103/PhysRevLett.106.142003 http://doi.org/10.1103/PhysRevLett.106.142003 Spontaneous electromagnetic superconductivity of vacuum in strong magnetic field: evidence from the Nambu-Jona-Lasinio model Y. Hidaka , A. Yamamoto , Charged vector mesons in a strong magnetic field . Phys. Rev. D 87 , 094502 ( 2013 ). doi: 10.1103/PhysRevD.87.094502 http://doi.org/10.1103/PhysRevD.87.094502 Charged vector mesons in a strong magnetic field H. Liu , L. Yu , M. Huang , Charged and neutral vector ρ mesons in a magnetic field . Phys. Rev. D 91 , 014017 ( 2015 ). doi: 10.1103/PhysRevD.91.014017 http://doi.org/10.1103/PhysRevD.91.014017 Charged and neutral vector ρ mesons in a magnetic field H. Liu , L. Yu , M. Huang , Survival of charged rho condensation at high temperature and density . Chin. Phys. C 40 , 023102 ( 2016 ). doi: 10.1088/1674-1137/40/2/023102 http://doi.org/10.1088/1674-1137/40/2/023102 Survival of charged rho condensation at high temperature and density G.Q. Cao , X.G. Huang , Electromagnetic triangle anomaly and neutral pion condensation in QCD vacuum . Phys. Lett. B 757 , 1 - 5 ( 2016 ). doi: 10.1016/j.physletb.2016.03.066 http://doi.org/10.1016/j.physletb.2016.03.066 Electromagnetic triangle anomaly and neutral pion condensation in QCD vacuum S.I. Braginskii , Transport processes in a plasma . in Reviews of Plasma Physics , Vol. 1 , Consultants Bureau , New York, NY , p. 205 ( 1965 ) Transport processes in a plasma E.M. Lifshitz , L.P. Pitaevskii , Transport processes in a plasma . Physcial Kinetics, Course of Theoretical Physics , Vol. 10 , Pergamon , New York , ( 1981 ) Transport processes in a plasma X.G. Huang , M. Huang , D.H. Rischke , et al ., Anisotropic hydrodynamics, bulk viscosities and R-modes of strange quark stars with strong magnetic fields . Phys. Rev.D 81 , 045015 ( 2010 ). doi: 10.1103/PhysRevD.81.045015 http://doi.org/10.1103/PhysRevD.81.045015 Anisotropic hydrodynamics, bulk viscosities and R-modes of strange quark stars with strong magnetic fields X.G. Huang , A. Sedrakian , D.H. Rischke , Kubo formulae for relativistic fluids in strong magnetic fields . Annals Phys . 326 , 3075 - 3094 ( 2011 ). doi: 10.1016/j.aop.2011.08.001 http://doi.org/10.1016/j.aop.2011.08.001 Kubo formulae for relativistic fluids in strong magnetic fields K. Tuchin , On viscous flow and azimuthal anisotropy of quark-gluon plasma in strong magnetic field . J. Phys. G 39 , 025010 ( 2012 ). doi: 10.1088/0954-3899/39/2/025010 http://doi.org/10.1088/0954-3899/39/2/025010 On viscous flow and azimuthal anisotropy of quark-gluon plasma in strong magnetic field S.I. Finazzo , R. Rougemont , H. Marrochio , et al ., Hydrodynamic transport coefficients for the non-conformal quark-gluon plasma from holography . JHEP , 02 , 051 ( 2015 ). doi: 10.1007/JHEP02(2015)051 http://doi.org/10.1007/JHEP02(2015)051 Hydrodynamic transport coefficients for the non-conformal quark-gluon plasma from holography K. Tuchin , Synchrotron radiation by fast fermions in heavy-ion collisions . Phys. Rev. C 82 , 034904 ( 2010 ). doi: 10.1103/PhysRevC.82.034904 http://doi.org/10.1103/PhysRevC.82.034904 Synchrotron radiation by fast fermions in heavy-ion collisions K. Tuchin , Photon decay in strong magnetic field in heavy-ion collisions . Phys. Rev. C 83 , 017901 ( 2011 ). doi: 10.1103/PhysRevC.83.017901 http://doi.org/10.1103/PhysRevC.83.017901 Photon decay in strong magnetic field in heavy-ion collisions K. Tuchin , Electromagnetic radiation by quark-gluon plasma in a magnetic field . Phys. Rev. C 87 , 024912 ( 2013 ). doi: 10.1103/PhysRevC.87.024912 http://doi.org/10.1103/PhysRevC.87.024912 Electromagnetic radiation by quark-gluon plasma in a magnetic field K. Tuchin , Particle production in strong electromagnetic fields in relativistic heavy-ion collisions . Adv. High Energy Phys . 2013 , 490495 ( 2013 ). doi: 10.1155/2013/490495 http://doi.org/10.1155/2013/490495 Particle production in strong electromagnetic fields in relativistic heavy-ion collisions K. Tuchin , Synchrotron radiation of vector bosons at relativistic colliders . Phys. Rev. D 91 , 033004 ( 2015 ). doi: 10.1103/PhysRevD.91.033004 http://doi.org/10.1103/PhysRevD.91.033004 Synchrotron radiation of vector bosons at relativistic colliders K. Tuchin , Role of magnetic field in photon excess in heavy ion collisions . Phys. Rev. C 91 , 014902 ( 2015 ). doi: 10.1103/PhysRevC.91.014902 http://doi.org/10.1103/PhysRevC.91.014902 Role of magnetic field in photon excess in heavy ion collisions G. Basar , D. Kharzeev , D. Kharzeev , et al ., Conformal anomaly as a source of soft photons in heavy ion collisions . Phys. Rev. Lett . 109 , 202303 ( 2012 ). doi: 10.1103/PhysRevLett.109.202303 http://doi.org/10.1103/PhysRevLett.109.202303 Conformal anomaly as a source of soft photons in heavy ion collisions K. Fukushima , K. Mameda , Wess-Zumino-Witten action and photons from the Chiral Magnetic Effect . Phys. Rev.D 86 071501 ( 2012 ). doi: 10.1103/PhysRevD.86.071501 http://doi.org/10.1103/PhysRevD.86.071501 Wess-Zumino-Witten action and photons from the Chiral Magnetic Effect K. Hattori , K. Itakura , S. Ozaki , Neutral-pion reactions induced by chiral anomaly in strong magnetic fields .( 2013 ). arXiv: 1305.7224 Neutral-pion reactions induced by chiral anomaly in strong magnetic fields A. Ayala , J.D. Castano-Yepes , C.A. Dominguez , et al ., Thermal photon production from gluon fusion induced by magnetic fields in relativistic heavy-ion collisions .( 2016 ). arXiv: 1604.02713 Thermal photon production from gluon fusion induced by magnetic fields in relativistic heavy-ion collisions B.G. Zakharov , Effect of magnetic field on the photon radiation from quark-gluon plasma in heavy ion collisions .( 2016 ). arXiv: 1609.04324 Effect of magnetic field on the photon radiation from quark-gluon plasma in heavy ion collisions D. Kharzeev , K. Landsteiner , A. Schmitt , et al . Strongly interacting matter in magnetic fields . Lect. Notes Phys . 871 , 1 - 624 ( 2013 ). doi: 10.1007/978-3-642-37305-3 http://doi.org/10.1007/978-3-642-37305-3 Strongly interacting matter in magnetic fields J.O. Andersen , W.R. Naylor , A. Tranberg , Phase diagram of QCD in a magnetic field: A review . Rev. Mod. Phys . 88 , 025001 ( 2016 ). doi: 10.1103/RevModPhys.88.025001 http://doi.org/10.1103/RevModPhys.88.025001 Phase diagram of QCD in a magnetic field: A review V.A. Miransky , I.A. Shovkovy , Quantum field theory in a magnetic field: From quantum chromodynamics to graphene and Dirac semimetals . Phys. Rept . 576 , 1 - 209 ( 2015 ). doi: 10.1016/j.physrep.2015.02.003 http://doi.org/10.1016/j.physrep.2015.02.003 Quantum field theory in a magnetic field: From quantum chromodynamics to graphene and Dirac semimetals X.G. Huang , Electromagnetic fields and anomalous transports in heavy-ion collisions — A pedagogical review . Rept. Prog. Phys . 79 , 076302 ( 2016 ). doi: 10.1088/0034-4885/79/7/076302 http://doi.org/10.1088/0034-4885/79/7/076302 Electromagnetic fields and anomalous transports in heavy-ion collisions — A pedagogical review S.A. Olausen , V.M. Kaspi , The McGill magnetar catalog . Astrophys. J. Suppl . 212 , 6 ( 2014 ). doi: 10.1088/0067-0049/212/1/6 http://doi.org/10.1088/0067-0049/212/1/6 The McGill magnetar catalog R. Turolla , S. Zane , A. Watts , Magnetars: the physics behind observations . A review. Rept. Prog. Phys . 78 , 116901 ( 2015 ). doi: 10.1088/0034-4885/78/11/116901 http://doi.org/10.1088/0034-4885/78/11/116901 Magnetars: the physics behind observations P.B. Arnold , G.D. Moore , L.G. Yaffe , Transport coefficients in high temperature gauge theories. 2. Beyond leading log . JHEP 05 , 051 ( 2003 ). doi: 10.1088/1126-6708/2003/05/051 http://doi.org/10.1088/1126-6708/2003/05/051 Transport coefficients in high temperature gauge theories. 2. Beyond leading log S. Gupta , The Electrical conductivity and soft photon emissivity of the QCD plasma . Phys. Lett. B 597 , 57 - 62 ( 2004 ). doi: 10.1016/j.physletb.2004.05.079 http://doi.org/10.1016/j.physletb.2004.05.079 The Electrical conductivity and soft photon emissivity of the QCD plasma G. Aarts , C. Allton , J. Foley , et al ., Spectral functions at small energies and the electrical conductivity in hot, quenched lattice QCD . Phys. Rev. Lett . 99 , 022002 ( 2007 ). doi: 10.1103/PhysRevLett.99.022002 http://doi.org/10.1103/PhysRevLett.99.022002 Spectral functions at small energies and the electrical conductivity in hot, quenched lattice QCD H. -T. Ding , A. Francis , O. Kaczmarek , et al ., Thermal dilepton rate and electrical conductivity: An analysis of vector current correlation functions in quenched lattice QCD . Phys. Rev. D 83 , 034504 ( 2011 ). doi: 10.1103/PhysRevD.83.034504 http://doi.org/10.1103/PhysRevD.83.034504 Thermal dilepton rate and electrical conductivity: An analysis of vector current correlation functions in quenched lattice QCD A. Francis , O. Kaczmarek , On the temperature dependence of the electrical conductivity in hot quenched lattice QCD . Prog. Part. Nucl. Phys . 67 , 212 - 217 ( 2012 ). doi: 10.1016/j.ppnp.2011.12.020 http://doi.org/10.1016/j.ppnp.2011.12.020 On the temperature dependence of the electrical conductivity in hot quenched lattice QCD H.-T. Ding , O. Kaczmarek , F. Meyer , Vector spectral functions and transport properties in quenched QCD . Proceedings, 32nd International Symposium on Lattice Field Theory (Lattice 2014) : Brookhaven, NY, USA , June 23-28, 2014 . Vector spectral functions and transport properties in quenched QCD B.B. Brandt , A. Francis , H.B. Meyer , et al ., Thermal Correlators in the ρ channel of two-flavor QCD . JHEP 03 , 100 ( 2013 ). doi: 10.1007/JHEP03(2013)100 http://doi.org/10.1007/JHEP03(2013)100 Thermal Correlators in the ρ channel of two-flavor QCD A. Amato , G. Aarts , C. Allton , et al ., Electrical conductivity of the quark-gluon plasma across the deconfinement transition . Phys. Rev. Lett . 111 , 172001 ( 2013 ). doi: 10.1103/PhysRevLett.111.172001 http://doi.org/10.1103/PhysRevLett.111.172001 Electrical conductivity of the quark-gluon plasma across the deconfinement transition G. Aarts , C. Allton , A. Amato , et al ., Electrical conductivity and charge diffusion in thermal QCD from the lattice . JHEP 02 , 186 ( 2015 ). doi: 10.1007/JHEP02(2015)186 http://doi.org/10.1007/JHEP02(2015)186 Electrical conductivity and charge diffusion in thermal QCD from the lattice H.-T. Ding , O. Kaczmarek , F. Meyer , Thermal dilepton rates and electrical conductivity of the QGP from the lattice . Phys. Rev. D 94 , 034504 ( 2016 ). doi: 10.1103/PhysRevD.94.034504 http://doi.org/10.1103/PhysRevD.94.034504 Thermal dilepton rates and electrical conductivity of the QGP from the lattice H.-T. Ding , F. Karsch , S. Mukherjee , Thermodynamics of strong-interaction matter from Lattice QCD . Int. J. Mod. Phys. E 24 , 1530007 ( 2015 ). doi: 10.1142/S0218301315300076 http://doi.org/10.1142/S0218301315300076 Thermodynamics of strong-interaction matter from Lattice QCD K. Tuchin , Time and space dependence of the electromagnetic field in relativistic heavy-ion collisions . Phys. Rev. C 88 , 024911 ( 2013 ). doi: 10.1103/PhysRevC.88.024911 http://doi.org/10.1103/PhysRevC.88.024911 Time and space dependence of the electromagnetic field in relativistic heavy-ion collisions K. Tuchin , Electromagnetic field and the chiral magnetic effect in the quark-gluon plasma . Phys. Rev. C 91 , 064902 ( 2015 ). doi: 10.1103/PhysRevC.91.064902 http://doi.org/10.1103/PhysRevC.91.064902 Electromagnetic field and the chiral magnetic effect in the quark-gluon plasma K. Tuchin , Electromagnetic fields in high energy heavy-ion collisions . Int. J. Mod. Phys. E 23 , 1430001 ( 2014 ). doi: 10.1142/S021830131430001X http://doi.org/10.1142/S021830131430001X Electromagnetic fields in high energy heavy-ion collisions K. Tuchin , Initial value problem for magnetic fields in heavy ion collisions . Int. J. Mod. Phys. Rev. C 93 , 014905 ( 2016 ). doi: 10.1103/PhysRevC.93.014905 http://doi.org/10.1103/PhysRevC.93.014905 Initial value problem for magnetic fields in heavy ion collisions U. Gursoy , D. Kharzeev , K. Rajagopal , Magnetohydrodynamics, charged currents and directed flow in heavy ion collisions . Phys. Rev. C 89 , 054905 ( 2014 ). doi: 10.1103/PhysRevC.89.054905 http://doi.org/10.1103/PhysRevC.89.054905 Magnetohydrodynamics, charged currents and directed flow in heavy ion collisions B. G. Zakharov , Electromagnetic response of quark?Cgluon plasma in heavy-ion collisions . Phys. Lett. B 737 , 262 - 266 ( 2014 ). doi: 10.1016/j.physletb.2014.08.068 http://doi.org/10.1016/j.physletb.2014.08.068 Electromagnetic response of quark?Cgluon plasma in heavy-ion collisions L. McLerran , V. Skokov , Comments about the electromagnetic field in heavy-ion collisions . Nucl. Phys. A 929 , 184 - 190 ( 2014 ). doi: 10.1016/j.nuclphysa.2014.05.008 http://doi.org/10.1016/j.nuclphysa.2014.05.008 Comments about the electromagnetic field in heavy-ion collisions Y. Hirono , M. Hongo , T. Hirano , Estimation of electric conductivity of the quark gluon plasma via asymmetric heavy-ion collisions . Phys. Rev. C 90 , 021903 ( 2014 ). doi: 10.1103/PhysRevC.90.021903 http://doi.org/10.1103/PhysRevC.90.021903 Estimation of electric conductivity of the quark gluon plasma via asymmetric heavy-ion collisions W.-T. Deng , X.-G. Huang , Electric fields and chiral magnetic effect in Cu+Au collisions . Phys. Lett. B 742 , 296 - 302 ( 2015 ). doi: 10.1016/j.physletb.2015.01.050 http://doi.org/10.1016/j.physletb.2015.01.050 Electric fields and chiral magnetic effect in Cu+Au collisions V. Voronyuk , V.D. Toneev , S.A. Voloshin , et al ., Charge-dependent directed flow in asymmetric nuclear collisions . Phys. Rev. C 90 , 064903 ( 2014 ). doi: 10.1103/PhysRevC.90.064903 http://doi.org/10.1103/PhysRevC.90.064903 Charge-dependent directed flow in asymmetric nuclear collisions V. Toneev , O. Rogachevsky , V. Voronyuk , Evidence for creation of strong electromagnetic fields in relativistic heavy-ion collisions . ( 2016 ). arXiv: 1604.06231 Evidence for creation of strong electromagnetic fields in relativistic heavy-ion collisions S. Chatterjee , P. Tribedy , Separation of flow from the chiral magnetic effect in U + U collisions using spectator asymmetry . Phys. Rev. C 92 , 011902 ( 2015 ). doi: 10.1103/PhysRevC.92.011902 http://doi.org/10.1103/PhysRevC.92.011902 Separation of flow from the chiral magnetic effect in U + U collisions using spectator asymmetry A. Rybicki , A. Szczurek , M. Klusek-Gawenda , Charge splitting of directed flow and charge-dependent effects in pion spectra in heavy ion collisions . Proceedings, 10th Workshop on Particle Correlations and Femtoscopy (WPCF 2014) : Gyongyos, Hungary , August 25-29, 2014 . Charge splitting of directed flow and charge-dependent effects in pion spectra in heavy ion collisions L. Adamczyk , et al ., Charge-dependent directed flow in Cu+Au collisions at <math id="M150"><mrow><msqrt><mrow><msub><mi>s</mi><mrow><msub><mrow/><mrow><mi>N</mi><mi>N</mi></mrow></msub></mrow></msub></mrow></msqrt></mrow></math> = 200 GeV . ( 2016 ). arXiv: 1608.04100 Charge-dependent directed flow in Cu+Au collisions at sNN = 200 GeV D.E. Kharzeev , Topologically induced local P and CP violation in QCD x QED . Annals Phys . 325 , 205 - 218 ( 2010 ). doi: 10.1016/j.aop.2009.11.002 http://doi.org/10.1016/j.aop.2009.11.002 Topologically induced local P and CP violation in QCD x QED D.E. Kharzeev , The chiral magnetic eEffect and anomaly-iInduced transport . Prog. Part. Nucl. Phys . 75 , 133 - 151 ( 2014 ). doi: 10.1016/j.ppnp.2014.01.002 http://doi.org/10.1016/j.ppnp.2014.01.002 The chiral magnetic eEffect and anomaly-iInduced transport K. Fukushima , Views of the chiral magnetic effect . Lect. Notes Phys . 871 , 241 - 259 ( 2013 ). doi: 10.1007/978-3-642-37305-3_9 http://doi.org/10.1007/978-3-642-37305-3_9 Views of the chiral magnetic effect D.E. Kharzeev , Topology, magnetic field, and strongly interacting matter . Ann. Rev. Nucl. Part. Sci . 65 , 193 - 214 ( 2015 ). doi: 10.1146/annurev-nucl-102313-025420 http://doi.org/10.1146/annurev-nucl-102313-025420 Topology, magnetic field, and strongly interacting matter J.F. Liao , Anomalous transport effects and possible environmental symmetry 'violation' in heavy-ion collisions . Pramana 84 , 901 - 926 ( 2015 ). doi: 10.1007/s12043-015-0984-x http://doi.org/10.1007/s12043-015-0984-x Anomalous transport effects and possible environmental symmetry 'violation' in heavy-ion collisions D.E. Kharzeev , J. Liao , S.A. Voloshin , et al ., Chiral magnetic and vortical effects in high-energy nuclear collisions? A status report . Prog. Part. Nucl. Phys . 88 , 1 - 28 ( 2016 ). doi: 10.1016/j.ppnp.2016.01.001 http://doi.org/10.1016/j.ppnp.2016.01.001 Chiral magnetic and vortical effects in high-energy nuclear collisions? A status report A. Vilenkin , Equilibrium parity violating current in a magnetic field . Phys. Rev. D 22 , 3080 - 3084 ( 1980 ). doi: 10.1103/PhysRevD.22.3080 http://doi.org/10.1103/PhysRevD.22.3080 Equilibrium parity violating current in a magnetic field Q. Li , D.E. Kharzeev , C. Zhang , et al ., Observation of the chiral magnetic effect in ZrTe5 . Nature Phys . 12 , 550 - 554 ( 2016 ). doi: 10.1038/nphys3648 http://doi.org/10.1038/nphys3648 Observation of the chiral magnetic effect in ZrTe5 J. Xiong , S.K. Kushwaha , T. Liang , et al ., Signature of the chiral anomaly in a Dirac semimetal: a current plume steered by a magnetic field . Science 350 , 413 ( 2015 ). Signature of the chiral anomaly in a Dirac semimetal: a current plume steered by a magnetic field F. Arnold , et al ., Negative magnetoresistance without well-defined chirality in the Weyl semimetal TaP . Nature Commun . 7 , 1615 ( 2016 ). doi: 10.1038/ncomms11615 http://doi.org/10.1038/ncomms11615 Negative magnetoresistance without well-defined chirality in the Weyl semimetal TaP X. Huang , et al ., Negative magnetoresistance without well-defined chirality in the Weyl semimetal TaP . Phys. Rev. X 5 , 031023 ( 2015 ). arXiv: 1506.06577 Negative magnetoresistance without well-defined chirality in the Weyl semimetal TaP M. A. Zubkov , Absence of equilibrium chiral magnetic effect . Phys. Rev. D 93 , 105036 ( 2016 ). doi: 10.1103/PhysRevD.93.105036 http://doi.org/10.1103/PhysRevD.93.105036 Absence of equilibrium chiral magnetic effect D. Kharzeev , Parity violation in hot QCD: Why it can happen, and how to look for it . Phys. Lett. B 633 , 260 - 264 ( 2006 ). doi: 10.1016/j.physletb.2005.11.075 http://doi.org/10.1016/j.physletb.2005.11.075 Parity violation in hot QCD: Why it can happen, and how to look for it I. Iatrakis , S. Lin , Y. Yin , The anomalous transport of axial charge: topological vs non-topological fluctuations . JHEP 09 , 030 ( 2015 ). doi: 10.1007/JHEP09(2015)030 http://doi.org/10.1007/JHEP09(2015)030 The anomalous transport of axial charge: topological vs non-topological fluctuations N. S. Manton , Topology in the Weinberg-Salam theory . Phys. Rev. D 28 , 2019 ( 1983 ). doi: 10.1103/PhysRevD.28.2019 http://doi.org/10.1103/PhysRevD.28.2019 Topology in the Weinberg-Salam theory F.R. Klinkhamer , N.S. Manton , A saddle point solution in the Weinberg-Salam theory . Phys. Rev. D 30 , 2212 ( 1984 ). doi: 10.1103/PhysRevD.30.2212 http://doi.org/10.1103/PhysRevD.30.2212 A saddle point solution in the Weinberg-Salam theory D. Kharzeev , A. Krasnitz , R. Venugopalan , Anomalous chirality fluctuations in the initial stage of heavy ion collisions and parity odd bubbles . Phys. Lett. B 545 , 298 - 306 ( 2002 ). doi: 10.1016/S0370-2693(02)02630-8 http://doi.org/10.1016/S0370-2693(02)02630-8 Anomalous chirality fluctuations in the initial stage of heavy ion collisions and parity odd bubbles T. Lappi , L. McLerran , Some features of the glasma . Nucl. Phys. A 772 , 200 - 212 ( 2006 ). doi: 10.1016/j.nuclphysa.2006.04.001 http://doi.org/10.1016/j.nuclphysa.2006.04.001 Some features of the glasma Y. Hirono , T. Hirano , D.E. Kharzeev , The chiral magnetic effect in heavy-ion collisions from event-by-event anomalous hydrodynamics . ( 2014 ). arXiv: 1412.0311 The chiral magnetic effect in heavy-ion collisions from event-by-event anomalous hydrodynamics S.A. Voloshin , Parity violation in hot QCD: How to detect it . Phys. Rev. C 70 , 057901 ( 2004 ). doi: 10.1103/PhysRevC.70.057901 http://doi.org/10.1103/PhysRevC.70.057901 Parity violation in hot QCD: How to detect it B. I. Abelev , et al ., 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 , et al ., 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 , et al ., Fluctuations of charge separation perpendicular to the event plane and local parity violation in <math id="M151"><mrow><msqrt><mrow><msub><mi>s</mi><mrow><mi>N</mi><mi>N</mi></mrow></msub></mrow></msqrt><mo>=</mo><mn>200</mn></mrow></math> GeV Au+Au collisions at the BNL Relativistic Heavy Ion Collider . Phys. Rev. C 88 , 064911 ( 2013 ). doi: 10.1103/PhysRevC.88.064911 http://doi.org/10.1103/PhysRevC.88.064911 Fluctuations of charge separation perpendicular to the event plane and local parity violation in sNN=200 GeV Au+Au collisions at the BNL Relativistic Heavy Ion Collider N.N. Ajitanand , R.A. Lacey , A. Taranenko , et al ., A New method for the experimental study of topological effects in the quark-gluon plasma . Phys. Rev. C 83 , 011901 ( 2011 ). doi: 10.1103/PhysRevC.83.011901 http://doi.org/10.1103/PhysRevC.83.011901 A New method for the experimental study of topological effects in the quark-gluon plasma B. Abelev , et al ., Charge separation relative to the reaction plane in Pb-Pb collisions at <math id="M152"><mrow><msqrt><mrow><msub><mi>s</mi><mrow><mi>N</mi><mi>N</mi></mrow></msub></mrow></msqrt><mo>=</mo><mn>2.76</mn></mrow></math> TeV . Phys. Rev. Lett . 110 , 012301 ( 2013 ). doi: 10.1103/PhysRevLett.110.012301 http://doi.org/10.1103/PhysRevLett.110.012301 Charge separation relative to the reaction plane in Pb-Pb collisions at sNN=2.76 TeV G. Wang , Search for chiral magnetic effects in high-energy nuclear collisions . Proceedings, 23rd International Conference on Ultrarelativistic Nucleus-Nucleus Collisions: Quark Matter 2012 (QM 2012) : Washington, DC, USA , August 13-18, 2012 . Search for chiral magnetic effects in high-energy nuclear collisions G. Wang , Search for chiral magnetic effects in high-energy nuclear collisions Nucl. Phys. A 904-905 , 248c - 255c ( 2013 ). doi: 10.1016/j.nuclphysa.2013.01.069 http://doi.org/10.1016/j.nuclphysa.2013.01.069 Search for chiral magnetic effects in high-energy nuclear collisions L. Adamczyk , et al ., 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 S. Pratt , Alternative contributions to the angular correlations observed at RHIC associated with parity fluctuations . ( 2010 ). arXiv: 1002.1758 Alternative contributions to the angular correlations observed at RHIC associated with parity fluctuations S. Pratt , S. Schlichting , S. Gavin , Effects of momentum conservation and flow on angular correlations at RHIC . Phys. Rev. C 84 , 024909 ( 2011 ). doi: 10.1103/PhysRevC.84.024909 http://doi.org/10.1103/PhysRevC.84.024909 Effects of momentum conservation and flow on angular correlations at RHIC A. Bzdak , V. Koch , J.F. Liao , Azimuthal correlations from transverse momentum conservation and possible local parity violation . Phys. Rev. C 83 , 014905 ( 2011 ). doi: 10.1103/PhysRevC.83.014905 http://doi.org/10.1103/PhysRevC.83.014905 Azimuthal correlations from transverse momentum conservation and possible local parity violation S. Schlichting , S. Pratt , Effects of charge conservation and flow on fluctuations of parity-odd observables ar RHIC . ( 2010 ) arXiv: 1005.5341 Effects of charge conservation and flow on fluctuations of parity-odd observables ar RHIC S. Schlichting , S. Pratt , Charge conservation at energies available at the BNL Relativistic Heavy Ion Collider and contributions to local parity violation observables . Phys. Rev. C 83 , 014913 ( 2011 ). doi: 10.1103/PhysRevC.83.014913 http://doi.org/10.1103/PhysRevC.83.014913 Charge conservation at energies available at the BNL Relativistic Heavy Ion Collider and contributions to local parity violation observables F. Wen , L. Wen , G. Wang , Procedure for removing flow backgrounds from the charge-separation observable perpendicular to the reaction plane in heavy-ion collisions . ( 2016 ). arXivi: 1608.03205 Procedure for removing flow backgrounds from the charge-separation observable perpendicular to the reaction plane in heavy-ion collisions F. Wang , J. Zhao , Challenges in flow background removal in search for the chiral magnetic effect . ( 2016 ). arXivi: 1608.06610 Challenges in flow background removal in search for the chiral magnetic effect S.A. Voloshin , Testing the chiral magnetic effect with central U+U collisions . Phys. Rev. Lett . 105 , 172301 ( 2010 ). doi: 10.1103/PhysRevLett.105.172301 http://doi.org/10.1103/PhysRevLett.105.172301 Testing the chiral magnetic effect with central U+U collisions L. Adamczyk , et al ., Azimuthal anisotropy in U+U and Au+Au collisions at RHIC . Phys. Rev. Lett . 115 , 222301 ( 2015 ). doi: 10.1103/PhysRevLett.115.222301 http://doi.org/10.1103/PhysRevLett.115.222301 Azimuthal anisotropy in U+U and Au+Au collisions at RHIC Q.Y. Shou , harge asymmetry dependency of π/K anisotropic flow in U+U <math id="M153"><mrow><msqrt><mrow><msub><mi>s</mi><mrow><mi>N</mi><mi>N</mi></mrow></msub></mrow></msqrt><mo>=</mo></mrow></math> 193 GeV and Au+Au <math id="M154"><mrow><msqrt><mrow><msub><mi>s</mi><mrow><mi>N</mi><mi>N</mi></mrow></msub></mrow></msqrt><mo>=</mo></mrow></math> 200 GeV collisions at STAR . J. Phys. Conf. Ser . 509 , 012033 ( 2014 ). doi: 10.1088/1742-6596/509/1/012033 http://doi.org/10.1088/1742-6596/509/1/012033 harge asymmetry dependency of π/K anisotropic flow in U+U sNN= 193 GeV and Au+Au sNN= 200 GeV collisions at STAR W.T. Deng , X.G. Huang , G.L. Ma , et al ., Test the chiral magnetic effect with isobaric collisions . ( 2016 ) arXiv: 1607.04697 Test the chiral magnetic effect with isobaric collisions V. Skokov , P. Sorensen , V. Koch , et al ., Chiral magnetic effect task force report . ( 2016 ) arXiv: 1608.00982 Chiral magnetic effect task force report V. Khachatryan , Observation of charge-dependent azimuthal correlations in pPb collisions and its implication for the search for the chiral magnetic effect . ( 2016 ). arXivi: 1610.00263 Observation of charge-dependent azimuthal correlations in pPb collisions and its implication for the search for the chiral magnetic effect G.M. Newman , D.T. Son , Response of strongly-interacting matter to magnetic field: Some exact results . Phys. Rev. D 73 , 045006 ( 2006 ). doi: 10.1103/PhysRevD.73.045006 http://doi.org/10.1103/PhysRevD.73.045006 Response of strongly-interacting matter to magnetic field: Some exact results O. Bergman , G. Lifschytz , M. Lippert , Magnetic properties of dense holographic QCD . Phys. Rev. D 79 , 105024 ( 2009 ). doi: 10.1103/PhysRevD.79.105024 http://doi.org/10.1103/PhysRevD.79.105024 Magnetic properties of dense holographic QCD E.V. Gorbar , V.A. Miransky , I.A. Shovkovy , Chiral asymmetry of the Fermi surface in dense relativistic matter in a magnetic field . Phys. Rev. C 80 , 032801 ( 2009 ). doi: 10.1103/PhysRevC.80.032801 http://doi.org/10.1103/PhysRevC.80.032801 Chiral asymmetry of the Fermi surface in dense relativistic matter in a magnetic field G. Basar , G.V. Dunne , D.E. Kharzeev , Chiral magnetic spiral . Phys. Rev. Lett . 104 , 232301 ( 2010 ). doi: 10.1103/PhysRevLett.104.232301 http://doi.org/10.1103/PhysRevLett.104.232301 Chiral magnetic spiral E.V. Gorbar , V.A. Miransky , I.A. Shovkovy , Chiral asymmetry and axial anomaly in magnetized relativistic matter . Phys. Lett. B 695 , 354 - 358 ( 2011 ). doi: 10.1016/j.physletb.2010.11.022 http://doi.org/10.1016/j.physletb.2010.11.022 Chiral asymmetry and axial anomaly in magnetized relativistic matter D.K. Hong , Anomalous currents in dense matter under a magneticield . Phys. Lett. B 699 , 305 - 308 ( 2011 ). doi: 10.1016/j.physletb.2011.04.010 http://doi.org/10.1016/j.physletb.2011.04.010 Anomalous currents in dense matter under a magneticield K. Landsteiner , E. Megias , F. Pena-Benitez , Gravitational anomaly and transport . Phys. Rev. Lett . 107 , 021601 ( 2011 ). doi: 10.1103/PhysRevLett.107.021601 http://doi.org/10.1103/PhysRevLett.107.021601 Gravitational anomaly and transport E. V. Gorbar , V. A. Miransky , I.A. Shovkovy , et al ., Radiative corrections to chiral separation effect in QED . Phys. Rev. D 88 , 025025 ( 2013 ). doi: 10.1103/PhysRevD.88.025025 http://doi.org/10.1103/PhysRevD.88.025025 Radiative corrections to chiral separation effect in QED N. Yamamoto , Generalized Bloch theorem and chiral transport phenomena . Phys. Rev. D 92 , 085011 ( 2015 ). doi: 10.1103/PhysRevD.92.085011 http://doi.org/10.1103/PhysRevD.92.085011 Generalized Bloch theorem and chiral transport phenomena X.G. Huang , Simulating chiral magnetic and separation effects with spin-orbit coupled atomic gases . Sci. Rep . 6 , 20601 ( 2016 ). doi: 10.1038/srep20601 http://doi.org/10.1038/srep20601 Simulating chiral magnetic and separation effects with spin-orbit coupled atomic gases M. Stephanov , H.-U. Yee , Y. Yin , Collective modes of chiral kinetic theory in a magnetic field . Phys. Rev. D 91 , 125014 ( 2015 ). doi: 10.1103/PhysRevD.91.125014 http://doi.org/10.1103/PhysRevD.91.125014 Collective modes of chiral kinetic theory in a magnetic field Y. Burnier , D.E. Kharzeev , J. Liao , et al ., From the chiral magnetic wave to the charge dependence of elliptic flow . ( 2012 ). arXiv: 1208.2537 From the chiral magnetic wave to the charge dependence of elliptic flow S.F. Taghavi , U.A. Wiedemann , Chiral magnetic wave in an expanding QCD fluid . Phys. Rev. C 91 , 024902 ( 2015 ). doi: 10.1103/PhysRevC.91.024902 http://doi.org/10.1103/PhysRevC.91.024902 Chiral magnetic wave in an expanding QCD fluid H.-U. Yee , Y. Yin , Realistic implementation of chiral magnetic wave in heavy ion collisions . Phys. Rev. C 89 , 044909 ( 2014 ). doi: 10.1103/PhysRevC.89.044909 http://doi.org/10.1103/PhysRevC.89.044909 Realistic implementation of chiral magnetic wave in heavy ion collisions H. Ke , Charge asymmetry dependency of π+/π- elliptic flow in Au + Au collisions at <math id="M155"><mrow><msqrt><mrow><msub><mi>s</mi><mrow><mi>N</mi><mi>N</mi></mrow></msub></mrow></msqrt></mrow></math> =200 GeV . J. Phys. Conf. Ser . 389 , 012035 ( 2012 ). doi: 10.1088/1742-6596/389/1/012035 http://doi.org/10.1088/1742-6596/389/1/012035 Charge asymmetry dependency of π+/π- elliptic flow in Au + Au collisions at sNN =200 GeV L. Adamczyk , et al ., 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 J. Adam , et al ., Charge-dependent flow and the search for the chiral magnetic wave in Pb-Pb collisions at <math id="M156"><mrow><msqrt><mrow><msub><mi>s</mi><mrow><mtext>N</mtext><mi>N</mi></mrow></msub></mrow></msqrt><mo>=</mo></mrow></math> 2.76 TeV . Phys. Rev. C 93 , 044903 ( 2016 ). doi: 10.1103/PhysRevC.93.044903 http://doi.org/10.1103/PhysRevC.93.044903 Charge-dependent flow and the search for the chiral magnetic wave in Pb-Pb collisions at sNN= 2.76 TeV M. Stephanov , H.-U. Yee , Charged elliptic flow at zero charge asymmetry . Phys. Rev. C 88 , 014908 ( 2013 ). doi: 10.1103/PhysRevC.88.014908 http://doi.org/10.1103/PhysRevC.88.014908 Charged elliptic flow at zero charge asymmetry J.C. Dunlop , M.A. Lisa , P. Sorensen , Constituent quark scaling violation due to baryon number transport . Phys. Rev. C 84 , 044914 ( 2011 ). doi: 10.1103/PhysRevC.84.044914 http://doi.org/10.1103/PhysRevC.84.044914 Constituent quark scaling violation due to baryon number transport J. Xu , L.W. Chen , C.M. Ko , et al ., Effects of hadronic potentials on elliptic flows in relativistic heavy ion collisions . Phys. Rev. C 85 , 041901 ( 2012 ). doi: 10.1103/PhysRevC.85.041901 http://doi.org/10.1103/PhysRevC.85.041901 Effects of hadronic potentials on elliptic flows in relativistic heavy ion collisions C.M. Ko , T. Song , F. Li , et al ., Partonic mean-field effects on matter and antimatter elliptic flows . Nucl. Phys. A 928 , 234 - 246 ( 2014 ). doi: 10.1016/j.nuclphysa.2014.05.016 http://doi.org/10.1016/j.nuclphysa.2014.05.016 Partonic mean-field effects on matter and antimatter elliptic flows A. Bzdak , P. Bozek , Contributions to the event-by-event charge asymmetry dependence for the elliptic flow of pi+ and pi- in heavy-ion collisions . Phys. Lett. B 726 , 239 - 243 ( 2013 ). doi: 10.1016/j.physletb.2013.08.003 http://doi.org/10.1016/j.physletb.2013.08.003 Contributions to the event-by-event charge asymmetry dependence for the elliptic flow of pi+ and pi- in heavy-ion collisions Y. Hatta , A. Monnai , B.W. Xiao , Elliptic flow difference of charged pions in heavy-ion collisions . Nucl. Phys. A 947 , 155 - 160 ( 2016 ). doi: 10.1016/j.nuclphysa.2015.12.009 http://doi.org/10.1016/j.nuclphysa.2015.12.009 Elliptic flow difference of charged pions in heavy-ion collisions X.G. Huang , J.F. Liao , Axial current generation from electric field: Chiral electric separation effect . Phys. Rev. Lett . 110 , 232302 ( 2013 ). doi: 10.1103/PhysRevLett.110.232302 http://doi.org/10.1103/PhysRevLett.110.232302 Axial current generation from electric field: Chiral electric separation effect Y. Jiang , X.G. Huang , J.F. Liao , Chiral electric separation effect in the quark-gluon plasma . Phys. Rev. D 91 , 045001 ( 2015 ). doi: 10.1103/PhysRevD.91.045001 http://doi.org/10.1103/PhysRevD.91.045001 Chiral electric separation effect in the quark-gluon plasma S. Pu , S.Y. Wu , D.L. Yang , Holographic chiral electric separation effect . Phys. Rev. D 89 , 085024 ( 2014 ). doi: 10.1103/PhysRevD.89.085024 http://doi.org/10.1103/PhysRevD.89.085024 Holographic chiral electric separation effect S. Pu , S.Y. Wu , D.L. Yang , Chiral hall effect and chiral electric waves . Phys. Rev. D 91 , 025011 ( 2015 ). doi: 10.1103/PhysRevD.91.025011 http://doi.org/10.1103/PhysRevD.91.025011 Chiral hall effect and chiral electric waves 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 D. Kharzeev , A. Zhitnitsky , Charge separation induced by P-odd bubbles in QCD matter . Nucl. Phys. A 797 , 67 - 79 ( 2007 ). doi: 10.1016/j.nuclphysa.2007.10.001 http://doi.org/10.1016/j.nuclphysa.2007.10.001 Charge separation induced by P-odd bubbles in QCD matter J. Erdmenger , M. Haack , M. Kaminski , et al ., Fluid dynamics of R-charged black holes . JHEP 01 , 055 ( 2009 ). doi: 10.1088/1126-6708/2009/01/055 http://doi.org/10.1088/1126-6708/2009/01/055 Fluid dynamics of R-charged black holes N. Banerjee , J. Bhattacharya , S. Bhattacharyya , et al ., Hydrodynamics from charged black branes . JHEP 01 , 094 ( 2011 ). doi: 10.1007/JHEP01(2011)094 http://doi.org/10.1007/JHEP01(2011)094 Hydrodynamics from charged black branes D.T. Son , P. Surowka , Hydrodynamics with triangle anomalies . Phys. Rev. Lett . 103 , 191601 ( 2009 ). doi: 10.1103/PhysRevLett.103.191601 http://doi.org/10.1103/PhysRevLett.103.191601 Hydrodynamics with triangle anomalies Y. Jiang , X.G. Huang , J.F. Liao , Chiral vortical wave and induced flavor charge transport in a rotating quark-gluon plasma . Phys. Rev. D 92 , 071501 ( 2015 ). doi: 10.1103/PhysRevD.92.07150 http://doi.org/10.1103/PhysRevD.92.07150 Chiral vortical wave and induced flavor charge transport in a rotating quark-gluon plasma L.P. Csernai , V.K. Magas , D.J. Wang , Flow vorticity in peripheral high energy heavy ion collisions . Phys. Rev. C 87 , 034906 ( 2013 ). doi: 10.1103/PhysRevC.87.034906 http://doi.org/10.1103/PhysRevC.87.034906 Flow vorticity in peripheral high energy heavy ion collisions F. Becattini , G. Inghirami , V. Rolando , et al ., A study of vorticity formation in high energy nuclear collisions . Eur. Phys. J. C 75 , 406 ( 2015 ). doi: 10.1140/epjc/s10052-015-3624-1 http://doi.org/10.1140/epjc/s10052-015-3624-1 A study of vorticity formation in high energy nuclear collisions Y. Jiang , Z.W. Lin , J.F. Liao , Rotating quark-gluon plasma in relativistic heavy ion collisions . ( 2016 ). arXiv: 1602.06580 Rotating quark-gluon plasma in relativistic heavy ion collisions W.T. Deng , X.G. Huang , Vorticity in heavy-ion collisions . Phys. Rev. C 93 , 064907 ( 2016 ). doi: 10.1103/PhysRevC.93.06490 http://doi.org/10.1103/PhysRevC.93.06490 Vorticity in heavy-ion collisions D.E. Kharzeev , D.T. Son , Testing the chiral magnetic and chiral vortical effects in heavy ion collisions . Phys. Rev. Lett . 106 , 062301 ( 2011 ). doi: 10.1103/PhysRevLett.106.062301 http://doi.org/10.1103/PhysRevLett.106.062301 Testing the chiral magnetic and chiral vortical effects in heavy ion collisions F. Zhao , Λ(KS0)?Ch± and Λ-p azimuthal correlations with respect to event plane and search for chiral magnetic and vortical effects . Nucl. Phys. A 931 , 746 - 751 ( 2014 ). doi: 10.1016/j.nuclphysa.2014.08.108 http://doi.org/10.1016/j.nuclphysa.2014.08.108 Λ(KS0)?Ch± and Λ-p azimuthal correlations with respect to event plane and search for chiral magnetic and vortical effects M. D’Elia , F. Negro , Chiral properties of strong interactions in a magnetic background . Phys. Rev.D 83 , 114028 ( 2011 ). doi: 10.1103/PhysRevD.83.114028 http://doi.org/10.1103/PhysRevD.83.114028 Chiral properties of strong interactions in a magnetic background G.S. Bali , F. Bruckmann , G. Endrodi , et al ., Magnetic field-induced gluonic (inverse) catalysis and pressure (an)isotropy in QCD . JHEP 04 , 130 ( 2013 ). doi: 10.1007/JHEP04(2013)130 http://doi.org/10.1007/JHEP04(2013)130 Magnetic field-induced gluonic (inverse) catalysis and pressure (an)isotropy in QCD A. Andronic , et al ., Heavy-flavour and quarkonium production in the LHC era: from proton?proton to heavy-ion collisions . Eur. Phys. J. C 76 , 107 ( 2016 ). doi: 10.1140/epjc/s10052-015-3819-5 http://doi.org/10.1140/epjc/s10052-015-3819-5 Heavy-flavour and quarkonium production in the LHC era: from proton?proton to heavy-ion collisions M.A. Andreichikov , V.D. Orlovsky , Yu. A. Simonov , Asymptotic freedom in strong magnetic fields . Phys. Rev. Lett . 110 , 162002 ( 2013 ). doi: 10.1103/PhysRevLett.110.162002 http://doi.org/10.1103/PhysRevLett.110.162002 Asymptotic freedom in strong magnetic fields M.A. Andreichikov , B. O. Kerbikov , V.D. Orlovsky , et al ., Meson spectrum in strong magnetic fields . Phys. Rev. D 87 , 094029 ( 2013 ). doi: 10.1103/PhysRevD.87.094029 http://doi.org/10.1103/PhysRevD.87.094029 Meson spectrum in strong magnetic fields B.B. Brandt , G. Bali , G. Endrodi , et al ., QCD spectroscopy and quark mass renormalisation in external magnetic fields with Wilson fermions . PoS , 265 ( 2016 ). arXiv: 1510.03899 doi: 10.1103/PhysRevD.87.094029 http://doi.org/10.1103/PhysRevD.87.094029 QCD spectroscopy and quark mass renormalisation in external magnetic fields with Wilson fermions E.V. Luschevskaya , O.E. Solovjeva , O.A. Kochetkov , et al ., Magnetic polarizabilities of light mesons in SU(3) lattice gauge theory . Nucl. Phys. B 898 , 627 - 643 ( 2015 ). doi: 10.1016/j.nuclphysb.2015.07.023 http://doi.org/10.1016/j.nuclphysb.2015.07.023 Magnetic polarizabilities of light mesons in SU(3) lattice gauge theory H. Taya , Hadron masses in strong magnetic fields . Phys. Rev. D 92 , 014038 ( 2015 ). doi: 10.1103/PhysRevD.92.014038 http://doi.org/10.1103/PhysRevD.92.014038 Hadron masses in strong magnetic fields N. Brambilla , P. Pietrulewicz , A. Vairo , Model-independent study of electric dipole transitions in quarkonium . Phys. Rev. D 85 , 094005 ( 2012 ). doi: 10.1103/PhysRevD.85.094005 http://doi.org/10.1103/PhysRevD.85.094005 Model-independent study of electric dipole transitions in quarkonium A. Pineda , J. Segovia , Improved determination of heavy quarkonium magnetic dipole transitions in potential nonrelativistic QCD . Phys. Rev. D 87 , 074024 ( 2013 ). doi: 10.1103/PhysRevD.87.074024 http://doi.org/10.1103/PhysRevD.87.074024 Improved determination of heavy quarkonium magnetic dipole transitions in potential nonrelativistic QCD E. Eichten , K. Gottfried , T. Kinoshita , et al ., Charmonium: The model . Phys. Rev. D 17 , 3090 ( 1978 ). doi: 10.1103/PhysRevD.17.3090 http://doi.org/10.1103/PhysRevD.17.3090 Charmonium: The model E. Eichten , K. Gottfried , T. Kinoshita , et al ., Charmonium: Comparison with experiment . Phys. Rev. D 21 , 203 ( 1980 ). doi: 10.1103/PhysRevD.21.203 http://doi.org/10.1103/PhysRevD.21.203 Charmonium: Comparison with experiment J.S. Schwinger , On gauge invariance and vacuum polarization . Phys. Rev . 82 , 664 - 679 ( 1951 ). doi: 10.1103/PhysRev.82.664 http://doi.org/10.1103/PhysRev.82.664 On gauge invariance and vacuum polarization A.V. Sadofyev , Y. Yin , The charmonium dissociation in an "ganomalous wind" . JHEP 01 , 052 ( 2016 ). doi: 10.1007/JHEP01(2016)052 http://doi.org/10.1007/JHEP01(2016)052 The charmonium dissociation in an "ganomalous wind" M.A. Shifman , A.I. Vainshtein , V.I. Zakharov , QCD and resonance physics. theoretical foundations . Nucl. Phys. B 147 , 385 - 447 ( 1979 ). doi: 10.1016/0550-3213(79)90022-1 http://doi.org/10.1016/0550-3213(79)90022-1 QCD and resonance physics. theoretical foundations M.A. Shifman , A.I. Vainshtein , V.I. Zakharov , QCD and resonance physics: Applications . Nucl. Phys. B 147 , 448 - 518 ( 1979 ). doi: 10.1016/0550-3213(79)90023-3 http://doi.org/10.1016/0550-3213(79)90023-3 QCD and resonance physics: Applications L.J. Reinders , H. Rubinstein , S. Yazaki , Hadron properties from QCD sum rules . Phys. Rept . 127 , 1 ( 1985 ). doi: 10.1016/0370-1573(85)90065-1 http://doi.org/10.1016/0370-1573(85)90065-1 Hadron properties from QCD sum rules M.A. Shifman , Snapshots of hadrons or the story of how the vacuum which are produced, live and die in the QCD vacuum . Prog. Theor. Phys. Suppl . 131 , 1 - 71 ( 1998 ). doi: 10.1143/PTPS.131.1 http://doi.org/10.1143/PTPS.131.1 Snapshots of hadrons or the story of how the vacuum which are produced, live and die in the QCD vacuum V.A. Novikov , M.A. Shifman , A.I. Vainshtein , et al ., Calculations in external fields in quantum chromodynamics . Technical review. Fortsch. Phys . 32 , 585 ( 1984 ). Calculations in external fields in quantum chromodynamics S. Ozaki , K. Itakura , Y. Kuramoto , Magnetically induced QCD kondo effect . ( 2015 ). arXiv: 1509.06966 Magnetically induced QCD kondo effect K. Hattori , K. Itakura , S. Ozaki , et al ., QCD Kondo effect: quark matter with heavy-flavor impurities . Phys. Rev. D 92 , 065003 ( 2015 ). doi: 10.1103/PhysRevD.92.065003 http://doi.org/10.1103/PhysRevD.92.065003 QCD Kondo effect: quark matter with heavy-flavor impurities S. Yasui , K. Suzuki , K. Itakura , Kondo phase diagram of quark matter . ( 2016 ). arXiv: 1604.07208 Kondo phase diagram of quark matter D.K. Hong , An effective field theory of QCD at high density . Phys. Lett. B 473 , 118 - 125 ( 2000 ). doi: 10.1016/S0370-2693(99)01472-0 http://doi.org/10.1016/S0370-2693(99)01472-0 An effective field theory of QCD at high density D.K. Hong , Aspects of high density effective theory in QCD . Nucl. Phys. B 582 , 451 - 476 ( 2000 ). doi: 10.1016/S0550-3213(00)00330-8 http://doi.org/10.1016/S0550-3213(00)00330-8 Aspects of high density effective theory in QCD T. Schefer , Hard loops, soft loops, and high density effective field theory . Nucl. Phys. A 728 , 251 - 271 ( 2003 ). doi: 10.1016/j.nuclphysa.2003.08.028 http://doi.org/10.1016/j.nuclphysa.2003.08.028 Hard loops, soft loops, and high density effective field theory D.K. Hong , Y. Kim , S.-J. Sin , RG analysis of magnetic catalysis in dynamical symmetry breaking . Phys. Rev. D 54 , 7879 - 7883 ( 1996 ). doi: 10.1103/PhysRevD.54.7879 http://doi.org/10.1103/PhysRevD.54.7879 RG analysis of magnetic catalysis in dynamical symmetry breaking D.K. Hong , Magnetic catalysis in quantum electrodynamics . Phys. Rev. D 57 , 3759 - 3762 ( 1998 ). doi: 10.1103/PhysRevD.57.3759 http://doi.org/10.1103/PhysRevD.57.3759 Magnetic catalysis in quantum electrodynamics G.D. Moore , D. Teaney , How much do heavy quarks thermalize in a heavy ion collision? Phys . Rev. C 71 , 064904 ( 2005 ). doi: 10.1103/PhysRevC.71.064904 http://doi.org/10.1103/PhysRevC.71.064904 How much do heavy quarks thermalize in a heavy ion collision? Phys L. Dolan , R. Jackiw , Symmetry behavior at finite temperature . Phys. Rev. D 9 , 3320 - 3341 ( 1974 ). doi: 10.1103/PhysRevD.9.3320 http://doi.org/10.1103/PhysRevD.9.3320 Symmetry behavior at finite temperature R. Baier , E. Pilon , On the axial anomaly at finite temperature in the Schwinger model . Z. Phys. C 52 , 339 - 342 ( 1991 ). doi: 10.1007/BF01560452 http://doi.org/10.1007/BF01560452 On the axial anomaly at finite temperature in the Schwinger model K. Fukushima , Magnetic-field induced screening effect and collective excitations . Phys. Rev. D 83 , 11150 ( 2011 ). doi: 10.1103/PhysRevD.83.111501 http://doi.org/10.1103/PhysRevD.83.111501 Magnetic-field induced screening effect and collective excitations K. Hattori , K. Itakura , Vacuum birefringence in strong magnetic fields: (I) Photon polarization tensor with all the Landau levels . Annals Phys . 330 , 23 - 54 ( 2013 ). doi: 10.1016/j.aop.2012.11.010 http://doi.org/10.1016/j.aop.2012.11.010 Vacuum birefringence in strong magnetic fields: (I) Photon polarization tensor with all the Landau levels K. Hattori , K. Itakura , Vacuum birefringence in strong magnetic fields: (II)Complex refractive index from the lowest Landau level . Annals Phys . 334 , 58 - 82 ( 2013 ). doi: 10.1016/j.aop.2013.03.016 http://doi.org/10.1016/j.aop.2013.03.016 Vacuum birefringence in strong magnetic fields: (II)Complex refractive index from the lowest Landau level A.V. Sadofyev , Y. Yin , Drag suppression in anomalous chiral media . Phys. Rev. D 93 , 125026 ( 2016 ). doi: 10.1103/PhysRevD.93.125026 http://doi.org/10.1103/PhysRevD.93.125026 Drag suppression in anomalous chiral media S. Caron-Huot , G.D. Moore , Heavy quark diffusion in perturbative QCD at next-to-leading order . Phys. Rev. Lett . 100 , 052301 ( 2008 ). doi: 10.1103/PhysRevLett.100.052301 http://doi.org/10.1103/PhysRevLett.100.052301 Heavy quark diffusion in perturbative QCD at next-to-leading order K. Rajagopal , A.V. Sadofyev , Chiral drag force . JHEP 10 , 018 ( 2015 ). doi: 10.1007/JHEP10(2015)018 http://doi.org/10.1007/JHEP10(2015)018 Chiral drag force M.A. Stephanov , H.-U. Yee , No-drag frame for anomalous chiral fluid . Phys. Rev. Lett . 116 , 122302 ( 2016 ). doi: 10.1103/PhysRevLett.116.122302 http://doi.org/10.1103/PhysRevLett.116.122302 No-drag frame for anomalous chiral fluid S.Y. Li , K.A. Mamo , H.-U. Yee , Jet quenching parameter of the quark-gluon plasma in a strong magnetic field: Perturbative QCD and AdS/CFT correspondence . Phys. Rev.D 94 085016 , ( 2016 ). doi: 10.1103/PhysRevD.94.085016 http://doi.org/10.1103/PhysRevD.94.085016 Jet quenching parameter of the quark-gluon plasma in a strong magnetic field: Perturbative QCD and AdS/CFT correspondence K. Hattori , S.Y. Li , D. Satow , et al ., Longitudinal conductivity in strong magnetic field in perturbative QCD: Complete Leading Order . 2016 , arXiv:1610.06839 Longitudinal conductivity in strong magnetic field in perturbative QCD: Complete Leading Order K. Hattori , and D. Satow , Electrical conductivity of Quark-Gluon plasma in strong magnetic fields . 2016 , arXiv:1610.06818 Electrical conductivity of Quark-Gluon plasma in strong magnetic fields Y.-S. Oh , S. Kim , S.H. Lee , Quarkonium hadron interactions in QCD . Phys. Rev. C 65 , 067901 ( 2002 ). doi: 10.1103/PhysRevC.65.067901 http://doi.org/10.1103/PhysRevC.65.067901 Quarkonium hadron interactions in QCD T. Song , S.H. Lee , Quarkonium-hadron interactions in perturbative QCD . Phys. Rev. D 72 , 034002 ( 2005 ). doi: 10.1103/PhysRevD.72.034002 http://doi.org/10.1103/PhysRevD.72.034002 Quarkonium-hadron interactions in perturbative QCD K.H. thanks Kei Suzuki for useful comments on the manuscript. 10.1007/s41365-016-0178-3 This work was supported by Shanghai Natural Science Foundation (No. 14ZR1403000), 1000 Young Talents Program of China, and the National Natural Science Foundation of China (No. 11535012). K.H. is also supported by China Postdoctoral Science Foundation under Grant No. 2016M590312 and is grateful to support from RIKEN-BNL Research Center. 2016-10-14 2016-11-16 2016-11-26 2017-01-11 2017-02-01 2021-04-11 © 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-016-0178-3 http://dx.doi.org/10.1007/s41365-016-0178-3 . 2017 Appendix: Mixing strengths from the Bethe-Salpeter amplitudes 2 It would be instructive to see the computation of the coupling strength at the three-point vertex among two quarkonia and an external magnetic field by using the Bethe-Salpeter amplitudes obtained in the ladder approximation and the heavy-quark limit [ 240 240 , 241 241 ]. This computation involves a typical technique for the perturbation theory in the presence of a magnetic field. For the static charmonium carrying a momentum q =(2 m - ε 0 , 0,0,0) with ≤ b being the binding energy, the Bethe-Salpeter amplitudes for η c and J/ψ are, respectively, given by <math display="block" id="M157"><mrow><msup><mi>Γ</mi><mn>5</mn></msup><mo stretchy="false">(</mo><mi>p</mi><mn>,</mn><mi>p</mi><mo>−</mo><mi>q</mi><mo stretchy="false">)</mo><mo>=</mo><mrow><mo>(</mo><mrow><msub><mi>ϵ</mi><mn>0</mn></msub><mo>+</mo><mfrac><mrow><msup><mi>p</mi><mn>2</mn></msup></mrow><mi>m</mi></mfrac></mrow><mo>)</mo></mrow><mtext>​</mtext><msqrt><mrow><mfrac><mrow><msub><mi>m</mi><mrow><mi>c</mi><mover accent="true"><mi>c</mi><mo>¯</mo></mover></mrow></msub></mrow><mrow><msub><mi>N</mi><mi>c</mi></msub></mrow></mfrac></mrow></msqrt><mtext> </mtext><msub><mi>ψ</mi><mrow><mn>1</mn><mi>S</mi></mrow></msub><mo stretchy="false">(</mo><mi>p</mi><mo stretchy="false">)</mo><mtext> </mtext><msub><mi>P</mi><mo>+</mo></msub><msup><mi>γ</mi><mn>5</mn></msup><msub><mi>P</mi><mo>−</mo></msub><mn>,</mn></mrow></math> <math display="block" id="M158"><mrow><msup><mi>Γ</mi><mi>μ</mi></msup><mo stretchy="false">(</mo><mi>p</mi><mn>,</mn><mi>p</mi><mo>−</mo><mi>q</mi><mo stretchy="false">)</mo><mo>=</mo><mrow><mo>(</mo><mrow><msub><mi>ϵ</mi><mn>0</mn></msub><mo>+</mo><mfrac><mrow><msup><mi>p</mi><mn>2</mn></msup></mrow><mi>m</mi></mfrac></mrow><mo>)</mo></mrow><mtext>​</mtext><msqrt><mrow><mfrac><mrow><msub><mi>m</mi><mrow><mi>c</mi><mover accent="true"><mi>c</mi><mo>¯</mo></mover></mrow></msub></mrow><mrow><msub><mi>N</mi><mi>c</mi></msub></mrow></mfrac></mrow></msqrt><mtext> </mtext><msub><mi>ψ</mi><mrow><mn>1</mn><mi>S</mi></mrow></msub><mo stretchy="false">(</mo><mi>p</mi><mo stretchy="false">)</mo><mtext> </mtext><msub><mi>P</mi><mo>+</mo></msub><msup><mi>γ</mi><mi>μ</mi></msup><msub><mi>P</mi><mo>−</mo></msub><mn>,</mn></mrow></math> where we have the projection operators <math id="M159"><mrow><msub><mi>P</mi><mo>±</mo></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo stretchy="false">(</mo><mn>1</mn><mo>±</mo><msup><mi>γ</mi><mn>0</mn></msup><mo stretchy="false">)</mo></mrow></math> and the ground-state wave function of the S-wave bound state ψ 1 S ( p ). The mass <math id="M160"><mrow><msub><mi>m</mi><mrow><mi>c</mi><mover accent="true"><mi>c</mi><mo>¯</mo></mover></mrow></msub></mrow></math> is those of η c and J/ψ , which are degenerated in the heavy-quark limit. The number of the color is N c =3. We show a calculation of a coupling strength in the mixing between η c and the longitudinal J/ψ from triangle diagrams ( Fig. 21 Fig. 21 ). Interactions between quarks and external magnetic fields are taken into account by employing the Fock-Schwinger gauge. In this gauge, the quark propagators with one and two insertions of constant external fields are expressed as [ 216 216 ] <math display="block" id="M161"><mtable columnalign="left"><mtr><mtd><msub><mi>S</mi><mn>1</mn></msub><mo stretchy="false">(</mo><mi>p</mi><mo stretchy="false">)</mo><mo>=</mo><mo>−</mo><mfrac><mi>i</mi><mn>4</mn></mfrac><msub><mi>Q</mi><mrow><mtext>e</mtext><mi>m</mi></mrow></msub><msub><mi>F</mi><mrow><mi>α</mi><mi>β</mi></mrow></msub><mfrac><mn>1</mn><mrow><msup><mrow><mo stretchy="false">(</mo><msup><mi>p</mi><mn>2</mn></msup><mo>−</mo><msup><mi>m</mi><mn>2</mn></msup><mo>+</mo><mi>i</mi><mi>ε</mi><mo stretchy="false">)</mo></mrow><mn>2</mn></msup></mrow></mfrac></mtd></mtr><mtr><mtd><mtext> </mtext><mo>×</mo><mrow><mo>{</mo><mrow><msup><mi>σ</mi><mrow><mi>α</mi><mi>β</mi></mrow></msup><mo stretchy="false">(</mo><mi>p</mi><mo>+</mo><mi>m</mi><mo stretchy="false">)</mo><mo>+</mo><mo stretchy="false">(</mo><mi>p</mi><mo>+</mo><mi>m</mi><mo stretchy="false">)</mo><msup><mi>σ</mi><mrow><mi>α</mi><mi>β</mi></mrow></msup></mrow><mo>}</mo></mrow><mtext> </mtext><mn>,</mn></mtd></mtr></mtable></math> <math display="block" id="M162"><mtable columnalign="left"><mtr><mtd><msub><mi>S</mi><mn>2</mn></msub><mo stretchy="false">(</mo><mi>p</mi><mo stretchy="false">)</mo><mo>=</mo><mo>−</mo><mfrac><mn>1</mn><mn>4</mn></mfrac><msubsup><mi>Q</mi><mrow><mtext>e</mtext><mi>m</mi></mrow><mn>2</mn></msubsup><msub><mi>F</mi><mrow><mi>α</mi><mi>β</mi></mrow></msub><msub><mi>F</mi><mrow><mi>μ</mi><mi>ν</mi></mrow></msub><mfrac><mn>1</mn><mrow><msup><mrow><mo stretchy="false">(</mo><msup><mi>p</mi><mn>2</mn></msup><mo>−</mo><msup><mi>m</mi><mn>2</mn></msup><mo>+</mo><mi>i</mi><mi>ε</mi><mo stretchy="false">)</mo></mrow><mn>5</mn></msup></mrow></mfrac></mtd></mtr><mtr><mtd><mtext> </mtext><mo>×</mo><mo stretchy="false">(</mo><mi>p</mi><mo>+</mo><mi>m</mi><mo stretchy="false">)</mo><mrow><mo>{</mo><mrow><msup><mi>f</mi><mrow><mi>α</mi><mi>β</mi><mi>μ</mi><mi>ν</mi></mrow></msup><mo>+</mo><msup><mi>f</mi><mrow><mi>α</mi><mi>μ</mi><mi>β</mi><mi>ν</mi></mrow></msup><mo>+</mo><msup><mi>f</mi><mrow><mi>α</mi><mi>μ</mi><mi>ν</mi><mi>β</mi></mrow></msup></mrow><mo>}</mo></mrow><mo stretchy="false">(</mo><mi>p</mi><mo>+</mo><mi>m</mi><mo stretchy="false">)</mo><mtext> </mtext><mn>,</mn></mtd></mtr></mtable></math> where Q em denotes an electromagnetic charge of a quark and the gamma matrix structures are given by <math display="block" id="M163"><mrow><msup><mi>σ</mi><mrow><mi>α</mi><mi>β</mi></mrow></msup><mtext> </mtext><mo>=</mo><mfrac><mi>i</mi><mn>2</mn></mfrac><mo stretchy="false">[</mo><msup><mi>γ</mi><mi>α</mi></msup><mn>,</mn><msup><mi>γ</mi><mi>β</mi></msup><mo stretchy="false">]</mo><mtext> </mtext><mn>,</mn></mrow></math> <math display="block" id="M164"><mrow><msup><mi>f</mi><mrow><mi>α</mi><mi>β</mi><mi>μ</mi><mi>ν</mi></mrow></msup><mo>=</mo><msup><mi>γ</mi><mi>α</mi></msup><mo stretchy="false">(</mo><mi>p</mi><mo>+</mo><mi>m</mi><mo stretchy="false">)</mo><msup><mi>γ</mi><mi>β</mi></msup><mo stretchy="false">(</mo><mi>p</mi><mo>+</mo><mi>m</mi><mo stretchy="false">)</mo><msup><mi>γ</mi><mi>μ</mi></msup><mo stretchy="false">(</mo><mi>p</mi><mo>+</mo><mi>m</mi><mo stretchy="false">)</mo><msup><mi>γ</mi><mi>ν</mi></msup><mtext> </mtext><mn>.</mn></mrow></math> The coupling strength can be read off from the sum of the two diagrams <math id="M165"><mrow><mi>i</mi><msup><mi mathvariant="script">M</mi><mi>μ</mi></msup><mo>=</mo><mi>i</mi><msubsup><mi mathvariant="script">M</mi><mi>a</mi><mi>μ</mi></msubsup><mo>+</mo><mi>i</mi><msubsup><mi mathvariant="script">M</mi><mi>b</mi><mi>μ</mi></msubsup></mrow></math> , where the each diagram is written down as <math display="block" id="M166"><mtable columnalign="left"><mtr><mtd><mi>i</mi><msubsup><mi mathvariant="script">M</mi><mi>a</mi><mi>μ</mi></msubsup><mo>=</mo><mo>−</mo><mtext>​</mtext><mtext>​</mtext><mstyle displaystyle="true"><mrow><mo>∫</mo></mrow></mstyle><mtext>​</mtext><mfrac><mrow><msup><mi>d</mi><mn>4</mn></msup><mi>p</mi></mrow><mrow><msup><mrow><mo stretchy="false">(</mo><mn>2</mn><mi>π</mi><mo stretchy="false">)</mo></mrow><mn>4</mn></msup></mrow></mfrac><mtext>T</mtext><mi>r</mi><mrow><mo>[</mo><mrow><mtext> </mtext><msubsup><mi>Γ</mi><mn>5</mn><mo>†</mo></msubsup><mo stretchy="false">(</mo><mi>p</mi><mo>−</mo><mi>q</mi><mn>,</mn><mi>p</mi><mo stretchy="false">)</mo><msub><mi>S</mi><mn>1</mn></msub><mo stretchy="false">(</mo><mi>p</mi><mo stretchy="false">)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mtext> </mtext><mrow><mrow><mo>×</mo><msup><mi>Γ</mi><mi>μ</mi></msup><mo stretchy="false">(</mo><mi>p</mi><mn>,</mn><mi>p</mi><mo>−</mo><mi>q</mi><mo stretchy="false">)</mo><msub><mi>S</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mi>p</mi><mo>−</mo><mi>q</mi><mo stretchy="false">)</mo><mtext> </mtext></mrow><mo>]</mo></mrow><mtext> </mtext><mn>,</mn></mtd></mtr></mtable></math> <math display="block" id="M167"><mtable columnalign="left"><mtr><mtd><mi>i</mi><msubsup><mi mathvariant="script">M</mi><mi>b</mi><mi>μ</mi></msubsup><mo>=</mo><mo>−</mo><mtext>​</mtext><mtext>​</mtext><mstyle displaystyle="true"><mrow><mo>∫</mo></mrow></mstyle><mtext>​</mtext><mfrac><mrow><msup><mi>d</mi><mn>4</mn></msup><mi>p</mi></mrow><mrow><msup><mrow><mo stretchy="false">(</mo><mn>2</mn><mi>π</mi><mo stretchy="false">)</mo></mrow><mn>4</mn></msup></mrow></mfrac><mtext>T</mtext><mi>r</mi><mrow><mo>[</mo><mrow><mtext> </mtext><msubsup><mi>Γ</mi><mn>5</mn><mo>†</mo></msubsup><mo stretchy="false">(</mo><mi>p</mi><mo>+</mo><mi>q</mi><mn>,</mn><mi>p</mi><mo stretchy="false">)</mo><msub><mi>S</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mi>p</mi><mo>+</mo><mi>q</mi><mo stretchy="false">)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mtext> </mtext><mrow><mrow><mo>×</mo><msup><mi>Γ</mi><mi>μ</mi></msup><mo stretchy="false">(</mo><mi>p</mi><mn>,</mn><mi>p</mi><mo>+</mo><mi>q</mi><mo stretchy="false">)</mo><msub><mi>S</mi><mn>1</mn></msub><mo stretchy="false">(</mo><mi>p</mi><mo stretchy="false">)</mo><mtext> </mtext></mrow><mo>]</mo></mrow><mtext> </mtext><mn>.</mn></mtd></mtr></mtable></math> In the leading order of the heavy-quark expansion, the p 0 -integral is easily performed, and we find <math display="block" id="M168"><mtable columnalign="left"><mtr><mtd><mi>i</mi><msubsup><mi mathvariant="script">M</mi><mi>a</mi><mi>μ</mi></msubsup><mo>=</mo><mi>i</mi><msubsup><mi mathvariant="script">M</mi><mi>b</mi><mi>μ</mi></msubsup></mtd></mtr><mtr><mtd><mo>=</mo><mn>2</mn><msub><mi>Q</mi><mrow><mtext>e</mtext><mi>m</mi></mrow></msub><msup><mover accent="true"><mi>F</mi><mo>˜</mo></mover><mrow><mn>0</mn><mi>μ</mi></mrow></msup><mstyle displaystyle="true"><mrow><mo>∫</mo></mrow></mstyle><mfrac><mrow><msup><mi>d</mi><mn>3</mn></msup><mtext>p</mtext></mrow><mrow><msup><mrow><mo stretchy="false">(</mo><mn>2</mn><mi>π</mi><mo stretchy="false">)</mo></mrow><mn>3</mn></msup></mrow></mfrac><msup><mrow><mo>|</mo><msub><mi>ψ</mi><mrow><mn>1</mn><mi>S</mi></mrow></msub><mo stretchy="false">(</mo><mi>p</mi><mo stretchy="false">)</mo><mo>|</mo></mrow><mn>2</mn></msup><mtext> </mtext><mn>.</mn></mtd></mtr></mtable></math> From the normalization of the wave function, <math display="block" id="M169"><mrow><mstyle displaystyle="true"><mrow><mo>∫</mo></mrow></mstyle><mfrac><mrow><msup><mi>d</mi><mn>3</mn></msup><mi>p</mi></mrow><mrow><msup><mrow><mo stretchy="false">(</mo><mn>2</mn><mi>π</mi><mo stretchy="false">)</mo></mrow><mn>3</mn></msup></mrow></mfrac><mo>|</mo><msub><mi>ψ</mi><mrow><mn>1</mn><mi>S</mi></mrow></msub><mo stretchy="false">(</mo><mi>p</mi><msup><mrow><mo stretchy="false">)</mo><mo>|</mo></mrow><mn>2</mn></msup><mo>=</mo><mn>1,</mn></mrow></math> the amplitude is independent of the wave functions, and the sum of two triangle diagrams is obtained as <math display="block" id="M170"><mrow><mi>i</mi><msup><mi mathvariant="script">M</mi><mi>μ</mi></msup><mo>=</mo><mn>4</mn><msub><mi>Q</mi><mrow><mtext>e</mtext><mi>m</mi></mrow></msub><msup><mover accent="true"><mi>F</mi><mo>˜</mo></mover><mrow><mn>0</mn><mi>μ</mi></mrow></msup><mtext> </mtext><mn>.</mn></mrow></math> By contracting with the polarization vector for the longitudinal (transverse) vector state ϵ μ = (0,0,0,1) ( <math id="M171"><mrow><msup><mover accent="true"><mi>ϵ</mi><mo>˜</mo></mover><mi>μ</mi></msup><mo>=</mo><mo stretchy="false">(</mo><mn>0,</mn><msub><mi>n</mi><mo>⊥</mo></msub><mn>,0</mn><mo stretchy="false">)</mo></mrow></math> ), the amplitude vanishes for the transverse modes as <math display="block" id="M172"><mrow><mi>i</mi><msub><mi mathvariant="script">M</mi><mi>μ</mi></msub><msup><mover accent="true"><mi>ϵ</mi><mo>˜</mo></mover><mi>μ</mi></msup><mo>=</mo><mn>0</mn><mtext> </mtext><mn>,</mn></mrow></math> while the longitudinal mode has a nonvanishing amplitude <math display="block" id="M173"><mrow><mi>i</mi><msub><mi mathvariant="script">M</mi><mi>μ</mi></msub><msup><mi>ϵ</mi><mi>μ</mi></msup><mo>=</mo><mn>4</mn><msub><mi>Q</mi><mrow><mtext>e</mtext><mi>m</mi></mrow></msub><msup><mover accent="true"><mi>F</mi><mo>˜</mo></mover><mrow><mn>0</mn><mi>μ</mi></mrow></msup><msub><mi>ϵ</mi><mi>μ</mi></msub><mo>=</mo><mn>4</mn><msub><mi>Q</mi><mrow><mtext>e</mtext><mi>m</mi></mrow></msub><mi>B</mi><mtext> </mtext><mn>.</mn></mrow></math> Therefore, the coupling strength in Eq. (25) is found to be <math display="block" id="M174"><mrow><msub><mi>g</mi><mrow><msub><mrow/><mrow><mtext>PV</mtext></mrow></msub></mrow></msub><mo>=</mo><mn>4</mn><msub><mi>Q</mi><mrow><mtext>em</mtext></mrow></msub><mtext> </mtext><mn>.</mn></mrow></math> The coupling strength depends only on the electric charge and is given by <math id="M175"><mrow><msub><mi>g</mi><mrow><msub><mrow/><mrow><mtext>PV</mtext></mrow></msub></mrow></msub><mo>=</mo><mn>8</mn><mo>/</mo><mn>3</mn><mo>≃</mo><mn>2.66</mn></mrow></math> ( <math id="M176"><mrow><msub><mi>g</mi><mrow><msub><mrow/><mrow><mtext>PV</mtext></mrow></msub></mrow></msub><mo>=</mo><mn>4</mn><mo>/</mo><mn>3</mn><mo>≃</mo><mn>1.33</mn></mrow></math> ) for the transition between η c and J/ψ ( η b and Υ ). This is consistent with the value obtained by fitting the measured radiative decay width [see below Eq. (33)] , but is slightly overestimated. The radiative decay widths in J / ψ → η c + γ and Υ→ η b + γ computed with the coupling strength (75) agree with the leading-order results by the potential Non-Relativistic QCD (pNRQCD) [ 208 208 , 209 209 ]. The overestimate can be improved with the inclusion of the subleading terms [ 208 208 , 209 209 ]. Extension to the open heavy flavors was done in Ref. [ 31 31 ]. Nuclear Science and Techniques 02 28