pith. sign in

arxiv: 2408.07937 · v3 · pith:6HBOD575new · submitted 2024-08-15 · ❄️ cond-mat.soft · cond-mat.stat-mech

Spectral Decomposition of Liquid Viscosity into Instantaneous Normal Modes

Long-Zhou Huang , Bingyu Cui , Min-Qiang Jiang , Matteo Baggioli , Yun-Jiang Wang This is my paper

Pith reviewed 2026-05-23 22:26 UTC · model grok-4.3

classification ❄️ cond-mat.soft cond-mat.stat-mech
keywords liquid viscosityinstantaneous normal modesmode-coupling temperaturenonaffine responseunstable localized modespotential energy landscapeviscous dynamics
0
0 comments X

The pith

Viscosity of liquids decomposes into sums of contributions from individual instantaneous normal modes, with unstable localized modes dominating above the mode-coupling temperature.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper applies nonaffine linear response theory to express the viscosity of metallic and model liquids as a sum over contributions from each instantaneous normal mode. Above the mode-coupling temperature, the unstable localized modes provide the main resistance to flow and serve as starting points for momentum diffusion. Below that temperature a crossover occurs and stable modes take over, matching the shift from saddle points to energy minima in the underlying landscape. The authors also give a quantitative relation that links the viscosity to the unstable modes in both the Arrhenius and non-Arrhenius regimes. The result supplies an atomic-scale accounting of what produces viscous drag.

Core claim

Viscosity equals a sum of individual instantaneous-normal-mode contributions obtained from nonaffine linear response; above T_MC the unstable localized modes control the value and act as precursors for diffusive momentum transport, while below T_MC stable modes dominate in agreement with the change from saddle-dominated to minima-dominated dynamics; a closed-form model connects viscosity to the unstable localized modes across both temperature regimes.

What carries the argument

Nonaffine linear-response decomposition of viscosity into instantaneous normal modes, isolating the role of unstable localized INMs.

If this is right

  • Viscosity can be predicted by counting or characterizing the unstable localized modes rather than running full flow simulations.
  • The dynamical crossover at T_MC accounts for the change between non-Arrhenius and Arrhenius temperature dependence of viscosity.
  • Transport properties become directly traceable to the topology of the potential-energy landscape (saddles versus minima).
  • The same mode decomposition applies to both metallic liquids and simple model liquids.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The decomposition could be used to engineer liquids whose viscosity is tuned by altering the stability of specific modes.
  • The same framework may clarify how the glass transition arrests flow once minima-dominated dynamics take over.
  • Application to molecular liquids such as water would test whether the unstable-mode mechanism is universal.

Load-bearing premise

The nonaffine linear response framework isolates viscosity contributions from each instantaneous normal mode without appreciable cross-mode interference or higher-order corrections.

What would settle it

A calculation in which the summed single-mode viscosity contributions differ measurably from the directly simulated total viscosity at the same state point, or in which the temperature of the ULINM-to-stable-mode crossover fails to coincide with the independently determined mode-coupling temperature.

Figures

Figures reproduced from arXiv: 2408.07937 by Bingyu Cui, Long-Zhou Huang, Matteo Baggioli, Min-Qiang Jiang, Yun-Jiang Wang.

