1、題目：計算機模擬MCM-4怫石分子篩吸附COf生能的研究完成日期：2013年6月計算機模擬MCM-41沸石分子篩吸附CO性能的研究摘要：MCM-41是M41S族中的典型代表，它具有六方有序的孔道結構，孔徑尺寸可隨合成時加入導向劑及合成件的不同在（1.510）nm之間變化，MCM-41孔徑均勻，具有較高的比表面積（1000m2/g）和大的吸附容量（0.7mL/g）,有利于有機分子的自由擴散。是優良的催化劑載體。這種新型的介孔分子篩具有穩定的骨架結構、孔道規則排列有序、孔徑分布窄等優點，填補了大分子催化材料的空缺。相較于老式的沸石分子篩材料，MCM-41介孔分子篩最大的優勢就在于：MCM-41介孔

2、分子篩的孔徑大大超出了常規沸石分子篩（孔徑小于1.5nm的孔徑范圍）。MCM-41被廣泛用于吸附、分離、催化等領域由于其理想的結構，可用于現有吸附理論的檢驗和發展。該課題的最大意義在于，可以將研究結果應用于汽車尾氣中CO去除，工業廢氣中CO的分離，煤礦坑道氣中CO的分離和去除以及CO氣體的濃縮等等。本文使用計算機，運用MaterialStudio4.3軟件建立了MCM-41沸石分子篩和CO分子模型，并采用巨正則蒙特卡羅（GCMC）分子模擬法模擬了在不同溫度和力場條件下，MCM-41沸石分子篩對CO的吸附量，并通過對數據的對比分析，得出了吸附量與溫度、壓力、力場之間的關系，為MCM-41沸石分子

3、篩在實際生產中的應用奠定了一定的基礎。關鍵詞：MCM-41沸石分子篩；巨正則蒙特卡羅（GCMC）;吸附；分子模擬ComputersimulationofMCM-41zeolitemolecularsieveadsorptionperformanceoftheCOAbstract:MCM-41isatypicalrepresentativeoftheM41Sfamily,Ithassixpartyorderedporestructure,MCM-41hasuniformaperture,highspecificsurfacearea(1000m2/g)andlargeradsorptioncap

4、acity(0.7ml/g),isconducivetothefreediffusionoforganicmolecules.Itsagoodcatalystcarrier.Thisnewtypeofmesoporousmolecularsievehasthestableskeletonstructure,channelrulesorderly,narrowporesizedistribution,etc;tocompensateforthelossofthemacromoleculecatalyticmaterials.Comparedwiththeoldmaterialofzeolitem

5、olecularsieves,MCM-41mesoporousmolecularsieveisthat:thebiggestadvantageofMCM-41mesoporousmolecularsieveapertureisgreatlybeyondtheconventionalzeolitemolecularsieves(poresizelessthanthediameterof1.5nmrange).Duetoitsidealstructure,MCM-41waswidelyusedinadsorption,separation,catalysisandotherfields.Inadd

6、ition,itcanbeusedfortheinspectionofexistingadsorptiontheoryanddevelopment.ThegreatestsenseofthesubjectisthattheresultscanbeusedintheautomobilecollectCO,COseparationinindustrialwastegas,thecoalminetunnelCOinthegasseparationandremovalandCOgasenrichment,andsoon.Inthispaper,usingthecomputerandtheMateria

7、lStudio4.3softwareestablishedtheMCM-41zeolitemolecularsieveandCOmoleculemodel,andusinggrandcanonicalmontecarlo(GCMC)molecularsimulationmethodissimulatedunderdifferentconditionsoftemperatureandforcefield,MCM-41zeolitemolecularsieveadsorptionofCO,andthroughthecomparisonandanalysisofdata,theadsorptionq

8、uantityisobtainedandtherelationshipbetweenthetemperature,pressure,forcefieldofMCM-41zeolitemolecularsieveintheactualproductionapplicationlaidacertainfoundation.Keywords:MCM-41zeolitemolecularsieve;GrandCanonicalEnsembleMonteCarlo(GCMC);Adsorption;Molecularsimulation一、緒論11.1 課題研究的背景11.2 課題研究的目的和意義21.

9、3 國內外相關課題的研究現狀31.4 選題依據及論文設想6二、MCM-41沸石分子篩82.1 MCM-41沸石分子篩的結構82.2 MCM-41沸石分子篩結構模型的建立92.3 MCM-41沸石分子篩的吸附原理122.4 MCM-41沸石分子篩的實驗室制備122.4.1 實驗儀器及藥品131 .實驗儀器132 .實驗藥品133 .4.2實驗部分132.5 MCM-41沸石分子篩的應用14三、計算軟件和理論計算方法173.1 計算軟件-MaterialStudio173.1.1 誕生背景173.1.2 軟件概況173.1.3 模塊簡介173.1.4 力場簡介193.2 理論計算方法-GCMC22

10、3.3 計算過程24四、結果與討論254.1 溫度對MCM-41沸石分子篩吸附性能的影響254.1.1 COMPASS力場下溫度對MCM-41沸石分子篩吸附性能的影響254.1.2 Universal力場下溫度對MCM-41沸石分子篩吸附性能的影響264.1.3 cvf力場下溫度對MCM-41沸石分子篩吸附性能的影響274.1.4 力場、溫度對MCM-41沸石分子篩吸附CO能力的影響284.2 壓力對MCM-41沸石分子篩吸附性能的影響294.2.1 COMPASS力場下，壓力對MCM-41沸石分子篩吸附性能的影響294.2.2 Universal力場下壓力對MCM-41沸石分子篩吸附性能的影

