湯夢(mèng)瑤，孫天一，李佳書(shū)，丁延，李進(jìn)龍

石墨烯與碳納米管氫吸附性能的分子模擬

湯夢(mèng)瑤，孫天一，李佳書(shū)，丁延，李進(jìn)龍

（常州大學(xué) 石油化工學(xué)院，江蘇 常州 213164）

在293.15 K、101.30 kPa的條件下，測(cè)定了氫氣在石墨烯、單壁碳納米管和多壁碳納米管中的儲(chǔ)氫密度，對(duì)比了同等條件下采用不同力場(chǎng)計(jì)算得到的數(shù)據(jù)，篩選了3種碳材料的最佳計(jì)算力場(chǎng)。在此基礎(chǔ)上，進(jìn)一步計(jì)算了3種碳材料在0～1 000.00 kPa、77.00～573.15 K條件下的儲(chǔ)氫密度。結(jié)果表明，Dreiding力場(chǎng)是計(jì)算石墨烯吸附儲(chǔ)氫密度的最佳力場(chǎng)，Universal力場(chǎng)是計(jì)算碳納米管吸附儲(chǔ)氫密度的最佳力場(chǎng)；在給定條件下，3種材料吸附儲(chǔ)氫能力強(qiáng)弱排序?yàn)槭?單壁碳納米管>多壁碳納米管，儲(chǔ)氫能力與材料的比表面積及其與氫氣的弱結(jié)合力緊密相關(guān)。該研究結(jié)果可為分子模擬碳材料吸附儲(chǔ)氫和儲(chǔ)氫材料設(shè)計(jì)提供數(shù)據(jù)和理論支撐。

石墨烯； 碳納米管； 吸附儲(chǔ)氫性能； 分子模擬

氫氣具有單位質(zhì)量能量密度大、環(huán)境友好、持續(xù)再生和大規(guī)模存儲(chǔ)等優(yōu)點(diǎn)，是理想的清潔能源。碳材料具有比表面積大、強(qiáng)度高、韌性好、制備難度低等優(yōu)點(diǎn)，因此在儲(chǔ)氫材料領(lǐng)域引起了研究者和工業(yè)界的廣泛興趣，并已通過(guò)實(shí)驗(yàn)和計(jì)算機(jī)模擬等方法開(kāi)展了深入研究[1?6]。

對(duì)碳納米管材料，A.C.Dillon等[7]曾報(bào)道單壁碳納米管在室溫條件下氫氣吸附量為10%；E.E.Dogan等[8]通過(guò)實(shí)驗(yàn)證明，碳納米管比在無(wú)超聲波和有超聲波介質(zhì)中合成的活性炭?jī)?chǔ)氫能力分別高300%和265%；Y.Li等[9]和I.K.Petrushenko等[10]采用計(jì)算模擬方法研究了氫氣在碳納米管中的吸附和傳輸機(jī)理，結(jié)果表明材料的氫吸附性能與材料結(jié)構(gòu)、分子作用密切相關(guān)，如過(guò)大的層間距會(huì)削弱碳?xì)涞南嗷プ饔?，從而?dǎo)致吸氫能力降低。對(duì)石墨烯材料，已有實(shí)驗(yàn)[11?14]表明在一定溫度和壓力下儲(chǔ)氫量可達(dá)3.1%，但總體氫吸附性能遜色于碳納米管。計(jì)算機(jī)模擬表明，石墨烯儲(chǔ)氫性能也與材料結(jié)構(gòu)、修飾基團(tuán)種類(lèi)以及吸附條件等密切相關(guān)，且低壓下氫氣在石墨烯上的吸附行為符合Freundlich吸附模型。近年來(lái)，為提升碳納米管和石墨烯的儲(chǔ)氫性能，對(duì)其進(jìn)行了各種摻雜或改性研究，如在材料表面或內(nèi)部嫁接Li[15?16]、Ti[17?18]、Ni[19?20]等金屬原子。

為獲得氫氣在碳基材料中的吸附機(jī)理并指導(dǎo)材料的設(shè)計(jì)和制備，基于不同力場(chǎng)的計(jì)算機(jī)模擬方法被廣泛用于探索氫氣與碳基材料的微觀作用及形成的微觀結(jié)構(gòu)，如密度泛函理論(DFT)[21]、蒙特卡羅方法(MC)[22]和分子動(dòng)力學(xué)(MD)[23]等，這些計(jì)算機(jī)模擬方法在碳納米管和石墨烯氫吸附研究中被廣泛應(yīng)用[24?26]。然而，因力場(chǎng)和求解方法的不同，計(jì)算機(jī)模擬方法存在計(jì)算方法多、計(jì)算標(biāo)準(zhǔn)無(wú)法統(tǒng)一、模型普適性待提高等不足，如多數(shù)DFT方法選擇有效描述色散作用[27]、MC體系樣本空間取樣量大等。因此，針對(duì)具體研究系統(tǒng)，需要根據(jù)其特性獲得合適的模擬方法，進(jìn)一步研究系統(tǒng)的微觀作用機(jī)理。

本研究以碳納米管和石墨烯吸附氫為研究對(duì)象，通過(guò)實(shí)驗(yàn)測(cè)定一定條件下不同碳基材料的氫飽和吸附量，再選擇不同力場(chǎng)模擬實(shí)驗(yàn)條件下不同材料的氫吸附性能，從而獲得每種碳材料吸氫分子模擬的最優(yōu)力場(chǎng)，并據(jù)此進(jìn)一步模擬不同溫度和壓力下氫氣在碳材料中的儲(chǔ)氫密度和吸附熱，探索氫氣在不同碳基材料中的吸附機(jī)理，為碳基儲(chǔ)氫材料的設(shè)計(jì)和制備提供參考。

1 實(shí)驗(yàn)部分

1.1 材料與儀器

石墨烯納米片(Graphene)，片徑為5～10 μm，厚度為3～10 nm，比表面積為31.66 m2/g，純度大于99.50%；單壁碳納米管(SWNT)，直徑為1～2 nm，長(zhǎng)度為5～30 μm，比表面積大于380.00 m2/g，純度大于95.00%；多壁碳納米管(MWNT)，直徑為5～15 nm，內(nèi)徑為2～5 nm，比表面積大于200.00 m2/g，純度大于95.00%。本實(shí)驗(yàn)相關(guān)碳材料均購(gòu)自南京先豐納米材料科技有限公司。

高性能吸附分析儀ASAP 2020 Plus，麥克默瑞提克(上海)儀器有限公司；真空管式爐BTF?1500C，安徽意貝克設(shè)備技術(shù)有限公司；電子分析天平ME2002(精度0.000 1 g)，梅特勒?托利多儀器有限公司。

1.2 實(shí)驗(yàn)方法

碳材料的吸附儲(chǔ)氫性能使用ASAP 2020 Plus高性能吸附分析儀進(jìn)行測(cè)量。將樣品放入儀器前需要對(duì)其進(jìn)行研磨處理，粉碎的顆粒需盡可能小，在真空和150 ℃的條件下進(jìn)行預(yù)處理，將樣品中的殘氣脫除，處理時(shí)間為6 h；然后，在20 ℃、0～101.30 kPa的條件下進(jìn)行氫氣吸附測(cè)試，獲得氫氣的吸附等溫線。

1.3 模擬方法

儲(chǔ)氫吸附的模擬計(jì)算基于巨正則蒙特卡洛法(GCMC)。為了確定3種碳材料合適的計(jì)算力場(chǎng)參數(shù)，選用了COMPASS、COMPASSⅡ、Dreiding、Universal四種不同力場(chǎng)，利用Fix Pressure和Adsorption Isotherm計(jì)算3種碳材料中的氫氣吸附量，以實(shí)驗(yàn)值為基準(zhǔn)，篩選獲得對(duì)應(yīng)碳材料的最優(yōu)吸附儲(chǔ)氫計(jì)算力場(chǎng)，由此進(jìn)一步研究不同材料的吸附儲(chǔ)氫性能。計(jì)算中，采用Metropolis抽樣方法；靜電力計(jì)算借助Ewald求和法；范德華力計(jì)算采用原子截?cái)喾?，截?cái)喟霃皆O(shè)為1.25 nm；精度均為中等精度。溫度分別為77.00、87.00、293.15、373.15、473.15、573.15 K，壓力為0～1 000.00 kPa。模擬計(jì)算總步數(shù)為107步，前5×106步用于平衡體系，后5×106步用于統(tǒng)計(jì)平均值。圖1給出了在293.15 K和1 000.00 kPa的條件下3種碳基材料的平衡吸附示意圖。圖1中點(diǎn)狀物為氫氣分子，在模擬中通過(guò)調(diào)控體系內(nèi)的氫氣分子數(shù)量達(dá)到控壓效果。

