{"id":"c66a4806-e9a4-463a-988a-e028f304996e","arxiv_id":"2412.14259","paper_version":4,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"In a consistent Einstein-Maxwell-scalar theory, black hole ringdown echoes seen in linear theory survive fully nonlinear radial evolution.","lead":"Black holes in a theory with scalar and electric hair can bounce incoming waves between two potential barriers, creating repeated echo pulses in their ringdown. This paper reports a fully nonlinear, consistent simulation showing that those echoes survive even when the perturbation is violent.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Nonlinear-echo claim rests on a single run whose remnant identification is tested only by a 4-mode fit; the extracted late-time signal shows no measured echo delay and the claimed final-state QNM fit may not uniquely certify echoes.","rationale":"The reader's verdict (CONDITIONAL) already captures the main structural weaknesses of the nonlinear-echo claim, and my stress test does not contradict the numerical data or the linear-level derivations. My specific concern is narrower and more concrete than the reader's weakest_assumption (which focused on strong hyperbolicity). I judge the well-posedness worry to be real but secondary: the EMS system is a GR-like two-derivative theory, the paper explicitly cites the expectation of strong hyperbolicity, and the 1+1 metric ansatz is fixed-gauge, so a loss of well-posedness would likely manifest as code crash or constraint blowup, which the paper partially monitors. The more load-bearing issue for the central claim is the evidential weight of a single tuned nonlinear run: the echoes are identified only visually in a semi-log plot of ∂tϕ, no echo delay is measured, the remnant identification is not independently tested, and the QNM fit is not accompanied by residuals or uncertainty quantification. These are standard checks for claims of new phenomena, and the paper does not provide them. Still, none of these gaps amounts to demonstrated internal inconsistency: the derivations are detailed, the linear and nonlinear codes are cross-checked, and the convergence tests in App. B are a real, if partial, form of support. A conditional acceptance requiring the additional checks is therefore the appropriate disposition, and my recommendation is UNCHANGED relative to the reader's verdict. I set agreement_with_reader to 'partial' because I agree with the overall CONDITIONAL verdict but identify a different load-bearing concern than the reader's well-posedness emphasis.","tokens_in":26769,"tokens_out":1903,"duration_ms":15678,"concrete_test":"Re-run the nonlinear case of Fig. 10 with two higher resolutions (e.g., Δr0/2 and Δr0/4) and, for each, extract ∂tϕ at R = 100 M0 and compute the time interval between successive echo peaks in the late-time signal; verify that Δt_echo is resolution-independent and matches the cavity round-trip time 2∫ dr*/v_g derived from the final background potential Vϕ. Additionally, perform a second independent fit of the same late-time signal without assuming the remnant is the static scalarized BH: fit the raw ∂tϕ(t, R) starting at t > 400 M0 with a free sum of damped sinusoids and check whether the best-fit frequencies agree with the QNMs of the claimed remnant at the 1% level, and report reduced chi-squared and fit residuals.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim, 'first example of echoes appearing in a consistent theory beyond a linearized analysis,' rests almost entirely on the single nonlinear run of Sec. V.C (Fig. 10, α = 0.6, q_initial = 1.0194, A = 0.03 M0). The paper reports that the remnant horizon area grows by ≈32, the final configuration is assumed to be the static scalarized BH with horizon radius Rh and charge Qh, and the late-time ∂tϕ at R = 100 M0 is fit by the first four QNMs of that assumed final BH. Three load-bearing gaps follow. (i) Echoes are asserted visually from scalloped structure in a semi-log plot of a quantity that is a time derivative of ϕ; no echo delay Δt_echo is measured or compared to the cavity round-trip time, so the 'echo' identification is not quantitatively established. (ii) The remnant identification is circular to a degree: the QNM template is built from the same static-solution family used to infer the remnant, and the fit starts at t > 200 M0 with only four damped sinusoids and no reported residuals or parameter uncertainties; a different final state (e.g., an oscillating or partially charged remnant) might fit equally well. (iii) The code crash noted in App. B for the nonlinear run at low resolution is acknowledged, but no convergence test of the actual echo feature itself is shown; only global constraint violations at t = 812 M0 are presented, and the one-sigma significance of the echoed bumps against numerical noise is not quantified.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies linear and nonlinear perturbations of spherically symmetric black holes in Einstein-Maxwell-scalar (EMS) theory with coupling F[phi] = exp(alpha phi^2). It derives the linearized radial and non-radial perturbation equations, computes quasinormal-mode frequencies with a frequency-domain shooting method, and evolves the linear perturbations in the time domain. In a parameter region where the scalarized black holes have a stable photon sphere and multipeaked effective potentials, the linear response shows repeated echo-like structures. The paper then presents fully nonlinear 1+1 simulations of scalar wave packets falling onto scalarized black holes. In the linear-amplitude regime the nonlinear code reproduces the echoes found by the linear code; in one high-amplitude run, in which the horizon area grows by a factor of about 32, it claims that echoes persist and that the late-time signal is well approximated by the first four quasinormal modes of the presumed final scalarized black hole. The abstract states that this is the first example of echoes appearing in a consistent theory beyond a linearized analysis.","tokens_in":27072,"tokens_out":3773,"duration_ms":36330,"significance":"If the nonlinear echo claim is upheld, the paper would be a valuable step beyond the large body of work on echoes in linearized and often ad hoc models: it would show the phenomenon in a concrete field-theory setting that has a variational principle and is amenable to full nonlinear evolution. The linear part of the paper is solid and provides useful material in itself: the perturbative equations are obtained from the action, the RN limit is recovered, the frequency-domain and time-domain linear codes agree, and the linear echo patterns are clearly exhibited. The nonlinear claim, however, is the main advertised result and currently rests on a single simulation whose echo identification is largely visual and whose remnant identification is tested only by a four-mode fit. These points need to be strengthened before the 'first beyond-linearized example' claim can be accepted.","major_comments":[{"comment":"The central claim that echoes survive in the nonlinear regime is based on a single run, and the echo identification is not quantitative. No echo delay is measured or compared with the cavity round-trip time implied by the effective potential, and no criterion is given that distinguishes the scalloped structure in the semi-log plot from a generic smooth late-time decay. I ask the authors to define the echo feature quantitatively (e.g., peak times and expected delay from the potential profile), to report the residuals of the four-mode fit, and ideally to repeat the nonlinear run with at least one different set of parameters to show that the feature is robust.","section":"Sec. V.C, Fig. 10"},{"comment":"The remnant identification is partly circular: the final state is assumed to be the static scalarized black hole with the measured horizon radius and charge, and the fit uses quasinormal frequencies of the same solution family. This is a consistency check but not an independent identification of the remnant. I request a comparison with at least one alternative model (for example, the RN quasinormal modes of the same mass and charge, or a non-static remnant) and a report of the fit residuals and parameter uncertainties, so that the reader can judge whether the data actually prefer the scalarized final state.","section":"Sec. V.C, Fig. 10"},{"comment":"The convergence tests in Appendix B establish fourth-order scaling of the constraint violations for the nonlinear run at t = 812 M0, but they do not test convergence of the extracted observable ∂tϕ at R = 100 M0 or of the echo bumps themselves. The lowest-resolution nonlinear run crashes, and no estimate of numerical noise in Fig. 10 is provided, so the significance of the echoed bumps relative to numerical error is unquantified. I ask for a resolution study of the echo feature itself and an error estimate for the fitted amplitudes and phases.","section":"Appendix B"},{"comment":"The framing of the nonlinear simulation as a 'consistent theory' rests on the well-posedness of the EMS system, but the paper only states an expectation of strong hyperbolicity for reasonable F[phi]. Since the nonlinear runs use very small horizon radii, a large amplitude wave packet, and a regime in which the final configuration is close to the critical charge, the authors should either provide a reference or a short numerical/analytic check of strong hyperbolicity in the regime they simulate, or explicitly qualify the claim so that it is not stronger than the evidence.","section":"Sec. I, bullet i; Sec. V.A"}],"minor_comments":[{"comment":"The heading 'Frequecy domain method' contains a typo and should read 'Frequency domain method'.","section":"Sec. IV.A"},{"comment":"Reference [65] contains the typo 'Einstein-Mawxell-dilaton theory'; it should read 'Einstein-Maxwell-dilaton theory'.","section":"Ref. [65]"},{"comment":"The y-axis labels in the text rendering of Figs. 9 and 10 appear garbled as 't (t, R = R) M0'; they should be typeset as ∂tϕ(t, R = 100 M0) with appropriate units.","section":"Figs. 9 and 10"},{"comment":"There is a small inconsistency in the reported start time of the nonlinear fit: the text in Fig. 10 says 't > 200 M0', while the discussion in Sec. V.C and the caption should be harmonized with the exact value used.","section":"Sec. V.C and Fig. 10"},{"comment":"The paper should state explicitly that in the fit of Eq. (90) only the amplitudes A_n and phases phi_n are fitted, while the frequencies are taken from the frequency-domain calculation of Sec. IV.C; this is clear from the text but an explicit sentence would avoid any impression that the frequencies are being tuned.","section":"Eq. (90)"}],"recommendation":"major_revision","confidential_remarks":"The linear perturbative analysis is careful and the code validation is convincing, but the headline nonlinear echo claim is currently supported by a single run with only a four-mode fit and no quantitative echo diagnostics. The requested additions are substantial but well within the scope of a revision: a quantitative echo definition, residual and uncertainty reporting, a convergence test of the extracted signal, and a well-posedness statement. I therefore recommend major revision rather than rejection; if the nonlinear evidence is not strengthened, the abstract's 'first example' claim should be downgraded."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth a serious referee. The new piece is the fully nonlinear 1+1 evolution showing echo-like structure survives in a consistent field theory. The supporting work is done well: the perturbative equations are derived from the action, the frequency-domain QNM code agrees with the RN limit, the time-domain code reproduces the same linear response, and convergence tests show clean fourth-order scaling. The linear sector alone, including the axial and polar systems, is a useful extension of the earlier scalarized-BH literature.\n\nSoft spots, in proportion. First, well-posedness is asserted as an expectation, not proved. That is a minor gap; the action is GR with a dielectric Maxwell term, so strong hyperbolicity is the reasonable default, but a sentence with a reference or a short argument would close it. Second, the nonlinear echo evidence is a single run, interpreted through a late-time fit to four QNMs of the assumed final state. The fit is shown but residuals and parameter uncertainties are not, and no echo delay is measured or compared against the cavity round-trip time. I do not think this is circular — the QNM frequencies come from a boundary-value problem, not from the time signal — but the remnant identification rests on that fit alone. Third, the 'first example beyond linearized analysis' claim is not established inside the paper. The cited nonlinear boson-star and ECO simulations may or may not contain echoes; the authors should say why those cases do not count, or soften the claim. Fourth, the low-resolution run crashed and the convergence test is performed at t=812 M0 rather than at the final time. This is minor because the crash occurs after the echo feature appears, but a convergence test of the echo feature itself would be more convincing.\n\nBottom line: this is a paper for the ringdown and beyond-GR community. The linear results are citable on their own; the nonlinear run is the novelty that makes it interesting. I would send it to peer review and ask for a revision that computes the echo delay, reports fit residuals and uncertainties, and tightens the literature claim. Conditional acceptance.","headline":"Worth refereeing: the linear work is careful and the nonlinear echo run is a genuine first, but the 'first' claim and the single-run identification need tightening.","tokens_in":27611,"tokens_out":4636,"would_cite":true,"duration_ms":43486,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["04.70.-s","04.30.-w"],"model":"deepseek-v4-flash","headline":"Black hole echoes survive fully nonlinear collapse in a consistent theory.","keywords":["black hole echoes","Einstein-Maxwell-scalar theory","quasinormal modes","scalarization","ringdown spectroscopy","photon sphere","nonlinear numerical relativity"],"falsifier":"Repeat the $A=0.03\\,M_0$ collapse with an independent, manifestly strongly hyperbolic formulation of the same theory; if the late-time signal extracted at $R=100\\,M_0$ no longer matches the four-QNM fit of the reconstructed final black hole, the claim that echoes persist in a consistent nonlinear theory fails.","tokens_in":26551,"feed_emoji":"🕳️","tokens_out":7318,"duration_ms":60297,"temperature":0.7,"pith_summary":"The paper tries to establish that gravitational-wave echoes—repeated pulses in a black hole's ringdown—are not an artifact of linearized toy models. In Einstein-Maxwell-scalar theory with a nonminimal electromagnetic-scalar coupling, hairy black holes can have a stable photon sphere, which turns the effective potential seen by perturbations into a multi-peaked cavity. The paper computes the quasinormal-mode spectrum of these backgrounds and shows by 1+1 nonlinear simulations that the echo pattern persists even when the perturbation is strong enough to grow the horizon area by a factor of about 32. This is claimed as the first example of echoes appearing in a consistent theory beyond a linearized analysis. If true, it means ringdown spectroscopy of extreme compact objects must take such cavity modes seriously in a fully nonlinear setting.","feed_headline":"Black hole echoes survive fully nonlinear collapse","feed_subtitle":"A 1+1 simulation shows cavity-trapped ringdown persists even when the horizon area grows by a factor of 32.","key_machinery":"The load-bearing object is the one-dimensional effective potential $V_\\phi$ in the linearized radial perturbation equation for the scalar field, written in tortoise coordinates; its wells and barriers turn the wave equation into a Schr\\\"odinger-like scattering problem whose cavity modes leak out through the barrier and arrive as echoes. For non-spherical perturbations, a coupled system of axial equations for the electromagnetic-led and gravitational-led channels plays the same role, with the coupling potential $V_{UH}$ controlling how echoes mix between channels. On the nonlinear side, the machinery is a 1+1 evolution in Painlev\\'e-Gullstrand-like coordinates with horizon excision, fourth-order finite differences, and Kreiss-Oliger dissipation, whose output is compared with a four-damped-sinusoid fit using frequencies computed in the linear frequency domain.","core_discovery":"The central claim is that in the EMS theory with coupling $F[\\phi]=e^{\\alpha\\phi^2}$, scalarized black holes near the critical charge develop an effective potential $V_\\phi$ with multiple maxima and a minimum, corresponding to a stable photon sphere. Such a potential traps low-frequency scalar perturbations in a cavity; they escape slowly by tunneling, producing repeated echo pulses in the time-domain response instead of a single prompt ringdown. The paper shows that a superposition of the first four quasinormal modes of the final black hole fits the late-time signal, and reports that fully nonlinear 1+1 evolutions, in which the horizon area grows by a factor of about 32, still display the echoes. The authors state that, to their knowledge, this is the first example of echoes in a consistent theory beyond a linearized analysis.","pith_inferences":["An implication the authors leave implicit is that echo morphology depends on how efficiently the horizon absorbs low-frequency radiation; one could test this by repeating the nonlinear run with a different coupling $F[\\phi]$ that produces a deeper or wider cavity and measuring how echo amplitude scales with horizon area growth.","The QNM-fit consistency between linear frequency-domain and nonlinear time-domain calculations suggests that the final state is well captured by the static scalarized solution; a sharper check would compare the full metric functions of the reconstructed final black hole with the simulated end state, not just the horizon radius and charge.","If echoes survive nonlinearity generically in EMS-like theories, post-merger gravitational-wave searches would need template banks containing cavity modes in addition to ordinary damped sinusoids, and the absence of such modes in observed ringdowns could constrain the coupling $\\alpha$."],"forward_implications":["Echoes are not wiped out by nonlinearities in this theory; a fully nonlinear spherical collapse onto a scalarized black hole retains the echo pattern seen in linear perturbation theory.","The late-time ringdown of a nonlinearly perturbed scalarized black hole is well described by a superposition of the first four quasinormal modes of the final static configuration, not the initial one.","Scalarized black holes with a stable photon sphere can have long-lived cavity modes with quality factor up to about 15, roughly four times the Schwarzschild value.","Axial $l=2$ perturbations also produce echoes, and the echo pattern depends on the relative amplitude of electromagnetic-led and gravitational-led initial data.","No nonlinear photon-sphere instability is found in the radial simulations, consistent with the theorem's hypotheses not applying to these black-hole backgrounds."],"supporting_citations":[{"why":"Supplies the scalarized black-hole solutions and the domain of existence used as backgrounds for the perturbation analysis.","marker":"[44]"},{"why":"Provides the master equation and effective potential for radial perturbations that the paper adopts and extends.","marker":"[45]"},{"why":"Introduced the echo phenomenon for horizonless compact objects that this paper extends to black holes.","marker":"[12]"},{"why":"Frames cavities in the effective potential as the origin of echoes and of modified ringdown signals.","marker":"[7]"},{"why":"Explains how low-frequency modes trapped in the cavity leak out by tunneling, producing repeated pulses.","marker":"[15]"},{"why":"Supplies the Painlevé-Gullstrand-like coordinate setup and numerical method used in the nonlinear 1+1 evolutions.","marker":"[41]"},{"why":"Motivates the theory's well-posedness and merger viability through prior binary simulations in the EMS model.","marker":"[50]"},{"why":"Identifies the nonlinear instability associated with stable photon spheres that the paper checks for in its radial simulations.","marker":"[17]"}],"fun_headline_variants":["First nonlinear black hole echoes in a consistent theory","Black hole echoes survive beyond linearized analysis","Echoes in black hole ringdown: nonlinearity isn't the end","Black hole cavity echoes persist in fully nonlinear evolutions","Nonlinear black hole echoes: a first in consistent gravity"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper's status as a consistent theory rests on the expectation, not a proof, that the Einstein-Maxwell-scalar equations remain well-posed for the chosen coupling and initial data; if they do not, the late-time echoes could be numerical artifacts rather than physical signals.","fun_headline_variants_meta":{"raw":{"variants":["First nonlinear black hole echoes in a consistent theory","Black hole echoes survive beyond linearized analysis","Echoes in black hole ringdown: nonlinearity isn't the end","Black hole cavity echoes persist in fully nonlinear evolutions","Nonlinear black hole echoes: a first in consistent gravity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000829,"raw_usage":{"total_tokens":3596,"prompt_tokens":894,"completion_tokens":2702,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":510,"completion_tokens_details":{"reasoning_tokens":2623}},"tokens_in":510,"tokens_out":2702,"duration_ms":20932,"temperature":1.0,"reasoning_tokens":2623,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T12:23:48.707978+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Repeat the $A=0.03\\,M_0$ collapse with an independent, manifestly strongly hyperbolic formulation of the same theory; if the late-time signal extracted at $R=100\\,M_0$ no longer matches the four-QNM fit of the reconstructed final black hole, the claim that echoes persist in a consistent nonlinear theory fails.","supporting_citations":[{"cited_title":"Nonperturbative gedanken experiments in Einstein-dilaton-Gauss-Bonnet gravity: nonlinear transitions and tests of the cosmic censorship beyond General Relativity","cited_arxiv_id":"2205.13007","evidence_quote":"Supplies the scalarized black-hole solutions and the domain of existence used as backgrounds for the perturbation analysis."}],"review_version":1}