11、響304.2.3 cvff力場下壓力對MCM-41沸石分子篩吸附性能的影響304.2.4 力場、壓力對MCM-41沸石分子篩吸附性能的影響314.3 結論32參考文獻33致謝35一、緒論1.1 課題研究的背景一氧化碳（carbonmonoxide,CQ純品為無色、無臭、無刺激性的氣體。分子量28.01,密度1.250g/L,冰點為-207C,沸點-190C。在水中的溶解度甚低。空氣混合爆炸極限為12.5%74%。一氧化碳進入人體之后會和血液中的血紅蛋白結合，進而使血紅蛋白不能與氧氣結合，從而引起機體組織出現缺氧，導致人體窒息死亡。因此一氧化碳具有毒性。一氧化碳是無色、無臭、無味的氣體，故易于忽

12、略而致中毒。常見于家庭居室通風差的情況下，煤爐產生的煤氣或液化氣管道漏氣或工業生產煤氣以及礦井中的一氧化碳吸入而致中毒。一氧化碳的吸附濃縮便具有了一定的意義。1992年,Kresge等人在Nature雜志上首次報道了一種名為MCM-41的有序介孔材料，它是一種新型的納米結構材料，具有孔道呈六方有序排列、大小均勻、孔徑可在270nm范圍內連續調節、比表面積大等特點。MCM-41合成區別于傳統Molecule分子篩合成的最大特點是所用模板劑不同，傳統沸石或Molecule分子篩的合成是以單個有機小Molecule分子或金屬離子為模板劑，以ZSM-5為例，所用典型模板劑為四丙基胺離子，晶體是通過酸鹽

13、在模板劑周圍的縮聚形成的。而MCM-41的合成則不同，它是以大Molecule分子表面活性劑為模板劑，模板劑的烷基鏈一般多于10個碳原子，關于它的形成目前已提出兩種機理，而且仍不斷進行改進及完善。與其它沸石材料相比，MCM-41的骨架鋁物種熱穩定性相對較差，在焙燒過程中，骨架鋁物種由骨架脫落成為非骨架鋁物種。1992年，美國Mobil公司的研究者首次合成了MCM-41S1系列硅基介孔分子篩，揭開了分子篩應用研究的新紀元。這種新型的介孔分子篩具有穩定的骨架結構、孔道規則排列有序、孔徑分布窄等優點，填補了大分子催化材料的空缺。MCM-41是M41S族中的典型代表，它具有六方有序的孔道結構，孔徑尺寸

14、可隨合成時加入導向劑及合成件的不同在1.510nm之間變化，MCM-41孔徑均勻，具有較高的比表面積(1000m2/g)和大的吸附容量(0.7mL/g),有利于有機分子的自由擴散。是優良的催化劑載體。MCM-41被廣泛用于吸附、分離、催化等領域由于其理想的結構，可用于現有吸附理論的檢驗和發展。有關MCM-41的吸附及各種表面性質3的理論和實驗的研究很多如N2、。2、Ar、SQ、CO2、環戊烷、環己烷、苯和水等各種吸附質在MCM-41中的吸附性質的測量及其吸附機理的研究等2區別了孔內流體達到毛細凝聚的臨界溫度Tcp與迂回滯線的臨界溫度，探討了Tcp與孔徑和吸附質的關系，并研究了水在MCM-41中

15、的相態、冷凝和熔解溫度與孔徑的關系等3-4以及氧在MCM-41中的相變。1.2 課題研究的目的和意義由于具有可調的孔徑，因此人們預計MCM-41分子篩將在多相催化、離子交換、吸附分離以及高等無機材料等領域有較高的工程應用及學術價值4。MCM-41不僅廣泛用于吸附、分離、催化等領域，亦被用于吸附理論的驗證與改進方面。一氧化碳可以作為氣體燃料、一氧化碳作為還原劑，高溫時能將許多金屬氧化物還原成金屬單質，因此常用于金屬的冶煉。然而，一氧化碳進入人體之后會和血液中的血紅蛋白結合，進而排擠血紅蛋白與氧氣的結合，從而出現缺氧，造成一氧中毒。為了避免CO中毒事故以及回收利用CO,所以在某些場合我們要收集CO

16、,或者對CO進行增濃處理等等；于是將MCM-41沸石分子篩和CO聯系在了一起，利用MCM-41沸石分子篩對CO進行吸附，為了研究MCM-41沸石分子篩對CO的吸附性能，本文采用MaterialStudio軟件進行模擬5計算。該課題的最大意義在于，可以將研究結果應用于汽車尾氣中CO去除，工業廢氣中CO的分離，煤礦坑道氣中CO的分離和去除以及CO氣體的濃縮等等。1.3 國內外相關課題的研究現狀由于MCM-41沸石分子篩具有良好的特性、優異的吸附性能，關于MCM-41沸石分子篩的計算機模擬研究早在1992年就開始了。在理論研究方面，Ravikovitch等6利用密度泛函理論研究了N2在MCM-41中

