{"id":"1f56e27e-a9f6-4961-a7f7-718ab9c9237f","arxiv_id":"1908.02296","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Neumann boundary conditions in anti-de Sitter space admit non-collapsing, multi-mode scalar solutions with arbitrarily small amplitude despite a fully resonant spectrum.","lead":"A numerical study finds that a scalar field in a boxed-in spacetime called anti-de Sitter space, with Neumann boundary conditions, can form non-collapsing patterns made of many waves even at very small amplitude, despite a resonant frequency spectrum. This challenges the common assumption that only single-wave patterns are stable in such resonant systems.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 'arbitrarily close to AdS' claim is inferred from a trend ending at ε=0.001, not from a convergent sequence; continuation to smaller ε and a grid-convergence check are needed.","rationale":"The paper's central claim is that Neumann boundary conditions admit non-collapsing, strongly multi-mode solutions arbitrarily close to AdS despite a resonant spectrum. The numerical construction is credible: pseudo-spectral BVP, Hamiltonian and energy-conservation checks, a time evolution at ε=0.06, and a Dirichlet limit that behaves as expected. The load-bearing step is the passage from computed solutions at ε≥0.001 to 'arbitrarily close'. Nothing in the paper provides a convergent sequence, a continuation argument, or a rigorous normal-form construction; Fig. 5's constant decay rate is a trend over roughly one order of magnitude. The linearized operator about AdS is degenerate at ε=0, so the Newton solver could fail or find spurious branches at smaller amplitudes, and the power-law tail in the normal-mode projection complicates the exponential decay-rate fit. The reader's weakest assumption identifies exactly this point, and I agree with it. The proposed continuation to smaller ε with resolution checks would settle whether the branch persists. This does not refute the result; it supports keeping the verdict conditional rather than accepting the 'arbitrarily close' phrasing as established.","tokens_in":9607,"tokens_out":20894,"duration_ms":244869,"concrete_test":"Compute the κ=0 double oscillator at ε = 3e-4, 1e-4 and 3e-5 with at least double the Fourier and Legendre resolutions used for ε=0.001, monitoring the PDE residual and the Hamiltonian constraint. Re-measure the exponential decay rate of |a_n| from the same fitting procedure. If the decay rate remains constant within fit error and residuals decrease to machine precision as resolution increases, the extrapolation to arbitrarily small ε is supported. If the solver fails to converge or the decay rate moves toward the two-mode value, the 'arbitrarily close' claim is not established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim in the abstract is that Neumann (κ=0) double oscillators exist arbitrarily close to AdS and are strongly multi-mode. The evidence is Fig. 4 at ε=0.001 and the plateau of the exponential decay rate in Fig. 5 for ε down to 0.001. That is an extrapolation, not a demonstrated limit. The paper reports no convergence study or residual tolerances for the boundary-value problem at these small amplitudes, and the normal-mode amplitudes used for the decay-rate fit develop a power-law tail at high n because the non-linear solutions are not even about the AdS boundary, which complicates the exponential fit. Near ε=0 the linearized problem about AdS has a large kernel (every normal-mode combination is a linear solution), so the Newton-Raphson branch could terminate or become numerical noise before reaching arbitrarily small ε. If the family does not extend, the claim reduces to existence at small but finite amplitude, which is weaker than 'arbitrarily close to AdS'. The discussion itself concedes that the origin of the Neumann solutions is not yet understood and that the perturbative cancellation of secular terms is only a hint.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper constructs quasi-periodic 'double oscillator' solutions for a massive complex scalar field (m^2 = -2) in spherically symmetric global AdS, with Robin boundary conditions parameterized by κ that interpolate between resonant Dirichlet (κ = 1) and resonant Neumann (κ = 0) boundary conditions. The central findings are: (i) two-mode equal-amplitude double oscillators exist for generic non-resonant Robin boundary conditions, and their maximum energy vanishes as κ → 1, consistent with the expected instability of two-mode data in the Dirichlet case; and (ii) surprisingly, at the Neumann endpoint κ = 0 the family does not vanish, and the solutions remain strongly multi-mode as the amplitude ε is decreased down to ε = 0.001, despite the fully resonant spectrum. The authors support the existence claim with a Newton-Raphson pseudo-spectral solver, energy and constraint checks at the 10^{-10} level, and a numerical time evolution at one amplitude (ε = 0.06) that shows no collapse up to t ≈ 4630. They also give a perturbative discussion suggesting that the Neumann solutions require seeding an infinite number of modes at leading order.","tokens_in":9791,"tokens_out":5183,"duration_ms":58594,"significance":"If the central claim is correct, the paper provides the first example in AdS of non-collapsing, strongly multi-mode initial data in a system with a fully resonant spectrum, contradicting the emerging picture that islands of stability in resonant cases are single-mode dominated. It also gives a quantitative map of how the space of double oscillators depends on the boundary-condition parameter κ, connecting non-resonant stability to the existence of multi-mode islands. The numerical evidence is substantial: the solver uses a consistent spectral discretization, the Hamiltonian constraint and energy are checked to 10^{-10}, and the paper demonstrates robustness of the existence curve to grid size and parametrization. The main weakness is that the 'arbitrarily close to AdS' statement rests on extrapolation of a trend ending at ε = 0.001, not on a demonstrated limit or a convergence study in ε.","major_comments":[{"comment":"The claim in the abstract and discussion that Neumann double oscillators exist 'arbitrarily close to AdS' is inferred from the plateau of the exponential decay rate in Fig. 5 and the spectrum in Fig. 4 down to ε = 0.001. The paper does not report solutions at smaller ε, does not state residual tolerances for the Newton-Raphson solver, and does not provide a convergence study in grid resolution for the small-ε branch. Since the linearized problem at ε = 0 has an infinite-dimensional kernel (every normal-mode combination is a solution), the numerical branch could terminate or become dominated by numerical noise before reaching arbitrarily small amplitude. Please provide either a convergent sequence of solutions with ε → 0 (with stated residuals and grid sizes) or explicitly weaken the central claim to existence for small but finite ε.","section":"Double Oscillators, Figures 4 and 5"},{"comment":"The statement that the Neumann solutions are 'not single-mode dominated' relies on the exponential decay rate of the normal-mode amplitudes |a_n|, but the n-range used for the linear-regression fit in Fig. 5 is not stated. The text itself notes that at very high n there is a power-law tail because the normal-mode basis is even about the boundary while the nonlinear solutions are only Neumann. If the fitting window includes part of this tail, the 'plateau' in the Neumann decay rate could be an artifact of the fit rather than a physical property of the solution. Please specify the fitted n-interval, show the decay rate as a function of the fit window, and demonstrate that the plateau is robust to the choice of window.","section":"Figures 4 and 5, normal-mode decay-rate fit"},{"comment":"The nonlinear-stability evidence is presented for a single evolution at ε = 0.06 with energy E = 0.24, evolved to t ≈ ε^{-3} ≈ 4630. The central claim of the paper concerns the limit ε → 0, so the absence of collapse at one moderate amplitude does not establish that the small-ε family is stable on the instability timescale. Moreover, the Discussion states the solutions are 'stable (until at least t ∼ 1/ε^2)', which is inconsistent with the actual evolution time reported (ε^{-3}); please clarify the intended stability timescale and, if the stability claim is to be part of the abstract-level conclusions, provide evolutions at smaller ε or soften the wording to an existence statement with a single finite-amplitude stability check.","section":"Numerical evolution"}],"minor_comments":[{"comment":"There are several typos: 'investigate the affect' should be 'investigate the effect', and the Figure 1 caption has 'mulit-oscillators' instead of 'multi-oscillators'.","section":"Introduction and Figure 1"},{"comment":"The text says the functions are 'periodic in time with period ω2', but since the Fourier expansion uses cos(k ω2 t), the quantity ω2 is an angular frequency and the period is 2π/ω2. Please correct this wording.","section":"Double Oscillators"},{"comment":"The y-axis label in Figure 2 appears garbled: 'ϵmax( e3 e1 )=0.1' should be replaced by a clear expression such as the maximum ε for which |(f_3^{(3)}, ê_3)/(f_3^{(1)}, ê_1)| < 0.1.","section":"Figure 2"},{"comment":"The phrase 'positive definite frequencies' should be 'positive frequencies'; positive definiteness is a property of operators or bilinear forms, not of frequencies.","section":"Perturbative analysis"},{"comment":"The paper would benefit from stating the grid sizes, number of Fourier modes, and residual/error tolerances used for the boundary-value problem, and from completing references [1], [2], and [49], which currently lack full bibliographic data.","section":"Numerical methods and references"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The punchline: this paper has a genuinely new numerical result—Neumann boundary conditions in AdS admit non-collapsing, strongly multi-mode scalar configurations that sit inside the islands of stability despite a fully resonant spectrum. If the result holds, it breaks the single-mode-dominated picture for resonant systems. I think it likely holds, but the 'arbitrarily close to AdS' part of the abstract is not as solid as the rest.\n\nWhat is actually new: prior work in the Dirichlet case only found single-mode-dominated double oscillators; here the Neumann double oscillators are multi-mode dominated, with the mode-amplitude decay rate approaching a constant as epsilon goes down to 0.001. That is a clear counterexample to the expectation. The Robin non-resonant case also behaves as expected, with the double oscillators disappearing as you approach Dirichlet. The numerics are credible: Newton-Raphson with pseudo-spectral discretization, energy and constraint checks at 10^-10, and a time evolution at epsilon=0.06 that stays periodic to t ~ 4630. The method is inherited from their prior work, and they cite it, which is appropriate. The perturbative discussion is honest, including the admission that the origin of the Neumann solutions is not yet understood.\n\nThe soft spots are in the extrapolation to zero amplitude. The smallest epsilon is 0.001; the claim 'arbitrarily close to AdS' is inferred from a trend in Fig. 5, not from a convergent sequence. No grid-convergence or residual tolerances are reported at these small amplitudes, and the normal-mode fit has a known power-law tail at high n because the nonlinear solutions are not even about the AdS boundary. Near epsilon=0, the linearized problem has a large kernel (every combination of normal modes is a solution), so the Newton branch could in principle terminate or become numerical noise before reaching zero. The authors themselves concede that the origin is not understood and that the perturbative secular-term cancellation is only a hint. Also, the time evolution at epsilon=0.06 (energy 0.24) is not small compared to the claimed small-epsilon limit, so nonlinear stability at small amplitudes is inferred, not directly tested.\n\nNone of this kills the paper. Even the finite-amplitude statement is important, and the multi-mode character is established at the smallest computed amplitude. But the abstract should either be softened or supported by a continuation study.\n\nThis paper is for the AdS stability and holography community. It deserves a serious referee; I'd send it to review and ask the authors to add a small-epsilon convergence study or qualify the claim. I would cite it.","headline":"A credible numerical discovery—Neumann AdS admits strongly multi-mode non-collapsing solutions near AdS—but the 'arbitrarily close' claim rests on extrapolation from epsilon=0.001, so treat that part as suggestive.","tokens_in":10292,"tokens_out":2995,"would_cite":true,"duration_ms":28134,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Fully resonant Neumann boundary conditions still admit stable, strongly multi-mode oscillators arbitrarily close to AdS, contradicting the single-mode-dominated picture of the islands of stability.","keywords":["anti-de Sitter instability","islands of stability","multi-oscillators","double-trace deformations","Robin boundary conditions","resonant spectrum","quasi-periodic solutions","black hole formation"],"falsifier":"Repeat the construction at $\\varepsilon=10^{-4}$ or smaller and check whether the branch persists with a nonzero exponential mode-decay rate; or evolve the $\\kappa=0$ double oscillator with $\\varepsilon=0.001$ for times much longer than $\\varepsilon^{-3}$ and watch whether $\\varphi\\varphi^*$ at the boundary remains periodic or develops drift toward horizon formation. Failure of either test would disqualify the 'arbitrarily close to AdS' claim.","tokens_in":9395,"feed_emoji":"🌊","tokens_out":7244,"duration_ms":70414,"temperature":0.7,"pith_summary":"The paper asks what role a resonant linear spectrum plays in the nonlinear instability of anti-de Sitter space (AdS), and whether breaking the resonance restores stability. It studies quasi-periodic, non-collapsing 'double-oscillator' solutions of Einstein–Klein–Gordon theory with a massive complex scalar field ($m^2=-2$), using Robin boundary conditions that interpolate between Dirichlet and Neumann. The authors find two-mode, equal-amplitude double oscillators of arbitrarily small energy for non-resonant Robin spectra, and find that these solutions vanish as the boundary condition approaches the fully resonant Dirichlet case. The surprise is the Neumann endpoint: although its spectrum is also fully resonant, strongly multi-mode non-collapsing solutions persist arbitrarily close to AdS and are not captured by two-mode perturbation theory. If correct, the common picture that islands of stability in resonant systems are single-mode dominated is incomplete.","feed_headline":"Fully resonant AdS still hosts stable multi-mode waves","feed_subtitle":"A resonant spectrum does not force collapse: two-mode oscillators persist at tiny amplitudes under Neumann conditions.","key_machinery":"The central object is the double-oscillator, a quasi-periodic solution oscillating on two frequencies, constructed by imposing a time-periodic Fourier ansatz and solving the Einstein–Klein–Gordon equations as a boundary value problem with a Newton–Raphson method and pseudo-spectral discretization. The boundary-condition family is the Robin condition $\\sin(\\pi\\kappa/2)\\varphi_1-\\cos(\\pi\\kappa/2)\\varphi_2=0$, with $\\kappa=1$ giving the fully resonant Dirichlet spectrum (frequencies $\\omega_n=2n$) and $\\kappa=0$ giving the fully resonant Neumann spectrum (frequencies $\\omega_n=2n+1$). Resonance is diagnosed through frequency quadruples satisfying $\\Delta\\omega_J=\\omega_{j_2}+\\omega_{j_3}-\\omega_{j_4}-\\omega_{j_1}=0$. The key diagnostic is the projection of a solution onto the normal modes; the exponential decay rate of those mode amplitudes distinguishes two-mode-dominated from strongly multi-mode data.","core_discovery":"The paper claims that the fully resonant Neumann spectrum admits islands of stability containing strongly multi-mode data arbitrarily close to AdS, and that these solutions are not single-mode dominated. Concretely, for a complex scalar field with $m^2=-2$, the authors construct double-oscillators, quasi-periodic solutions on two frequencies, with equal amplitude imposed on the first two normal modes. For Neumann boundary conditions ($\\kappa=0$) these solutions are found down to amplitude $\\varepsilon=0.001$, and their mode spectra show an exponential fall-off whose decay rate approaches a nonzero constant as $\\varepsilon\\to 0$, so the small-amplitude limit does not resemble two-mode data. A numerical evolution of a representative solution with $\\varepsilon=0.06$ shows periodic $\\varphi\\varphi^*$ at the boundary and no horizon formation up to $t\\approx\\varepsilon^{-3}$, which the authors take as evidence of nonlinear stability. This contradicts the expectation that, in a fully resonant system, stable solutions must be single-mode dominated.","pith_inferences":["The Neumann result suggests that it is the structure of the secular terms, not merely the presence of a resonance, that decides whether multi-mode data collapse; a resonant spectrum may be necessary for instability but is not sufficient.","If these bulk solutions have a CFT dual, a double-trace-deformed boundary theory would admit long-lived, non-thermal states at arbitrarily small energy, which is relevant for how double-trace deformations affect thermalization.","A natural testable extension is to scan other scalar masses and boundary-condition families; the flat-space Dirichlet-box analogue considered in the paper does not reproduce the Neumann plateau, suggesting the AdS boundary structure is essential.","The observation that extra secular terms shrink when more modes are added at $O(\\varepsilon)$ hints at a resummation or non-perturbative description in which the Neumann double oscillator is a coherent multi-mode state; constructing such a description would make the claim testable beyond numerics."],"forward_implications":["For non-resonant Robin spectra, islands of stability include strongly multi-mode data of arbitrarily small energy, supporting the idea that breaking resonance restores stability.","As the boundary condition approaches the fully resonant Dirichlet case, the maximum energy of two-mode, equal-amplitude double oscillators tends to zero, consistent with such data collapsing.","In the Neumann case, despite a fully resonant spectrum, strongly multi-mode double oscillators exist arbitrarily close to AdS and are not single-mode dominated.","These Neumann double oscillators are not described by two-mode perturbation theory: their higher-mode amplitudes decay exponentially with a rate that stays nonzero as $\\varepsilon\\to 0$, so a perturbative description would require infinitely many modes at leading order.","Numerical evolution shows a Neumann double oscillator with $\\varepsilon=0.06$ remains periodic and non-collapsing up to $t\\approx\\varepsilon^{-3}$, evidence that the solutions are nonlinearly stable at least on that timescale."],"supporting_citations":[{"why":"Supplies the original numerical instability result and the $\\varepsilon^{-2}$ collapse timescale that motivates the search for non-collapsing data.","marker":"[3]"},{"why":"Supplies the analytical argument that a resonant spectrum may be a necessary condition for the AdS instability.","marker":"[6]"},{"why":"Supports the expectation that non-resonant spectra restore stability, which the Robin-case results confirm.","marker":"[36]"},{"why":"Introduces multi-oscillators as a tool to chart islands of stability and reports single-mode-dominated solutions in the massless Dirichlet case, the baseline the Neumann result is compared against.","marker":"[42]"},{"why":"Maps Robin boundary conditions to double-trace deformations, grounding the boundary-condition family in the dual CFT.","marker":"[48]"},{"why":"Defines the general multi-oscillator construction from which the double-oscillator solutions are taken.","marker":"[49]"}],"fun_headline_variants":["Neumann AdS defies resonant instability with stable two-mode waves","Stable multi-mode waves survive full resonance in AdS","Resonant AdS not doomed: Neumann admits stable multi-mode islands","Beyond single-mode: AdS stability islands under full resonance","Contrary to expectations, resonant AdS hosts stable multi-mode data"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the numerical family of Neumann double oscillators continues to exist as the amplitude $\\varepsilon$ tends to zero; the smallest computed case is $\\varepsilon=0.001$, so if the branch stops at some positive amplitude, the claim weakens from 'arbitrarily close to AdS' to existence at finite small amplitudes.","fun_headline_variants_meta":{"raw":{"variants":["Neumann AdS defies resonant instability with stable two-mode waves","Stable multi-mode waves survive full resonance in AdS","Resonant AdS not doomed: Neumann admits stable multi-mode islands","Beyond single-mode: AdS stability islands under full resonance","Contrary to expectations, resonant AdS hosts stable multi-mode data"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00022,"raw_usage":{"total_tokens":1412,"prompt_tokens":873,"completion_tokens":539,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":489,"completion_tokens_details":{"reasoning_tokens":451}},"tokens_in":489,"tokens_out":539,"duration_ms":5575,"temperature":1.0,"reasoning_tokens":451,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:48:16.040709+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Repeat the construction at $\\varepsilon=10^{-4}$ or smaller and check whether the branch persists with a nonzero exponential mode-decay rate; or evolve the $\\kappa=0$ double oscillator with $\\varepsilon=0.001$ for times much longer than $\\varepsilon^{-3}$ and watch whether $\\varphi\\varphi^*$ at the boundary remains periodic or develops drift toward horizon formation. Failure of either test would disqualify the 'arbitrarily close to AdS' claim.","supporting_citations":[{"cited_title":"Necessary conditions for an AdS-type instability","cited_arxiv_id":"1509.00232","evidence_quote":"Supports the expectation that non-resonant spectra restore stability, which the Robin-case results confirm."},{"cited_title":"Charting the AdS Islands of Stability with Multi-oscillators?","cited_arxiv_id":"1803.02830","evidence_quote":"Introduces multi-oscillators as a tool to chart islands of stability and reports single-mode-dominated solutions in the massless Dirichlet case, the baseline the Neumann result is compared against."}],"review_version":1}