엔탱글드 폴리머 크리프 거동의 원시 사슬 네트워크 시뮬레이션

초록

본 연구는 다중 사슬 슬립링크 모델을 기반으로 한 원시 사슬 네트워크 시뮬레이션을 변형하여 응력 제어 크리프 실험을 모사한다. 시뮬레이션 결과를 엔탱글드 폴리부타디엔 용액의 실험 데이터와 비교했을 때, 비선형 영역(큰 응력)에서 정성적으로 일치함을 확인하였다. 전이 과정에서 사슬의 스트레치와 방향 텐서, 그리고 슬립링크(엔탱글먼트) 밀도의 변화를 추적함으로써, 응력 제어 모드에서는 전단 속도 제어에 비해 분자 회전의 일관성이 방해받아 변형 유도 탈엔탱글먼트가 상대적으로 억제된다는 메커니즘을 제시한다.

상세 분석

이 논문은 기존의 원시 사슬 네트워크(Primitive Chain Network, PCN) 모델을 크리프 실험에 적용하기 위해 두 가지 주요 변형을 도입한다. 첫째, 슬립링크 위치와 사슬 길이의 라그랑지안 운동 방정식을 응력 제어 조건에 맞게 피드백 제어 루프와 결합함으로써, 외부 전단 응력이 목표값에 도달하도록 시뮬레이션 배경 흐름의 전단 속도를 자동 조정한다. 둘째, 슬립링크 생성·소멸 메커니즘을 구현하여 사슬 말단이 다른 사슬에 접촉할 때 새로운 엔탱글먼트를 형성하고, 말단이 슬립링크를 떠날 때 기존 엔탱글먼트를 소멸시키는 동적 토폴로지 변화를 허용한다. 이러한 구현은 기존의 전단 속도 제어 시뮬레이션보다 복잡하지만, 실제 실험에서 흔히 사용되는 응력 제어 크리프 조건을 직접 모사할 수 있게 한다.

시뮬레이션 파라미터는 실험 데이터와의 정량적 매핑을 위해 평균 사슬 길이, 열에너지, 사슬 이완 시간 등을 비차원화하였다. 특히, 실험에 사용된 폴리부타디엔 용액의 플래토 모듈러스와 엔탱글먼트 분자량을 기반으로 평균 엔탱글먼트 수와 사슬당 평균 쿠흔 세그먼트 수를 추정하였다. 이 과정에서 얻어진 비차원화된 탄성계수 W와 시간 스케일 τ는 선형 점탄성 영역에서 실험과 좋은 일치를 보였으나, 고주파 영역에서는 라우스 모드와 유리 모드가 모델에 포함되지 않아 차이가 발생한다.

비선형 크리프 거동을 분석하기 위해 전단 응력과 전단 속도의 시간 전개를 비교하였다. 전단 속도 제어 시뮬레이션은 실험과 유사한 점탄성 오버슈트와 전단 속도 언더슈트를 재현했으며, 이는 사슬의 방향 텐서(p*)가 급격히 증가하고 이후 사슬 스트레치(m*)가 뒤따라오는 전형적인 “스트레치‑오리엔테이션 디커플링” 메커니즘과 일치한다. 반면, 응력 제어 시뮬레이션에서는 전단 속도 언더슈트가 존재하지만, 전단 속도 자체는 실험보다 빠르게 정상 상태에 도달한다. 이는 응력 일정 조건 하에서 사슬 회전(코히런트 텀블링)이 방해받아 방향 텐서가 급격히 변하지 않으며, 결과적으로 엔탱글먼트 감소 속도가 전단 속도 제어보다 느려지는 현상으로 해석된다.

슬립링크 생존 비율(s₁)의 시간 전개를 통해 변형 유도 탈엔탱글먼트를 정량화하였다. 높은 전단 응력 혹은 전단 속도에서는 s₁이 급격히 감소하지만, 동일한 정상 상태 전단 응력에서도 전단 속도 제어 경우가 응력 제어 경우보다 s₁ 감소가 더 빠르다. 이는 코히런트 텀블링이 활발히 일어날 때 사슬이 서로 교차하고 슬립링크가 파괴되는 확률이 높아지기 때문이다. 따라서, 응력 제어 모드에서는 텀블링의 일관성이 깨져 탈엔탱글먼트가 억제된다는 결론을 도출한다.

전체적으로 이 연구는 원시 사슬 네트워크 모델을 응력 제어 크리프 실험에 성공적으로 적용함으로써, 비선형 크리프 영역에서의 미시적 메커니즘—특히 엔탱글먼트 밀도 변화와 사슬 회전의 상관관계—을 밝히는 데 기여한다. 또한, 전단 속도와 전단 응력이라는 두 제어 방식이 동일한 정상 상태 전단 응력을 달성하더라도 내부 토폴로지와 동역학에 서로 다른 영향을 미친다는 중요한 물리적 통찰을 제공한다.