17、的吸附等溫線并表征了MCM-41的孔徑分布。吳志強、劉志平7采用巨正則系綜MonteCarlo模擬了CO/H2在碳納米狹縫孔內的吸附和分離機理。研究結果表明：混合物中H2的吸附量高于與其分壓相同壓力下純H2的吸附量，而CO則與之相反；止匕外，還詳細研究了壓力、溫度和混合氣體組成對吸附量和平衡分離因子的影響，發現平衡分離因子隨壓力降低而提高，而低壓下尤其明顯。Maddox等8-9采用巨正則系綜MonteCarlo(GCMC)方法研究了純流體及二元混合物在MCM-41中的吸附性質。王巍10采用基于MonteCarlo取樣方法的分子力學模擬方法，研究了丙烯和4種C4單烯炫異構體在絲光沸石內的物理吸附

18、，得到了單組分烯炫在絲光沸石上的等壓吸附性能曲線，以及等比例多組分體系和不等比例多組分體系的等壓吸附曲線。侯廷軍11等采用分子動力學模擬的方法研究了純硅MCM-22型分子篩(ITQ-1)中苯分子的吸附行為。表明在較低吸附值的情況下，分子篩骨架的柔性對苯分子吸附和擴散并沒有產生大的影響。劉蓿、張現仁12提出了一個混合勢模型，利用巨正則系綜蒙特卡羅方法分別研究了一氧化碳、甲烷及乙烷在MCM-41中的吸附等溫線。本課題采用的研究方法是蒙特卡羅法，蒙特卡羅法可分為巨正則系綜蒙特卡羅法(GrandCanonicalEnsembleMonteCarlo,GCMC)和構型偏倚蒙特卡羅法(Configurat

19、ional-biasMonteCarlo,CBMC)。利用蒙特卡羅法可得到氣體小分子在催化劑中的吸附性質，包括吸附位信息、吸附量和吸附能。GCMC方法最早被用于主體流體的模擬，后來推廣到限定空間。巨正則系綜MonteCarlo方法廣泛用于研究吸附質的吸附相平衡的研究，它不僅可以模擬分子的靜態結構和動態行為（如氫鍵的締合與解締、吸附、擴散等），也可以模擬微孔中吸附質的化學勢、溫度保持恒定時的平衡及預測分子在微孔固體（如分子篩）中的吸附性質，并且可以計算吸附等溫線、結合位、吸附熱、擴散途徑及分子選擇性等。MonteCarlo模擬已廣泛應用于研究分子篩的吸附性能及分子篩內吸附質的動態分布。曾勇平等1

20、3采用GCMC方法模擬了嚷吩-苯二元組分和曦吩-苯-正己烷三元組分在MFI和MOR沸石中的吸附分離性能。孫曉巖等14采用巨正則MonteCarlo方法研究了不同溫度、不同吸附方式下純硅MCM-22型分子篩ITQ-1上苯與丙烯分子的吸附行為。呂玲紅等15用巨正則系綜MonteCarlo與構型偏倚相結合的方法模擬了烷炫、異構烷炫及其多組分混合物在分子篩上的吸附與分離行為，得到了不同混合物的吸附選擇性。張舟等16用巨正則系綜蒙特卡羅及構型偏移蒙特卡羅的方法模擬了300K時，CiC7直鏈烷炷分子在MFI沸石分子篩中的吸附行為，并用目前常用的幾種吸附理論（單雙點Langmuir吸附等溫式、擬化學近似吸附

21、等溫式及Gibbs吸附等溫式）對之進行擬和。在分子篩的吸附擴散研究中，分子模擬技術已成為必不可少的研究手段。其中，MonteCarlo方法和分子動力學模擬的組合方法更是得到了廣泛的應用，并取得良好的效果。1.4 選題依據及論文設想通過查閱文獻發現，研究計算機模擬MCM-41沸石分子篩吸附CO性能的文獻很少，這給我的研究帶來了些許的困難，但通過對前人研究成果的閱讀和思考，初步形成了自己的研究思路：首先運用MaterialStudio軟件建立MCM-41沸石分子篩，使用Sorption模塊進行GCMC（巨正則系綜MonteCarlo）模擬；GCMC模擬的對象是一溫度T、體積V、化學位M亙定的開放體

22、系，而這些參數在實驗室吸附測定過程中也是恒定的，因此模擬結果可以與實驗結果直接進行比較。采用GCMC模擬時，只有固相被模擬，氣相被當作一個假想的溫度和化學勢已知的粒子庫。模擬過程中，分子篩微孔（模擬盒子）內的氣體分子與體相流體（粒子庫）中的氣體分子之間進行不斷的質量和能量交換，只要成功交換的氣體分子數足夠多，孔內與體相氣體的化學勢就會相等，平均吸附量就可以直接由計算得到。模擬過程中，吸附劑（MCM-41）與吸附質（CO）、吸附質與吸附質之間相互作用時，范德華相互作用截斷距離為1.25nm,靜電相互作用采用Ewald加和法處理，相互作用參數來自于COMPASS力場，采用周期性邊界條件。在每次GC

