Figures:

                                                                      Figure 1

Fig 1:  Sketch of a heavy-ion collision. The impact parameter and the beam direction are along the $x$ and $z$ axes, respectively. The $x$-$z$ plane is called the reaction plane. Participating nucleons in the overlap region create a hot and dense medium of  quark-gluon plasma. Spectator nuclear fragments generate strong electromagnetic fields.        

                                                                                                                                                    
                                                                                 
Fig 2: Schematic side view of a heavy-ion collision. The dashed lines represent the motion of quarks due to the QGP expansion.The black solid lines indicate the Hall effect, where the magnetic field exerts the Lorentz force on moving charges. The impact of the Coulomb field from spectators and the Faraday field induced by the fast decay of the magnetic field is marked by red lines.

Fig 3: Illustration of different contributions to $\Delta(dv_1/dy)$ between protons and anti-protons versus centrality.  Panel (a) sketches the transported-quark effect. Panel (b) depicts the electromagnetic-field 

Fig 4: Panel (a): the ionization energy loss $dE/dx$ of charged particles as a function of momentum in TPC. The solid lines are the corresponding Bichsel function. Panel (b): the mass-squared as a function of momentum in TOF. The dashed boxed areas indicate the pion, kaon and proton candidates.

                                                                                             
 
Fig 5: $v_1$ for protons and anti-protons as a function of rapidity in (a) Au+Au collisions at $\sqrt{s_{NN}}=$ 200 GeV, (b) isobar (Ru+Ru and Zr+Zr) collisions at $\sqrt{s_{NN}}=$ 200 GeV, and (c) Au+Au collisions at $\sqrt{s_{NN}}=$ 27 GeV in the centrality interval of 50--80\%. Protons and anti-protons are marked with solid and open circles, respectively. Panels (d), (e) and (f) show $\Delta v_1 \equiv v_1^{p} - v_1^{\bar{p}}$ versus rapidity. The $d\Delta v_1/d\mathsf{y}$ values are obtained with linear fits (solid lines). Systematic uncertainties are indicated with shaded bands. 

Conclusions: