Relaxation dynamics and long-time tails explain shear-induced diffusion of soft athermal particles near jamming

이 논문은 재밍 (jamming) 전이 근처에 있는 2 차원 연성 비열적 (athermal) 입자 시스템에서 전단 (shear) 에 의해 유도된 확산 (shear-induced diffusion) 의 거동을 수치 시뮬레이션을 통해 연구한 것입니다. 저자들은 확산 계수 $D$와 입자 속도의 자기 상관 함수 (velocity auto-correlation function) 및 이완 시간 (relaxation time) 사이의 관계를 그린 - 쿠보 (Green-Kubo, GK) 공식을 통해 규명하고, 재밍 전이 근처에서의 임계 스케일링 (critical scaling) 과 장시간 꼬리 (long-time tail) 현상을 설명했습니다.

다음은 논문의 상세한 기술적 요약입니다.

1. 연구 배경 및 문제 제기 (Problem)

• 배경: 거품, 유화액, 콜로이드 현탁액, 입자성 물질 등 연성 비열적 입자 시스템은 재밍 전이 (jamming transition) 에서 임계적인 거동을 보입니다. 기존 연구들은 전단 유도 확산 계수 $D$가 입자 직경 $d_0$, 전단율 $\dot{\gamma}$, 그리고 시스템의 특징 길이 척도 (collective motion size, $\xi$) 와 어떻게 관련되는지에 초점을 맞추어 왔습니다.
• 문제점: 열적 시스템 (고전적 유체, 유리 등) 의 수송 특성은 이완 역학 (relaxation dynamics) 으로 설명되지만, 비열적 입자 시스템의 전단 유도 확산에서 이완 역학의 역할은 잘 알려져 있지 않았습니다. 특히, 재밍 전이 근처에서 확산 계수가 시스템 크기에 따라 발산하는 현상과 속도 상관 함수의 장시간 꼬리 (long-time tail) 사이의 인과 관계를 명확히 규명할 필요가 있었습니다.

2. 방법론 (Methodology)

• 수치 시뮬레이션: 2 차원 연성 비열적 입자 시스템에 대한 분자 동역학 (MD) 시뮬레이션을 수행했습니다.
  • 시스템: 직경이 $d_0$와 $1.4d_0$인 두 종류의 입자를 1:1 비율로 혼합하여 총 $N=8192$개의 입자로 구성했습니다.
  • 상호작용: 접촉 시 선형 탄성력 ($f_{ij}^e = k\delta_{ij}n_{ij}$) 과 감쇠력 ($f_i^d = -\zeta\{v_i - u(r_i)\}$) 을 적용하여 과감쇠 (overdamped) 역학을 모델링했습니다.
  • 조건: Lees-Edwards 경계 조건을 사용하여 전단 흐름을 구현했습니다. 충진률 ($\phi$) 은 재밍 임계값 ($\phi_J \approx 0.8433$) 주변 ($0.8 \le \phi \le 0.9$) 에서, 전단율 ($\dot{\gamma}$) 은 $10^{-7} \le \dot{\gamma}t_0 \le 10^{-2}$ 범위에서 변화시켰습니다.
• 분석 도구:
  • 확산 계수 ($D$): 횡방향 평균 제곱 변위 (MSD) 를 통해 정의 ($D = \lim_{\tau \to \infty} \Delta y^2 / 2\tau$).
  • 그린 - 쿠보 (GK) 공식: $D = \langle v_y^2 \rangle \int_0^\infty C(\tau) d\tau$. 여기서 $\langle v_y^2 \rangle$는 평균 제곱 속도, $C(\tau)$는 정규화된 속도 자기 상관 함수입니다.
  • 스케일링 분석: 다양한 물리량 ($D$, $\langle v_y^2 \rangle$, 이완 시간 $\tau^*$) 에 대한 임계 지수 (critical exponents) 를 추정하고 데이터 붕괴 (data collapse) 를 확인했습니다.

3. 주요 기여 및 결과 (Key Contributions & Results)

A. 평균 제곱 속도와 이완 시간의 임계 스케일링

• 평균 제곱 속도 ($\langle v_y^2 \rangle$): 재밍 근처에서 $\langle v_y^2 \rangle$는 임계 스케일링을 따릅니다.
  • 뉴턴 유체 영역 ($\phi < \phi_J$): $\langle v_y^2 \rangle \sim \dot{\gamma}^2$ (점도 $\eta$가 발산하는 것과 일치).
  • 준정적 영역 ($\phi > \phi_J$): 항복 응력 ($\sigma_Y$) 으로 인해 $\langle v_y^2 \rangle \sim \dot{\gamma}$로 선형적으로 비례합니다.
  • 스케일링 함수: $\langle v_y^2 \rangle / |\Delta\phi|^\psi = F(\dot{\gamma}/|\Delta\phi|^\lambda)$ 형태의 스케일링 관계가 성립하며, 지수 $\psi \approx 5.39$, $\lambda \approx 3.87$로 추정되었습니다.
