胆甾相液晶COMSOL仿真指南：介电张量设置与布拉格反射模拟

做光学仿真的人，迟早会碰液晶；做液晶仿真的人，迟早会碰胆甾相。这种材料在器件里太常见了——反射式显示、激光护目镜、圆偏振片、温度传感器，背后都靠它那套“挑偏振、挑波长”的反射机制。但真要在 COMSOL 里把胆甾相液晶的光反射与光透射算清楚，头一回动手的同行，十个有八个会卡在“介电张量到底该怎么旋转”和“为什么仿真出来的布拉格峰和理论公式对不上”这两件事上。

这篇内容我按自己实际跑通的全流程来讲：从胆甾相液晶的螺旋结构开始，解释它的选择性反射从哪来，然后落到 COMSOL 的建模方案、材料设置、端口与边界条件、网格与求解，最后附上我做参数扫描和排查问题时的经验。无论你是刚接触液晶光学的新人，还是已经会用 COMSOL 但想验证手性光子结构的老手，这套思路基本可以直接照搬。

1. 胆甾相液晶到底“甾”在哪：先弄懂你算的是什么

1.1 螺旋结构与光学周期

胆甾相液晶其实是向列相液晶加了一条“手性约束”：分子在每一小层内平行排列，但相邻层之间的指向矢会绕螺旋轴连续旋转。从侧面看过去，就是一根根分子方向像拧麻花一样绕着轴线转。

这里有一个必须正面说清的概念：螺旋的螺距 p 有两种定义习惯。一种是把“指向矢旋转 360° 所经过的距离”称为螺距，另一种则把“旋转 180° 的距离”称为螺距。你在查文献时会看到两种写法都有，千万别直接用错公式。

在 COMSOL 建模时，我的建议是把螺距参数设为 p 代表旋转 360° 的距离，也就是指向矢

n(x) = (cos(2πx/p), sin(2πx/p), 0)

这里的 x 是螺旋轴方向。注意一个关键物理事实：液晶分子有一个 n 和 -n 等价的对称性，所以介电张量沿 x 方向实际变化的周期是 p/2，而不是 p。这个 p/2 就是光学上真正“看见”的周期。布拉格反射的中心波长

λ₀ = n_avg · p

而不是 2·n_avg·(p/2) 的简单叠加结果，具体推导后面会讲。如果你一开始就用光学周期 p/2 去换算波长，很容易把整个反射带的位置算偏。

1.2 选择性布拉格反射的三个关键参数

胆甾相液晶对光的反射并不是普通反射镜那种“什么都弹回去”，它有很强的选择性，被人称为选择性布拉格反射。工程上关心的核心参数有三个。

中心波长：正入射时，满足布拉格条件的波长近似为

λ₀ = n_avg · p

其中 n_avg 是液晶平均折射率。对于由寻常光折射率 n_o 和非寻常光折射率 n_e 组成的双折射材料，常见近似取 n_avg ≈ (n_e + n_o)/2。波长在这个位置附近的光会强烈反射。

反射带宽：反带宽度与双折射直接相关，近似为

Δλ ≈ (n_e - n_o) · p

双折射越大，反射带越宽。这就是为什么通常在仿真里把 n_e 和 n_o 拉开差距后，能明显看到反射谱的宽谷变宽。如果要实现宽带反射，基本思路就是提高双折射或者叠加多个螺距不同的结构。

反射率与层数：反射率不是一上来就接近 100% 的，它取决于螺旋周期数 N。周期数越少，反射率越低，带内可能只有几十个百分点的反射；周期数足够多时，带内反射率趋于饱和，接近 100%。在 COMSOL 里这一点特别直观：把模型长度从 5 个螺距拉到 20 个螺距，反射谱的深度会肉眼可见地变化。

1.3 反射和透射的偏振密码

胆甾相液晶最迷人的地方是它对圆偏振有“挑食”行为。

当一束光沿着螺旋轴传播时，只有与螺旋手性匹配的那一种圆偏振光会被布拉格反射，另一种圆偏振光基本不受影响，直接透射过去。比如右手螺旋的胆甾相液晶反射右旋圆偏振光，左旋圆偏振光则透射。在反射带内，反射光的旋向与入射光相同；在反射带外，两种圆偏振都能透射。

