王越婷，賈如月，劉棟梁，邱俊，陳國偉，姜姍

(1.山東科技大學(xué) 安全與環(huán)境工程學(xué)院，山東 青島 266590；2.山東科技大學(xué) 化學(xué)與生物工程學(xué)院，山東 青島 266590)

膨潤土又名“微晶高嶺石”或“膠嶺石”，是一種含水的層狀鋁硅酸鹽礦物[1]，蒙脫石是其主要礦物成分，具有納米級的二八面體(TOT)結(jié)構(gòu)，即中間為鋁氧八面體，上下兩層為硅氧四面體。由于異價類質(zhì)同象替代的原因，蒙脫石上下粘土片層常帶有永久性的負(fù)電荷，可以吸附環(huán)境中的各類陽離子進入層間。若在其結(jié)構(gòu)層空間中插入有機改性劑，親水性的蒙脫石片層將會改為疏水性，對有機物分子的親和力也會相應(yīng)增加[2]，進而可以應(yīng)用于高分子聚合物、藥物緩釋材料的制備及有機污染物的吸附等眾多領(lǐng)域。

在該研究方向，以往研究者大多借用X射線衍射儀(XRD)、掃描電鏡(SEM)、透射電鏡(TEM)及傅里葉紅外光譜儀(FTIR)等多種分析儀器研究并推測有機蒙脫石的層間微環(huán)境特征，雖然取得了較大的進展，但由于粘土片層結(jié)構(gòu)的屏蔽作用，這些分析手段往往不能有效地得出原子以及分子尺度上蒙脫石礦物結(jié)構(gòu)片層間各類微環(huán)境特征和微觀結(jié)構(gòu)的信息[3]。而計算機分子模擬技術(shù)作為一種既可進行理論研究又可與實驗研究相互印證的微觀實驗方法，不僅可以建立模擬各物質(zhì)本身的分子結(jié)構(gòu)，對材料性能進行預(yù)測；而且也可以模擬物質(zhì)發(fā)生反應(yīng)時的動態(tài)變化過程，對宏觀觀察到的實驗現(xiàn)象進行解釋。因此其在材料科學(xué)、生命科學(xué)及化工科學(xué)等多個領(lǐng)域已經(jīng)得到了廣泛應(yīng)用，也愈來愈成為許多研究者進行科學(xué)研究的重要技術(shù)手段[4]。

1 分子動力學(xué)模擬方法概述

計算機分子模擬技術(shù)主要包括分子力學(xué)模擬技術(shù)(MM)、分子動力學(xué)模擬技術(shù)(MD)和蒙特卡洛模擬方法(MC)三種，而分子動力學(xué)模擬則是現(xiàn)階段對于有機蒙脫石研究中的應(yīng)用最為廣泛且重要的一種方法[5]。

分子動力學(xué)，簡稱MD模擬計算，起源于20世紀(jì)六七十年代，距今已有50多年的歷史。隨著計算機技術(shù)的飛速發(fā)展，分子動力學(xué)模擬技術(shù)已被廣泛應(yīng)用于物理、化學(xué)、醫(yī)藥、環(huán)境保護等眾多領(lǐng)域。

分子動力學(xué)模擬基本原理是將體系中的各分子(或者原子、離子)看作牛頓粒子[6]，在指定的系統(tǒng)內(nèi)，通過給定所有粒子一定的初始條件，利用求解牛頓運動方程的方法，往復(fù)循環(huán)得出模擬時間內(nèi)各粒子的運動軌跡，最后根據(jù)結(jié)果分析得到整個體系的力學(xué)、熱力學(xué)和動力學(xué)性質(zhì)，其計算容量一般以 5 000 個原子、100 ns為宜[7]。因此，與傳統(tǒng)研究方法相比，分子動力學(xué)模擬在研究微觀時間或空間尺度下物質(zhì)結(jié)構(gòu)與微環(huán)境的變化趨勢、機理時，具有無可比擬的優(yōu)勢。

分子動力學(xué)模擬一般包括以下幾個步驟[8-11]：①第一步也是整個模擬過程中最為重要的一環(huán)：為模擬研究建立恰當(dāng)?shù)哪Ｐ?，并對其進行幾何或者能量優(yōu)化；②確定系統(tǒng)的初始條件，包括模擬的力場類型、溫度、壓強、積分步長及時間設(shè)置等；③系統(tǒng)逐步趨于穩(wěn)定的模擬過程，在設(shè)定的初始條件下運行模擬至整個系統(tǒng)逐漸平衡；④根據(jù)模擬目的，選擇有針對性的參數(shù)文件(例如徑向分布函數(shù)、密度分布函數(shù)、自擴散系數(shù)等)進行分析，得出結(jié)論。

