Oscillon decay via parametric resonance: the case of three-point scalar interactions

Basic Information

• Paper ID: 2511.03501
• Title: Oscillon decay via parametric resonance: the case of three-point scalar interactions
• Author: Siyao Li (Institute of Science Tokyo & IBS Center for Theoretical Physics)
• Categories: hep-ph, astro-ph.CO, gr-qc, hep-th
• Publication Date: November 5, 2025 (arXiv v1)
• Paper Link: https://arxiv.org/abs/2511.03501

Abstract

This paper investigates the decay dynamics of oscillons through interactions with external scalar fields. To test the robustness of the parametric resonance decay mechanism previously discovered in four-point couplings, the authors extend their analysis to three-point interactions $g_3\phi\chi^2$. By computing Floquet exponents for the external field $\chi$ in the oscillon background, the authors analyze how instability bands depend on coupling constants and oscillon shape. Numerical simulations of the two-field system reveal that, similar to the four-point coupling case, parametric resonance may terminate before complete oscillon destruction, leaving a smaller oscillon that decays only through perturbative channels. This demonstrates that partial oscillon decay via parametric resonance is a universal phenomenon for oscillon-scalar couplings, qualitatively independent of the specific interaction form, while the instability band shape, parameter dependence, and precise critical energy depend on the particular coupling.

Research Background and Motivation

1. Core Problem

Oscillons are localized, oscillating soliton solutions in nonlinear potentials of real scalar fields, capable of long-term existence due to suppressed radiation decay channels. During the preheating phase following inflation, coherent oscillations of the inflaton field can lead to oscillon formation through parametric resonance. Understanding how oscillons couple to external fields and decay is crucial for the cosmological reheating process.

2. Problem Significance

• Cosmological Importance: Oscillons play important roles in post-inflationary cosmic evolution, reheating processes, and dark matter formation scenarios
• Theoretical Value: Coupling of oscillons to external fields can modify decay rates through particle production, with parametric resonance enabling rapid decay
• Physical Mechanism: Unlike uniform inflationary backgrounds, the spatial localization of oscillons can suppress resonance, as produced particles can escape from the oscillon

3. Existing Research Limitations

In prior work Ref. 1, the authors studied four-point couplings $g_4\phi^2\chi^2$ and found that resonance terminates before complete oscillon destruction when the oscillon energy drops below a critical value. The key question is: To what extent does this result depend on the specific interaction form?

4. Research Motivation

This paper addresses this by studying three-point interactions $g_3\phi\chi^2$ to answer:

• Which decay dynamics features are universal?
• Which features undergo qualitative changes with different interaction forms?
• Is partial oscillon decay a universal phenomenon?

Core Contributions

1. Extended Floquet Analysis: First systematic analysis of the instability band structure for the external field $\chi$ under three-point coupling $g_3\phi\chi^2$, computing Floquet exponent dependence on coupling constants and oscillon shape
2. Universality Revelation: Demonstrates that partial oscillon decay via parametric resonance is a universal phenomenon, qualitatively independent of specific coupling form (three-point or four-point)
3. Coupling-Dependent Features Identification: Clarifies that instability band shape, parameter dependence relationships, and critical oscillon energy $E_0^{\text{osc}}$ depend on the specific coupling
4. Complete Numerical Simulations: Performs two-field dynamics simulations for both $g_3 \neq 0, g_4 = 0$ and $g_3 \neq 0, g_4 \neq 0$ cases, validating theoretical predictions
5. Cosmological Applications: Provides deeper physical insights for understanding the role of post-inflationary oscillons in the reheating process

Methodology Details

Task Definition

Input:

• Self-interaction potential of oscillon field $\phi$: $V(\phi) = \frac{1}{2}m_\phi^2\phi^2 - \lambda\phi^4 + \lambda_6\phi^6$
• Mass of external scalar field $\chi$: $m_\chi$
• Coupling parameters: $g_3, g_4$

Output:

• Floquet exponents $\mu$ of external field (determining instability bands)
• Critical coupling $(g_3)_0$ and critical oscillon energy $E_0^{\text{osc}}$
• Complete decay dynamics evolution

Constraints:

• Oscillon stability condition: $\min[2V(\psi)/\psi^2] < \omega^2 < m_\phi^2$
• Energy conservation and spherical symmetry

Model Architecture

1. Two-Field Lagrangian

$\mathcal{L} = \frac{1}{2}\partial_\mu\phi\partial^\mu\phi - V(\phi) + \frac{1}{2}\partial_\mu\chi\partial^\mu\chi - \frac{1}{2}m_\chi^2\chi^2 - g_3\phi\chi^2 - g_4\phi^2\chi^2$

2. Equations of Motion