23、MC模擬過程中，吸附劑的每一個構型根據能量變化運用Metropolis運算規則接受或拒絕，其中分子交換被接受概率為39%、分子構象異構化被接受概率為20%、分子轉動被接受概率為20%、分子平動被接受概率為20%。為使體系達到真正平衡，計算進行了2X105個MonteCarlo步，其中前1M05個MonteCarlo步用來達到平衡，后1M05個MonteCarlo步用來統計平均。整個模擬計算采用AccelrysInc.的MaterialStudio4.3軟件完成。通過計算機模擬得出不同溫度和壓力下MCM-41沸石分子篩對CO吸附量的具體數值；最后對數值進行分析，得出最佳吸附溫度和最佳吸附壓力。2

24、2二、MCM-41沸石分子篩2.1 MCM-41沸石分子篩的結構MCM-41是M41S族中的典型代表，它具有六方有序的孔道結構，孔徑尺寸可隨合成時加入導向劑及合成件的不同在（1.510）nm之間變化，MCM-41孔徑均勻。是優良的催化劑載體。孔徑約3.5nm,晶格參數約為4.5nm,壁厚約1nm;壁表面積約1000m2,g-1;孔體積約為1mL.g-1;顆粒具有規則外形，且可在微米尺度內保持高度的孔道有序性。圖2.1MCM-41沸石分子篩MCM-41分子篩有幾個顯著特點：（1）與其它介孔材料相比，孔徑分布狹窄；（2）孔徑大小可通過改變表面活性劑的鏈長來調節；（3）具有較高的熱穩定性和水熱穩定性

25、；（4）孔道排列有序。其中MCM-41分子篩的孔道呈六方有序排列，是一致的平行軌道，穩定性高，因而引起人們的更多關注。由于微孔材料和介孔材料具有較大的內表面積，在作為催化劑和吸附劑上有了相當大的應用，但包括一些晶態的沸石在內，目前已知最大孔徑不超過1.4nm,如一些金屬磷酸鹽（1.01.2nm）和黃硫鐵礦（1.4nm）。作為一種理想的催化劑，要求其可逆吸附量大、孔徑分布窄、催化活性高、疏水性好、水熱穩定性好，所以就要尋求更好的分子篩材料來滿足現代工業的需求。與其它介孔材料相比，MCM-41分子篩是一種性能極為優良的分子篩，它的出現給分子篩領域帶來新活力。MCMW1分子篩屬于一維孔道體系結構，其

26、孔徑均勻，具有高比表面積（1000m2/g）和大吸附容量（0.7mL/g）的特點，比沸石和磷鋁酸鹽等微孔材料更有利于有機分子的快速擴散，這使得它能為大分子尤其是石油化工過程中重油有機分子進行擇型反應提供無可比擬的有利空間和有效酸性活性中心，可根據需要調節孔徑和酸性濃度、強度，在渣油催化裂化、重油加氫、潤滑油加氫、烷基化、烯炫聚合、CO2-CH4的分離等酸催化領域和石油化工的分離過程中具有相當大的潛在價值。2.2 MCM-41沸石分子篩結構模型的建立MaterialStudio里面自帶了很多分子篩，但是像MCM-41這樣的分子篩,MaterialStudio里面沒有，那就得自己建了，以下是建立M

27、CM-41分子篩的主要步驟：1 .建立MCM-41沸石分子篩模型1. FileImportMaterialstudio4.3/share/structures/glasses/SiO221A.3d.msio2. BuildSymmetrySupercell輸入A,B,C數值（比如說3*7*2）按一CreatesupercellII3. EditAtomSelectionSelectbyproperty在下拉菜單中選擇一RadialDistanceIIWithin-10A（孔的半徑）-Relativeto-location7輸入坐標（如21.42,21.42,0）,Selectionmode選擇

28、第一個select右鍵delete每次在z方向推進2，這樣在慢慢推進后就會挖出一個圓孔)，無定形SiO2的單胞的邊長是一定的，那么此時3*7*2的超晶胞的各個邊的邊長也是確定的，超晶胞的邊長分別為64.18458,149.76402,42.78972,按照幾何關系，可以計算出每個孔道的圓心位置，中間圓心開始的位置分別是：(32.0923,18.7205,0)(32.092,56.1615,0)(32.0922,93.6025,0)、(32.0923,131.0435,0);左邊上半圓以及四分之一圓孔開始的圓心位置為(0,0,0)、(0,37.44105,0)(0,74.88201,0)、(0,

29、112.323015,0)、(0,149.76402,0);右邊的半圓和四分之一圓孔的開始圓心位置也可以按照此種方法寫出。4. 點擊一AdjustHydrogen的(H)后的Si補齊電荷。5. EditAtomSelectionSeleCtyproperty下拉菜單中選擇Element7選擇SiSelectionmode下選擇第一個點擊SelectSelectbyproperty下拉菜單中選擇connectedSelectionmode選擇第一個點擊SelectSeleclbyproperty下拉菜單中選擇element選擇HHSelectionmode下選擇第二個-點擊Select7Modi

30、fyModifyelementoxygen點擊一AdjusHydrogenI典做這一步的原因是Si不能直接連著H,只能連著O。6.EditFindpatterns找出孔壁上連有3個-OH的Si,將該結構刪除,再補足H原子，完成后的MCM-41結構如圖1所示：圖1MCM-4的四孔模型2.使用Dmol3模塊量化計算簇模型的ESP電荷得到了MCM-41的結構，下一步就要進行量化計算電荷，才能進行GCMC和MD模擬。計算電荷時需要在MCM-41中切取簇模型才行，要切取簇模型的話，就需要破壞骨架，Buildsymmetryfnon-periodicstructure,然后你就可以隨意的刪除原子了，得到你