圖1　293.15 K、1 000.00 kPa的條件下3種碳基材料的氫氣吸附平衡示意圖

2 結(jié)果與討論

2.1 吸附儲(chǔ)氫實(shí)驗(yàn)結(jié)果

圖2為石墨烯納米片、單壁碳納米管和多壁碳納米管在293.15 K、101.30 kPa下的氫氣吸附等溫線。由圖2可見(jiàn)，3種碳基材料的儲(chǔ)氫密度隨著氫氣壓力線性增加，其吸附等溫線不屬于任何一種常見(jiàn)的吸附等溫線類(lèi)型，但與C.Liu[28]和L.P.Ma[29]等報(bào)道的相同；在293.15 K和101.30 kPa的條件下，石墨烯、單壁碳納米管和多壁碳納米管的儲(chǔ)氫密度分別為0.05%、0.10%和0.06%，與G.G.Tibbetts[30]和G.Srinivas[31]等報(bào)道的基本一致。單壁碳納米管的儲(chǔ)氫能力最強(qiáng)，多壁碳納米管次之，石墨烯的儲(chǔ)氫能力最弱，這與它們的比表面積正相關(guān)，但比表面積與儲(chǔ)氫能力并無(wú)線性關(guān)聯(lián)，如多壁碳納米管的比表面積較石墨烯大很多，儲(chǔ)氫密度差值卻較小。

圖2　293.15 K、101.30 kPa的條件下3種碳基材料的氫氣吸附等溫線

2.2 模擬力場(chǎng)篩選

圖3為293.15 K、101.30 kPa的條件下石墨烯、單壁碳納米管和多壁碳納米管儲(chǔ)氫吸附的實(shí)驗(yàn)值和模擬值對(duì)比曲線。由圖3（a）可知，在293.15 K、101.30 kPa的條件下，石墨烯對(duì)氫氣的等溫吸附實(shí)驗(yàn)值小于不同力場(chǎng)計(jì)算的模擬值，且等溫吸附線同實(shí)驗(yàn)結(jié)果一樣近似直線；Universal力場(chǎng)的計(jì)算結(jié)果相對(duì)偏大，COMPASS和COMPASSⅡ力場(chǎng)計(jì)算的吸附等溫線幾乎重合，兩者差別很?。籇reiding力場(chǎng)相對(duì)最接近實(shí)驗(yàn)值。因此，可認(rèn)為Dreiding力場(chǎng)為計(jì)算石墨烯儲(chǔ)氫吸附的最優(yōu)力場(chǎng)。

由圖3（b）和（c）可知，對(duì)單壁碳納米管體系，COMPASS、COMPASSⅡ和Dreiding力場(chǎng)計(jì)算的儲(chǔ)氫密度過(guò)小，與實(shí)驗(yàn)值嚴(yán)重不符，而Universal力場(chǎng)的計(jì)算結(jié)果與實(shí)驗(yàn)值的偏差最??；對(duì)多壁碳納米管體系，Universal力場(chǎng)的計(jì)算結(jié)果與實(shí)驗(yàn)值最為接近。因此，Universal力場(chǎng)可作為計(jì)算單壁和多壁碳納米管的最優(yōu)力場(chǎng)。

圖3　293.15 K、101.30 kPa的條件下3種碳基材料儲(chǔ)氫吸附的實(shí)驗(yàn)值和模擬值對(duì)比曲線