这句话听起来简单，但做仿真时很容易踩坑：如果用两个正交线偏振分别入射，会看到两者的响应几乎一样，很难直接感受到“手性选择”。只有当你在后处理里把两个正交线偏振的复振幅按 90° 相位差组合起来，或者说干脆用圆偏振模式激励，才能把胆甾相液晶的偏振特性真实还原出来。

提示：圆偏振光在这个结构里才是“本征模式”，线偏振光进入后会被劈裂成两个圆偏振分量，在带内经历完全不同的命运。所以建模阶段就要想清楚，你最终要提取的是线偏振响应还是圆偏振响应，这两种结果不能混着看。

2. COMSOL 建模的整体思路与接口选型

2.1 为什么选 COMSOL 而不是手写传输矩阵

很多人问过我一个问题：胆甾相液晶的一维光学性质，用传输矩阵或者商业薄膜设计软件几行代码就能算，为什么还要上 COMSOL？

我的回答是：如果你的目标只是拿到一维正入射的反射率、透射率曲线，确实不需要 COMSOL。但如果你接下来要处理斜入射、非均匀螺距、缺陷层、外加电场导致指向矢重排、多层膜叠加、手征结构与其他微纳结构的耦合，传输矩阵的推导成本会迅速变大，边界条件也容易写错。COMSOL 的优势是你可以在一个物理接口里统一处理这些复杂情况，尤其是场分布的可视化——想看带内电磁场如何在螺旋结构中指数衰减，只有 FEM 场图最直观。

而且 COMSOL 里做完一个基础模型后，参数化扫描、优化模块、与热场或结构力学场耦合都是现成的。我的做法是：先用 COMSOL 验证光路，再决定是否要进一步做多物理场。这套基础模型的价值比手写几个公式要大得多。

2.2 选哪个物理接口：ewfd 与频域求解

在 COMSOL 中，做这类光学仿真使用的是波动光学模块里的“电磁波，频域”接口，也就是物理接口简写 ewfd。它求解的是频域亥姆霍兹方程组，能正确处理各项异性介电张量、周期边界和多端口问题。

不要选“射线光学”，那是几何光学范畴，看不到干涉和布拉格效应。也不要选“电磁波，波束包络”，虽然它适合长距离传播，但对亚波长螺旋结构没什么优势。ewfd 接口配合频域研究，足够覆盖正入射、斜入射和参数扫描。

研究类型方面，新建研究时选择“频域”，然后在扫描设置里定义波长或频率范围。我习惯用波长作为扫描变量，因为液晶光学通常是以“xx nm 的反射带”来描述的，直接扫波长更符合阅读习惯。COMSOL 里可以通过表达式把波长转成频率，例如 freq = c_const/λ，这里的 c_const 是 COMSOL 内置真空光速常量。

2.3 模型维度与单胞策略

胆甾相液晶的光学结构，严格说是三维的：分子指向矢在面内旋转，光沿面内某一方向传播。但结构在垂直于螺旋轴的平面内是均匀的，这就给了我们很大的简化空间。

有两个主流方案。

方案一是纯三维单胞：取一小段螺旋结构，横向尺寸很小，只相当于螺距的十分之一左右，侧面加周期性边界条件。这个方案配置简单，电场可以自由设置为任意偏振方向，最适合做圆偏振分析和斜入射。计算量也非常小，因为横向网格只划分三五个单元就够了。

方案二是二维模型：在光传播方向和面内一个偏振方向上建立几何。缺点是 COMSOL 二维波动光学默认只能处理一个面外偏振分量和面内偏振分量，如果你想同时得到两个正交偏振的响应，要么建两个接口并行计算，要么还是要回到三维。所以我个人建议，除非你只关心单一偏振，否则直接用三维单胞更省心，别贪二维的这点计算量。

注意：三维单胞的横向尺寸一定要远小于工作波长，否则会产生不希望的横向模式。通常取 p/10 到 p/20 就够，加上周期性边界条件后，物理上等价于无限大的均匀薄膜。横向网格太粗虽然省时间，但对后处理和端口模式计算都不利。

3. 实操：从几何到材料的完整配置

3.1 几何搭建与周期性边界

