[论文解读] Soft modes in jammed hard spheres (I): Mean field theory of the isostatic transition

Giorgio Parisi|arXiv (Cornell University)|Jan 17, 2014

Advanced Thermodynamics and Statistical Mechanics被引用 7

一句话总结

本文为受 jammed 硬球中软模态发展了一套平均场理论，从理论上证明在等价刚性转变点（配位数 z = 2D）处，振动态密度在低频区域表现出 ∝ λ⁻¹/² 的平方根发散，表明存在额外的软模态。在无限维极限下，得到精确解，证实了 λ⁻¹/² 奇异性以及在 ω = 0 附近存在准局域态。

ABSTRACT

In this paper we consider different models for soft modes in jammed hard spheres. We show how one can construct mean field models that can be solved analytically. The analytic solution of these models displays an excess of low energy soft modes that become more and more localized by decreasing the energy. A simple solution of these models is found in the infinite dimensional limit.

研究动机与目标

• 通过可解析求解的平均场模型，理解 jammed 硬球中低能振动模态过剩的起源。
• 研究在机械稳定性仅勉强满足的等价刚性转变点（z = 2D）处的临界行为。
• 确定软模态是否在低能区域趋于局域化，以及这如何影响谱密度和预解函数行为。
• 推导无限维极限下的精确结果，并评估密度谱中 λ⁻¹/² 奇异性在该极限下的鲁棒性。

提出的方法

• 构建一个基于伪 Wishart 矩阵的平均场模型，其中相互作用矩阵 J 定义为 J_{i,k} = ∑_ν x^ν_i x^ν_k，其中向量 x^ν 为高斯分布。
• 分析在 N → ∞ 极限下矩阵 J 的谱密度 ρ(λ)，表明其服从 Marcenko-Pastur 分布，当 α = M/N = α_c 时，在 λ = 0 处出现 λ⁻¹/² 发散。
• 使用预解函数 R(λ) = (λ - J)⁻¹，并定义 I(λ) = lim_{ε→0} ⟨(Im R_{i,i}(λ+iε))²⟩ / ⟨Im R_{i,i}(λ+iε)⟩² 以探测局域化相变。
• 应用腔方法和在 D → ∞ 极限下预解函数的自洽方程，推导出 s = -s q(-sλ, γ) + 1，其中 s = 1 + m，且在 γ = 2 时 s ∝ λ⁻¹/²。
• 考虑一个投影预解函数 Q(λ) 以消除平动零模态，从而能够分析对单极力的响应，并探测在 λ = 0 附近的大幅涨落。
• 将结果与 DD 模型的数值数据进行比较，发现 ω_c ∝ (z - z_c)^ν 的关系中 ν ≈ 0.94，符合良好。

实验结果

研究问题

• RQ1 在 jammed 硬球的平均场模型中，低能振动模态过剩的起源是什么？
• RQ2 在等价刚性点（z = 2D）处，谱密度 ρ(λ) 在 λ = 0 附近的特性如何？
• RQ3 系统是否在 ω = 0 处表现出局域化相变？如果是，它在预解函数中如何体现？
• RQ4 动力学矩阵的本征态在低能区域如何行为？是否存在局域化趋势？
• RQ5 在无限维极限下，预解函数及其对角元的特性如何？

主要发现

• 在等价刚性点（α = α_c）处，谱密度 ρ(λ) 以 λ⁻¹/² 发散，导致在 ω = 0 附近呈现平坦的态密度 D(ω) ∝ ω⁰。
• ρ(λ) 在 λ = 0 处的发散具有鲁棒性，且在 D → ∞ 极限下依然存在，与 Marcenko-Pastur 分布中观察到的相同 λ⁻¹/² 奇异性一致。
• 预解函数对角元在 λ → 0 时表现出发散涨落 I(λ)，标志在零频率附近出现准局域扩展态。
• 在无限维极限下，自洽方程给出在 γ = 2 时 s ∝ λ⁻¹/²，从而从理论上确认了 λ⁻¹/² 行为。
• 对负 λ 的投影预解函数 Q(λ) 表现出与小正 λ 时 Im(R(λ)) 类似的大幅涨落，表明在 ω = 0 附近响应增强。
• 该模型预测 ω_c ∝ (z - z_c)^ν，其中 ν ≈ 0.94，与 DD 模型的数值数据高度一致。