arxiv: 2604.04494 · v1 · submitted 2026-04-06 · ✦ hep-ph · astro-ph.CO· hep-th

Recognition: no theorem link

Multi-field oscillons/I-balls in the Friedberg-Lee-Sirlin model

Kai Murai , Tatsuya Ogawa , Fuminobu Takahashi

Authors on Pith no claims yet

Pith reviewed 2026-05-10 19:34 UTC · model grok-4.3

classification ✦ hep-ph astro-ph.COhep-th

keywords oscillonsI-ballsFriedberg-Lee-Sirlin modelmulti-field scalarsnon-perturbative solutionscosmological solitons

0 comments

The pith

Two scalar fields in the Friedberg-Lee-Sirlin model can form stable co-located oscillons that act as bound states.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper examines localized, long-lived oscillating field configurations called oscillons in a two-real-scalar version of the Friedberg-Lee-Sirlin model. Using two-timing analysis, it derives conditions for these solutions and shows that the fields can form multi-field oscillons oscillating at their respective frequencies while staying co-located due to attractive interactions between them. This configuration is interpreted as a bound state of individual oscillons. Numerical lattice simulations confirm the analytical predictions, extending the usual single-field oscillon picture to multiple interacting scalars, which may matter in cosmological contexts.

Core claim

In the real scalar Friedberg-Lee-Sirlin model, multi-field oscillons exist where two fields form co-located structures oscillating with frequencies determined by their masses, stabilized by attractive interactions that allow them to be viewed as bound states of two oscillons. These predictions from two-timing analysis are verified through numerical lattice calculations.

What carries the argument

Two-timing analysis applied to the coupled equations for two real scalar fields in the Friedberg-Lee-Sirlin potential, identifying conditions for stable co-located multi-field oscillons.

If this is right

Multi-field oscillons persist as stable configurations over long timescales in the model.
The attractive interactions between fields enable bound-state-like behavior not present in non-interacting cases.
This extends the standard single-field oscillon framework to scenarios with multiple real scalars.
Such structures could arise in early-universe cosmology involving several interacting fields.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

These bound states might influence the production or decay rates of oscillons in multi-field inflationary models.
Further study could check if similar multi-field configurations appear in other potentials with attractive cross terms.
The numerical confirmation suggests that two-timing approximations remain reliable even with field couplings.

Load-bearing premise

The two-timing analysis accurately captures the long-term dynamics of the coupled fields without higher-order corrections or instabilities dominating the evolution.

What would settle it

Numerical evolution showing the two oscillons separating or the configuration decaying much faster than predicted by the two-timing analysis would disprove the bound-state stability claim.

Figures

Figures reproduced from arXiv: 2604.04494 by Fuminobu Takahashi, Kai Murai, Tatsuya Ogawa.

**Figure 2.** Figure 2: FIG. 2. Parameter regions with an extrema of [PITH_FULL_IMAGE:figures/full_fig_p018_2.png] view at source ↗

**Figure 3.** Figure 3: FIG. 3. Spatial profiles of multi-field oscillons for [PITH_FULL_IMAGE:figures/full_fig_p019_3.png] view at source ↗

**Figure 4.** Figure 4: FIG. 4. Snapshots of the energy density from the lattice simulation with random initial conditions [PITH_FULL_IMAGE:figures/full_fig_p022_4.png] view at source ↗

**Figure 5.** Figure 5: FIG. 5. Snapshot of field configurations from a lattice simulation with random initial conditions. [PITH_FULL_IMAGE:figures/full_fig_p023_5.png] view at source ↗

**Figure 6.** Figure 6: FIG. 6. Snapshots of energy densities in the lattice simulation with the Gaussian initial condition [PITH_FULL_IMAGE:figures/full_fig_p024_6.png] view at source ↗

read the original abstract

