We study fermion mass correction to chiral kinetic equations in electromagnetic fields. Different from the chiral limit where fermion number density is the only independent distribution, the number and spin densities are coupled to each other for massive fermion systems. To the first order in h, we derived the quantum correction to the classical on-shell condition and the Boltzmann-type transport equations. To the linear order in the fermion mass, the mass correction does not change the structure of the chiral kinetic equations and behaves like additional collision terms. While the mass correction exists already at classical level in general electromagnetic fields, it is only a first order quantum correction in the study of chiral magnetic effect.

We study neutral and charged meson properties in the magnetic field. Taking the bosonization method in a two-flavor Nambu-Jona-Lasinio model, we derive effective meson Lagrangian density with minimal coupling to the magnetic field, by employing derivative expansion for both the meson fields and Schwinger phases. We extract from the effective Lagrangian density the meson curvature, pole and screening masses. As the only Goldstone mode, the neutral pion controls the thermodynamics of the system and propagates the long range quark interaction. The magnetic field breaks down the space symmetry, and the quark interaction region changes from a sphere in vacuum to a ellipsoid in magnetic field.

As the core ingredient for spin polarization, the equilibrium spin distribution function that eliminates the collision terms is derived from the detailed balance principle. The kinetic theory for interacting fermionic systems is applied to the Nambu–Jona-Lasinio model at quark level. Under the semi-classical expansion with respect to $$\hbar $$ ħ , the kinetic equations for the vector and axial-vector distribution functions are obtained with collision terms. For an initially unpolarized system, spin polarization can be generated at the first order of $$\hbar $$ ħ from the coupling between the vector and axial-vector charges. Different from the classical transport theory, the collision terms in a quantum theory vanish only in global equilibrium with Killing condition.

We investigate the critical behavior of pion superfluidity in the frame of functional renormalization group. By solving the flow equation in the SU(2) linear sigma model with pion superfluidity at finite temperature and isospin density, and making comparison with the fixed point analysis of the flow diagram in a general O(N ) model at zero temperature and density but with continuous dimension, the pion superfluidity is a second order phase transition subject to a O(2) universality class with a dimension crossover from d eff = 4 to d eff = 3. The usually used large N expansion gives a temperature independent critical exponent β and agrees with the renormalization group only at zero temperature.