一、单轴液晶的偏振特性
Polarization properties of a uniaxial liquid crystal
在这里插入图片描述
在这个例子中,我们将模拟一个简单的单轴 LCD 层。 有关更复杂的 LC 方向,请参阅扭曲向列 LCD。
1.1 模拟设置
由液晶取向引起的各向异性导致沿液晶指向矢偏振的光以与沿垂直方向偏振的光不同的速度传播。 这种双折射导致入射光和透过 LC 层的光之间的偏振发生变化。
为了在 FDTD 中模拟这种效果,我们将使用各向异性材料和 LC 网格属性对象,它允许我们为 LC 层指定任意方向。 在这种情况下,我们所要做的就是在 LC 属性对象中设置与所需 LC 方向相对应的角度“theta”和“phi”。 在 LCD.fsp 中,我们从在 x-y 平面上以 45 度偏振的平面波开始,并使用偏振椭圆分析组(polarization ellipse analysis group)来研究透射光的偏振变化。 偏振椭圆分析对象可以在对象库( object library)中找到,在远场投影(far field projections)部分。
1.1.1 各向异性材料
1.1.2 LC 网格属性对象
位置在Attributes
参数编辑
液晶指向为y轴
1.1.3 平面波光源
45度偏振
1.1.4 仿真区域
大小是个大长条,mesh accuracy设置为6,很高的精度
边界条件:x,y方向上周期边界条件,z上完全吸收、无反射。
除此之外,还设置了两个mesh,上下各一个。
1.1.5 偏振椭圆分析组
1.1.5.1 位置
1.1.5.2 Setup
varibables:
Script:
############################################### Polarization ellipse# This script sets up the monitor size 该脚本设置监视器大小## Input properties# x,y span: span of the monitor## Tags: far field polarization ellipse grating# 远场偏振椭圆光栅# Copyright Lumerical Solutions Inc############################################### simplify variable names by removing spacesx_span = %x span%;y_span = %y span%;select("field");set("x",0);set("y",0);set("z",0);set("x span",x_span);set("y span",y_span);
1.1.5.3 Analysis
Variables:
4个参数:
f:感兴趣的频率
Gn、Gm:感兴趣的光栅阶次 (Gn,Gm)
make plots: 1 表示绘图,0 不绘制
2个输出属性:
major_angle:偏振椭圆长轴的角度
比率:长轴长度除以椭圆短轴长度
Script:
############################################### Polarization ellipse# This script calculates the polarization of a grating order. 此脚本计算光栅阶的极化。# The result is plotted in terms of the polarization ellipse. 结果以偏振椭圆的形式绘制。## Input properties# f: frequency of interest# Gn,Gm: grating order of interest (Gn,Gm)# make plots: 1 to make plots, 0 otherwise## Output properties # major_angle: angle of major axis of polarization ellipse# ratio: length of major axis divided by length of minor axis of ellipse## Tags: far field polarization ellipse grating# 远场偏振椭圆光栅# Copyright Lumerical Solutions Inc############################################### simplify input variable names by removing spacesmake_plots = %make plots%;# find frequency closest to that specified# by the user# 找到最接近指定频率的频率# 由用户f_vec=getdata("field","f"); fi=find(f_vec,f); # Do grating calculation, and separate# into each field component# 进行光栅计算,并分离进入每个场分量G_vec = gratingpolar("field",fi);# 返回球座标的向量场信息n = gratingn("field",fi); # 返回第一维级的级数m = gratingm("field",fi); # 返回第二维级的级数# n和m都是0ux= gratingu1("field",fi);uy= gratingu2("field",fi);Gtheta = pinch(G_vec,3,2); # 相当于只取了第二种值的所有信息,也就是thetaGphi = pinch(G_vec,3,3);G=abs(Gtheta)^2+abs(Gphi)^2;# grating order strength# find index of selected grating order 查找所选光栅阶次的索引ni = find(n,Gn);mi = find(m,Gm);# convert spherical coordinates polarization into # S,P polarization and select grating order# 将球坐标偏振转换为 S,P 偏振并选择光栅阶数Gs=Gphi(ni,mi);Gp=Gtheta(ni,mi);# Calculate polarization ellipse 计算极化椭圆kappa = linspace(0,360,400)*pi/180;Gs = real( Gs*exp(1i*kappa) );Gp = real( Gp*exp(1i*kappa) );# measure parameters from ellipse 从椭圆测量参数diameter = sqrt( (Gs)^2+(Gp)^2 );major_axis = max( diameter);minor_axis = min( diameter );ratio = major_axis/minor_axis;pos_major = find(diameter,max(diameter));major_angle = atan2(Gs(pos_major),Gp(pos_major))*180/pi;if(major_angle < -90) {major_angle = major_angle+180; }if(major_angle > 90) {major_angle = major_angle-180; }# 把角度限制在-90-90之间Gs=Gs/major_axis;Gp=Gp/major_axis;theta0 = acos(sqrt(1-ux(ni)^2-uy(mi)^2))*180/pi;phi0 = atan2(uy(mi),ux(ni))*180/pi; if (make_plots) {# plot grating order strenghtsimage(n,m,G ,"Grating number n","Grating number m","Grating order strength");# plot ellipse, and ellipse parametersplot(Gp,Gs,"P","S","polarization ellipse (N,M)=("+num2str(Gn)+","+num2str(Gm)+")" +endl+"(theta,phi)=("+num2str(theta0)+","+num2str(phi0)+")" +endl+"axis ratio = " + num2str(ratio) +endl+"major axix angle = "+num2str(major_angle) );setplot("x min",-1); # 设置图形属性setplot("x max", 1);setplot("y min",-1);setplot("y max", 1);}
gratingpolar及相关函数
linspace
索引用kappa(1),或者kappa(1:6),kappa(1:2:6).
li*kappa的效果是把kappa的数据写到虚部,比如1变成0+1i
atan2
atan2()是对象限敏感的,因此需要两个参数
image
endl
1.1.6 sweep
1.2 结果
我们将查看两种 LC 方向的结果:第一种是 LC 导向器指向 z 轴(theta = 0,phi = 0),第二种是导向器指向 y 轴(theta = 90° , phi = 0)。 在第一种情况下,沿 x 和 y 方向偏振的光仍将以相同的速度传播,因此我们看不到输出光的偏振有任何变化。 在第二种情况下,在 y 方向上偏振的光将以与 x 不同的速度传播,我们将看到由于两个分量之间的相位差异而导致的偏振变化。 该相位差与 LC 层的厚度成正比:
分析组绘图结果:
为了确定将偏振改变 90 度(即从 45 度到 -45 度)所需的最佳厚度,我们建立了一个参数扫描项目,将 LCD 层的厚度从 2 到 3 微米扫描,并跟踪椭圆率ellipticity (即长轴和短轴之间的比率)作为这种变化的函数。 参数sweep的结果如下图(在sweep项目上右击,选择Visualize-ratio可以看到)。我们可以看到最佳厚度在2.75um左右。
这也是对于Δθ=π\Delta\theta=\piΔθ=π所期待的结果,因为,h=λ2Δn=2.75μmh=\frac{\lambda}{2\Delta n}=2.75\mu mh=2Δnλ=2.75μm.
PS:45度到-45度,难道不是应该90度,π2\frac{\pi}{2}2π吗
还有这里这个Δn=0.1\Delta n=0.1Δn=0.1,查了查应该是双折射率差,1.6-1.5=0.1
二、螺旋扭曲液晶的偏振特性
在本例中,我们将对扭曲向列 LCD (TN-LCD) 进行建模。 对于简单的单轴 LCD,请参阅简单 LCD。
2.1 模拟设置
当 TN-LCD 上没有施加电压时,液晶以扭曲配置取向,并且 LC 导向器从 LC 层的底面到顶部均匀地旋转 90 度。由于光以平行于底部导向器(本例中的 x 轴)的偏振进入 LC,光的线性偏振将大致跟随导向器的旋转,因此透射的光将在 y 方向上偏振.当施加电压时,外电场迫使 LC 指向矢变得均匀且平行于传播方向,光的偏振态不会发生变化。
为了模拟 FDTD 中 LC 方向的 90 度扭曲,我们将使用各向异性材料与 LC 导入网格属性相结合。脚本 LCD_twist.lsf 设置螺旋扭曲配置中的 LC 方向,如文件 LCD_twist.fsp 中所示。我们从与底部 LC 导向器(沿 x 轴)具有相同偏振的平面波源开始,并使用 2D DFT 监视器跟踪光在传播通过 TN-LCD 时的偏振变化。
2.1.1 材料设置
各向异性材料和网格属性设置
这个LC attribute跟上一个案例的设置截面又不一样了
2.1.2 光源
550nm输入
2.1.3 脚本 LCD_twist.lsf
# delete existing attributeselect('LC attribute');delete; # 删除所选择的实体对象# define x/y/zx = 0;y = 0;z = linspace(0e-6,5e-6,100);X = meshgrid3dx(x,y,z);Y = meshgrid3dy(x,y,z);Z = meshgrid3dz(x,y,z);n = matrix(length(x),length(y),length(z),3); # 初始化一个矩阵# define the orientation function 定义方向函数n(1:length(x),1:length(y),1:length(z),1) = cos(Z*pi*1e5);n(1:length(x),1:length(y),1:length(z),2) = sin(Z*pi*1e5);n(1:length(x),1:length(y),1:length(z),3) = 0;# create dataset containing orientation vectors and position parameters# # 创建包含方向向量和位置参数的数据集LC=rectilineardataset("LC",x,y,z);LC.addattribute("u",n);# add LC import grid attribute 添加 LC 导入网格属性addgridattribute("lc orientation",LC); # lc orientation在attribute下拉栏里setnamed("LC attribute","nz",50); # set resolution 设置分辨率
2.1.3.1 meshgrid3dx
2.1.3.2 rectilineardataset
2.1.3.3 addattribute
2.1.3.4 addgridattribute
2.2 结果
在这里,我们可以通过 TE 分数定义极化(沿 z 的每个点):
∣Ex∣2∣Ex∣2+∣Ey∣2\frac{|Ex|^2}{|Ex|^2+|Ey|^2}∣Ex∣2+∣Ey∣2∣Ex∣2
脚本 LCD_twist_analysis.lsf 绘制了 TN-LCD 层中的“TE 分数与 z”:
实际操作时,TN-LCD层会夹在两个偏光片之间(偏光差为90度)。 在没有外部场的情况下,光线将在 LCD 层中发生 90 度的偏振变化,因此将被允许通过顶部偏振器(因此处于“开启”状态)。 当施加外场时,晶体将与场对齐,因此入射光的偏振不会发生变化,顶部偏光片将不允许任何光通过(因此“关闭” 状态)。
2.2.1 field的结果
Ex随着z的增加,而逐渐减小
Ey随着z的增加,而逐渐增大
2.2.2 脚本 LCD_twist_analysis.lsf
closeall;clear;E = getresult('field','E');Ex = E.Ex;Ey = E.Ey;z = E.z;# track the TE fraction in side the rotator 跟踪旋转器侧的 TE 分数TE_frac = matrix(length(z));for (i=1:length(z)){Ex_slice = pinch(Ex,3,i);Ey_slice = pinch(Ey,3,i);TE_frac(i) = abs(Ex_slice)^2/(abs(Ex_slice)^2+abs(Ey_slice)^2);}# plot resultsplot(z*1e6,TE_frac,'z (um)','TE Fraction');
如果觉得《Lumerical官方案例 FDTD时域有限差分法仿真学习(五)——液晶(liquid crystal)》对你有帮助,请点赞、收藏,并留下你的观点哦!