Energy Barriers for Reversible Chain Scission and Healing under Tension with Displacement Control

기술 요약: 변위 제어 하의 인장 상태에서 가역적 사슬 절단 및 회복을 위한 에너지 장벽

문제 정의
고분자 사슬의 절단(scission)은 엘라스토머나 하이드로젤과 같은 연성 재료의 파괴를 지배하는 근본적인 메커니즘이다. 단일 분자 힘 분광법 실험은 사슬 절단이 확률적이며 로딩 속도와 환경 조건에 의존한다는 것을 입증해 왔으나, 전통적인 이론적 기술들은 대개 힘 제어(force-controlled) 조건만을 가정한다. 일정한 힘이 가해지는 상황에서 고분자 사슬의 자유 에너지 지형은 일반적으로 하나의 국소 최솟값(온전한 상태)과 하나의 전이 상태만을 가지며, 이는 사슬 절단을 비가역적으로 만들어 회복을 불가능하게 한다. 그러나 연성 재료에서 관찰되는 자가 치유 현상은 사슬 절단이 가역적일 수 있음을 시사한다. 본 논문은 힘 대신 고분자 사슬의 종단 간 거리(end-to-end distance)가 규정되는 변위 제어(displacement-controlled) 조건 하에서 가역적인 사슬 절단과 회복이 가능한지 조사함으로써 이러한 이론적 공백을 다룬다.

방법론
저자들은 인장 하의 단일 고분자 사슬의 열역학 및 역학을 분석하기 위해 "절단 가능한 자유 결합 사슬(breakable freely-jointed chain, FJC)" 모델을 개발하였다.

• 모델 정식화: 사슬은 $n+1$개의 동일한 링크(쿠른 세그먼트)로 구성된다. 모델은 하나의 임의의 링크가 해리되는 동안 나머지 $n$개의 링크는 강직한 세그먼트로 행동한다고 가정한다. 절단 가능한 링크에는 결합 신장 및 해리 에너지를 고려하기 위해 레너드-존스(Lennard-Jones, LJ) 퍼텐셜이 할당된다.
• 열역학적 프레임워크: 시스템은 온도 $T$인 열적 저장소로 취급된다. 자유 에너지 $A(x; y)$는 내부 변수 $x$(절단 가능한 링크의 길이)와 제어 매개변수 $y$(규정된 종단 간 거리)의 함수로 유도된다. 이 유도 과정은 절단 가능한 링크의 퍼텐셜 에너지와 강직한 사슬 파편들의 엔트로피 기여를 통합한다.
• 통계적 근사: 강직한 파편들에 대해 두 가지 통계적 처리가 사용된다:
  1. 가우시안 근사(Gaussian Approximation): 작은 신장 범위에서 유효하며, 닫힌 형태의 자유 에너지 함수를 제공한다.
  2. 비가우시안(Kuhn-Grun) 근사: 유한한 신장성을 고려하기 위해 역 라플란지 함수(inverse Langevin function)를 활용하며, 경로 길이 근처에서 더 정확한 묘사를 제공한다.
• 역학적 분석: 유도된 자유 에너지 지형을 사용하여 저자들은 사슬 절단($E_s$)과 회복($E_h$)을 위한 에너지 장벽을 계산한다. 저자들은 사슬의 생존 확률을 시간에 따라 예측하기 위해 1차 속도 방정식(아레니우스 법칙)에 기반한 역학적 접근 방식을 채택한다. 또한, 사슬 절단과 회복의 확률적 특성을 시각화하고 분석적 속도 방정식을 검증하기 위해 키네틱 몬테카를로(Kinetic Monte Carlo, KMC) 시뮬레이션을 수행한다.

주요 기여 및 결과

1. 변위 제어 하의 가역성: 본 연구는 변위 제어 조건 하에서 자유 에너지 함수가 특정 종단 간 거리 범위 내에서 두 개의 국소 최솟값(온전한 상태와 해리된 상태)과 이를 구분하는 국소 최댓값(전이 상태)을 나타냄을 보여준다. 이러한 이중 최솟값 지형은 힘 제어 모델에서는 불가능했던 사슬 절단($E_s$)과 회방($E_h$)에 대한 별도의 에너지 장벽 정의를 가능하게 한다.
2. 비선형 에너지 장벽:
   • 절단 장벽 ($E_s$): 변위가 0일 때 결합 에너지 근처에서 시작하여, 변위가 증가함에 따라 약 이차 함수 형태로 감소하며, 임계 변위($y_c \approx n \cdot l$)에서 소멸한다.
   • 회복 장벽 ($E_h$): 낮은 수준($k_B T$와 유사)에서 시작하여 변위에 따라 약 이차 함수 형태로 증가한다. 결정적으로, $E_h$는 사슬 길이에 매우 의존적(엔트로피적 성질)인 반면 결합 에너지에는 거의 무관하지만, $E_s$는 링크의 결합 에너지에 의해 지배된다.
3. 속도 의존적 강도: 모델은 고분자 사슬의 파단력이 확률적이며 속도 의존적임을 예측한다.
   • 극도로 낮은 인장 속도에서 시스템은 평형에 도달하며, 평균 파단력은 평형 확률 분포에 의해 결정된다. 이 하한값은 이론적인 공유 결합 강도보다 2~3 자릿수(orders of magnitude) 더 낮게 나타난다.
   • 높은 인장 속도에서 사슬은 강직한 막대처럼 행동하며 파단력이 증가하지만, 분자 동역학 시뮬레이션에서 일반적인 속도에서도 에너지 장벽 없이 즉각적으로 결합을 끊는 데 필요한 이론적 상한값보다는 여전히 현저히 낮다.
4. 임계 변위: 모델은 온전한 상태의 최솟값이 사라져 즉각적인 사슬 절단으로 이어지는 임계 종단 간 거리를 식별한다. 이 임계 거리는 결합 에너지에는 민감하지 않으나 사슬 길이에 따라 스케일링된다.

의의 및 주장
본 논문은 기존 모델들이 힘 제어 조건을 가정함으로써 가역적인 사슬 절단과 회복의 가능성을 가려왔다고 주장한다. 변위 제어로 전환함으로써, 저자들은 변위가 임계치 미만일 때 회복이 에너지적으로 유리함을 보여줌으로써 고분자 네트워크의 자가 치유에 대한 열역학적 근거를 제공한다.

이러한 발견의 의의는 연성 재료에서 관찰되는 속도 의존적 파괴 및 자가 치유 행동을 설명할 수 있는 잠재력에 있다. 저자들은 이러한 통계적 및 역학적 메커니즘을 거시적인 고분자 네트워크 모델에 통합함으로써, 특히 가역적 절단과 회피가 역할을 하는 반복 하중 조건에서의 균열 저항 및 손상 진화에 대한 이해를 높일 수 있다고 제안한다. 본 연구는 고분자 네트워크의 거시적 강도가 단순히 결합 강도의 함수가 아니라, 인장 속도와 사슬 구조에 의해 강력하게 영향을 받는 사슬 절단의 통계적 역학과 회복에 대한 엔트로피적 장벽에 의해 결정된다는 점을 확립한다.