1. 引言
液晶是一种介于流动液体和固态晶体之间的物质状态。它一方面像流体一样可以自由移动,另一方面又呈现出与晶体类似的各向异性,这一独特的物理性质通常来源于液晶分子几何上的各向异性,如棒状、盘状和香蕉形状等 [1]。液晶具有几个相,例如低温下的向列相和高温下的无序相。在向列相中棒状分子没有位置次序,但它们会自我对齐,具有长程定向次序,稳定的向列相和无序相可以共存。不同相之间的相变引发了各种有趣的数学问题。刻画液晶的相变现象以及动力学规律最有力的工具就是选取适当的数学模型,物理学家和数学家先后建立了各种数学模型。经典的液晶模型有Oseen-Frank模型和Ericksen模型,但是它们在刻画液晶时也存在不足之处,例如Oseen-Frank模型只能刻画液晶点缺陷,Ericksen模型只能刻画液晶分子的单轴分布。针对这些问题,Landau-de Gennes从统计物理角度出发,对向列型液晶提出了能全面描述液晶物理现象的Q-张量模型,称为Landau-de Gennes理论 [2]。
在Landau-de Gennes理论框架下,Park等人从分析的角度研究了无序相–向列相界面的静力学。考虑弹性常数为小参数,费明稳等人利用匹配渐近展开方法,从非惯性的Beris-Edwards模型 [3] 形式推导出描述无序相–向列相相变的尖锐界面极限模型。对于耦合流体费明稳等人形式推导了非惯性Beris-Edwards模型的无序相–向列相尖锐界面极限模型 [4];对于不耦合流体费明稳等人严格证明了Landau-de Gennes张量模型(Q-张量形式的梯度流)的无序相–向列相尖锐界面极限。费明稳等人还严格证明了非惯性Beris-Edwards模型与尖锐界面极限模型之间光滑解的收敛关系 [5]。有几种动力学Q-张量模型描述向列型液晶的流动,这些理论通过封闭近似例如 [6] [7] 得到分子动力学理论,或直接通过变分方法得到,如Beris-Edwards模型和Qian-Sheng模型 [8]。
描述液晶中无序相–向列相界面问题存在两种模型框架:尖锐界面模型和相场模型。尖锐界面模型涉及到在分离无序相和向列相的移动界面上求解具有匹配边界条件的控制微分方程 [9] [10];相场模型 [11] [12] 即在相场变量(序参量) Q中用有限宽度的光滑过渡区来模拟无序相–向列相界面,在无序相和向列相区域对应于不同的序参量Q值,在过渡区序参量在两个平衡值之间不断变化。本文采用相场模型,通过匹配的渐近展开方法,从Landau-de Gennes 理论出发,在相变外区域作外展开,相变内区域作内展开,推导出无序相–向列相两相共存时不耦合流体情形下的尖锐界面极限。
2. Landau-de Gennes模型
Landau-de Gennes理论能够相当全面地描述介质的局部行为,它解释了更复杂的液晶现象,如线缺陷和双轴构型。该理论采用了一个对称迹零的矩阵(序参量)
来表征液晶分子的排序行为。物理上
可以理解为描述棒状液晶分子指向的概率分布函数f的对称迹零二阶矩:
其中
表示在点x方向与m平行的分子的微观分布,如若张量
的所有特征值都为零,则称张量为无序的。如若张量
有两个相等的非零特征值,则称张量为单轴的。如若张量
有三个不同的特征值,则称张量为双轴的。
在没有边界约束和外场的情况下,Landau-de Gennes自由能有如下形式:
(1)
其中
是取决于材料和温度的非负参数,而
是与材料相关的弹性常数,
是描述无序相–向列相相变的体积能密度,而弹性能量密度
描述空间非均匀性。有关详细的介绍可以参考 [2] [13]。
3. 体积能的临界点和线性化算子
如果
,则矩阵
称为体积能
的临界点,从 [14] [15] 中可以看到临界点有
以下特征。
命题3.1对于一些
,
。其中
或者s是
的解,即
或
此外,当
时,临界点
是稳定的。
给定一个临界点
,
附近的
的线性化算子
定义为
通过直接计算得出
(2)
线性化算子
的核空间是
的一个二维子空间,可以定义为
对任意
,其中
。
4. 非线性波映照流方程的尖锐界面极限
假设存在一个宽度为
的过渡区域,分隔两个区域
(无序相和向列相区域),其中
为向列相区域,
为无序相区域。设
为以过渡区域为中心的光滑表面,
为到尖锐界面的符号距离。非线性波映照流方程为:
(3)
其中f和Q分别定义为
4.1. 外展开
在
区域对
进行展开,得到
(4)
通过泰勒展开得到
(5)
其中
取决于
。
将(4)和(5)代入(3),令
项对应相等,得到方程:
(6)
(7)
(8)
结合(6)和命题3.1,得到
对于一些
,
。令
,方程(7)表明
,即
那么
,通过(2)得到
(9)
通过(8)和(9),得到
即对于任意
,
(10)
结合 [16] 中的命题2.3,简单地计算得到
其中分子场h通过Oseen-Frank自由能
定义为
其中
表示的弹性系数为
上述结合(10)计算得
其中
,
。特别地,如果
,
。有
(11)
此外,当系数
,
时,(11)可以被简化为波映照热流:
4.2. 内展开
在界面附近的内部区域,我们引入了收缩变量
,并对(3)的解进行内部展开:
(12)
其中
是到界面
的符号距离,且
。通过使用泰勒展开式和(12),得到
其中
取决于
。通过计算有
其中
因此在新的变量
下,方程(3)变为
(13)
其中
假设展开式
(14)
其中
定义在
上,
定义在
的一个邻域中的一个邻域中,方程
等价于
将(12)和(14)代入(13),令
项对应相等,得到方程:
(15)
(16)
假设
时,有
注意到
,令(15)与
进行缩并并对z作从
到
的积分,得到
将相容性条件代入(16)中,得到
(17)
其中
对于
,有
,因此方程(17)变为
(18)
接下来考虑特殊情况
时,方程(15)和(16)可以简化为
(19)
(20)
我们定义
方程(19)描述了过渡区附近Q的首阶方程:
方程(19)存在一个单轴解
,且
由以下ODE确定:
(21)
方程(21)的解为
然后有
意味着函数
独立于
,则方程(18)可以简化为
(22)
定义
,观察到对于任意
,
然后通过(20)和(22)得到
意味着对于任意
,
此外,我们知道
。因此在尖锐界面上,n满足Neumann条件:
其中
是尖锐界面的单位法向量。
5. 结语
本文通过使用匹配的渐近展开方法,针对液晶无序相-向列相相变问题,对于不耦合流体推导出非线性波映照流方程的尖锐界面极限:在向列相区域
,在尖锐界面上
,尖锐界面的演化由平均曲率流
确定。
NOTES
*通讯作者。