2025-11-30T12:52:18.148453

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

Li
We investigate the decay dynamics of oscillons through interactions with an external scalar field. To examine how robust the decay dynamics of oscillons via parametric resonance we previously found in Li et al. 2025 are to the specific form of the coupling, we extend the analysis to include a three-point interaction $g_3ϕχ^2$. We compute the Floquet exponents of the external field $χ$ under an oscillating oscillon background and analyze how the instability bands depend on the coupling constants and the oscillon shapes. Numerical simulations of the two-field system show that, similar to the four-point case, the parametric resonance may cease before the oscillon is destroyed, leaving a smaller oscillon that decays only perturbatively. This indicates that the partial decay of oscillons through parametric resonance is a generic phenomenon of oscillon-scalar couplings, qualitatively insensitive to the specific interaction form, while the shape of instability bands, parameter dependence, and the precise critical oscillon energies depend on the specific coupling. Our findings provide further insights into the decay dynamics of oscillons and their potential role in the post-inflationary reheating process.
academic

Decaimiento de oscilones mediante resonancia paramétrica: el caso de interacciones escalares de tres puntos

Información Básica

  • ID del Artículo: 2511.03501
  • Título: Oscillon decay via parametric resonance: the case of three-point scalar interactions
  • Autor: Siyao Li (Institute of Science Tokyo & IBS Center for Theoretical Physics)
  • Clasificación: hep-ph, astro-ph.CO, gr-qc, hep-th
  • Fecha de Publicación: 5 de noviembre de 2025 (arXiv v1)
  • Enlace del Artículo: https://arxiv.org/abs/2511.03501

Resumen

Este artículo investiga la dinámica de decaimiento de oscilones (oscillons) mediante interacción con campos escalares externos. Para examinar la robustez del mecanismo de decaimiento por resonancia paramétrica previamente encontrado en acoplamientos de cuatro puntos, los autores extienden el análisis a interacciones de tres puntos g3ϕχ2g_3\phi\chi^2. Mediante el cálculo de exponentes de Floquet del campo externo χ\chi en el fondo del oscilón, se analiza cómo las bandas de inestabilidad dependen de la constante de acoplamiento y la forma del oscilón. Las simulaciones numéricas del sistema de dos campos muestran que, de manera similar al caso de acoplamiento de cuatro puntos, la resonancia paramétrica puede terminar antes de la destrucción completa del oscilón, dejando un oscilón más pequeño que decae únicamente de manera perturbativa. Esto sugiere que el decaimiento parcial de oscilones mediante resonancia paramétrica es un fenómeno universal en acoplamientos oscilón-escalar, que cualitativamente no depende de la forma de interacción específica, aunque la forma de las bandas de inestabilidad, la dependencia de parámetros y la energía crítica exacta dependen del acoplamiento específico.

Contexto de Investigación y Motivación

1. Problema Central

Los oscilones son soluciones de solitones oscilatoria localizados en potenciales no lineales de campos escalares reales, que pueden existir durante largos períodos porque sus canales de decaimiento por radiación están suprimidos. Durante la fase de precalentamiento (preheating) posterior a la inflación, las oscilaciones coherentes del campo inflacionario pueden conducir a la formación de oscilones mediante resonancia paramétrica. Comprender cómo los oscilones se acoplan con campos externos y decaen es crucial para el proceso de recalentamiento cosmológico posterior a la inflación.

2. Importancia del Problema

  • Significado Cosmológico: Los oscilones juegan un papel importante en la evolución del universo posterior a la inflación, en procesos de recalentamiento y en la formación de materia oscura
  • Valor Teórico: El acoplamiento de oscilones con campos externos puede modificar su tasa de decaimiento mediante producción de partículas, particularmente la resonancia paramétrica puede conducir a decaimiento rápido
  • Mecanismo Físico: A diferencia del fondo inflacionario uniforme, la localización espacial del oscilón puede suprimir la resonancia, ya que las partículas producidas pueden escapar del oscilón

