Misalignment dynamics of Scalar Condensates with Yukawa coupling: Particle and Entropy Production

Misalignment dynamics, the non-equilibrium evolution of a scalar (or pseudoscalar) condensate in a potential landscape, broadly describes a solution to the strong CP problem, a mechanism for cold dark matter production and (pre) reheating post inflation. Often, radiative corrections are included phenomenologically by replacing the potential by the effective potential which is a \emph{static quantity}, its usefulness is restricted to (near) equilibrium situations. We study the misalignment dynamics of a scalar condensate Yukawa coupled to $N_f$ fermions in a manifestly energy conserving, fully renormalized Hamiltonian framework. A large $N_f$ limit allows us to focus on the fermion degrees of freedom, which yield a negative contribution to the effective potential, a radiatively induced instability and ultraviolet divergent field renormalization. We introduce an adiabatic basis and an adiabatic expansion that embodies the derivative expansion in the effective action, the zeroth order is identified with the effective potential, higher orders account for the derivative expansion including field renormalization and describe profuse particle production. Energy conserving dynamics leads to the conjecture of emergent asymptotic highly excited stationary states with a distribution function $n_k(\infty)\propto 1/k^6$ and an extensive entropy which is identified with an entanglement entropy. Subtle aspects of renormalization associated with the initial value problem are analyzed and resolved. Possible new manifestations of asymptotic spontaneous symmetry breaking (SSB) as a consequence of the dynamics even in absence of tree level (SSB), and cosmological inferences are discussed.

💡 Research Summary

The paper presents a comprehensive, energy‑conserving Hamiltonian treatment of the non‑equilibrium evolution of a homogeneous scalar (or pseudoscalar) condensate that is Yukawa‑coupled to (N_f) fermion species. By taking the large‑(N_f) limit the authors isolate the fermionic one‑loop effects while suppressing bosonic quantum fluctuations. First, they recover the static one‑loop effective potential for a constant background field: the fermion Dirac sea contributes a negative term (-\frac{N_f}{8\pi^2}m_f^4(\phi)\ln(m_f^2/\mu^2)), producing a Coleman–Weinberg‑type maximum and signalling a radiatively induced instability even when the tree‑level potential has no symmetry‑breaking minima.

Moving to the dynamical case, the exact Heisenberg equations for the scalar field and its conjugate momentum are combined with the Dirac equations for the fermions, whose mass (m_f(t)=y\langle\phi(t)\rangle) becomes time‑dependent. The authors introduce an adiabatic basis: instantaneous fermionic mode functions (u_k^{\pm}(t)) diagonalize the quadratic fermion Hamiltonian at each moment. A Bogoliubov transformation to this basis defines “adiabatic particles”. Expanding the mode functions in an adiabatic (derivative) series, the zeroth order reproduces the static effective potential, while higher orders generate derivative operators such as (\partial_\mu\phi,\partial^\mu\phi). These higher‑order terms are precisely the wave‑function renormalization pieces that appear in the derivative expansion of the full effective action.

Renormalization is performed with dimensional regularization and minimal subtraction. The fermion loop yields a logarithmically divergent field‑strength renormalization (\delta Z_\phi = \frac{N_f y^2}{12\pi^2}\ln(\Lambda/\mu)), which is absorbed into a redefinition of the kinetic term. Consequently, the renormalized equations of motion contain finite, time‑dependent mass and coupling counterterms that run with the renormalization scale, guaranteeing UV finiteness and exact energy conservation.

The time‑dependent background drives copious fermion‑pair production. By solving the mode equations (analytically in limiting regimes and numerically otherwise) the authors find that the occupation numbers of the adiabatic particles settle at late times into a non‑thermal power‑law spectrum
\