紧聚焦光场中的磁场分量如何计算？

聚焦是光学中比较常见的现象，当用来聚焦的透镜具有较高的数值孔径时，称为紧聚焦​。与弱聚焦相比，紧聚焦存在一些特殊的现象，如纵向场的产生​。关于紧聚焦场在matlab中的计算之前已经介绍过很多（参考：论文复现--紧聚焦场的横向自旋，论文复现-双共轴纵向偏振涡旋，论文复现及其扩展-焦点的三维移动及任意轨迹拼接），然而，我们通常仅仅关注其电矢量分布，很少关注磁矢量分布​。但实际上，光场是电磁场，故而同时存在电矢量和磁矢量。此外，在很多情况下也需要知道磁矢量，如计算能流、光场手性等​。

下面我们来介绍如何基于matlab计算紧聚焦条件下的磁矢量分布​。为了说明结果的正确性，我们选择一篇论文进行对比​。论文信息如下​：

这篇论文正好是计算光场手性，但这不是本期介绍的重点，这里我们主要关注紧聚焦条件下的磁矢量分布​。先来看原文中计算的电矢量和磁矢量分布

原文计算了两种情况，即入射场分别具有-1和+1的轨道角动量，偏振为x方向线偏振，聚焦的数值孔径为0.9​。这里我们首先给出计算紧聚焦场电矢量的表达式为​：

其中和分别表示入射场的径向和旋向分量，任何偏振态都可以通过坐标变化转化为径向和旋向分量的和​。例如对于x方向线偏振，有

对于磁矢量，由于其与电矢量和波矢相互垂直，则有

其中

然后用替换积分中的，则得到紧聚焦下的磁矢量分布​。以拓扑荷为-1时为例（对应原文图左边），电矢量和磁矢量分布为

图1 电矢量分布​。

图2 磁矢量分布​。

可以看出，分布与原文基本一致，但存在两处差别​。第一，我们计算出的分布需要逆时针旋转180度才能与原文一致，这是由于旋向坐标与原文不一致，从原文图(a)的入射场相位分布可以看出，其旋向坐标相对正常坐标本身就沿逆时针旋转了180度，故而我们的计算结果实际上是与​原图一致的。第二，磁场y分量的相位分布有少许差别，这可能是由入射场E0的取值所导致的，但这并不影响结果的正确性​。

详情可关注公众号“光学师兄的日常”。