Cohesion-induced hysteresis and breakdown of marginal stability in jammed granular materials

이 논문은 이산 요소 시뮬레이션과 유효 매질 이론을 통해 점성 입자 응집체에서 인력 상호작용이 한계 안정성을 붕괴시켜 전단 탄성률의 히스테리시스를 유발하고 압력을 유일한 상태 변수로 만드는 기존 가정을 무효화한다는 것을 규명합니다.

핵심

이 논문은 **"찌그러진 모래알갱이와 끈적한 모래알갱이의 비밀"**에 대한 이야기입니다. 과학자들이 왜 어떤 모래는 눌렀다 놓으면 원래대로 돌아오지 않고, 또 어떤 모래는 눌렀을 때와 놓았을 때의 상태가 완전히 다르게 행동하는지 그 이유를 찾아냈습니다.

간단한 비유로 설명해 드릴게요.

1. 두 가지 모래 알갱이: "미운 친구" vs "친구"

우리가 상상하는 모래 알갱이는 보통 **미운 친구 (반발력)**입니다. 서로 닿으면 밀어내려고 하죠. 이 모래들을 통에 넣고 꽉꽉 누르면 (압축), 단단해집니다. 하지만 누르는 힘을 빼면 (이완), 다시 원래대로 돌아갑니다. 이 경우, "얼마나 눌렀느냐 (압력)"만 알면 그 상태가 딱 정해집니다.

하지만 현실의 모래, 특히 젖은 모래나 끈적이는 물질에는 **친구 (인력)**가 있습니다. 서로 붙어 있으려는 성질이죠. 이 논문은 바로 이 "끈적이는 모래"를 연구했습니다.

2. 놀라운 발견: "기억력"이 생긴 모래

과학자들은 끈적이는 모래를 실험실에서 다음과 같이 실험했습니다.

1. 압축: 모래를 점점 더 꽉 누릅니다.
2. 이완: 다시 힘을 빼면서 원래대로 풀어줍니다.

그런데 놀라운 일이 벌어졌습니다. 반발력만 있는 모래는 누를 때와 놓을 때의 상태가 압력만 같으면 똑같았습니다. 하지만 끈적이는 모래는 달랐습니다.

• 비유: 마치 진흙 공을 생각해보세요.
  • 반발력 모래는 공기 풍선 같습니다. 누르면 작아지고, 손을 떼면 원래 크기로 돌아갑니다.
  • 끈적이는 모래는 점토나 달걀 흰자 같습니다. 누르고 나면 모양이 변해버리고, 힘을 빼도 원래대로 돌아오지 않습니다.

가장 중요한 점은, 끈적이는 모래는 "압력"이라는 숫자만으로는 상태를 설명할 수 없다는 것입니다. 같은 압력이라도, "누가 먼저 눌렀는지 (압축 과정)"와 "누가 먼저 놓았는지 (이완 과정)"에 따라 단단함 (전단 탄성률) 이 완전히 다릅니다. 마치 기억력이 있는 것처럼 말이죠.

3. 왜 이런 일이 일어날까? "마진 (가장자리) 의 붕괴"

과학자들은 이 현상을 설명하기 위해 **'유효 매질 이론 (EMT)'**이라는 수학적 도구를 사용했습니다. 이를 쉽게 비유하자면, 복잡한 모래 알갱이들의 관계를 **'스프링으로 연결된 거미줄'**로 상상하는 것입니다.

• 반발력 모래 (마진 안정성):
  이 모래들은 서로 밀어내면서 딱딱하게 서 있습니다. 마치 마진 (가장자리) 에 서 있는 줄다리기처럼, 한쪽이 조금만 힘을 빼도 무너질 것 같은 불안정한 상태입니다. 이 상태에서는 압력과 연결된 스프링의 개수 (접촉 수) 가 딱딱하게 연결되어 있어서, 압력만 알면 단단함도 정확히 예측할 수 있습니다.

• 끈적이는 모래 (마진 붕괴):
  끈적이는 모래는 서로 붙어 있으려고 (인력) 합니다. 이 때문에 "마진 (가장자리) 에 서 있는 줄다리기" 상태가 깨집니다. 서로 붙어있기 때문에, 압력이 0 이 되어도 (힘을 완전히 뺐을 때) 서로 떨어지지 않고 단단하게 붙어 있는 상태를 유지합니다.

  핵심: 끈적이는 모래는 압력이 0 이 되어도 단단할 수 있습니다. 반발력 모래는 압력이 0 이 되면 흐물흐물해지지만, 끈적이는 모래는 서로 붙어있기 때문에 여전히 단단합니다. 이것이 바로 "마진 안정성의 붕괴"입니다.

4. 결론: 왜 이 연구가 중요할까?

