[논문 리뷰] Mean-field density of states of a small-world model and a jammed soft spheres model

M. Pernici|arXiv (Cornell University)|2020. 01. 08.

Material Dynamics and Properties인용 수 5

한 줄 요약

이 논문은 조임점 근처 소프트 스피어 시스템의 진동 밀도 상태(DOS)를 위한 평균장 랜덤 블록 행렬 모델을 제안한다. 이 모델은 랜덤 정규 그래프와 소월드 단절을 조합하여 구성되며, $ d=3 $ 에서 $ t = Z/d $ 파라미터를 가진 이동된 Kesten-McKay DOS가 보존 피크와 저주파수 플랫폼을 정확히 기록함을 보여준다. 이는 $ Z=7,8 $ 에서 분자 동역학 시뮬레이션과 일치하며, 무작위 탄성 네트워크 모델보다 개선된 이소스태틱 플랫폼을 재현한다.

ABSTRACT

We consider a class of random block matrix models in $d$ dimensions, $d \\ge 1$, motivated by the study of the vibrational density of states (DOS) of soft spheres near the isostatic point. The contact networks of average degree $Z = z_0 + \\zeta$ are represented by random $z_0$-regular graphs (only the circle graph in $d=1$ with $z_0=2$) to which Erd\\"os-Renyi graphs having a small average degree $\\zeta$ are superimposed. In the case $d=1$, for $\\zeta$ small the shifted Kesten-McKay DOS with parameter $Z$ is a mean-field solution for the DOS. Numerical simulations in the $z_0=2$ model, which is the $k=1$ Newman-Watts small-world model, and in the $z_0=3$ model lead us to conjecture that for $\\zeta \ o 0$ the cumulative function of the DOS converges uniformly to that of the shifted Kesten-McKay DOS, in an interval $[0, \\omega_0]$, with $\\omega_0 < \\sqrt{z_0-1} + 1$. For $2 \\le d \\le 4$, we introduce a cutoff parameter $K_d \\le 0.5$ modeling sphere repulsion. The case $K_d=0$ is the random elastic network case, with the DOS close to the Marchenko-Pastur DOS with parameter $t=\\frac{Z}{d}$. For $K_d$ large the DOS is close for small $\\omega$ to the shifted Kesten-McKay DOS with parameter $t=\\frac{Z}{d}$; in the isostatic case the DOS has around $\\omega=0$ the expected plateau. The boson peak frequency in $d=3$ with $K_3$ large is close to the one found in molecular dynamics simulations for $Z=7$ and $8$.

연구 동기 및 목표

• 소프트 스피어 시스템의 진동 밀도 상태(DOS)를 랜덤 블록 행렬 접근법을 사용하여 조임점 전이 근처에서 모델링하기.
• 평균장 근사법을 통해 무질서한 고체에서 보존 피크와 저주파수 플랫폼의 기원을 이해하기.
• 구형 체적 효과와 상호작용을 모델링하기 위해 각도 차단 파라미터 $ K_d $ 를 도입하여 무작위 탄성 네트워크 모델을 개선하기.
• 이동된 Kesten-McKay DOS가 하이퍼스태틱 및 이소스태틱 영역에서 DOS를 타당한 평균장 기술로 제공하는지 테스트하기.
• 모든 $ \zeta \to 0 $ 조건에서 DOS가 이동된 Kesten-McKay 형태로 균일 수렴하는 조건을 규명하기.

제안 방법

• 평균 연결수 $ \zeta $ 를 가진 에르되시-레니 그래프를 $ z_0 $-정규 그래프에 중첩시켜 평균 연결수 $ Z = z_0 + \zeta $ 를 가진 소월드 네트워크를 구성한다.
• 구형 체적 효과를 표현하기 위해 접촉 단위 벡터 간의 각도 상관관계를 모델링하기 위해 $ K_d \leq 0.5 $ 의 차단 파라미터를 도입한다.
• 시스템의 헤시안 행렬을 사용하여 랜덤 블록 행렬의 스펙트럼 분석을 통해 진동 DOS를 계산한다.
• $ 2 \leq d \leq 4 $ 에서 캐비티 방법과 정확한 고유값 분해를 적용하여 DOS를 시뮬레이션하고, 분석적 형태와 비교한다.
• 모의 DOS를 이동된 Kesten-McKay(SKM) 및 Marchenko-Pastur(MP) 분포와 비교하며, 특히 작은 $ \omega $ 에서의 행동을 분석한다.
• 저주파수 행동과 보존 피크 위치를 분석하기 위해 효과적 매질 근사(Efficient Medium Approximation, EMA)와 스펙트럼 분포 모멘트를 활용한다.

실험 결과

연구 질문

• RQ1 소규모 세계 모델 $ z_0=2 $ 의 누적 DOS 는 $ \zeta \to 0 $ 일 때 이동된 Kesten-McKay DOS 로 균일하게 수렴하는가?
• RQ2 이동된 Kesten-McKay DOS 가 파라미터 $ t = Z/d $ 를 사용하여 조임점 근처 소프트 스피어 시스템의 진동 DOS 를 정확히 기술할 수 있는가, 특히 $ d=3 $ 에서는?
• RQ3 각도 차단 $ K_d $ 를 도입할 경우 저주파수 DOS 와 플랫폼 또는 갭 존재에 어떤 영향을 미치는가?
• RQ4 분자 동역학 시뮬레이션과의 정량적 일치도는 $ d=3 $ 에서 $ Z=7,8,9 $ 인 경우에 대해 어떻게 되는가?
• RQ5 모델은 $ \omega=0 $ 에서 이소스태틱 플랫폼을 재현하는가? 그리고 시뮬레이션 데이터와 가장 잘 일치하는 $ K_3 $ 값은 무엇인가?

주요 결과

• $ d=1 $ 에서는 파라미터 $ Z $ 를 가진 이동된 Kesten-McKay DOS 가 평균장 해로서 유효하며, $ \zeta \to 0 $ 일 때 구간 $[0, \omega_0]$ 에서 누적 DOS 가 이에 균일하게 수렴한다. 여기서 $ \omega_0 < \sqrt{z_0 - 1} + 1 $ 이다.
• Newman-Watts 소월드 모델인 $ z_0=2 $ 모델에서 수치 시뮬레이션은 $ \zeta \to 0 $ 극한에서 DOS 가 이동된 Kesten-McKay(SKM) DOS 로 균일하게 수렴한다는 추측을 지지한다.
• $ 2 \leq d \leq 4 $ 에서 큰 $ K_d $ 를 사용할 경우, 작은 $ \omega $ 에서 DOS 는 파라미터 $ t = Z/d $ 를 가진 SKM DOS 와 매우 유사하며, 특히 이소스태틱 경우에 잘 맞는다.
• $ d=3 $ 에서 SKM DOS 가 예측한 보존 피크 주파수는 분자 동역학 시뮬레이션과 잘 일치하며, $ Z=7 $ 과 $ 8 $ 에서는 오차가 4% 이내이다.
• $ \omega=0 $ 에서 이소스태틱 플랫폼은 $ K_3 \approx 0.3 $ 로 잘 재현되며, 분자 동역학 시뮬레이션에서 관측된 높이와 일치한다.
• 큰 $ K_d $ 를 사용할 경우 $ \omega=0 $ 근처에 갭이 나타나지만, SKM DOS 의 경우보다 작다. 따라서 $ D(\omega) \sim \omega^4 $ 는 성립하지 않으며, 이는 $ K_d=0 $ (무작위 탄성 네트워크) 경우와 다름을 보여준다.