3. Limitaciones de Investigaciones Previas

En trabajos anteriores Ref. 1, los autores estudiaron el caso de acoplamiento de cuatro puntos g4ϕ2χ2g_4\phi^2\chi^2, encontrando que cuando la energía del oscilón cae por debajo de un valor crítico, la resonancia termina antes de la destrucción completa del oscilón. Sin embargo, la pregunta clave es: ¿En qué medida este resultado depende de la forma específica de la interacción?

4. Motivación de la Investigación

Este artículo responde mediante el estudio de interacciones de tres puntos g3ϕχ2g_3\phi\chi^2:

  • ¿Cuáles características de la dinámica de decaimiento son universales?
  • ¿Cuáles características cambian cualitativamente con cambios en la forma de interacción?
  • ¿Es el fenómeno de decaimiento parcial del oscilón una ley universal?

Contribuciones Principales

  1. Análisis de Floquet Extendido: Análisis sistemático por primera vez de la estructura de bandas de inestabilidad del campo externo χ\chi bajo acoplamiento de tres puntos g3ϕχ2g_3\phi\chi^2, calculando la dependencia de los exponentes de Floquet en la constante de acoplamiento y la forma del oscilón
  2. Revelación de Universalidad: Demostración de que el decaimiento parcial de oscilones mediante resonancia paramétrica es un fenómeno universal, cualitativamente independiente de la forma de acoplamiento específica (tres o cuatro puntos)
  3. Identificación de Dependencia de Interacción: Clarificación de que la forma de las bandas de inestabilidad, las relaciones de dependencia de parámetros y la energía crítica del oscilón E0oscE_0^{\text{osc}} dependen del acoplamiento específico
  4. Simulaciones Numéricas Completas: Realización de simulaciones de dinámica de dos campos incluyendo casos con g30,g4=0g_3 \neq 0, g_4 = 0 y g30,g40g_3 \neq 0, g_4 \neq 0, verificando predicciones teóricas
  5. Aplicaciones Cosmológicas: Proporciona una comprensión física más profunda del papel de los oscilones en el recalentamiento posterior a la inflación

Explicación Detallada de Métodos

Definición de Tareas

Entrada:

  • Potencial de autointeracción del campo oscilón ϕ\phi: V(ϕ)=12mϕ2ϕ2λϕ4+λ6ϕ6V(\phi) = \frac{1}{2}m_\phi^2\phi^2 - \lambda\phi^4 + \lambda_6\phi^6
  • Masa del campo escalar externo χ\chi: mχm_\chi
  • Parámetros de acoplamiento: g3,g4g_3, g_4

Salida:

  • Exponentes de Floquet del campo externo μ\mu (determinan bandas de inestabilidad)
  • Acoplamiento crítico (g3)0(g_3)_0 y energía crítica del oscilón E0oscE_0^{\text{osc}}
  • Evolución dinámica completa del decaimiento

Condiciones de Restricción:

  • Condición de estabilidad del oscilón: min[2V(ψ)/ψ2]<ω2<mϕ2\min[2V(\psi)/\psi^2] < \omega^2 < m_\phi^2
  • Conservación de energía y simetría esférica

Arquitectura del Modelo

1. Lagrangiano de Dos Campos

L=12μϕμϕV(ϕ)+12μχμχ12mχ2χ2g3ϕχ2g4ϕ2χ2\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. Ecuaciones de Movimiento

ϕ¨2ϕ+mϕ2ϕ+dVnldϕ+g3χ2+2g4ϕχ2=0\ddot{\phi} - \nabla^2\phi + m_\phi^2\phi + \frac{dV_{nl}}{d\phi} + g_3\chi^2 + 2g_4\phi\chi^2 = 0χ¨2χ+mχ2χ+2g3ϕχ+2g4ϕ2χ=0\ddot{\chi} - \nabla^2\chi + m_\chi^2\chi + 2g_3\phi\chi + 2g_4\phi^2\chi = 0

