\section{Introduction}
\label{}
Non-central heavy ion collisions produce extremely strong magnetic fields~\cite{Kharzeev:2007jp,Skokov:2009qp,Tuchin:2010vs}. For such strong fields, electromagnetic interactions significantly affect the evolution of the produced quark-gluon plasma (QGP). In particular, the interplay between magnetic fields and chiral anomaly has been predicted to lead to interesting phenomena including the chiral magnetic effect~\cite{kharzeev2006parity,Kharzeev:2007jp,Fukushima:2008xe}
and the chiral magnetic wave~\cite{Kharzeev:2010gd,Burnier:2011bf}. These quantum chiral effects, if observed, would reveal the existence of topological charge fluctuations in the QGP that are similar to the topological fluctuations in the primordial electroweak plasma responsible for the baryon asymmetry of the Universe. 

This motivation drives the need for establishing the presence and the magnitude of electromagnetic fields in the QGP. Here we report on the results of our studies \cite{Gursoy:2014aka, gursoy2018charge}  of the effects of electromagnetic fields on the collective flow of charged hadrons. 
Our treatment of electromagnetic effects is perturbative: we solve the Maxwell equations to evaluate the currents induced by electromagnetism on top of the conventional hydrodynamical evolution (no back-reaction is included). 

\section{Electromagnetic fields in heavy ion collisions: a qualitative discussion}
\begin{figure}[t]
 \begin{center}
\includegraphics[scale=0.35]{flowschemaColorv4}
 \end{center}
 \vspace{-0.3in}
 \caption{Schematic illustration of how the electromagnetic fields in a heavy ion collision result
in a directed flow of electric charge, $\Delta v_1$, see text for explanation. From \cite{gursoy2018charge}.}
\label{figschema}
\end{figure}


\begin{figure}[t!]
 \begin{center}
\includegraphics[scale=0.28]{E_field_tau_1_eta_0}
\includegraphics[scale=0.29]{B_field_tau_1_eta_0}
 \end{center}
 \caption[]{The electric (left) and magnetic (right) fields in the transverse plane at $z=0$  in the lab frame 
 at a proper time $\tau=1$~fm/c after a Pb+Pb collision
 with 20-30\% centrality (impact parameters in the range 6.24~fm~$<b<$~9.05~fm) and with a collision energy $\sqrt{s} =2.76$~ATeV.   From \cite{gursoy2018charge}. }
\label{figEB}
\end{figure}

As illustrated in Fig.~\ref{figschema}, there are three distinct electromagnetic effects 
on charged components of the fluid, resulting in a sideways current:
\begin{enumerate}
\item {\it Faraday:} as the magnetic field decreases in time, Faraday's law implies the induction of an electric field and, since the QGP possesses mobile charges, an electric current. We denote this electric field by $\vec E_F$.  Since $\vec E_F$ curls around the (decreasing in time) $\vec B$ that points in the $y$-direction, the sideways component of $E_F$ points in opposite directions at opposite rapidity, see Fig.~\ref{figschema}.
\item {\it Lorentz:} since the QGP fluid exhibits a strong longitudinal flow velocity $\vec{v}_\mathrm{flow}$ denoted by $\vec u$ in Fig.~\ref{figschema} pointing along the beam direction (and hence perpendicular to $\vec B$), the Lorentz force exerts a sideways push on charged particles in opposite directions at opposite rapidity. Equivalently, upon boosting to the local fluid rest frame, the lab frame $\vec B$ yields a fluid frame $\vec E$ in which the effects on the charged components of the fluid are equivalent to the effects of the Lorentz force in the lab frame. We denote this electric field by $\vec E_L$.  Both $\vec E_F$ and $\vec E_L$ originate from the magnetic field.
\item {\it Coulomb:} The positive spectators exert an electric force on the charged plasma produced in the collision, which again points in opposite directions at opposite rapidity; the corresponding 
electric field is denoted by $\vec E_C$. 
\end{enumerate}
All three of these electric fields --- and the electric currents induced by them --- have opposite directions at positive and negative rapidity. As shown in Fig.~\ref{figschema}, 
$\vec E_F$ and  $\vec E_C$ have the same sign, whereas  $\vec E_L$ opposes them. Therefore, 
 the sign of the total rapidity-odd, charge-odd, directed flow charge dependence $\Delta v_1$ that results 
 from these electric fields
 depends on whether  $\vec E_F + \vec E_C$ or $\vec E_L$ is dominant. The relative magnitude of these fields is quite sensitive to the space-time evolution of the fluid. 



