🔷電場如何驅動電子，電子又如何改變分子與化學反應？

🔷從碰撞截面、EEDF到流體模型，拆解電漿體模擬底層原理

在半導體製造、電漿體刻蝕、薄膜沉積以及電漿體催化等技術背後，都繞不開一個核心問題：我們究竟該如何把肉眼不可見的電子運動、分子碰撞和化學反應“算出來”？

這次小編帶大家一起學習這份課程，從電子碰撞、E/N、EEDF、振動激發與電漿體化學出發，一路深入到 BOLSIG+、ZDPlasKin、0D/1D/2D流體模型以及多時間尺度建模。

在小編看來，這份資料真正值得關注的地方，是它把“電子如何獲得能量—如何激發和解離分子—如何產生自由基—最終如何影響宏觀反應”這一整套電漿體建模邏輯串聯了起來。接下來，我們就沿著這條路線，看看複雜的電漿體世界究竟是如何被數學模型一步步“還原”出來的。

一、這份材料真正想說明什麼？

這是一份比較系統的電漿體建模課程資料。整份文件不是單純介紹某一種模擬軟體，而是試圖建立一套從微觀電子碰撞 → 電漿體化學 → 宏觀放電與流體 → 多時間尺度加速 → 電漿體催化多尺度耦合的完整建模框架。

二、電漿體建模的起點，是搞清楚“電子如何把能量傳遞給分子”

文件前半部分首先從電漿體化學（Plasma Chemistry）講起。在非平衡電漿體中，電子是驅動整個反應網路的關鍵角色，而電子碰撞反應並不簡單由氣體溫度決定，而是高度依賴電子能量分佈函數 EEDF、電子—分子碰撞截面以及約化電場 E/N。換句話說，E/N 決定電子獲得多少能量，而電子能量又決定那些化學反應通道能夠被打開。

電子與 N₂、O₂、CH₄、CO₂ 等分子碰撞後，會產生多種不同過程，包括電子激發、振動激發、解離、電離以及電子附著。這些過程進一步產生激發態分子、自由基、正負離子等高活性粒子，使得原本需要較高溫度才能進行的化學反應可以在非平衡條件下發生。例如電子碰撞 CH₄ 可以形成 CH₃、CH₂、CH 等活性碎片，而 O₂ 則可以形成 O、O(¹D) 等活性氧物種。

因此，這一部分真正想說明的是，電漿體不是簡單地把氣體“加熱”，而是在利用高能電子選擇性地把電能注入分子的電子、振動和化學自由度，從而重新設計化學反應路徑。

三、振動激發是非平衡電漿體化學中非常重要的一條“隱藏通道”

文件花了相當大的篇幅討論 Vibrational Excitation（振動激發）。電子碰撞可以把 N₂、O₂、CH₄、CO₂ 等分子激發到不同振動態，而這些振動態並不是簡單的中間產物，而是一個重要的能量記憶體和傳遞通道。

振動激發後的分子主要經歷三類過程：V-T relaxation（振動—平動能量交換）、V-V relaxation（不同分子振動態之間的能量交換）以及振動激發分子參與的化學反應。尤其是後者，可以提高反應速率、打開新的反應路徑並改變體系的熱釋放。

文件進一步利用 Fridman–Macheret α-Model解釋振動能為什麼能夠促進化學反應：振動能相當於提前給反應物“墊付”了一部分跨越勢壘所需要的能量，因此能夠有效降低反應活化能。這種作用對於強吸熱反應尤其明顯，而對於本身沒有明顯活化勢壘的放熱反應，作用則相對較弱。

這也是整份課程一個非常重要的物理思想，真正理解非熱電漿體，不能只看電子溫度和氣體溫度，還必須追蹤電子激發態、振動態、自由基和離子之間複雜的非平衡能量傳遞。

四、從微觀反應走向0D模型：BOLSIG+、ZDPlasKin與CHEMKIN構成核心計算鏈條

理解了碰撞截面和反應機理之後，課程進一步進入真正的電漿體化學建模。

這裡介紹了一套非常典型的建模體系：首先利用 BOLSIG+ 求解電子 Boltzmann 方程，獲得電子能量分佈函數，並進一步計算電子輸運參數以及電子碰撞反應速率；隨後使用 ZDPlasKin 建立0D電漿體動力學模型，計算各種粒子濃度以及氣體溫度隨時間的變化。ZDPlasKin本身就是針對複雜非熱電漿體化學體系設計的零維動力學求解器。

但只計算 plasma chemistry 還不夠，因為電漿體產生的自由基、激發態最終還要進入傳統燃燒化學網路。因此課程進一步把 ZDPlasKin 與 CHEMKIN 耦合：

BOLSIG+ → EEDF/反應速率係數 → ZDPlasKin → 電漿體物種演化 → CHEMKIN → 燃燒反應網路。

