{"id":"8653f745-30df-4935-8212-96afb7ecaadd","arxiv_id":"2412.12353","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"At equal R-charge chemical potentials, the holographic AdS5-Reissner-Nordström black brane has a negative R-charge diffusion coefficient below a critical μ/T, making the low-temperature phase unstable.","lead":"This paper shows that the AdS-Reissner-Nordström black brane, the holographic description of N=4 super-Yang-Mills theory at equal R-charge chemical potentials, is dynamically unstable at low temperature. The instability appears as negative R-charge diffusion, so the low-temperature phase must be described by a different gravity solution.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The dynamical-instability claim rests on quasinormal-mode equations (7)-(8) that are asserted without derivation; if those equations omit couplings or the boundary conditions are mis-specified, the reported sign flip in the diffusion coefficient could be an artifact.","rationale":"The thermodynamic part of the paper is internally consistent and well supported: the equation of state (2) follows from known STU thermodynamics, and the Hessian eigenvalues in §3 demonstrably change sign at κ=1. The hydrodynamic argument in Appendix B is also coherent: for a symmetric charge matrix, pair-difference fluctuations obey diffusion with D=(σ11−σ12)/(χ11−χ12), and a negative susceptibility eigenvalue with positive conductivity implies D<0. What is not independently established is the dynamical content. The quasinormal-mode equations (7)–(8) are stated without derivation, the decoupling of the difference sector is assumed, and the numerical solution lacks convergence data. The phrase 'one can show' is the weakest link because the paper's headline claim about the low-temperature phase goes beyond thermodynamics: it asserts that AdS-RN is dynamically unstable. That claim would fail if the correct fluctuation system contained additional couplings or if the boundary conditions implemented a different holographic prescription. This is a missing-support concern rather than a demonstrated contradiction, so a conditional verdict is appropriate. I am not raising disagreement with the consensus view or questioning the authors' integrity; the concern is internal to the paper's presentation. A concrete independent check of eqs. (7)–(8) and a second numerical implementation would resolve it.","tokens_in":11223,"tokens_out":4570,"duration_ms":45823,"concrete_test":"Independently linearize the action (9) around the equal-κ background, fix δA^a_u=0, impose the scalar-field constraint, and project the fluctuation equations onto the subspace with ΣδA=0 and Σδs=0; verify that the resulting decoupled system matches eqs. (7)–(8) identically and that no residual coupling to δg_{μν}, ΣδA, or the remaining scalar combination contributes at linear order. Then solve the same boundary-value problem with an independent numerical method, such as Chebyshev pseudospectral collocation with a Newton root-finder, at κ=1.1 and q=0.05, and check that Im(ω) becomes positive and that D(κ)=−lim_{q→0} Im(ω)/q^2 follows the reported linear (1−κ) behavior near κ=1. Optionally, compute σ11−σ12 from the Kubo formula in the same background to confirm that it is positive and finite as κ→1.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim that AdS5-RN is not the low-temperature phase depends on the numerical result that the R-charge diffusion coefficient D becomes negative for κ>1. That result comes directly from solving eqs. (7)–(8), which are introduced with 'one can show' and no derivation from the action (9). The decoupling of the difference combinations E^a_z − (1/3)ΣE^b_z and s^a − (1/3)Σs^b from the total charge mode, the metric fluctuations, and the remaining scalar combination is asserted rather than demonstrated. The boundary treatment is also stated, not derived: for E_z the exponent-1 mode is chosen at u=0, and for s the logarithmic branch coefficient A is set to zero. If the correct holographic prescription selects a different linear combination, or if couplings to the total-charge or sound-channel sectors survive, the quasinormal-mode crossing at κ=1 could be an artifact of an inconsistent truncation. The thermodynamic instability itself is solid—the Hessian eigenvalues in §3 change sign at κ=1—and the hydrodynamic relation D=(σ11−σ12)/(χ11−χ12) is coherent. However, the headline claim additionally requires the dynamical mode to exist as computed, and the sign of σ11−σ12 is imported from ref [14] rather than computed here. These are missing derivations and checks, not demonstrated contradictions, so the paper should be accepted only after the QNM system is independently verified.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyzes the thermodynamic and hydrodynamic stability of strongly coupled N=4 supersymmetric Yang-Mills theory at finite R-charge chemical potentials, using the five-dimensional STU black brane as the holographic dual. Starting from the known equation of state, the authors derive explicit thermodynamic stability conditions and show that the Hessian matrix of the energy density develops negative eigenvalues when the STU parameter κ exceeds 1 along the equal-chemical-potential line. They then argue, using a general hydrodynamic analysis for multiple conserved charges, that this thermodynamic instability should be accompanied by a dynamical instability in the charge-diffusion sector. The central holographic computation is a numerical quasinormal-mode analysis of the equal-charge background: eqs. (7)-(8) are solved for the decoupled fluctuations, and a diffusive mode is found to cross into the upper complex frequency half-plane for κ>1, corresponding to a negative R-charge diffusion coefficient. The paper concludes that the low-temperature phase of N=4 SYM with equal chemical potentials is not described by the AdS-Reissner-Nordström black brane, because this background is unstable to fluctuations that are invisible in a minimal Einstein-Maxwell truncation but present in the STU model. Appendix B develops the hydrodynamic derivation of the diffusion mode and the relation between D and the conductivity and susceptibility matrices.","tokens_in":11571,"tokens_out":4476,"duration_ms":43477,"significance":"If the central numerical claim is correct, this is a significant result: it provides a concrete example in which a consistent supergravity truncation displays a dynamical instability that is invisible in the minimal charged black brane, and it demonstrates in a holographic setting the general relation between thermodynamic Hessian eigenvalues and hydrodynamic instabilities. The thermodynamic part is clean and explicit: the stability conditions (4), the Hessian eigenvalues, and the susceptibility eigenvalues are derived from a closed-form equation of state, and the hydrodynamic derivation in Appendix B is coherent. The paper also gives a falsifiable prediction—the sign flip of the diffusion coefficient at κ=1—and it is honest about the external input from ref. [14] for the equation of state. However, the headline dynamical-instability claim rests on a numerical quasinormal-mode calculation whose equations are asserted rather than derived, whose boundary conditions are stated without full justification, and whose numerical implementation is not documented. The result is therefore plausible but not yet independently verifiable from the manuscript as written.","major_comments":[{"comment":"The quasinormal-mode equations are introduced with the phrase 'one can show' and no derivation is given. The central claim that the diffusion coefficient becomes negative for κ>1 depends entirely on these equations being the correct linearized equations for the decoupled fluctuations E^a_z and s^a. In particular, the manuscript does not demonstrate that these combinations decouple from the total-charge mode, from the metric fluctuations, and from the remaining scalar combination. Please provide a derivation in an appendix or a precise reference to a source where the decoupling and the explicit equations are obtained. Without this, the numerical results in Fig. 2 are not independently checkable.","section":"Section 4, Eqs. (7)-(8)"},{"comment":"The treatment of the boundary conditions at u=0 is load-bearing. For E_z the exponent-1 mode is chosen, and for s the logarithmic branch coefficient A is set to zero in the expansion s = A u log u + B u + ... . The manuscript invokes the 'standard holographic recipe', but for coincident indicial exponents the choice of boundary condition is not uniquely fixed by the usual normalizability criterion. A different linear combination of the logarithmic and power-law branches would change the quasinormal spectrum and could remove or alter the reported sign flip. The authors should justify this boundary condition, for example by deriving it from the poles of the dual retarded correlator or by an independent variational principle.","section":"Section 4, boundary conditions at u=0"},{"comment":"No numerical method is described, and no convergence tests, resolution checks, or error estimates are reported. The right panel of Fig. 2 shows D(κ) as a smooth curve near κ=1, but no data points are displayed and the fit underlying the linear dependence is not specified. Since the paper's main conclusion is the crossing of the diffusive mode into the upper half-plane at κ≈1, the numerical procedure needs to be reproducible: the discretization scheme, the way the diffusive mode is identified among the quasinormal spectrum, and the extrapolation to q→0 should all be stated explicitly.","section":"Section 4, Fig. 2 and numerical procedure"},{"comment":"The hydrodynamic prediction D = (σ11−σ12)/(χ11−χ12) is used together with the statement that σ11−σ12 stays finite and non-zero as κ→1, so that the sign of D is determined by the thermodynamic instability. The sign of σ11−σ12 is imported from ref. [14] and is not computed in this paper. Please state explicitly that this sign is taken as an external input, and either verify it by an independent calculation or clearly mark it as a reliance on ref. [14]. As written, the agreement between the numerical D(κ) and the hydrodynamic formula is presented as stronger evidence than the paper's own computations support.","section":"Section 4 and Appendix B, Eq. (19)"}],"minor_comments":[{"comment":"The name 'Behrnd' should be 'Behrndt' to match the reference list and the original paper.","section":"Abstract and Introduction"},{"comment":"The sentence 'The eigenvalues diverge at κ=1, suggesting that the corresponding diffusion coefficients vanish at κ=1' is potentially misleading: the divergence of χ11−χ12 implies D→0 only if σ11−σ12 remains finite, and the divergence of the susceptibility itself is a thermodynamic signal. Please rephrase to make the logic explicit.","section":"Section 3, text after Eq. (5)"},{"comment":"The right panel lacks axis labels and data points; please show the actual numerical values, the quoted linear fit, and, if possible, an estimate of the numerical error.","section":"Fig. 2"},{"comment":"The relation 'G5 = πL3/2Nc2' is ambiguous; it should be written as G5 = π L^3/(2 N_c^2) to avoid confusion.","section":"Appendix A, text before Eq. (9)"},{"comment":"Ref. [39] is cited as 'To appear, (2025)' with no further information; if the companion paper is available, please provide an arXiv identifier or a more complete citation.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the thermodynamic analysis is solid. The main issue is that the quasinormal-mode calculation, which carries the dynamical-instability conclusion, is not sufficiently documented in the manuscript: the equations are asserted, the boundary conditions are stated without derivation, and the numerical method is absent. These are fixable within the manuscript's scope by adding an appendix with the derivation and by reporting the numerical details, so I do not think rejection is warranted. I would also suggest that the authors explicitly delineate which inputs come from ref. [14] and which are computed here, since one of the authors of the present paper is also an author of ref. [14]. This is not a concern about correctness, but it would help the reader assess the independence of the hydrodynamic check."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Interesting paper, but the headline dynamical instability is asserted rather than demonstrated. The thermodynamic part is clean and solid; the QNM part needs work before I'd trust it.\n\nWhat's new: previous work established thermodynamic instability for generic STU black branes and hydrodynamic instability for a single chemical potential. This paper shows that even at the symmetric point—equal chemical potentials—the AdS5-RN brane has unstable R-charge diffusion modes. That's a genuinely new result, and if right, it says the simplest charged brane is not the low-temperature description of N=4 SYM with equal chemical potentials. The Hessian analysis in §3 is the best part: the eigenvalue crossing at κ=1 is explicit, the eigenvectors correctly identify the difference channels, and the stability region in Fig 1 is clear. The hydrodynamic argument in Appendix B is also coherent: for equal charges, the differences decouple and D=(σ11−σ12)/(χ11−χ12). With positive σ, a negative susceptibility eigenvalue gives D<0. That's a clean conceptual link between thermodynamics and dynamics.\n\nThe soft spots are real. Eqs. (7)-(8) are the linearized equations for the decoupled difference combinations, but they are introduced with 'one can show' and no derivation. The decoupling from the total-charge mode, the metric, and the remaining scalar is asserted. This is the load-bearing part of the dynamical claim, so the omission matters. The numerics are also undocumented: no method, no convergence checks, no code. Figure 2 shows a mode crossing at κ=1, but I can't verify it. And the sign of σ11−σ12 is imported from Son-Starinets rather than computed here. That's a reputable source, but a direct check would close the argument.\n\nNone of this is fatal. The thermodynamic instability stands on its own, and the hydrodynamic argument makes the dynamical result expected. But as written, the paper states the QNM claim rather than proves it. A serious referee should ask for the derivation of (7)-(8) and for numerical details before the headline is accepted. This is a paper for holographers working on finite-density N=4 SYM; the thermodynamic part is a solid contribution, the dynamical part is a plausible conjecture with insufficient evidence. Send to peer review with those requests.","headline":"Thermodynamic instability claim is solid; the dynamical QNM claim needs derivation and numerics before I'd trust it.","tokens_in":12061,"tokens_out":2897,"would_cite":false,"duration_ms":28467,"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 low-temperature phase of N=4 SYM with three equal R-charge chemical potentials is not the AdS-Reissner-Nordström black brane; that background is thermodynamically and dynamically unstable.","keywords":["N=4 supersymmetric Yang-Mills","holographic duality","R-charge chemical potentials","thermodynamic stability","hydrodynamic instability","quasinormal modes","negative diffusion coefficient","AdS-Reissner-Nordström black brane"],"falsifier":"A direct derivation of equations (7) and (8) from the full linearized STU equations of motion, followed by an independent numerical solution of the complete coupled system without the decoupling ansatz, would settle the matter; if no quasinormal mode crosses into the upper half-plane for $\\kappa>1$, the instability claim fails. Alternatively, computing $\\sigma_{11}-\\sigma_{12}$ from first principles in this background and checking whether it changes sign before $\\kappa=1$ would test the hydrodynamic relation that turns thermodynamic instability into negative diffusion.","tokens_in":11033,"feed_emoji":"📉","tokens_out":11092,"duration_ms":84656,"temperature":0.7,"pith_summary":"N=4 supersymmetric Yang-Mills theory at strong coupling and finite R-charge density is studied through its five-dimensional gravitational dual. The paper argues that the charged black brane describing the equilibrium state becomes thermodynamically unstable at low temperature, and that for equal chemical potentials this thermodynamic instability is accompanied by a dynamical one: the R-charge diffusion coefficient turns negative and diffusive quasinormal modes move into the upper half-plane. Consequently, the low-temperature phase with equal chemical potentials is not described by the AdS-Reissner-Nordström black brane. The result matters because it shows that the standard holographic description of cold, dense SYM plasma breaks down, and it gives a concrete instance of the general relation between Hessian eigenvalues and hydrodynamic instabilities.","feed_headline":"Low-temperature N=4 SYM with equal R-charge densities is unstable","feed_subtitle":"R-charge diffusion goes negative, so the AdS-Reissner-Nordström black brane fails at low T.","key_machinery":"The five-dimensional STU black brane solution—a consistent truncation of IIB supergravity on S5 containing the metric, three U(1) gauge fields, and two neutral scalars—is the holographic background. Three pieces carry the argument: the thermodynamic Hessian $H^{\\epsilon}_{ij}=\\partial^2\\epsilon/\\partial y_i\\partial y_j$ and its eigenvectors, which for equal $\\kappa$ locate the critical point $\\kappa=1$; the decoupled fluctuation combinations $E^a_z$ and $s^a$, obtained by subtracting the average over the three charges, which reduce the linearized equations to the coupled ODEs (7) and (8); and the hydrodynamic formula $D=(\\sigma_{11}-\\sigma_{12})/(\\chi_{11}-\\chi_{12})$ for the diffusion of charge-difference modes. Imposing infalling boundary conditions at the horizon and normalizable ones at infinity (with the scalar log term set to zero), the numerical quasinormal mode solution gives a diffusive branch whose imaginary part becomes positive for $\\kappa>1$.","core_discovery":"The central claim is that the STU black brane with equal $\\kappa$ (all three chemical potentials equal) has a thermodynamic Hessian whose two coincident eigenvalues pass through zero at $\\kappa=1$ ($\\mu/2\\pi T=\\sqrt{2}$), making the state thermodynamically unstable for $\\kappa>1$. Relativistic hydrodynamics with three conserved charges predicts that the charge-difference fluctuations $\\delta(n_1-n_2)$, $\\delta(n_2-n_3)$, $\\delta(n_3-n_1)$ diffuse with coefficient $D=(\\sigma_{11}-\\sigma_{12})/(\\chi_{11}-\\chi_{12})$. Because $\\chi_{11}-\\chi_{12}$ changes sign at $\\kappa=1$ while $\\sigma_{11}-\\sigma_{12}$ stays positive, $D$ becomes negative. The paper verifies this by numerically solving the decoupled quasinormal mode equations for the relevant sound-channel fluctuations and finds that the diffusive mode crosses into the upper complex half-plane for $\\kappa>1$. Hence the AdS5-RN background is unstable, and the low-temperature phase of N=4 SYM with equal chemical potentials is not described by this background.","pith_inferences":["The instability at $\\kappa=1$ may mark the onset of a phase with broken U(1)^3 symmetry or spatially modulated order, similar to hairy-black-hole phases in other holographic models; the paper does not construct such a phase.","Because the diffusive mode returns to the lower half-plane at larger wavenumber, the effective long-wavelength description breaks down near the onset; higher-order gradient terms could select a finite-wavelength instability.","The Hessian-eigenvalue diagnostic used here could be applied to other strongly coupled theories with multiple conserved charges to predict which hydrodynamic modes go unstable before performing a full quasinormal mode analysis.","If the scalar boundary condition $A=0$ is the wrong choice, the quasinormal mode spectrum could shift; a comparison against an independent Kubo-formula computation of the diffusion constant would test this."],"forward_implications":["For equal chemical potentials, cooling N=4 SYM at finite density inevitably hits an instability before reaching low temperature, so the grand-canonical phase diagram must either end or pass through a phase transition.","The AdS5-RN black brane, commonly used to model strongly coupled plasma at finite density, is not the correct dual for the cold phase even in the symmetric three-charge configuration.","The Hessian eigenvectors identify the unstable modes as the charge-difference fluctuations $\\delta(n_1-n_2)$, $\\delta(n_2-n_3)$, and $\\delta(n_3-n_1)$; these are the fluctuations that acquire a negative diffusion coefficient.","The same thermodynamic-to-hydrodynamic correspondence is expected to hold for unequal chemical potentials, including the single-charge case where a hydrodynamic instability was seen earlier."],"supporting_citations":[{"why":"Supplies the five-dimensional STU black brane background with three gauge fields and two scalars on which the entire analysis is built.","marker":"[6]"},{"why":"Provides the hydrodynamics of R-charged black holes, including the conductivity matrix whose combination $\\sigma_{11}-\\sigma_{12}$ is imported to fix the sign of the diffusion coefficient.","marker":"[14]"},{"why":"Establishes the correspondence between thermodynamic and dynamical instabilities for charged AdS black holes, the template for the analysis here.","marker":"[12]"},{"why":"Finds the hydrodynamic instability for a single chemical potential, the result this paper extends to the equal-charge case.","marker":"[15]"},{"why":"Supplies the thermodynamic stability criteria via positive-definiteness of the energy/entropy Hessian.","marker":"[27]"},{"why":"Gives the quasinormal mode formalism and boundary conditions used to compute the holographic poles.","marker":"[31]"},{"why":"Provides the Minkowski-space AdS/CFT recipe that fixes the infalling boundary condition at the horizon.","marker":"[32]"}],"fun_headline_variants":["Equal chemical potentials trigger N=4 SYM instability at low T","N=4 SYM with equal chemical potentials: unstable at low T","Low-T N=4 SYM with equal charges: AdS-RN fails","N=4 SYM at equal charge density turns unstable at low temperature","AdS-RN black brane fails for N=4 SYM with equal charges at low T"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The linearized fluctuation equations (7) and (8), which the paper states with 'one can show' rather than deriving, are the correct decoupled equations for the R-charge and scalar modes, and the boundary condition that removes the logarithmic term of the scalar at the AdS boundary is the correct holographic choice.","fun_headline_variants_meta":{"raw":{"variants":["Equal chemical potentials trigger N=4 SYM instability at low T","N=4 SYM with equal chemical potentials: unstable at low T","Low-T N=4 SYM with equal charges: AdS-RN fails","N=4 SYM at equal charge density turns unstable at low temperature","AdS-RN black brane fails for N=4 SYM with equal charges at low T"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000885,"raw_usage":{"total_tokens":3859,"prompt_tokens":1018,"completion_tokens":2841,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":634,"completion_tokens_details":{"reasoning_tokens":2752}},"tokens_in":634,"tokens_out":2841,"duration_ms":16894,"temperature":1.0,"reasoning_tokens":2752,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T14:09:59.333814+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct derivation of equations (7) and (8) from the full linearized STU equations of motion, followed by an independent numerical solution of the complete coupled system without the decoupling ansatz, would settle the matter; if no quasinormal mode crosses into the upper half-plane for $\\kappa>1$, the instability claim fails. Alternatively, computing $\\sigma_{11}-\\sigma_{12}$ from first principles in this background and checking whether it changes sign before $\\kappa=1$ would test the hydrodynamic relation that turns thermodynamic instability into negative diffusion.","supporting_citations":[{"cited_title":"Callen, Thermodynamics and introduction to thermostatistics , Wiley (1985)","cited_arxiv_id":null,"evidence_quote":"Supplies the thermodynamic stability criteria via positive-definiteness of the energy/entropy Hessian."}],"review_version":1}