$\ddot{\phi} - \nabla^2\phi + m_\phi^2\phi + \frac{dV_{nl}}{d\phi} + g_3\chi^2 + 2g_4\phi\chi^2 = 0$$\ddot{\chi} - \nabla^2\chi + m_\chi^2\chi + 2g_3\phi\chi + 2g_4\phi^2\chi = 0$

3. Oscillon Approximation

Employing single-frequency approximation: $\phi(t,r) \approx 2\psi(r)\cos(\omega t)$, where $\psi(r)$ is a localized profile function satisfying the boundary value problem: $\frac{d^2\psi}{dr^2} + \frac{2}{r}\frac{d\psi}{dr} - \left[(m_\phi^2-\omega^2)\psi + \frac{1}{2}\frac{\partial V_{eff}}{\partial\psi}\right] = 0$

Boundary conditions: $\frac{d\psi}{dr}|_{r=0} = 0$, $\psi(r\to\infty) \to 0$

Technical Innovations

1. Floquet Analysis Method

Uniform Approximation: Taking the oscillon center amplitude $\psi_0 = \psi(r=0)$, the equation of motion for $\chi$ simplifies to the standard Mathieu equation: $\chi_k'' + (A_k + 2q\cos(2z))\chi_k = 0$ where:

• $A_k = \frac{4(k^2+m_\chi^2)}{\omega^2}$
• $q = \frac{8g_3\psi_0}{\omega^2}$
• $z = \omega t/2$

Key Finding: Unlike four-point coupling, the parameters $A_k$ and $q$ are independent for three-point coupling, and the equation is symmetric about $g_3=0$ (attractive and repulsive interactions show identical behavior).

2. Narrow Band Resonance Conditions

First Narrow Band (when $m_\chi/m_\phi < 0.5$, allowing $\phi \to \chi\chi$): $\mu_{max} \approx |q|/2 = \frac{4\psi_0}{\omega^2}|g_3|, \quad |g_3| \ll 1$

Second Narrow Band (when $m_\chi/m_\phi > 0.5$, corresponding to $\phi\phi \to \chi\chi$): $\mu_{max} \approx q^2/16 = \frac{4\psi_0^2}{\omega^4}g_3^2, \quad |g_3| \ll 1$

3. Finite-Size Effects and Critical Coupling

Accounting for oscillon spatial localization, the particle escape rate is: $\Gamma_{escape} \sim \frac{v_\chi}{R} \sim \frac{1}{R}\sqrt{1 - \frac{4}{n^2}\left(\frac{m_\chi}{m_\phi}\right)^2}$

The critical coupling $(g_3)_0$ is determined by the condition $\Re(\mu) \simeq \Gamma_{escape}$: $|(g_3)_0| \sim \frac{\omega}{2\psi_0 R}\sqrt{1-4\left(\frac{m_\chi}{m_\phi}\right)^2}, \quad m_\chi \lesssim 0.5m_\phi$

This is the core innovation of the paper: establishing a quantitative relationship between critical coupling and oscillon geometric parameters.

4. Numerical Treatment of Non-Uniform Background

• Evolve $\chi$ field on fixed oscillon background $\phi(t,r) = 2\psi(r)\cos(\omega t)$
• Use fourth-order symplectic integrator for time evolution
• Apply adiabatic damping boundary conditions to eliminate non-physical reflections
• Extract Floquet exponents by fitting $E_\chi \propto e^{\Re(\mu)\omega t}$

Experimental Setup

Numerical Parameters

• Simulation Box Size: $r_{box} = 192 m_\phi^{-1}$
• Grid Points: 3072
• Time Step: $\Delta t = 0.01 m_\phi^{-1}$
• Monitoring Radius: $r_{max} = 30 m_\phi^{-1}$
• Spatial Derivatives: Fourth-order central differences
• Time Integration: Fourth-order symplectic integrator

Dimensionless Variables

All calculations use dimensionless variables: $\tilde{x}^\mu = m_\phi x^\mu, \quad \tilde{\phi} = \frac{\sqrt{\lambda}\phi}{m_\phi}, \quad \tilde{\chi} = \frac{\sqrt{\lambda}\chi}{m_\phi}$$\tilde{m}_\chi = \frac{m_\chi}{m_\phi}, \quad \tilde{g}_3 = \frac{g_3}{m_\phi\sqrt{\lambda}}, \quad \tilde{g}_4 = \frac{g_4}{\lambda}$

The only free parameter is $m_\phi^2\lambda_6/\lambda^2 = 0.8$ (adopted in this work)

Initial Conditions

Fixed Background Simulation: $\chi(0,r) = \chi_0 e^{-r^2/R_\chi^2}, \quad \dot{\chi}(0,r) = 0$ Taking $\chi_0 = 0.1$, $R_\chi = 6$