2 分子動力學(xué)模擬在有機蒙脫石合成中的應(yīng)用

現(xiàn)階段利用分子動力學(xué)模擬技術(shù)研究有機蒙脫石的成果主要體現(xiàn)在兩個方面，一是從三維角度直觀地展示表面活性劑在蒙脫石納米級層間空間內(nèi)所引起的結(jié)構(gòu)及環(huán)境變化；二是為宏觀觀察到的現(xiàn)象及性質(zhì)提供分子或原子水平上的深入理解[12]。

2.1 分子動力學(xué)模擬在有機蒙脫石層間微結(jié)構(gòu)和分子環(huán)境影響因素方面的研究

2.1.1 有機改性劑的種類和用量 對于不同種類的表面活性劑，由于其改性蒙脫石的機理不同，所以會導(dǎo)致改性蒙脫石的層間微結(jié)構(gòu)有很大差別。例如Liao等[13]根據(jù)Materials Studio(MS)軟件建立了一系列與XRD測得層間距數(shù)值相對應(yīng)的超晶胞有機蒙脫石模型，在Forcite模塊NVT系綜下模擬至整個系統(tǒng)平衡后，直觀探究了三類荷電性不同的表面活性劑連續(xù)改性引起蒙脫石結(jié)構(gòu)變化的機理：陽離子表面活性劑十六烷基三甲基溴化銨(CTAB)通過陽離子交換反應(yīng)取代蒙脫石層間金屬陽離子，之后在蒙脫石層間域內(nèi)呈雙層平鋪排列使蒙脫石層間距由1.26 nm增大到2.10 nm；而陰離子表面活性劑十二烷基硫酸鈉(SDS)插入后，由于其在蒙脫石層間域內(nèi)呈垂直排列促使有機蒙脫石的層間距驟增至 5.20 nm；但與前兩者不同的是非離子表面活性劑氫化蓖麻油(HCO)改性后的有機蒙脫石層間距反而從5.20 nm減小到4.10 nm，出現(xiàn)這種情況的原因可能是之前CTA+和SDS在蒙脫石層間域內(nèi)排列不規(guī)則、不穩(wěn)定所造成的，當(dāng)HCO參與改性后，三者在蒙脫石層間域內(nèi)的排列變得穩(wěn)定有序。Meng[14]和Martina等[15]也采用了根據(jù)XRD數(shù)據(jù)建立有機蒙脫石超晶胞模型的方法做了相似的研究，其中Meng等模擬了十六烷基三甲基氯化銨(CTAC)和十二烷基硫酸鈉(SDS)對蒙脫石的共改性過程，發(fā)現(xiàn)SDS由于自身帶有負(fù)電性的原因不能單獨插入到蒙脫石層間，而是與CTA+形成離子對進入到蒙脫石層間，并且兩種表面活性劑改性的順序不同，得到的有機蒙脫石層間距也不同；在兩性表面活性劑和陽離子表面活性劑改性蒙脫石的結(jié)構(gòu)差異方面，Martina等通過將實驗研究與建立相應(yīng)的有機蒙脫石模型相結(jié)合的方法發(fā)現(xiàn)磷脂酰膽堿分子(DSPC)可直接插入到蒙脫石層間域內(nèi)，而且分子動力學(xué)模擬的排列結(jié)果(圖1和圖2)顯示磷脂酰膽堿分子(DSPC)比十八烷基三甲基溴化銨分子(ODTMA)對蒙脫石層間距的影響更大。

圖1 ODTMA在蒙脫石層間域內(nèi)的排列方式Fig.1 The arrangement of ODTMA in the interlayer of montmorillonite a.雙層排列，d=1.90 nm；b.偽三層排列，d=2.10 nm

圖2 DSPC在蒙脫石層間域內(nèi)垂直排列(d=5.20 nm)Fig.2 The vertical arrangement of DSPC in the interlayer of montmorillonite (d=5.20 nm)