這樣就可以同時追蹤電子、激發態、自由基、溫度以及燃燒反應的演化。課程甚至討論了利用納秒脈衝放電和DC放電組合，分別控制電子激發與振動激發，從而實現對化學反應選擇性的控制。

五、真正困難的是跨尺度：從0D化學模型升級到1D/2D流體模型

當空間分佈不能忽略以後，問題迅速複雜化。

課程指出，電漿體輔助燃燒本質上是一個典型的多物理場、多長度尺度、多時間尺度、多組分強剛性問題。其中同時存在光子、電子、電子激發態、振動態、電磁場、聲波、衝擊波、點火波以及燃燒波；空間尺度則從 Debye 長度和鞘層一直跨越到反應區、混合層以及宏觀流場；時間尺度甚至可以橫跨 ps—ms。

因此課程進一步引入 1D/2D Fluid Model。此時需要同時求解 Poisson方程 + 帶電粒子輸運方程 + 電子能量方程 + 化學反應方程 + 質量守恆 + 動量守恆 + 能量守恆。

例如1D模型中，電子輸運參數和電子碰撞反應速率仍然可以由 BOLSIG 提供，但模型開始顯式描述電勢、電場、電子密度、離子密度、流體速度以及總能量等空間分佈。

到了2D模型，則進一步模擬納秒放電中 streamer（流注）傳播、電子密度分佈、局部電場、自由基生成以及氣體加熱。課程展示的結果表明，電壓脈衝結束以後電子密度會迅速下降，但 O、H、OH 等活性物種仍可以在後放電階段通過解離淬滅等過程繼續產生。

這意味著放電結束 ≠ 電漿體化學結束。後放電階段的自由基化學和能量釋放同樣可能決定最終的點火與燃燒行為。

六、多時間尺度演算法，是解決大型電漿體化學模型計算量爆炸的關鍵

隨著模型從幾十個反應擴展到幾百甚至上千個反應，一個嚴重問題出現了：不同物種的特徵時間尺度可能相差十幾個數量級。

有些電子激發態可能在皮秒或納秒尺度變化，而穩定分子和宏觀燃燒過程可能需要微秒甚至毫秒。因此，如果整個模型統一使用最快過程要求的時間步長，計算成本會極其巨大。

課程因此重點介紹 MTS（Multi-Time Scale）與 HMTS（Hybrid Multi-Time Scale）方法：根據不同物種的特徵時間，把反應體系劃分為 Fast / Medial / Slow Groups，不同組使用不同時間步長進行計算。

進一步又引入 DAC（Dynamic Adaptive Chemistry）以及 CO-DAC（Correlated DAC），利用時間和空間相關性動態減少需要即時計算的化學反應網路。這樣就把傳統的“所有位置、所有時間、所有反應全部計算”，轉變為只在真正需要的時候，計算真正重要的反應。

這實際上是整份課程在數值計算層面的核心思想——模型不僅要精準，還必須通過降階、自適應和多尺度演算法讓複雜電漿體模擬變得可計算。

七、最終目標是建立從原子尺度到反應器尺度的“多尺度電漿體數字模型”

整份文件最終指向的是 Multi-scale Plasma Modeling。

給出的總體架構圖實際上可以看作整份課程的“總路線圖”：最底層是 Atomic-scale Model（原子尺度電漿體—催化劑相互作用）；往上是 Plasma Chemistry Model（0D電漿體化學）；再進一步進入催化劑孔隙內部的 3D Micro-scale Model；最後上升到整個反應器的 3D Macro-scale Model。這些模型之間通過電場、溫度、物種通量、關鍵粒子濃度以及簡化化學機制不斷交換資訊，最終形成統一的 Integrated Model for Plasma Catalysis。

所以，這份課程真正想表達的，並不是“如何使用某一個電漿體軟體”，而是一套非常完整的建模哲學 電子碰撞截面 → EEDF → 電子碰撞反應 → 激發/解離/電離/附著 → 振動能量傳遞 → Plasma Chemistry → 0D動力學 → 1D/2D流體模型 → 多時間尺度演算法 → 多尺度電漿體/催化模型。

最終希望把電子尺度、分子尺度、化學反應尺度、放電尺度、流體尺度以及反應器尺度真正連接起來，從而預測“外部電場如何通過電子驅動微觀化學，最終改變宏觀反應過程”。

小編總結

系統建立了從電子碰撞與非平衡電漿體化學，到0D/1D/2D流體模擬，再到多時間尺度與電漿體催化多尺度耦合的完整建模框架，其核心目標是打通“電場—電子—激發態/自由基—化學反應—傳熱傳質—宏觀反應器”之間的跨尺度聯絡，實現對複雜電漿體過程的可預測模擬。(芯聯匯)