{"id":"9e5154ec-ebb1-4a29-b5ae-292bb2c6970b","arxiv_id":"2509.06523","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Finite-size quantum critical systems thermalize in three regimes: recurrent non-equilibrium steady states, long-lived confined shock waves, or oscillatory decay, with near-complete energy swapping between the two sides.","lead":"Two holographic simulations follow two finite-size quantum systems at different temperatures that are joined at a perfect interface and then evolve as one isolated system. The simulations find three distinct thermalization patterns depending on system size and energy imbalance, including repeated formation and collapse of steady energy flows and a near-complete swap of the two sides' energies.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Brane no-matter condition (SI1 Eq. S35) is not enforced for the tanh initial profile; a residual boundary stress could change the observed recurrent dynamics.","rationale":"Agree with the reader: the weakest assumption is the consistency of the tensionless EOW brane condition for the chosen initial data. This is the only premise that is both load-bearing and explicitly acknowledged as unproved in the manuscript. The three-regime classification and the near-complete energy-swapping claim are statements about the isolated BCFT dual to the tensionless brane; if the brane actually carries a spurious stress tensor, the numerics realize a different boundary condition and the holographic interpretation of the results is not justified. The mirror-reflection class (S37) provides an exact consistency proof, but Eq. (4) only approaches it as α→0; the authors do not quantify the violation or check it against numerical error. A quantitative check of the S35 residual is the single experiment that would settle the concern. Other issues (schematic phase-diagram boundaries, missing code/data, Fig. 2 caption typo) are real but secondary: they affect confidence in the breadth or reproducibility of the claim, not its internal soundness. The recommended verdict remains conditional: the test could upgrade to accept if the residual is at truncation level, or downgrade to unverified if it is not. I therefore leave the reader's verdict unchanged.","tokens_in":14966,"tokens_out":8743,"duration_ms":93378,"concrete_test":"From the stored numerical data for the recurrent-NESS run (Δ=0.02, L̄=100), evaluate the left-hand side of Eq. (S35) — or equivalently the brane stress components \\hat T_uu, \\hat T_tu, \\hat T_yy from Eqs. (S23)-(S25) — on the grid at the brane location, sampling several times over the full simulation (e.g., t/L=0, 20, 100, 200). Normalize by the local energy density and compare with the truncation errors reported in SI5 (Figs. S7-S8). If the normalized residual stays at or below the truncation error and does not grow over the recurrence cycles, the concern is resolved; if it is larger or trends upward, the no-matter condition is not enforced and the central three-regime claim is not established for the stated theory.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The manuscript's own SI1 (Eq. S35 and following text) states that the no-matter brane condition \\hat T_uu=0 is 'not guaranteed' for generic initial data, and that a proof for the actual initial data would require a full nonlinear analysis. The only class shown to satisfy it is the mirror-reflection class (Eq. S37), requiring all odd x-derivatives of s1 and B to vanish on the brane. The initial profile (4), a tanh-smoothed step with α=0.07, is not in this class: at x=±L, derivatives of s1 are proportional to sech^2(1/α) and are exponentially small but not zero. The evolution algorithm (SI2) solves the bulk equations and imposes the derived boundary conditions (S38), but the no-matter condition S35 is an extra constraint that is not separately enforced or monitored. If its residual is not at the level of the truncation error, the simulations describe a brane carrying a small spurious stress tensor, i.e., a different BCFT with boundary degrees of freedom, and the clean mapping to the isolated tensionless-brane system used to interpret the three regimes is not established. Since the recurrent-NESS claim relies on repeated boundary reflections over t≈200L, even a small violation could in principle accumulate and change the revival amplitude. The paper does not report any check of the S35 residual.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper uses holographic duality (AdS/BCFT) to simulate the real-time thermalization of a finite-size strongly coupled CFT prepared with a spatially varying energy density. The setup is an asymptotically AdS_4 spacetime with two dynamical, tensionless end-of-the-world branes, dual to an isolated strip CFT. The initial state is a tanh-smoothed step in energy density (Eq. 4). By numerically solving the bulk Einstein equations, the authors identify three thermalization regimes controlled by two dimensionless parameters: the energy imbalance Δ and the scaled system size L̄. For large L̄ and small Δ, a recurrent NESS forms and collapses, with near-mirror revivals at t=2L/c_s. For large Δ, a long-lived confined shock wave persists, while for small L̄, the system undergoes oscillatory decay without persistent structures. Late-time decay is compared with linear quasi-normal modes of the final black brane, including nonlinear forced-oscillation corrections. The supplementary information contains the metric ansatz, boundary-condition derivation, numerical scheme, perturbation analysis, and error estimates.","tokens_in":15324,"tokens_out":7114,"duration_ms":87040,"significance":"If the results hold, they constitute a useful extension of holographic studies of NESS from infinite systems to finite, isolated systems, with falsifiable predictions: a three-regime phase diagram in (Δ, L̄), a revival time 2L/c_s, and late-time decay governed by black-brane QNMs. The numerical machinery is a notable strength: the paper reports two independent error estimators at the 10^-6–10^-4 level (Fig. S7–S8), energy conservation to 10^-9–10^-6 (Fig. S3), and entropy monotonicity (Fig. S2). The late-time comparison with linear QNMs appears to be an external, parameter-free benchmark, which is valuable. The main risk to the central claim is the consistency of the initial data with the tensionless-brane no-matter condition, as detailed below.","major_comments":[{"comment":"The no-matter brane condition \\hat T_uu=0 is not guaranteed for the initial data actually used. The paper itself states that for generic initial conditions (S35) requires a full nonlinear proof, and the only proven class is the mirror-reflection class (S37). The tanh-smoothed step (4) with α=0.07 is not exactly in this class: solving s1 from the energy profile gives odd x-derivatives at x=±L that are exponentially small but nonzero (∼sech^2(1/α)). The evolution scheme imposes (S38) but does not monitor the residual of (S35). Since the recurrent-NESS claim relies on many boundary reflections over t≈200L, a small spurious boundary stress could accumulate and change revival amplitudes. Please report the residual of (S35) on Q during the runs, or repeat the key simulations with initial data constructed to lie exactly in the class (S37), and show that the three regimes persist. This is load-b","section":"SI1, Eqs. (S35)–(S41); main text Sec. II, Eq. (4)"},{"comment":"The phase diagram is obtained for a single interface sharpness α=0.07. The statement that this is 'sufficiently small' to not qualitatively alter the dynamics is not supported by a convergence study in α. Because α also controls the magnitude of the boundary odd derivatives of s1, this is intertwined with the consistency issue above. Please provide at least one representative trajectory (e.g., the recurrent-NESS case) for a smaller and a larger α (say 0.03 and 0.1), showing that the revival amplitude and the qualitative phase boundaries are stable.","section":"Sec. II, Eq. (4) and Fig. 2"}],"minor_comments":[{"comment":"The text says the energy current for oscillatory decay is shown in the 'bottom-right panel of Fig. 2', but the figure caption identifies the bottom-right panel as recurrent NESS and the top-left panel as oscillatory decay. Please correct this cross-reference.","section":"Sec. IV (Oscillatory decay)"},{"comment":"Typo: 'preform' should be 'perform'.","section":"SI4"},{"comment":"Typo: 'impost' should be 'impose'.","section":"SI1"},{"comment":"Eq. (5) introduces constants τ and ω without specifying their origin. If they are taken from the QNM spectrum of the final black brane, state this explicitly in the main text; if they are fits, clarify what is being fitted. This would make the predictive content of the late-time comparison precise.","section":"Eq. (5)"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Finite-size effects in holographic thermalization are the real new thing here. The three regimes — recurrent NESS, confined shock, oscillatory decay — and the near-complete energy swap do not reduce to the known infinite-size NESS results, and the late-time match to QNMs of the final black brane is a genuine external benchmark with no fitted constants. The numerics look careful: two independent error estimators, energy conservation, entropy monotonicity. Credit where due: this is a serious computation.\n\nThe main soft spot is the brane boundary condition. The supplement itself states that the no-matter condition on the tensionless brane (S35) is not guaranteed for general initial data; only the mirror-reflection class (S37) is proven to satisfy it. The tanh-smoothed profile only approximates that class, with odd x-derivatives of s1 exponentially small but not zero. The code imposes the derived conditions (S38) but does not separately enforce or monitor S35. Over t ~ 200L a small residual boundary stress could accumulate and change the revival amplitude. The paper reports no check of S35. This is not a fatal objection by itself — it could well be that the residual is at truncation level — but it is a specific, answerable one. The fix is easy: compute the S35 residual from the numerical data or argue that the tanh profile is effectively in the basin. Until that's shown, the clean mapping to an isolated tensionless-brane BCFT is not established.\n\nThe phase diagram rests on few parameter points and the boundaries are schematic; that's acceptable in a letter, but it makes \"three regimes\" a set of examples rather than a mapped phase diagram. The Fig. 2 caption inconsistency (which run is recurrent NESS) is minor. The absence of released code or data is a negative for this kind of numerical claim.\n\nWho this is for: people working on holographic thermalization, AdS/BCFT, or finite-size quantum critical transport. It deserves a serious referee; the referee should ask for the S35 residual check and a commitment to release data. I would not cite it in my own work until that's resolved.","headline":"New finite-size holographic thermalization phenomena, but the no-matter brane constraint for the actual initial data is unproven and needs a direct numerical check.","tokens_in":15803,"tokens_out":2332,"would_cite":false,"duration_ms":26602,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81T40","83C57","83-08"],"pacs":[],"model":"deepseek-v4-flash","headline":"A holographic simulation shows finite-size quantum critical systems thermalize through three distinct regimes controlled by energy imbalance and scaled size, with wave reflection causing near-complete energy swapping.","keywords":["thermalization","holographic duality","quantum critical systems","non-equilibrium steady states","shock waves","quasi-normal modes","end-of-the-world brane","finite-size systems"],"falsifier":"Run the same thermalization with an initial profile that exactly satisfies the odd-derivative mirror-reflection conditions—for example, one constructed from even-parity spatial modes—against the tanh-smoothed step; if the recurrence peak at t=2L/c_s or the extracted quasi-normal frequencies shift measurably, the approximate boundary compatibility is contributing to the claimed dynamics.","tokens_in":14826,"feed_emoji":"🌊","tokens_out":5989,"duration_ms":65666,"temperature":0.7,"pith_summary":"The paper addresses how an isolated, finite-size strongly coupled conformal system approaches equilibrium when two halves at different temperatures are suddenly joined. Using holographic duality, it simulates the full gravitational dynamics in a box and finds that thermalization is not diffusive: it is a wave phenomenon controlled by the energy imbalance Δ and the rescaled size L̄. Depending on those two numbers, the system either cycles through repeated formation and collapse of a non-equilibrium steady state, sustains a long-lived confined shock wave with boundary reflections, or bypasses structure entirely and decays in synchronized oscillations. The paper claims this wave-dominated transport plus boundary reflection can nearly swap the energies of the two halves, and that the late-time decay is quantitatively the linear quasi-normal ringdown of the final equilibrium black hole with nonlinear corrections.","feed_headline":"Finite-size critical systems thermalize in three distinct ways","feed_subtitle":"Holographic simulations show energy swapping via reflected waves, with behavior set by size and energy imbalance.","key_machinery":"The central object is a holographic black brane in an asymptotically anti-de Sitter spacetime with a tensionless end-of-the-world brane at each spatial end; the end-of-the-world brane—a bulk hypersurface anchored to the field-theory boundary—encodes the isolated-system condition of zero energy current at the boundary. Real-time evolution uses a characteristic numerical scheme in Eddington-Finkelstein coordinates, with the apparent horizon area supplying entropy production. Late-time decay is analyzed by Fourier-decomposing the energy current and comparing the modes to quasi-normal modes of the final equilibrium black brane, including forced-oscillation corrections when a dominant mode drives","core_discovery":"The paper claims that when two finite-size conformal field theories at different energy densities are joined at a perfectly transmitting interface, thermalization is controlled by two dimensionless parameters: the energy imbalance Δ=(E_R−E_L)/E_L and the scaled size L̄=L((E_R+E_L)/2)^{1/3}. In the holographic dual—a black brane confined between two tensionless end-of-the-world branes—the early contact region forms a non-equilibrium steady state with shock and rarefaction fronts moving at the sound speed. The finite boundaries reflect these fronts. For large L̄ and small Δ, the reflection produces near-mirror recurrences of the initial state at t=2L/c_s, with the NESS re-forming and collapsin","pith_inferences":["A natural but untested extension: in two spatial dimensions the same phase diagram should persist with curved shock fronts, since the paper states that generalization to other dimensions is straightforward but does not demonstrate it.","If the energy-swap mechanism survives weak external coupling, timed separation of the two subsystems after roughly one sound-crossing time could serve as a cooling or energy-harvesting protocol; the paper mentions cooling applications without detailing an engineered cycle.","The recurrence time t=2L/c_s is a concrete, parameter-free prediction: ultracold-atom or nanowire realizations of conformal systems could look for oscillatory energy exchange at that period, which would support the holographic picture and whose absence would challenge it.","The forced-oscillation crossover suggests that mode amplitudes, not just frequencies, retain memory of the initial energy gap; varying Δ in small systems could expose this memory effect in the late-time ringdown."],"forward_implications":["For large systems with small energy imbalance, thermalization is not monotonic: the system cycles through repeated formation and collapse of a non-equilibrium steady state, with near-mirror revivals every sound-crossing time 2L/c_s.","For large energy imbalance, shock waves persist through many boundary reflections and dominate late-time energy transport, while rarefaction waves quickly homogenize.","For sufficiently small systems, dissipation dominates: no sustained NESS or shock structure forms, and the system enters rapid synchronized oscillatory decay.","The late-time decay is described by linear quasi-normal modes of the final equilibrium black brane, with higher modes entering a forced-oscillation regime where the fundamental mode acts as the driver.","Wave-propagated energy transfer plus boundary reflection allows near-complete energy swapping between subsystems after separation; total entropy still increases, so the second law is not violated."],"supporting_citations":[{"why":"Supplies the conformal-field-theory construction of steady energy flow in non-equilibrium states, the baseline for the NESS concept.","marker":"[13]"},{"why":"Establishes NESS formation in infinite quantum critical systems, the infinite-size case that the finite-size analysis extends.","marker":"[14]"},{"why":"Identifies shock waves, rarefaction waves, and NESS in quantum critical systems, providing the wave-front mechanism used here.","marker":"[15]"},{"why":"Gives black-brane steady states, the holographic description of NESS the paper builds on.","marker":"[17]"},{"why":"Previous numerical demonstration of NESS formation in 3+1 dimensions, a direct comparison point for the present simulations.","marker":"[18]"},{"why":"Supplies the characteristic numerical evolution method used to solve the bulk gravitational equations in real time.","marker":"[19]"},{"why":"Develops aspects of the holographic boundary construction with end-of-the-world branes, the setup for finite-size systems.","marker":"[23]"},{"why":"Introduces the holographic dual of boundary conformal field theory, the source of the end-of-the-world brane description.","marker":"[25]"}],"fun_headline_variants":["Holographic view reveals three thermalization paths in finite-size critical systems","Size and energy gap dictate how quantum critical systems thermalize","Reflected waves enable energy swapping in finite-size critical systems","Thermalization in finite-size critical systems: three regimes from holography"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The reliable long-time evolution assumes the boundary stays a perfect, matter-free wall, but that is proven only for specially symmetric starting states; the smoothed-step initial profile used here obeys that condition only approximately.","fun_headline_variants_meta":{"raw":{"variants":["Holographic view reveals three thermalization paths in finite-size critical systems","Size and energy gap dictate how quantum critical systems thermalize","Reflected waves enable energy swapping in finite-size critical systems","Thermalization in finite-size critical systems: three regimes from holography"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000609,"raw_usage":{"total_tokens":2657,"prompt_tokens":713,"completion_tokens":1944,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":457,"completion_tokens_details":{"reasoning_tokens":1879}},"tokens_in":457,"tokens_out":1944,"duration_ms":14003,"temperature":1.0,"reasoning_tokens":1879,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T23:29:21.699613+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same thermalization with an initial profile that exactly satisfies the odd-derivative mirror-reflection conditions—for example, one constructed from even-parity spatial modes—against the tanh-smoothed step; if the recurrence peak at t=2L/c_s or the extracted quasi-normal frequencies shift measurably, the approximate boundary compatibility is contributing to the claimed dynamics.","supporting_citations":[{"cited_title":"Bernard and B","cited_arxiv_id":null,"evidence_quote":"Supplies the conformal-field-theory construction of steady energy flow in non-equilibrium states, the baseline for the NESS concept."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes NESS formation in infinite quantum critical systems, the infinite-size case that the finite-size analysis extends."},{"cited_title":"Lucas, K","cited_arxiv_id":null,"evidence_quote":"Identifies shock waves, rarefaction waves, and NESS in quantum critical systems, providing the wave-front mechanism used here."},{"cited_title":"Amado and A","cited_arxiv_id":null,"evidence_quote":"Gives black-brane steady states, the holographic description of NESS the paper builds on."},{"cited_title":"Non-equilibrium steady state formation in 3+1 dimensions","cited_arxiv_id":"2103.10435","evidence_quote":"Previous numerical demonstration of NESS formation in 3+1 dimensions, a direct comparison point for the present simulations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the characteristic numerical evolution method used to solve the bulk gravitational equations in real time."},{"cited_title":"Fujita, T","cited_arxiv_id":null,"evidence_quote":"Develops aspects of the holographic boundary construction with end-of-the-world branes, the setup for finite-size systems."},{"cited_title":"Takayanagi, Holographic Dual of BCFT, Phys","cited_arxiv_id":null,"evidence_quote":"Introduces the holographic dual of boundary conformal field theory, the source of the end-of-the-world brane description."}],"review_version":1}