• 이완 시간 ($\tau^*$): 속도 상관 함수를 늘어난 지수 함수 (stretched exponential, $C(\tau) = \exp[-(\tau/\tau^*)^\alpha]$) 로 피팅하여 이완 시간을 추출했습니다.
  • 재밍 아래 ($\phi < \phi_J$) 또는 전단율이 충분히 큰 경우, $\tau^*$는 임계 스케일링을 따릅니다 ($\tau^* \sim |\Delta\phi|^{-\epsilon} \dot{\gamma}^{-1}$ 등).
  • 지수 $\epsilon \approx 2.57$로 추정되었습니다.

B. 장시간 꼬리 (Long-time Tail) 와 확산 계수의 발산

• 장시간 꼬리의 발견: 시스템이 재밍 위 ($\phi > \phi_J$) 이고 전단율이 매우 작은 준정적 (quasi-static) 영역 ($\dot{\gamma}t_0 \le 10^{-6}$) 에 있을 때, 속도 자기 상관 함수 $C(\tau)$는 늘어난 지수 함수로 설명되지 않습니다. 대신 멱함수 법칙 (power-law) 인 $C(\tau) \sim (\dot{\gamma}\tau)^{-1}$ 형태의 장시간 꼬리가 관찰됩니다.
• GK 공식의 발산: $C(\tau) \sim \tau^{-1}$일 때, GK 공식의 시간 적분 $\int_0^\infty C(\tau) d\tau$가 발산합니다. 이는 열역학적 극한 ($L \to \infty$) 에서 확산 계수 $D$가 발산함을 의미하며, 이는 기존 연구에서 보고된 유한 크기 스케일링 ($D/\dot{\gamma} \sim L^\beta, \beta \ge 1$) 과 일치합니다.
• 2 차원 vs 3 차원: 이는 고전적 유체에서 2 차원일 때 GK 관계가 장시간 꼬리에 의해 깨지지만 3 차원에서는 그렇지 않은 현상과 유사한 비평형 비열적 시스템의 특징을 보여줍니다.

C. 축소된 GK 공식과 임계 지수 스케일링 관계

• 축소된 GK 공식: $C(\tau)$가 늘어난 지수 함수로 잘 설명되는 영역에서, GK 공식은 $D = \frac{\tau^* \langle v_y^2 \rangle}{\alpha} \Gamma(1/\alpha)$로 단순화됩니다.
• 임계 지수 관계식: 이 축소된 공식에 각 물리량의 임계 스케일링을 대입하여 다음 관계를 유도했습니다:
  $$\nu + \psi - \epsilon = \lambda$$
  여기서 $\nu$는 확산 계수의 지수, $\psi$는 속도 지수, $\epsilon$은 이완 시간 지수, $\lambda$는 전단율 스케일링 지수입니다.
  • 계산된 값: $\nu + \psi - \epsilon \approx 3.76$, $\lambda \approx 3.87$로, 수치적으로 매우 잘 일치합니다. 이는 서로 다른 물리량 사이의 임계 지수가 보편성 (universality) 을 공유함을 시사합니다.

4. 의의 및 결론 (Significance)

• 이론적 통찰: 기존 연구들이 확산 계수를 공간적 길이 척도 ($\xi, L$) 와 연결했던 것과 달리, 본 연구는 시간 척도 (이완 시간) 와 장시간 꼬리를 통해 확산 계수의 거동과 발산 메커니즘을 설명했습니다.
• 재밍 전이 이해: 재밍 전이 근처에서 비열적 입자 시스템의 수송 특성이 임계 현상 (critical phenomena) 으로 잘 설명될 수 있음을 입증했습니다.
• 장시간 꼬리의 중요성: 준정적 전단 흐름에서 장시간 꼬리가 발생하여 확산 계수가 발산한다는 점은, 2 차원 비열적 시스템의 유한 크기 효과를 이해하는 데 핵심적인 요소임을 강조했습니다.
• 스케일링 관계의 제안: 다양한 물리량의 임계 지수를 연결하는 새로운 스케일링 관계식 ($\nu + \psi - \epsilon = \lambda$) 을 제안하여, 향후 이론적 배경을 규명하는 데 기여할 수 있는 기반을 마련했습니다.

요약하자면, 이 논문은 그린 - 쿠보 공식을 활용하여 재밍 근처의 비열적 입자 시스템에서 속도 상관 함수의 장시간 꼬리가 확산 계수의 발산을 유발한다는 메커니즘을 규명하고, 이를 통해 다양한 물리량의 임계 지수 간의 보편적인 스케일링 관계를 제시한 중요한 연구입니다.