直接建立一个长方体，定义全局参数：

几何的长宽高就设置为：

• 螺旋轴方向长度 L = N·p
• 横向尺寸 a = p/10
• 另一个横向尺寸同样为 a

模型的 x 方向从 0 到 L，光沿 +x 入射。y 和 z 方向分别取 0 到 a。

两侧与光传播方向垂直的面是端口面：x = 0 是入射/反射端口，x = L 是透射端口。四个平行于 x 轴的侧面都加上周期性边界条件。

具体操作是：在 ewfd 接口下添加“周期性条件”，类型选“周期性”，手动指定源边界和目标边界。由于这里入射光垂直照射，横向波矢分量为 0，Floquet 波矢就是 0，等价于普通的周期性连接。如果是斜入射，则需要在这里填上对应的控制波矢表达式。

3.2 材料参数：旋转介电张量怎么设

这是整个模型里最容易被卡住的一步。COMSOL 默认的材料“折射率”往往是各向同性的，而胆甾相液晶必须用各向异性介电张量描述。

我的做法是新建一个“空材料”，然后在“相对介电常数”里选择“各向异性”或直接按矩阵形式填写相对介电常数的分量。定义方法如下。

液晶主轴折射率对应介电常数：

ε_e = n_e²
ε_o = n_o²

在第 3.1 节的几何坐标里，假设分子指向矢在 xy 平面内旋转：

n_x = cos(θ)
n_y = sin(θ)
θ = 2πx/p

相对介电张量的分量可以写成：

ε_xx = ε_o + (ε_e - ε_o) · cos²(θ)
ε_yy = ε_o + (ε_e - ε_o) · sin²(θ)
ε_xy = ε_yx = (ε_e - ε_o) · cos(θ) · sin(θ)
ε_zz = ε_o

其余分量取 0。

这里最关键的是理解为什么 ε 分量会包含 cos²、sin² 和 cos·sin 的混合项。胆甾相液晶的介电张量是对称张量，主轴与分子指向矢对齐。当你把它从液晶主轴坐标系转换到实验室坐标系时，就会出现旋转张量的二次项。如果直接只用主对角分量而不加交叉项，电场会被错误地解耦，反射谱会完全失真。

实操提醒：在 COMSOL 材料设置里，把上述表达式中的 x 变量写成坐标系变量名。三维模型默认是 x、y、z，如果你的模型旋转了坐标系，就要注意变量名是否仍然是 x。做材料表达式时我习惯先把坐标系列出来检查一遍，省得到时候出现“为什么材料旋转不对”的排查难题。

3.3 端口设置与圆偏振入射

在 ewfd 接口中添加两个“端口”边界条件，分别施加到 x = 0 和 x = L 的两个边界。

端口 1 在入射面，设置为“波激励”，用于发射入射波并提取反射 S 参数。端口 2 在透射面，设置为常规端口，用于提取透射 S 参数。

对于正入射的平面波端口，端口类型通常选“平面波”。如果只是做线偏振验证，可以把偏振方向分别设为 y 方向和 z 方向，各算一次。

如果希望直接得到圆偏振响应，我推荐使用“用户定义”端口模式。在端口设置里把模式类型选为“用户定义”，然后在模式电场里输入复数矢量。对于沿 +x 传播的右旋圆偏振光，一种可行写法是：在端口面上设定电场矢量为 (0, 1, j)，其中 j 是 COMSOL 中的虚数单位。左旋圆偏振则设为 (0, 1, -j)。具体取值还要看你的坐标定义和传播方向，如果发现手性反了，取共轭即可。

但这里有一个经验之谈：头一回做时，我建议先用两个线偏振端口分别计算，把结果解出来后再去合成圆偏振。因为线偏振响应最容易对理论公式、最容易排查网格问题。等确认基础模型没问题，再改成用户定义圆偏振模式做最终计算，两条路可以互相验证。

3.4 网格划分的教训

网格划分是这类周期螺旋结构最容易翻车的地方，也是最容易用计算量换教训的地方。

螺旋轴方向也就是 x 方向，必须解析介电张量的旋转周期。由于介电张量实际振荡周期是 p/2，这个周期内至少要有足够的单元数。我的经验是最少 10 到 20 个网格单元。假如 p = 400 nm，那么每 20 nm 至少要有一个网格，如果取 10 个周期，x 方向网格总数就在 200 到 400 左右。这个数量级对三维单胞来说毫无压力。

