-
空气中的各种悬浮颗粒附着在镜面上形成表面污染物,这些颗粒污染物可近似视作球形粒子,进而可以用米氏理论来对镜面污染物的吸收特性进行定量分析。污染粒子的吸收主要来源为光照射在粒子上引起粒子内部电场的调制。如图1所示,单色平面光照射在颗粒上,入射光除了发生散射外,同时也会有部分能量在颗粒内部形成内电场,从而产生吸收特性。
由Mie理论[13]可知,单位振幅的入射光场照射颗粒体内后形成的内电场Eint平方的期望值可表达为:
$$\left\langle {{{\left| {{E_{{\rm{int}}}}} \right|}^2}} \right\rangle = \frac{1}{4}\sum\limits_{n = 1}^\infty {\left( {{m_{{n}}}{{\left| {{c_{{n}}}} \right|}^2} + {n_{{n}}}{{\left| {{d_{{n}}}} \right|}^2}} \right)} $$ (1) 式中:参量cn、dn为Mie吸收系数;mn和nn为与散射角函数相关的量,可分别表示为:
$$ {c_{{n}}}{\rm{ = }}\frac{{{j_{{n}}}\left( x \right)[xh_{{n}}^{(1)}(x)]' - h_{m{n}}^{(1)}(x)[x{j_{{n}}}(mx)]'}}{{{j_{{n}}}(mx)[xh_{{n}}^{(1)}(x)]' - {h_{{n}}}^{(1)}(x)[mx{j_{{n}}}(mx)]'}} $$ (2) $$ {d_{{n}}} = \frac{{m{j_{{n}}}\left( x \right)[xh_{{n}}^{(1)}(x)]' - mh_{{n}}^{(1)}(x)[x{j_{{n}}}(mx)]'}}{{{m^2}{j_{{n}}}(mx)[xh_{{n}}^{(1)}(x)]' - {h_{{n}}}^{(1)}(x)[mx{j_{{n}}}(mx)]'}} $$ (3) $${m_{{n}}} = 2\left( {2n + 1} \right){\left| {{j_{{n}}}\left( z \right)} \right|^2}$$ (4) $${n_{{n}}} = 2n\left( {2n + 1} \right)\left. {\left\{ {\left( {n + 1} \right){{\left| {\frac{{{j_{{n}}}\left( z \right)}}{z}} \right|}^2} + {{\left| {\frac{{\left( {z{j_{{n}}}\left( z \right)} \right)}}{z}} \right|}^2}} \right.} \right\}$$ (5) 式中:x=ka/2用于表征颗粒的尺寸,为无量纲量,a为球形颗粒直径;m为相对折射率;z=x或z=mx。
利用公式(1)可计算单个颗粒在单位振幅下吸收值随粒子尺寸变化。图2给出了不同尺寸的污染颗粒内部的电场强度随粒子径向的变化。从图2可看出,粒子内部的电场强度随着颗粒径向尺寸的增大而减小,且内电场在颗粒中心处达到最大,而在边缘处最小;在材料折射率不变的情况下,不同尺寸颗粒内电场分布随径向尺寸的增大而明显减小,且在靠近边缘处的电场强度变化幅度会变小,并逐渐趋于平缓。
图 2 不同尺寸污染颗粒的内电场强度
Figure 2. Internal electric field strength of contaminated particles with different sizes
进一步对公式(1)做球面积分,则可获得单个颗粒的吸收截面为:
$${\sigma _{{\rm{abs}}}} = k\varepsilon ''{\int_V {\left| {{{{E}}_{{\rm{int}}}}} \right|} ^2}{\rm{d}}V$$ (6) 式中:k=2π/λ为入射光波数,λ为入射光波长;ε″为该颗粒的相对介电常数的虚部(颗粒介质相对环境介质)。
由参考文献[14]可知,吸收截面与吸收的关系为:Qabs=4σabs/(πa2),进而可得到单个颗粒对光的吸收为:
$${Q_{{\rm{abs}}}} = \frac{{4\text{π} \varepsilon ''}}{{{x^2}}}\int_0^x {\left\langle {{{\left| {{{{E}}_{{\rm{int}}}}} \right|}^2}} \right\rangle } {x'^2}{\rm{d}}x'$$ (7) 图3给出了在不同折射率实部及虚部情况下,污染颗粒吸收值随粒子尺寸的变化,如图3所示。
图 3 不同折射率(a)实部和(b)虚部污染颗粒的吸收
Figure 3. Absorption of polluted particles with different refractive index (a) real part and (b) virtual part
从图3(a)可以看出,随着折射率实部的增大,颗粒吸收曲线第一个峰值向左移动,且第一个峰值不断升高,曲线峰值的变化越来越大;从图3(b)可以看出,折射率虚部的变化只会导致曲线峰值的变化。并且,随着虚部的增大,其峰值逐渐减小,在靠近粒子大尺寸处,吸收变化逐渐趋于平缓,几乎不随粒子折射率虚部的变化而变化。
-
在实际应用中,通常使用表面洁净度来表示光学元件的污染情况。根据美国军用标准MIL-STD-209D,镜面的污染程度可用表面洁净度来表征。表面洁净度用一个数值描述,记作X1,表示每平方英尺镜面上大于该数值的颗粒数目刚好为1。该标准表明,镜面上粒子直径与其数量分布近似满足:
$$\log {N_{\rm{a}}} = C \cdot \left[ {{{\left( {\log {X_1}} \right)}^2} - {{\left( {\log X} \right)}^2}} \right]$$ (8) 式中:X为颗粒直径;Na为元件表面每平方英尺上颗粒直径大于等于X的粒子数量(1 ft=0.305 m);C是归一化常量,对于清洁表面,C的值近似为0.926,对于未清洁表面,C的值近似为0.383。系数C对于清洁过和未清洁情况存在差异的主要原因在于对于未清洁的表面,系数0.383可能与观察结果更好地吻合,也就是说,来自空气的颗粒沉降产生了朝向大颗粒的尺寸分布,而表面清洁作用则在去除大颗粒方面比小颗粒更有效。因此,当使用MIL-STD指定清洁度水平时,将其限制在暴露于沉降物后已清洁的表面是很重要的。
图4给出了不同洁净度光学镜面的污染颗粒分布。分析图4可知,污染颗粒的尺寸越大,附着在镜面上的数量就越少。进一步地,C系数的选取对粒子数目影响也较大,在相同表面洁净度时,清洁和未清洁镜面的小尺寸颗粒数目大致相差2个数量级。然而,随着颗粒直径增大,清洁与未清洁的镜面污染颗粒之间的差距逐渐减小。
图 4 不同洁净度光学镜面的污染颗粒分布
Figure 4. Distribution of polluted particles on optical mirrors with different cleanliness
根据公式(7)和(8),镜面的颗粒污染物吸收Aabs可写为:
$${A_{{\rm{abs}}}} = \sum\limits_{X = 5}^\infty {[{N_{\rm{a}}}(X) - {N_{\rm{a}}}(X + 1)]} {Q_{{\rm{abs}}}}(X)$$ (9) 式中:X为大于5的整数。将镜面各尺寸污染颗粒吸收值与镜面吸收总量的比值记作Ra(n),即:
$${R_{\rm{a}}}(n) = \frac{{\left[ {{N_{\rm{a}}}(n) - {N_{\rm{a}}}(n + 1)} \right]{Q_{{\rm{abs}}}}(n)}}{{\sum\limits_{n = 1}^{{X_1}} {\left[ {{N_{\rm{a}}}(n) - {N_{\rm{a}}}(n + 1)} \right]{Q_{{\rm{abs}}}}(n)} }}$$ (10) 式中:Na(n)–Na(n+1)表示在n~(n+1) μm尺寸范围内镜面污染颗粒的数量。图5给出了不同洁净度镜面的Ra值随颗粒尺寸的变化曲线。由图5可知,随着颗粒尺寸的增大,其吸收系数呈下降趋势,经归一化处理后,其变化趋势不随表面洁净度等级的变化而变化。
图 5 不同洁净度时镜面各尺寸颗粒吸收值与颗粒吸收总值的比值
Figure 5. Ratio of absorptive value of particles of different sizes to total absorptive value of particles in different cleanliness
图6和图7分别给出了不同入射波长和灰尘不同折射率时,镜面不同尺寸灰尘的光吸收与镜面灰尘吸收总量的比值。其中,镜面洁净度X1=300,常数C=0.926。不同入射波长下镜面灰尘的折射率,如表1[15]所示。
由图6、7可知,镜面颗粒吸收主要集中在小尺寸的颗粒,随着颗粒尺寸增大,其吸收量呈下降趋势。
进一步从图7可以明显看出,当入射波长和表面洁净度一定的情况下,元件表面吸收几乎只与折射率虚部有关,虚部越大尺寸对吸收的影响越大,吸收比值变化梯度也越大。
图 6 不同波长下镜面各尺寸灰尘的吸收值与灰尘收总量的比值
Figure 6. Ratio of absorptive value of specular dust to total dust collection at different wavelengths
图 7 不同折射率下镜面各尺寸污染颗粒的吸收与镜面颗粒吸收总量的比值
Figure 7. Ratio of polluted particle absorption to total absorption of specular particle with different refractive index
表 1 不同波长下的灰尘折射率
Table 1. Refractive index of dust at different wavelengths
Item Value λ/μm 0.325 0.633 1.15 3.39 10.6 n 1.53+0.01i 1.53+0.000 5i 1.5+0.001i 1.5+0.02i 1.7+0.2i -
空气洁净度是影响镜面洁净度的重要因素之一,为此需要明确不同空气洁净等级下,污染颗粒附着在镜面上的规律。美国联邦标准(FED-STD-209D)规定了空气洁净度等级标准:每立方英尺或每立方米空气中,直径大于或等于0.5 μm粒子的数目。空气中对的颗粒尺寸分布满足如下关系:
$$\log {N_{\rm{c}}} = - 2.197\log X - 0.662 + \log {X_{\rm{c}}}$$ (11) 式中:X为颗粒的直径;Xc是空气洁净度等级;Nc是每立方英尺中颗粒直径大于等于X的粒子数量。基于大量关于粒子沉降速率的实验,Hamberg在文献中总结得出[16]:当洁净室的空气洁净度已知时,粒子在空气中的沉降规律可表示为:
$$\frac{N}{t} = \rho N_{\rm{c}}^{0.773}$$ (12) 式中:N为每平方英尺表面上直径大于5 μm的颗粒数量;t为元件曝露在洁净室空气中的时间,单位为天;ρ为洁净室相关因数,一般情况下ρ=2 851。结合公式(8)、(11)和(12)可得到镜面洁净度Xl与曝露时间t和空气洁净度Xc之间的关系,即:
$$\log {X_{\rm{l}}} = {\left[ {1.08\left( {\log \rho + \log t + 0.773\log {X_{\rm{c}}} - 1.248} \right)} \right]^{1/2}}$$ (13) 由公式(13)可知,当元件工作面朝向和曝露在空气中的时间给定时,元件表面洁净度随着空气洁净度的增大而增大。
将公式(7)与(13)结合,可计算空气洁净度对镜面颗粒污染物吸收的影响,如图8所示(图中1 ppm=10−6)。
图 8 空气洁净度对镜面颗粒污染物吸收的影响
Figure 8. Effects of air cleanliness on specular particulate pollutant absorption
由图8可知,当曝露时间一定时,镜面吸收会随着空气洁净度等级的增大而增大,并且对元件表面吸收产生影响的主要因素是其折射率的虚部。
需要指出的是,Hamberg所考虑的镜面是水平放置时,但实际工作中的光学元件放置方式多种多样。由Tribble AC[17]等的研究可以得知,当光学元件垂直于水平面放置时,其工作表面上所附着的粒子数量大约是工作面水平向上放置时的十分之一,当元件的工作面是水平向下放置时,其工作表面附着的粒子数量则大约为工作面水平向上放置时的百分之一。图9给出了不同的放置条件下,元件表面吸收随时间的变化。
分析图9可知,在空气洁净度一定的情况下,曝露时间越长,元件表面洁净度也越大,进而造成镜面吸收的增大。进一步比较不同波长的情况表明,颗粒吸收主要受折射率虚部的影响。此外,元件不同的放置方式也会对吸收产生很大影响。
在实际工作中,光学元件表面洁净度必须控制在某一范围,才能保证系统的性能保持基本不变,而当表面洁净度超过这个洁净度等级时,系统的工作性能将会迅速降低,此时系统所能容忍的洁净等级就是污染容限。当空气洁净度一定时,空气洁净度越大达到污染容限所用的时间也越短,进而,可以推测出达到污染容限所要的时间,并以此作为对元件进行清洗及维护所要求时间容限提供理论基础。
表2~4给出了入射光波长为0.633 μm和1.15 μm时,不同空气洁净度、工作面放置方式以及放置时间对元件表面吸收的影响。
表 2 不同空气洁净度等级吸收系数
Table 2. Absorption efficiency of different air cleanliness grades
Air class 1 000 10 000 100 000 Absorption(0.633 μm) 2.80×10−8 1.92×10−7 1.31×10−6 Absorption (1.15 μm) 3.46×10−8 2.37×10−7 1.62×10−6 表 3 不同朝向时的吸收系数
Table 3. Absorption efficient for different placement
Placement Horizontal downward surface Verticle surface Horizontal upward surface Absorption
(0.633 μm)1.94×10−10 2.33×10−9 2.80×10−8 Absorption
(1.15 μm)2.39×10−10 2.88×10−9 3.46×10−8 表 4 不同曝露时间下的吸收系数
Table 4. Absorption efficiency at exposure times
Exposure time 1 d 10 d 20 d Absorption(0.633 μm) 2.80×10−8 3.37×10−7 7.13×10−7 Absorption(1.15 μm) 3.46×10−8 4.16×10−7 8.80×10−7 由表2~4可以看出,空气洁净度等级越高,工作面朝向上,放置在空气中的时间越长时,空气中的颗粒沉积在元件表面的数量也越多,这样会致使元件表面吸收系数越大。其中,光学元件的放置方式对表面吸收特性影响显著,使得元件表面吸收系数有量级上的差异。此外,不同入射波长对吸收的影响也有较大差异,分析其主要原因是不同波长时的污染颗粒折射率虚部的差异造成的。光学镜面污染颗粒的吸收属于界面微弱吸收,其吸收值相对较低。在参考文献[6]中,鲁江涛等通过热透镜法测量得到的界面吸收系数为10−7,论文表2~4中结果与其相近,说明模拟结果符合实际。需要指出的是,光学薄膜的吸收已达到10−5量级甚至更小,则此时镜面污染颗粒的吸收影响不可忽视。
Analysis of absorption characteristics of polluted particles on optical mirrors
-
摘要: 针对实际工作中的光学元件表面粒子污染,以米氏理论为基础,分析了在风沙地区空气中污染颗粒(以SiO2为主)对入射激光能量的吸收效应。此外,进一步分析空气洁净度、元件工作面的朝向以及放置时间等对光学元件表面吸收特性的影响。结果表明,空气洁净度等级越高、元件工作面朝上、放置时间越久,光学元件的表面吸收能力越强。研究成果可为实际应用中光学系统的镜面污染控制提供理论参考。Abstract: Aiming at the particle contamination on the surface of optical components in practical work, based on Mie theory, the absorption of incident laser energy by polluted particles (mainly SiO2) in the air in sandy areas was discussed. In addition, the effects of air cleanliness, orientation of working face and placement time on the surface absorption characteristics of optical components were further analyzed. The results show that the higher the air cleanliness level is, the optical component face up and the longer exposure time are, would make the surface absorption of the optical component stronger. The research results can provide guidance for pollution control of optical system in practical applications.
-
Key words:
- absorption /
- optical elements /
- Mie theory /
- cleanliness
-
表 1 不同波长下的灰尘折射率
Table 1. Refractive index of dust at different wavelengths
Item Value λ/μm 0.325 0.633 1.15 3.39 10.6 n 1.53+0.01i 1.53+0.000 5i 1.5+0.001i 1.5+0.02i 1.7+0.2i 表 2 不同空气洁净度等级吸收系数
Table 2. Absorption efficiency of different air cleanliness grades
Air class 1 000 10 000 100 000 Absorption(0.633 μm) 2.80×10−8 1.92×10−7 1.31×10−6 Absorption (1.15 μm) 3.46×10−8 2.37×10−7 1.62×10−6 表 3 不同朝向时的吸收系数
Table 3. Absorption efficient for different placement
Placement Horizontal downward surface Verticle surface Horizontal upward surface Absorption
(0.633 μm)1.94×10−10 2.33×10−9 2.80×10−8 Absorption
(1.15 μm)2.39×10−10 2.88×10−9 3.46×10−8 表 4 不同曝露时间下的吸收系数
Table 4. Absorption efficiency at exposure times
Exposure time 1 d 10 d 20 d Absorption(0.633 μm) 2.80×10−8 3.37×10−7 7.13×10−7 Absorption(1.15 μm) 3.46×10−8 4.16×10−7 8.80×10−7 -
[1] 肖静, 张彬. 光学元件污染对红外光学系统信噪比的影响分析[J]. 红外与激光工程, 2012, 41(4): 1010−1016. doi: 10.3969/j.issn.1007-2276.2012.04.033 Xiao Jing, Zhang Bin. Influence of the optical components contamination on the signal to noise ratio in infrared optical systems [J]. Infrared and Laser Engineering, 2012, 41(4): 1010−1016. (in Chinese) doi: 10.3969/j.issn.1007-2276.2012.04.033 [2] 孙腾飞, 张骏, 吕海兵, 等. 光学镜面污染对激光传输特性的影响[J]. 红外与激光工程, 2014, 43(5): 1444−1448. doi: 10.3969/j.issn.1007-2276.2014.05.017 Sun Tengfei, Zhang Jun, Lv Haibing, et al. Influence of optical mirror pollution on laser transmission characteristics [J]. Infrared and Laser Engineering, 2014, 43(5): 1444−1448. (in Chinese) doi: 10.3969/j.issn.1007-2276.2014.05.017 [3] 张晓荣. 光学元件吸收损耗的高灵敏度检测技术研究[D]. 成都: 中国科学院光电技术研究所, 2015. Zhang Xiaorong. Study on highly sensitive testing method for weak absorptance of optical componrnts[D]. Chengdu: Institute of Optics and Electronics, Chinese Academy of Sciences, 2015. [4] 苏锴隆. CO2激光器的光学元件: 吸收是主要作用[J]. 激光与光电子学进展, 1984, 36(3): 19−24. Su Kailong. Optical components of CO2 lasers: absorption is the main role [J]. Laser & Optoelectronics Progress, 1984, 36(3): 19−24. (in Chinese) [5] 邱荣. 强激光诱导光学元件损伤的研究[D]. 绵阳: 中国工程物理研究院, 2013. Qiu Rong. Damage of optical elements induced by high power laser[D]. Mianyang: China Academy of Engineering Physics, 2013. (in Chinese) [6] 鲁江涛, 程鑫彬, 沈正祥, 等. 单层膜体吸收与界面吸收研究[J]. 物理学报, 2011, 60(4): 047802. doi: 10.7498/aps.60.047802 Lu Jiangtao, Cheng Xingbin, Shen Zhenxiang, et al. Volume and interface absorptions of single layer [J]. Journal of Physics, 2011, 60(4): 047802. (in Chinese) doi: 10.7498/aps.60.047802 [7] Brown R G W. Absorption and Scattering of Light by Small Particles[M]. New York: Wiley, 1983. [8] Jennings S G, Pinnick R G, Auvermann H J. Effects of particulate complex refractive index and particle size distribution variations on atmospheric extinction and absorption for visible through middle ir wavelengths [J]. Applied Optics, 1978, 17(24): 3922−3929. doi: 10.1364/AO.17.003922 [9] Genin F, Kozlowski M, Brusasco R. Catastrophic failure of contaminated fused silica optics at 355 nm[C]// 2nd Annual International Conference on Solid-State Lasers for Applications to Inertial Confinement Fusion, 1996: 1-9. [10] Bloembergen N. Role of cracks, pores, and absorbing inclusions on laser induced damage threshold at surfaces of transparent dielectrics [J]. Applied Optics, 1973, 12(4): 661−664. doi: 10.1364/AO.12.000661 [11] 王毅, 许乔, 柴立群, 等. 熔石英表面划痕附近电磁场分布模拟分析[J]. 强激光与粒子束, 2005, 17(1): 67−70. Wang Yi, Xu Qian, Chai Liqun, et al. Simulation and analysis of electromagnetic field distributing around fused silica scratch [J]. High Power Laser and Particle Beams, 2005, 17(1): 67−70. (in Chinese) [12] 章春来, 王治国, 向霞, 等. 熔石英后表面坑点型划痕对光场调制的近场模拟[J]. 物理学报, 2012, 61(11): 114210. doi: 10.7498/aps.61.114210 Zhang Chunlai, Wang Zhiguo, Xiang Xia, et al. Simulation of field intensification induced by pit-shaped crack on fused silica rear-surface [J]. Acta Physica Sinica, 2012, 61(11): 114210. (in Chinese) doi: 10.7498/aps.61.114210 [13] Mätzler C. MATLAB functions for Mie scattering and absorption[D]. Bern, Switzerland: Institute of Applied Physics, University of Bern, 2002. [14] Miroshnichenko A E, Tribelsky M I. Ultimate absorption in light scattering by a finite obstacle [J]. Physical Review Letters, 2018, 120(3): 033902. doi: 10.1103/PhysRevLett.120.033902 [15] Dittman M G. Contamination scatter functions for stray-light analysis[C]// Proceedings of SPIE - The International Society for Optical Engineering, 2002, 4774: 99-100. [16] Hamberg O, Shon E M. Particle size distribution on surfaces in clean rooms[R]. California: Vehicle Engineering Division of the Group of the Aerospace Corporation, 1984. [17] Tribble A C, Boyadjian B, Davis J, et al. Contamination control engineering design guidelines for the aerospace community[C]// SPIE's 1996 International Symposium on Optical Science, Engineering, and Instrumentation. International Society for Optics and Photonics, 1996. [18] 施泽明, 倪师军, 张成江. 成都市近地表大气尘的地球化学特征[J]. 地球与环境, 2004, 32(3): 53−58. Shi Zeming, Ni Shijun, Zhang Chengjiang. The Geochemical characteristics of air dust near the ground in Chengdu [J]. Earth and Environment, 2004, 32(3): 53−58. (in Chinese)