Preparation history can outweigh descriptor definition in frictionless bidisperse packings below jamming

기술 요약: 마찰 없는 이분산 충전물(Bidisperse Packings)의 준비 이력

문제 제기
알갱이 충전물의 구조적 및 구성적 기술은 크기 비(size ratio), 미세 입자 분율(fine fraction), 밀집도(crowding)와 같은 거시적 상태 변수를 고정하면 충전물의 형태와 그로부터 유도된 모든 기술자(descriptor)가 유일하게 결정된다는 가정에 의존한다. 이러한 가정은 준비 이력이 중요하다고 알려진 잼 상태(jammed states)에서는 잘 확립되어 있으나, 잼 점(jamming point)보다 훨씬 낮은 수준의 비잼(unjammed), 마찰 없는 이분산 충전물에 대해서는 검증되지 않은 상태이다. 구체적으로, 조립 경로(준비 프로토콜)가 기술자 자체의 임의적인 파라미터(예: 연결성을 정의하는 간극 임계값)를 변화시키는 것보다 접촉 기반 기술자를 더 크게 변화시키는지 여부는 불분명하다. 본 연구는 하부 잼 영역의 미세 네트워크 연결성 맥락에서 "준비 이력이 기술자 정의보다 더 큰 영향을 미칠 수 있는지"를 조사한다.

방법론
저자들은 주기적 입방 셀 내의 마찰 없는 이분산 구형 충전물에 대한 수치 시뮬레이션을 사용하였다. 연구는 모든 상태 변수—크기 비($\lambda$), 미세 부피 분율($\xi$), 고체 부피 분율($\phi$), 밀집도($s = \phi/\phi_m$)—를 제어하면서 오직 조립 프로토콜만을 변화시켰다. 두 가지 서로 다른 준비 패밀리가 비교되었다:

1. 성장 준비형 (G): 구들을 무작위로 균일하게 배치하여 명목 크기로 성장시킨 후, 중첩을 제거하기 위해 집단 재배열 기법을 적용하여 생성된 충전물.
2. 흔들림 평형형 (S): G 상태로부터 "흔들림-이완(shake–relax)" 사이클을 적용하여 유도된 충전물. 이 사이클에서 입자 위치는 가우시안 섭동(요동을 시뮬레이션함)을 받은 후 이완되며, 이때 고체 부피 분율은 엄격히 보존된다. 이 과정은 시스템을 압축하지 않으면서 성장 프로토콜에 의해 남겨진 근접 쌍(near-contact pairs)을 배출한다.

주요 분석 대상인 미세 네트워크 연결성은 연속체 퍼콜레이션(continuum percolation)을 통해 정의되었다: 두 미세 구는 표면 간극이 임계값 $t$(미세 직경의 비율) 미만일 때 연결된 것으로 간да. 연결성 지표는 최대 클러스터 분율($P_\infty^f$)이다. 본 연구는 다음을 체계적으로 비교하였다:

• 임계값 민감도 ($S_t$): 간극 임계값 $t$를 10배 범위(2% ~ 20%)로 변화시킴에 따라 발생하는 퍼콜레이션 개시점($\xi_g$)의 이동.
• 준비 민감도 ($S_P$): 고정된 참조 임계값($t=5\%$)에서 G 및 S 패밀리 간의 $\xi_g$ 이동.
• 이력 좌표 ($\chi_h$): G와 S 참조 상태 간의 근접 쌍(간극 < 미세 직경의 2%)의 접촉 분율을 정규화한 연속 변수로, 궤적 내의 연결성을 예측하는 데 사용됨.

실험은 크기 비 $\lambda = 4, 6, 8$ 및 밀집도 $s = 0.60, 0.65, 0.70$에서 수행되었다.

주요 결과