31、所期望的簇模型了；從基本道理上來說，簇模型是隨機選擇的，但是為了計算的準確性，最好切取多個簇模型，然后計算各自簇模型中不同原子的電荷，再次求各簇模型的特定原子的電荷平均；因此可以切取如圖2所示的那樣的簇模型，然后使用Dmol3模塊計算電荷。圖2MCM-41中所切取的簇模型對于某個特定的簇模型，Dmol3中有三種電荷表示方式，在接下來的MC和MD計算中，ESP電荷描述方式是最好的，因此在開始計算的時候，一定要勾選ModuleDmolACalculationPropertiesPopulationanalysis（前面打勾）fESPcharges面打勾），其他的設置按照文獻來，然后點擊一run1鈕

32、，開始量化計算；計算結束后，在Dmol3的結果文件（左邊Project中擴展名為.outmol的文件）中找到ESP電荷,使用Excel文件分列和排序，求解不同原子的電荷數值平均，然后再把不同簇模型所得到的某個特定原子的平均電荷再次求平均，這個平均電荷將要進行下一步的GCMC計算。現在，模型也搭好了，電荷也計算好了，下面要做的就是平衡電荷，使用的是MaterialStudio中的Sorption模塊，在MCM-41的四孔模型中，點擊某種原子，按住Alt鍵，鼠標雙擊，在Propertiesexplorer里面找到charge那一項，輸入剛才量化計算得到的那個原子的電荷，每個原子的電荷賦予過程都是一

33、樣的，然后就得到了一個真實的帶電體系；但是此時由于簇模型和整體模型的電荷差異，也許此時整體模型有少許的電荷過剩，在Modifycharges里面關注整體四孔模型的帶電剩余，把此電荷剩余除以四孔模型整體的原子數量（這個在propertiesFilterSymmetsystem里面可以看到）,得到一個數值，然后把四孔模型中每個種類的原子電荷加上或者減去這個數值，四孔模型的電荷就平衡了。2.3 MCM-41沸石分子篩的吸附原理。2、N2、H2和CO等高純度氣體的生產方法有很多種，分子篩以其具有高吸附選擇性和高吸附容量的特性在這些工藝中起關鍵作用。分子篩在深冷生產高純度也和CO工藝中是不可缺少的一部分

34、。在典型的深冷法工藝中，由天然氣重整得到的合成氣（主要成分為CO和HJ首先通過一個吸附床，除去大部分的CO2,隨后，飽含水汽的合成氣再通過分子篩吸附床，吸水和CO2的含量脫除到遠小于1ppm的水平；保證了在以后的工序的分離過程中，管道不會因水和CO2的結冰而堵塞。使用高吸附容量和長壽命的分子篩吸附劑，還有助于降低設備投資和操作費用。2.4 MCM-41沸石分子篩的實驗室制備目前合成中孔分子篩主要采用水熱合成法。室溫合成16,17、微波合成18、濕膠焙燒法19、相轉變法20及在非水體系21中合成也有少量報道。水熱合成的一般過程為：（1）生成比較柔順、松散的表面活性劑和無機物種的復合產物；（2）水

35、熱處理提高無機物種的縮聚程度，提高復合產物結構的穩定性；（3）焙燒或溶劑抽提除掉復合產物中的表面活性劑后得到類似液晶結構的無機多孔骨架，即中孔分子篩。2.4.1 實驗儀器及藥品1 .實驗儀器DF-101S集熱恒溫加熱磁力攪拌加熱器、電子天平AL-204、SX2-12-16馬氟爐、DGB/20-002臺式干燥箱、TSA-900型原子吸收分光光度計、800B離心機、SHD-2000型循環水多用真空泵2 .實驗藥品十六烷基三甲基澳化氨（CTAB）、正硅酸乙酯（TEOS）、氨水、無水乙醇、氫氧化鈉（NaOH）、硫酸鎘2.4.2 實驗部分方案一按照TEOS（正硅酸乙酯）：NH3H2O:CTAB:H2O=

36、1:1.64:0.15:126的摩爾比稱取十六烷基三甲基澳化鏤2.4011g,置于250mL的燒杯中，加入120mL去離子水，在30c恒溫下攪拌至澄清，加入8mLNH3H2O（32%）攪拌5min,再緩慢滴加10mL正硅酸乙酯，在30c恒溫下攪拌12h,減壓過濾，用去離子水洗滌5次，再用乙醇洗1次,室溫干燥得到MCM-41原粉。將原粉放入馬弗爐中在550c下，焙燒6h,除去表面活性劑，最后得到中孔孔分子篩MCM-41。方案二按照TEOS:CTAB:NaOH:H2O=1：0.6:0.24：100的摩爾比稱取十六烷基三甲基澳化鏤和氫氧化鈉。將氫氧化鈉溶于去離子水中，控制水溫在36c左右，緩慢攪拌下