3. Aproximación del Oscilón

Se adopta la aproximación de frecuencia única: ϕ(t,r)2ψ(r)cos(ωt)\phi(t,r) \approx 2\psi(r)\cos(\omega t), donde ψ(r)\psi(r) es una función de perfil localizado que satisface el problema de valores en la frontera: d2ψdr2+2rdψdr[(mϕ2ω2)ψ+12Veffψ]=0\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

Condiciones de frontera: dψdrr=0=0\frac{d\psi}{dr}|_{r=0} = 0, ψ(r)0\psi(r\to\infty) \to 0

Puntos de Innovación Técnica

1. Método de Análisis de Floquet

Aproximación Uniforme: Tomando la amplitud en el centro del oscilón ψ0=ψ(r=0)\psi_0 = \psi(r=0), la ecuación de movimiento de χ\chi se simplifica a la ecuación de Mathieu estándar: χk+(Ak+2qcos(2z))χk=0\chi_k'' + (A_k + 2q\cos(2z))\chi_k = 0 donde:

  • Ak=4(k2+mχ2)ω2A_k = \frac{4(k^2+m_\chi^2)}{\omega^2}
  • q=8g3ψ0ω2q = \frac{8g_3\psi_0}{\omega^2}
  • z=ωt/2z = \omega t/2

Hallazgo Clave: A diferencia del acoplamiento de cuatro puntos, los parámetros AkA_k y qq en el acoplamiento de tres puntos son independientes, y la ecuación es simétrica respecto a g3=0g_3=0 (los comportamientos de interacción atractiva y repulsiva son idénticos).

2. Condiciones de Resonancia de Banda Estrecha

Primera Banda Estrecha (cuando mχ/mϕ<0.5m_\chi/m_\phi < 0.5, permitiendo ϕχχ\phi \to \chi\chi): μmaxq/2=4ψ0ω2g3,g31\mu_{max} \approx |q|/2 = \frac{4\psi_0}{\omega^2}|g_3|, \quad |g_3| \ll 1

Segunda Banda Estrecha (cuando mχ/mϕ>0.5m_\chi/m_\phi > 0.5, correspondiente a ϕϕχχ\phi\phi \to \chi\chi): μmaxq2/16=4ψ02ω4g32,g31\mu_{max} \approx q^2/16 = \frac{4\psi_0^2}{\omega^4}g_3^2, \quad |g_3| \ll 1

3. Efectos de Tamaño Finito y Acoplamiento Crítico

Considerando la localización espacial del oscilón, la tasa de escape de partículas es: ΓescapevχR1R14n2(mχmϕ)2\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}

El acoplamiento crítico (g3)0(g_3)_0 se determina por la condición (μ)Γescape\Re(\mu) \simeq \Gamma_{escape}: (g3)0ω2ψ0R14(mχmϕ)2,mχ0.5mϕ|(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

Esta es la innovación central del artículo: establece una relación cuantitativa entre el acoplamiento crítico y los parámetros geométricos del oscilón.

4. Tratamiento Numérico de Fondo No Uniforme

  • Evolución de campo χ\chi en fondo oscilón fijo ϕ(t,r)=2ψ(r)cos(ωt)\phi(t,r) = 2\psi(r)\cos(\omega t)
  • Uso de método de integración simpléctico de cuarto orden para evolución temporal
  • Adopción de condiciones de frontera de amortiguamiento adiabático para eliminar reflexiones no físicas
  • Extracción de exponentes de Floquet mediante ajuste de Eχe(μ)ωtE_\chi \propto e^{\Re(\mu)\omega t}

Configuración Experimental

Parámetros Numéricos

  • Tamaño de la caja de simulación: rbox=192mϕ1r_{box} = 192 m_\phi^{-1}
  • Número de puntos de malla: 3072
  • Paso de tiempo: Δt=0.01mϕ1\Delta t = 0.01 m_\phi^{-1}
  • Radio de monitoreo: rmax=30mϕ1r_{max} = 30 m_\phi^{-1}
  • Derivadas espaciales: Diferencia central de cuarto orden
  • Integración temporal: Integrador simpléctico de cuarto orden

Adimensionalización

Todos los cálculos utilizan variables adimensionales: x~μ=mϕxμ,ϕ~=λϕmϕ,χ~=λχmϕ\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}m~χ=mχmϕ,g~3=g3mϕλ,g~4=g4λ\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}