2.3 吸附儲(chǔ)氫性能

圖4為不同壓力和溫度下石墨烯、單壁碳納米管和多壁碳納米管的氫氣吸附等溫線。由圖4可知，石墨烯在較低壓力下吸附儲(chǔ)氫能力最低，隨著壓力增加，無(wú)論在77.00 K或293.15 K條件下，其儲(chǔ)氫密度均線性增加，且當(dāng)壓力達(dá)到某一數(shù)值后，其吸附儲(chǔ)氫能力均高于碳納米管材料，在1 000.00 kPa的壓力下，溫度為77.00、293.15 K時(shí)對(duì)應(yīng)儲(chǔ)氫密度分別為1.79%和0.48%，遠(yuǎn)高于單壁碳納米管和多壁碳納米管；對(duì)碳納米管，隨著壓力的變化，不同溫度條件下吸附等溫線變化規(guī)律稍有不同，溫度較高時(shí)，儲(chǔ)氫密度在低壓下增加較快，隨壓力增加其增勢(shì)放緩，而在研究的壓力范圍和低溫條件下儲(chǔ)氫密度變化緩慢。

圖4　77.00、293.15 K的條件下氫氣的吸附等溫線

2.4 溫度對(duì)儲(chǔ)氫性能的影響

圖5為不同溫度下氫氣在石墨烯、單壁碳納米管、多壁碳納米管中的吸附等溫線。由圖5（a）可以看出，不同溫度下石墨烯的氫氣吸附等溫線都隨壓力呈線性增長(zhǎng)；隨著溫度上升，等溫線斜率逐漸減小，且減小趨勢(shì)逐漸減緩。由圖5（a）還可以看出，77.00、87.00 K的低溫有利于石墨烯儲(chǔ)氫，在這兩個(gè)溫度條件下增大壓力時(shí)石墨烯儲(chǔ)氫能力的提升幅度較大；在1 000.00 kPa下，溫度從77.00 K升高至87.00 K，石墨烯最終的儲(chǔ)氫密度從1.79%降低至1.59%；溫度升高至293.15 K時(shí)，儲(chǔ)氫密度顯著減小，但隨著溫度的進(jìn)一步升高，儲(chǔ)氫密度的下降并不明顯，在1 000.00 kPa的壓力下，溫度為293.15、373.15、473.15、573.15 K時(shí)儲(chǔ)氫密度分別為0.48%、0.38%、0.30%和0.25%。綜上，溫度越高，升高溫度對(duì)儲(chǔ)氫能力的負(fù)面影響越小。

圖5　不同溫度下氫氣在3種碳基材料中的吸附等溫線

由圖5（b）可以看出，在77.00、87.00 K的溫度下，單壁碳納米管在壓力為10.00 kPa時(shí)已吸附相當(dāng)數(shù)量的氫氣，隨著體系加壓，氫氣吸附密度增長(zhǎng)緩慢。P.Benard[32]等認(rèn)為，在低壓下物理吸附主要由表面結(jié)合能決定，在高壓下主要由比表面積決定，因此造成這種現(xiàn)象的原因可能是碳納米管比表面積有限，在這兩個(gè)超低溫條件下，即使壓力不大，氫氣也已經(jīng)接近飽和吸附狀態(tài)，并且在低溫下氫氣和碳納米管的結(jié)合作用可能較弱。

由圖5（c）可以看出，在77.00 K時(shí)的儲(chǔ)氫密度比87.00 K時(shí)大，約高0.17%；溫度從293.15 K升高至573.15 K，儲(chǔ)氫密度下降幅度則變小，其總體變化規(guī)律與石墨烯一致。根據(jù)圖5（a）和圖5（b）觀察到的吸附等溫線變化情況，推測(cè)單壁碳納米管吸附等溫線符合Ι類(lèi)等溫線[33?34]。為此，借助Langmuir模型和Freundlich模型對(duì)吸附等溫線進(jìn)行了擬合，模型具體公式為：

