本文介紹了模擬光在均勻介質(zhì)中傳播的四種快速而嚴格的方法。結果表明，在自由空間傳播中，對光滑強相位項的解析處理在減少計算量方面是非常有效的。因此，在不限制快速傅里葉變換算法應用的情況下，我們重新設計了平面波角譜（SPW）算子來(lái)處理線(xiàn)性、球形和一般光滑相位項。特別是對于非傍軸場(chǎng)傳播，所提出的技術(shù)可以顯著(zhù)地減少所需的采樣點(diǎn)數量。數值結果表明了新方法的有效性和準確性。
一.文章介紹
光學(xué)建模與設計是研究與開(kāi)發(fā)中極其重要的一部分。由于人們對高質(zhì)量光學(xué)系統（包括衍射光學(xué)和微光學(xué)、散射物體和部分相干源）的需求日益增加，基于幾何光學(xué)和物理光學(xué)相結合的模擬方法，即場(chǎng)追跡變得非常重要。這種模擬技術(shù)的一個(gè)重要部分是諧波場(chǎng)在均勻介質(zhì)中的傳播。然而，能夠快速、準確地模擬一般光場(chǎng)在自由空間中的傳播仍然是一項具有挑戰性的任務(wù)。常用的算法只能做到快速或者只是準確。
在本文中，我們沒(méi)有進(jìn)一步的物理近似，介紹了四種新的算法，基于平面波（SPW）算子的角譜，有效地計算包含平滑但強相位項的非傍軸矢量光場(chǎng)的傳播。根據光滑相位項的形狀，可以使用不同的傳播算子。它們的共同點(diǎn)是避免了光滑相位項指數函數的采樣。相反，平滑相位項是解析處理的，只需對殘差進(jìn)行采樣即可執行傳播操作；因此，稱(chēng)為半解析傳播技術(shù)。
首先，在第二節中我們給出一個(gè)問(wèn)題的描述并引入數學(xué)符號。然后，在第3節中，我們考慮了一個(gè)球面相位項，Mansuripur[6]為此引入了一種嚴格的技術(shù)，稱(chēng)為使用快速傅里葉變換（FFT）的擴展菲涅耳衍射積分。在本節中，通過(guò)應用Van der Avoort等人最初使用的數值合適的拋物線(xiàn)擬合技術(shù)改進(jìn)了該概念。在另一種情況下[7]，詳細討論了擴展菲涅耳算子在數值上可行的參數空間。此外，我們還介紹了擴展的菲涅耳算符的快速反演方法，用于快速計算非傍軸場(chǎng)到焦點(diǎn)區域的傳播。
在第四節中，我們描述了一個(gè)用于光場(chǎng)快速傳播的半解析SPW算子，它包含一個(gè)光滑的線(xiàn)性相位項。該方法基于線(xiàn)性相位項和橫向偏移量的解析處理。之后，我們將這兩種技術(shù)結合起來(lái)，得到了一個(gè)數值有效的半解析SPW算子，它能夠同時(shí)解析地處理線(xiàn)性和球形相位項。
最后，在第6節中，我們通過(guò)將光場(chǎng)分解成具有平滑線(xiàn)性相位項的子光場(chǎng)，將半解析SPW算子概念推廣到平滑相位的通用形狀。在目標平面上，所有傳播子光場(chǎng)被相干地相加，其中解析已知的平滑線(xiàn)性相位項以數值有效的方式使用第7節中介紹的逆拋物面分解技術(shù)（PDT）進(jìn)行處理。數值結果證明了新的傳播方法的有效性和準確性。所有的模擬都是用光學(xué)軟件VirtualLab完成的。
二.均勻介質(zhì)中的場(chǎng)追跡
在光場(chǎng)追跡法中，光在線(xiàn)性、均勻和各向同性介質(zhì)中快速而精確的傳播是由諧波場(chǎng)的概念處理的。結果表明，任何電磁場(chǎng)都可以分解為一組諧波場(chǎng)[8,9]。在空間頻率域中，以特定角頻率ω0振蕩的單次諧波場(chǎng)定義為
 (1)
用位置向量  和角頻率ω分別表示。請注意，下列理論是完全矢量的，因為在式（1）中，諧波場(chǎng)分量代表三個(gè)電場(chǎng)分量和三個(gè)磁場(chǎng)分量，由于計算效率高，常用的諧波傳播技術(shù)基于FFT算法[10]。一種嚴格的傳播技術(shù)是SPW算子[5]，其中各諧波場(chǎng)分量的復振幅在與傳播方向正交的平面邊界上，通過(guò)傅里葉變換（FT）分解成一組平面波
 (2)
 是初始平面邊界上的橫向位置向量，是  對應的空間頻率矢量。用  表示的平面波通過(guò)與傳播因子相乘，在距離z上傳播
 (3)
 表示折射率為n的均勻介質(zhì)中的波數，c為光的真空速度。最后，利用逆傅里葉變換將所有平面波疊加，從而得到SPW傳播算子，
 (4)
從物理角度來(lái)看，SPW算子對任何傳播距離z和任何空間頻率矢量k[5]都是有效的。然而，對于長(cháng)的傳播距離，采樣公式（4）的數值工作量太大。對于非傍軸場(chǎng)，它包含高頻分量，數值工作量將變得更高。圖1示意性地示出了由于快速增長(cháng)的數值工作量而導致的SPW算子的有限范圍。
一篇文獻中報道了兩個(gè)近似公式（4），以克服這一局限性。在這兩種情況下，即菲涅耳和遠場(chǎng)積分，使用近似來(lái)分離球面相位項與積分的數值計算，并且球面相位項是被解析地處理的。這就大大減少了數值計算的工作量。然而，由于這些近似，兩種解決方案的適用性受到限制。菲涅耳積分[11]使用空間頻率分量的泰勒展開(kāi)。用這種方法，將式（3）的球面相位函數替換為拋物線(xiàn)相位函數，從而得到式（4）中逆傅里葉變換積分的半解析解。如圖1所示，該概念僅適用于具有低空間頻率的傍軸場(chǎng)。對于長(cháng)傳播距離z，可以應用遠場(chǎng)積分[5]，其中公式（4）的逆傅里葉變換積分可以用固定相位的想法來(lái)半解析地求解。
圖1：場(chǎng)追跡中常用的基于FFT的傳播算子的有效范圍的示意圖�；疑珔^域表示在數值上可行的參數空間。嚴格的SPW算子在長(cháng)傳播距離或大空間頻率下很難突破數值工作量。由于物理近似的性質(zhì)，遠場(chǎng)和菲涅耳傳播算子具有有限的有效范圍。
由于普通的自由空間傳播算子使用范圍的限制，先進(jìn)的傳播技術(shù)的發(fā)展對于非傍軸場(chǎng)追跡是必不可少的。這些技術(shù)應覆蓋圖1中剩余的白參數空間，并在精度和數值計算之間取得最佳折衷。文獻[2,3,12-14]中有幾種方法。然而，這個(gè)問(wèn)題的一般解決辦法仍不清楚。在實(shí)踐中，非傍軸諧波場(chǎng)振幅出現了采樣問(wèn)題，其中部分包含強而光滑的相位項  數學(xué)表達式為
 (5)
|