而對于同種類但尺寸大小不同的表面活性劑，Koteja等[16]首先模擬了在NPT系綜下達到恒溫恒壓時，不同鏈長、構(gòu)型的表面活性劑插入蒙脫石超晶胞后所形成的最優(yōu)結(jié)構(gòu)，結(jié)果表明對于表面活性劑碳鏈長度，蒙脫石層間距隨碳鏈長度的增大而增大；之后，又模擬了有機蒙脫石對水分子的吸附能變化，發(fā)現(xiàn)對于偶氮苯活性劑，雖然反式偶氮苯尺寸較大，但順式偶氮苯會加強蒙脫石對層間水的吸附從而使其層間距增大。Miroslav等[17]采用Cerius2軟件在UNIVERSAL立場下分別研究了兩種碳原子數(shù)相同但基團不同的表面活性劑十六烷基三甲基銨陽離子(CTA+)和十六烷基吡啶(CP)對蒙脫石層間結(jié)構(gòu)及能量變化的影響，發(fā)現(xiàn)與CTA-MMT(層間距d=1.81 nm)相比，CP-MMT改性陽離子與蒙脫石硅氧烷表面之間的靜電力作用較弱，吸附能也較低，因此層間距較大，為2.06 nm，模擬得到的層間距數(shù)據(jù)也與實驗結(jié)果較為吻合。

另一方面，為探究表面活性劑用量不同時蒙脫石層間結(jié)構(gòu)的變化規(guī)律，Lü等[18]根據(jù)蒙脫石對苯海拉明(DPH)的吸附實驗數(shù)據(jù)，用MS軟件建立了與之對應(yīng)的不同負(fù)載量的有機蒙脫石模型，經(jīng)過研究先幾何優(yōu)化再NVT系綜運行后的模擬結(jié)果，發(fā)現(xiàn)當(dāng)其吸附量在0.2～0.3 CEC時，DPH分子在蒙脫石層間域內(nèi)呈單層平臥排列，改性蒙脫石層間距d=1.50 nm；而當(dāng)DPH吸附量為在0.6～0.7 CEC時，DPH分子在蒙脫石層間域中呈雙層排布，蒙脫石d值增大到1.70 nm；Li等[19]以新型Gemini表面活性劑丙基雙十六烷基二甲基陽離子為模擬對象，采用MS中Discover模塊研究了其在蒙脫石層間域內(nèi)排列結(jié)構(gòu)隨用量的變化，發(fā)現(xiàn)當(dāng)表面活性劑用量從0.2 CEC增大到1.0 CEC時，改性蒙脫石層間距從1.53 nm逐漸增大到3.76 nm，而當(dāng)表面活性劑用量從1.0 CEC增大到1.5 CEC時，改性蒙脫石層間距不再發(fā)生變化，在此過程中表面活性劑在蒙脫石層間域內(nèi)逐漸從雙層排列向三層排列、石蠟型傾斜雙層排列轉(zhuǎn)變。Schampera等[20]在比較十六烷基三甲基溴化銨(HDTMA)和十六烷基吡啶(HDPy)對蒙脫石層間結(jié)構(gòu)的影響時，對蒙脫石、有機陽離子和氯離子分別采用了LAMMPS模塊下粘土礦物專用CLAYFF力場和OPLS-AA力場，發(fā)現(xiàn)在低負(fù)載量(<1.0 CEC)的條件下，兩種有機陽離子大致呈單層排列于蒙脫石層間域內(nèi)，在高負(fù)載量(>1.0 CEC)條件下，有機陽離子在蒙脫石層間域內(nèi)呈雙層排列結(jié)構(gòu)，但HDTMA+在蒙脫石層間域內(nèi)形成的有機片層密度卻明顯大于HDPy+所形成的，與采用UNIVERSAL力場模擬的結(jié)果相比，這一模擬結(jié)果則更接近實際情況。除研究蒙脫石層間域內(nèi)有機分子自身的排列變化規(guī)律外，利用MD模擬技術(shù)也可以直觀地觀察表面活性劑分子對蒙脫石層間微環(huán)境的影響，例如Zhao等[21]在采用LAMMPS模塊中CLAYFF力場及SPC水模型的條件下，通過NPT系綜及NVE系綜模擬研究十六烷基三甲基銨陽離子改性有機蒙脫石的動力學(xué)變化時，發(fā)現(xiàn)疏水性有機陽離子與蒙脫石片層的相互作用明顯強于無機金屬離子與蒙脫石片層的相互作用，層間有機陽離子插入量越多，則蒙脫石表面疏水作用越強，靠近蒙脫石片層的水分子就越少；該模擬結(jié)果也從側(cè)面證明了上述不同類型的力場能夠達到很好的契合效果。