表1列出了不同溫度下氫氣在單壁碳納米管中的Langmuir和Freundlich吸附擬合參數(shù)。由表1可知，低溫下Langmuir模型擬合結(jié)果較差，77.00、87.00 K時(shí)回歸系數(shù)僅為0.544 10、0.779 28，而Freundlich模型擬合效果較好；溫度為293.15～573.15 K時(shí)，除573.15 K時(shí)Freundlich模型擬合不收斂，兩個(gè)模型的關(guān)聯(lián)結(jié)果均令人滿意，但相較于Freundlich模型，Langmuir模型擬合結(jié)果更好。

表1　不同溫度下氫氣在單壁碳納米管中的Langmuir和Freundlich吸附擬合參數(shù)

表2　不同溫度下氫氣在多壁碳納米管中的Langmuir和Freundlich吸附擬合參數(shù)

2.5 吸附熱變化

為定量評(píng)估石墨烯、單壁碳納米管和多壁碳納米管與氫氣分子的相互作用，分別計(jì)算了77.00 K、1 000.00 kPa和293.15 K、1 000.00 kPa條件下氫氣的等量吸附熱，結(jié)果見(jiàn)圖6。由圖6（a）可知，在77.00 K、1 000.00 kPa的條件下，3種碳基材料對(duì)氫氣的吸附熱均隨儲(chǔ)氫密度增大而減小，且吸附熱遵循單壁碳納米管(16～18 kJ/mol) >多壁碳納米管(9～11 kJ/mol)>石墨烯(1～3 kJ/mol)的規(guī)律。由此推斷，此條件下的吸附均為物理吸附[29]，且碳材料與氫氣分子結(jié)合力較弱；此外，氫氣在石墨烯上的吸附值為零時(shí)，石墨烯對(duì)氫氣的吸附熱約為3 kJ/mol，這與氫氣在碳材料上物理吸附的最小值相對(duì)應(yīng)[36]。由圖6（a）還可以看出，隨著儲(chǔ)氫密度增加，單壁碳納米管的吸附熱下降幅度最大，表明樣品表面存在多個(gè)結(jié)合能不同的吸附點(diǎn)。上述結(jié)果表明，為了提高石墨烯和碳納米管的儲(chǔ)氫能力，需要對(duì)其結(jié)構(gòu)和成分進(jìn)行修飾，從而在吸附過(guò)程中保持高比表面積。

由圖6（b）可知，293.15 K、1 000.00 kPa下石墨烯對(duì)氫氣的吸附熱明顯大于77.15 K、1 000.00 kPa下的吸附熱，而碳納米管體系在兩種條件下對(duì)氫氣的吸附熱相差不大，這與氫分子的熱運(yùn)動(dòng)隨溫度的升高而加劇有關(guān)。吸附熱是從游離態(tài)到吸附態(tài)所釋放的熱量，在吸附過(guò)程中，隨著氫分子向吸附材料表面靠近，分子熱運(yùn)動(dòng)速度逐步降低，同時(shí)勢(shì)能轉(zhuǎn)化為內(nèi)能而放熱，從而釋放出熱量；氫分子在293.15 K下的熱運(yùn)動(dòng)比77.00 K下要?jiǎng)×业枚?，?dāng)這些氫分子吸附在碳材料表面時(shí)，其動(dòng)能被轉(zhuǎn)化為熱量釋放出來(lái)，故吸附熱隨之增大。同時(shí)，在293.15 K下，材料表面氫氣覆蓋率較低，吸附熱的變化幅度比較小，但仍呈現(xiàn)隨儲(chǔ)氫密度增大而減小的規(guī)律。

圖6　壓力為1 000.00 kPa時(shí)、溫度為77.00 K和293.15 K的條件下氫氣的等量吸附熱

3 結(jié) 論

（1）借助GCMC分子模擬方法，計(jì)算了石墨烯、單壁碳納米管和多壁碳納米管對(duì)氫氣的吸附性能，并基于等溫吸附實(shí)驗(yàn)結(jié)果獲得了模擬不同體系的最佳力場(chǎng)：Dreiding力場(chǎng)計(jì)算石墨烯體系偏差最小，而Universal力場(chǎng)計(jì)算碳納米管系統(tǒng)偏差最小。