37、加入模板劑十六烷基三甲基澳化鏤(CTAB)至溶液澄清；待溶液澄清后，在劇烈攪拌下緩慢滴入正硅酸乙酯，繼續攪拌2h。接著用0.3molL-1的HCl溶液調節pH值至(J910,將反應后所得懸濁液裝入聚四氟乙烯襯墊的不銹鋼自升壓反應釜中，然后置于110c的烘箱中靜置晶化3天，取出產物抽濾、洗滌，在80c下干燥得到MCM-41原粉。將原粉放入馬弗爐中在550c下，焙燒6h,除去表面活性劑，最后得到中孔孔分子篩MCM-41。2.5 MCM-41沸石分子篩的應用MCM-41的有序介孔材料，它是一種新型的納米結構材料。與其它沸石材料相比，MCM-41的骨架鋁物種熱穩定性相對較差，在焙燒過程中，骨架鋁物種由

38、骨架脫落成為非骨架鋁物種。它作為一種特殊的介孔分子篩，在化工領域、環境治理和保護領域、生物醫藥領域等眾多領域，均有重要應用。而MCM-41的負載產品則具有更多新的作用。M41S具有有序孔道排列和狹窄的孔分布、大比表面積和高熱穩定性，是一種良好的催化劑載體，但硅基介孔材料自身無催化活性，因此需要采取摻雜、改性等多種方法引入催化活性中心。1.環境保護領域TiO2-MCM-41降解有機污染物高級氧化對徹底廣泛的降解有機污染物和細菌類物質是一個有效的措施。而TiO2作為催化劑來高級氧化氣體污染物和水污染物是近幾年來所有已發表的化工文獻中的最有效的措施。但是TiO2本身對于催化反應來講又有些局限性。它不

39、具備孔道，所以呈現出較低的吸附性能，尤其是對于非極性的化合物來講尤其明顯。無孔道，降低了它的吸附能力，然而，吸附步驟是對于反應來說極其重要的一個前期步驟。MCM-41的孔徑在1.872.25nm之間，將TiO2負載到MCM-41上，較之裸露的TiO2極大提高了該催化劑的反應活性，因為負載之后的TiO2的比表面積得到了提高，從而使有機污染物在催化劑表面聚集，同時催化劑表面的電子空穴得以減少。與此相關的實驗已取得一些成果，雖然實際效果遠非理想，但科研人員仍在堅持不懈地鉆研。2.材料開發領域一一對氧氣敏感的光致激發光材料此種材料對一種新型的銅復合材料Cu(POP)(PTI)BF4,此材料由銅化合物與

40、無機多孔分子篩組成一一在一種特定溶劑中溶解銅化合物，然后加入一定量MCM-41分子篩制成溶膠，在室溫下攪拌、過濾，然后用溶劑洗滌三次，在空氣中風干。這種復合物感光材料色白，在紫外燈照射下的最大波長為500nm,與目前廣被使用的高價惰性金屬復合材料相比，此感光符合材料造價低，制作簡單，敏感性高，毒性低，非常適用于氧探測器。3 .自身性狀改性純硅MCM-41具有典型的六方結構，具有很高的比表面積，接枝改性(通過嫁接新的基團，使原來的物質具有新的性質)獲得的Al-MCM-41能夠保持較好的規整性和很高的比表面積，同時獲得了一定的酸性；水熱合成的Al-MCM-41結晶度稍差，但酸量要比同硅鋁比的接枝改

41、性的Al-MCM-41要高出很多。4 .有機氧化領域(1)含有雜原子的多孔的MCM-41沸石分子篩作為雙氧水氧化環己烷的催化劑。Cr-MCM-41展示出了最好的催化行為。環己烷的轉化率和雙氧水的利用率都受到反應溫度、雙氧水濃度、雙氧水添加方式的影響。環己烷的轉化率是33.7%,反應對環己酮和環己醇的選擇性是97.2%,環己酮與環己酉I的分子比率是1.99。當溫度處于70OC,環己烷與雙氧水的摩爾比是1:2并且分兩次加入雙氧水時，雙氧水的轉化率是37.1%。(2)Al-MCM-41負載到鉗催化劑上，用以催化甲烷的燃燒。這項工作目的是闡明處理過的鉗催化劑在甲烷燃燒過程中的表現行為。將鉗負載到有序的

42、介孔分子篩Al-MCM-41(Si:Ai摩爾比75:1),作為分子支架。與沒有處理過的催化劑相比，Pa/Al-MCM-41顯示出較高的活性，和較好的穩定性。三、計算軟件和理論計算方法3.1 計算軟件-MaterialStudio3.1.1 誕生背景美國Accelrys公司的前身為四家世界領先的科學軟件公司一一美國MolecularSimulationsInc.(MSI)公司、GeneticsComputerGroup(GCG)公司、英國SynopsysScientific系統公司以及OxfordMolecularGroup(OMG)公司，由這四家軟件公司于2001年6月1日合并組建的Accel

43、rys公司，是目前全球范圍內唯一能夠提供分子模擬、材料設計以及化學信息學和生物信息學全面解決方案和相關服務的軟件供應商。Accelrys材料科學軟件產品提供了全面完善的模擬環境，可以幫助研究者構建、顯示和分析分子、固體及表面的結構模型，并研究、預測材料的相關性質。Accelrys的軟件是高度模塊化的集成產品，用戶可以自由定制、購買自己的軟件系統，以滿足研究工作的不同需要。Accelrys軟件用于材料科學研究的主要產品包括運行于UNIX工作站系統上的Cerius2軟件，以及全新開發的基于PC平臺的MaterialsStudio軟件。Accelrys材料科學軟件被廣泛應用于石化、化工、制藥、食品、

44、石油、電子、汽車和航空航天等工業及教育研究部門，在上述領域中具有較大影響的世界各主要跨國公司及著名研究機構幾乎都是Accelrys產品的用戶。3.1.2 軟件概況MaterialsStudio是專門為材料科學領域研究者開發的一款可運行在PC上的模擬軟件。它可以幫助我們解決當今化學、材料工業中的一系列重要問題。支持Windows98、2000、NT、Unix以及Linux等多種操作平臺的MaterialsStudio使化學及材料科學的研究者們能更方便地建立三維結構模型，并對各種晶體、無定型以及高分子材料的性質及相關過程進行深入的研究。多種先進算法的綜合應用使MaterialsStudio成為一個

45、強有力的模擬工具。無論構型優化、性質預測和X射線衍射分析，以及復雜的動力學模擬和量子力學計算，我們都可以通過一些簡單易學的操作來得到切實可靠的數據。MaterialsStudio軟件采用靈活的Client-Server結構。其核心模塊Visualizer運行于客戶端PC,支持的操作系統包括Windows98、2000、NT；計算模塊(如Discover,Amorphous,DMol3,CASTEP等)運行于服務器端，支持的系統包括Windows2000、NT、SGIIRIX以及RedHatLinux。浮動許可(FloatingLicense)機制允許用戶將計算作業提交到網絡上的任何一臺服務器上

46、，并將結果返回到客戶端進行分析，從而最大限度地利用了網絡資源。MaterialsStudio軟件能使研究者達到與世界一流研究部門相一致的材料模擬的能力。模擬的內容包括了催化劑、聚合物、固體及表面、晶體與衍射、化學反應等材料和化學研究領域的主要課題。3.1.3 模塊簡介MaterialsStudio采用了大家非常熟悉的Microsoft標準用戶界面，允許用戶通過各種控制面板直接對計算參數和計算結果進行設置和分析。目前，MaterialsStudio軟件包括如下功能模塊：MaterialsVisualizer:提供了搭建分子、晶體及高分子材料結構模型所需要的所有工具，可以操作、觀察及分析結構模型，

47、處理圖表、表格或文本等形式的數據，并提供軟件的基本環境和分析工具以及支持MaterialsStudio的其他產品。是MaterialsStudio產品系列的核心模塊。Discover:MaterialsStudio的分子力學計算引擎。使用多種分子力學和動力學方法，以仔細推導的力場作為基礎，可準確地計算出最低能量構型、分子體系的結構和動力學軌跡等。COMPASS:支持對凝聚態材料進行原子水平模擬的功能強大的力場。是第一個由凝聚態性質以及孤立分子的各種從頭算和經驗數據等參數化并經驗證的從頭算力場。可以在很大的溫度、壓力范圍內精確地預測孤立體系或凝聚態體系中各種分子的結構、構象、振動以及熱物理性質。

48、DMo獨特的密度泛函（DFT）量子力學程序，是唯一的可以模擬氣相、溶液、表面及固體等過程及性質的商業化量子力學程序，應用于化學、材料、化工、固體物理等許多領域。可用于研究均相催化、多相催化、分子反應、分子結構等，也可預測溶解度、蒸氣壓、配分函數、熔解熱、混合熱等性質。CASTEP:先進的量子力學程序，廣泛應用于陶瓷、半導體、金屬等多種材料，可研究：晶體材料的性質（半導體、陶瓷、金屬、分子篩等）、表面和表面重構的性質、表面化學、電子結構（能帶及態密度）、晶體的光學性質、點缺陷性質（如空位、間隙或取代摻雜）、擴展缺陷（晶粒間界、位錯）、體系的三維電荷密度及波函數等。3.1.4 力場簡介分子模擬的基

49、礎，是準確計算原子之間的相互作用，包括組成同一分子的原子之間的成鍵相互作用，和不同分子間的范德華相互作用，有的分子間還有氫鍵相互作用。描述原子間的這些相互作用，有兩種方式，一個是通過量子化學計算，另外一種方式就是采用分子力場計算。1、分子力場的來源我們知道，量子化學計算分子結構和原子、分子間相互作用比較準確，但是很慢；而采用分子力場計算就會很快，因為分子力場并不計算電子相互作用，它是對分子結構的一種簡化模型，所以計算很快。在這個模型中，它把組成分子的原子看成是由彈簧連接起來的球，然后用簡單的數學函數來描述球與球之間的相互作用。比如，氫分子，看做有彈簧鏈接的兩個球的話，可以用胡克定律描述兩個氫原

50、子間的能量：E=k(b-b0)2。其中，b表木兩氫原子間距離，b0表不平衡時原子間距，k為鍵能系數，bo和K稱為力場參數。更復雜一點可以用四次方表達：E=Ki(b-bo)2+K2(b-b0)3+K3(b-b0)4,更多的參數可以獲得對成鍵分子的更精確的描述。這是描述成鍵作用，不成鍵的原子間的相互作用則采用Legendre-Jones®數，或者Bukingham函數描述。從上面可以看出來，力場用簡單的數學函數描述原子間作用，稱為分子力場，又叫分子力學力場。采用分子力場的分子模擬稱為經典分子模擬。這是相對于采用量子力學計算的分子模擬來說的。分子力學對分子結構和原子間相互作用描述的是否準確

51、，主要依賴于所用的參數。而這些參數通常擬合自實驗數據，或者量子化學結果。它屬于經驗描述，顯然品質要低一些，但是由于計算速度快，適合于描述上千個乃至百萬個原子的模擬，在這些情況下，我們無法采用量子力學計算，因此，只能采用經典模擬。2、不同分子力場間的區別分子力場有很多，材料領域常用的CFF,MMFF,COMPASS等。一個力場通常包括三個部分：原子類型，勢函數，和力場參數。也就是說不同的力場，他們的函數形式可能不一樣，或者函數形式一樣而力場參數不一樣。其中，最關鍵的差別取決于分子力學模型，比如有的力場考慮氫鍵，有氫鍵函數；有的考慮極化，有極化函數。其次，分子力場參數都是擬合特定分子的數據而生成的

52、，比如，面向生物模擬的力場選擇生物領域的分子模擬得到參數，而材料的，則側重選擇材料方面的分子。這些被擬合的分子成為訓練基。1)COMPASS力場COMPASS力場是第一個把以往分別處理的有機分子體系的力場與無機分子體系的力場統一的分子力場。COMPASS力場能夠模擬小分子與高分子，一些金屬離子、金屬氧化物與金屬。在處理有機與無機體系時，采用分類別處理的方式，不同的體系采用不同的模型，即使對于兩類體系的混合，仍然能夠采用合理的模型描述。2)cvff力場cvff力場全名為一致性價力場(consistantvalenceforcefield),最初以生化分子為主，適應于計算氨基酸、水及含各種官能團的