横向 y、z 方向只需要 3 到 5 个单元，因为结构在横向上均匀，场分布也很平缓。千万别学普通三维结构的习惯，把横向网格划得很细，那是白费计算量。网格单元选择映射网格或者扫掠网格都很适合这个几何：x 方向加密，y、z 方向保持少量单元即可。

注意：如果你在反射带的边缘或者高折射率差场景下计算，默认的二阶单元可能出现轻微数值冗余，建议对 x 方向使用更细的网格，或者改用更高质量的网格序列。反射谱出现非物理毛刺时，第一反应不是调求解器，而是加密网格。

4. 跑通核心仿真：反射率与透射率曲线、扫描与验证

4.1 频率扫描与 S 参数提取

在频域研究设置中，把扫描参数设为波长 λ，范围建议覆盖理论中心波长的两侧。以 p = 400 nm、n_avg = 1.6 为例，中心波长大约在 640 nm，那么扫描范围设在 500 nm 到 800 nm 就足够看到完整的反射带。

计算完成后，提取全局变量。COMSOL 中端口 S 参数可以直接用：

• 反射率 R = abs(ewfd.S11)^2
• 透射率 T = abs(ewfd.S21)^2

注意 S11 和 S21 本身是复数，在绘图或全局计算中不要忘了取模平方。很多人刚开始画出的曲线是复数实部，看起来一团混乱，就是因为没有取 abs。

用一维绘图组画这两条曲线，横轴设为波长，纵轴分别为 R 和 T，就能看到典型的胆甾相液晶光谱：反射带内 R 接近 1、T 接近 0；反射带外 R 接近 0、T 接近 1。由于模型没有添加材料损耗，带外透射率应该很接近 100%，带内反射率随着周期数增加快速逼近 100%。

4.2 偏振分解与手性验证

如果用的是线偏振端口，下一步需要做偏振分解。方法是在两个正交线偏振计算中提取 S11 和 S21，再在后处理中构造圆偏振响应。

具体做法是：对入射圆偏振光，反射的同旋分量通过两个正交线偏振的反射复振幅求组合。例如：

R_right = abs(S11_y + j·S11_z)^2 / 2

这里 j 是虚数单位，系数 1/2 来自圆偏振的归一化。如果计算结果正确，你会发现入射右旋圆偏振和入射左旋圆偏振得到的反射谱完全不同：一个在带内强烈反射，另一个几乎全透。

这一步是我最喜欢做的验证，因为它直接体现了胆甾相液晶的手性响应。如果两种圆偏振算出来一样，那就说明你的材料介电张量交叉项设错了，或者应该检查网格方向是否在旋转过程中被数值抹平。

4.3 参数化扫描：螺距、双折射和层数的影响

基础模型跑通后，参数化扫描才是 COMSOL 真正的价值所在。

螺距扫描：把 p 从 350 nm 扫到 450 nm，步长 10 nm。你会看到反射带中心波长跟着 p 线性移动，这也是胆甾相液晶温度传感器的物理基础——温度改变螺距，光谱位置就变化。扫描结果可以做一个二维图：横轴波长、纵轴螺距、颜色为反射率，布拉格反射带会形成一条清晰的斜亮带。

双折射扫描：把 n_e 从 1.6 扫到 1.8，保持 n_o 不变。反射带宽随 Δn = n_e - n_o 变宽，中心波长略有移动。这个特性对反射式显示器的色域设计很有用，想加宽反射带就把双折射做大。

层数扫描：把 N 从 2 扫到 30。反射带深度先快速增长，后期饱和。饱和点对应的周期数取决于双折射，Δn 越小，需要的周期数越多。比如 n_e - n_o = 0.2 时，大概 10 个周期就能接近饱和；Δn 只有 0.05 时，可能要几十个周期才能看到明显的强反射。

每完成一组扫描，我都建议顺手画一张反射带深度随 N 的曲线，然后和理论估计做对比。误差在几个百分点以内说明模型没问题，误差很大时优先检查端口相位参考面和材料设置。

5. 常见问题与排查实录