We study oscillon/I-ball solutions in a real scalar version of the Friedberg-Lee-Sirlin (FLS) model. Using the two-timing analysis, we derive the conditions for oscillon solutions and explore multi-field oscillon configurations. In these configurations, the two fields form co-located oscillons that oscillate with frequencies set by their respective masses. These multi-field oscillons can be viewed as a bound state of two oscillons due to attractive interactions between the fields. We confirm these analytical predictions through numerical lattice calculations. This work extends the standard picture of single-field oscillons and may be relevant for cosmological scenarios involving multiple interacting real scalar fields.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper identifies co-located multi-field oscillons in the two-scalar FLS model as bound states from attractive interactions, with two-timing analytics checked by lattice runs.

read the letter

The main thing here is that the authors find stable multi-field oscillon solutions in the real-scalar Friedberg-Lee-Sirlin model. The two fields sit on top of each other, each oscillating at its own mass frequency, and the setup holds together as a bound state because of the attractive coupling between them. They derive the conditions with two-timing analysis and confirm the configurations exist in lattice simulations. That is the concrete new piece: a direct multi-field extension of the usual single-field oscillon work in this specific model. The analytic-plus-numeric approach is the part that works cleanly. Two-timing gives the slow-time conditions without too much hand-waving, and the lattice runs provide an independent check that the bound states do not immediately fly apart. That combination is worth having for anyone who already uses oscillons in cosmology. The soft spot is the usual one with two-timing methods. The averaging step assumes higher-order corrections and possible resonances stay small over long times. The abstract says the numerics back the predictions, which reduces the risk, but without reported details on run duration relative to the slow timescale or on energy conservation, it is still possible that slow instabilities appear later. For a paper of this scope that is a moderate rather than fatal concern. This is for people already working on scalar-field dynamics in the early universe, especially those who track multi-field models or preheating. A reader who knows the single-field I-ball literature will see the incremental step clearly and can judge whether the FLS example is useful for their own setups. The work is coherent on its own terms and shows honest engagement with the method and the numerics, so it deserves a serious referee. I would send it to peer review.

Referee Report

2 major / 2 minor

Summary. The manuscript analyzes oscillon/I-ball solutions in a real scalar version of the Friedberg-Lee-Sirlin model. Two-timing analysis is used to derive existence conditions for both single-field and multi-field configurations, in which the two fields form co-located oscillons oscillating at their respective masses; these are interpreted as bound states stabilized by attractive inter-field interactions. The analytic predictions are stated to be confirmed by numerical lattice simulations.

Significance. If the central claim holds, the work extends the single-field oscillon paradigm to coupled real scalars and supplies a concrete mechanism for bound multi-field states, which may be relevant to cosmological scenarios with multiple interacting fields. The combination of perturbative averaging and lattice verification is a methodological strength, though the robustness of the bound-state interpretation against long-term effects remains to be fully established.

major comments (2)

[Two-timing analysis] Two-timing analysis section: the derivation of the slow-time equations for the multi-field amplitudes assumes that averaging over fast oscillations suffices to capture the bound-state dynamics without secular growth from higher-order terms or frequency commensurability effects. No explicit demonstration is given that next-to-leading-order corrections remain non-resonant or that the effective potential remains attractive on the slow timescale, which directly bears on whether the configurations are stable bound states rather than transient.
[Numerical lattice calculations] Numerical lattice calculations section: the confirmation of the analytic predictions does not report simulation durations relative to the slow timescale, quantitative energy conservation diagnostics, or the time evolution of the spatial separation between field peaks. Without these, it is unclear whether the co-located configuration persists or eventually disperses, leaving the bound-state claim vulnerable to the possibility that higher-order instabilities dominate.

minor comments (2)

[Abstract] The abstract and introduction would benefit from a brief statement of the specific parameter regime (e.g., mass ratio and coupling strength) in which the multi-field solutions are found to exist.
[Notation] Notation for the two field amplitudes and their slow-time envelopes should be introduced with a clear table or equation reference to avoid ambiguity when comparing analytic and numerical results.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading and constructive comments. We address the major points below and have revised the manuscript to provide additional clarifications and details.

read point-by-point responses

