Analysing gelation transition through fractional viscoelasticity and Mittag-Leffler-Prabhakar function

기술 요약: 분수 차수 점탄성 및 미텔-레플레르-프라바카르 함수를 통한 겔화 전이 분석

문제 정의
점탄성 액체가 점탄성 고체로(또는 그 반대로) 변하는 과정인 겔화 전이는 스케일 불변의 멱법칙(power-law) 점탄성을 특징으로 하는 임계 겔 상태에 의해 특징지어진다. 분수 미분(fractional calculus), 특히 스프링-팟(spring-pot) 요소는 이러한 멱법칙 거동을 설명하는 자연스러운 수학적 언어를 제공하지만, 기존 문헌에는 상당한 한계가 존재한다. 기존의 분수 점탄성 모델은 주로 임계 겔 지점에 엄격하게 국한되어 적용되거나, 겔 형성 전(sol) 및 겔 형성 후(gel) 영역에 독립적으로 적용되어 왔다. 결정적으로, 기존의 프레임워크는 시스템이 임계점을 통과할 때 동적 탄성률 및 그 1차 도함수가 가교 밀도($p$)에 대해 연속적임을 명시적으로 강제하지 않는다. 이러한 누락은 대칭 조건($\kappa_S = \kappa_G = \kappa$)과 하이퍼스케일링 관계($n = z/(z+s)$)를 경험적 관찰이 아닌 이론적 필연성으로서 엄밀하게 유도하는 것을 방해한다. 또한, 표준적인 분수 역학 조립체는 허용 가능한 전체 임계 지수 범위에서 물리적 제약 조건(예: 완화의 단조성)을 충족하지 못하는 경우가 많다.

방법론
저자들은 겔화 전이에 걸친 선형 점탄성 특성의 진화를 모델링하기 위해 분수 점탄성과 3-매개변수 미텔-레플레르-프라바카르(Mittag-Leffler-Prabhakar) 함수에 기반한 통합 이론적 프레임워크를 개발하였다. 방법론은 다음과 같다:

1. 제약 조건 정식화: 모델은 내재적 열역학적 제약 조건(인과율, 비음의 탄성률, 단조성)과 임계 겔 이론에서 발견되는 특정 스케일링 관계를 충족하도록 구성된다. 가장 중요한 것은, 이 프레임워크는 임계 겔 지점($p_c$)에서 복소 탄성률 $G^*(\omega)$와 그 $p$에 대한 1차 도함수의 연속성을 강제한다는 점이다.
2. Pre-Gel(Sol) 모델링:
   • 모델 1: 두 개의 스프링-팟과 하나의 대시팟(dashpot)이 직렬로 연결된 조립체이다. 이는 2-매개변수 미텔-레플레르 함수로 표현되는 완화 계수를 산출한다. 그러나 이 모델은 $n + \kappa < 1$인 시스템으로 제한되며, 임계점 근처($\varepsilon \to 0$)에서만 점근적으로 정확하다.
   • 모델 2: 3-매개변수 미텔-레플레르-프라바카르 함수를 사용한 일반화 모델이다. 이 모델은 $n + \kappa < 1$ 제한을 제거하여 모든 $n, \kappa \in (0,1)$에 대해 물리적으로 허용 가능한 기술을 가능하게 한다. 이 모델은 프라바카르 지수 $\gamma_S$가 점도 발산 지수 $s$와 같음을 유일하게 결정한다.
3. Post-Gel(Gel) 모델링:
   • 모델 1: 후크 법칙을 따르는 스프링(평형 탄성률 제공)과 스프링-팟 배열이 병렬로 연결된 조립체이다.
   • 모델 2: 프라바카르 함수의 절단된 급수 전개이다. 이는 특정 계수 치환 하에서 모델 1과 대수적으로 동등함을 보이지만, 더 제한적이고 물리적으로 구조화된 매개변수 집합을 제공한다.
