Self-consistent field theory of semiflexible nematics: Density-nematic coupling, anisotropic elasticity, and defect core sizes

기술 요약: 반강성 네마틱의 자기 일관적 장 이론(Self-consistent field theory)

문제 제기
소분자 열역학적 액정의 선형 응답는 잘 이해되어 있는 반면, 고분자, 미셀 응집체, 생체 필라멘트 등을 포함하는 라이오트로픽(lyotropic) 액정의 거동은 상당한 이론적 및 실험적 과제를 안겨준다. 핵심적인 문제는 라이오트로픽 시스템에서 관찰되는 큰 탄성 이방성의 기원과 사슬 밀도와 네마틱 질서 사이의 결합이다. 온사거(Onsager) 모델에 주로 기반한 기존의 미시적 이론들은 공간적으로 변화하는 반강성 네마틱의 탄성, 위상 결함, 그리고 밀도-질서 결합을 동시에 다루는 데 어려움을 겪고 있다. 또한, 거시적 연속체 이론(Oseen–Frank 및 Landau–de Gennes)은 일반적으로 미시적 분자 통계와의 명시적인 연결 없이 현상론적 계수를 도입하며, 이로 인해 탄성 이방성과 결함 코어 구조의 미시적 기원은 탐구되지 않은 채 남아 있다.

방법론
저자들은 반강성 네마틱을 웜라이크 체인(wormlike chains)으로 모델링하는 자기 일관적 장 이론(SCFT) 프레임워크를 사용한다. 시스템은 다음을 포함하는 해밀토니안(Hamiltonian)으로 정의된다:

1. 굽힘 에너지(Bending energy): 지속 길이($l_p$)와 윤곽 길이($L$)에 의해 지배됨.
2. 마이어-소프(Maier–Saupe, MS) 상호작용: 네마틱 질서를 유도하는 비등방적 배향 상호작용($u_2$).
3. 등방적 배제 부피(Isotropic excluded volume): 라이오트로픽 조건을 나타내는 용매 매개 상호작용($u_0$, 플로리-헉킨스 파라미터를 통해 용매 효과를 암시적으로 처리함).

연구는 세 단계로 진행된다:

• 균일 상태 분석: 저자들은 사슬 세그먼트 밀도($\rho$)와 단축 질서 파라미터($S$)에 의존하는 유효 자유 에너지 범함수를 유도한다. 이를 통해 등방성-네마틱 공존(binodal) 및 안정성 한계(spinodal) 경계를 결정한다.
• 선형 응답 분석: 약한 외부 장을 도입하여 밀도, 단축 질서 및 디렉터(director)의 섭동에 대한 시스템의 응답을 조사한다. 이를 통해 역량성도 행렬(inverse susceptibility matrix)을 얻으며, 여기서 프랭크 탄성 상수($K_1, K_2, K_3$)와 상관 길이를 추출한다.
• 결함 구조 시뮬레이션: +1/2 디스클리네이션(disclination)에 대해 전체 SCFT 방정식을 수치적으로 해결하여 디렉터 프로파일과 코어 크기를 분석하고, 강직한 사슬 한계와 반강성 영역을 비교한다.

주요 기여 및 결과

• 밀도-네마틱 결합 및 안정성:
  유도된 유효 자유 에너지는 밀도($\rho$)와 질서도($S$) 사이의 강력한 결합을 드러낸다. 저자들은 강한 밀도-네마틱 결합이 네마틱 상태를 불안정하게 만들 수 있음을 입증한다. 스피노달 불안정성의 위치는 배제 부피와 네마틱 상호작용 강도의 비율($u_0/u_2$)에 따라 결정된다. 이와 대조적으로, 디렉터 왜곡은 밀도 및 질서 변화와 고차 항에서만 결합하며, 선형 응답 영역에서는 효과적으로 디커플링된 상태로 남는다.