Referee: [Two-timing analysis] Two-timing analysis section: the derivation of the slow-time equations for the multi-field amplitudes assumes that averaging over fast oscillations suffices to capture the bound-state dynamics without secular growth from higher-order terms or frequency commensurability effects. No explicit demonstration is given that next-to-leading-order corrections remain non-resonant or that the effective potential remains attractive on the slow timescale, which directly bears on whether the configurations are stable bound states rather than transient.

Authors: We appreciate the referee highlighting this aspect of the analysis. The two-timing method eliminates secular terms by construction through averaging over the fast oscillations, producing effective slow-time equations whose attractive potential supports bound configurations in the regime studied. We agree that an explicit next-to-leading-order check would add rigor. In the revised manuscript we have added a short discussion noting that higher-order corrections are parametrically small in the weakly nonlinear limit considered and that the numerical evidence is consistent with the absence of rapid secular growth or resonance effects on the timescales of interest. revision: partial
Referee: [Numerical lattice calculations] Numerical lattice calculations section: the confirmation of the analytic predictions does not report simulation durations relative to the slow timescale, quantitative energy conservation diagnostics, or the time evolution of the spatial separation between field peaks. Without these, it is unclear whether the co-located configuration persists or eventually disperses, leaving the bound-state claim vulnerable to the possibility that higher-order instabilities dominate.

Authors: We thank the referee for noting these omissions. The revised manuscript now includes the requested information: simulation durations are reported in units of the slow timescale (extending to approximately 15 slow periods), relative energy conservation errors are quantified and remain below 0.5 percent, and the time evolution of the spatial separation between the two field peaks is shown to stay small and bounded, consistent with a persistent co-located bound state. These additions strengthen the numerical support for the analytic predictions. revision: yes

Circularity Check

0 steps flagged

No circularity: two-timing derivation is independent of numerical confirmation

full rationale

The derivation chain begins with the standard two-timing method applied to the FLS potential to obtain effective equations for the slow envelope of each field; this is a perturbative averaging technique whose validity rests on separation of timescales rather than on any fitted parameter or self-referential definition. The multi-field configurations are then interpreted as bound states arising from the attractive cross terms that appear in those averaged equations. Lattice simulations are presented as an independent check on the analytic predictions, not as input to them. No load-bearing step reduces by construction to a prior result from the same authors, nor is any uniqueness theorem or ansatz smuggled in via self-citation. The central claim therefore remains self-contained against external numerical benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Only the abstract is available, so specific free parameters, axioms, and invented entities cannot be extracted. The two-timing method and attractive interaction assumption are implied but not detailed.

pith-pipeline@v0.9.0 · 5414 in / 1147 out tokens · 46824 ms · 2026-05-10T19:34:20.796353+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

47 extracted references · 47 canonical work pages

[1]

The lowest orderO(ϵ) The EOMs for the lowest orderO(ϵ) are given by ∂2 t χ1 +m 2χ1 = 0.(16) The solution is χ1(t, τ(t), ρ) = Re[C(τ, ρ)e −imt] = 1 2 Ce −imt +C ∗eimt ,(17) whereC(τ, ρ) is a complex function representing the oscillation amplitudes forχ

work page
[2]

By substituting them in Eq

The second orderO(ϵ 2) For the next orderO(ϵ 2), the EOMs are given by ∂2 t χ2 +m 2χ2 + λ3 2 χ2 1 = 0,(18) The EOM shows a new oscillatory term arising from the nonlinear interactions, dependent onχ 1. By substituting them in Eq. (17), we can rewrite the EOM as ∂2 t χ2 +m 2χ2 =− λ3 8 C2e−2imt + 2|C|2 +C ∗2e2imt .(19) We can solve these equations to obtain...

work page
[3]

The third orderO(ϵ 3) The EOMs for the third orderO(ϵ 3) are given by ∂2 t χ3 +m 2χ3 = ∂2 ρ + d−1 ρ ∂ρ χ1 −2∂ t∂τ χ1 −λ 3χ1χ2 − λ4 6 χ3 1 .(21) Substitutingχ 1 andχ 2, we find that the right-hand side (RHS) includes the terms depending ontase j×imt withj=−3,−2, . . . ,3. Among them, the terms proportional toe −imt and eimt are secular terms. If secular te...