53、分子體系。其后，經過不斷的強化，cvff力場可適用于計算多肽、蛋白質與大量的有機分子。此力場以計算系統的結構與結合能最為準確，亦可提供合理的構型能與振動頻率。3 )Dreiding力場Dreiding力場是Mayo、Olafson等人于1990年開發的一種較好的、具有各種用途的力場，它的最大優點在于有很強的預測能力，相對于那些為十分有限的體系提供較高精確度的特殊力場，Dreiding力場允許合理地預測大量的體系，包括那些含有新元素化合類型的體系，以及沒有或很少實驗數據的體系，也能較好地進行能量、結構的預測和動力學計算。Dreiding力場在預測依賴分子內相互作用的聚集態材料的性質比UFF力場更

54、好，但它未對整個周期表進行參數化，因而不及UFF力場應用廣泛。Dreiding力場沒有像UFF力場和ESFF力場那樣的-參數產生器U,但它的參數也是由一種基于規則的方法而產生。Dreiding是目前最新的版本，也是高分子材料模擬中最常用的力場。4 )UFF力場UFF(UniversalForceField)力場是由Rappe及其研究小組在1992年開發的優秀通用力場它對整個元素周期表進行了參數化，所有力場參數是由一系列基于元素、雜化及化合連接的規則產生的，它是適用范圍最廣泛的計算方法，對元素周期表內的所有元素適用，因此一般結構都可以采用這個算法。UFF力場被證明對許多體系有效，因而被推薦用來預