El único parámetro libre es mϕ2λ6/λ2=0.8m_\phi^2\lambda_6/\lambda^2 = 0.8 (adoptado en este artículo)

Condiciones Iniciales

Simulación de Fondo Fijo: χ(0,r)=χ0er2/Rχ2,χ˙(0,r)=0\chi(0,r) = \chi_0 e^{-r^2/R_\chi^2}, \quad \dot{\chi}(0,r) = 0 Se toma χ0=0.1\chi_0 = 0.1, Rχ=6R_\chi = 6

Simulación de Evolución de Dos Campos: ϕ(0,r)=2ψ(r),χ(0,r)=χ~0χˉ(r)\phi(0,r) = 2\psi(r), \quad \chi(0,r) = \tilde{\chi}_0 \bar{\chi}(r)ϕ˙(0,r)=χ˙(0,r)=0\dot{\phi}(0,r) = \dot{\chi}(0,r) = 0 Se toma χ~0=0.01\tilde{\chi}_0 = 0.01

Cálculo de Energía

Energía promediada en tiempo (ventana de promedio Tave=20100mϕ1T_{ave} = 20-100 m_\phi^{-1}): Eχ(t)=1Tavett+Tavedt0rmaxdr4πr2[12χ˙2+12(rχ)2+12mχ2χ2]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]

Resultados Experimentales

Resultados Principales

1. Características del Espectro de Floquet (Figura 2)

Para oscilón con ω/mϕ=0.8\omega/m_\phi = 0.8:

  • mχ/mϕ=0.3m_\chi/m_\phi = 0.3: Primera banda estrecha dominante, correspondiente al canal de decaimiento ϕχχ\phi \to \chi\chi
  • mχ/mϕ=0.6m_\chi/m_\phi = 0.6: Segunda banda estrecha dominante (ϕϕχχ\phi\phi \to \chi\chi), porque el decaimiento de una sola partícula está prohibido por conservación de energía

Observaciones clave:

  • Ausencia de bandas de taquiones (Ak0A_k \geq 0 siempre)
  • Simetría respecto a g3=0g_3=0
  • Mayor mχm_\chi desplaza el espectro de Floquet hacia abajo a lo largo del eje kk

2. Tasa de Crecimiento Máxima (Figura 3)

μmaxmax((μ))\mu_{max} \equiv \max(\Re(\mu))

  • Región de acoplamiento pequeño: dependencia lineal μmaxg3\mu_{max} \propto |g_3| cuando mχ/mϕ=0.3m_\chi/m_\phi = 0.3
  • Caso de masa grande: dependencia cuadrática μmaxg32\mu_{max} \propto g_3^2 cuando mχ/mϕ=0.6m_\chi/m_\phi = 0.6
  • Coincidencia exacta con predicciones teóricas (ecuaciones 4.7 y 4.8)

3. Supresión de Tamaño Finito (Figura 4)

Las simulaciones numéricas del fondo oscilón no uniforme muestran:

  • Supresión significativa de la tasa de crecimiento en comparación con el caso uniforme (expansión de región azul estable)
  • Existencia de acoplamiento crítico (g3)0|(g_3)_0|:
    • mχ/mϕ=0.3m_\chi/m_\phi = 0.3: (g3)0/(mϕλ)0.02|(g_3)_0|/(m_\phi\sqrt{\lambda}) \approx 0.02
    • mχ/mϕ=0.6m_\chi/m_\phi = 0.6: (g3)0/(mϕλ)0.1|(g_3)_0|/(m_\phi\sqrt{\lambda}) \approx 0.1
  • Concordancia cuantitativa con estimaciones teóricas de la ecuación (4.14)