5.1 反射谱中心波长与理论对不上

这是最常遇到的问题，几乎每个新手都会遇到一次。反射带中心不是理论 λ₀ = n_avg·p 的位置，而是差了一倍或者是某个奇怪偏移。

排查顺序有三步。

第一确认螺距定义：检查你参数里 p 到底是“指向矢转 360°”还是“指向矢转 180°”。很多教程里把 p 写成半周期，套公式时就会差两倍。第二检查介电张量表达式：看 θ = 2πx/p 是否正确，如果误写成了 θ = 2πx/(p/2)，实际旋转速度翻倍，布拉格波长也会翻倍。第三检查网格是否解析足够：如果每个光学周期只有两三个网格，反射带会明显偏移甚至消失。先加密到每周期 20 个网格再判断。

5.2 透射率大于反射率导致能量不守恒

无损耗模型里，理论上 R + T = 1。如果你算出 R + T 明显偏离 1，比如只有 0.85，那问题通常出在端口或边界上。

最常见的原因是周期边界与端口面没匹配好，导致部分能量从侧面漏走。检查周期性条件是否完整施加到四个侧面，并且源边界和目标边界的映射方向一致。其次，检查端口模式是否与企业号的参考面错位。如果实际模式与入射模式不完全正交，S 参数可能会出现串扰。

经验：我在单胞横向尺寸设置过大时，也遇到过能量守恒偏差。横向尺寸一旦达到与波长同数量级，单胞里会出现高阶横向模式，端口假设的平面波模式就不再是完备描述，这时把横向尺寸缩回 p/10 量级即可。

5.3 圆偏振响应不对称

当我们构造左旋和右旋圆偏振入射时，发现两者的透射谱几乎一样，完全没有手性选择。这个问题的根源基本都是材料介电张量里 xy 交叉项写错了或者手性方向取反。

检查方式：在 COMSOL 后处理里画某一个平面的介电张量主方向云图，看指向矢是否确实沿 x 方向旋转，并且旋转方向是否符合预期。如果 θ = +2πx/p 和 θ = -2πx/p 的仿真结果只是把明暗带翻转，那就说明结构本身是手性的，只是初始手性定义反了。此时把 θ 加个负号重新扫描即可。

5.4 反射带边缘出现毛刺或振荡

出现这个现象通常有两种原因：网格不够密，或者扫描步长太大。带边缘附近场梯度很大，如果网格分辨率不足，数值色散会带来虚假振荡。先加密网格，再把扫描步长从 5 nm 缩到 1 nm。如果毛刺仍存在，检查是否在带边缘波长附近出现了数值共振，这可以通过观察场分布图来辅助判断。

6. 一些值得继续做的事

基础的正入射模型跑通后，这个模型的价值才开始释放。我建议你按自己的方向挑几条路继续做。

斜入射是需要经常用到的扩展。把周期性边界条件的 Floquet 波矢从 0 改成与入射角度相关的表达式，同时调整端口的入射方向和偏振，就可以得到角度依赖的反射谱。胆甾相液晶在斜入射时布拉格波长会蓝移，反射带宽也会变化，对宽角器件设计很关键。

缺陷态和可调结构也是热门方向。在螺旋结构中间插入一层各向同性介质，会在禁带里打开一个窄透射通道，对应缺陷态。这个模型对于窄带滤光和传感器设计特别有用。如果再把液晶材料替换成电调谐材料，或者耦合一个结构力学场让螺距变化，就能直接仿真可调滤波器。

另外，如果你手头有真实的反射率测量数据，强烈建议做一步对比验证：先用 COMSOL 算理想结构，再根据工艺误差加入随机螺距抖动，看反射谱如何退化。这个容差分析在工程上非常实用，也是手写传输矩阵很难快速完成的。

从我个人经验来说，把胆甾相液晶的 COMSOL 模型吃透以后，再看其他手性光子结构都会感觉顺畅很多，因为介电张量旋转、布拉格条件、偏振选择这些核心概念是共通的。刚开始做的时候遇到波长对不上、偏振不对称的问题，不必怀疑软件，多数情况都是定义和细节没对齐。把每一步的设置变量统一记清楚，这套模型就能成为你后续各种液晶光学设计的基础工具。