（2）在77.00 K、1 000.00 kPa的條件下，吸附儲(chǔ)氫能力排序?yàn)槭?單壁碳納米管>多壁碳納米管，對(duì)應(yīng)的最大儲(chǔ)氫密度分別為1.79%、1.02%和0.85%；在293.15 K、1 000.00 kPa的條件下，其排序?yàn)槭?多壁碳納米管>單壁碳納米管，對(duì)應(yīng)的最大儲(chǔ)氫密度分別為0.48%、0.28%和0.21%。

（3）氫氣在石墨烯中的吸附等溫線呈線性，符合一次函數(shù)關(guān)系；氫氣在碳納米管中的吸附符合Ι類(lèi)等溫線類(lèi)型，符合Langmuir和Freundlich模型。

（4）氫氣的吸附熱隨儲(chǔ)氫密度增大而減小，吸附熱排序?yàn)閱伪谔技{米管>多壁碳納米管>石墨烯。氫氣與碳材料的弱結(jié)合及氫氣覆蓋率的升高，導(dǎo)致可用比表面積減少，這極大地限制了材料的物理吸附儲(chǔ)氫性能。因此，需要對(duì)其進(jìn)行結(jié)構(gòu)和成分修飾，以保持固有的高比表面積，并增強(qiáng)與氫氣的結(jié)合，從而提高材料的吸附儲(chǔ)氫能力。

［1］ Zhang R F，Zhang Y Y，Wei F.Horizontally aligned carbon nanotube arrays： Growth mechanism，controlled synthesis，characterization，properties and applications［J］.Chemical Society Reviews，2017，46（12）：3661?3715.

［2］ Lee C G，Wei X D，Kysar J W，et al.Measurement of the elastic properties and intrinsic strength of monolayer graphene［J］.Science，2008，321（5887）：385?388.

［3］ Geim A K，Novoselov K S.The rise of graphene［J］.Nature Materials，2007，6（3）：183?191.

［4］ Lyu J Z，Kudiiarov V，Lider A.An overview of the recent progress in modifications of carbon nanotubes for hydrogen adsorption［J］.Nanomaterials，2020，10（2）：255?281.

［5］ Salehabadi A，Umar M F，Ahmad A，et al.Carbon‐based nanocomposites in solid‐state hydrogen storage technology：An overview［J］.International Journal of Energy Research，2020，44（14）：11044?11058.

［6］ 張紹鵬，建偉偉，馬丹竹，等.碳基及碳基復(fù)合材料吸附劑對(duì)VOCs吸附性能研究進(jìn)展［J］.遼寧石油化工大學(xué)學(xué)報(bào)，2021，41（1）：30?36.

Zhang S P，Jian W W，Ma D Z，et al.Research progress of adsorption properties of carbon adsorbents for VOCs［J］.Journal of Liaoning Petrochemical University，2021，41（1）：30?36.

［7］ Dillon A C，Jones K M，Bekkedahl T A，et al.Storage of hydrogen in single?walled carbon nanotubes［J］.Nature，1997，386（6623）：377?379.

［8］ Dogan E E，Tokcan P，Kizilduman B K.Storage of hydrogen in activated carbons and carbon nanotubes［J］.Advances in Materials Science，2018，18（4）：5?16.

［9］ Li Y，Liu H P.Grand canonical monte carlo simulation on the hydrogen storage behaviors of the cup?stacked carbon nanotubes at room temperature［J］.International Journal of Hydrogen Energy，2021，46（9）：6623?6631.

［10］ Petrushenko I K.DFT calculations of hydrogen adsorption inside single?walled carbon nanotubes［J］.Advances in Materials Science and Engineering，2018（1）：1?6.

［11］ Hassan I A，Ramadan H S，Saleh M A，et al.Hydrogen storage technologies for stationary and mobile applications：Review，analysis and perspectives［J］.Renewable and Sustainable Energy Reviews，2021，149：111311.

［12］ Zivar D，Kumar S，F(xiàn)oroozesh J.Underground hydrogen storage：A comprehensive review［J］.International Journal of Hydrogen Energy，2021，46（45）：23436?23462.

［13］ Jain V，Kandasubramanian B.Functionalized graphene materials for hydrogen storage［J］.Journal of Materials Science，2020，55（5）：1865?1903.

［14］ Ghosh A，Subrahmanyam K S，Krishna K S，et al.Uptake of H2and CO2by graphene［J］.Journal of Physical Chemistry C，2008，112（40）：15704?15707.

［15］ Alvarez?Zapatero P，Herrero A，Lebon A，et al.Ab initio study of lithium decoration of popgraphene and hydrogen storage capacity of the hybrid nanostructure［J］.International Journal of Hydrogen Energy，2021，46（29）：15724?15737.

［16］ Anikina E，Banerjee A，Beskachko V，et al.Li?functionalized carbon nanotubes for hydrogen storage：Importance of size effects［J］.Acs Applied Nano Materials，2019，2（5）：3021?3030.

［17］ Huo Y J，Zhang Y，Wang C B，et al.Boron?doping effect on the enhanced hydrogen storage of titanium?decorated porous graphene：A first?principles study［J］.International Journal of Hydrogen Energy，2021，46（80）：40301?40311.

［18］ Rather S.Hydrogen uptake of Ti?decorated multiwalled carbon nanotube composites［J］.International Journal of Hydrogen Energy，2021，46（34）：17793?17801.

［19］ Singh S B，De M.Improved hydrogen uptake of metal modified reduced and exfoliated graphene oxide［J］.Journal of Materials Research，2021，36（15）：3109?3120.

［20］ Lu X，Chu W，Sun W J，et al.Enhancement of hydrogen sorption on metal（Ni，Rh，Pd） functionalized carbon nanotubes：A DFT study［J］.Chemical Research in Chinese Universities，2017，33（3）：422?429.

［21］ Liu P P，Zhang Y F，Xu X J，et al.Ti decorated B8 as a potential hydrogen storage material：A DFT study with van der Waals corrections［J］.Chemical Physics Letters，2021，765：138277.

［22］ Srivastava S，Shet S P，Priya S S，et al.Molecular simulation of copper based metal?organic framework （Cu?MOF） for hydrogen adsorption［J］.International Journal of Hydrogen Energy，2022，47（35）：15820?15831.

［23］ Kag D，Luhadiya N，Patil N D，et al.Strain and defect engineering of graphene for hydrogen storage via atomistic modelling［J］.International Journal of Hydrogen Energy，2021，46（43）：22599?22610.

［24］ Froudakis G E .Hydrogen interaction with carbon nanotubes：A review of ab initio studies［J］.Journal of Physics Condensed Matter，2002，14（17）：453?465.

［25］ Lamari D F，Malbrunot P，Tartaglia G P.Review of hydrogen storage by adsorption in carbon nanotubes［J］.International Journal of Hydrogen Energy，2002，27（2）：193?202.

［26］ Dicko M，Seydou M，Lamari F D，et al.Hydrogen adsorption on graphane：An estimate using ab?initio interaction［J］.International Journal of Hydrogen Energy，2017，42（15）：10057?10063.

［27］ Rayón V M，Cabria I.Assessment of density functional approximations for N2and CO2physisorption on benzene and graphene［J］.Journal of Computational Chemistry，2022，43（21）：1403?1419.

［28］ Liu C，F(xiàn)an Y Y，Liu M，et al.Hydrogen storage in single?walled carbon nanotubes at room temperature［J］.Science，1999，286（5442）：1127?1129.

［29］ Ma L P，Wu Z S，Li J，et al.Hydrogen adsorption behavior of graphene above critical temperature［J］.International Journal of Hydrogen Energy，2009，34（5）：2329?2332.

［30］ Tibbetts G G，Meisner G P，Olk C H.Hydrogen storage capacity of carbon nanotubes，filaments，and vapor?grown fibers［J］.Carbon，2001，39（15）：2291?2301.

［31］ Srinivas G，Zhu Y W，Piner R，et al.Synthesis of graphene?like nanosheets and their hydrogen adsorption capacity［J］.Carbon，2010，48（3）：630?635.

［32］ Benard P，Chahine R.Storage of hydrogen by physisorption on carbon and nanostructured materials［J］.Scripta Materialia，2007，56（10）：803?808.

［33］ 鄧佳，張奇，王棟，等.壓力驅(qū)動(dòng)條件下頁(yè)巖微納米孔隙CO2/CH4競(jìng)爭(zhēng)吸附特性［J］.東北石油大學(xué)學(xué)報(bào)，2021，45（5）：109?116.

Deng J，Zhang Q，Wang D，et al. Competitive adsorption characteristics of CO2/CH4in micro?nano pores by pressure?driven for shale reservoir［J］.Journal of Northeast Petroleum University，2021，45（5）：109?116.

［34］ 高良軍，葉佳君，陳彩勝，等.貝殼粉海洋油污吸附材料的改性實(shí)驗(yàn)［J］.油氣儲(chǔ)運(yùn)，2018，37（3）：343?347.

Gao L J，Ye J J，Chen C S，et al.Modification experiment of shell powder marine oil adsorption material［J］.Oil & Gas Storage and Transportation，2018，37（3）：343?347.

［35］ 朱孟爽，羅根祥，許凌子，等.Mg?Al?LDH對(duì)塔河瀝青質(zhì)及其亞組分的吸附研究［J］.石油化工高等學(xué)校學(xué)報(bào)，2022，35（3）：10?15.

Zhu M S，Luo G X，Xu L Z，et al.Study on the adsorption of asphaltenes and their sub?components by Mg?Al?LDH in Tahe［J］.Journal of Petrochemical Universities，2022，35（3）：10?15.

［36］ 楊宵，劉晶，胡建波.氫氣在共價(jià)有機(jī)骨架材料中的吸附機(jī)理［J］.化工學(xué)報(bào)，2015，66（7）：2540?2546.

Yang X，Liu J，Hu J B.Adsorption mechanism of H2on covalent organic frameworks［J］.CIESC Journal，2015，66（7）：2540?2546.

Molecular Simulation of Hydrogen Adsorption Performance in Graphene and Carbon Nanotubes

Tang Mengyao， Sun Tianyi， Li Jiashu， Ding Yan， Li Jinlong

（School of Petrochemical Engineering，Changzhou University，Changzhou Jiangsu 213164，China）

The adsorption data of hydrogen in graphene,single?walled carbon nanotubes,and multi?walled carbon nanotubes at 293.15 K and 101.30 kPa were measured experimentally.Through the comparison of the data calculated by different force fields under the same conditions, the optimal force fields of the three carbon materials for calculation were selected.On this basis,the hydrogen adsorption data of the three carbon materials at 0～1 000.00 kPa and 77.00～573.15 K were calculated.The results show that the Dreiding force field is the optimal force field to calculate the hydrogen adsorption and storage in graphene,and the Universal force field is the best force field to calculate that in carbon nanotubes.Under given conditions,the hydrogen adsorption and storage capacity of the three materials is ranked as follows:graphene single?walled carbon nanotubes multi?walled carbon nanotubes.The hydrogen storage capacity is closely related to the specific surface area of the material and its weak binding force with hydrogen.The results can provide data and theoretical support for the molecular simulation of hydrogen adsorption and storage in carbon materials as well as hydrogen storage material design.

Graphene； Carbon nanotubes； Adsorption hydrogen storage performance； Molecular simulation

1006?396X（2022）04?0010?08

2022?06?13

2022?08?02

國(guó)家自然科學(xué)基金項(xiàng)目（22078026, 21878025）。

湯夢(mèng)瑤（1998?），女，碩士研究生，從事分子熱力學(xué)與藥物活性分離方面的研究；E?mail:870019653@qq.com。

李進(jìn)龍（1976?），男，博士，教授，從事分子熱力學(xué)實(shí)驗(yàn)與理論方面的研究；E?mail:lijinlong@cczu.edu.cn。

TQ021.2

A

10.3969/j.issn.1006?396X.2022.04.002

（編輯 閆玉玲）