• 탄성 이방성 및 유연성:
  본 연구는 사슬의 유연성($L/l_p$)이 프랭크 탄성 상수를 어떻게 변화시키는지 정량화한다.

  • 강직한 사슬: $K_3$(bend)는 정렬 강도가 증가함에 따라 급격히 증가하는 반면, $K_1$(splay)과 $K_2$(twist)는 유사한 수준을 유지한다.
  • 반강성 사슬: 유연성이 증가하면 트위스트($K_2$)와 벤드($K_3$) 탄성이 강력하게 억제되는 반면, 스플레이($K_1$) 탄성은 상대적으로 약하게 영향을 받는다.
  • 교차(Crossover): 이러한 차별적 억제는 유연성이 증가함에 따라 벤드 지배형에서 스플레이 지배형 탄성으로의 교차를 초래하며, 이는 $K_1/K_3$ 비율의 역전($<1$에서 $>1$로)으로 입증된다.

• 상관 길이 및 이방성:
  선형 응답 분석은 역량도 계수의 대각 방향 구배 계수가 디렉터와 평행한 방향의 상관 길이가 더 짧음($K_{\parallel} < K_{\perp}$)을 시사함에도 불구하고, 밀도와 질서 사이의 결합이 유효 계수를 재규격화한다는 것을 밝혀냈다. 이는 디렉터와 평행한 방향의 상관 길이가 실제로 더 큰($\xi_{\parallel} > \xi_{\perp}$) 이방성을 초래한다. 이 발견은 등방성-네마틱 계면 및 결함 코어 구조에 대한 실험적 관찰과 일치한다.

• 결함 코어 크기:
  +1/2 디스클리네이션에 대한 수치 시뮬레이션은 코어 크기가 이방성임을 보여주며, 밀도 구배가 디렉터와 평행한 방향을 따라 반지름이 더 크다. 본 연구는 특징적인 결함 코어 크기가 평형 상관 길이($\xi_Q$)에 비례함을 확립하며, 이는 란다우-드 젠(Landau–de Gennes) 스케일링과 일치한다. 저자들은 라이오트로픽 크로믹스가 종종 비정상적으로 큰 코어를 보이는 점(잠재적으로 응집체 절단 때문)을 언급하면서도, 다분산성을 제외한 현재의 모델이 밀도 요동에 의해 구동되는 코어 크기와 상관 길이 사이의 스케일링 관계를 포착하고 있음을 명시한다.

• 스피노달 거동:
  시스템이 네마틱 스피노달(변수 $u_0$에 의해 제어됨)에 접근함에 따라, 상관 길이는 $\xi \sim (u_0 - u_{0,i})^{-1/2}$의 스케일링을 따르며 발산한다. 그러나 결함 구성에서는 국부적인 등방성 도메인(tactoid)의 형성이 스피노달에 도달하기 전에 네마틱 질서를 억제함으로써 이 발산을 선제적으로 차단한다.

의의
본 논문은 반강성 네마틱에서 분자 유연성, 밀도-네마틱 질서 결합, 그리고 탄성 이방성을 연결하는 통합된 미시적 프레임워크를 제공한다. SCFT 형식론 내에서 사슬 통계로부터 프랭크 탄성 상수를 직접 유도함으로써, 이 연구는 라이오트로픽 시스템에서 관찰되는 큰 탄성 이방성의 미시적 기원을 규명한다. 유효 계수를 도출한 결과는 결함 구조를 설명하는 데 있어 연속체 이론(예: Landau–de Gennes)의 사용을 정당화하는 동시에, 탄성 응답과 안정성 한계를 정확하게 예측하기 위해 밀도 요동과 사슬 유연성을 포함하는 것이 필수적임을 강조한다. 저자들은 본 프레임워크가 자기 조립된 다분산 사슬 및 가역적 응집체로의 향후 확장을 위한 기초를 제공한다고 결론짓지만, 현재의 이론적 분석 범위를 넘어서는 구체적인 새로운 실험 프로토콜이나 응용 분야를 제안하지는 않는다.