We review the properties of quarkonia under strong magnetic fields. The main phenomena are (i) mixing between different spin eigenstates, (ii) quark Landau levels and deformation of wave function, (iii) modification of $bar{Q}Q$ potential, and (iv) t
he motional Stark effect. For theoretical approaches, we review (i) constituent quark models, (ii) effective Lagrangians, (iii) QCD sum rules, and (iv) holographic approaches.
We investigate time evolution of $S$-wave charmonium populations under a time-dependent homogeneous magnetic field and evaluate survival probabilities of the low-lying charmonia to the goal of estimating the magnetic field strength at heavy-ion colli
sions. Our approach implements mixing between different spin eigenstates and transitions to radially excited states. We show that the survival probabilities can change even by an extremely short magnetic field. Furthermore, we find that the survival probabilities depend on the initial spin states. We propose the sum of the survival probabilities over spin partners as an observable insensitive to the initial states. We also find that the sum can be approximately given as a function of $sigma B_0^2$ with a duration time $sigma$ and the maximum strength of the magnetic field $B_0$.
We study the radiative (E1 and M1) decays of P-wave quarkonia in a strong magnetic field based on the Lagrangian of potential nonrelativistic QCD. To investigate their properties, we implement a polarized wave function basis justified in the Paschen-
Back limit. In a magnetic field stronger than the spin-orbit coupling, the wave functions of the P-wave quarkonia are drastically deformed by the Hadronic Paschen-Back effect. Such deformation leads to the anisotropy of the direction of decays from the P-wave quarkonia. The analytic formulas for the radiative decay widths in the nonrelativistic limit are shown, and the qualitative decay properties are discussed.
