焦宏飛，汲小川，張錦龍，程鑫彬，王占山
（同濟大學物理科學與工程學院精密光學工程技術研究所，先進微結構材料教育部重點實驗室，上海市數字光學前沿科學研究基地，上海市全光譜高性能光學薄膜器件與應用專業技術服務平臺，上海200092）
摘要：高光譜性能光學薄膜是國家重大光學工程、光電子產業的基石。為了提高光學薄膜的光學性能并開發精確制備技術，現有研究主要圍繞設計理念、監控技術和薄膜材料等方面展開。目前，薄膜技術已經取得了長足進步，能實現高光譜性能光學薄膜的魯棒性設計，多種基于高光譜性能光學薄膜的精確制備技術相繼被提出。從薄膜設計、精確制備技術以及薄膜材料幾方面出發，本文對現代高光譜性能光學薄膜研究進行了回顧和討論，并對高性能光學薄膜潛在的挑戰和進一步研究方向進行了展望。
關鍵詞：光學薄膜及器件；魯棒性設計；光學監控；薄膜制備；計算制造
高光譜性能光學薄膜技術是現代光學和光電子系統的基石。現代大科學裝置的發展對光學薄膜器件的光譜性能提出了越來越高的要求，然而，這種光學薄膜的材料生長及加工難度極大、成本非常高。如用于引力波探測裝置的光學干涉鏡［1］、空間光學系統的高質量薄膜［2］、用于激光聚變裝置、激光對抗系統中的薄膜器件［3］等。同時，在消費電子、光纖通信、生物醫學和激光制造等領域［45］，薄膜器件的制造效率是影響其成本、競爭力的重要因素。因此，現代光學薄膜制造技術面臨著超高性能、高效率的迫切需求。
近幾十年來，在重大光學工程、光學產業需求的牽引下，光學薄膜技術取得了長足的進步。在薄膜設計方面，基于20世紀50年代FlorinAbelès提出的薄膜光譜系數的矩陣計算方法［6］，人們提出了GradualEvolution［7-8］，Flip-FlopTechnique［9-10］，NeedleOptimization［11-13］，GeneticAlgorithms［14-16］等一系列薄膜設計算法。得益于現代薄膜設計方法和商業軟件［17-19］的發展，設計人員可以完成光譜性能要求苛刻的復雜膜系設計，因此找到全局最優的解析解顯得不再那么迫切。然而，在設計階段考慮誤差敏感性，提高設計的魯棒性是一個不容忽視的問題［20-22］。一種普遍的魯棒性設計方法是通過在價值函數中增加懲罰項來降低誤差靈敏度［23］。隨著薄膜沉積過程研究的深入，在更全面的誤差因素上建立的計算制造實驗有望幫助設計者篩選出最佳的薄膜設計［24］。然而，如何有效進行極端性能光學薄膜的魯棒性設計，仍然是一個亟待解決的問題。
如果說光學薄膜的魯棒性設計是現代高性能光學薄膜制造技術的關鍵問題，那么薄膜厚度的精確制備則是生產高性能光學薄膜的核心問題。電子束蒸發［25-26］、離子束濺射［27-28］等沉積技術的發展使薄膜制造的穩定性不斷提升，為制備高性能光學薄膜提供了技術保障。鍍膜生產中，石英晶振監控被廣泛應用于各種沉積工藝［29］。在一些高能沉積工藝條件下，允許通過時間監控薄膜厚度來制備色散薄膜［30］。然而，由于光學監控對薄膜光學厚度的監測精度更高［31-32］，光學監控依然是研究熱點。單波長監控方法與光學薄膜有著一樣悠久的歷史。目前，單波長監控方式已經衍生出各種不同的監控策略，如極值點監控、差值監控、擺動值監控等［33-34］。在制備多腔窄帶濾光薄膜時，單波長極值點監控策略表現出極強的誤差自補償效應［35］。最新的研究表明，在使用其他單波長監控策略制備非規整的光學薄膜時也存在誤差自補償效應［36］。現代計算機處理數據能力的巨大提升以及光電二極管和CCD探測器的發展，為寬光譜監控技術的應用開辟了新的道路。與單波長監控相比，寬光譜監控技術具有對測試數據隨機誤差靈敏度低的優點［37］，同時由于其具有寬帶光譜反演分析的特性，國內外多家研究所已經開展了基于寬光譜監控技術的薄膜制造策略的研究。寬光譜監控中存在某種誤差自補償效應［38-39］，可減小監控誤差對薄膜光譜性能的負面影響。在近30多年的研究中，人們逐漸認識到，光學薄膜生產中監控策略的選擇很大程度上取決于薄膜的種類［40-41］。
總的來說，現代高光譜性能光學薄膜的精確制備技術面臨著光學薄膜的魯棒性設計和監控策略的正確應用兩方面問題。當前，高性能光學薄膜的高效、確定性制備仍然是光學薄膜領域的難點問題和研究熱點。
本文回顧了光學薄膜這一領域的近期發展，首先介紹光學薄膜的魯棒性設計的發展，其次介紹基于高光譜性能光學薄膜精確監控的確定性制造技術，然后將介紹混合膜領域以及rugate薄膜的相關生產，最后對高光譜性能光學薄膜鍍膜技術進行了總結和展望。
1高性能光學薄膜的設計
以目前的薄膜設計軟件和專業知識，當以光譜性能為唯一標準時，找到薄膜設計問題的解決方案是相對簡單的。然而，在設計階段找到更簡單且對制造誤差容忍度更高的設計方案則更具實際意義。制造誤差通常包括由沉積設備引起的隨機誤差、系統誤差、監控誤差以及材料參數誤差等因素，因此，在設計階段考慮并降低誤差敏感性并非易事。
1.1光學薄膜的脫敏設計
解決光學薄膜的魯棒性問題最常見的做法是先以光譜性能要求生成設計，并檢查其誤差靈敏度。通常脫敏設計的簡單方法是，抑制設計中的超薄層、優化材料、甚至逆轉設計，將敏感層放在最后。早期，加拿大國家研究院的Do?browolski提出一種在生成設計階段降低靈敏度的方法，即在評價函數中添加懲罰項［42］，表達式如下：
MF=MF0+α∑twi∂MF∂xi，
其中：MF0是不受擾動的評價函數，wi是表示結構參數xi的預期誤差水平的權重，α是用于權衡光譜性能和誤差的可調常數（0<α<1）。應用此方法設計的增益平坦濾光薄膜降低了對隨機厚度誤差的靈敏度，制備該濾光片的成功率達到98%［43］。
2003年，Tikhonravov等在此基礎上提出一種變體［45］，將二階項添加到泰勒級數中，簡化后的二階表達式如下：
MF=MF0+a2∑t(wi∂MF∂xi)2.
此方法已經在商業薄膜設計軟件中實現，并在2011年成功應用于色散薄膜的脫敏設計中［46］。
圖1（a）增益平坦濾光薄膜的魯棒性設計［44］；（b）色散薄膜的脫敏設計［45］
此外，光學薄膜的設計結構與薄膜材料的選擇也會直接影響多層薄膜的制備過程及難易程度。當光線傾斜入射光學薄膜時不可避免地會產生偏振效應，即S偏振光與P偏振光特性產生分離，光學系統性能下降，因此需要通過光學薄膜的設計實現對偏振的調控。中國科學院上海技術物理研究所針對航天任務的應用，設計了多種偏振效應可調的分色鏡和增透膜。蔡清元等提出了一種基于介質-金屬-介質膜堆的消偏振濾光薄膜的簡單設計方法［46］。段微波等利用光在金屬膜中傳播的傍軸傾向特性，選用金屬銀和Ta2O5，SiO2作為光學薄膜的材料設計用于空間環境下量子通信的保偏反射膜［47］。基于金屬-介質設計的光學薄膜對厚度誤差敏感性低，對膜厚控制要求低，然而金屬吸收率高，光譜效率低等問題限制它在強激光中的應用。因此，多層介質薄膜成為應用于強激光環境下的最佳選擇。Baumeister使用駐波比技術來匹配反射帶，設計了一個45°的消偏振帶通濾波器，然而這種設計多達7種材料［48］，給制備帶來了困難。Thelen發現了另一種基于失諧F-P腔的新方法［49］。然而，這種濾光薄膜通常需要不規整的膜層結構實現高透過率，制作過程復雜繁瑣，同時由于薄膜結構的低容差和高敏感性，受監控誤差和誤差累積效應的影響，通帶效率會降低，不利于薄膜的高效制備。同濟大學焦宏飛等提出了一種消偏振長波通濾光薄膜的設計初始結構，適當優化薄膜結構后降低了層厚度誤差對光譜特性的敏感度［50］。圖2所示為兩種設計結構與沉積實驗的光譜曲線對比，實驗結果表明，所提出的消偏振設計結構對誤差敏感性更低。
圖2新型消偏振結構和基于F-P腔優化結構在45°入射的設計和制備光譜［50］
1.2基于監控特性的魯棒性設計
在進行薄膜的魯棒性設計時，除了考慮沉積設備造成的隨機或系統誤差外，由薄膜生長檢測引起的厚度誤差也是不容忽視的重要因素。這種厚度誤差量級很大程度上取決于監控方式以及薄膜種類。因此其模擬更復雜［51］，模擬中將沉積過程中的厚度誤差與監測信號、沉積速率和擋板延遲等沉積參數關聯，由此模擬程序生成的設計比純數值結果更加可靠。早在1972年，Macleod等就提出了一種將設計與光學監控模擬結合生成的魯棒性設計［52］。然而30年后，Tikhonravov等才給出關于這種魯棒性設計的物理解釋，并基于此提出一種波分復用濾光薄膜的自動化設計方法，再結合單波長拐點監控策略能夠獲得十分穩定的生產結果［53］。Trubetskov等提出通過適當約束如厚度、層數等條件，生成眾多可行的設計，隨后進行模擬沉積實驗，從中篩選出最佳設計［54］。人們將這一步驟稱為計算制造實驗［55-56］。該方法不僅有助于選擇最佳設計，而且為設計的工程反演和采用合適的膜厚監控策略提供了重要信息。
1.3負濾光片的設計
負濾光片可選擇性地抑制波段，在較短波長和較長波長處具有良好的透過率。Bovard等提出連續調制折射率的方法設計負濾光片［57］，此方法可以實現優異的光譜特性［58-59］。然而，梯度折射率濾光薄膜的實現和高效制備仍然是一個難題，無法普及［60］。采用兩種材料設計負濾光片是目前主流的方法，Thelen使用等效層概念開發了一種在抑制帶兩側具有相對平滑的高透過率區域的負濾光片設計方法［61］。Yong基于天線理論的類比，描述了一種改善透過率的窄帶負濾光片設計方法［62］。然而，受材料折射率的限制，兩種材料中有非常薄的層，很難精確制備，需要采用離子束濺射技術，這導致制備成本高，殘余應力大等問題。
同濟大學張錦龍等提出使用常規四分之一波膜堆的二次諧波的高反區域來設計負濾光片［63］，如圖3所示，通過變跡法調制厚度可以有效抑制通帶的旁瓣，使用這種方法設計的負濾光片膜層厚度適中，能夠在沉積過程中進行精確監控。
圖3變跡法調制負濾光片的膜層結構及設計光譜［63］
1.4抑制非均質性半波孔的設計方法
除了鍍膜中可能產生的厚度誤差，材料參數的偏差也會嚴重影響薄膜的光譜特性。如電子束蒸發工藝制備的HfO2薄膜折射率具有非均質性，導致短波通濾光片在二倍頻的導納發生顯著變化，破壞原有的匹配，產生半波孔現象［64］。一些研究指出，改變鍍膜沉積參數甚至沉積工藝有利于得到均勻薄膜［65-66］。然而，為了保證最佳的激光損傷特性，薄膜材料和沉積工藝參數無法改變。同濟大學焦宏飛等提出了兩種抑制半波孔影響的設計方法，一種是偏移中心波長法（見圖4（a）），即通過回避半波孔的辦法來設計膜系，這種設計方法相對直接，結構簡單，缺點是透射帶寬較窄，且對沉積精度要求較高。另一種方法是通過導納補償法抑制半波孔（見圖4（b）），即將非均質性為正（0.5LH0.5L）和非均質性為負（HL）的兩種膜堆組合，合成后的新膜堆（0.5LH0.5LHL）在半波處與材料的非均質性無關。
圖4（a）偏移中心波長法得到的倍頻分離膜；（b）導納補償法得到的倍頻分離膜［67］
2高性能光學薄膜的精確制備技術
薄膜厚度的精確制備是高效率生產高光譜性能光學薄膜的核心。隨著光學薄膜理論的發展，監控方式也從早期的目視法進化為技術成熟的單波長監控、寬光譜監控等光學監控技術。光學監控中，不同的監控技術和監控策略會引入不同特性的厚度誤差。因此，了解各類監控的優缺點是選擇最佳生產策略的關鍵。
2.1單波長監控
單波長監控系統目前廣泛配備于各類真空鍍膜設備中，單波長監控策略也衍生出極值監控策略、差值監控策略和擺動值監控策略等［68］。直接單波長拐點監控策略的一大優點是在制備由1/4波長厚度或其厚度整數倍組成的窄帶濾光薄膜時表現出極強的誤差自補償效應。因此，這一監控技術被廣泛應用于波分復用中濾光片的大量生產中［35］。Tikhonravov等為這種誤差自補償效應提供了一種物理解釋［53］，但這種監控策略在制備非1/4監控波長光學厚度的多層薄膜時沒有明顯的優勢。因此，非規整的多層薄膜通常采用單波長擺動值監控策略。李正中等提出了一種選擇敏感監控波長的方法，可以對確定的光學薄膜設計生成監控波長序列，這種監控策略可以有效削弱監控誤差對光學特性的負面影響［69］。張錦龍等在間接單波長監控制備高性能薄膜的研制中，通過設計-制造-反演分析的制造流程，獲得了制造過程中薄膜誤差的演化規律，進而改進單波長監控序列，成功制備了具有數十層結構的高性能光學薄膜［70］，如圖5所示。
間接單波長監控一個明顯的缺點是必須要抓取準確的Tooling值以重新計算薄膜厚度。直接單波長監控則可以避免這類問題，但是其缺點是層厚度誤差會隨著監控層數的增加而累積，即誤差累積效應［71］。Zöller等提出一種結合間接和直接監控優點的綜合監控方式，將原固定的直接監控片改進成監控芯片切換器，這種切換器一次可安裝4個監控片，從而抑止誤差累積效應的發展［72］。另一項研究表明，采用非拐點監控策略制備非規整多層薄膜時存在誤差自補償效應，通過短波通濾光薄膜的計算制造實驗證明，擺動值監控策略比電平監控策略的誤差自補償效應更強［73］。
圖5（a）修正Tooling值后制備光譜與理論高度吻合；（b）由逆向工程反演確定低折射率層的相對誤差［70］
2.2寬光譜監控
沉積過程中監測生長薄膜的光譜透過率的技術發展可以追溯到20世紀60年代初，Vidal等［74］將半自動單色儀與沉積設備耦合，對置于夾具中心的基板的光譜特性進行重復測試，但當時計算機性能較低，只能分析測試光譜的部分數據點。隨著計算機的發展，1994年Tilsch等［75］首次將光譜測試應用在離子束濺射系統中。由于緊湊型分光光度計和高速CCD相機的發展，在線寬光譜監控系統已經集成到先進的鍍膜設備中［76］。圖6為配有直接寬光譜監控系統的離子束濺射設備示意圖。
圖6配有直接寬光譜監控系統的離子束濺射設備示意圖
近年來，寬光譜監控技術成為研究熱點，主要原因是寬光譜監控很好地解決了非規整膜系的監控難題［77-78］。與單波長監控相比，寬光譜監控具有對測量數據隨機誤差靈敏度低的優點［79］。除此之外，直接寬光譜監控還存在兩個重要特征，一種是厚度誤差隨層數的增加而累積［80］，在對一種40層短波通濾光片的模擬沉積中發現，隨著膜層數的增加，厚度誤差水平逐漸升高，后面層的相對厚度誤差是前幾層的近10倍［68］。另一種重要的特征是厚度誤差的相關性，先前層的監控誤差會影響監控信號，導致當前層的誤差依賴于先前層的誤差，厚度誤差相關性的積極作用是產生厚度誤差自補償效應［81］。Pelletier等在對非規整的多層膜進行計算制造時首次發現寬光譜監控中誤差的自補償效應［82］。2017年，Zhupanov等證明寬光譜監控中存在極強的誤差自補償效應［83］。由ZrO2和SiO2兩種材料設計一種布儒斯特角偏振濾光薄膜，由于ZrO2薄膜的折射率存在不穩定性，因此ZrO2薄膜折射率的變化也會引起較大的厚度誤差。實驗結果顯示，厚度誤差隨著沉積誤差數量的增加，可以清晰地觀察到誤差累積效應，部分膜層的相對厚度誤差高達16%，如此大的厚度誤差會導致光譜性能完全失效。然而，由于監控過程將層厚度誤差關聯在一起，從而提供了非常強的誤差自補償效應，最終成功制備了偏振片。
后續研究指出，不同類型的光學薄膜都可能存在誤差自補償效應［84］。這種效應不僅取決于光學薄膜的類型，而且取決于光學薄膜設計的具體選擇。因此，對光譜特性接近的薄膜設計，進行預生產分析，選擇最合適寬光譜監控策略，以提升復雜的光學薄膜制備成功率。莫斯科國立大學計算研究中心認識到由光學監控引起的各膜層厚度誤差之間存在相關性，厚度誤差相關性可能是引起誤差自補償效應的直接原因［85］，并提出了一種評估厚度誤差相關性強度的方法［86］。
圖7（a）偏振片生產過程中的厚度誤差；（b）此誤差下S和P偏振光的透過率曲線（實線），沒有厚度誤差的S和P偏振光的透過率曲線（虛線）；5種厚度誤差誤差不相關且平均誤差與a相同的設計；（c）S偏振光的透過率曲線；（d）P偏振光的透過率曲線［83］
同濟大學團隊則通過計算制造與真實沉積實驗，驗證了基于諧振腔結構的超陡度二向色鏡在直接寬光譜監控條件下具有較強的誤差自補償效應。如圖8所示，采用直接寬光譜監控制備的高敏感性超陡度二向色鏡的設計細節、模擬和測試結果，最終制得陡度為8nm，透過效率為95.01%，反射效率為95.52%的超陡度二向色鏡。除此之外，我們通過厚度誤差相關與非相關沉積實驗比較，證明了厚度誤差相關性在誤差自補償效應中的積極作用［87］，如圖9所示。
圖8（a）超陡度二向色鏡的設計；（b）膜層敏感性分析；（c）基于模擬計算預測誤差自補償效應強度；（d）寬光譜監控制備二向色鏡的理論和實際結果對比［87］
圖9（a）沉積膜層的層厚度相對誤差（RUN1厚度誤差相關，RUN2厚度誤差非相關）；（b）厚度誤差相關沉積實驗結果；（c）厚度誤差非相關沉積實驗結果
2.3其他監控方式
光學薄膜的光譜性能要求越來越高，傳統的光學監控方式已無法成為現代光學薄膜最有效的監控手段。例如，色散鏡是超快激光系統中控制色散的關鍵元件，其色散帶寬和調控能力都被嚴格要求。傳統的光學監控方式僅能通過記錄的透過率/反射率信息優化膜層的光學參數或厚度，無法提供足夠的自由度來完全補償GDD中發生的偏差。德國漢諾威激光中心將邁克爾遜干涉儀與IBS鍍膜設備相結合，實現了在沉積過程中對光學薄膜的群延遲色散進行原位測量［88］。圖10為相位測試系統的示意圖。
圖10相位測試系統的示意圖［88］
除此之外，混合監控策略是制備高光譜性能光學薄膜的重要手段。在制備復雜的多層膜時，部分膜層采用光學監控，其他層的沉積則使用非光學監控方式進行，通常包括石英晶振監控或者時間監控。德國MaxPlank研究所在磁控濺射方面的研究表明，采用時間監控方式在制備應用于超短脈沖激光器的色散薄膜方面具有突出的可靠性［89］。中國科學院上海技術物理研究所采用單波長與時間監控混合模式，成功研制出應用于空間儀器，效率在70%左右的亞納米帶寬的光學薄膜器件［90］。
同濟大學團隊系統研究了石英監控誤差的來源，明確系統誤差是造成薄膜光學性能下降的主要原因，并通過修正系統誤差實現了超寬帶增透薄膜制造［91］。隨后，我們應用雙離子束濺射技術，采用寬光譜監控與時間監控策略，成功制備兩種折射率材料光學厚度比接近10∶1，超薄層為15nm的四通道負濾光薄膜。還結合單波長監控與時間監控，實現了寬帶截止的窄帶濾光片的精確制備。值得一提的是，通過對行星系統掩膜板的高精度修正，210mm口徑內的薄膜均勻性可控制在0.1%以內。
圖11（a）應用石英晶振監控技術制備寬帶增透薄膜［91］；（b）應用混合監控模式成功制備四通道負濾光片；（c）210mm大口徑窄帶濾光元件；（d）寬帶外截止的超窄帶帶濾光片（210mm口徑內中心波長偏差小于0.1%）
3混合膜及rugate薄膜制備技術
由一系列具有恒定光學性質的離散層組成的層結構，只能代表一種特殊的光學薄膜。因此，光學薄膜可以理解為折射率在深度上任意連續變化的結構。基于此，可以實現一種理想的抗反射薄膜結構，其折射率從基板不斷下降到環境的折射率。漸變性設計除了可以實現優異的光譜性能外，在溫度和激光穩定性方面更具優勢［92］。這類具有連續變化折射率的光學薄膜被稱為Rugate薄膜。
3.1混合膜的制備
針對Rugate薄膜，科研人員對已有的光學薄膜制備技術開展了大量的研究，以確定合適的生產策略。漸變折射率的技術實現，通常采用共沉積兩種材料混合膜或兩種反應氣體組分的工藝。日本NTT應用電子實驗室將兩種反應氣體的混合，利用濺射技術制備了具有確定X組分的（SiO2）X（Si3N4）1-X混合膜，并將其成功應用于減反射薄膜的制備中［93］。德國漢諾威激光中心應用兩個獨立的電子束蒸發源，實現了對兩種材料沉積速率的獨立控制，大大提升了材料混合比例的可控性［92］。此外，離子束共濺射技術也應用于混合膜和Rugate薄膜的生產中［94］。圖12給出了共濺射和共蒸發工藝制備混合膜示意圖。
圖12離子束共濺射［94］及電子束共蒸發工藝制備混合膜示意圖
現有研究表明，利用共蒸發、共濺射等制備工藝在折射率薄膜中摻雜非晶、低吸收、寬帶隙的低折射率材料形成混合膜，不僅可以有效抑制薄膜結晶，降低薄膜的吸收，還可以拓寬薄膜的能帶隙，提高薄膜的本征損傷閾值［95-96］。Tokas［97］等報道了EB-Hf1-xSixO2混合膜在不同SiO2組分下的形態演變規律，認為Hf1-xSixO2混合膜不僅可以抑制薄膜的結晶，優化薄膜的微觀結構，從而有效降低紫外反射薄膜的散射損耗，還可以拓寬薄膜的能帶隙。Jensen［98］等通過實驗證明IBSHf1-xSixO2混合膜可以有效降低HfO2薄膜吸收，提升薄膜在納秒和飛秒脈沖作用下的激光損傷閾值。Lappschies等通過比較納米層與共濺射層制備的厚度一致的單層膜，證明共濺射制備的三元氧化物的帶隙向UV區偏移了近20nm，如圖13所示。這種帶隙向短波偏移的特性為提高紫外光譜范圍高功率激光薄膜的研制開辟了道路［99］。
圖13由TiO2和SiO2制備的納米層與共濺射層的透射光譜與純TiO2薄膜的消光系數比較［99］
同濟大學團隊系統地研究了SiO2摻雜含量和退火溫度對Hf1-xSixO2混合膜的微觀結構特性的影響［100］，在共蒸發工藝實現了對SiO2摻雜比例的精確控制。圖14展示了一系列不同SiO2組分Hf1-xSixO2混合膜的光學參數。此外，同濟大學還建立了Hf1-xSixO2薄膜表面形貌、表面粗糙度以及薄膜散射損耗之間的耦合關系，明確了薄膜結構缺陷、化學計量比失調缺陷與薄膜吸收之間的聯系，成功制備了低吸收、低散射損耗和高光譜效率的高性能反射薄膜。
圖14采用Optilayer擬合得到的不同SiO2組分Hf1-xSixO2混合膜的折射率［100］及消光系數
3.2 Rugate濾光薄膜
基于混合膜制備工藝的研究基礎，迄今為止已經形成幾種有效的沉積方法來制備梯度折射率層。俄勒岡州立大學團隊采用等離子體增強化學氣相沉積方法，以SiH4，N2和N2O為反應氣體制備了折射率從1.48~2.05呈正弦變化的ru?gate濾光薄膜［101］。Tang等提出用兩級高速反應濺射法制備由TaxSiyOz非均質薄膜組成的rugate濾光薄膜，最終制備的負rugate濾光薄膜具有低吸收和非晶性質，并且光譜特性優良，如圖15所示［102］。
圖15兩級高速反應濺射法制備的窄帶負Rugate濾光片的光譜特性［102］
德國漢諾威激光中心將離子束共濺射技術與寬光譜監控技術完美結合，對混合材料厚度進行精確控制，Rugate薄膜的理論光譜與制備結果具有良好的一致性。圖16為中心波長為800nm，帶寬40nm，0~55°的全向抗反射薄膜。
圖16測試與設計的159層全向抗反射薄膜光譜及全向抗反射薄膜800nm處的角譜［94］
同濟大學團隊提出了一種改進的傅里葉變換合成方法，該方法首先固定有效控制折射率范圍的二次Q位相，通過迭代法改變光譜函數，提高了傅里葉變換設計方法的通用性，保證了設計結果的準確性和可制備性［103］。圖17展示了在限制折射率的條件下，經過170次迭代理想光譜與合成光譜之間基本重合。
圖17類屋反射器的折射率分布和反射光譜［103］
此外，同濟大學團隊提出通過傾斜沉積控制薄膜孔隙率，從而生長出具有正弦波折射率特征的Rugate薄膜的工藝［104］。圖18展示了利用傾斜沉積技術實現多阻帶濾光薄膜的設計，這種方法的優勢是僅需要單一材料就可以實現多阻帶濾光薄膜制備，不過傾斜薄膜的控制精度有待進一步提升。
圖18多阻帶Rugate薄膜的設計和透射光譜［104］
4結論與展望
高性能光學薄膜的高效、確定性制備一直是光學薄膜領域的難點問題和研究熱點。近幾十年來，在重大光學工程、光學產業需求的牽引下，無論是高性能光學薄膜的設計理念還是復雜薄膜結構的精確制備技術都取得了長足的進步和發展。本文回顧了用于高性能光學薄膜的魯棒性設計方法，對制造誤差的脫敏性設計是一個具有實際意義的重要課題。早期，薄膜脫敏性設計是通過添加懲罰項來約束薄膜光學性能對隨機誤差的靈敏度，然而這種設計缺少對沉積過程中誤差因素的考慮。目前，考慮真實沉積誤差（監控誤差、速率波動、擋板延遲等其他相關參數）對目標光譜性能影響的優化設計已經實現。
高性能光學薄膜的精確制備離不開光學監控的發展，到目前為止，各種各樣的光學監控策略被提出，然而沒有一種通用的策略可以在所有情況下提供精確的厚度監控。本文介紹了監控策略的優缺點，可以幫助做出正確的選擇，對于某些類型的光學薄膜，監控策略的選擇與光學薄膜的設計相關，例如多腔窄帶濾光薄膜。基于此理念，計算制造實驗在現代光學薄膜生產中發揮著越來越重要的作用。計算制造實驗能夠反映出誤差累積、厚度誤差相關性、誤差自補償等一系列信息，這些信息不僅有利于薄膜魯棒性設計，而且為定制合適的膜厚監控策略提供了重要參考，最大限度地減少研發和迭代的時間。
折射率隨深度任意連續變化的Rugate薄膜的制備是另一挑戰。混合膜作為一種折射率可控的特殊材料，可以用于制備Rugate薄膜，為光學薄膜的設計提供了新的自由度。在離子束共濺射工藝基礎上，我們成功將寬光譜監控應用于rugate薄膜的制備。
在未來，薄膜設計和監控方式的協同優化會成為高光譜性能薄膜的基礎，光學監控策略和數據分析算法的進步將實現更高精度薄膜的監控，納米復合材料、微結構材料等新型薄膜材料生長技術也會推動高光譜性能薄膜的發展。
參考文獻：
［1］ ACCADIA T，ACERNESE F，ANTonUCCI F，et al. Performance of the Virgo interferometer longi?tudinal control system during the second science Run［J］. Astroparticle Physics，2011，34（7）：521-527.
［2］ SOLANKI S K， DEL TORO INIESTA J C，WOCH J，et al. The polarimetric and helioseismicimager for Solar orbiter：so/PHI［J］. Proceedings of the International Astronomical unio，2014，10（S305）：108-113.
［3］ SAKHAROV V K. Model of lock-in in a ring laser and a semiconductor laser gyro［J］. Technical Physics，2011，56（8）：1135-1141.
［4］ TIKHonRAVOV A，KOCHIKOV I，SHARAP OVA S. Broad band optical monitoring in the pro? duction of gain flattening filters for telecommunication applications［J］. Moscow University Physics Bulletin，2019，74（2）：160-164.
［5］ VENGHAUS H. Wavelength Filters in Fibre Optics［M］. Springer，2006，123：71-116.
［6］ BORN M，WOLF E. Principles of Optics；Electromagnetic Theory of Propagation，Interference and Diffraction of Light［M］. 2d rev. ed. New York：Macmillan，1964.
［7］ DOBROWOLSKI J A. Completely automatic syn?thesis of optical thin film systems［J］. Applied Op?tics，1965，4（8）：937-946.
［8］ DOBROWOLSKI J A. Versatile computer pro?gram for absorbing optical thin film systems［J］. Ap?plied Optics，1981，20（1）：74-81.
［9］ SOUTHWELL W H. Coating design using very thin high- and low-index layers［J］. Applied Optics，1985，24（4）：457-460.
［10］ DOBROWOLSKI J A. Comparison of the Fourier transform and flip-flop thin-film synthesis methods ［J］. Applied Optics，1986，25（12）：1966-1972.
［11］ TIKHonRAVOV A V，BAUMEISTER P W，POPOV K V. Phase properties of multilayers［J］.Applied Optics，1997，36（19）：4382-4392.
［12］ TIKHonRAVOV A V，TRUBETSKOV M K，DEBELL G W. Optical coating design approaches based on the needle optimization technique［J］. Applied Optics，2007，46（5）：704-710.
［13］ SULLIVAN B T，DOBROWOLSKI J A. Imple?mentation of a numerical needle method for thinfilm design［J］. Applied Optics，1996，35（28）：5484-5492.
［14］ BINDA P D，ZOCCHI F E. Genetic algorithm op?timization of X-ray multilayer coatings［C］. Opti?cal Science and Technology，the SPIE 49th Annu?al Meeting. Proc SPIE 5536，Advances in Compu?tational Methods for X-Ray and Neutron Optics，Denver，Colorado，USA. 2004，5536：97-108.
［15］ AL-MARZOUG S M，HODGSON R J. Optimization of multilayer mirrors at 13. 4 nm with more than two materials［J］. Applied Optics，2008，47（12）：2155-2160.
［16］ WANG Z，KAISER N，LEQUIME M，et al. Non-periodic multilayer coatings in EUV，soft xray and x-ray range［J］. International Society for Optics and Photonics，2008，7101：710110.
［17］ WILLEY R R. Practical Design and Production of Optical Thin Films［M］. New York： M. Dekker，1996.
［18］ OptiLayer. com［EB/OL］. http：//www. optilay?er. com
［19］ BAUMEISTER P W. Optical Coating Technology［M］. SPIE Press，2004.
［20］ BAUMEISTER P. evaluation of the solutions for two design problems presented at the 1998 optical interference coatings conference［J］. Applied Optics，2000，39（13）：2230-2234.
［21］ THELEN A，TILSCH M，TIKHonRAVOV AV，et al. Topical Meeting on Optical interference Coatings（OIC’2001）：design contest results［J］.Applied Optics，2002，41（16）：3022-3038.
［22］ TILSCH M，HENDRIX K，VERLY P. Optical interference coatings design contest 2004［J］. Ap?plied Optics，2006，45（7）：1544-1554.
［23］ PERVAK V， TIKHonRAVOV A V， TRUBETSKOV M K，et al. 1. 5-octave chirped mirror for pulse compression down to sub-3 fs［J］. Applied Physics B，2007，87（1）：5-12.
［24］ TIKHonRAVOV A V，TRUBETSKOV M K. Modern design tools and a new paradigm in optical coating design［J］. Applied Optics， 2012， 51（30）：7319-7332.
［25］ HARPER J M E，CUOMO J J，KAUFMAN HR. Technology and applications of broad ?beam ion sources used in sputtering. Part II. Applications［J］. Journal of Vacuum Science and Technology，1982，21（3）：737-756.
［26］ KAUFMAN H R. Developments in broad-beam，ion source technology and applications［J］. Journal Of Vacuum Science and Technology，1982，21（3）：1982，764-767.
［27］ HSU J C，LEE C C. Single- and dual-ion-beam sputter deposition of titanium oxide films［J］. Applied Optics，1998，37（7）：1171-1176.
［28］ LEE C C，HSU J C，WONG D H. Low loss niobium oxides film deposited by ion beam sputter deposition［J］. Optical and Quantum Electronics，2000，32（3）：327-337.
［29］ FRAZIER G A，GLOSSER R. Phase diagrams of thin films of the palladium hydrogen system using a quartz crystal thickness monitor［J］. Journal of Physics D：Applied Physics，1979，12（10）：L113-L115.
［30］ LIU H X，XIONG S M，LI L H，et al. Variation of the deposition rate during ion beam sputter deposition of optical thin films［J］. Thin Solid Films，2005，484（1/2）：170-173.
［31］ MACLEOD H A. Monitoring of optical coatings［J］. Applied Optics，1981，20（1）：82-89.
［32］ FURMAN S，TIKHonRAVOV A V. Basics of optics of multilayer systems［J］. 1996.
［33］ BAUMEISTER P W. Optical Coating Technology［M］. SPIE Press，2004.
［34］ ZHAO F T. Monitoring of periodic multilayers by the level method［J］. Applied Optics，1985，24（20）：3339-3342.
［35］ MACLEOD H A. Turning value monitoring of narrow-band all-dielectric thin-film optical filters［J］. Optica Acta：International Journal of Optics，1972，19（1）：1-28.
［36］ TIKHonRAVOV A V， LAGUTINA A A，LAGUTIN I S，et al. Self-compensation of errors in optical coating production with monochromatic monitoring［J］. Optics Express，2021，29（26）：44275.
［37］ TIKHonRAVOV A V，TRUBETSKOV M K，KOKAREV M A，et al. Effect of systematic errors in spectral photometric data on the accuracy of determination of optical parameters of dielectric thin films［J］. Applied Optics，2002，41（13）：2555-2560.
［38］ TIKHonRAVOV A V，TRUBETSKOV M K，AMOTCHKINA T V. Does broadband optical monitoring provide an error self-compensation mechanism？ ［C］. Optical Interference Coatings 2010，Tucson，Arizona United States，6-11 June 2010，2010：TuC3.
［39］ TIKHonRAVOV A，KOCHIKOV I，YAGO LA A. Mathematical investigation of the error selfcompensation mechanism in optical coating technology［J］. Inverse Problems in Science and Engineering，2018，26（8）：1214-1229.
［40］ TIKHonRAVOV A V， KOCHIKOV I V，YAGOLA A G. Investigation of the error selfcompensation effect associated with direct broad band monitoring of coating production［J］. Optics Express，2018，26（19）：24964-24972.
［41］ TIKHonRAVOV A V， KOCHIKOV I V，YAGOLA A G. Error self-compensation mechanism in the optical coating production with direct broad band monitoring［J］. Optics Express，2017，25（22）：27225-27233.
［42］ DOBROWOLSKI J A. Versatile computer program for absorbing optical thin film systems［J］.Applied Optics，1981，20（1）：74-81.
［43］ DOBROWOLSKI J A，DALACU D，LI L，etal. Fifty years of optical interference coatings at the national research council of Canada［J］. Optics and Photonics News，2007，18（6）：24.
［44］ TIKHonRAVOV A V，TRUBETSKOV M K，AMOTCHKINA T V，et al. Application of advanced optimization concepts to the design of high quality optical coatings［C］. Proceedings of SPIEThe International Society for Optical Engineering，2003，4829：1061-1062.
［45］ PERVAK V，TRUBETSKOV M K，TIKHON RAVOV A V. Robust synthesis of dispersive mirrors［C］. SPIE Optical Systems Design. Proc SPIE 8168，Advances in Optical Thin Films IV，Marseille，France. 2011，8168：95-103.
［46］ CAI Q Y，LUO H H，ZHENG Y X，et al. Design of non-polarizing cut-off filters based on dielectric-metal-dielectric stacks ［J］. Optics Express，2013，21（16）：19163-19172.
［47］ 余德明，段微波，李大琪，等 . 偏振和相位調控反射 鏡 的 設 計 與 制 備［J］. 光 學 學 報 ，2020，40（15）：222-226.YU D M，DUAN W B，LI D Q，et al. Design and fabrication of polarization- and phase-modulated mirror［J］. Acta Optica Sinica，2020，40（15）：222-226.（in Chinese）
［48］ BAUMEISTER P. Bandpass design-applications to nonnormal incidence［J］. Applied Optics，1992，31（4）：504-512.
［49］ THELEN A. Nonpolarizing edge filters［J］. Journal of the Optical Society of America，1981，71（3）：309-314.
［50］ JIAO H F，NIU X S，ZHANG X M，et al. Design and fabrication of a superior nonpolarizing longwavelength pass edge filter applied in laser beam combining technology［J］. Applied Optics，2020，59（5）：A162-A166.
［51］ TIKHonRAVOV A V，TRUBETSKOV M K.Computational manufacturing as a bridge between design and production［J］. Applied Optics，2005，44（32）：6877-6884.
［52］ MACLEOD H A. Turning value monitoring of narrow-band all-dielectric thin-film optical filters［J］. Optica Acta：International Journal of Optics，1972，19（1）：1-28.
［53］ TIKHonRAVOV A V，TRUBETSKOV M K. Automated design and sensitivity analysis of wavelengh-division multiplexing filters［J］. Applied Optics，2002，41（16）：3176-3182.
［54］ TIKHonRAVOV A V，TRUBETSKOV M K. Modern design tools and a new paradigm in optical coating design［J］. Applied Optics， 2012， 51（30）：7319-7332.
［55］ AMOTCHKINA T V， SCHLICHTING S，EHLERS H，et al. Computational manufacturing as a tool for the selecion of the most manufacturable design［J］. Applied Optics，2012，51（36）：
8677-8686.
［56］ WILBRANDT S，STENZEL O，KAISER N. All-oxide broadband antireflection coatings by plasma ion assisted deposition： design， simulation，manufacturing and re-optimization［J］. Optics Express，2010，18（19）：19732-19742.
［57］ BOVARD B G. Rugate filter design：the modified Fourier transform technique［J］. Applied Optics，1990，29（1）：24-30.
［58］ SOUTHWELL W H. Using apodization functions to reduce sidelobes in rugate filters［J］. Applied Optics，1989，28（23）：5091-5094.
［59］ TIKHonRAVOV A， TRUBETSKOV M，AMOCHKINA T，et al. New optimization algorithm for the synthesis of rugate optical coatings［C］. Optical Systems Design 2005. Proc SPIE 5963，Advances in Optical Thin Films II，Jena，Germany. 2005，5963：18-26.
［60］ LAPPSCHIES M，GÖRTZ B，RISTAU D. Application of optical broadband monitoring to quasirugate filters by ion-beam sputtering［J］. Applied Optics，2006，45（7）：1502-1506.
［61］ THELEN A. Design of optical minus filters［J］.Journal of the Optical Society of America，1971，61（3）：365-369.
［62］ YOUNG L. Multilayer interference filters with narrow stop bands［J］. Applied Optics，1967，6（2）：297-315.
［63］ ZHANG J L，XIE Y J，CHENG X B，et al. Thinfilm thickness-modulated designs for optical minus filter［J］. Applied Optics，2013，52（23）：5788-5793.
［64］ MACLEOD H A. Half wave holes［C］. Leaks and Other Problems in Proceedings of the 39th Annual Technical conference of the Society of Vacuum Coaters. Society of Vacuum Coaters， 1996：193-198.
［65］ ALVISI M，DE TOMASI F，PERRONE M R，et al. Laser damage dependence on structural and optical properties of ion-assisted HfO2 thin films［J］. Thin Solid Films，2001，396（1/2）：44-52.
［66］ WILBRANDT S，STENZEL O，KAISER N.All-optical in situ analysis of PIAD deposition processes［C］. Advances in Optical Thin Films III，SPIE Proceedings. Glasgow，Scotland，United Kingdom. SPIE，2008，7101（3）：71010D.
［67］ JIAO H F，CHENG X B，BAO G H，et al.Study of HfO2/SiO2 dichroic laser mirrors with refractive index inhomogeneity［J］. Applied Optics，2014，53（4）：A56-A61.
［68］ TIKHonRAVOV A V，TRUBETSKOV M K，AMOTCHKINA T V. Statistical approach to choosing a strategy of monochromatic monitoring of optical coating production［J］. Applied Optics，2006，45（30）：7863-7870.
［69］ LEE C C，WU K，KUO C C，et al. Improvement of the optical coating process by cutting layers with sensitive monitor wavelengths［J］. Optics Express，2005，13（13）：4854-4861.
［70］ ZHANG J L，TIKHonRAVOV A V，LIU Y L，et al. Design，production and reverse engineering of ultra-steep hot mirrors［J］. Optics Express，2014，22（11）：13448-13453.
［71］ STECKELMACHER W. Thin-film optical filters，3rd edition［J］. Vacuum，2002，66（1）：91-92.
［72］ ZOELLER A， HAGEDORN H， WEINRICHW，et al. Testglass changer for direct optical monitoring［C］. SPIE Optical Systems Design. ProcSPIE 8168，Advances in Optical Thin Films IV，Marseille，France. 2011，8168：375-380.
［73］ TIKHonRAVOV A V， LAGUTINA A A，LAGUTIN I S，et al. Self-compensation of errors in optical coating production with monochromatic monitoring［J］. Optics Express，2021，29（26）：44275-44282.
［74］ VIDAL B，FORNIER A，PELLETIER E. Optical monitoring of nonquarterwave multilayer filters［J］. Applied Optics，1978，17（7）：1038-1047.
［75］ TILSCH M，SCHEUER V，STAUB J，et al. Direct optical monitoring instrument with a double detection system for the control of multilayer systems from the visible to the near infrared［C］. Proc SPIE 2253，Optical Interference Coatings，1994，2253：414-422.
［76］ RISTAU D，GUENSTER S. Optimization of optical coatings for the UV/VUV-range［C］. SPIE Proceedings， Advanced Characterization Techniques for Optics，Semiconductors，and Nanotechnologies. San Diego，California，USA. SPIE，2003，5188：80-95.
［77］ VIDAL B，FORNIER A，PELLETIER E. Wideband optical monitoring of nonquarterwave multilayer filters［J］. Applied Optics，1979，18（22）：3851-3856.
［78］ LI L，YEN Y H. Wideband monitoring and measuring system for optical coatings［J］. Applied Optics，1989，28（14）：2889-2894.
［79］ TIKHonRAVOV A V，TRUBETSKOV M K，AMOTCHKINA T V. Investigation of the effect of accumulation of thickness errors in optical coating production by broadband optical monitoring ［J］. Applied Optics，2006，45（27）：7026-7034.
［80］ MACLEOD H A. Thin-Film Optical Filters［M］. CRC Press，2001.
［81］ TIKHonRAVOV A V，TRUBETSKOV M K，AMOTCHKINA T V. Investigation of the error self-compensation effect associated with broadband optical monitoring［J］. Applied Optics，2011，50（9）：C111-C116.
［82］ VIDAL B， PELLETIER E. Nonquarterwave multilayer filters：optical monitoring with a minicomputer allowing correction of thickness errors［J］. Applied Optics，1979，18（22）：3857-3862.
［83］ ZHUPANOV V，KOZLOV I，FEDOSEEV V，et al. Production of Brewster angle thin film polarizers using a ZrO2/SiO2 pair of materials［J］. Applied Optics，2017，56（4）：C30-C34.
［84］ TIKHonRAVOV A V，TRUBETSKOV M K，AMOTCHKINA T V. Investigation of the error self-compensation effect associated with broadband optical monitoring［J］. Applied Optics，2011，50（9）：C111-C116.
［85］ TIKHonRAVOV A V， KOCHIKOV I V，YAGOLA A G. Investigation of the error selfcompensation effect associated with direct broad band monitoring of coating production［J］. Optics Express，2018，26（19）：24964-24972.
［86］ TIKHonRAVOV A，KOCHIKOV I，SHARAPOVA S，et al. Optical monitoring of coating production：correlation of errors and errors self-com pensation ［C］. SPIE Optical Systems Design.Proc SPIE 11872， Advances in Optical Thin Films VII，online Only. 2021，11872：96-102.
［87］ JI X C，ZHANG J L，JIAO H F，et al. Production of ultra-steep dichroic filters with broad band optical monitoring［J］. Optics Express，2022，30（13）：22501-22511.
［88］ SCHLICHTING S，WILLEMSEN T，EHLERSH，et al. Fourier-transform spectral interferometry for in situ group delay dispersion monitoring of thin film coating processes［J］. Optics Express，2016，24（20）：22516-22527.
［89］ PERVAK V， TIKHonRAVOV A V， TRU BETSKOV M K，et al. Band filters：two-material technology versus rugate ［J］. Applied Optics，2007，46（8）：1190-1193.
［90］ 王凱旋，陳剛，劉定權，等.綠光波段60pm超窄帶濾光片的研制［J］.中國光學，2022，15（1）：119-131.
［91］ DONG S Y，JIAO H F，BAO G H，et al. Origin and compensation of deposition errors in a broadband antireflection coating prepared using quartz crystal monitoring［J］. Thin Solid Films，2018，660：54-58.
［92］ RISTAU D，AKHTAR S，EBERT J. Laser Induced Damage of Dielectric Phase Retardation Mirrors at 1. 064 Microns［M］. Laser Induced Damage in Optical Materials：1986. 100 Barr Harbor Drive，PO Box C700，West Conshohocken，PA 19428-2959：ASTM International，2009：300-300-13.
［93］ KATAGIRI Y，UKITA H. Ion beam sputtered（SiO（2））（x）（Si（3）N（4））（1-x） antireflection coatings on laser facets produced using O（2）-N（2）discharges［J］. Applied Optics，1990，29（34）：5074-5079.
［94］ BECKER H， TonOVA D， SUNDERMANNM，et al. Design and realization of advanced multiindex systems［J］. Applied Optics，2014，53（4）：A88-A95.
［95］ LAPPSCHIES M，GÖRTZ B，RISTAU D. Application of optical broadband monitoring to quasirugate filters by ion-beam sputtering［J］. Applied Optics，2006，45（7）：1502-1506.
［96］ MANGOTE B，GALLAIS L，COMMANDRÉM，et al. Femtosecond laser damage resistance of oxide and mixture oxide optical coatings［J］. Optics Letters，2012，37（9）：1478-1480.
［97］ TOKAS R B，SAHOO N K，THAKUR S，et al. A comparative morphological study of electron beam co-deposited binary optical thin films of HfO2：SiO2 and ZrO2：SiO2［J］. Current Applied
Physics，2008，8（5）：589-602.
［98］ JENSEN L O，MENDE M，BLASCHKE H，et al. Investigations on SiO2/HfO2 mixtures for nanosecond and femtosecond pulses［C］. Laser Damage Symposium XLII：Annual Symposium on Optical Materials for High Power Lasers. Proc SPIE 7842，Laser-Induced Damage in Optical Materials： 2010， Boulder， Colorado， USA. 2010，7842：68-77.
［99］ LAPPSCHIES M，JUPÉ M，RISTAU D. Extension of ion beam sputtered oxide mixtures into the UV spectral range［C］. Optical Interference Coatings，OSA Technical Digest，2007：TuA7.
［100］ JIAO H F，NIU X S，ZHANG J L，et al. Hf1-xSixO2 nanocomposite coatings prepared by ionassisted Co-evaporation process for low-loss and high-LIDT optics［J］. Materials（basel，Switzerland），2021，14（10）：2606.
［101］ LIM S，RYU J H，WAGER J F，et al. Rugate filters grown by plasma-enhanced chemical vapor deposition［J］. Thin Solid Films，1994，245（1/2）：141-145.
［102］ TANG Q，MATSUDA H，KIKUCHI K，et al. Fabrication and characteristics of rugate filters deposited by the TSH reactive sputtering method ［J］. Journal of Vacuum Science & Technology A： Vacuum， Surfaces， and Films， 1998， 16（6）：3384-3388.
［103］ CHENG X B，FAN B，DOBROWOLSKI J A，et al. Gradient-index optical filter synthesis with controllable and predictable refractive index profiles［J］. Optics Express，2008，16（4）：2315-2321.
［104］ ZHU Y M，JIAO H F. Rugate filter with multichannel grown by glancing angle deposition［J］.Optik，2012，123（16）：1501-1503.
玻纖含量對長玻纖...
鈣鈦礦薄膜的均勻...
用于光伏板靜電除...
聚砜醫療干粉吸入...