4. Análisis Espectral (Figura 5)

El espectro de Fourier del valor central del campo χ\chi muestra:

  • Picos principales ubicados en ωk=k2+mχ2=nω/2\omega_k = \sqrt{k^2+m_\chi^2} = n\omega/2 (resonancia paramétrica de banda estrecha nn)
  • Picos secundarios en (n±1)ω(n\pm1)\omega (efectos de modulación de frecuencia del término de acoplamiento)

Simulaciones de Dinámica de Dos Campos (Figuras 7-8)

Caso 1: Destrucción Completa del Oscilón

Parámetros: mχ/mϕ=0.3m_\chi/m_\phi = 0.3, g3/(mϕλ)=0.3g_3/(m_\phi\sqrt{\lambda}) = 0.3, ω/mϕ=0.8\omega/m_\phi = 0.8 inicial

Resultados:

  • Crecimiento exponencial rápido del campo χ\chi (t<50mϕ1t < 50 m_\phi^{-1})
  • Transferencia rápida de energía del oscilón al campo χ\chi
  • Porque E0oscEdeathE_0^{osc} \approx E_{death}, el oscilón se destruye antes de que termine la resonancia
  • Finalmente Eχ/Etot>0.95E_\chi/E_{tot} > 0.95 (transferencia de energía casi completa)

Caso 2: Terminación Prematura de Resonancia

Parámetros: mχ/mϕ=0.6m_\chi/m_\phi = 0.6, g3/(mϕλ)=0.7g_3/(m_\phi\sqrt{\lambda}) = 0.7, ω/mϕ=0.8\omega/m_\phi = 0.8 inicial

Resultados:

  • Breve período de crecimiento exponencial (t<30mϕ1t < 30 m_\phi^{-1})
  • Campo χ\chi sale de la banda de inestabilidad, el crecimiento se detiene
  • Oscilón residual continúa existiendo, decayendo únicamente de manera perturbativa
  • Eϕ/EtotE_\phi/E_{tot} se estabiliza en aproximadamente 0.4 (energía parcialmente retenida)

Caso 3: Acoplamiento Mixto (Figura 8)

Parámetros: mχ/mϕ=0.3m_\chi/m_\phi = 0.3, g3/(mϕλ)=0.3g_3/(m_\phi\sqrt{\lambda}) = 0.3, g4/λ=0.2g_4/\lambda = 0.2

Hallazgo Clave:

  • g4g_4 no nulo modifica los exponentes de Floquet
  • E0oscE_0^{osc} disminuye, la resonancia termina antes de la destrucción del oscilón
  • Demuestra la universalidad del decaimiento parcial en casos de acoplamiento múltiple

Resumen de Resultados Cuantitativos

Combinación de ParámetrosRelación de Energía CríticaModo de DecaimientoDistribución Final de Energía
mχ=0.3mϕm_\chi=0.3m_\phi, g3=0.3g_3=0.3E0oscEdeathE_0^{osc} \lesssim E_{death}Decaimiento completoEχ/Etot>95%E_\chi/E_{tot} > 95\%
mχ=0.6mϕm_\chi=0.6m_\phi, g3=0.7g_3=0.7E0oscEdeathE_0^{osc} \gtrsim E_{death}Decaimiento parcialEϕ/Etot40%E_\phi/E_{tot} \sim 40\%
Acoplamiento mixtoControlableDecaimiento parcialDepende de g3,g4g_3, g_4

Hallazgos Experimentales

  1. Confirmación de Universalidad: Los acoplamientos de tres y cuatro puntos exhiben patrones de comportamiento cualitativamente idénticos en decaimiento
  2. Escala de Energía Crítica: La tasa de decaimiento del oscilón sigue una ley por segmentos:
undefined