蒙特卡洛（MC）模擬研究加氫石油樹脂結(jié)晶過程，是通過隨機抽樣與統(tǒng)計力學(xué)方法，在分子層面模擬樹脂分子鏈從無規(guī)線團向有序晶態(tài)轉(zhuǎn)變的成核、生長全過程，核心流程與關(guān)鍵要點如下：

一、模擬體系構(gòu)建

選取加氫石油樹脂的典型分子結(jié)構(gòu)（如氫化茚滿、氫化萘滿及其烷基衍生物），構(gòu)建含20~50個重復(fù)單元的分子鏈模型，采用旋轉(zhuǎn)異構(gòu)態(tài)描述鏈段柔性?；跇渲瑢嶋H密度（0.90~0.95g/cm³）搭建三維周期性立方模擬盒，初始構(gòu)象設(shè)為無規(guī)線團，分子鏈隨機分布。相互作用勢函數(shù)選用Lennard-Jones（L-J）勢描述鏈段間范德華力，含極性基團的改性樹脂需額外引入庫侖勢描述氫鍵作用，同時用諧振勢固定鍵長與鍵角。

二、核心模擬算法

采用鏈段位移、構(gòu)象翻轉(zhuǎn)等抽樣方式驅(qū)動分子鏈運動，通過Metropolis準(zhǔn)則判定構(gòu)象是否接受——能量降低的構(gòu)象直接接受，能量升高的構(gòu)象以概率P=exp(-△E/k_BT) 接受。通過鍵取向序參數(shù)q_4與密度序參數(shù)判定結(jié)晶狀態(tài)：q_4≥0.6且區(qū)域尺寸大于臨界晶核（5~10個鏈段）時，定義為穩(wěn)定晶核；晶核形成后，周圍無序鏈段通過構(gòu)象調(diào)整附著于晶核表面，實現(xiàn)晶粒生長。

三、模擬流程

先在樹脂熔融溫度以上（T_m+50℃）進行NVT系綜平衡模擬（10^6~10^7步），消除初始構(gòu)象干擾；再通過線性降溫（速率10^4~10^6K/s）將體系冷卻至目標(biāo)溫度，降溫過程中持續(xù)抽樣并統(tǒng)計有序參數(shù)、晶粒尺寸等數(shù)據(jù)；降溫結(jié)束后繼續(xù)模擬，直至結(jié)晶達到穩(wěn)定狀態(tài)。

四、關(guān)鍵影響因素與結(jié)果分析

模擬可定量揭示分子結(jié)構(gòu)與外界條件對結(jié)晶的調(diào)控規(guī)律：分子鏈支鏈越少、規(guī)整度越高，結(jié)晶度越高；長鏈樹脂結(jié)晶誘導(dǎo)期長但晶粒生長速率快。冷卻速率越快，晶粒尺寸越小，可實現(xiàn)樹脂微晶化；含羥基、酯基的改性樹脂，分子間氫鍵可促進鏈段有序排列，顯著提升結(jié)晶度與晶體穩(wěn)定性。最終通過統(tǒng)計不同時刻的結(jié)晶度、晶粒尺寸分布，繪制結(jié)晶動力學(xué)曲線，結(jié)合差示掃描量熱儀（DSC）、X射線衍射（XRD）等實驗數(shù)據(jù)驗證模擬結(jié)果。

該方法可彌補實驗難以實時觀測微觀結(jié)晶過程的不足，為加氫石油樹脂的分子結(jié)構(gòu)設(shè)計與加工工藝優(yōu)化提供理論依據(jù)。

本文來源：河南向榮石油化工有限公司 http://m.jieqite.com/