55、測那些未被特殊力場覆蓋的體系。通常UFF力場對分子結構預測比Dreiding力場好，但對依賴于分子內相互作用的性質預測則差一些。UFF力場另一特點是有一個整數產生器U,它能結合原子的參數產生力場參數，因此力場所需的任何原子化合類型參數都能被產生出來。UFF力場已發展到與電荷平衡法相結合，用UFF力場進行的計算極力主張使用電荷平衡計算，并為許多體系所證實。3、分子力場的局限由于力場參數是擬合訓練基分子得到的，那么這些參數用于計算其它分子準確性降低，這叫分子力場的遷移性問題。遷移性問題還包括狀態遷移性問題，就是說所擬合的實驗數據是常溫常壓下測量的，然后你模擬的可能是高溫高壓下的，那么分子力學的準確

56、性也會降低。這些都屬于分子力場的局限性。5 .2理論計算方法-GCMCGCMC模擬的對象是一溫度T、體積V、化學位區恒定的開放體系，而這些參數在實驗室吸附測定過程中也是恒定的，因此模擬結果可以與實驗結果直接進行比較。采用GCMC模擬時，只有固相被模擬，氣相被當作一個假想的溫度和化學勢已知的粒子庫22。模擬過程中，分子篩微孔（模擬盒子）內的氣體分子與體相流體（粒子庫）中的氣體分子之間進行不斷的質量和能量交換，只要成功交換的氣體分子數足夠多，孔內與體相氣體的化學勢就會相等，平均吸附量就可以直接由計算得到。模擬過程中，吸附劑（MCM-41）與吸附質（CH,、沖或CO?）、吸附質與吸附質之間相互作用時，范德華相互作用截斷距離為1.25nm,靜電相互作用采用Ewald加和法處理，相互作用參數來自于COMPASS力場，采用周期性邊界條件。在每次GCMC模擬過程中，吸附劑的每一個構型根據能量變化運用Metropolis運算規則接受或拒絕，其中分子交換被接受概率為39%、分子構象異構化被接受概率為20%、分子轉動被接受概率為20%、分子平動被接受概率為20%。為使體系達到真正平衡，計算進行了2M05個MonteCarlo步，其中前1X105個MonteCarlo步用來達到平衡，后1X105個MonteCarlo步用來統計平均。整個模擬計算采用AccelrysInc.的MaterialStudio