브라운 운동을 포함한 별형 폴리머 네트워크 파단의 사이클 랭크 마스터 커브

초록

본 연구는 별형 프리핸 네트워크에 브라운 운동을 도입하고 다양한 스트레인 레이트에서 단일축 신장을 수행하였다. 사이클 랭크를 기준으로 파단 시 변형 및 응력은 마스터 커브를 형성했으며, 이 커브는 스트레인 레이트에 따라 이동한다. 스트레인 레이트가 낮아질수록 에너지 최소화 시 얻은 결과에 수렴한다는 점을 확인하였다.

상세 분석

본 논문은 기존의 에너지 최소화 기반 네트워크 파단 연구가 열적 요동을 무시한다는 한계를 극복하고자, 별형 프리핸(phantom) 사슬 네트워크에 브라운 동역학을 적용하였다. 입자 위치는 Langevin 방정식으로 기술되며, 비선형 FENE 형태의 스프링으로 결합을 모델링한다. 반응 거리 C 이내에서 두 말단이 만나면 확률 D 로 새로운 스프링이 형성되어 네트워크가 구축되고, 전환율 E/ 로 전환 정도를 제어한다. 구축된 네트워크는 축방향 스트레인 레이트 F* 로 변형되며, 결합 길이가 임계값 9/ 를 초과하면 파단된다.

시뮬레이션은 스트레인 레이트를 10⁻³, 10⁻⁴, 10⁻⁵ (단위: τ⁻¹) 등 3가지 수준으로 변화시켜 수행하였다. 주요 관측값은 진응력 σ와 진변형 ε의 피크값을 각각 Y₄, F₄ 로 정의하고, 파단 전까지 끊어진 사슬 비율 E₄₅ 를 추적하였다. 결과는 다음과 같다. (1) 브라운 운동이 존재함에도 불구하고, 다양한 기능도 f와 전환율 E/ 에 대해 Y₄와 F₄는 사이클 랭크 ^(=N_c) 를 축으로 하는 마스터 커브를 형성한다. (2) 이 마스터 커브는 스트레인 레이트에 따라 좌우되며, 레이트가 감소할수록 에너지 최소화 시 얻은 마스터 커브에 근접한다. (3) 응력 피크 이후에도 완전한 응력 소멸이 일어나지 않으며, 이는 남아 있는 클러스터가 드래그 힘을 제공하기 때문이다. (4) 파단 전후의 구조적 이질성은 에너지 최소화 경우와 달리 크게 나타나지 않으며, 이는 열적 진동이 파단 과정을 연속적으로 퍼뜨리기 때문이다.

또한, 파워‑law 형태의 스케일링 지수(σ₄∝^‑0.16, ε₄∝^0.42)는 스트레인 레이트와 무관하게 일정했지만, 이전 에너지 최소화 연구에서 보고된 –½, 1/3 와는 차이를 보였다. 이는 본 연구에서 열적 분해율 9/ 를 보다 보수적으로 설정했기 때문으로 해석된다. 시스템 크기 의존성도 확인했으며, 큰 시스템에서는 에너지 최소화 결과에 수렴하기 위해 더 낮은 스트레인 레이트가 필요함을 제시한다.

결론적으로, 사이클 랭크는 열적 요동과 동역학적 파라미터가 포함된 실제 네트워크에서도 파단 특성을 통합적으로 설명하는 핵심 구조적 지표임을 재확인하였다. 다만, 열적 결합 파단, 배제 부피, 얽힘, 사슬 강성 등 실험적 현실을 반영하기 위한 추가 모델링이 필요하다.