4. 연속성 강제: Pre-gel 모델과 Post-gel 모델은 동시에 풀린다. 임계 겔 지점에서의 동적 탄성률의 도함수가 일치해야 한다는 요구 조건은 양측의 완화 스케일링 지수가 동일($\kappa_S = \kappa_G = \kappa$)하고, 주요 차수 이탈 계수가 같아야 함을 강제한다.
5. 실험적 검증: 제안된 모델은 다음 두 가지 뚜렷한 실험 데이터셋에 의해 검증된다:
   • 솔-겔 전이를 겪는 화학적 가교 폴리디메틸실록산(PDMS)의 시간 영역 완화 계수 데이터($n=0.45$).
   • 겔-솔 전이를 겪는 열가역적 폴리비닐 알코올(PVA) 용액의 주파수 영역 동적 탄성률 데이터($n=0.77$).

주요 기여 및 결과

• 통합 프레임워크: 본 논문은 연속성 제약을 통해 pre-gel 및 post-gel 영역을 엄밀하게 연결하여, 전이를 별개의 상태가 아닌 단일한 연속 과정으로 취급하는 최초의 분수 점탄성 프레임워크를 제시한다.
• 하이퍼스케일링의 이론적 필연성: 동적 탄성률의 1차 도함수 연속성을 강제함으로써, 저자들은 대칭 조건($\kappa_S = \kappa_G$)과 하이퍼스케일링 관계($n = z/(z+s)$)가 경험적 우연이 아니라 물리적 일관성을 위해 요구되는 이론적 필연성임을 입증한다.
• 프라바카르 함수의 우수성: 프라바카르 기반의 Pre-Gel 모델 2는 표준적인 스프링-팟/대시팟 조립체보다 우수함이 입증되었다. 이 모델은 점도 발산 지수 $s$를 프라바카르 지수 $\gamma_S = s$를 통해 내재적으로 포착하며, 점근적 극한에서 사라지는 대시팟 요소를 필요로 하지 않는다. 이는 프라바카르 함수가 기계적 아날로그보다 pre-gel 완화 역학을 나타내는 데 더 자연스러운 표현일 수 있음을 시사한다.
• 모델 불가지론적 지표(Model-Agnostic Fingerprint): 본 분석은 임계 겔 상태를 위한 주파수 독립적인 모델 불가지론적 점탄성 지표, 즉 $(\partial G'' / \partial G')_{\varepsilon \to 0} = \tan((n-\kappa)\pi/2)$를 식별한다. 이 양은 임계 지수 $n$과 $\kappa$에 의해 유일하게 결정된다.
• 실험적 일치: 모델은 화학적 가교(비가역적) 및 물리적 결합(가역적) 시스템 모두에 대한 실험 데이터를 성공적으로 재현한다. 피팅은 임계점 근근처에서 매우 우수하며, 시스템이 솔 또는 겔 영역으로 멀어짐에 따라 편차가 증가하는데, 이는 유도된 식의 점근적 성질과 일치한다.

의의
본 논문의 주요 의의는 분수 점탄성 모델링과 겔화의 근본적인 스케일링 법칙 사이의 단절을 해결하는 데 있다. 연속성을 엄격한 물리적 제약으로 통합함으로써, 본 연구는 하이퍼스케일링 관계를 경험적 관찰에서 이론적 명령으로 격상시킨다. 프라바카르 지수 $\gamma$를 측정 가능한 물리량(점도 지수 $s$와 동일)으로 식별한 것은 수학적 피팅 매개변수와 클러스터 크기 분포 및 완화 스펙트럼 사이의 근본적인 물리 사이의 새로운 연결 고리를 제공한다. 마지막으로, 모델 불가지론적인 주파수 독립적 지표 $(\partial G'' / \partial G')$의 식별은 겔화의 화학적 성질이나 경로에 관계없이 임계 겔 상태를 특징짓는 강력하고 보편적인 기준을 제공한다.