{"id":"fc9517d6-f9a4-4bef-8ba7-e901e09e27b6","arxiv_id":"2412.14755","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"For a boosted strongly coupled plasma in d=2,3,4, holographic subregion complexity increases with temperature, velocity, and region size, and diverges as the Lorentz factor squared, gamma^2, when velocity approaches light speed.","lead":"This paper computes a quantum information measure, subregion complexity, for a strongly coupled plasma moving at constant velocity using holographic duality. It finds that complexity grows with temperature, velocity, and subregion size, and diverges as the velocity approaches the speed of light.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The γ² universality is asserted for arbitrary l but all scaling fits use l=1, so the central claim is unsupported.","rationale":"The central advertised result is not merely that HSC increases with T, v, and l, which the plotted curves plausibly support, but that the divergence as v→1 is universally γ² for arbitrary T and l in d=2,3,4. The presented fits only address l=1 and a small set of temperatures, and the authors explicitly note that high temperatures cannot be reached near v→1. The l-dependence plots are taken at fixed v and therefore contain no information about the divergence exponent. Without a scaling analysis at multiple l values, the claimed universality in l is unsupported. The reader's weakest assumption about Eq. (10) is also real: the PDE is un-derived, no code or data accompany the paper, and the only convergence figure is not a mesh-convergence study. That compounds the concern, but the missing l-scaling analysis alone means the headline universality cannot be accepted as stated. A conditional verdict is appropriate pending a scaling analysis at multiple l and a verified numerical pipeline, so the reader's conditional verdict remains unchanged.","tokens_in":12750,"tokens_out":9098,"duration_ms":81412,"concrete_test":"For d=2, T=1/(6π), compute C(v) at l=0.5, 1.0, and 2.0 with the same numerical pipeline, and fit each curve to a γ^p + b with p as a free parameter. If the best-fit p differs by more than the fit uncertainty across these l values, the claimed universality in subregion length fails; if p≈2 for all three, the l-universality claim survives, though a mesh-convergence study would still be needed to trust the pipeline.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's headline result is a dimension-independent γ² divergence for arbitrary temperature and subregion length. The numerical evidence for that scaling is restricted to l=1 in every fitting plot: d=2 is fitted at T=1/6π and T=1/60π (Fig. 2), and d=3,4 are each fitted at a single temperature T=1/4π (Fig. 5). No fit of C(v) as a function of l is reported anywhere. The l-dependence plots (Figs. 4, 7, 8) are all at fixed v, so they test monotonicity but not the v→1 divergence exponent. The text also states that high temperatures cannot be probed near v→1, so the 'arbitrary T' part of the claim is not demonstrated either. In addition, Eq. (10) is presented without derivation, no code or data are provided, and the single convergence illustration is not a mesh-convergence study. These facts combined mean the central universal claim is an extrapolation from a narrow parameter set, and the scaling-in-l assertion is not backed by any direct numerical test.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript computes holographic subregion complexity (HSC) for a strip-like boundary region in a boosted black-brane background in d = 2, 3, 4. The authors numerically solve a nonlinear PDE for the bulk profile t(x,z), using the HRT surface from [42] as boundary data, and then evaluate the enclosed volume. They report three qualitative results: HSC increases with plasma temperature T, boost velocity v, and strip length l; and as v approaches 1, the subtracted complexity C-C0 diverges as the Lorentz factor squared, gamma^2, which they tentatively call universal across dimensions, temperatures, and subregion lengths. The fits use the form a(T)gamma^2+b(T) and have small reported relative errors.","tokens_in":13005,"tokens_out":6856,"duration_ms":60616,"significance":"If the gamma^2 divergence and the monotonicity results are correct, the paper would provide the first all-temperature, all-velocity numerical treatment of holographic subregion complexity in this stationary moving-plasma background, going beyond the perturbative calculations of [39,40]. The numerical computation of a genuinely nonlinear PDE with HRT boundary data is nontrivial, and the small fit errors and the consistency with the perturbative gamma^2 result of [40] are genuine strengths. However, the paper's headline universality claim currently rests on a narrow set of fits and on an unproven identification of the fitted constant b with numerical error, so the significance of the paper depends on how much of that claim can be supported by additional evidence or by a more restricted statement.","major_comments":[{"comment":"The central claim that the gamma^2 divergence holds for arbitrary subregion length l is not supported by the data: every velocity-scaling fit is performed at l=1. Figures 4, 7, and 8 vary l only at fixed v, which tests monotonicity but cannot determine the v->1 divergence exponent. To support the abstract and the conclusions, the authors should either extract the divergence exponent from data at several values of l (for example by fitting C as a function of gamma at each l and showing the coefficient scales consistently), or explicitly restrict the universality claim to l=1.","section":"Numerical results, d=2 and d=3,4; Figs. 2 and 5"},{"comment":"The claim of universality for arbitrary temperature is similarly not demonstrated. In d=2 the fits use only T=1/(6 pi) and T=1/(60 pi); in d=3 and d=4 each scaling fit is at a single temperature T=1/(4 pi). Moreover, Appendix A concedes that velocities close to one cannot be probed at high temperatures, so the v->1 exponent is not checked across a meaningful temperature range. The authors should either restrict the claim to the temperatures actually probed or perform fits at additional temperatures and show the gamma^2 coefficient and exponent are stable.","section":"Numerical results; Figs. 2, 5, 6 and Appendix A"},{"comment":"Equation (10), the second-order nonlinear PDE that is the core of the paper's numerical method, is stated without derivation. Since the entire numerical computation and all subsequent claims depend on this equation, the authors should provide its derivation from the volume functional (8)-(9), or at least a clear reference where it is derived. They should also comment on the existence, uniqueness, and single-valuedness of the solution t(x,z) for the HRT boundary data, since the text asserts single-valuedness without discussion.","section":"Holographic subregion complexity, Eq. (10)"},{"comment":"The convergence evidence is not a mesh-convergence study. Figure 10 shows the number of unrelaxed grid nodes as a function of iteration count for one example at one set of parameters; it does not show how the computed complexity C depends on the mesh sizes Delta x and Delta z, nor does it quantify discretization error. This matters because the authors attribute the fitted constant b to numerical precision. I ask for a systematic grid-convergence test for the subtracted complexity, for example repeating a representative case at cell sizes 10^-3, 10^-4, and 10^-5, and reporting the resulting values of C and of the fitted b. Availability of the code and data would also substantially strengthen reproducibility.","section":"Appendix A and Fig. 10"},{"comment":"The treatment of the constant b in the fit C_bar = a(T)gamma^2 + b(T) is not justified. The authors state that a and b are of the same order and that both are much smaller than C0, and they then conclude that b is due to numerical error. The comparison with C0 is irrelevant for this conclusion: b is comparable to a (for example, for T=1/(6 pi), a=0.9260 and b=-0.8898), so over the accessible velocity range the data are fit by a constant plus a gamma^2 term. The asymptotic v->1 exponent would still be gamma^2 if b is constant, but the stronger statement C proportional to a(T)gamma^2 requires b to vanish in the continuum limit. The authors should test this by increasing the numerical precision and by a mesh-convergence study, rather than inferring it from a and b being small compared with C0.","section":"Numerical results, d=2; Fig. 2 and text after it"}],"minor_comments":[{"comment":"There are several typos and reference formatting issues: 'Hubney' should be 'Hubeny' in the first occurrence and in reference [3]; Appendix A contains 'differntial'; reference [14] contains 'Infromation'; several arXiv identifiers are missing their archive prefix, e.g. [35] should read arXiv:hep-ph/0607062.","section":"Introduction"},{"comment":"The notation is potentially confusing because C is used both for the unsubtracted complexity in Eq. (4) and for the subtracted quantity in Eq. (11a). The authors should introduce distinct symbols for the two quantities and use them consistently in the figures and captions.","section":"Equation (11) and surrounding text"},{"comment":"The quantity Cv is used for dC/dv but is not defined before the figure; please define it explicitly and, if possible, state the numerical differentiation method used to obtain the slope.","section":"Fig. 2, right panel"},{"comment":"The red and blue curves in Fig. 1 are described as the plus and minus branches but the caption does not say which color corresponds to which branch; in addition, the discretization details in Appendix A state the maximum cell size is O(10^-3) but do not specify how nx, nz, and the tolerance were chosen for each run.","section":"Fig. 1 and Appendix A"},{"comment":"The conclusions repeat the universality claim for 'all temperatures and subregion lengths', but the preceding numerical section explicitly notes that high temperatures cannot be probed near v->1. The wording should be reconciled with the actual data coverage so that the claim does not overstate the numerical evidence.","section":"Conclusions"}],"recommendation":"major_revision","confidential_remarks":"The paper makes an interesting numerical contribution, and I do not see a definitive error in the monotonicity results. However, the headline universality claim is currently an extrapolation from fits at l=1 and at very few temperatures, and the treatment of the fitted constant b as a numerical artifact is not established. I would be willing to accept a revision that either substantially expands the numerical evidence for the gamma^2 divergence across l and T, or carefully rewrites the central claim to match the actually probed parameter range. The lack of a derivation for Eq. (10) and the absence of a mesh-convergence study are additional issues that need to be addressed in the revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my read. The paper does something legitimately new: it computes holographic subregion complexity for a boosted black brane non-perturbatively across d=2,3,4, for a range of temperatures and velocities, and shows monotonic growth in T, v, and strip length. The fits are impressively tight (sub-percent relative errors), and the consistency with the perturbative gamma^2 result from [40] is a real point in its favor. The numerical work is not trivial, and the authors are honest about the numerical challenges near v=1 and about citing the earlier perturbative calculations.\n\nThe soft spots are real but not fatal. The headline claim—universal gamma^2 divergence for arbitrary T and l—is not backed by the reported fits. Every scaling fit uses l=1, and only one or two temperatures per dimension. The l-dependence plots are all at fixed v, so they test monotonicity, not the divergence exponent. The paper's own footnote admits that high temperatures cannot be probed near v=1, so the 'arbitrary T' part is also not demonstrated. Calling this universal is an extrapolation. I think the stress-test note is correct on this point.\n\nTwo technical gaps also bother me. Equation (10) is stated without derivation, and for a numerical paper that is a missing load-bearing piece: the reader cannot check that the PDE is the correct Euler-Lagrange equation. And there is no mesh-convergence study, no code or data release, just one example solution and an iteration-count plot. That is not enough to establish numerical reliability, especially when the central claim is a fitted scaling law with 'b' discarded as numerical error.\n\nThat said, the monotonic trends are very likely correct, and the gamma^2 scaling is externally supported by [40]. This is a solid numerical study overreaching in its conclusions, not a flawed one. The paper deserves a serious referee, but the referee should push for a derivation or at least a clear justification of (10), a mesh-convergence analysis, and fits at multiple l and T values before the universality claim is accepted. I would bring it to a reading group as an example of how numerical hysteresis can look convincing while missing the stated domain of validity.","headline":"A genuinely new numerical map of subregion complexity in boosted plasmas, but the universal gamma^2 claim is fitted from l=1 data only and is broader than the evidence.","tokens_in":13505,"tokens_out":1493,"would_cite":true,"duration_ms":13284,"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":"In a moving strongly coupled plasma, holographic subregion complexity grows with temperature, velocity, and strip length, and diverges as the square of the Lorentz factor as velocity approaches light speed, in d=2,3,4.","keywords":["holographic subregion complexity","boosted black brane","moving strongly coupled plasma","Lorentz factor divergence","HRT surface","finite difference method","strip subregion","AdS/CFT"],"falsifier":"Compute the same subtracted complexity at velocities closer to 1 (for example $v=0.999$) with a finer mesh and a convergence study, and extract the divergence exponent from a log-log plot of $\\mathcal{C}$ versus $\\gamma$ at fixed $T$ and $l$. If the exponent is not 2, or if the fitted function requires higher powers of $\\gamma$, the claimed universal $\\gamma^2$ divergence fails. An analytic derivation of the leading velocity dependence from the volume functional that yields a different Lorentz-factor power would also settle the question.","tokens_in":12568,"feed_emoji":"⚡","tokens_out":14698,"duration_ms":101524,"temperature":0.7,"pith_summary":"The paper studies holographic subregion complexity in a strongly coupled plasma moving at constant velocity, using the holographic proposal that identifies the complexity of a mixed state with the volume of the codimension-one hypersurface enclosed by the HRT surface (the covariant extremal surface used in holographic entanglement entropy). For a strip-shaped boundary region aligned with the boost, the authors solve the resulting extremal-volume equation numerically by finite differences in boundary dimensions $d=2,3,4$, without restricting temperature or velocity. Their central finding is that holographic subregion complexity increases with temperature, velocity, and subregion length, and that as the velocity approaches the speed of light the complexity diverges as the square of the Lorentz factor, $\\gamma^2$, for arbitrary temperature and subregion length. They conclude that this $\\gamma^2$ divergence is a universal feature of strip subregions in a moving plasma across the three dimensions studied.","feed_headline":"Holographic complexity diverges like γ² in hot moving plasma","feed_subtitle":"Numerics in d=2,3,4 give the same Lorentz-factor-squared divergence for strips at any T and l.","key_machinery":"The central object is the volume functional $V=\\int dz\\,dx\\,\\mathcal{V}$ for the region enclosed by the HRT surface in the boosted black brane background, with integrand (9) and the associated Euler-Lagrange equation (10). The load-bearing computational step is the finite-difference solution of this nonlinear PDE for $t(x,z)$, with boundary values taken from the HRT profile $(x_\\pm(z), t_\\pm(z))$; the iteration relaxes all grid nodes to a tolerance of order $10^{-6}$. The Lorentz factor enters through the metric function $g(z)=\\gamma^2(z/z_h)^d$, so the $\\gamma^2$ divergence of the complexity as $v\\to1$ tracks the divergence of this metric coefficient. The subtracted quantities $\\mathcal{C}$ and $\\hat{\\mathcal{C}}$ compare the complexity against the static-plasma and pure-AdS baselines, and the fitted form $a(T)\\gamma^2+b(T)$ is what converts the numerical curves into the claimed universal scaling law.","core_discovery":"The paper establishes numerically that for a strip-like subregion of length $l$ and infinite width in a $d$-dimensional strongly coupled plasma moving with constant velocity $v$ parallel to the strip's short edge, the holographic subregion complexity $C_A = V_{\\gamma_A}/(8\\pi G_N)$ is an increasing function of temperature $T$, velocity $v$, and subregion length $l$ for $d=2,3,4$. Working with the boosted black brane metric and using the HRT surface (the covariant extremal surface used in holographic entanglement entropy) as boundary data, the authors solve the nonlinear partial differential equation (10) for the extremal volume slice $t(x,z)$ by finite differences. Subtracting either the static-plasma value or the pure-AdS value gives the quantities $\\mathcal{C}$ and $\\hat{\\mathcal{C}}$ defined in (11), both of which are positive and increase with all three parameters. Fitting the velocity dependence shows that as $v\\to 1$ the complexity diverges as $a(T)\\gamma^2+b(T)$, with $a,b\\ll C_0$, so the divergence is characterized by the Lorentz factor squared $\\gamma^2$; the same behavior is found in $d=2,3,4$ for arbitrary temperature and subregion length. The paper reads this as evidence for a universal $\\gamma^2$ divergence for strip subregions when the plasma moves parallel to the short edge.","pith_inferences":["The $\\gamma^2$ scaling may reflect a kinematic boost effect rather than a dynamical property of the plasma; one way to test this is to repeat the calculation with the strip at an oblique angle to the boost direction and check whether the Lorentz-factor power changes.","The same finite-difference pipeline could be applied to other mixed-state observables, such as the entanglement wedge cross-section or mutual information, to see whether they share the $\\gamma^2$ tail; if they do, the divergence would be a general property of the boosted geometry rather than of the specific volume functional.","The fitted form $a(T)\\gamma^2+b(T)$ with $b\\ll C_0$ suggests that in the strict large-$\\gamma$ regime the subtracted complexity is controlled by a single temperature-dependent coefficient $a(T)$; deriving $a(T)$ analytically would turn the numerical fit into a closed-form prediction."],"forward_implications":["Holographic subregion complexity of a strip region in a moving plasma is larger than in a static plasma at the same temperature, so specifying the mixed state requires more information once the plasma moves.","As the plasma velocity approaches the speed of light, the subtracted complexity diverges like $\\gamma^2$ for fixed $T$ and $l$ in $d=2,3,4$, making the divergence a candidate universal signature of the boost.","For small subregion lengths ($T_l\\gg T$) the complexity approaches the pure-AdS value, so short-distance probes are insensitive to both temperature and velocity.","Higher boundary dimension increases the absolute value of the complexity and makes the differences between velocities and temperatures more pronounced at larger subregion lengths.","At low temperature the velocity dependence is suppressed because the metric coefficients reduce to $f(z)\\to1$ and $g(z)\\to0$, recovering the static vacuum limit."],"supporting_citations":[{"why":"Proposes identifying holographic subregion complexity with the volume of the codimension-one hypersurface enclosed by the HRT surface; this is the definition the paper computes.","marker":"[19]"},{"why":"Defines the covariant HRT surface that bounds the volume whose size is interpreted as subregion complexity.","marker":"[3]"},{"why":"Supplies the numerical HRT surface profiles for the boosted black brane background, which serve as boundary conditions for the volume PDE (10).","marker":"[42]"},{"why":"Earlier perturbative computation of holographic subregion complexity in a moving plasma reported $\\Delta C\\propto\\gamma^2$ with a vanishing constant term for $d=2$; the paper's all-parameter numerical result extends and matches that scaling.","marker":"[40]"},{"why":"Precedent for solving a time-dependent extremal-volume PDE with boundary data from the HRT surface, motivating the paper's numerical approach.","marker":"[45]"},{"why":"Supplies the finite-difference method used to solve the nonlinear PDE for the extremal volume slice.","marker":"[44]"}],"fun_headline_variants":["Moving plasma sends holographic complexity up as γ²","Complexity of moving plasma diverges universally like γ²","Velocity drives holographic subregion complexity to γ² blow-up","Holographic complexity in hot moving plasma: γ² divergence","Universal γ² divergence in moving plasma complexity"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that equation (10) is the correct equation for the surface whose volume defines the complexity, and that it has a unique smooth solution matching the HRT surface; the paper states this without derivation, checks only one sample solution, and gives no mesh-convergence study.","fun_headline_variants_meta":{"raw":{"variants":["Moving plasma sends holographic complexity up as γ²","Complexity of moving plasma diverges universally like γ²","Velocity drives holographic subregion complexity to γ² blow-up","Holographic complexity in hot moving plasma: γ² divergence","Universal γ² divergence in moving plasma complexity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000413,"raw_usage":{"total_tokens":2159,"prompt_tokens":995,"completion_tokens":1164,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":611,"completion_tokens_details":{"reasoning_tokens":1085}},"tokens_in":611,"tokens_out":1164,"duration_ms":7969,"temperature":1.0,"reasoning_tokens":1085,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T11:55:46.172121+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the same subtracted complexity at velocities closer to 1 (for example $v=0.999$) with a finer mesh and a convergence study, and extract the divergence exponent from a log-log plot of $\\mathcal{C}$ versus $\\gamma$ at fixed $T$ and $l$. If the exponent is not 2, or if the fitted function requires higher powers of $\\gamma$, the claimed universal $\\gamma^2$ divergence fails. An analytic derivation of the leading velocity dependence from the volume functional that yields a different Lorentz-factor power would also settle the question.","supporting_citations":[{"cited_title":"Numerical Techniques in Electromagnetics,","cited_arxiv_id":null,"evidence_quote":"Supplies the finite-difference method used to solve the nonlinear PDE for the extremal volume slice."}],"review_version":1}