{"id":"2b58bf49-38e5-4772-8e61-8429541da61f","arxiv_id":"2510.12198","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A driven holographic black hole exhibits turbulence whose horizon has fractal dimension D≈2.65 and whose dual fluid energy spectrum scales as k^{−1.79}.","lead":"This paper simulates a turbulent black hole in a holographic model of gravity, finding the horizon has a fractal shape with dimension about 2.65. The result hints that black holes and turbulent fluids may share a common, universal kind of disorder.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The quoted horizon fractal dimension is measured from the gauge-dependent function z_H(x,y); the authors concede no covariant definition exists, so D≈2.65 may not be a property of the black hole itself.","rationale":"The reader's weakest assumption flagged the quasi-steady-state issue as primary and the coordinate dependence as secondary. I agree the non-stationarity (v=4000–10000 with no large-scale friction) is a genuine fragility, but I find the coordinate/gauge dependence more load-bearing: it directly undercuts the interpretation of D as the fractal dimension of the horizon, independent of how long the run is. The authors themselves acknowledge this in Sec. VI, so this is an admitted missing piece rather than a speculative concern. The numerical scheme is described in detail and the spectrum results are plausible, and I do not see grounds for rejection; the correct disposition is a conditional acceptance requiring either a covariant estimator or a gauge-invariance test. Hence the reader's CONDITIONAL verdict stands, but the emphasis of the condition should be on the geometric definition of D, not only on the time-averaging window.","tokens_in":22136,"tokens_out":14249,"duration_ms":132342,"concrete_test":"Recompute the madogram dimension from the saved apparent-horizon data of the run behind Fig. 9 after transforming the radial coordinate to an affine-null coordinate (e.g., the gauge of Ref. [66]), keeping the same boundary coordinates and the same transect/fitting procedure. If the resulting D differs from 2.65 by more than the quoted ±0.02, the reported fractal dimension is gauge-dependent and the central claim is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that D≈2.65 is the fractal dimension of the turbulent black hole (Eq. 1, Sec. V). The measurement uses the madogram of the apparent-horizon location z_H(v,x,y), where z is the Bondi-Sachs areal-radius coordinate (Sec. II.B, App. B1). The function z_H is not a geometric invariant: residual coordinate freedom in the Bondi-Sachs gauge (radial redefinitions satisfying the areal condition, and spatial reparameterizations preserving g_vv=g_vi=0) changes the graph whose roughness the madogram quantifies. There is no demonstration that the inferred Hölder exponent, and hence D=2.65±0.02, survives such a change. Sec. VI explicitly states that 'since the fractal dimension of the turbulent black hole is a geometric quantity, its definition should be covariantly defined and be further investigated.' This is an admission that the present number is not established as a geometric property of the horizon. The agreement with Ref. [49] does not resolve the issue because that work also used a coordinate-based estimate. If D is gauge-dependent, the headline claim—a first nonlinear estimate of the horizon fractal dimension—is premature.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper constructs driven, fully nonlinear turbulent black holes in asymptotically AdS4 by solving the Einstein-scalar system in a Bondi–Sachs null foliation, with a random boundary scalar source serving as the external driving force. It reports an inverse energy cascade in the dual (2+1)-dimensional compressible fluid, a time-averaged total energy spectrum E(k) ~ k^{-1.79±0.03} in an inertial range k∈(10,65), compressible and incompressible components scaling as k^{-1.80±0.03} and k^{-1.99±0.03}, and a horizon fractal dimension D=2.65±0.02 obtained from the madogram of the apparent-horizon location z_H(v,x,y). The authors claim this is the first estimate of the fractal dimension from fully nonlinear driven black-hole dynamics and note agreement with the boundary-fluid result of Ref. [49].","tokens_in":22499,"tokens_out":4649,"duration_ms":43181,"significance":"If correct, the result is significant: it provides a first fully nonlinear, driven-evolution estimate of a turbulent horizon's fractal character, going beyond earlier derivative-expansion and force-free constructions, and it connects a bulk geometric measurement to a boundary spectral exponent. The paper also contains a useful technical contribution: a Bondi–Sachs scheme in which the evolution equations for the spatial metric are decoupled through an SO(2) rotation, allowing efficient GPU-based evolution. The numerical setup is described in unusual detail in Appendices A and B, and the reported exponents are given with time-averaging uncertainties. However, the central quantitative claims currently rest on a short inertial range and on a coordinate-dependent definition of the horizon's fractal dimension, both of which the manuscript itself partially acknowledges.","major_comments":[{"comment":"The inertial range used for the power-law fits is k∈(10,65), a factor of only 6.5 (about 0.8 decades), with the driving scale at kf≈100 and the grid resolution quoted only as '330 or 512 Fourier modes' in each spatial direction. No convergence or resolution study is reported (no variation of Nx, Ny, Nz, or δv; no comparison of results between 330 and 512 modes; no estimate of discretization error). The fitted exponent −1.79±0.03 and the fractal dimension D=2.65±0.02 are therefore not robustly separated from systematic effects. The authors should provide a resolution test (at least two well-separated resolutions) and a discussion of whether the short range and the proximity of kf to the dissipation scale render the exponents fitting-range dependent.","section":"Sec. IV.B, Figs. 6–7; App. B1"},{"comment":"The headline claim D≈2.65 is a 'fractal dimension of the turbulent black hole,' but the measurement is made on the graph z_H(v,x,y), the location of the apparent horizon in a particular Bondi–Sachs areal-radius gauge. The manuscript itself states in Sec. VI that 'since the fractal dimension of the turbulent black hole is a geometric quantity, its definition should be covariantly defined and be further investigated.' This is an admission that the present number has not been shown to be a geometric invariant. Because residual gauge freedom in the Bondi–Sachs construction (radial redefinitions consistent with the areal condition, and spatial diffeomorphisms preserving the gauge conditions) changes the function z_H, the madogram and hence D can change. The agreement with Ref. [49] does not remove this concern, since that work also used a coordinate-based estimator. At minimum, the authors sh","section":"Sec. V and Sec. VI"},{"comment":"The time average over v=4000–10000 is said to describe a quasi-steady turbulent state, but the system has no large-scale friction and the authors note that energy condensation will eventually occur. The mean kinetic energy in Fig. 5 appears to fluctuate but the paper does not demonstrate that the chosen window is stationary, e.g., by showing that the fitted exponent is stable under shifting the window, splitting the interval in half, or checking that the low-k energy content has saturated. Since the central exponents are obtained by averaging over this window, the possibility remains that the reported values are time-dependent transients of the inverse cascade rather than steady-state scaling exponents. A quantitative stationarity test should be added.","section":"Sec. IV.B, Figs. 5 and 7"},{"comment":"The apparent horizon is located by solving the nonlinear elliptic equation (B15) with a Newton–Krylov method, but no numerical tests are reported for this solver: no convergence of the Newton iteration, no comparison of the horizon location with an independent method, and no sensitivity of the fractal dimension to the horizon-finding tolerance or to the interpolation procedure. Since D=2.65±0.02 is estimated from z_H, the accuracy of z_H directly affects the central claim. Please report at least a basic validation of the horizon solver (e.g., convergence of H under mesh refinement and residual reduction).","section":"App. B2"}],"minor_comments":[{"comment":"The definition of Π_C appears to contain a typo: the term −(1/2) f ∂_z B should likely be −(1/2) f ∂_z C, consistent with the structure of (16) and with S_ΠC in App. B3. Please check and correct.","section":"Eq. (17)"},{"comment":"The text says 'Serval profiles' (should be 'several'), and the figure would benefit from a legend identifying the fitted ranges and exponents. The k^{-5} range is mentioned in the text but not clearly discussed; a sentence on its origin and lifetime would help.","section":"Fig. 2 and Sec. III.B"},{"comment":"The phrase '330 or 512 Fourier modes on each spatial xi direction' is ambiguous. Which resolution was used for the decaying runs and which for the driven runs? Please specify the actual grid for each simulation.","section":"App. B1"},{"comment":"For the madogram estimate, please state the exact range of separations r used in the linear fit, the number of transects and time snapshots used in the average, and whether the x- and y-direction transects are averaged with equal weight. The current Fig. 10 refers only to a shaded range that is said to correspond to Fig. 8.","section":"Sec. V"},{"comment":"Reference [50] is cited as a Springer book chapter; please give the full bibliographic details. Also correct 'GMRES' misspelled as 'GRMES' in App. B2.","section":"Refs."}],"recommendation":"major_revision","confidential_remarks":"The paper presents a potentially important numerical breakthrough, and the technical scheme is described in commendable detail. My main reservation is that the two headline numbers—the spectral exponent and the fractal dimension—are not yet supported with the numerical evidence a reader needs: no resolution study, a sub-decade inertial range, and a gauge-dependent horizon observable whose geometric status the authors concede. These are fixable within the manuscript's scope by adding numerical tests and substantially softening or qualifying the geometric claims. I therefore recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nShort version: this is a serious numerical relativity paper that does something new—drives a holographic black hole into a turbulent state with a scalar source, evolves the full Einstein-scalar system in Bondi–Sachs gauge, and reports a horizon fractal dimension D≈2.65. The scheme is described in unusual detail, and the decoupling of the h_ij evolution equations via an SO(2) rotation is a genuinely useful technical contribution.\n\nBut the headline claim is softer than it looks. D is extracted from the madogram of the apparent-horizon location z_H(v,x,y), and z_H is a coordinate-dependent function. The authors themselves write in Sec. VI that a covariant definition “should be further investigated.” That is an admission that the number is not established as a geometric property of the black hole. The agreement with Ref. [49] is suggestive but not resolving: that work used the same coordinate-based estimator on a fluid/gravity approximate metric. I’m not saying the result is wrong; I’m saying the headline is premature until someone checks gauge invariance or produces a covariant estimator.\n\nSecond soft spot: the averaging window v=4000–10000 is called quasi-steady, but the system has no infrared friction and the authors expect eventual energy condensation. The fitted exponent −1.79±0.03 is also taken over a short inertial range—about a decade and a half (k=10–65, forcing at k≈100). That is enough for a first look, not for a universal claim.\n\nThird, no code, no data, and no convergence or resolution study. For a numerical paper of this type, that is a real gap, though the detailed appendix partly compensates and the reported cost (five days on a GPU) suggests reproducibility is feasible.\n\nOn the plus side, the spectrum results—compressible component dominating, E_c~k^{−1.80}, E_i~k^{−1.99}—are consistent with earlier compressible fluid simulations and are physically interesting. The paper is honest about its own limitations, which is more than many are.\n\nBottom line: worth a serious referee, but the referee should ask for (1) a gauge-invariance check or a covariant definition of the fractal dimension, (2) a convergence study, and (3) evidence that the averaging window is not a pre-condensation transient. If those get fixed, the paper becomes a solid contribution; if not, the spectrum analysis still stands on its own.","headline":"A genuinely new numerical experiment with a useful Bondi–Sachs scheme, but the headline fractal dimension is measured in a gauge-dependent way the authors concede, so treat the number as provisional.","tokens_in":22947,"tokens_out":2979,"would_cite":true,"duration_ms":28691,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper claims that a holographically driven black hole horizon becomes a turbulent surface with fractal dimension D≈2.65, the first such estimate from full nonlinear evolution.","keywords":["holographic turbulence","AdS/CFT","Bondi-Sachs formalism","fractal dimension","inverse energy cascade","compressible turbulence","black hole horizon","numerical relativity"],"falsifier":"Repeat the same driving protocol with explicit large-scale drag or with a much longer run and re-fit the spectrum in successive windows: if the exponent drifts outside ±0.03 or D moves beyond ±0.02, the quoted values are transient. Alternatively, recompute D from a covariant intrinsic horizon quantity instead of the coordinate surface z=z_H: a significant change would mark the number as a gauge artifact.","tokens_in":22011,"feed_emoji":"🌀","tokens_out":6232,"duration_ms":54050,"temperature":0.7,"pith_summary":"The paper seeks to establish that a black hole horizon driven by a scalar source in a holographic spacetime becomes a genuinely turbulent surface, with a fractal dimension D≈2.65±0.02 and a boundary energy spectrum E(k)∼k^{−1.79±0.03}. This is the first estimate of the horizon fractal dimension obtained from fully nonlinear, driven black hole dynamics, and its agreement with earlier boundary-fluid simulations points toward a universal rough-geometry signature of turbulence. The deviation from the classical −5/3 scaling is attributed to the compressible character of the flow produced by scalar driving, and the cascade exponents match those seen in weakly coupled compressible two-dimensional fluid simulations. A sympathetic reader would take the paper as evidence that holography gives a concrete working definition of turbulent horizon roughness and a practical numerical route to measure it.","feed_headline":"Black hole horizon is a 2.65-dimensional fractal","feed_subtitle":"Driven AdS black hole simulations match boundary-fluid turbulence, hinting that horizon roughness is universal.","key_machinery":"The Bondi-Sachs null-foliation of the four-dimensional asymptotically anti-de Sitter spacetime: on each ingoing-null hypersurface the field equations reduce to nested radial integrations, and the two evolution equations for the shear modes B and C are decoupled by a field redefinition and an SO(2) rotation whose angle is K=∫ dz ∂_z B sinh C. This makes the radial operators time-independent, so spectral evolution can be accelerated on GPUs. The fractal dimension is then read off the apparent horizon with the madogram γ_1(r)=½⟨|z_H(x+r)−z_H(x)|⟩, using the scaling γ_1∝r^{2−D} and D=D_transect+1.","core_discovery":"The authors start from a static planar anti-de Sitter black hole and turn on a massive scalar field at the boundary whose value is updated by random, periodic white noise. Solving the full nonlinear gravitational-scalar system in a null-hypersurface gauge, they drive the boundary fluid into a quasi-steady two-dimensional turbulent state with an inverse energy cascade. Time-averaged fits over v=4000–10000 give E(k)∼k^{−1.79±0.03} in k∈(10,65); a Helmholtz decomposition yields compressible and incompressible components scaling as k^{−1.80±0.03} and k^{−1.99±0.03}, with the compressible part dominating. On the gravity side, the apparent horizon location z_H(v,x,y) is rough; applying the madogra","pith_inferences":["If D≈2.65 is an invariant property of the turbulent horizon, then horizon fluctuations carry a roughness signature that could in principle be probed through horizon-membrane observables in the dual plasma.","A decisive test is to switch to a divergence-free driving protocol: the paper's incompressible exponent near −2 predicts the fractal dimension should drop toward ≈2.58, distinguishing forcing-dependent roughness from a universal number.","Because z_H is coordinate-dependent, recomputing the dimension from a covariant intrinsic horizon quantity would show whether D is a geometric property or a gauge artifact.","The quoted exponents are averaged over a window before energy condensation; rerunning with large-scale friction or a longer averaging window would test whether −1.79 and 2.65 are true steady-state values."],"forward_implications":["The horizon's fractal dimension can be extracted from a direct, fully nonlinear bulk evolution, not only from boundary-fluid perturbation theory.","Scalar (compressible) driving yields a total cascade steeper than the classical −5/3 scaling, consistent with weakly coupled compressible fluid simulations.","The measured D≈2.65 is compatible with the earlier boundary-fluid result, reinforcing the idea of a universal fractal dimension for turbulent horizons.","The decoupled Bondi-Sachs scheme widens the accessible inertial range for holographic turbulence at moderate computational cost.","At higher resolution below the driving scale, the incompressible branch near k^{−2} should become visible and should give a second, slightly smaller fractal dimension."],"fun_headline_variants":["Black hole horizon roughness measured: D=2.65","Turbulent horizon fractal dimension hits 2.65","Horizon turbulence matches compressible fluid cascade","AdS black hole turbulence: k^-1.79 spectrum","Fractal dimension 2.65 for turbulent horizon"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The results presume that v∈[4000,10000] is a statistically stationary turbulent window, but with no large-scale friction the system is still transferring energy to the largest scales, so the fitted exponents and D may be time-dependent; the authors also note that the horizon location z_H is coordinate-dependent.","fun_headline_variants_meta":{"raw":{"variants":["Black hole horizon roughness measured: D=2.65","Turbulent horizon fractal dimension hits 2.65","Horizon turbulence matches compressible fluid cascade","AdS black hole turbulence: k^-1.79 spectrum","Fractal dimension 2.65 for turbulent horizon"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000165,"raw_usage":{"total_tokens":1088,"prompt_tokens":748,"completion_tokens":340,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":492,"completion_tokens_details":{"reasoning_tokens":275}},"tokens_in":492,"tokens_out":340,"duration_ms":2956,"temperature":1.0,"reasoning_tokens":275,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T09:59:22.231575+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Repeat the same driving protocol with explicit large-scale drag or with a much longer run and re-fit the spectrum in successive windows: if the exponent drifts outside ±0.03 or D moves beyond ±0.02, the quoted values are transient. Alternatively, recompute D from a covariant intrinsic horizon quantity instead of the coordinate surface z=z_H: a significant change would mark the number as a gauge artifact.","supporting_citations":[],"review_version":1}