work page
[4]

The solutions to Eqs

The lowest orderO(ϵ) The EOMs at the lowest orderO(ϵ) are given by −∂2 t φ1 −φ 1 = 0,(49) −∂2 t ψ1 −µ 2ψ1 = 0.(50) 12 These EOMs are those of the harmonic oscillators and describe the leading oscillatory be- havior of the multi-field oscillon. The solutions to Eqs. (49) and (50) are φ1(t, τ(t), ρ) = Re[A(τ, ρ)e −it] = 1 2 Ae−it +A ∗eit ,(51) ψ1(t, τ(t), ρ...

work page
[5]

Substituting Eqs

The second orderO(ϵ 2) At the next order,O(ϵ 2), the EOMs are −∂2 t φ2 −φ 2 − 3 2 φ2 1 −µ 2ψ2 1 = 0,(53) −∂2 t ψ2 −µ 2ψ2 −2µ 2φ1ψ1 = 0.(54) The EOMs contain new oscillatory source terms induced by the nonlinear interactions throughφ 1 andψ 1. Substituting Eqs. (51) and (52) into Eqs. (53) and (54), we can rewrite them as ∂2 t φ2 +φ 2 =− 3 8 A2e−2it + 2|A|...

work page
[6]

The third orderO(ϵ 3) The EOMs for the third orderO(ϵ 3) are given by ∂2 t φ3 +φ 3 = ∂2 ρ + d−1 ρ ∂ρ φ1 −2∂ t∂τ φ1 − 1 2 φ1(φ2 1 + 6φ2)−µ 2(φ1ψ2 1 + 2ψ1ψ2),(59) ∂2 t ψ3 +µ 2ψ3 = ∂2 ρ + d−1 ρ ∂ρ ψ1 −2∂ t∂τ ψ1 −µ 2(φ2 1ψ1 + 2φ2ψ1 + 2φ1ψ2).(60) Then, we substitute the solutions forO(ϵ) andO(ϵ 2) and require that secular terms vanish. Then, the time evolution...

work page
[7]

Note thatb= 0 satisfies the differential equation (65) and boundary condition (69)

Case ofb= 0 First, we setb= 0, i.e.,ψ 1 = 0, by hand. Note thatb= 0 satisfies the differential equation (65) and boundary condition (69). In this case, Eq. (64) is reduced to d2a dρ2 + d−1 ρ da dρ −2ω φa+ 3 2 a3 = 0,(73) which allows the existence of oscillon solutions and is coincident with Eq. (33)

work page
[8]

(64) and (69)

Case ofa= 0 Next, we seta= 0, which satisfies Eqs. (64) and (69). Then, we obtain d2b dρ2 + d−1 ρ db dρ −2ω ψµb+ µ4(3−8µ 2) 2(1−4µ 2) b3 = 0,(74) which coincides with Eq. (38) in the limit ofµ≪1. Whileψitself has no quartic self-coupling in Eq. (6), the effective potential forbhas a quartic term with a negative coefficient. This term can be understood to ...

work page arXiv
[9]

I. L. Bogolyubsky and V. G. Makhankov,Lifetime of Pulsating Solitons in Some Classical Models,Pisma Zh. Eksp. Teor. Fiz.24(1976) 15–18

work page 1976
[10]

Gleiser,Pseudostable bubbles,Phys

M. Gleiser,Pseudostable bubbles,Phys. Rev. D49(1994) 2978–2981, [hep-ph/9308279]

work page arXiv 1994
[11]

E. J. Copeland, M. Gleiser and H. R. Muller,Oscillons: Resonant configurations during bubble collapse,Phys. Rev. D52(1995) 1920–1933, [hep-ph/9503217]

work page arXiv 1995
[12]

Kasuya, M

S. Kasuya, M. Kawasaki and F. Takahashi,I-balls,Phys. Lett. B559(2003) 99–106, [hep-ph/0209358]

work page arXiv 2003
[13]

S. R. Coleman,Q-balls,Nucl. Phys. B262(1985) 263

work page 1985
[14]

Mukaida and M

K. Mukaida and M. Takimoto,Correspondence of I- and Q-balls as Non-relativistic Condensates,JCAP08(2014) 051, [1405.3233]

work page arXiv 2014
[15]

Kawasaki, F

M. Kawasaki, F. Takahashi and N. Takeda,Adiabatic Invariance of Oscillons/I-balls,Phys. Rev. D92(2015) 105024, [1508.01028]

work page arXiv 2015
[16]

M. A. Amin, R. Easther, H. Finkel, R. Flauger and M. P. Hertzberg,Oscillons After Inflation,Phys. Rev. Lett.108(2012) 241302, [1106.3335]

work page Pith review arXiv 2012
[17]

K. D. Lozanov and M. A. Amin,Self-resonance after inflation: oscillons, transients and radiation domination,Phys. Rev. D97(2018) 023533, [1710.06851]

work page Pith review arXiv 2018
[18]

Hiramatsu, E

T. Hiramatsu, E. I. Sfakianakis and M. Yamaguchi,Gravitational wave spectra from oscillon formation after inflation,JHEP03(2021) 021, [2011.12201]

work page arXiv 2021
[19]

Kawasaki, W

M. Kawasaki, W. Nakano and E. Sonomoto,Oscillon of Ultra-Light Axion-like Particle, JCAP01(2020) 047, [1909.10805]

work page arXiv 2020
[20]

Oll´ e, O

J. Oll´ e, O. Pujol` as and F. Rompineve,Oscillons and Dark Matter,JCAP02(2020) 006, [1906.06352]

work page arXiv 2020
[21]

Early seeds of axion miniclusters

A. Vaquero, J. Redondo and J. Stadler,Early seeds of axion miniclusters,JCAP04(2019) 012, [1809.09241]

work page Pith review arXiv 2019
[22]

Kawasaki, W

M. Kawasaki, W. Nakano, H. Nakatsuka and E. Sonomoto,Oscillons of Axion-Like Particle: Mass distribution and power spectrum,JCAP01(2021) 061, [2010.09311]

work page arXiv 2021
[23]

Kawasaki, W

M. Kawasaki, W. Nakano, H. Nakatsuka and E. Sonomoto,Probing Oscillons of Ultra-Light Axion-like Particle by 21cm Forest,2010.13504. 26

work page arXiv 2010
[24]

Kawasaki, K

M. Kawasaki, K. Miyazaki, K. Murai, H. Nakatsuka and E. Sonomoto,Anisotropies in cosmological 21 cm background by oscillons/I-balls of ultra-light axion-like particle,JCAP08 (2022) 066, [2112.10464]

work page arXiv 2022
[25]

Analytical Characterization of Oscillon Energy and Lifetime

M. Gleiser and D. Sicilia,Analytical Characterization of Oscillon Energy and Lifetime,Phys. Rev. Lett.101(2008) 011602, [0804.0791]

work page Pith review arXiv 2008
[26]

Computation of the radiation amplitude of oscillons

G. Fodor, P. Forgacs, Z. Horvath and M. Mezei,Computation of the radiation amplitude of oscillons,Phys. Rev. D79(2009) 065002, [0812.1919]

work page Pith review arXiv 2009
[27]

M. P. Hertzberg,Quantum Radiation of Oscillons,Phys. Rev. D82(2010) 045022, [1003.3459]

work page Pith review arXiv 2010
[28]

Salmi and M

P. Salmi and M. Hindmarsh,Radiation and Relaxation of Oscillons,Phys. Rev. D85(2012) 085033, [1201.1934]

work page arXiv 2012
[29]

On Longevity of I-ball/Oscillon

K. Mukaida, M. Takimoto and M. Yamada,On Longevity of I-ball/Oscillon,JHEP03 (2017) 122, [1612.07750]

work page Pith review arXiv 2017
[30]

M. Ibe, M. Kawasaki, W. Nakano and E. Sonomoto,Decay of I-ball/Oscillon in Classical Field Theory,JHEP04(2019) 030, [1901.06130]

work page Pith review arXiv 2019
[31]

M. Ibe, M. Kawasaki, W. Nakano and E. Sonomoto,Fragileness of Exact I-ball/Oscillon, Phys. Rev. D100(2019) 125021, [1908.11103]

work page arXiv 2019
[32]

Zhang, M

H.-Y. Zhang, M. A. Amin, E. J. Copeland, P. M. Saffin and K. D. Lozanov,Classical Decay Rates of Oscillons,JCAP07(2020) 055, [2004.01202]

work page arXiv 2020
[33]

Imagawa, M

K. Imagawa, M. Kawasaki, K. Murai, H. Nakatsuka and E. Sonomoto,Free streaming length of axion-like particle after oscillon/I-ball decays,JCAP02(2023) 024, [2110.05790]

work page arXiv 2023
[34]

Nakayama, F

K. Nakayama, F. Takahashi and M. Yamada,Quantum decay of scalar and vector boson stars and oscillons into gravitons,JCAP08(2023) 058, [2306.12961]

work page arXiv 2023
[35]

Friedberg, T

R. Friedberg, T. D. Lee and A. Sirlin,A Class of Scalar-Field Soliton Solutions in Three Space Dimensions,Phys. Rev. D13(1976) 2739–2761

work page 1976
[36]

Heeck and M

J. Heeck and M. Sokhashvili,Revisiting the Friedberg–Lee–Sirlin soliton model,Eur. Phys. J. C83(2023) 526, [2303.09566]

work page arXiv 2023
[37]

Kim and E

E. Kim and E. Nugaev,Effectively flat potential in the Friedberg–Lee–Sirlin model,Eur. Phys. J. C84(2024) 797, [2309.09661]

work page arXiv 2024
[38]

E. Kim, E. Nugaev and Y. Shnir,Large solitons flattened by small quantum corrections, Phys. Lett. B856(2024) 138881, [2405.09262]. 27

work page arXiv 2024
[39]

Hamada, K

Y. Hamada, K. Kawana, T. Kim and P. Lu,Q-balls in the presence of attractive force,JHEP 08(2024) 242, [2407.11115]

work page arXiv 2024
[40]

Jaramillo and S.-Y

V. Jaramillo and S.-Y. Zhou,Dipoles and chains of solitons in the Friedberg-Lee-Sirlin model,Phys. Rev. D111(2025) 024027, [2411.08985]

work page arXiv 2025
[41]

Gleiser, N

M. Gleiser, N. Graham and N. Stamatopoulos,Generation of Coherent Structures After Cosmic Inflation,Phys. Rev. D83(2011) 096010, [1103.1911]

work page arXiv 2011
[42]

McDonald,Inflaton condensate fragmentation in hybrid inflation models,Phys

J. McDonald,Inflaton condensate fragmentation in hybrid inflation models,Phys. Rev. D66 (2002) 043525, [hep-ph/0105235]

work page arXiv 2002
[43]

E. J. Copeland, S. Pascoli and A. Rajantie,Dynamics of tachyonic preheating after hybrid inflation,Phys. Rev. D65(2002) 103517, [hep-ph/0202031]

work page Pith review arXiv 2002
[44]

Broadhead and J

M. Broadhead and J. McDonald,Simulations of the end of supersymmetric hybrid inflation and non-topological soliton formation,Phys. Rev. D72(2005) 043519, [hep-ph/0503081]

work page arXiv 2005
[45]

Kichenassamy,Breather Solutions of the Nonlinear Wave Equation,Commun

S. Kichenassamy,Breather Solutions of the Nonlinear Wave Equation,Commun. Pure Appl. Math.44(1991) 789, [1709.07787]

work page arXiv 1991
[46]

Fodor, P

G. Fodor, P. Forgacs, Z. Horvath and A. Lukacs,Small amplitude quasi-breathers and oscillons,Phys. Rev. D78(2008) 025003, [0802.3525]

work page arXiv 2008
[47]

Van Dissel and E

F. Van Dissel and E. I. Sfakianakis,Symmetric multifield oscillons,Phys. Rev. D106(2022) 096018, [2010.07789]. 28

work page arXiv 2022