이 연구는 **"끈적이는 힘 (인력) 이 물질의 성질을 어떻게 바꾸는지"**를 수학적으로 증명했습니다.

• 기존 생각: "압력만 알면 그 물질의 단단함을 알 수 있다."
• 새로운 발견: "끈적이는 물질은 **어떻게 눌렀는지 (역사)**를 기억합니다. 같은 압력이라도, 눌렀다가 놓은 상태와 처음부터 놓아둔 상태는 완전히 다릅니다."

일상생활에서의 예시:

• 건조한 모래 (반발력): 모래성 쌓을 때, 모래를 꽉 누르면 단단해지지만, 손을 떼면 흐트러집니다.
• 젖은 모래 (끈적임): 젖은 모래로 성을 쌓으면, 손을 떼어도 성이 그대로 유지됩니다. 심지어 너무 많이 물을 부어 (압력을 줄여) 성이 무너져도, 모래 알갱이들이 서로 붙어있어 완전히 흐물흐물해지지 않고 덩어리를 유지합니다.

이 논문은 바로 그 젖은 모래 (끈적이는 입자) 가 왜 이렇게 '기억력'을 가지고 행동하는지에 대한 물리학적 원리를 찾아낸 것입니다. 이는 지진, 눈사태, 혹은 산업용 분말 처리 등 다양한 분야에서 끈적이는 입자 물질의 거동을 예측하는 데 큰 도움이 될 것입니다.

기술 요약

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

• 배경: 무질서한 고체 (입자 집합체, 콜로이드, 거품 등) 는 임계 부피 분율 ($\phi_J$) 을 넘으면 '잠김 (Jamming)' 전이를 겪어 고체처럼 거동합니다. 잠김 근처에서 기계적 특성 (전단 탄성률 $G$ 등) 은 부피 분율과의 거리 ($\phi - \phi_J$) 가 아닌, 압력 ($p$) 이나 배위수 ($Z$) 와 같은 상태 변수로 설명될 때 보편적인 스케일링 법칙을 따릅니다.
• 기존 지식: 순수하게 반발력 (repulsive) 만 작용하는 입자의 경우, 잠김 상태는 '마진 안정성 (Marginal Stability)'을 가지며, 이는 압력 $p$와 배위수 $Z$가 고유한 관계를 갖는다는 것을 의미합니다. 따라서 $p$나 $Z$로 특성을 표현하면 준비 과정 (압축/이완) 에 따른 히스테리시스가 사라집니다.
• 문제점: 실제 입자계 (습한 모래 등) 에는 모세관력이나 젖음 현상과 같은 인력 (attractive force, 응집력) 이 존재합니다. 그러나 응집력이 있는 시스템에서 동일한 이론적 설명이 적용되는지, 그리고 압력이나 배위수로 특성을 표현할 때 히스테리시스가 사라지는지에 대한 명확한 이론적 틀과 체계적인 분석이 부족했습니다.
• 핵심 질문: 응집력이 있는 입자계에서도 압력 ($p$) 이 유일한 상태 변수로 작용하여 히스테리시스를 제거할 수 있는가?

2. 연구 방법론 (Methodology)

• 시뮬레이션 (DEM):
  • 마찰이 없는 단분산 입자 시스템을 3 차원 주기적 경계 조건 하에서 모델링했습니다.
  • 입자 간 힘은 반발 코어, 단거리 인력 꼬리, 컷오프 거리를 가진 조각별 선형 힘 법칙 (piecewise linear force law) 으로 정의했습니다.
  • 부피 분율 ($\phi$) 을 제어하여 압축 (compression) 과 이완 (decompression) 과정을 반복 수행하며, 기계적 평형 상태에서의 압력 ($p$), 전단 탄성률 ($G$), 배위수 ($Z$) 를 측정했습니다.
  • 인력 강도 ($\alpha$) 를 0 (반발력만) 과 0.001 (응집력 있음) 로 설정하여 비교했습니다.
• 이론적 접근 (유효 매질 이론, EMT):
  • 응집 상호작용을 포함하도록 유효 매질 이론 (Effective Medium Theory, EMT) 을 확장했습니다.
  • 시스템의 기계적 응답을 결정하는 헤시안 행렬 (Hessian matrix) 을 분석하여, 인력이 시스템의 안정성 (마진 안정성) 에 미치는 영향을 이론적으로 규명했습니다.

3. 주요 결과 및 발견 (Key Results)

가. 응집력에 의한 히스테리시스의 지속성