根據(jù)以上研究內(nèi)容可以看出，分子動力學(xué)模擬主要依托MS和Cerius2兩款軟件中的Forcite模塊或LAMMPS模塊實現(xiàn)，在研究時有UNIVERSAL力場和CLAYFF力場等力場可供選擇，不過就有機蒙脫石而言，CLAYFF力場模擬結(jié)果更接近實際情況。分子動力學(xué)模擬技術(shù)與常規(guī)實驗手段息息相關(guān)，但又有所區(qū)別，例如有機蒙脫石超晶胞模型的建立需要根據(jù)基礎(chǔ)XRD實驗數(shù)據(jù)，但利用該模擬技術(shù)可直接觀察表面活性劑分子在有機蒙脫石內(nèi)的排列狀態(tài)及引起的能量和環(huán)境變化，所獲得的結(jié)果與常規(guī)測試手段相比則更為直觀。

模擬結(jié)論顯示不同種類的表面活性劑對蒙脫石的改性機理不同，陽離子、非離子表面活性劑均可直接對蒙脫石進行插層改性，而陰離子表面活性劑卻不能單獨對蒙脫石進行改性；隨著表面活性劑鏈長、用量的增加，蒙脫石的層間距會逐漸增大、層間疏水作用也會逐漸增強；在碳鏈長度相同，但基團不同或構(gòu)型不同的情況下，表面活性劑與蒙脫石片層間的相互作用力則對改性蒙脫石的結(jié)構(gòu)起主要影響。

2.1.2 晶體化學(xué)性質(zhì) 晶體化學(xué)性質(zhì)也被認(rèn)為是影響蒙脫石有機改性過程的關(guān)鍵因素之一[3]，而利用分子動力學(xué)模擬技術(shù)可以很好地從分子/原子視角研究其對蒙脫石層間結(jié)構(gòu)和微環(huán)境的影響。例如在探究層間金屬陽離子對蒙脫石有機化的影響時，Wu等[22]采用分子動力學(xué)模擬技術(shù)觀察研究有機蒙脫石的層間微結(jié)構(gòu)時，發(fā)現(xiàn)在Na+、K+、Mg2+、Ca2+和Fe3+中，水化作用強、與蒙脫石片層間相互作用弱的金屬離子(Ca2+、Fe3+)更容易被有機陽離子/分子取代。對于改性后的蒙脫石結(jié)構(gòu)和熱力學(xué)性質(zhì)的變化，Anoukou等[23]主要以八面體Mg2+取代Al3+為蒙脫石層電荷的來源，進而對構(gòu)建的蒙脫石超晶胞模型作研究，發(fā)現(xiàn)當(dāng)蒙脫石層間陽離子交換容量由68 meq/100 g增大至91 meq/100 g時，有機改性劑在層間域內(nèi)的插入量會隨之增多，導(dǎo)致其層間距和體積都會相應(yīng)增大。而Shen等[24]分別將CLAYFF力場、CVFF力場和SPC水的模型應(yīng)用于蒙脫石、四甲基銨陽離子(TMA+)和水分子，然后以有機蒙脫石(TMA-MMT)為基礎(chǔ)研究了不同層電荷條件下蒙脫石的水化性能，結(jié)果表明隨著層電荷的增加，水分子逐漸被硅氧烷表面吸引，這不僅會導(dǎo)致氫鍵數(shù)目增多，而且也會使水化能增大，從而減少水分子的流動性，這與Fu等[25]模擬的結(jié)論一致。

綜合上述模擬結(jié)果，在經(jīng)過按照實驗?zāi)康慕?、幾何?yōu)化至動力學(xué)模擬后，研究者可以直觀地觀察到有機蒙脫石的層間微環(huán)境結(jié)構(gòu)及分子運動狀態(tài)，且這種模擬研究方法既高效、簡便、經(jīng)濟，又可探究XRD、FTIR、SEM等常規(guī)研究手段觸及不到的領(lǐng)域，具有傳統(tǒng)實驗方法無可比擬的優(yōu)勢。分子動力學(xué)模擬結(jié)果表明蒙脫石的有機化改性難易程度與層間金屬離子種類有很大關(guān)系，且蒙脫石的層電荷越高，其對有機陽離子的吸附能力越高以及改性后的蒙脫石層間距就越大，層間水分子的流動性則就越差。

