강제 고정 경계조건에서의 액정 거동: 오센‑프랭크 한계와 소형 영역 이론의 고차 보정

초록

본 논문은 2차원 제한된 영역에서 Rapini‑Papoular 표면 에너지를 유지한 채 강제 고정(Strong anchoring) 경계조건을 재검토한다. 작은 영역(ε→0)에서는 경계 데이터의 평균이 주도해 무향성 용융이 발생할 수 있음을 보이고, 큰 영역(ε→∞)에서는 표면 고정이 O(1/ε) 수준의 교정항을 만든다는 점을 Dirichlet 조건이 예측하는 O(1/ε²)와 대비한다. 수치 실험은 정사각형·원형 웰에서 결함 형태의 차이를 명확히 보여준다.

상세 분석

이 연구는 Landau‑de Gennes 이론에 Rapini‑Papoular 형태의 표면 에너지를 명시적으로 포함함으로써 전통적인 Dirichlet 강제 고정 모델의 한계를 극복한다. 핵심은 외삽 길이 l_ex와 코히런스 길이 ξ가 동일한 차원에서 비교되는 스케일링(γ=ξ/l_ex=O(1))을 채택한 점이다. 이 스케일링은 실제 나노·마이크로 규모의 제한된 도메인에서 l_ex가 ξ보다 크게 되거나 비슷해지는 상황을 반영한다.

소형 영역 한계(ε=h/ξ→0)에서는 전반적인 해가 상수이며, 그 상수값은 경계 텐서 Q_b의 평균값으로 고정된다. 이는 경계 조건이 서로 충돌하는 대칭적 기하학(예: 정다각형의 접선 고정)에서 평균이 0이 되어 Q=0, 즉 등방성 용융이 일어나는 메커니즘을 제공한다. 첫 차수 교정은 Neumann 형태의 포아송 방정식으로, 경계에서의 Robin 조건(∂_ν q = -εγ(q−q_b))이 내부에 전달되는 방식으로 나타난다. 고차 항들은 일련의 포아송 방정식과 적분 제약조건을 통해 순차적으로 해결되며, 각 차수마다 상수항이 고유하게 결정된다.

반면 대형 영역 한계(ε→∞)에서는 외부(벌크)와 내부(결함) 영역으로 나뉘는 특이 섭동 문제가 발생한다. 외부 영역에서는 Oseen‑Frank 이론이 지배하지만, 표면 고정이 존재하면 지배적인 교정항이 O(1/ε) 규모로 나타난다. 이는 Dirichlet 고정이 예측하는 O(1/ε²)와 근본적으로 다른 스케일이며, 결함 코어 주변에서의 director 필드가 더 부드럽게 변한다는 물리적 의미를 가진다. 특히, Robin‑type 경계조건은 결함 코어를 완화시켜 실제 실험에서 관찰되는 매끄러운 결함 구조와 일치한다.

수치 시뮬레이션은 정사각형과 원형 웰에 접선 고정을 적용한 경우를 다룬다. Dirichlet 조건을 사용하면 코너에서 급격한 방향 전이가 발생해 강한 점 결함이 형성되지만, Rapini‑Papoular 기반 Robin 조건에서는 결함이 경계에 걸쳐 연속적인 변화를 보이며 코어 크기가 ξ에 비례한다. 또한, 작은 영역에서는 평균 p가 0이 되는 경우 전역적으로 등방성 상태가 나타나는 반면, 큰 영역에서는 전체 위상 전위(위상 차수 d)에 따라 최소 결함 수와 전하가 결정된다.

결론적으로, 표면 에너지를 명시적으로 포함한 강제 고정 모델은 기존 Dirichlet 기반 Oseen‑Frank 근사보다 더 정확한 고차 교정을 제공하고, 특히 나노·마이크로 스케일의 제한된 기하학에서 결함 형성 및 안정성을 예측하는 데 필수적이다.