1. 준비가 기술자 정의를 압도함 (성장 준비형): 테스트된 모든 크기 비에서, 조립 경로를 변경함으로써 발생하는 퍼콜레이션 개시점의 이동($S_P$)은 성장 준비형 패밀리에 대해 10배의 임계값 변화에 의한 이동($S_t$)보다 컸다.

   • $\lambda = 4, s=0.70$에서 준비에 의한 이동은 $\Delta \xi_g = 0.334$였던 반면, 10배의 임계값 변화는 성장 준비형의 개시점을 $0.086$만큼만 이동시켰다.
   • 성장 준비형 분모를 기준으로 할 때, $\lambda = 4$에서 임계값 1 e-fold당 준비 민감도와 임계값 민감도의 비율($R = S_P/S_t$)은 약 8.9였다.
   • $\lambda = 6$ 및 $\lambda = 8$에서도 $R$은 1보다 유의미하게 큰 값(각각 3.3과 9의 하한값)을 유지하였으며, 이는 준비 이력이 성장 준비형 임계값 민감도에 비해 지배적인 요인임을 나타낸다.
   • 미묘한 차이: 그러나 $\lambda = 4$에서 특정 5–10% 임계값 대역에 대해 흔들림 평형형 패밀리의 임계값 민감도와 비교했을 때, 비율 $R$은 1 미만($R = 0.83$ per e-fold)으로 떨어졌다. 이는 임계값 민감도 자체가 준비 경로에 따라 달라질 수 있음을 의미하며, 모든 준비 패밀리에서 각 민감도를 확인해야 함을 시사한다.

2. 구조적 차이의 국소화: G와 S 충전물 사이의 구조적 차이는 매우 국소적이다. 가장 큰 불일치는 미세 직경의 2%보다 작은 간극을 가진 쌍의 집단에서 나타난다. 성장 프로토콜은 이러한 쌍들을 상당 부분 접촉(또는 근접 접촉) 상태로 남겨두는 반면, 흔들림-이완 과정은 이 집단을 배출한다. 이 좁은 구간($> 2\%$ 간극)을 넘어서면 두 패밀리의 쌍 분포는 수렴한다.

3. 예측 도구로서의 이력 좌표: 범주적 프로토콜 레이블(G vs. S)은 연속적인 이력 좌표 $\chi_h$로 대체될 수 있다. 고정된 조건(고정된 $\lambda, \xi, s$) 내에서 $\chi_h$는 최대 클러스터 분율($P_\infty^f$)을 강력하게 예측하며, 상태 전용 베이스라인에 비해 표본 외 예측 오차를 78% 감소시킨다. 그러나 이 관계는 서로 다른 밀집도 수준에 대해 보편적이지 않다. $s=0.70$에서 피팅된 모델은 $s=0.60$에서의 거동을 예측하는 데 실패한다.

4. 기계적 결과: 매칭된 조건($\lambda=4, \xi=0.35, s=0.70$)에서, 성장 준비 상태는 동일한 상태 변수와 입자 수를 가졌음에도 불구하고 흔들림 평형 상태보다 전단 강성이 약 6.83% 더 높게 나타났다. 이는 이력 의존적인 형태가 유효 기계적 특성에 직접적인 영향을 미친다는 것을 확인시켜 준다.

의의 및 주장
본 논문은 비잼 상태의 접촉 기반 기술자의 경우, 조립 경로가 기술자의 내부 파라미터(간극 임계값)의 임의적인 선택보다 기술자 값의 더 중요한 결정 요인이라고 주장한다(성장 준비형 임계값 민감도와 비교했을 때).

• 상태 변수 모델에 대한 도전: 저자들은 형태 정보가 포함된 유효 특성 폐쇄식(예: 강성을 상태 변수와 기술자 $M$의 관계로 나타내는 식)이 조성과 밀집도만으로는 유일하게 규정될 수 없다고 주장한다. 기술자 $M$ 자체도 이력에 의존하므로, 구성 모델에는 이력 좌표(토질 역학의 입자 배열 텐서와 유사한 개념)라는 네 번째 인자가 필요하다.
• 방법론적 테스트: 본 연구는 접촉 기술자가 상태 변수 모델에 통합되기 전에, 온셋 이동(준비 vs. 임계값)을 직접 비교하는 것이 필수적인 테스트라고 제안한다. 결정적으로, 임계값 민감도는 준비 경로에 따라 달라지므로(앞서 $\lambda=4$의 흔들림 평형 패밀리 사례에서 확인됨), 모든 준비 패밀리에 대해 두 민감도를 모두 확인해야 한다.
• 범위: 본 연구 결과는 잼 한계 근처의 마찰 없는 정적인 충전물에 특화되어 있다. 저자들은 근접 접촉 쌍의 상관관계가 잼 한계에서 보편적이라는 점은 알고 있으나, 여기서 관찰된 진폭의 차이(연결성 변화를 유발하는)는 그러한 점근적 보편성이 기대되지 않는 밀집도($s \le 0.70$)에서 발생한다고 언급하였다.

결론적으로, 접촉 기술자를 사용하여 알갱이 시스템의 상태를 정의하기 전에는, 두 가지 민감도(준비 및 임계값) 모두를 모든 준비 패밀리에 대해 확인해야 한다.