2.1.3 水分子的含量 有機蒙脫石層間域中通常含有不同含量的水分子，利用分子動力學(xué)模擬技術(shù)可以對不同水分子含量的影響進行最大程度的仿真研究。魏景明等[26]模擬了十六烷基三甲基溴化銨改性蒙脫石的分子空間分布特征，發(fā)現(xiàn)1.0 CEC改性的有機蒙脫石在含水量低的條件下，表面活性劑分子主要在層間域內(nèi)呈雙層平臥排列，而隨著含水量的升高，表面活性劑在層間域內(nèi)會逐漸從雙層分布向均勻分布轉(zhuǎn)變；2.0 CEC改性的有機蒙脫石在含水量低的情況下，表面活性劑分子在層間域內(nèi)呈雙層傾斜分布，但極性基團靠近四面體片層，在含水量高的條件下，表面活性劑分子整體會向?qū)娱g域中心移動，此時水分子則占據(jù)靠近硅氧烷表面的位置。

綜上所述，現(xiàn)階段分子動力學(xué)模擬技術(shù)主要應(yīng)用于研究表面活性劑、層電荷及水分子含量等因素對有機蒙脫石晶體結(jié)構(gòu)及層間域內(nèi)客體分子的影響，得到的結(jié)果對于補充有機蒙脫石的基礎(chǔ)理論及指導(dǎo)實踐工作都具有重要意義。不過，對于一些外部條件如溫度、壓強、pH等環(huán)境因素以及一些更為復(fù)雜的情況如類質(zhì)同象取代位置的不同等晶體化學(xué)因素對有機蒙脫石層間結(jié)構(gòu)和微環(huán)境等產(chǎn)生的影響則鮮有報道。

2.2 分子動力學(xué)模擬在有機蒙脫石應(yīng)用性能方面的研究

聚合物/粘土納米復(fù)合材料(PCNs)因具有穩(wěn)定、耐熱、阻隔及阻燃[31]等眾多優(yōu)異特性，一直是有機蒙脫石應(yīng)用領(lǐng)域的熱點問題，而近年來由于分子動力學(xué)模擬技術(shù)的發(fā)展，已逐漸有學(xué)者利用該技術(shù)從微觀角度對復(fù)合物的結(jié)構(gòu)及性質(zhì)進行了探索研究。如Wang等[32]利用分子動力學(xué)模擬手段研究了半纖維素木聚糖/蒙脫石納米復(fù)合材料的膨脹和穩(wěn)定機理，發(fā)現(xiàn)含水量較低時，木聚糖通過在不影響材料結(jié)構(gòu)的情況下促進粘土片層形成納米級孔隙，進而抑制蒙脫石片層的膨脹；而在含水量高的條件下，木聚糖會由平鋪構(gòu)象轉(zhuǎn)變?yōu)槭D型構(gòu)象，與粘土片層相連接致使蒙脫石層間距增大。Anoukou等[33]分別將CLAYFF力場和CHARMM力場應(yīng)用于蒙脫石模型和MT2EtOH、尼龍模型的構(gòu)建中，之后先經(jīng)NVE系綜模擬使體系能量達到最小化后，再經(jīng)NPT系綜模擬尼龍合成聚合物粘土納米復(fù)合材料的結(jié)構(gòu)的動力學(xué)變化，發(fā)現(xiàn)與聚合物相比，有機改性劑與蒙脫石片層間的吸引作用更為強烈，而聚合物鏈中單體的數(shù)量則對復(fù)合材料的層間距和彈性常數(shù)沒有明顯影響。Scocchi等[34]先以Discover模塊中COMPASS力場建立了二牛脂基二甲基氯化銨(M2(C18)2)和牛脂基三甲基氯化銨(M3C18)改性蒙脫石的最優(yōu)模型，后在NVT系綜下模擬觀察了微觀形態(tài)下研究了尼龍聚合有機蒙脫石體系的相互作用規(guī)律，結(jié)果表明表面活性劑內(nèi)的極性基團被吸引到硅氧烷表面，非極性的尾鏈則平鋪于蒙脫石層間域內(nèi)，但與雙鏈M2(C18)2不同的是，單鏈M3C18對蒙脫石片層表面的占據(jù)度較低一些。