Two-Field Evolution Simulation: $\phi(0,r) = 2\psi(r), \quad \chi(0,r) = \tilde{\chi}_0 \bar{\chi}(r)$$\dot{\phi}(0,r) = \dot{\chi}(0,r) = 0$ Taking $\tilde{\chi}_0 = 0.01$

Energy Calculation

Time-averaged energy (averaging window $T_{ave} = 20-100 m_\phi^{-1}$): $E_\chi(t) = \frac{1}{T_{ave}}\int_t^{t+T_{ave}}dt \int_0^{r_{max}}dr\, 4\pi r^2\left[\frac{1}{2}\dot{\chi}^2 + \frac{1}{2}(\partial_r\chi)^2 + \frac{1}{2}m_\chi^2\chi^2\right]$

Experimental Results

Main Results

1. Floquet Spectrum Characteristics (Figure 2)

For an oscillon with $\omega/m_\phi = 0.8$:

• $m_\chi/m_\phi = 0.3$: First narrow band dominates, corresponding to $\phi \to \chi\chi$ decay channel
• $m_\chi/m_\phi = 0.6$: Second narrow band dominates ($\phi\phi \to \chi\chi$), as single-particle decay is kinematically forbidden

Key Observations:

• No tachyonic bands ($A_k \geq 0$ always holds)
• Symmetric about $g_3=0$
• Larger $m_\chi$ shifts Floquet diagram downward along $k$ axis

2. Maximum Growth Rate (Figure 3)

$\mu_{max} \equiv \max(\Re(\mu))$

• Small coupling region: Linear dependence $\mu_{max} \propto |g_3|$ for $m_\chi/m_\phi = 0.3$
• Large mass case: Quadratic dependence $\mu_{max} \propto g_3^2$ for $m_\chi/m_\phi = 0.6$
• Exactly matches theoretical predictions (Eqs. 4.7 and 4.8)

3. Finite-Size Suppression (Figure 4)

Numerical simulations with non-uniform oscillon background show:

• Significantly suppressed growth rates compared to uniform case (blue stable region expands)
• Existence of critical coupling $|(g_3)_0|$:
  • $m_\chi/m_\phi = 0.3$: $|(g_3)_0|/(m_\phi\sqrt{\lambda}) \approx 0.02$
  • $m_\chi/m_\phi = 0.6$: $|(g_3)_0|/(m_\phi\sqrt{\lambda}) \approx 0.1$
• Quantitatively consistent with theoretical estimate (Eq. 4.14)

4. Spectral Analysis (Figure 5)

Fourier spectrum of $\chi$ field center value shows:

• Main peaks at $\omega_k = \sqrt{k^2+m_\chi^2} = n\omega/2$ (n-th narrow band parametric resonance)
• Secondary peaks at $(n\pm1)\omega$ (frequency modulation effects from coupling terms)

Two-Field Dynamics Simulations (Figures 7-8)

Case 1: Complete Oscillon Destruction

Parameters: $m_\chi/m_\phi = 0.3$, $g_3/(m_\phi\sqrt{\lambda}) = 0.3$, initial $\omega/m_\phi = 0.8$

Results:

• Rapid exponential growth of $\chi$ field ($t < 50 m_\phi^{-1}$)
• Swift energy transfer from oscillon to $\chi$ field
• Since $E_0^{osc} \approx E_{death}$, oscillon is destroyed before resonance termination
• Final $E_\chi/E_{tot} > 0.95$ (nearly complete energy transfer)

Case 2: Premature Resonance Termination

Parameters: $m_\chi/m_\phi = 0.6$, $g_3/(m_\phi\sqrt{\lambda}) = 0.7$, initial $\omega/m_\phi = 0.8$

Results:

• Brief exponential growth phase ($t < 30 m_\phi^{-1}$)
• $\chi$ field exits instability band, growth stops
• Residual oscillon persists, decaying only through perturbative channels
• $E_\phi/E_{tot}$ stabilizes at approximately 0.4 (partial energy retention)

Case 3: Mixed Coupling (Figure 8)

Parameters: $m_\chi/m_\phi = 0.3$, $g_3/(m_\phi\sqrt{\lambda}) = 0.3$, $g_4/\lambda = 0.2$

Key Finding:

• Non-zero $g_4$ modifies Floquet exponents
• $E_0^{osc}$ decreases, resonance terminates before oscillon destruction
• Demonstrates universality of partial decay in multi-coupling scenarios

Quantitative Results Summary

Experimental Findings

1. Universality Confirmation: Three-point and four-point couplings exhibit qualitatively identical decay behavior patterns
2. Critical Energy Scaling: Oscillon decay rate follows piecewise law: