1．介紹

動態(tài)多模分析的目的是進(jìn)行激光多模和激光調(diào)Q運(yùn)轉(zhuǎn)分析。激光腔內(nèi)橫模結(jié)構(gòu)近似為HG和LG模式。HG和LG模式是不同本征頻率對應(yīng)的正交特征函數(shù)，我們假設(shè)模式之間的橫模振蕩互不干擾，因此模式之間的短時(shí)干涉影響可以忽略。基于這個(gè)假設(shè)，起振模式中的反轉(zhuǎn)粒子數(shù)密度和光子數(shù)是由下面的以時(shí)間為變量的3D速率方程描述：

方程2-3中的參數(shù)如下

方程2-3用于四能級激光系統(tǒng)。準(zhǔn)三能級系統(tǒng)的多模分析還在研究中。三能級（泵浦能級）和能級2（激光上能級）之間，能級1（激光下能級）和能級（基態(tài)能級）之間的快速衰減速率已被假定。

在前面的章節(jié)中有詳細(xì)的關(guān)于計(jì)算激光輸出功率，Q開關(guān)運(yùn)轉(zhuǎn)和光闌影響的參數(shù)和方程之間的數(shù)學(xué)關(guān)系的介紹。或者，可以點(diǎn)擊LASCAD主窗口主菜單的“Help DMA Code”。

接下來有一個(gè)關(guān)于DMA編碼的指南。它顯示了怎樣比較合理地定義DMA GUI中單個(gè)輸入?yún)?shù)來模擬CW多模操作，Q開關(guān)運(yùn)轉(zhuǎn)和光闌的影響。

2．激光器連續(xù)輸出時(shí)輸出功率，模式競爭，和光束質(zhì)量的模擬

要使用DMA編碼，需要在腔內(nèi)插入一個(gè)熱透鏡晶體。可以參照教程1的說明，準(zhǔn)備一個(gè)端面泵浦的晶體。簡便地，我們可以直接打開tutorial-1.lcd文件激活腔結(jié)構(gòu)，該文件可以在LASCAD的子目錄“tutorials”中找到。在我們設(shè)置好FEA編碼和在模式腔內(nèi)插入晶體后，在LASCAD主菜單中選擇 “Dynamic Multimode Analysis”，打開DMA窗口。在該窗口中，點(diǎn)擊Open GUI for DMA，打開“Dynamic Multimode Analysis”。該窗口有5個(gè)標(biāo)簽，如下：

2.1 高斯模式選項(xiàng)

點(diǎn)選該選項(xiàng)并選擇“type of Gaussian modes”來近似激光模式結(jié)構(gòu)。如果模式結(jié)構(gòu)是象散的，選擇Hermite-Gaussian模式。在旋轉(zhuǎn)對稱腔結(jié)構(gòu)設(shè)計(jì)時(shí)，選擇LG模式更好。但是，如果橫模數(shù)量較多，則推薦采用Hermite-Gaussian模式，高階LG在數(shù)值上更復(fù)雜。模式類型選擇會自動地基于激光模式結(jié)構(gòu)的象散選擇。“Maximum transverse mode order”定義在x和y方向的最大傳輸模式階數(shù)N_max。但是，隨著N_max的增加，由于總模式數(shù)

因此M也會跟著變大，將耗費(fèi)一定的計(jì)算時(shí)間。

如果N_max已經(jīng)定義，我們需要設(shè)置“Number of grid points in x- and y- direction”足夠大來解決高階傳輸模式的傳輸強(qiáng)度振幅分布。否則，該數(shù)值和熱FEA中使用的橫向格點(diǎn)大小一樣。熱FEA大小仍然保留“Number of grid points in z-direction”。

“Stretch factor in x- and y-direction related to beam diameter”的定義嚴(yán)格依賴于N_max和在Help中描述的泵浦光分布情況。如果N_max=0，我們就需要將伸展因子設(shè)置為2來考慮基模和吸收泵浦光功率分布的疊加情況。

2.2 速率方程選項(xiàng)

雖然該標(biāo)簽的定義已經(jīng)在LASCAD的“Laser Power Output”窗口定義過，但是還是有必要在該標(biāo)簽中再定義一次。與時(shí)間無關(guān)的激光功率輸出結(jié)果相比，該窗口處的定義需要與別的輸出一致。

2.3 CW運(yùn)轉(zhuǎn)選項(xiàng)

“Time resolution”的默認(rèn)值10ns，與普通激光結(jié)構(gòu)相符。但是，該條目和“Time period used for simulation”應(yīng)該受計(jì)算結(jié)果控制。

2.4 光束質(zhì)量

光束質(zhì)量因子

和

是根據(jù)Siegman and Townsend用如下方程計(jì)算的

2.5 結(jié)果

DMA GUI中點(diǎn)擊“Calculate”，彈出DMA Calculation窗口顯示計(jì)算進(jìn)程。在窗口顯示的信息中，如“Power output for TEM00 mode using CW time independent recursion formula” 用于與相關(guān)計(jì)算的對比，在Laser Power Output窗口打開后開始計(jì)算。

在計(jì)算末，顯示了最后四分之一模擬時(shí)間的單個(gè)模式平均功率輸出。為了減少輸出尖峰的影響，平均時(shí)間必須嚴(yán)格與模擬時(shí)間的四分之一相同，如圖1。如果尖峰擴(kuò)展到最后四分之一模擬時(shí)間，我們推薦增加“CW Operation”標(biāo)簽中的“Time period of simulation”。接下來，如果沒有其他具體的參數(shù)指定，物理量按照上述方式進(jìn)行平均。

如果按照教程1中的條件設(shè)置腔結(jié)構(gòu)，并保持DMA GUI選項(xiàng)中的默認(rèn)設(shè)置（除了將Nmax設(shè)置為3），我們可以得到表1中的結(jié)果。

Mode (0,3): 1.23341

Mode (0,2): 0.598921

Mode (0,1): 0.560508

Mode (1,2): 0.550806

Mode (3,0): 0.545464

Mode (1,3): 0.491664

Mode (0,0): 0.456887

Mode (3,1): 0.417439

Mode (1,1): 0.413295

Mode (1,0): 0.346157

Mode (3,3): 0.323292

Mode (2,2): 0.321322

Mode (2,1): 0.305183

Mode (2,0): 0.160791

Mode (3,2): 0.149285

Mode (2,3): 7.35643e-040

表1 單個(gè)模式功率輸出教程1中腔結(jié)構(gòu)

該結(jié)果顯示教程1中象散對吸收泵浦功率密度的影響，在x-z平面呈高斯分布，在y-z平面是指數(shù)為10的超高斯分布。如表中所示，模式（m ,n）和模式（n, m）的功率大小不同，例如，Mode(0,3)功率為1.31W，Mode(3,0)的功率僅有0.57W。結(jié)果，x方向和y方向的光束質(zhì)量

，

不同。總輸出功率為6.87W。

在關(guān)閉DMA和LASCAD之后，我們也可以通過打開DMA相關(guān)目錄下的output.txt文件查看數(shù)值結(jié)果。

點(diǎn)擊DMA GUI的“Show Results”按鈕，打開DMA觀察器查看結(jié)果。在觀察窗口下方的下拉框中可以選擇重要的2D和3D圖。

圖1 輸出功率隨時(shí)間變化

圖1顯示，隨時(shí)間變化的輸出功率圖。計(jì)算開始時(shí)假定反轉(zhuǎn)粒子數(shù)密度

。因?yàn)檫@不同于熱平衡條件，所以一開始會看到尖峰，隨著時(shí)間的增加輸出功率逐漸變小，最后會趨向于常數(shù)。基于腔結(jié)構(gòu)，我們有必要增加“Time period used for simulation”來實(shí)現(xiàn)收斂。

單個(gè)模式隨時(shí)間變化的輸出功率也可以單獨(dú)顯示。

同時(shí)也可以顯示隨時(shí)間變化的輸出功率二維圖，沿腔軸TM00模式光斑圖。 “Beam profile in the crystal” 三維圖，基于單個(gè)模式對總輸出功率的貢獻(xiàn)，顯示了它們之間的疊加。基于教程1的腔結(jié)構(gòu)的例子，如圖2。該圖顯示了教程1中吸收泵浦光功率密度的象散分布影響。可見，強(qiáng)度分布有很明顯的象散，y方向的高階橫模比較突出。

圖2 晶體中光束截面分布

3 Q開關(guān)運(yùn)轉(zhuǎn)模擬

在“Pumping”下拉框中可以選擇兩個(gè)不同模式的泵浦，CW泵浦和脈沖泵浦。脈沖泵浦可用于近似被動Q開關(guān)模型，見第3.3節(jié)。

3.1 CW泵浦

對于連續(xù)波長泵浦，由重復(fù)固定頻率引發(fā)的預(yù)定義脈沖也可以計(jì)算。在速率方程中引入較高的人為腔損耗可以抑制裝入時(shí)的激光震蕩，可以在“Q-switch induced loss during load phase”框中定義。我們通常采用該參數(shù)的默認(rèn)值0.8。因?yàn)樵谳d入期沒有受激輻射發(fā)生，此時(shí)產(chǎn)生粒子束反轉(zhuǎn)。

如果打開周期大于0，人工調(diào)Q損耗不會立即減少，但是會在已定義的開放時(shí)間不斷地降低到普通腔損耗。但是，該參數(shù)對脈沖能量和形狀只有非常微小的影響。

因此，一個(gè)合適脈沖周期的定義很重要。該值不代表脈沖物理寬度，僅定義了脈沖計(jì)算的時(shí)域。在脈沖產(chǎn)生時(shí)，粒子束反轉(zhuǎn)和光子密度會瞬間變化，很有必要在脈沖周期內(nèi)定義一個(gè)數(shù)值較大的時(shí)步來獲得較好的離散。因?yàn)榍唤Y(jié)構(gòu)會導(dǎo)致脈沖生成被延遲，可能需要我們設(shè)置脈沖周期比脈沖帶寬大，來阻止脈沖展寬成弛豫周期。

加載時(shí)間“Load period”+打開時(shí)間“opening period”+脈沖時(shí)間“pulse period”必須小于脈沖重復(fù)周期“pulse repetition period”。在計(jì)算多個(gè)脈沖時(shí)，剩余時(shí)間就是松散周期，是脈沖周期和新的導(dǎo)入時(shí)間之間的緩沖區(qū)。在松散周期中的時(shí)步可以很小，因?yàn)榱Ｗ邮崔D(zhuǎn)和光子密度在此時(shí)幾乎不變化。

3.2 CW泵浦結(jié)果

再一次結(jié)合教程1中的腔結(jié)構(gòu)與DMA GUI標(biāo)簽中的默認(rèn)設(shè)置，除了設(shè)置基模（N_max=0），計(jì)算結(jié)果顯示了3個(gè)脈沖中的最后一個(gè)。

Power output averaged over pulse repetition period [W] = 2.92419 

Pulse energy [mJ] = 0.194946 

Pulse width (FWHM) [ns] = 5.15 

Average Beam Quality M² in x-direction = 1 

Average Beam Quality M² in y-direction = 1 

Peak power output 32216.1 [W] at time 0.000198894 [s]

輸出脈沖的最大值提供了光學(xué)元件損耗的重要信息。

像“DMA Viewer Help”描述的，在放大“Power output over time”的二維圖后，可見脈沖外形。圖3顯示了一個(gè)典型的脈沖形狀的例子。

圖3 CW泵浦典型脈沖外形

3.3 脈沖泵浦

除了“CW pumping”，“Pulsed synchronous with Q-switch”也可以選作“Mode of pumping”。該輸入可用于創(chuàng)建被動Q開關(guān)的近似模型。在脈沖泵浦的情況下，假設(shè)泵浦的起始時(shí)間與導(dǎo)入時(shí)間同步。在連續(xù)波泵浦中，脈沖泵浦頻率明顯要比脈沖重復(fù)頻率小很多，但是兩種脈沖的時(shí)間間隔很大。因此，脈沖之間沒有重疊發(fā)生，我們限制脈沖數(shù)為1。計(jì)算的時(shí)間周期以脈沖周期截止。

然而，脈沖重復(fù)頻率必須由以下因素來定義。目前的LASCAD版本使用一個(gè)時(shí)間獨(dú)立的代碼進(jìn)行熱分析。在熱分析時(shí)，如果采用脈沖泵浦，我們需要采用泵浦功率平均時(shí)間。如果泵浦頻率大于60HZ，近似精度就足夠了。在FEA窗口定義的值是用泵浦功率平均時(shí)間除以“Pulse repetition frequency”與“Load period”乘積，采用該計(jì)算可以得到脈沖泵浦的泵浦功率。因此，兩個(gè)條目都必須在選項(xiàng)卡Q-switch中定義。低泵浦頻率和單泵浦的與時(shí)間相關(guān)的熱FEA還在研究中。

因?yàn)槊}沖泵浦很容易實(shí)現(xiàn)很高的粒子束反轉(zhuǎn)，相比于連續(xù)波泵浦，輸出鏡的反射率低一些。我們繼續(xù)使用教程1中的腔面結(jié)構(gòu)，Nmax=0，速率方程標(biāo)簽下的輸出鏡的反射率為0.5，調(diào)Q條目的設(shè)置也一樣，得到如圖4的結(jié)果：

Pulse 0: M²x = 1 M²y =1 

Power output averaged over pulse repetition period [W] = 71.8257 

Pulse energy [mJ] = 14.3687 

Pulse width (FWHM) [ns] = 0.495 

Average Beam Quality M² in x-direction = 1 

Average Beam Quality M² in y-direction = 1 

Peak power output 2.3572e+007 [W] at time 0.000200002 [s] 

Power averaged over whole simulation time [W]: 71.8257, 

with whole simulation time [s]: 0.00020005. 

Power output of individual modes averaged over computation period [W]: 

Mode (0,0): 71.8257 

Power output summed up over all considered modes = 71.8257

圖4 脈沖調(diào)Q參數(shù)設(shè)置

圖5 脈沖泵浦脈沖外形

采用圖4的輸出參數(shù)，得到圖5所示的脈沖形狀。

4 光闌影響模擬

“Apertures”標(biāo)簽下的條目可以設(shè)置高斯光闌和硬邊光闌，另外，也可以是硬邊、高斯或者超高斯的平面鏡。

4.1 硬邊光闌和平面鏡

硬邊光闌或者平面鏡可以是橢圓的，矩形的。在“Type of Aperture”的下拉框中選擇相關(guān)條目。

如果孔徑的位置小于0，光闌代表了輸出平面鏡（即右端面鏡），在x，y方向的光闌大小不等。

采用教程1中的腔面結(jié)構(gòu)，N_max=3，半徑大小為0.2mm的圓形端面鏡，運(yùn)行連續(xù)波計(jì)算。表2顯示的結(jié)果與表1不同，因?yàn)檩敵銎矫骁R的半徑大小有限。只有模式(0,1)，（1,0）和（0,0）對輸出功率起作用。相應(yīng)地，與2.5節(jié)中沒有使用光闌的結(jié)果相比，輸出功率從6.9W減少到4W。但是，光束質(zhì)量得到改善，

。

Mode (0,1): 1.83459

Mode (1,0): 1.42093

Mode (0,0): 0.721171

Mode (1,1): 5.85802e-042

Mode (2,0): 4.07249e-042

Mode (0,2): 2.36935e-042

Mode (2,1): 9.62273e-043

Mode (1,2): 7.98643e-043

Mode (3,0): 6.69612e-043

Mode (0,3): 4.63281e-043

Mode (2,2): 3.62786e-043

Mode (3,1): 3.49377e-043

Mode (1,3): 2.9093e-043

Mode (3,2): 2.12117e-043

Mode (2,3): 2.08619e-043

Mode (3,3): 1.53504e-043

表2 最后四分之一模擬時(shí)間的單個(gè)模式平均功率輸出

采用教程1中的腔結(jié)構(gòu)，另外，輸出平面鏡的半徑大小為0.2m

4.2 超高斯輸出平面鏡

因?yàn)锳BCD矩陣算法中沒有考慮超高斯光闌，所以超高斯分布內(nèi)部就采用指數(shù)SG的高斯近似，如圖6所示，SG=4。該近似是合理的，因?yàn)橹挥邪霃叫∮?/p>

的模式才會對總的模式外形造成影響。獲得的高斯擬合用在ABCD矩陣算法中。

圖6 SG=4的超高斯反射率分布的高斯擬合

接著，我們繼續(xù)采用教程1的腔結(jié)構(gòu)，Nmax=3，其他與10.4.2.4，10.4.4.1節(jié)的設(shè)置相同。運(yùn)行連續(xù)波計(jì)算，我們得到表3的結(jié)果。我們可以看到，由于超高斯輸出鏡，僅有00模式存在。相應(yīng)地，光束截面分布顯示了著名的高斯基模外形。令人驚奇的是，該模式的功率為4.22W。該大小與2.4節(jié)得到結(jié)果相比，明顯高出了很多，4.2節(jié)中所有模式的總功率為6.87W。

與硬邊光闌的結(jié)果相比，硬邊光闌條件下的輸出功率為4W，而超高斯光闌的輸出功率為4.22W，然而，x，y方向的光束質(zhì)量改善到1。

Mode (0,0): 4.22301

Mode (0,1): 3.01251e-041

Mode (1,0): 1.1211e-041

Mode (0,2): 8.25512e-042

Mode (1,1): 6.72075e-042

Mode (0,3): 4.86947e-042

Mode (1,2): 4.73229e-042

Mode (2,0): 4.64064e-042

Mode (2,1): 4.03505e-042

Mode (1,3): 3.73545e-042

Mode (2,2): 3.60176e-042

Mode (2,3): 3.24972e-042

Mode (3,0): 3.24888e-042

Mode (3,1): 3.21808e-042

Mode (3,2): 3.16512e-042

Mode (3,3): 3.05088e-042

表3 最后四分之一受激時(shí)間的單個(gè)模式的功率輸出

基于超高斯輸出平面鏡，教程1的腔結(jié)構(gòu)的計(jì)算

5．結(jié)論

上面章節(jié)的所有計(jì)算證實(shí)了DMA編碼用于多模競爭的動態(tài)分析和Q開關(guān)都是一個(gè)非常好用的工具。該工具可以計(jì)算與時(shí)間相關(guān)一系列單個(gè)橫模的相互作用，光束質(zhì)量，激光輸出功率，橫模結(jié)構(gòu)的強(qiáng)度分布，連續(xù)波和脈沖泵浦的脈沖外形，脈沖波峰強(qiáng)度等。因此它給激光工程師在優(yōu)化激光腔結(jié)構(gòu)時(shí)提供了重要的信息和幫助。