• 반발력 입자 ($\alpha=0$): 압축 및 이완 과정에서 $G$와 $p$의 관계를 보면, $p$를 상태 변수로 사용할 때 히스테리시스가 사라지고 $G \propto p^{1/2}$의 스케일링 법칙을 따릅니다.
• 응집력 입자 ($\alpha>0$):
  • 압력-부피 분율 관계: 이완 과정에서 음의 압력 ($p<0$) 상태에서도 입자 집합체가 유지되며, $p=0$이 되는 부피 분율 ($\phi_{p=0}$) 이 압축 시와 다릅니다.
  • 전단 탄성률 ($G$) 의 히스테리시스: $G$를 $p$의 함수로 그렸을 때, 반발력 입자와 달리 히스테리시스가 명확하게 관찰됩니다. 특히 이완 과정에서 $p \to 0$일 때에도 $G$가 유한한 값을 유지합니다 (음의 압력에서도 강성이 존재함).
  • 결론: 응집력이 있는 시스템에서는 압력 ($p$) 만으로는 기계적 상태를 완전히 설명할 수 없으며, 준비 과정의 이력 (history) 이 시스템의 강성에 영향을 미칩니다.

나. 마진 안정성 (Marginal Stability) 의 붕괴

• 헤시안 분석: 응집력 입자의 경우, 이론적으로 불안정한 접촉 (국소적으로 불안정한 영역, $s_{ij}=-1$) 이 존재할 수 있으나, 시뮬레이션 결과 기계적 평형 상태에서는 이러한 접촉이 제거되어 모든 접촉이 안정적 ($s_{ij}=1$) 인 것으로 나타났습니다.
• EMT 예측: 헤시안 구조가 반발력 입자와 동일하므로 EMT 공식은 동일하게 적용되지만, 마진 안정성 조건 ($p = p_c(Z)$) 이 깨집니다.
  • 반발력 입자: $p = p_c(Z)$ (압력과 배위수가 1:1 대응).
  • 응집력 입자: $p < p_c(Z)$ (압력이 임계값보다 낮음).
• 과잉 강성 (Excess Rigidity): 마진 안정성이 깨지면서 시스템은 '과잉 강성'을 가지게 됩니다. 즉, 배위수 $Z$가 같더라도 압력 $p$가 낮을 때 더 큰 전단 탄성률 $G$를 갖게 됩니다.

다. 일반화된 스케일링 법칙 검증

• EMT 를 통해 유도된 새로운 스케일링 법칙을 검증했습니다:
  $$\frac{G - G_0(Z)}{G_0(Z)} = \sqrt{\frac{p_c(Z) - p}{p_c(Z)}}$$
  • 여기서 $G_0(Z)$는 마진 안정성 하에서의 기본 전단 탄성률, $p_c(Z)$는 임계 압력입니다.
• 시뮬레이션 데이터는 이 식을 따르는 것으로 확인되었으며, 이는 응집력 입자에서 관찰된 히스테리시스가 마진 안정성의 붕괴에서 기인함을 정량적으로 증명했습니다.

4. 주요 기여 및 의의 (Significance)

1. 히스테리시스의 물리적 기원 규명: 미세한 입자 간 힘 법칙 자체에 히스테리시스가 없더라도, 응집력 (인력) 이 존재할 경우 기계적 접촉 네트워크의 진화 과정에서 거시적인 히스테리시스가 발생함을 밝혔습니다.
2. 상태 변수의 한계 제시: 반발력 입자계에서는 압력 ($p$) 이 유효한 상태 변수였으나, 응집력 입자계에서는 압력만으로는 상태를 결정할 수 없으며, 배위수 ($Z$) 가 압축 이력을 기억하는 추가적인 구조 변수로 작용함을 보였습니다.
3. 마진 안정성 붕괴의 보편성: 인력 상호작용이 잠김 (Jamming) 근처의 안정성 지형을 근본적으로 변화시켜 마진 안정성 조건을 위반하고 과잉 강성을 생성한다는 것을 이론 (EMT) 과 시뮬레이션으로 동시 입증했습니다.
4. 이론적 확장: 기존에 반발력 입자계에만 적용되던 유효 매질 이론 (EMT) 을 응집력 시스템으로 성공적으로 확장하여, 실제 습한 모래나 콜로이드 등 다양한 응집성 입자계의 기계적 거동을 설명하는 통일된 틀을 제시했습니다.

5. 결론

이 연구는 응집력이 있는 잠김 입자계에서 관찰되는 강한 기계적 히스테리시스가 미세한 힘 법칙의 비가역성이 아니라, 응집력에 의한 마진 안정성의 붕괴에서 비롯됨을 규명했습니다. 이는 압력이나 배위수와 같은 전통적인 상태 변수만으로는 응집성 물질의 기계적 특성을 완전히 설명할 수 없으며, 시스템의 이력 의존성을 이해하기 위해 새로운 이론적 접근이 필요함을 시사합니다.