目前，利用分子動力學(xué)模擬研究蒙脫石作醫(yī)用藥物載體的可行性亦頗受關(guān)注。Sun等[35]以 Forcite 模塊模擬氯化十六烷基三甲基咪唑為改性劑時，其改性沸石、伊利石和蒙脫石后對氯霉素(CP)的吸附差異，發(fā)現(xiàn)改性沸石和伊利石對CP的吸附位點僅局限于外表面，而改性蒙脫石對CP的吸附位點除了外表面之外，在層間域內(nèi)也大量存在，且對CP的吸附量隨層電荷的增大而增大。Iannuccelli等[36]以分子動力學(xué)手段表征了蒙脫石負(fù)載醫(yī)用藥物艮他霉素(GM)的層間結(jié)構(gòu)變化，結(jié)果表明當(dāng)GM添加量為0.79 CEC時，其在蒙脫石層間域內(nèi)呈單層平臥排列，改性蒙脫石的層間距由原來的 1.49 nm 縮小到1.42 nm；而Wu等[37]以CASTEP模塊模擬研究蒙脫石對環(huán)丙沙星(CIP)的吸附性時，發(fā)現(xiàn)蒙脫石層電荷越高，其對CIP的靜電引力作用就越大，CIP越不容易從蒙脫石層間域內(nèi)釋放，這些都從模擬的角度證明了蒙脫石用作藥物緩釋材料的可行性。

在應(yīng)用性能方面，利用分子動力學(xué)模擬技術(shù)可以從微觀視角研究有機蒙脫石吸附污染物、負(fù)載醫(yī)用藥物及合成聚合物/粘土納米復(fù)合材料的機理，并對這些實驗現(xiàn)象加以解釋，模擬結(jié)果表明，利用表面活性劑改性的蒙脫石不僅對一些有機物具有很強的吸附能力，而且自身物理化學(xué)性質(zhì)也會有很大提高，可以為有機蒙脫石的大規(guī)模工業(yè)應(yīng)用提供理論指導(dǎo)和技術(shù)支持。不過以上多表現(xiàn)為利用分子動力學(xué)模擬技術(shù)對有機蒙脫石吸附性能的研究，而對其流變性、凝膠性等性能的探索有待進一步展開。

3 結(jié)束語

本文綜述了利用分子動力學(xué)模擬技術(shù)研究表面活性劑、蒙脫石晶體化學(xué)性質(zhì)等因素對有機蒙脫石合成的影響和有機蒙脫石在不同應(yīng)用領(lǐng)域的進展情況。研究發(fā)現(xiàn)，表面活性劑的電荷、尺寸、用量以及蒙脫石的層電荷等因素都對有機蒙脫石的晶體結(jié)構(gòu)和層間微環(huán)境變化有重要的影響，且有機蒙脫石在吸附環(huán)境污染物、負(fù)載引用藥物及合成高分子材料方面也有很好的應(yīng)用前景。

另外，由于分子動力學(xué)模擬技術(shù)可通過計算機軟件直接完成，因此具有方便、經(jīng)濟成本低及安全等優(yōu)點。且與常規(guī)XRD及FTIR等測試技術(shù)不同的是，通過研究分子動力學(xué)模擬技術(shù)運行后的有機蒙脫石模型，研究者可直接觀察到該體系層的間距變化、所含官能團種類、表面活性劑分子在有機蒙脫石內(nèi)部的排列狀態(tài)及其所導(dǎo)致的能量變化等現(xiàn)象。

雖然分子動力學(xué)模擬技術(shù)在研究微觀尺度下的有機蒙脫石時有很大優(yōu)勢，但在該研究領(lǐng)域仍存在一些不足，例如有些研究者沒有充分考慮到蒙脫石類質(zhì)同象取代位置的不同對有機蒙脫石的影響、建立的模型較為簡化等，因此提出如下建議：

(1)在模擬初始建模階段，應(yīng)該充分考慮蒙脫石層間環(huán)境條件、層間電荷來源等因素，以使建立的蒙脫石模型更接近實際。

(2)在有機蒙脫石的合成和應(yīng)用模擬研究中，應(yīng)該綜合考慮模擬體系的溫度、壓強、pH等外部環(huán)境因素的變化對有機蒙脫石層間結(jié)構(gòu)以及吸附離子賦存形式等的影響機制，這樣得到的模擬分析結(jié)果更具有實際意義。

相信隨著計算機運算能力的提高、模擬軟件的更新以及模擬方法的不斷完善，分子動力學(xué)模擬技術(shù)將在很大程度上促進有機蒙脫石的研究與應(yīng)用。