双层扭曲光子晶体辐射场拓扑光力； - 复现：2022子刊Science Advances； - 结果关键词：光子晶体，能带，远场偏振椭圆分布 (偏振场)，连续域束缚态 (BIC)，一幅图绘制品质因子和能带的3D图，光力的麦克斯韦张量面积分求法，坡印廷矢量 - 软件：comsol，matlab (也可以不用，仅comsol内出图) - 备注：这篇大文章光力的结果建立在能带和远场偏振基础上，因此价格上考虑累积。
最近在研究领域，2022 年子刊 Science Advances 上关于双层扭曲光子晶体辐射场拓扑光力的研究引起了广泛关注，今天就来和大家聊聊如何复现这一有趣的研究成果。
一、研究关键要素
1. 光子晶体与能带
光子晶体是一种人造的周期性结构材料，它能像半导体控制电子运动一样，控制光子的传播。在这个研究里，光子晶体的能带结构是关键。能带描述了光子在晶体中允许存在的能量范围。通过 Comsol 软件，我们可以对光子晶体进行建模，比如创建一个简单的二维光子晶体结构，以下是 Comsol 里的部分建模代码示例（伪代码）：
% 创建二维几何模型 model = createpde('electromagnetic','spherical'); geometryFromEdges(model,phc_geometry); % phc_geometry 为定义光子晶体几何形状的函数 % 设置材料属性 setmaterial(model,'all','ElectromagneticMaterial','Silver');在上述代码中，我们首先创建了一个电磁学相关的 PDE 模型，然后通过geometryFromEdges函数定义光子晶体的几何形状，最后设置材料为银。这只是简单示意，实际建模会复杂很多，不同的光子晶体结构和材料属性对能带影响巨大。
2. 远场偏振椭圆分布（偏振场）
远场偏振椭圆分布描绘了光在传播到远处时电场矢量的振动方向和幅度变化情况。它与光子晶体的结构和能带相互关联。在 Comsol 里模拟远场偏振时，我们需要设置合适的边界条件和场监视器来观察远场的偏振特性。例如：
% 设置远场边界条件 addboundary(model,'farfield','ElectromagneticFarField','PEC'); % 定义场监视器观察远场偏振 fem = createpde('electromagnetic','spherical'); addboundary(fem,'farfield','ElectromagneticFarField','PEC'); monitor = addmonitor(fem,'PolarizationMonitor');通过这些设置，我们能在模拟过程中捕捉到远场的偏振信息，进而分析其椭圆分布特征。
3. 连续域束缚态（BIC）
连续域束缚态是一种独特的光学现象，在这种状态下，光被限制在一个开放的连续谱区域而不发生辐射损耗。在光子晶体中，BIC 的存在对光力等特性有重要影响。在研究中，我们需要通过调整光子晶体的结构参数，如晶格常数、柱体半径等来寻找 BIC 状态。在 Comsol 模拟里，这意味着不断尝试不同的几何参数输入，观察能带结构和场分布的变化，以确定 BIC 出现的条件。
4. 绘制品质因子和能带的 3D 图
品质因子衡量了光学谐振腔中存储能量与耗散能量的比率，与 BIC 密切相关。我们可以通过 Comsol 模拟得到的数据，在 Matlab 中绘制品质因子和能带的 3D 图（当然也可以仅在 Comsol 内出图）。以下是 Matlab 绘制 3D 图的简单代码示例：
% 假设 data 是从 Comsol 导出的包含品质因子、能带和一个变量的数据 x = data(:,1); % 变量 y = data(:,2); % 能带 z = data(:,3); % 品质因子 figure; scatter3(x,y,z); xlabel('Variable'); ylabel('Band Structure'); zlabel('Quality Factor'); title('3D Plot of Quality Factor and Band Structure');这段代码将从 Comsol 获得的数据进行 3D 散点图绘制，直观展示品质因子和能带之间的关系。
5. 光力的麦克斯韦张量面积分求法与坡印廷矢量
光力的计算在这个研究里至关重要，这里采用麦克斯韦张量面积分的方法。麦克斯韦张量包含了电磁场的应力信息，通过对其在物体表面进行面积分可以得到光力。坡印廷矢量则描述了电磁能量的流动方向和大小，它在光力计算中也起到关键作用。在 Comsol 中，我们可以利用其内置的物理场接口和后处理功能来计算麦克斯韦张量和坡印廷矢量。例如：
% 获取麦克斯韦张量 Maxwell_tensor = evalexpression(model,'MaxwellTensor'); % 获取坡印廷矢量 Poynting_vector = evalexpression(model,'PoyntingVector');通过这些表达式，我们能得到所需的物理量，进而完成光力的计算。
二、软件选择与协作
本次研究主要用到 Comsol 和 Matlab，当然仅使用 Comsol 内出图也是可行的。Comsol 强大的多物理场建模能力可以精确模拟光子晶体的各种物理特性，而 Matlab 在数据处理和绘图方面有着出色的表现。两者结合能更高效地完成从建模到结果呈现的整个流程。
三、总结与价格考虑
这篇大文章中光力的结果建立在能带和远场偏振基础上，整个研究过程涉及多个复杂物理量的计算和分析，在价格上考虑累积也是合理的。通过一步步复现这些研究内容，不仅能深入理解双层扭曲光子晶体辐射场拓扑光力的物理机制，也能提升我们在光子学领域的研究和模拟能力。希望感兴趣的朋友可以一起探索，在这个有趣的研究方向上共同进步。
