REVIEW 4 major objections 5 minor 4 cited by
Instability in ${\cal N}=4$ supersymmetric Yang-Mills theory at finite density
T0 review · 4 major / 5 minor · reviewed 2026-08-11 · deepseek-v4-flash
Pith's one-line read 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.
desk verdict Thermodynamic instability claim is solid; the dynamical QNM claim needs derivation and numerics before I'd trust it. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
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$.
What would settle it
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.
Extended reading notes
Core claim
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.
Load-bearing premise
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.
Editorial extensions
If this is right
- 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.
Reading between the lines
- 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.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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.
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 (4)
- [Section 4, Eqs. (7)-(8)] 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 4, boundary conditions at u=0] 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 4, Fig. 2 and numerical procedure] 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 4 and Appendix B, Eq. (19)] 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.
minor comments (5)
- [Abstract and Introduction] The name 'Behrnd' should be 'Behrndt' to match the reference list and the original paper.
- [Section 3, text after Eq. (5)] 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.
- [Fig. 2] 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.
- [Appendix A, text before Eq. (9)] The relation 'G5 = πL3/2Nc2' is ambiguous; it should be written as G5 = π L^3/(2 N_c^2) to avoid confusion.
- [References] 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.
Circularity Check
No significant circularity: the thermodynamic-instability input is a standard equation of state, and the quasinormal-mode calculation independently confirms the hydrodynamic prediction rather than reducing to it by construction.
full rationale
The derivation chain is not circular. The equation of state (2) is taken from standard black-hole thermodynamics, with refs. [30] and [14] cited as sources; it is an input, not a quantity derived from the instability being claimed. The Hessian analysis in Sec. 3 is a direct computation from this equation of state, and the eigenvector directions are used only to suggest which fluctuations to examine. The central dynamical result comes from solving the quasinormal-mode equations (7) and (8), which are stated with the phrase "one can show" but are presented as equations of motion following from the action (9); the numerical solution then finds the diffusive mode crossing into the upper half-plane for kappa>1 without imposing that outcome. The hydrodynamic relation D = (sigma_11 - sigma_12)/(chi_11 - chi_12) is derived in Appendix B and used as a prediction, not as a fitted result, and the QNM calculation agrees with it independently. The citation of ref. [14], by one of the present authors, for the equation of state and for the sign of sigma_11 - sigma_12 is a citation to a published, independently checkable calculation, so it is real evidence rather than circular reliance. The "one can show" assertion and the imported sign of sigma_11 - sigma_12 are rigor gaps or external inputs, not reductions of the claim to its own assumptions; no step makes the conclusion equivalent to its premise by definition.
Assumptions & free parameters
assumptions (6)
- domain assumption The STU model is a consistent truncation of type IIB supergravity on S5 and is the correct holographic dual of N=4 SYM at finite R-charge density.
- domain assumption Equation of state (2) correctly describes the thermodynamics of the STU black brane.
- domain assumption The hydrodynamic constitutive relations in Appendix B, including the Landau-Lifshitz frame and the diffusion operator D_ab = χ_ab ∂t − σ_ab ∂^2, apply to the R-charge currents.
- domain assumption In the state with μa=μ, the susceptibility and conductivity matrices have equal diagonal and equal off-diagonal entries, and all charge densities are equal.
- ad hoc to paper The fluctuation equations (7) and (8) correctly describe decoupled perturbations, and the scalar boundary condition A=0 at u=0 is correct.
- ad hoc to paper The numerical integration used to find the quasinormal frequencies is accurate and the diffusion mode is correctly identified.
Cite this review
Pith. "Pith review of Instability in ${\cal N}=4$ supersymmetric Yang-Mills theory at finite density." pith.science (2026). https://pith.science/paper/LQU6P53G
@misc{pith2026241212353,
author = {Pith},
title = {Pith review of: Instability in $\cal N=4$ supersymmetric Yang-Mills theory at finite density},
year = {2026},
howpublished = {\url{https://pith.science/paper/LQU6P53G}},
note = {Machine review of arXiv:2412.12353}
}
abstract
Equilibrium states of ${\cal N}=4$ supersymmetric Yang-Mills theory can be characterized by the temperature and three chemical potentials, corresponding to the ${\rm U}(1)^3$ subgroup of the $R$-symmetry group. We investigate the phase diagram of the theory at strong coupling, in the grand canonical ensemble in flat space, using its holographic description via the five-dimensional model of Behrnd, Cveti\v{c}, and Sabra. The bulk action includes the metric, three Abelian gauge fields, and two neutral scalar fields. The equilibrium state described by the charged black brane is always thermodynamically unstable at low temperature. Relativistic hydrodynamics with multiple conserved charges predicts that thermodynamic instability is accompanied by a dynamical instability, with the eigenvalues and eigenvectors of the corresponding Hessian playing a key role in identifying the unstable modes. We explicitly demonstrate this for three equal chemical potentials, finding unstable quasinormal modes that describe $R$-charge diffusion. Consequently, the low-temperature phase of ${\cal N}=4$ supersymmetric Yang-Mills theory with equal chemical potentials is not described by the AdS-Reissner-Nordstr\"om black brane.
Figures
Forward citations
Cited by 4 Pith papers
-
Chiral anomaly from (anomalous) spin hydrodynamics
Spinning D3-brane hydrodynamics in 10d maps to anomalous chiral hydrodynamics in 4d, identifying the R-current anomaly as a 10d gravitational anomaly.
-
Instability of Baryonic Black Branes
Baryonic black branes in the conifold gauge theory are shown to be dynamically unstable below T/µ_B = 0.2770(5), with an accompanying subdominant ordered phase extending to high temperatures.
-
Quark-Antiquark Potential as a Probe for Holographic Phase Transitions
A holographic Wilson loop calculation shows the quark-antiquark potential has a crossing at the phase transition point ψ=1 for two AdS5 black hole phases.
-
The ordered phase of charged N=4 SYM plasma
A new ordered phase of charged N=4 SYM plasma, with a dimension-2 condensate O2 proportional to T^2 at high temperature, is constructed holographically and dominates the microcanonical ensemble.
Reference graph
Works this paper leans on
-
[14]
D.T. Son and A.O. Starinets, Hydrodynamics of R-charged black holes , JHEP 03 (2006) 052 [hep-th/0601157]
arXiv 2006
-
[1]
Maldacena, The Large N limit of superconformal field theories and supergravity , Adv
J.M. Maldacena, The Large N limit of superconformal field theories and supergravity , Adv. Theor. Math. Phys. 2 (1998) 231 [ hep-th/9711200]
arXiv 1998
-
[2]
S.S. Gubser, I.R. Klebanov and A.M. Polyakov, Gauge theory correlators from noncritical string theory , Phys. Lett. B 428 (1998) 105 [ hep-th/9802109]
arXiv 1998
-
[3]
Witten, Anti-de Sitter space and holography , Adv
E. Witten, Anti-de Sitter space and holography , Adv. Theor. Math. Phys. 2 (1998) 253 [hep-th/9802150]
arXiv 1998
-
[4]
D. Yamada and L.G. Yaffe, Phase diagram of N = 4 super-Yang-Mills theory with R-symmetry chemical potentials , JHEP 09 (2006) 027 [ hep-th/0602074]
arXiv 2006
-
[5]
O. Aharony, S.S. Gubser, J.M. Maldacena, H. Ooguri and Y. Oz, Large N field theories, string theory and gravity , Phys. Rept. 323 (2000) 183 [ hep-th/9905111]. 10
arXiv 2000
-
[6]
K. Behrndt, M. Cvetic and W.A. Sabra, Nonextreme black holes of five-dimensional N=2 AdS supergravity, Nucl. Phys. B 553 (1999) 317 [ hep-th/9810227]
arXiv 1999
-
[7]
M. Cvetic, M.J. Duff, P. Hoxha, J.T. Liu, H. Lu, J.X. Lu et al., Embedding AdS black holes in ten-dimensions and eleven-dimensions , Nucl. Phys. B 558 (1999) 96 [hep-th/9903214]
arXiv 1999
Show all 41 references
-
[8]
Cvetic, H
M. Cvetic, H. Lu, C.N. Pope, A. Sadrzadeh and T.A. Tran, Consistent SO(6) reduction of type IIB supergravity on S5, Nucl. Phys. B 586 (2000) 275 [ hep-th/0003103]
2000 arXiv
-
[9]
Cvetic and S.S
M. Cvetic and S.S. Gubser, Phases of R charged black holes, spinning branes and strongly coupled gauge theories , JHEP 04 (1999) 024 [ hep-th/9902195]
1999 arXiv
-
[10]
Cvetic and S.S
M. Cvetic and S.S. Gubser, Thermodynamic stability and phases of general spinning branes, JHEP 07 (1999) 010 [ hep-th/9903132]
1999 arXiv
-
[11]
Cai and K.-S
R.-G. Cai and K.-S. Soh, Critical behavior in the rotating D-branes , Mod. Phys. Lett. A 14 (1999) 1895 [ hep-th/9812121]
1999 arXiv
-
[12]
Gubser and I
S.S. Gubser and I. Mitra, Instability of charged black holes in Anti-de Sitter space , Clay Math. Proc. 1 (2002) 221 [ hep-th/0009126]
2002 arXiv
-
[13]
Gubser and I
S.S. Gubser and I. Mitra, The Evolution of unstable black holes in anti-de Sitter space , JHEP 08 (2001) 018 [ hep-th/0011127]
2001 arXiv
-
[15]
Buchel, Critical phenomena in N = 4 SYM plasma , Nucl
A. Buchel, Critical phenomena in N = 4 SYM plasma , Nucl. Phys. B 841 (2010) 59 [1005.0819]
2010 arXiv
-
[16]
Chamblin, R
A. Chamblin, R. Emparan, C.V. Johnson and R.C. Myers, Charged AdS black holes and catastrophic holography, Phys. Rev. D 60 (1999) 064018 [ hep-th/9902170]
1999 arXiv
-
[17]
Gubser, Breaking an Abelian gauge symmetry near a black hole horizon , Phys
S.S. Gubser, Breaking an Abelian gauge symmetry near a black hole horizon , Phys. Rev. D 78 (2008) 065034 [ 0801.2977]
2008 arXiv
-
[18]
O.J.C. Dias, P. Figueras, R. Monteiro, J.E. Santos and R. Emparan, Instability and new phases of higher-dimensional rotating black holes , Phys. Rev. D 80 (2009) 111701 [0907.2248]
2009 arXiv
-
[19]
Bhattacharyya, S
S. Bhattacharyya, S. Minwalla and K. Papadodimas, Small Hairy Black Holes in AdS5 × S5, JHEP 11 (2011) 035 [ 1005.1287]
2011 arXiv
-
[20]
O.J.C. Dias, P. Figueras, R. Monteiro, H.S. Reall and J.E. Santos, An instability of higher-dimensional rotating black holes , JHEP 05 (2010) 076 [ 1001.4527]
2010 arXiv
-
[21]
O.J.C. Dias, P. Figueras, R. Monteiro and J.E. Santos, Ultraspinning instability of anti-de Sitter black holes , JHEP 12 (2010) 067 [ 1011.0996]
2010 arXiv
-
[22]
Markeviciute and J.E
J. Markeviciute and J.E. Santos, Hairy black holes in AdS5 × S5, JHEP 06 (2016) 096 [1602.03893]. 11
2016 arXiv
-
[23]
Yaffe, Large N phase transitions and the fate of small Schwarzschild-AdS black holes, Phys
L.G. Yaffe, Large N phase transitions and the fate of small Schwarzschild-AdS black holes, Phys. Rev. D 97 (2018) 026010 [ 1710.06455]
2018 arXiv
-
[24]
O.J.C. Dias, P. Mitra and J.E. Santos, New phases of N = 4 SYM at finite chemical potential, JHEP 05 (2023) 053 [ 2207.07134]
2023 arXiv
-
[25]
S. Choi, D. Jain, S. Kim, V. Krishna, E. Lee, S. Minwalla et al., Dual Dressed Black Holes as the end point of the Charged Superradiant instability in N = 4 Yang Mills, 2409.18178
-
[26]
O.J.C. Dias, P. Mitra and J.E. Santos, Charged Rotating Hairy Black Holes in AdS5 × S5: Unveiling their Secrets , 2411.18712
-
[27]
Callen, Thermodynamics and introduction to thermostatistics , Wiley (1985)
H.B. Callen, Thermodynamics and introduction to thermostatistics , Wiley (1985)
1985
-
[28]
Kovtun, Lectures on hydrodynamic fluctuations in relativistic theories , J
P. Kovtun, Lectures on hydrodynamic fluctuations in relativistic theories , J. Phys. A45 (2012) 473001 [ 1205.5040]
2012 arXiv
-
[29]
Buchel, A Holographic perspective on Gubser-Mitra conjecture, Nucl
A. Buchel, A Holographic perspective on Gubser-Mitra conjecture, Nucl. Phys. B 731 (2005) 109 [ hep-th/0507275]
2005 arXiv
-
[30]
Harmark and N.A
T. Harmark and N.A. Obers, Thermodynamics of spinning branes and their dual field theories, JHEP 01 (2000) 008 [ hep-th/9910036]
2000 arXiv
-
[31]
Kovtun and A.O
P.K. Kovtun and A.O. Starinets, Quasinormal modes and holography , Phys. Rev. D72 (2005) 086009 [ hep-th/0506184]
2005 arXiv
-
[32]
Son and A.O
D.T. Son and A.O. Starinets, Minkowski space correlators in AdS / CFT correspondence: Recipe and applications, JHEP 09 (2002) 042 [ hep-th/0205051]
2002 arXiv
-
[33]
Herzog, Lectures on Holographic Superfluidity and Superconductivity, J
C.P. Herzog, Lectures on Holographic Superfluidity and Superconductivity, J. Phys. A 42 (2009) 343001 [ 0904.1975]
2009 arXiv
-
[34]
Friess, S.S
J.J. Friess, S.S. Gubser and I. Mitra, Counter-examples to the correlated stability conjecture, Phys. Rev. D 72 (2005) 104019 [ hep-th/0508220]
2005 arXiv
-
[35]
Buchel and C
A. Buchel and C. Pagnutti, Correlated stability conjecture revisited, Phys. Lett. B 697 (2011) 168 [ 1010.5748]
2011 arXiv
-
[36]
Buchel and A
A. Buchel and A. Patrushev, Can the correlated stability conjecture be saved? , JHEP 06 (2011) 090 [ 1102.5331]
2011 arXiv
-
[37]
Buchel, Fate of the conformal order , Phys
A. Buchel, Fate of the conformal order , Phys. Rev. D 103 (2021) 026008 [ 2011.11509]
2021 arXiv
-
[38]
Gout´ eraux and E
B. Gout´ eraux and E. Mefford,Linear dynamical stability and the laws of thermodynamics, 2407.07939
-
[39]
Gladden, V
L. Gladden, V. Ivo, P.K. Kovtun and A.O. Starinets, To appear, (2025)
2025
-
[40]
Freedman, S.D
D.Z. Freedman, S.D. Mathur, A. Matusis and L. Rastelli, Correlation functions in the CFT(d) / AdS(d+1) correspondence, Nucl. Phys. B 546 (1999) 96 [ hep-th/9804058]
1999 arXiv
-
[41]
Landau and E.M
L.D. Landau and E.M. Lifshitz, Fluid Mechanics, Pergamon Press (1987). 12
1987
Reviewed August 11, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.