다중입자 크기분포가 전이와 진동모드에 미치는 영향

초록

2차원에서 파워‑법 분포를 갖는 다중입자(λ=2~20) 연질 입자들을 MD로 압축한 결과, 폴리디스퍼시티는 접촉력과 배위수 분포, 그리고 임계 포장밀도 φc를 크게 변화시키지만, 압력·탄성계수·진동밀도(V​DOS)와 같은 거시적 기계적·진동적 특성은 Δz(초과배위수)만을 매개로 하는 기존의 임계 스케일링에 그대로 따른다.

상세 분석

본 연구는 2차원에서 파워‑법(지수 ν=3)으로 샘플링된 반경 분포를 갖는 연질 입자들을 N=2048개의 입자로 구성하고, 크기비 λ=Rmax/Rmin을 2에서 20까지 변화시켜 전이와 진동 특성을 정량적으로 조사하였다. 입자 간 상호작용은 선형 스프링 f = k δ n으로 모델링했으며, FIRE 알고리즘을 이용해 에너지 최소화를 수행해 정적 포장 상태를 얻었다.

첫 번째 주요 결과는 폴리디스퍼시티가 접촉력 분포 P(f)와 배위수 분포 P(z)를 현저히 넓힌다는 점이다. λ가 작을 때 P(f)는 평균 h_fi에 대해 가우시안 형태를 보이지만, λ가 10 이상으로 커지면 꼬리가 지수적으로 감쇠하는 형태로 전환된다. 이는 기존에 응력 이방성, 차원 증가, 비구형 입자 등으로 보고된 전이 현상과 동일한 메커니즘으로 해석될 수 있다. 배위수 분포는 λ가 커짐에 따라 파워‑법 P(z) ∝ z⁻⁴·²의 꼬리를 형성하고, 절단점 z*는 λ⁰·⁷⁴에 비례한다. 이는 작은 입자가 큰 입자 사이의 빈 공간을 메우면서 지역적인 과밀을 유발함을 의미한다.

두 번째로, 임계 포장밀도 φ_c는 λ에 따라 상승한다. φ_c − φ*_c ∝ (λ − 1)^0·32 로 피팅되었으며, φ*_c≈0.81은 단일입자(λ=1) 경우와 일치한다. 이는 폴리디스퍼시티가 작은 입자를 통해 공간 효율을 높여 전체적인 체적 점유율을 증가시키는 효과를 반영한다.

반면, 평균 배위수 ⟨z⟩, 압력 p, 초과배위수 Δz와 같은 거시적 양은 λ에 무관하게 기존의 임계 스케일링을 따른다. 구체적으로 p/k ∝ (φ − φ_c) ∝ Δz², Δz ∝ (φ − φ_c)¹ᐟ² 로서, λ가 변해도 기울기가 동일하게 유지된다. 탄성계수 G와 B 역시 Δz에 대한 함수 형태가 변하지 않는다. G/k ∝ Δz, B/k → 상수, 그리고 G/B ∝ Δz는 폴리디스퍼시티와 무관하게 동일한 임계 거동을 보인다. 이는 평균 배위수만이 선형 탄성 거동을 결정한다는 점을 강조한다.

진동 특성 측면에서는 동적 행렬을 대각화해 VDOS D(ω)를 구했으며, ω* ∝ Δz인 특성 주파수가 존재한다는 점에서 단일·이중입자 시스템과 동일한 플래토 형태의 VDOS를 확인했다. λ를 변화시켜도 VDOS 곡선 자체는 거의 겹쳐지며, 중간 주파수 영역에서의 참여비율 Pr(ω)는 λ가 증가할수록 약간 감소하지만 전체적인 모드 구조는 변하지 않는다. 따라서 폴리디스퍼시티는 미시적 분포를 바꾸지만, 거리‑전이(Δz)만이 기계적·진동적 거동을 지배한다는 결론에 도달한다.

이러한 결과는 폴리디스퍼시티가 큰 공학·지질 시스템(예: 지진단층의 파워‑법 입자 분포, 아폴로니안 포장)에서도 전이점 근처의 거시적 물성은 입자 크기 분포와 무관하게 거리‑전이 변수에 의해 결정된다는 중요한 시사점을 제공한다.