\section{Electromagnetic fields in heavy ion collisions: numerical results}

Our numerical results were obtained, as explained in \cite{gursoy2018charge}, by describing the evolution of the QGP fluid within the {\tt iEBE-VISHNU} framework \cite{Shen:2014vra}, and by solving Maxwell equations on top of this hydrodynamical background.  The 
hydrodynamic calculation 
assumes longitudinal boost-invariance and begins at $\tau_0 = 0.4$ fm/$c$. We evolve the relativistic viscous
hydrodynamic equations for a fluid with an equation of state based upon lattice QCD calculations \cite{Huovinen:2009yb}. The electromagnetic fields generated by the charges and currents evolve according to the Maxwell equations, in which we assume a constant conductivity. 

Our calculation yields the momentum distribution for hadrons with different charge, from which we  evaluate  the difference between the flow for positively and negatively charged hadrons 
$\Delta v_n \equiv v_n(h^+) - v_n(h^-)$ 
shown in Fig.~\ref{fig5} for the top RHIC and LHC energies.

\begin{figure*}[t!]
  \includegraphics[width=0.6\linewidth]{EM_induced_flow_ecm_dependence.pdf}   
  \caption{The collision energy dependence of the electromagnetically induced charge-odd contributions to flow observables. The difference of particle mean $p_T$ and $v_n$ between $\pi^+$ and $\pi^-$ are plotted as a function of particle rapidity for collisions at the top RHIC energy of 200 GeV and at two LHC collision energies; from \cite{gursoy2018charge}..}
  \label{fig5}
\end{figure*}


\section{Conclusions}

The data from RHIC \cite{Adam:2019wnk} and LHC \cite{Acharya:2019ijj} show that the observed $\Delta v_1$ is similar in magnitude to the one that we have predicted, but has an {\it opposite sign}. This 
indicates that the electromagnetic fields of expected magnitude are indeed present in the QGP, but our treatment of magnetohydrodynamics of the plasma has to be improved. This can be done along two directions: i) improve the computation within the perturbative framework, e.g. by including the time-dependent conductivity; and ii) perform a full magnetohydrodynamical (MHD) simulation. The recent pilot MHD simulation \cite{Inghirami:2019mkc} however indicates the same sign of $\Delta v_1$ as our perturbative result. Clearly, more work is needed to improve our understanding of this important issue.

\section{Acknowledgements}
This work was supported in part by the Netherlands Organisation for Scientific Research (NWO) under VIDI grant 680-47-518, the Delta Institute for Theoretical Physics (D-ITP) funded by the Dutch Ministry of Education, Culture and Science (OCW), the Scientific and Technological Research Council of Turkey (TUBITAK), the Office of Nuclear Physics of the U.S.~Department of Energy under Contract Numbers  DE-SC0011090, DE-FG-88ER40388 and DE-AC02-98CH10886, and the Natural Sciences and Engineering Research Council of Canada. KR gratefully acknowledges the hospitality of the CERN Theory Group. CS gratefully acknowledges a Goldhaber Distinguished Fellowship from Brookhaven Science Associates. Computations were made in part on the supercomputer Guillimin from McGill University, managed by Calcul Qu\'ebec and Compute Canada. The operation of this supercomputer is funded by the Canada Foundation for Innovation (CFI), NanoQu\'ebec, RMGA and the Fonds de recherche du Qu\'ebec -- Nature et technologies (FRQ-NT). UG is grateful for the hospitality of the Bo\~gazi\c ci University and the Mimar Sinan University in Istanbul. We gratefully acknowledge helpful discussions with Gang Chen, Ulrich Heinz, Jacopo Margutti, Raimond Snellings, Sergey Voloshin and Fuqiang Wang. 





\bibliographystyle{elsarticle-num}
