{"id":"a4f129a9-08c5-4b7b-8d2f-cfe2de721b82","arxiv_id":"2508.03519","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"For RN-AdS black holes in quintessence, Lyapunov exponents of unstable circular null and timelike geodesics decrease with the quintessence parameter, reach a finite cutoff, and at phase transitions jump with a claimed |T-Tc|^(1/2) scaling.","lead":"Charged anti-de Sitter black holes surrounded by quintessence dark energy are shown to have orbital-stability exponents that shrink as the quintessence strength rises and disappear beyond a finite cutoff. The paper connects jumps in these exponents at phase transitions to van der Waals critical behavior, proposing them as a dynamical thermometer for black hole thermodynamics.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The derivation of the central δ=1/2 scaling is invalid as written: Eq. (28) assumes a nonzero (∂²T/∂r²_+)_c while Eq. (11) defines the critical point by that derivative vanishing. The advertised critical exponent is therefore unsupported by the paper's own mathematics.","rationale":"The paper has two advertised outputs: (i) quintessence suppresses Lyapunov exponents and introduces a finite cutoff for unstable null circular orbits, and (ii) the discontinuity Δλ at the first-order phase transition scales as |T_p − T_c|^{1/2}, matching van der Waals behavior. Output (i) follows from the standard effective-potential analysis in Section III and is credible: for ω = −1/2 the quintessence term modifies f/r² in a way that can remove the potential maximum for sufficiently large b or r_+, and the contours in Figs. 5 and 10 are consistent with that picture. The central quantitative claim (ii) is where the paper fails. Section IV's derivation uses Eq. (28) with a nonzero (∂²T/∂r²_+)_c. However, Eq. (11) defines the critical point precisely by the vanishing of that derivative. Consequently Eq. (29) divides by zero and Eqs. (30)–(31) do not follow. A correct mean-field treatment requires expanding T(r_+, b) around (r_{+,c}, b_c) including the linear-in-δb shift, the cross term ∂²T/∂r_+∂b, and the cubic term ∂³T/∂r_+³, then imposing Maxwell/equal-free-energy coexistence; such a derivation is absent. The numerical fit in Fig. 14 is not a substitute: no error bars are given, and the functional form Δλ̃ = k sqrt(T̃−1) is fitted to data whose b-values are not stated. The malformed Eq. (14) for Q_c is a further correctness defect but is secondary: the numerical critical values in Eq. (16) appear internally consistent with Eqs. (12)–(13). Because the 1/2 exponent is the paper's headline result, the current text does not support it. The concrete Maxwell-construction test would tell whether the claim can be repaired; based on standard mean-field behavior and the cited literature on RN-AdS, I expect the exponent would survive, but that does not make the present derivation acceptable. No code or data were provided, and several typographical errors prevent full verification of the figures.","tokens_in":14433,"tokens_out":23825,"duration_ms":257811,"concrete_test":"Perform an independent numerical Maxwell construction near criticality for Q = 0.1: choose b = b_c (1 − 10^{−n}) for n = 3,...,7, solve T(r_s) = T(r_l) = T_p together with equal free energy F(r_s) = F(r_l), compute Δλ = λ_l − λ_s using the same Lyapunov formulas (Eqs. 19 and 23), and fit log Δλ versus log|T_p − T_c|. If the fitted slope is 1/2 to within, say, 1%, the underlying δ = 1/2 claim survives but the derivation in Section IV must still be rewritten with a valid joint expansion; if the slope deviates significantly, the advertised central claim is false.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's headline quantitative result is the critical exponent δ = 1/2 for the Lyapunov-exponent discontinuity Δλ = λ_l − λ_s (Eqs. 30–31). That result depends entirely on Eq. (28), which expands T(r_+) near the critical point as T(r_+) = T_c + (1/2)(∂²T/∂r²_+)_c (r_+ − r_{+,c})² + O(r_+ − r_{+,c})³. This expansion is only valid if (∂²T/∂r²_+)_c ≠ 0. But the paper's own criticality condition, Eq. (11), is ∂T/∂r_+ = 0 and ∂²T/∂r²_+ = 0 at the critical point. Hence the quadratic coefficient in Eq. (28) vanishes, the denominator in Eq. (29) is zero, and Eqs. (30)–(31) do not follow. A correct derivation would need a joint expansion in r_+ − r_{+,c} and b − b_c (or T − T_c), retaining the cubic term in T(r_+) and the cross term ∂²T/∂r_+∂b, and imposing the equal-free-energy Maxwell construction for coexistence. The numerical fit in Fig. 14 and Eq. (32) cannot substitute for this missing derivation, especially without error bars or stated b-values for the fitted points. The qualitative statements about quintessence suppressing λ and producing a finite cutoff are plausible and supported by the standard geodesic analysis in Section III, but the central critical-exponent claim is unsupported as written. A secondary defect is the malformed Eq. (14) for Q_c, which is imaginary for the plotted case ω = −1/2; this further undermines confidence in the formal derivations, though the numerical critical values in Eq. (16) appear internally consistent with Eqs. (12)–(13).","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies four-dimensional Reissner-Nordström-AdS black holes surrounded by a quintessence field and analyzes their thermodynamic phase structure and geodesic Lyapunov exponents. It derives thermodynamic quantities, locates critical points, computes Lyapunov exponents for null and timelike circular orbits, and observes that increasing the quintessence parameter suppresses the exponents and introduces a finite cutoff absent in RN-AdS. It then defines the discontinuity Δλ at the first-order transition as an order parameter and claims it scales as |T_p − T_c|^{1/2}, i.e., critical exponent 1/2.","tokens_in":14859,"tokens_out":4357,"duration_ms":49648,"significance":"If the central scaling claim were fully supported, the paper would provide an interesting dynamical probe of black-hole phase transitions and a dark-energy-dependent cutoff. The qualitative sections are based on standard geodesic equations and contain useful numerical surveys. However, the derivation of the headline critical exponent rests on an inconsistent Taylor expansion, and the analytic formula for the critical charge is not valid for the plotted case. The finite-cutoff observation is the most credible novel part, but it is qualitative. Overall the paper does not yet substantiate its main quantitative claim.","major_comments":[{"comment":"The expansion T(r_+) = T_c + (1/2)(∂²T/∂r_+²)_c (r_+ − r_{+,c})² + O(r_+ − r_{+,c})³ assumes (∂²T/∂r_+²)_c ≠ 0. This is contradicted by the paper's own criticality condition in Eq. (11), which states ∂T/∂r_+ = 0 and ∂²T/∂r_+² = 0 at the critical point. The quadratic coefficient in Eq. (28) therefore vanishes, the denominator in Eq. (29) is zero, and Eqs. (30)–(31) do not follow. A valid derivation would require a joint expansion in r_+ − r_{+,c} and b − b_c (or T − T_c) with the Maxwell construction; the manuscript provides no such argument.","section":"Section IV, Eq. (28)"},{"comment":"For the case ω = −1/2 used throughout the numerical analysis, the factor √(−3 + 9ω) is imaginary, yet Fig. 1 displays real critical values. Either the formula contains a typographical or algebraic error or the plotted curve is not generated from Eq. (14). Since the critical parameters underpin the phase diagram and the subsequent critical-exponent analysis, this error needs to be fixed.","section":"Section II, Eq. (14)"},{"comment":"The numerical fit Δλ̃ = k√(T̃−1) is presented without the underlying data points, error bars, or the values of b and Q used; the caption only says 'points near the critical point.' This fit cannot serve as independent evidence for δ = 1/2, especially when the analytical derivation is invalid.","section":"Section IV, Fig. 14 and Eq. (32)"}],"minor_comments":[{"comment":"The title and abstract contain typos: 'Ads' should be 'AdS' and 'These work suggest' should be 'These works suggest'.","section":"Title and Abstract"},{"comment":"The sentence 'where where L = L/l is scaled...' contains a duplicated 'where' and should be corrected.","section":"Section III.A"},{"comment":"The statement in Fig. 5b that 'for larger b exceeding 1.2, λ vanishes entirely' is made without a quantitative criterion for the disappearance of unstable circular orbits; the threshold should be defined precisely.","section":"Section III.A"},{"comment":"References [48] and [51] are identical, and reference [47] lacks publication details; check and consolidate the bibliography.","section":"References"},{"comment":"The conclusion repeats the claim that Δλ yields a critical exponent of 1/2, but this claim depends on the invalid derivation in Section IV and should be revisited.","section":"Section V"}],"recommendation":"reject","confidential_remarks":"The central quantitative claim is unsupported as written because Eq. (28) contradicts Eq. (11). The qualitative finite-cutoff result is interesting but does not by itself justify publication. I also note that the critical-exponent result overlaps with Refs. [42,58,59] from overlapping author groups, and the new element of quintessence does not appear to change the derivation's structure. The editor may wish to check the novelty and citation practices."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this morning. The genuinely new thing is the finite cutoff: in the RN-qAdS spacetime, null-geodesic Lyapunov exponents from unstable circular photon orbits vanish identically for sufficiently large r+ or b, unlike pure RN-AdS where they persist asymptotically. That follows directly from the effective potential losing its maximum, and the λ-T plots with multivalued SBH/IBH/LBH branches mirroring the free-energy swallowtail are a clean, plausible demonstration. The qualitative physics—quintessence suppresses geodesic instability and can shut off chaotic photon orbits entirely—looks right and worth having.\n\nBut the centerpiece, the claim that Δλ = λ_l − λ_s is an order parameter with critical exponent 1/2, is not supported by the paper's own math. Eq. (28) expands T(r+) near criticality keeping the quadratic term with coefficient (∂²T/∂r²_+)_c, while the paper's own criticality condition, Eq. (11), is precisely ∂²T/∂r²_+ = 0 at the critical point. So the quadratic coefficient vanishes, Eq. (29) divides by zero, and the square-root scaling does not follow. A correct derivation would require a joint expansion in r+ − r_c and b − b_c (or Q − Q_c), together with the Maxwell equal-area construction for coexistence. The numerical fit in Fig. 14, with no error bars and no stated b-values for the fitted points, cannot substitute for that derivation.\n\nSmaller issues add to the impression of haste: Eq. (14) gives an imaginary Q_c for ω = −1/2, the exact case plotted in Fig. 1; Eq. (32) writes the fit as k√(T̃ − 1) while the axis label says Δλ̃ = 2.19... T̃ − 1 (likely a typo); and no code or data are provided. None of these affect the qualitative result, but they reinforce that the manuscript needs careful revision.\n\nWho is this for? Anyone working on Lyapunov exponents as thermodynamic probes or on quintessence black hole phase structure. The finite-cutoff observation is a genuine contribution, and the branch structure is pedagogically useful. The critical-exponent claim is the kind of thing a referee could catch; I think the paper deserves a serious referee, not a desk reject, because the qualitative content is solid and the quantitative problem is fixable in principle—though as written the central advertised result does not stand.\n\nMy recommendation: send to peer review, require major revision, and ask the authors to either supply the correct joint expansion or remove the δ = 1/2 claim entirely.","headline":"The finite-cutoff observation is genuinely new, but the advertised δ=1/2 scaling is unsupported because Eq. (28) contradicts the paper's own criticality condition Eq. (11).","tokens_in":15427,"tokens_out":3268,"would_cite":false,"duration_ms":34937,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Lyapunov exponents of unstable circular orbits in quintessence-RN-AdS black holes jump at the small/large phase transition, with the jump scaling as the square root of the temperature gap.","keywords":["Lyapunov exponent","black hole thermodynamics","phase transition","quintessence","RN-AdS black hole","critical exponent","geodesic chaos","photon sphere"],"falsifier":"Numerically compute the coexisting small- and large-black-hole branches from the equal-free-energy condition for $b$ close to $b_c$, evaluate $\\Delta\\lambda = \\lambda_l - \\lambda_s$ on those branches, and fit $\\log \\Delta\\lambda$ against $\\log |T_p - T_c|$; if the fitted slope is not $1/2$, the claimed critical exponent is not supported.","tokens_in":14194,"feed_emoji":"🕳️","tokens_out":13616,"duration_ms":126972,"temperature":0.7,"pith_summary":"The paper attempts to establish that a purely dynamical quantity--the Lyapunov exponent measuring how quickly nearby geodesics diverge around a black hole--can act as an order parameter for thermodynamic phase transitions. For charged AdS black holes surrounded by a quintessence field, it argues that the Lyapunov exponent of unstable circular photon and massive-particle orbits becomes multivalued below the critical point, jumps discontinuously at the first-order transition temperature, and has a jump $\\Delta\\lambda$ that vanishes at the critical point as $|T_p - T_c|^{1/2}$, matching the van der Waals exponent. It further argues that quintessence suppresses geodesic chaos and, unlike pure RN-AdS, imposes a finite cutoff at which unstable circular orbits cease to exist and $\\lambda$ drops to zero. If correct, this connects geodesic instability to black hole thermodynamics and to the presence of dark energy.","feed_headline":"Black-hole chaos probe sees phase transition with exponent 1/2","feed_subtitle":"Lyapunov jumps at the small/large transition with a 1/2 exponent; quintessence adds a cutoff absent in RN-AdS.","key_machinery":"The central object is the Lyapunov exponent $\\lambda$ of unstable circular geodesic orbits, defined from the second derivative of the effective potential $V_{\\mathrm{eff}}(r)$ at the unstable radius $r_c$: for null geodesics $\\lambda = \\sqrt{-r_c^2 f(r_c) V_{\\mathrm{eff}}''(r_c)/(2L^2)}$, with an analogous expression for timelike geodesics. The paper uses this exponent as a dynamical mirror of the black hole's thermodynamic phases, with the jump $\\Delta\\lambda = \\lambda_l - \\lambda_s$ between the small- and large-black-hole branches serving as the order parameter. The metric ingredient is the quintessence-deformed Reissner-Nordström-AdS metric function $f(\\tilde r) = 1 - 2\\tilde M/\\tilde r + \\tilde Q^2/\\tilde r^2 + \\tilde r^2/l^2 - \\tilde b/\\tilde r^{3\\omega+1}$, whose normalization factor $\\tilde b$ suppresses $\\lambda$ and eventually removes the unstable circular orbit altogether.","core_discovery":"The paper's central claim is that for the Reissner-Nordström-AdS black hole enveloped by quintessence (RN-qAdS), the Lyapunov exponent $\\lambda$ of unstable circular orbits--both null (photon) and timelike (massive particle) geodesics--reproduces the thermodynamic phase structure of the black hole. Below the critical quintessence parameter $b_c$ or charge $Q_c$, the $\\lambda$ versus Hawking temperature $T$ curves show three branches (small, intermediate, large black holes), and at the phase transition temperature $T_p$ where free energies of small and large black holes are equal, $\\lambda$ jumps discontinuously. At criticality the jump $\\Delta\\lambda = \\lambda_l - \\lambda_s$ vanishes, and the paper derives $\\Delta\\lambda \\sim |T_p - T_c|^{1/2}$, identifying the jump as an order parameter with critical exponent $1/2$, consistent with the van der Waals fluid. A second claim is that quintessence introduces a finite cutoff: for sufficiently large $b$ or $r_+$, the effective potential no longer has an unstable circular-orbit extremum, so $\\lambda$ vanishes identically, a behavior that does not occur in pure RN-AdS spacetimes.","pith_inferences":["Extension beyond the paper: if the cutoff is confirmed, black hole shadow or photon-ring instability measurements could in principle constrain the quintessence normalization factor $b$, because the temperature or horizon radius at which $\\lambda$ vanishes depends on $b$.","Extension beyond the paper: the same $\\Delta\\lambda$ order-parameter construction could be applied to other dark-energy or dark-matter modified black hole spacetimes, such as perfect fluid dark matter backgrounds, to test whether the critical exponent is universal or model-dependent.","Extension beyond the paper: a two-variable expansion in $(r_+ - r_{+,c})$ and $(b - b_c)$, rather than the paper's one-variable expansion, would clarify whether the $1/2$ exponent survives at the true critical point."],"forward_implications":["The multivalued $\\lambda$--$T$ curves mirror the free-energy swallowtail, so a purely dynamical measurement of geodesic instability can locate the first-order small/large black hole transition and the second-order critical point.","The exponent $1/2$ ties the Lyapunov jump to van der Waals universality, extending the same critical behavior already reported for RN-AdS, Gauss-Bonnet AdS, and Born-Infeld AdS black holes to the quintessence case.","For $b > b_c$ or $Q > Q_c$ the $\\lambda$--$T$ curve is monotonic and single-valued, giving a dynamical signature that the black hole is in a single stable phase with no transition.","The finite cutoff where $\\lambda$ vanishes provides a testable distinction between quintessence and pure RN-AdS black holes, since the latter retains nonzero photon-sphere Lyapunov exponents asymptotically."],"supporting_citations":[{"why":"Establishes the method of using Lyapunov-exponent discontinuities to probe small/large black hole phase transitions in RN-AdS, which this paper extends to quintessence.","marker":"[42]"},{"why":"Provides the timelike Lyapunov exponent formula and the prior critical-exponent-1/2 result for Gauss-Bonnet AdS black holes used as a consistency check.","marker":"[58]"},{"why":"Reports the same 1/2 exponent for Born-Infeld AdS black holes, one of the results this paper's claim is said to be consistent with.","marker":"[59]"},{"why":"Gives the metric for black holes surrounded by quintessence that underlies Eq. (4).","marker":"[18]"},{"why":"Supplies the action and thermodynamic framework for charged RN-AdS black holes surrounded by quintessence used in Section II.","marker":"[57]"},{"why":"Provides the barotropic perfect-fluid Lagrangian used to model the quintessence contribution in the action.","marker":"[22]"},{"why":"Establishes the upper bound on the quintessence parameter for Hawking-Page phase transitions, used to delimit the phase-transition region.","marker":"[56]"},{"why":"Introduces the van der Waals-like small/large black hole phase transition in charged AdS spacetimes that the quintessence case is compared against.","marker":"[7]"}],"fun_headline_variants":["Lyapunov jump orders black-hole phase transition with 1/2 exponent","Black-hole chaos probe: Lyapunov jump gives exponent 1/2","Quintessence sets Lyapunov cutoff; jump marks black-hole phase transition","Lyapunov jump as order parameter: black-hole criticality matches van der Waals"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The 1/2 exponent rests on assuming the quadratic term in the expansion of $T(r_+)$ dominates near the critical point, yet the paper's criticality conditions set the coefficient of that quadratic term to zero.","fun_headline_variants_meta":{"raw":{"variants":["Lyapunov jump orders black-hole phase transition with 1/2 exponent","Black-hole chaos probe: Lyapunov jump gives exponent 1/2","Quintessence sets Lyapunov cutoff; jump marks black-hole phase transition","Lyapunov jump as order parameter: black-hole criticality matches van der Waals"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001248,"raw_usage":{"total_tokens":5143,"prompt_tokens":992,"completion_tokens":4151,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":608,"completion_tokens_details":{"reasoning_tokens":4066}},"tokens_in":608,"tokens_out":4151,"duration_ms":38401,"temperature":1.0,"reasoning_tokens":4066,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T04:26:04.562888+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Numerically compute the coexisting small- and large-black-hole branches from the equal-free-energy condition for $b$ close to $b_c$, evaluate $\\Delta\\lambda = \\lambda_l - \\lambda_s$ on those branches, and fit $\\log \\Delta\\lambda$ against $\\log |T_p - T_c|$; if the fitted slope is not $1/2$, the claimed critical exponent is not supported.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the method of using Lyapunov-exponent discontinuities to probe small/large black hole phase transitions in RN-AdS, which this paper extends to quintessence."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the timelike Lyapunov exponent formula and the prior critical-exponent-1/2 result for Gauss-Bonnet AdS black holes used as a consistency check."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports the same 1/2 exponent for Born-Infeld AdS black holes, one of the results this paper's claim is said to be consistent with."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the metric for black holes surrounded by quintessence that underlies Eq. (4)."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the action and thermodynamic framework for charged RN-AdS black holes surrounded by quintessence used in Section II."},{"cited_title":"Minazzoli and T","cited_arxiv_id":null,"evidence_quote":"Provides the barotropic perfect-fluid Lagrangian used to model the quintessence contribution in the action."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the upper bound on the quintessence parameter for Hawking-Page phase transitions, used to delimit the phase-transition region."},{"cited_title":"Chamblin, R","cited_arxiv_id":null,"evidence_quote":"Introduces the van der Waals-like small/large black hole phase transition in charged AdS spacetimes that the quintessence case is compared against."}],"review_version":1}