用于模拟波光学元件的非参节性高斯光束公式

2018年6月26日

在以前的博客文章中,我们讨论了近去高斯梁公式。今天,我们将讨论从Comsol®软件5.3A版本可用的高斯光束的更准确的配方。基于平面波膨胀的该配方比传统的携带式配方更准确地处理非帕型高斯光束。

高斯光束

众所周知高斯梁公式仅适用于近后高斯梁。近亲意味着光束主要沿光轴传播。有几篇论文在定量意义上谈论了近亲性(请参阅参考。1)。

粗略地说,如果梁腰大小接近波长,则梁以更高的角度传播到焦点。因此,近端假设分解,配方不再准确。为了减轻此问题,为您提供更通用和准确的高斯梁的配方,我们引入了非副束束的配方。在用户界面中,这称为飞机波膨胀

该方法基于平面波的角光谱((参考。2),有时被称为角光谱法((参考。3)。

平面波的角光谱

让我们简要介绍2D中的近去式高斯光束公式(为了更好的视觉效果和理解)。

我们从麦克斯韦(Maxwell)的方程式开始,假设有时间谐波场,从中,我们从中获得了以下helmholtz的方程式,用于波长的平面电场。\ lambda对于我们选择极化:

\ left(\ frac {\ partial^2} {\ partial x^2} + \ frac {\ partial^2} {\ partial y^2} + k^2 \ right)e_z = 0,)e_z = 0,

在哪里k = 2 \ pi/\ lambda

平面波的角光谱基于以下简单事实:随意的满足上述Helmholtz方程的场可以表示为以下平面波膨胀:

e_z(x,y)= \ int_ {k_x^2 +k_y^2 = k^2} a(k_y)e^{i(k_x x +k_y y)} dk_y,

在哪里一个(k_y)是一个任意功能。

整合路径是半径圆k真正的k_xk_y。(用于复杂k_xk_y,集成域延伸至复杂的平面。)功能一个(k_y)称为角谱函数。一个人可以证明这一点E_Z通过直接替换来满足Helmholtz的方程式。

现在,我们知道这种公式总是为Helmholtz的方程式提供精确的解决方案,让我们尝试以视觉上理解它。从约束,k_x^2+k_y^2 = k^2,我们可以设置k_x = k cos(\ varphi)k_y = k sin(\ varphi)并将上述方程式重写为:

e_z(x,y)= \ int _ { - \ pi/2}^{\ pi/2} a(\ varphi)e^{ik(x \ cos \ cos \ cos \ varphi +y \ sin \ sin \ varphi)}

上述公式的含义是,它构造波浪作为总和或整体,由许多在各个方向传播的波组成,所有波浪都具有相同的波数k。这在下图中显示。

平面波的角光谱的插图。
平面波角光谱的可视化。

回到k_y表示,当使用此公式实际解决问题时,您要做的就是找到角光谱函数一个(k_y)满足边界条件。通过假设横向场的轮廓(垂直于传播方向,即光轴)也是高斯形状(请参阅参考。4),一个人可以得出a(k_y)= \ exp(-k_y^2 / w_0^2), 在哪里W_0是频谱宽度。

通过更多的数学操纵,我们在光谱宽度之间得到了关系W_0和梁腰半径。例如,对于缓慢的高斯梁,角光谱是狭窄的。另一方面,平面波是角光谱函数是Delta函数的极端情况。对于快速高斯梁,角光谱更宽,反之亦然。

这是对非参节性高斯梁的基本理论的快速摘要。要回顾我们到目前为止所显示的内容,让我们通过使用极性坐标再次重写公式,x = r \ cos \ theta,\ y = r \ sin \ theta

e_z(r,\ theta)= \ int _ { - \ pi/2}^{\ pi/2} a(\ varphi)e^{ikr \ cos(\ theta- \ varphi)} d \ varphi。

这是出生和狼的表述(参考。2)在他们的书中使用。

3D公式更为复杂,并且由于极化而看起来不同,但是基本思想与上述参考文献中相同。根据您是否认为逃生的波浪,它看起来也不同。平面波扩展方法波光学模块RF模块,尽管基于角谱理论,但适用于数值计算。

平面波扩展:设置和结果

让我们比较新功能,飞机波膨胀,具有先前可用的功能,近似近似。这设置覆盖这两种方法的窗口如下所示。

电磁波的屏幕截图,comsol多物理学中的频域设置。
飞机波膨胀功能设置。

有了新功能,如果您有两个选择自动的设置并不能使您令人满意的近似值:

  1. 波矢量计数
  2. 最大横波数

第一个选项确定离散水平的数量,具体取决于您要表示高斯光束的罚款。飞机波越多,得到的越好。第二个选项与上一个方程中的积分界有关。IE。,- \ pi/2 \ le \ varphi \ le \ pi/2。这种整体界限可以是最大的\ pi/2对于最小的斑点大小,对于较慢的光束,可能会更浅,具体取决于高斯梁的速度。您需要更多具有较大横向波数的倾斜平面波来代表更快(更集中)的光束。

以下结果比较了斑点半径为的情况的两个公式\ lambda/2,这是相当非参与性的。与上一篇博客文章一样,模拟是通过分散的场公式和域被完美匹配的层(PML)包围。这样,散射场表示来自精确的Helmholtz解决方案的误差。

下面的左图显示了新功能,而右侧的图像显示了近似近似。顶部图像显示了计算的高斯光束背景字段的规范,ewfd.ebz,虽然底部图像显示散射的场范围,但ewfd.Relez,代表确切的Helmholtz解决方案的误差。显然,在非选项方法中,Helmholtz解决方案的误差大大减少。

波光学仿真结果显示了计算的高斯光束背景场和散射场标准的规范。
平面波的角光谱与近似公式之间的比较。

总结说

我们已经讨论了使用新的平面波膨胀选项的非允许高斯梁的近似方法的理论和结果。请记住,该公式非常准确,但仍然是假设下的近似值。首先,我们对焦平面中的场形状做出了假设。其次,我们假设逃生场为零。如果您对字段耦合到一些感兴趣纳米结构在快速高斯梁中,在焦点区域附近,您可能需要计算evanscent场。您可以在以前的博客文章中了解更多有关非参与性高斯光束的逃生成分的信息:非参与式高斯光束的渐渐成分

下一步

通过单击下面的按钮:

注意:此功能也可以在RF模块

参考

  1. P. vaveliuk,“激光束中近似近似的极限”,光学字母,卷。32,不。8,2007。
  2. M. Born和E. Wolf,光学原理,ed。7,剑桥大学出版社,1999年。
  3. J. W. Goodman,傅立叶光学元件
  4. G. P. Agrawal和M. Lax,“超出近似近似之外的自由空间波传播”,物理。修订版A。27,第1693–1695页,1983年。

评论(18)

发表评论
乐动体育app
加载...
穆罕默德·穆罕默德(Mohammed Mohammed)
穆罕默德·穆罕默德(Mohammed Mohammed)
2018年7月17日

亲爱的米苏山
非常感兴趣的话题。我有一个问题。
为什么lambda等于500nm并在comsol中用作频率计算的默认值(f = c_const/500 [nm])?

瑞安·弗里曼(Ryan Freeman)
瑞安·弗里曼(Ryan Freeman)
2019年2月25日

是否有人能够使用此方法来创建专注于衍射极限的高斯光束?

我目前正在尝试使用设置来创建一个紧密的横梁,但是所得梁的宽度约为5-10倍,太大了。

是否还有有关此方法或新示例的文档?

谢谢,
-ryan

Mizuyama Yosuke
Mizuyama Yosuke
2019年2月25日

亲爱的瑞安,

感谢您阅读我的博客。如果您没有得到专注的光束,请尝试增加我在帖子中写的横向波数。否则,请将您的模型发送给support@comsol.com
谢谢你!

此致,
Yosuke

迈克尔·托马斯
迈克尔·托马斯
2019年2月25日

你好Yosuke,

我有同样的问题。

我假设2Pi/lambda是支持的最大横向波数,对应于最大焦点。

问题可能会试图将光束聚焦到薄膜(30nm)上吗?

出色的博客文章Yosuke。为新功能感到兴奋。

-Michael

Mizuyama Yosuke
Mizuyama Yosuke
2019年2月25日

亲爱的迈克尔,
感谢您的评论。
请向我们发送您的文件,以便我们发现问题。
同时,我将研究功能。
此致,
Yosuke

Mizuyama Yosuke
Mizuyama Yosuke
2019年2月27日

亲爱的瑞安和迈克尔,
我模拟了半lambda腰。我认为在这里展示我对一个参数案例的观察是一件好事。这几乎接近衍射有限的光束尺寸。对于如此高的非律性,您需要修改默认设置。如果将“ Wave Vector分布类型”更改为“用户定义”,则将看到两个参数“ Wave Vector Count”,NK和“最大横向波数”,Kt,Max。请阅读博客内容以获取更多详细信息。我选择了nk = 401和kt,max = ewfd.k0/2用于获得平滑结果。对于衍射有限的光束,最大横向波数不一定是最高的,即EWFD.K0,因为从如此高角度的贡献也很小,即使对于接近衍射有限的梁。如果您选择KT,则最大为ewfd.k0,则NK还不够。您需要增加NK以保持相同的离散水平。
到目前为止,您需要这种考虑到非常高的非律束。
我希望这些信息对你有用。
此致,
Yosuke

尤里安格郑
尤里安格郑
2020年1月5日

尊敬的Yosuke Mizuyama,
我想让高斯光束在金属膜上入射,并使高斯束的腰部落在金属膜上,我应该如何设置它?高斯梁是从空气介质到金属介质表面的。I can’t use your method to achieve the purpose, I want to consult you.Do I need to set the incident electromagnetic field formula?Can you give me a reference formula?
如果您有时间,请帮助解决此问题。
谢谢,
-cheng

Mizuyama Yosuke
Mizuyama Yosuke
2020年1月6日 Comsol员工

亲爱的尤利安,

感谢您阅读我的博客。
我认为您可以使用此博客文章(第二个方程式)中给出的公式。
将金属表面放在x = 0。
将激发设置在左输入边界,例如x = -x0。
通过将x = -x0替换为y的函数,根据公式计算激发字段。
您应该能够激发一个非顺式高斯光束边界场,并将其聚焦在x = 0处的金属表面上。
但是,几天后,请再次阅读此博客。我在第二个公式中发现了一个错误。
我现在会纠正它。更新页面需要几天的时间。
感谢您让我再次阅读并找到错误!

此致,
Yosuke

Mizuyama Yosuke
Mizuyama Yosuke
2020年1月6日 Comsol员工

亲爱的尤利安,

现在进行了更正,并且页面已经实时。
角光谱函数A(ky)是EXP(-KY^2/W0^2),其中W0是一个代表梁的速度的常数。
祝你好运!

此致,

Yosuke

凯凯
凯凯
2020年8月6日

尊敬的Yosuke Mizuyama,
我想让高斯光束倾斜地出现在光栅上,并使高斯束的腰部落在光栅上,我该如何设置?

妙妙
妙妙
2020年4月24日

嗨,Yosuke博士
很抱歉再次打扰您。
我使用“平面波膨胀”方法设置了紧密浓缩的高斯光束,但是当我将波矢量计数增加到100时,发生了错误,日志文件显示“ java.lang.stackoverflowerror”。我对此一无所知,你能帮我吗?
最好的祝愿
真挚地,
Miao Peng

Mizuyama Yosuke
Mizuyama Yosuke
2020年4月24日 Comsol员工

嗨,Miao,
那不会发生在我身上。
请将模型文件发送到support@comsol.com
此致,
Yosuke

Mohammad Hosein Khosravi
Mohammad Hosein Khosravi
2020年5月26日

尊敬的Yosuke Mizuyama,

感谢您有用和有用的文章。
是否可以定义聚焦梁从此开始的端口?
当我使用这种方法模拟水中高度聚焦的高斯光束的传播时,我通过增加传播长度,焦点尺寸开始增加而发现,为什么会发生这种情况?我认为默认情况下,光束从真空环境开始,并进行了定义的培养基,对吗?
我希望波从玻璃培养基传播到水中培养基。如何改变波浪的初始介质?
真挚地,
穆罕默德·霍森(Mohammad Hosein)

Mizuyama Yosuke
Mizuyama Yosuke
2020年5月27日 Comsol员工

亲爱的穆罕默德,
感谢您阅读我的博客。
该端口功能应用于横向本本元素字段。
因此,高斯光束不是港口的合适梁。
材料特性包含在波矢量中,因为k = n*k0在公式中。
此致,
Yosuke
相反,您可以在散射边界条件或匹配的边界条件下作为用户定义的字段实现此目标。

Mohammad Hosein Khosravi
Mohammad Hosein Khosravi
2020年6月9日

感谢您的回复,亲爱的Yosuke Mizuyama,

您知道我如何在客观镜头的光束半径和Na之间建立联系吗?

真挚地,
穆罕默德

Mizuyama Yosuke
Mizuyama Yosuke
2020年6月10日 Comsol员工

亲爱的穆罕默德,
对于近距离高斯梁,雷利范围外的差异角约为theta = lambda/(n pi w0)。Na的定义是Na = N Sin(Theta)〜n Theta = lambda /(pi w0)。这是近期案件的关系。

Tau Rasethuntsa
Tau Rasethuntsa
2021年5月14日

Zubair Buch
Zubair Buch
2021年5月30日

我想知道是否有一种方法可以计算高斯光束的光功率?此外,我们如何才能改变我们在非掌式公式中产生的光束中插入的光功能的值?试图改变电子场(v/m)的价值似乎在这里没有任何功能,而不是飞机波,对吗?
谢谢,祖拜尔

探索comsol乐动体育赛事播报博客