Figure 1
Figure 1. Figure 1: FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2 [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4 [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6 [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗
read the original abstract

Viscosity, the resistance of a liquid to flow, is driven by atomic-scale friction but its microscopic origin remains poorly understood. We use a theoretical framework based on nonaffine linear response to decompose the viscosity of metallic and model liquids into contributions from individual instantaneous normal modes (INMs). Our approach reveals excellent agreement with simulations and exposes the specific excitations that govern viscous dynamics. Above the mode-coupling temperature ($T_{\text{MC}}$), viscosity is controlled by unstable localized INMs (ULINMs), which act as precursors for diffusive momentum transport. Below $T_{\text{MC}}$, we find a dynamical crossover where stable modes govern viscosity, a behavior consistent with a transition in the potential energy landscape from saddle-dominated to minima-dominated dynamics. We also propose a quantitative model connecting viscosity with ULINMs in both Arrhenius and non-Arrhenius regimes. This work provides a spectral decomposition of liquid viscosity, identifying the atomic modes responsible for it and opening a path to predict it from elementary excitations.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript develops a nonaffine linear-response framework to decompose the shear viscosity of metallic and model liquids into additive contributions from individual instantaneous normal modes (INMs). It reports that above the mode-coupling temperature T_MC viscosity is dominated by unstable localized INMs (ULINMs), which serve as precursors to diffusive momentum transport, while below T_MC stable modes take over, consistent with a saddle-to-minima transition in the potential-energy landscape. A quantitative model linking viscosity to ULINMs is proposed for both Arrhenius and non-Arrhenius regimes, with claimed excellent agreement to molecular-dynamics simulations.

Significance. If the per-mode decomposition is free of significant cross terms, the work supplies a concrete spectral picture of viscous flow in terms of elementary excitations and a direct link to the potential-energy landscape. The explicit identification of ULINMs as the controlling modes above T_MC, together with the proposed quantitative model, would constitute a substantive advance in the microscopic theory of liquid transport.

major comments (2)
  1. [Theoretical framework] Theoretical framework section: the nonaffine linear-response expression for viscosity is projected onto individual INM eigenvectors to obtain per-mode contributions. It is not shown that this projection commutes with the nonaffine correction or that the underlying stress autocorrelation is free of bilinear or higher-order cross-mode couplings; if such interference exists, the attribution of viscosity control specifically to ULINMs (and the claimed dynamical crossover at T_MC) does not follow.
  2. [Results] Results section (mode classification and quantitative model): the manuscript states excellent agreement with simulations and presents a quantitative model connecting viscosity to ULINMs, yet provides no explicit error bars, data-exclusion criteria, or quantitative definition of “localized” and “unstable” for the INM classification. Without these, the load-bearing claim that ULINMs govern viscosity above T_MC cannot be independently verified.
minor comments (2)
  1. [Abstract] The abstract and introduction use the acronym ULINM without an initial definition; a parenthetical expansion on first use would improve readability.
  2. [Figures] Figure captions should state the precise temperature range, system size, and ensemble used for each data set to allow direct comparison with the claimed agreement.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments, which help clarify key aspects of the work. We address each major comment below and will revise the manuscript to incorporate the requested clarifications and details.

read point-by-point responses

  1. Referee: [Theoretical framework] Theoretical framework section: the nonaffine linear-response expression for viscosity is projected onto individual INM eigenvectors to obtain per-mode contributions. It is not shown that this projection commutes with the nonaffine correction or that the underlying stress autocorrelation is free of bilinear or higher-order cross-mode couplings; if such interference exists, the attribution of viscosity control specifically to ULINMs (and the claimed dynamical crossover at T_MC) does not follow.

    Authors: We thank the referee for highlighting this point. The INMs form a complete orthogonal basis of eigenvectors of the instantaneous Hessian. The nonaffine correction enters the stress operator linearly, and the projection of the viscosity expression onto individual modes therefore yields additive contributions with vanishing cross terms by orthogonality. Nevertheless, we acknowledge that an explicit demonstration of the absence of bilinear couplings was not included. In the revised manuscript we will add a short derivation in the Theoretical framework section proving that cross-mode terms integrate to zero due to eigenvector orthogonality and the structure of the nonaffine stress, thereby confirming that the per-mode decomposition is exact within linear response and that the attribution to ULINMs remains valid. revision: yes

  2. Referee: [Results] Results section (mode classification and quantitative model): the manuscript states excellent agreement with simulations and presents a quantitative model connecting viscosity to ULINMs, yet provides no explicit error bars, data-exclusion criteria, or quantitative definition of “localized” and “unstable” for the INM classification. Without these, the load-bearing claim that ULINMs govern viscosity above T_MC cannot be independently verified.

    Authors: We agree that these elements are essential for independent verification. In the revised Results section we will: (i) add explicit statistical error bars to all reported viscosities and per-mode contributions, (ii) provide a quantitative definition of “unstable” (eigenvalue < 0) and “localized” (participation ratio below a stated numerical threshold, with the threshold value and justification given), and (iii) specify the data-exclusion criteria (e.g., number of independent configurations, equilibration protocol, and outlier rejection rules). These additions will make the classification and the claimed agreement with simulations fully reproducible. revision: yes

Circularity Check

0 steps flagged

No circularity; derivation uses external nonaffine framework without reduction to fitted inputs or self-citations

full rationale

The abstract and provided context describe a decomposition of viscosity via nonaffine linear response projected onto INM eigenvectors, with a proposed quantitative model linking viscosity to ULINMs. No equations are shown that define a parameter from data and then rename its output as an independent prediction. No self-citation chains or uniqueness theorems from the same authors are invoked as load-bearing. The central claim remains a forward application of an established response framework to identify mode contributions, which is self-contained against external benchmarks and does not reduce by construction to its inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 1 invented entities

The central claim rests on the validity of nonaffine linear response theory for viscosity and on the operational definition of ULINMs as unstable localized modes that act as diffusive precursors. No free parameters are explicitly named in the abstract. The framework introduces ULINMs as a functional entity whose independent falsifiability outside the viscosity fit is not stated.

axioms (1)
  • domain assumption Nonaffine linear response theory provides an exact decomposition of viscosity into instantaneous normal mode contributions.
    The entire spectral decomposition is based on this framework as stated in the abstract.
invented entities (1)
  • ULINMs (unstable localized instantaneous normal modes) no independent evidence
    purpose: To identify the specific excitations that control viscosity above T_MC and serve as precursors for diffusive momentum transport.
    Introduced in the abstract as the dominant contributors; no independent evidence (e.g., predicted observable outside viscosity) is provided.

pith-pipeline@v0.9.0 · 5715 in / 1614 out tokens · 22244 ms · 2026-05-23T22:26:56.829405+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

69 extracted references · 69 canonical work pages

  1. [1]

    D. S. Viswanath, T. K. Ghosh, D. H. Prasad, N. V. Dutt, and K. Y. Rani, Viscosity of liquids: theory, estimation, experiment, and data (Springer Science & Business Me- dia, 2007)

    work page 2007

  2. [2]

    L. D. Landau and E. M. Lifshitz, Fluid Mechanics (Sec- ond Edition) (Pergamon, 1987)

    work page 1987

  3. [3]

    J. M. Garcia, E. Chambers, Z. Matta, and M. Clark, Dysphagia 20, 325 (2005)

    work page 2005

  4. [4]

    Polini and A

    M. Polini and A. K. Geim, Phys. Today 73, 28 (2020)

    work page 2020

  5. [5]

    Shuryak, Prog

    E. Shuryak, Prog. Part. Nucl. Phys. 53, 273 (2004). 6

    work page 2004

  6. [6]

    M. S. Green, J. Chem. Phys. 22, 398 (1954)

    work page 1954

  7. [7]

    R. Kubo, J. Phys. Soc. Jpn. 12, 570 (1957)

    work page 1957

  8. [8]

    B. J. Alder and T. E. Wainwright, Phys. Rev. Lett. 18, 988 (1967)

    work page 1967

  9. [9]

    G. G. Stokes, Mathematical and Physical Papers , Cam- bridge Library Collection - Mathematics (Cambridge University Press, 2009)

    work page 2009

  10. [10]

    T. A. Broadbent, A. B. Einstein, H. Hertz, H. L. Dryden, F. P. Murnaghan, and H. Bateman, The Mathematical Gazette 41, 231 (1928)

    work page 1928

  11. [11]

    S. G. Moustafa, A. J. Schultz, and J. F. Douglas, J. Chem. Phys. 160, 024114 (2024)

    work page 2024

  12. [12]

    Trachenko, Theory of liquids: from excitations to ther- modynamics (Cambridge University Press, 2023)

    K. Trachenko, Theory of liquids: from excitations to ther- modynamics (Cambridge University Press, 2023)

    work page 2023

  13. [13]

    Moon, Heat Carriers in Liquids: An Introduction (Springer, 2024)

    J. Moon, Heat Carriers in Liquids: An Introduction (Springer, 2024)

    work page 2024

  14. [14]

    J. Moon, L. Lindsay, and T. Egami, Phys. Rev. E 108, 014601 (2023)

    work page 2023

  15. [15]

    Zwanzig, Phys

    R. Zwanzig, Phys. Rev. 156, 190 (1967)

    work page 1967

  16. [16]

    J. C. Maxwell, Phil. Trans. R. Soc. , 49 (1867)

    work page

  17. [17]

    Keyes, J

    T. Keyes, J. Phys. Chem. 101, 2921 (1997)

    work page 1997

  18. [18]

    Keyes, Normal Mode Analysis: Theory and Applica- tions to Biological and Chemical Systems (Chapman and Hall/CRC, 2005) Chap

    T. Keyes, Normal Mode Analysis: Theory and Applica- tions to Biological and Chemical Systems (Chapman and Hall/CRC, 2005) Chap. The Relation Between Unstable Instantaneous Normal Modes and Diffusion

    work page 2005

  19. [19]

    Rahman, M

    A. Rahman, M. Mandell, and J. McTague, J. Chem. Phys. 64, 1564 (1976)

    work page 1976

  20. [20]

    R. M. Stratt, Acc. Chem. Res. 28, 201 (1995)

    work page 1995

  21. [21]

    Keyes, J

    T. Keyes, J. Chem. Phys. 101, 5081 (1994)

    work page 1994

  22. [22]

    Keyes, J

    T. Keyes, J. Chem. Phys. 103, 9810 (1995)

    work page 1995

  23. [23]

    S. D. Bembenek and B. B. Laird, Phys. Rev. Lett. 74, 936 (1995)

    work page 1995

  24. [24]

    V. I. Clapa, T. Kottos, and F. W. Starr, J. Chem. Phys. 136, 144504 (2012)

    work page 2012

  25. [25]

    J. D. Gezelter, E. Rabani, and B. J. Berne, J. Chem. Phys. 107, 4618 (1997)

    work page 1997

  26. [26]

    Sciortino and P

    F. Sciortino and P. Tartaglia, Phys. Rev. Lett. 78, 2385 (1997)

    work page 1997

  27. [27]

    La Nave, A

    E. La Nave, A. Scala, F. W. Starr, F. Sciortino, and H. E. Stanley, Phys. Rev. Lett. 84, 4605 (2000)

    work page 2000

  28. [28]

    La Nave, A

    E. La Nave, A. Scala, F. W. Starr, H. E. Stanley, and F. Sciortino, Phys. Rev. E 64, 036102 (2001)

    work page 2001

  29. [29]

    Chowdhary and T

    J. Chowdhary and T. Keyes, Phys. Rev. E 65, 026125 (2002)

    work page 2002

  30. [30]

    Li and T

    W.-X. Li and T. Keyes, J. Chem. Phys. 107, 7275 (1997)

    work page 1997

  31. [31]

    Stratt, Int

    R. Stratt, Int. J. Thermophys 18, 899–907 (1997)

    work page 1997

  32. [32]

    Zaccone, Phys

    A. Zaccone, Phys. Rev. E 108, 044101 (2023)

    work page 2023

  33. [33]

    Lemaˆ ıtre and C

    A. Lemaˆ ıtre and C. Maloney, J. Stat. Phys. 123, 415 (2006)

    work page 2006

  34. [34]

    Cui and A

    B. Cui and A. Zaccone, Phys. Rev. E 97, 060102 (2018)

    work page 2018

  35. [35]

    B. Cui, J. Yang, J. Qiao, M. Jiang, L. Dai, Y.-J. Wang, and A. Zaccone, Phys. Rev. B 96, 094203 (2017)

    work page 2017

  36. [36]

    Carof, R

    A. Carof, R. Vuilleumier, and B. Rotenberg, J. Chem. Phys. 140, 124103 (2014)

    work page 2014

  37. [37]

    Lesnicki, R

    D. Lesnicki, R. Vuilleumier, A. Carof, and B. Rotenberg, Phys. Rev. Lett. 116, 147804 (2016)

    work page 2016

  38. [38]

    G. Jung, M. Hanke, and F. Schmid, J. Chem. Theory Comput. 13 6, 2481 (2017)

    work page 2017

  39. [39]

    M. I. Mendelev, Y. Sun, F. Zhang, C. Z. Wang, and K. M. Ho, J. Chem. Phys. 151, 214502 (2019)

    work page 2019

  40. [40]

    Kob and H

    W. Kob and H. C. Andersen, Phys. Rev. E 51, 4626 (1995)

    work page 1995

  41. [41]

    Kob and H

    W. Kob and H. C. Andersen, Phys. Rev. E 52, 4134 (1995)

    work page 1995

  42. [42]

    Zaccone and M

    A. Zaccone and M. Baggioli, PNAS. 118, e2022303118 (2021)

    work page 2021

  43. [43]

    Z¨ urcher and T

    U. Z¨ urcher and T. Keyes, Phys. Rev. E55, 6917 (1997)

    work page 1997

  44. [44]

    Stamper, D

    C. Stamper, D. Cortie, Z. Yue, X. Wang, and D. Yu, J. Phys. Chem. Lett. 13, 3105 (2022)

    work page 2022

  45. [45]

    S. Jin, X. Fan, C. Stamper, R. A. Mole, Y. Yu, L. Hong, D. Yu, and M. Baggioli, Scientific Reports 14, 18805 (2024)

    work page 2024

  46. [46]

    Vogel, Phys

    H. Vogel, Phys. Z 22, 645 (1921)

    work page 1921

  47. [47]

    G. S. Fulcher, J. Am. Ceram. Soc. 8, 339 (1925)

    work page 1925

  48. [48]

    Tammann and W

    G. Tammann and W. Hesse, Z. Anorg. Allg. Chem 156, 245 (1926)

    work page 1926

  49. [49]

    S. D. Bembenek and B. B. Laird, J. Chem. Phys. 104, 5199 (1996)

    work page 1996

  50. [50]

    Zhang, J

    W. Zhang, J. F. Douglas, and F. W. Starr, J. Chem. Phys. 151, 184904 (2019)

    work page 2019

  51. [51]

    Widmer-Cooper, H

    A. Widmer-Cooper, H. Perry, P. Harrowell, and D. R. Reichman, Nat. Phys. 4, 711 (2008)

    work page 2008

  52. [52]

    Anomalies

    H. Mizuno and A. Ikeda, in Low-Temperature Thermal and Vibrational Properties of Disordered Solids: A Half- Century of Universal “Anomalies” of Glasses (World Sci- entific, 2023) pp. 375–433

    work page 2023

  53. [53]

    Angelani, R

    L. Angelani, R. Di Leonardo, G. Ruocco, A. Scala, and F. Sciortino, Phys. Rev. Lett. 85, 5356 (2000)

    work page 2000

  54. [54]

    Broderix, K

    K. Broderix, K. K. Bhattacharya, A. Cavagna, A. Zip- pelius, and I. Giardina, Phys. Rev. Lett. 85, 5360 (2000)

    work page 2000

  55. [55]

    Adam and J

    G. Adam and J. H. Gibbs, J. Chem. Phys.43, 139 (1965)

    work page 1965

  56. [56]

    I. M. Douglass, J. C. Dyre, and L. Costigliola, Phys. Rev. Lett. 133, 068001 (2024)

    work page 2024

  57. [57]

    Cao and Y.-J

    L.-L. Cao and Y.-J. Wang, J. Phys. Chem. Lett. 15, 811 (2024)

    work page 2024

  58. [58]

    Bouchaud and G

    J.-P. Bouchaud and G. Biroli, J. Chem. Phys. 121, 7347 (2004)

    work page 2004

  59. [59]

    Ozawa, C

    M. Ozawa, C. Scalliet, A. Ninarello, and L. Berthier, J. Chem. Phys. 151 (2019)

    work page 2019

  60. [60]

    L. D. Landau and E. M. Lifshitz, Statistical Physics: Vol- ume 5 , Vol. 5 (Pergamon, 1980)

    work page 1980

  61. [61]

    A. S. Keys, L. O. Hedges, J. P. Garrahan, S. C. Glotzer, and D. Chandler, Phys. Rev. X 1, 021013 (2011)

    work page 2011

  62. [62]

    M. P. Ciamarra, W. Ji, and M. Wyart, PNAS. 121, e2400611121 (2024)

    work page 2024

  63. [63]

    E. d. C. Andrade, Nature 125, 309 (1930)

    work page 1930

  64. [64]

    M. Mohr, R. K. Wunderlich, S. Koch, P. K. Galenko, A. K. Gangopadhyay, K. F. Kelton, J.-Z. Jiang, and H.- J. Fecht, Microgravity. Sci. Tec. 31, 177 (2019)

    work page 2019

  65. [65]

    Mauro, M

    N. Mauro, M. Blodgett, M. Johnson, A. Vogt, and K. Kelton, Nat. Commun. 5, 4616 (2014)

    work page 2014

  66. [66]

    Nos´ e, J

    S. Nos´ e, J. Chem. Phys.81, 511 (1984)

    work page 1984

  67. [67]

    https://www.lammps.org/

    work page

  68. [68]

    Hill, Proc

    R. Hill, Proc. Phys. Soc. A 65, 349 (1952)

    work page 1952

  69. [69]

    Gottwald, S

    F. Gottwald, S. Karsten, S. D. Ivanov, and O. K¨ uhn, J. Chem. Phys. 142, 244110 (2015). END MATTER Appendix A: Extended analysis Cu-Zr metallic liquid - In the main text, we have shown the results for the vis- cosity η computed using the Green-Kubo formalism, Eq. (6). Alternatively, the viscosity can be obtained from non-equilibrium molecular dynamics us...

    work page 2015