Modulation of Molecular Sensing Properties of Graphdiyne Based on 3d Impurities

CHEN Xi , ZHANG Shengli

1 Department of Applied Physics, School of Physics and Optoelectronic Engineering, Xidian University, Xi’an 710071, P. R. China.

2 Department of Applied Physics, School of Science, Xi’An Jiaotong University, Xi’an 710049, P. R. China.

3 MOE Key Laboratory for Non-equilibrium Synthesis and Modulation of Condensed Matter, Xi’an Jiaotong University, Xi’an 710049, P. R. China.

Abstract: In recent years, the successful preparation of single-layer graphene, MoS2, and other two-dimensional materials has started a new era of two-dimensional materials. The potential applications of twodimensional materials in emerging electronics have drawn widespread attention. Two-dimensional carbon materials, with their unique properties, have become the research hotspot of condensed matter physics, nanoelectronics, and biological medicine. The remarkable success in preparing graphene provides additional possibilities for developing sensitive biodevices and medicine systems. However, graphene is gapless and thus is unsuitable for building nanoelectronic devices or biosensors due to the too low on/off current ratio. More than 20 years ago, graphyne and its family (viz. graphdiyne, graphyne-3, etc.), as hypothetical C allotropes, were theoretically predicted to be semiconductors with a layered structure. Recently, graphdiyne was successfully synthesized on the surface of copper via a cross-coupling reaction using hexaethynylbenzene. Graphdiyne, as a new two-dimensional carbon material with semiconductor properties and a unique porous structure, is more advantageous than graphene for nanoelectronic and biosensing applications. As the first discovered semiconducting two-dimensional carbon material, with independent intellectual property rights in China,graphdiyne has great research significance. Compared with graphene, graphdiyne has a unique structure with larger pores composed of high π-conjugated acetylenic bonds, which may facilitate strong adsorption to biomolecules. Therefore,further research is needed to reveal how the physical properties of graphdiyne can be modulated effectively to meet the requirements of practical applications. The interaction between biological molecules and materials is an important subject of research in condensed matter physics and materials science. Detailed understanding of the interactions between graphdiyne and small molecules may facilitate the development of advanced biological applications such as biosensors for the detection of biomolecules and living cells, drug delivery systems, and cell imaging technologies. In sensitive analysis,the ultimate goal is to achieve reliable detection of trace amounts of molecules. In this work, first-principles calculations were employed to investigate the electronic structure of graphdiyne nanoribbons and the adsorption of graphdiyne to small molecules. To improve the chemical response of graphdiyne to single molecules, we considered modifying graphdiyne by doping 3d transition metal atoms. We chose Sc and Ti, which have the largest adsorption energies on graphdiyne, and studied the room-temperature stabilities of Sc- and Ti-doped graphdiyne and the possibility of using Sc- and Ti-doped graphdiyne as materials for molecular sensing. Finally, we investigated the interaction between graphdiyne and amino acid molecules and discovered that the dispersion force plays a large role in the interaction. The influence of amino acids on the electronic transport properties of graphdiyne was also studied, and the potential applications of graphdiyne to biosensors were investigated.

Key Words: Graphdiyne; Two-dimensional material; Electronic transport; Nanodevice

1 引言

碳元素是人類最早認(rèn)識(shí)的元素之一，也是自然界成鍵結(jié)構(gòu)最豐富的元素1。碳元素可構(gòu)成許多具有特殊性質(zhì)的材料2–6。在古代人們便發(fā)現(xiàn)了碳的同素異形體，如石墨和金剛石。隨著近二十年來實(shí)驗(yàn)技術(shù)的發(fā)展，人們成功制備了許多不同結(jié)構(gòu)的碳納米結(jié)構(gòu)，對(duì)材料技術(shù)的進(jìn)步具有里程碑意義。含碳材料由于其具有的獨(dú)特電子特性，在某些方面有許多先進(jìn)應(yīng)用，如太陽能電池，有機(jī)發(fā)光二極管，場(chǎng)效應(yīng)管，化學(xué)傳感器等。

2004年，Geim等人以單晶石墨為原料，得到了單原子層的石墨烯7,8。此后，人們發(fā)現(xiàn)了越來越多的二維材料。二維材料電子特性的潛在應(yīng)用引起了廣泛的關(guān)注。由于二維單原子層具有特殊的結(jié)構(gòu)和獨(dú)特的物理性質(zhì)8–11，因此一直是眾多領(lǐng)域中科學(xué)工作者研究的熱點(diǎn)。然而，零帶隙的石墨烯不適用于制備晶體管等電子器件。早在數(shù)十年前，就有人從理論上預(yù)測(cè)出一類由sp和sp2雜化的碳原子構(gòu)成碳的同素異形體—石墨炔12,13。石墨炔具有平面層狀結(jié)構(gòu)，是由碳六元環(huán)和炔鍵構(gòu)成的半導(dǎo)體。近年來，中國科學(xué)院化學(xué)研究所的李玉良院士課題組在實(shí)驗(yàn)上首次在銅表面成功制備出γ-石墨二炔(屬于石墨炔中的一種，以下簡稱石墨二炔)14,15。石墨二炔是第一種具有半導(dǎo)體性質(zhì)的碳的同素異形體，帶隙與硅接近，比石墨烯更適用于制備納米電子器件15–17，是新一代場(chǎng)效應(yīng)管材料的候選者。

石墨二炔是我國獨(dú)立研發(fā)并具有自主知識(shí)產(chǎn)權(quán)的新型二維碳材料。在技術(shù)發(fā)展迅猛的今天，傳統(tǒng)電子學(xué)器件不斷地小型化，在10 nm尺寸以下即將面受量子效應(yīng)瓶頸。因此，人們迫切需要找到能取代傳統(tǒng)硅基材料的新型納米半導(dǎo)體材料。在納米科技發(fā)展日新月異的今天，我們亟需發(fā)展基于石墨二炔的二維電子學(xué)及生物器件。其中所面臨的主要問題有二：如何有效調(diào)控石墨二炔物理性質(zhì)使其滿足實(shí)際應(yīng)用的需求；如何根據(jù)石墨二炔的特性提出電子學(xué)器件的設(shè)計(jì)方案。石墨二炔具有獨(dú)特的由 π共軛炔鍵構(gòu)成的大孔網(wǎng)狀結(jié)構(gòu)，可被應(yīng)用于氣體分離膜、儲(chǔ)能材料、電池電極材料、太陽能電池材料18–20等領(lǐng)域。此外，石墨二炔亦對(duì)生物分子有吸附作用?；谑驳莫?dú)特優(yōu)勢(shì)，我們考慮石墨二炔將應(yīng)用到生物器件領(lǐng)域。由于氨基酸是生命體基本單元，我們需掌握氨基酸及其它生物小分子與石墨二炔表面相互作用規(guī)律及其物理機(jī)制，這是設(shè)計(jì)基于石墨二炔的納米生物傳感器的基礎(chǔ)。

分子檢測(cè)的理想目標(biāo)為達(dá)到終極的檢測(cè)靈敏度，即能夠?qū)蝹€(gè)分子進(jìn)行解析。迄今為止，人們已經(jīng)發(fā)展了大量的檢測(cè)技術(shù)手段，如掃描隧道顯微鏡、原子力顯微鏡21和熒光關(guān)聯(lián)光譜22等等?，F(xiàn)代的掃描隧道顯微鏡技術(shù)，已能夠在單分子水平上操作化學(xué)合成23。熒光關(guān)聯(lián)光譜主要用于檢測(cè)聚合物折疊24、大分子構(gòu)象的改變25及較快的酶動(dòng)力學(xué)過程26,27?，F(xiàn)今，人們更多地把注意力集中在小尺寸的分子檢測(cè)器上，并要求具備高靈敏度及易用性。二維材料的誕生給氣體檢測(cè)器的制備開拓了更多的可能性。二維材料具有非常大的比表面積，高的靈敏度。氣體分子在二維材料表面吸附后，起到了電子施主或受主的作用，能引起二維材料載流子濃度的改變，從而形成電信號(hào)響應(yīng)。人們已經(jīng)利用石墨烯的霍爾效應(yīng)來檢測(cè)吸附氣體分子造成的載流子濃度改變28，并解析出單個(gè)分子的信號(hào)?；谑┑膱?chǎng)效應(yīng)管也開發(fā)并用于檢測(cè)單分子的電信號(hào)29–32。為了提高檢測(cè)的靈敏度，人們考慮各種對(duì)二維材料表面進(jìn)行修飾的方法。過渡金屬原子往往可作為小分子的吸附中心。而石墨烯表面的π鍵是飽和的，具有化學(xué)惰性；因此金屬原子在石墨烯表面的吸附能低于2 eV，遷移勢(shì)壘低于0.8 eV33,34。較低的吸附能和遷移勢(shì)壘將使得金屬原子在室溫下聚集。相比之下，石墨二炔比石墨烯更容易吸附金屬原子，使之成為潛在的分子吸附材料。這是石墨二炔在分子檢測(cè)材料中的優(yōu)勢(shì)。

在本文中，我們首先研究了石墨二炔納米帶的電子輸運(yùn)性質(zhì)。之后，考慮利用石墨二炔的特殊大孔碳網(wǎng)狀結(jié)構(gòu)作為小分子傳感器材料。為了提高石墨二炔的響應(yīng)，在3d過渡金屬中，我們挑選出了在石墨二炔表面吸附能最大的Sc、Ti原子，并確定了Sc、Ti摻雜石墨二炔在室溫下的穩(wěn)定性。我們選擇HCHO作為典型的有代表性的小分子，從能帶、載流子濃度等方面探討了Sc、Ti摻雜石墨二炔對(duì)HCHO分子的響應(yīng)。我們還研究了氨基酸與石墨二炔的相互作用及物理機(jī)制，得出氨基酸對(duì)石墨二炔光吸收和電輸運(yùn)性質(zhì)影響的規(guī)律。本研究為今后將石墨二炔用于新型生物分子電子學(xué)及生物學(xué)檢測(cè)器件，在生物電子學(xué)及藥物載體方面的應(yīng)用提供重要理論指導(dǎo)。

2 理論模型與計(jì)算方法

在本文中，我們主要使用SIESTA軟件包35，采用密度泛函理論計(jì)算石墨二炔的電子性質(zhì)。在計(jì)算中采用Perdew-Burke-Ernzerhof (PBE)36的廣義梯度近似(GGA)方法來描寫電子的交換關(guān)聯(lián)作用，引入了包含色散力的PBE-D2修正37，以及模守恒的Troullier-Martins贗勢(shì)38。實(shí)空間積分格點(diǎn)截?cái)嗳?50 Ry。使用了極化的雙ζ基組對(duì)波函數(shù)進(jìn)行展開。對(duì)二維的石墨二炔單層，沿與其平面垂直的方向加上了2.0 nm的真空層。進(jìn)行計(jì)算前，我們先將各種結(jié)構(gòu)進(jìn)行馳豫，直到所有原子受力皆小于0.1 eV·nm?1。在分子動(dòng)力學(xué)模擬中，我們使用了Verlet算法對(duì)原子的牛頓方程進(jìn)行積分，時(shí)間步長為1 fs，采用Nose熱浴控制體系溫度。

考慮到PBE泛函往往低估帶隙，因此我們?cè)谀軒в?jì)算中采用雜化的 Heyd-Scuseria-Ernzerhof(HSE06)泛函39,40，配合 VASP 軟件41–43在 PAW 方法下進(jìn)行計(jì)算。平面波能量截?cái)酁?00 eV。在2 ×2 × 1的Γ中心Monkhorst-Pack采樣下選取k點(diǎn)進(jìn)行自洽計(jì)算，而后在布里淵區(qū)高對(duì)稱路徑上取30個(gè)k點(diǎn)計(jì)算能帶。

在介觀輸運(yùn)計(jì)算中，為了縮短計(jì)算時(shí)間，我們使用了極化的單ζ基組。納米帶k點(diǎn)用1 × 1 × 100進(jìn)行采樣。在系統(tǒng)兩端施加偏壓后，用非平衡格林函數(shù)方法44(TRANSIESTA軟件包45)計(jì)算電子的透射率，然后根據(jù) Landauer-Buttiker公式計(jì)算電流。

3 石墨二炔納米帶電子結(jié)構(gòu)與電子輸運(yùn)性質(zhì)

圖1 (a) 2 × 2的石墨二炔原子結(jié)構(gòu)，其中虛線范圍表示原胞；(b) PBE和HSE06泛函給出的2 × 2石墨二炔能帶(左)和態(tài)密度(右)Fig. 1 (a) The structure of graphdiyne 2 × 2 supercell, primitive cell is denoted by dashed lines;

石墨二炔是由sp雜化的碳原子鏈與sp2雜化的碳六圓環(huán)組成的二維平面碳網(wǎng)。從圖1a中我們看出，石墨二炔為二維平面網(wǎng)絡(luò)狀的碳結(jié)構(gòu)，其中碳碳三鍵(sp雜化)是構(gòu)成石墨二炔結(jié)構(gòu)中十分重要的連接單元。由于它不會(huì)受順反異構(gòu)的變化影響，因此可以一直保持線性的結(jié)構(gòu)。其次，它具有更小的空間位阻，這很有利于將sp雜化的碳連接到sp2雜化或者sp雜化的碳原子中心上。另外，碳碳三鍵與苯環(huán)之間形成了離域π鍵，因此其高度共軛、碳富集的有機(jī)分子結(jié)構(gòu)具有優(yōu)良的靈活可調(diào)性。雖然石墨烯具有比硅更快的電子移動(dòng)速率，但是石墨烯并不具有明顯的能帶隙，不適合被應(yīng)用于場(chǎng)應(yīng)效晶體管。而石墨二炔具有直接帶隙，其值約為1 eV。相比石墨烯，石墨二炔更適合被應(yīng)用于場(chǎng)效應(yīng)晶體管。值得注意的是，石墨二炔的帶隙值與硅的(1.1 eV)很相近。由于以硅為基礎(chǔ)的科學(xué)技術(shù)擁有其自身的局限性，石墨二炔將有可能代替硅來制成更快更小的電子學(xué)器件。而且石墨二炔直接帶隙的存在，促進(jìn)了光—電的高效轉(zhuǎn)換，有助于其在光電學(xué)器件上的應(yīng)用。這些優(yōu)點(diǎn)就使得石墨二炔很可能替代硅成為電子學(xué)器件的優(yōu)異材料。

石墨二炔的原胞如圖1a中的菱形虛線所示，菱形的邊長即為石墨二炔的晶格常數(shù)a0。我們的計(jì)算得到a0= 0.950 nm，與文獻(xiàn)46的計(jì)算值a0=0.948 nm吻合。能帶計(jì)算表明，石墨二炔具有直接帶隙。在后續(xù)的計(jì)算中(見第4、5節(jié))，為了研究系小分子吸附，我們?nèi)×? × 2晶胞作為襯底模型進(jìn)行能帶計(jì)算。圖1b左給出了PBE和HSE06泛函所計(jì)算得到的能帶及態(tài)密度。我們通過PBE泛函計(jì)算出的帶隙為0.44 eV，與文獻(xiàn)值0.46 eV46相近。HSE06泛函計(jì)算出的帶隙為0.83 eV，與文獻(xiàn)值(0.88 eV)47接近。這里需指出的是，普通的密度泛函方法計(jì)算的帶隙總是偏小，PBE計(jì)算值低于帶隙實(shí)驗(yàn)值。在文獻(xiàn)47中，準(zhǔn)粒子(GW)準(zhǔn)粒子近似給出的帶隙1.10 eV最接近實(shí)驗(yàn)值。GW準(zhǔn)粒子近似對(duì)帶隙的修正將近1.5倍，這是由低維材料中較強(qiáng)的電子間庫倫相互作用所引起。圖1b右展示了PBE和HSE06泛函所給出的石墨二炔態(tài)密度。

為了研究石墨二炔作為介觀器件材料的電子輸運(yùn)特性，我們考慮以平直的石墨二炔納米帶作為理論研究模型。在進(jìn)行計(jì)算之前，先對(duì)石墨二炔納米帶的結(jié)構(gòu)進(jìn)行分類。我們考慮了扶手椅型(AGDNR)和鋸齒型石墨二炔納米帶(ZGDNR)。各種結(jié)構(gòu)、各種寬度和各種邊緣形狀的石墨二炔納米帶結(jié)構(gòu)如圖2a所示。n-AGDNR邊緣平直，帶有六元環(huán)。這里的整數(shù)n表示其相對(duì)寬度(橫向碳原子排數(shù))。ZGDNR邊緣為突起狀，類似于鋸齒型石墨烯納米帶。但ZGDNR突起比鋸齒型石墨烯納米帶大得多。ZGDNR分為兩類，其中n-ZGDNR兩邊緣是對(duì)稱的，而n’-ZGDNR兩邊緣的鋸齒位置是錯(cuò)開的。納米帶沿z方向具有周期性，我們選擇基本重復(fù)單元(圖2a黑線框)作為原胞進(jìn)行幾何結(jié)構(gòu)優(yōu)化。計(jì)算結(jié)果表明，圖2a中各種納米帶的碳六圓環(huán)間距a與石墨二炔晶格常數(shù)a0差別在2%以內(nèi)，這說明石墨二炔的納米帶結(jié)構(gòu)與層狀石墨二炔接近。

圖2 (a)不同寬度的扶手椅型石墨二炔納米帶和兩種類型的鋸齒型石墨二炔納米帶，分別表示為n-AGDNR、n-ZGDNR及n’-ZGDNR，其中黑線框標(biāo)出每條納米帶的晶胞；(b)鐵磁態(tài)(FM)的3-ZGDNR與3’-ZGDNR的能帶圖，邊緣態(tài)能帶用藍(lán)色箭頭標(biāo)出；(c)鐵磁態(tài)(FM) 3-ZGDNR的電荷密度圖及其邊緣態(tài)軌道態(tài)密度(PDOS)圖Fig. 2 (a) Armchair and two different types of graphdiyne nanoribbons with different widths, which are called

從能帶結(jié)構(gòu)可以看出，各種不同寬度的石墨二炔納米帶均為直接帶隙半導(dǎo)體，帶隙出現(xiàn)在Γ點(diǎn)。此處為了節(jié)省計(jì)算時(shí)間，用PBE泛函計(jì)算帶隙，僅用于定性比較。2-，3-，4-和5-AGDNR的帶隙分別為0.79，0.70，0.65和0.61 eV。由于量子限域效應(yīng)，石墨二炔納米帶寬度越小，帶隙越大。隨著石墨二炔納米帶寬度的增大，帶隙下降并越來越接近無限大石墨二炔平面的PBE帶隙值(0.47 eV)。n-ZGDNR (表1)和n’-ZGDNR (表2)的帶隙亦表現(xiàn)出此特性。

我們考慮了石墨二炔納米帶的各種磁性狀態(tài)。計(jì)算結(jié)果表明，各種AGDNR均不具有磁性。而 ZGDNR具有磁性，其鐵磁(FM)態(tài)能量比相應(yīng)的非磁態(tài)(NM)能量低(表 1和表 2分別給出了 n-ZGDNR和n’-ZGDNR的帶隙)。FM態(tài)的磁矩主要分布在 ZGDNR邊緣突起部分的原子上。對(duì)各種不同形狀、不同寬度的n-ZGDNR，F(xiàn)M態(tài)比NM態(tài)的每晶胞能量低1.01–2.12 eV。對(duì)n’-ZGDNR，F(xiàn)M態(tài)比NM態(tài)的每晶胞能量低2.00–2.04 eV。此外，我們還考慮了兩種反鐵磁(AFM)態(tài)。第一種情況為，相鄰ZGDNR晶胞磁矩反平行排列(AFM1)。結(jié)果表明對(duì)同種ZGDNR，AFM1態(tài)的能量比FM態(tài)高約0.01–0.05 eV。第二種情況為，納米帶兩邊緣磁矩反平行排列的情形(AFM2)。結(jié)果表明對(duì)同種ZGDNR，AFM2態(tài)的能量比同種結(jié)構(gòu)的FM態(tài)高約0.08–0.10 eV。綜上所述，F(xiàn)M態(tài)為ZGDNR的最穩(wěn)定狀態(tài)。

表1 n-ZGDNR的帶隙(PBE泛函)、FM態(tài)與NM態(tài)能量差(EFM ? ENM)、FM態(tài)晶胞磁矩(μ)Table 1 Band gap, total energy difference between FM and NM state (EFM ? ENM) and the magnetic moment of FM unit cell (μ) of n-ZGDNRs at the level of PBE.

表2 n’-ZGDNR的帶隙(PBE)、FM態(tài)與NM態(tài)能量差(EFM ? ENM)、FM態(tài)晶胞磁矩(μ)Table 2 Band gap, total energy difference between FM and NM state (EFM ? ENM) and the magnetic moment of FM unit cell (μ) of n’-ZGDNRs.

對(duì)各種形狀、各種寬度的ZGDNR，由于其鐵磁基態(tài)特性，能帶均表現(xiàn)出自旋極化特征。作為例子，圖2b顯示了FM態(tài)的3-ZGDNR和3’-ZGDNR的能帶。在能帶圖中，能夠明顯看到若干特殊的呈現(xiàn)出幾乎平直E–k關(guān)系的能帶(藍(lán)色箭頭)。我們猜測(cè)這些電子態(tài)為局域電子軌道。為了揭示這些平直能帶的本質(zhì)，我們計(jì)算了它們的電荷密度分布，及其軌道投影態(tài)密度(PDOS) (圖2c)。從電荷分布可以看出，其電子云主要局域在鋸齒型邊緣的突起處，表明這些局域態(tài)為納米帶的邊緣態(tài)。由于ZGDNR的邊緣突起彼此間隔距離較遠(yuǎn)(比鋸齒型石墨烯邊緣突起間距大)，各個(gè)鋸齒尖端的邊緣態(tài)的相互作用較弱。所以這些局域態(tài)相對(duì)獨(dú)立，形成了這些平直帶。從其 PDOS可以看出，邊緣態(tài)主要由2s軌道及面內(nèi)2py+ 2pz軌道組成。

為了研究石墨二炔納米帶的介觀電子輸運(yùn)性質(zhì)，我們構(gòu)造了圖3a所示的模型。以納米帶為兩端電極(模擬時(shí)取原胞加上周期邊界條件)、3倍原胞長度作為散射區(qū)進(jìn)行計(jì)算。在兩端施加電壓 Vb后，計(jì)算通過體系的電流I。我們也嘗試過選取更長的散射區(qū)，結(jié)果 I–Vb曲線無明顯變化，因此 3倍原胞長度的散射區(qū)是足夠大的。對(duì)各種納米帶，I–Vb曲線呈現(xiàn)出半導(dǎo)體特征。圖 3b展示了 2–4-AGDNR的I–Vb曲線。當(dāng)eVb大于帶隙時(shí)電流I為幾到十幾μA量級(jí)，當(dāng)eVb小于帶隙時(shí)電流I極小(小于0.01 μA)。從圖3b可以看出，AGDNR的寬度越大，帶隙越小，導(dǎo)通電壓越小。當(dāng)eVb大于帶隙時(shí)，AGDNR導(dǎo)通，I–Vb曲線近似為直線。在導(dǎo)通區(qū)域，2–4-AGDNR的電導(dǎo)約為0.77–1.31 G0。4-AGDNR的電導(dǎo)最大。對(duì)各種 n-ZGDNR和 n’-ZGDNR，I–Vb曲線也呈現(xiàn)出半導(dǎo)體特征。圖3c呈現(xiàn)了 n-ZGDNR 的 I–Vb曲線。圖 3d呈現(xiàn)了 n’-ZGDNR的 I–Vb曲線。這些曲線均具備半導(dǎo)體特征，與 AGDNR類似?？偟膩砜?，納米帶寬度越寬，導(dǎo)通后的導(dǎo)電能力越強(qiáng)。對(duì) 3-ZGDNR、3’-ZGDNR、4-ZGDNR、4’-ZGDNR，導(dǎo)通段的電導(dǎo)約為0.26–0.60 G0。值得一提的是，F(xiàn)M的ZGDNR產(chǎn)生自旋極化電流。自旋極化率定義為P = | I↑?I↓|/(I↑+ I↓)，其中 I↑和 I↓分別為自旋向上和自旋向下電流。在導(dǎo)通段，3-ZGDNR、3’-ZGDNR、4-GDNR和 4’-ZGDNR的電流自旋極化率為 23%–64%。

4 Sc、Ti摻雜石墨二炔對(duì)HCHO分子的吸附與響應(yīng)

圖3 (a)石墨二炔納米帶介觀輸運(yùn)計(jì)算模型，(b) 2–4-AGDNR的I–Vb曲線，(c) 2–4-ZGDNR的I–Vb 曲線，(d) 2’–4’-ZGDNR 的 I–Vb 曲線Fig. 3 (a) The calculation model of electronic transport in graphdiyne nanoribbons, (b) I–Vb curves of 2–4-AGDNRs,(c) I–Vb curves of 2–4-ZGDNRs, (d) I–Vb curves of 2’–4’-ZGDNRs.

接下來，我們考慮將石墨二炔作為小分子傳感器的可能性。欲作為傳感器材料，首先要求此種材料對(duì)小分子具有一定的吸附能力。目前的研究表明，石墨炔的飽和炔鍵結(jié)構(gòu)對(duì)一些小分子(如H2、O2、CO、N2、CH4)的吸附較弱48,49，其相互作用主要是微弱的排斥作用。為了提高石墨二炔的吸附能力，我們考慮在其大孔中加入3d金屬原子。

這里選擇HCHO分子作為典型的研究對(duì)象。為了計(jì)算HCHO分子在石墨二炔上的吸附，我們選取各種不同初始位置及不同的HCHO分子取向來進(jìn)行結(jié)構(gòu)優(yōu)化，尋找最穩(wěn)定的吸附構(gòu)型(即吸附能最大的構(gòu)型)。模擬中采用2 × 2的石墨二炔晶胞，在2 × 2 × 1的Γ中心Monkhorst-Pack采樣下選取k點(diǎn)進(jìn)行自洽計(jì)算。結(jié)果表明，當(dāng)HCHO分子位于三角形炔鏈中心處(圖 4a)時(shí)為在石墨二炔表面最穩(wěn)定的吸附構(gòu)型，吸附能為 Ead= 0.43 eV(表3)。此外，若HCHO分子位于石墨二炔的碳六圓環(huán)上，則吸附能為負(fù)值，為不穩(wěn)定位置。相比之下，我們也獲得了HCHO在石墨烯表面的最穩(wěn)定吸附位置，相應(yīng)的最大吸附能是Ead= 0.28 eV，明顯小于在石墨二炔表面的吸附能。作為參考，我們以相同的方法也搜索了CO和C6H6分子在石墨二炔和石墨烯上的最穩(wěn)定吸附構(gòu)型。CO分子在石墨二炔和石墨烯上的最大吸附能Ead(CO)分別是0.55和0.28 eV，C6H6分子在石墨二炔和石墨烯上的最大吸附能Ead(C6H6)分別是0.81和0.75 eV (表3)。結(jié)果表明，兩種分子在石墨二炔上的吸附能都比在石墨烯上大。與石墨烯中的烯鍵相比，石墨二炔中炔鏈中額外的π電子可能會(huì)導(dǎo)致吸附增強(qiáng)。

圖4 (a) HCHO分子在石墨二炔上吸附的側(cè)視及俯視圖；(b) C、H、X表示Sc、Ti原子三種可能的吸附位置Fig. 4 (a) Side and top views of HCHO adsorption on graphdiyne; (b) three possible adsorption sites of Sc and Ti atoms shown by C, H and X.

表3 HCHO、CO、C6H6分子在石墨二炔和石墨烯上的吸附能EadTable 3 The adsorption energy Ead of HCHO,CO and C6H6 molecule on graphdiyne and graphene.

為了將3d金屬原子摻雜的石墨二炔作為分子傳感材料，我們首先要確定其在室溫下的穩(wěn)定性。對(duì)3d金屬原子，我們考慮了在石墨二炔表面三種可能的吸附位置，如圖4b所示。分別為由炔鏈組成的三角形大孔中心X位置、大孔角上C位置和六圓環(huán)中心H位置。幾何優(yōu)化結(jié)果表明，C位置和H位置是穩(wěn)定的吸附位置。而X位置不穩(wěn)定，因?yàn)樵趦?yōu)化過程中，Sc、Ti原子會(huì)自發(fā)遷移至 C位置。對(duì)于ScTi原子在C位置的吸附能Ead(C)大于H位置的吸附能Ead(H)，如表4所示。我們也考慮了其他的3d過渡金屬原子的吸附，發(fā)現(xiàn)Sc、Ti原子在石墨二炔表面的吸附能大于或接近 V、Cr、Mn，F(xiàn)e、Co、Ni吸附原子在石墨二炔表面的吸附能50。從穩(wěn)定性角度看，Sc、Ti原子是3d過渡金屬元素中最好的選擇。值得一提的是，Sc、Ti吸附原子在石墨二炔表面的吸附能，比在石墨烯表面吸附能(Ead＜ 2 eV)33,34,51大得多。因此，石墨二炔比石墨烯更適宜做過渡金屬原子的基底。

接下來，為了解Sc、Ti原子在石墨二炔表面的熱穩(wěn)定性，我們計(jì)算了它們?cè)谑脖砻孢w移的最小能量路徑以及勢(shì)壘E0。如圖4b中紅箭頭所示，Sc、Ti原子可從C位置，先跳到最近鄰的H位置再遷移至另一個(gè)C位置。也可以如圖4b中藍(lán)箭頭所示，越過炔鏈直接遷移至相鄰的C位置。根據(jù)計(jì)算結(jié)果，C→H→C的勢(shì)壘比C→C勢(shì)壘大兩倍，因此越過炔鏈的C→C遷移是Sc、Ti原子遷移的最小能量路徑。此外，Sc、Ti原子也可能會(huì)在炔鏈三角形大孔內(nèi)遷移。如果越過炔鏈的C→C勢(shì)壘足夠大，Sc、Ti原子只能停留在同一個(gè)炔鏈三角形內(nèi)部。計(jì)算得到的 Sc和 Ti原子的C→C擴(kuò)散勢(shì)壘 E0分別為 2.19和 1.74 eV (見表4)。一般來說，室溫下擴(kuò)散勢(shì)壘大于0.8 eV的反應(yīng)是難以進(jìn)行的。也就是說，在石墨二炔表面，足夠高的C→C擴(kuò)散勢(shì)壘能阻止Sc、Ti原子聚集。所以我們可以推斷，在室溫下Sc、Ti原子能很好地穩(wěn)定分散在石墨二炔表面。

表4 Sc、Ti原子分別在石墨二炔的C位置吸附能Ead(C)和H位置的吸附能Ead(H)，在C與C位置間的擴(kuò)散勢(shì)壘E0Table 4 The adsorption energy of Sc and Ti atoms on the C and H sites of graphdiyne (Ead(C) and Ead(H),respectively), and the migration barrier E0 between neighboring C sites.

研究HCHO在Sc、Ti摻雜石墨二炔上的吸附情況，是了解其作為分子傳感器的先決條件。我們隨機(jī)選取HCHO在摻Sc石墨二炔和摻Ti石墨二炔表面上各種初始的幾何構(gòu)型，通過幾何優(yōu)化來找到最穩(wěn)定的吸附構(gòu)型。石墨二炔-Sc和石墨二炔-Ti都表現(xiàn)出與HCHO分子的較強(qiáng)的相互作用。圖5展示出石墨二炔-Sc·HCHO和石墨二炔-Ti·HCHO 的最穩(wěn)定吸附構(gòu)型。在石墨二炔-Sc·HCHO構(gòu)型中，HCHO分子中的O原子與 Sc原子成鍵。在石墨二炔-Ti·HCHO構(gòu)型中，HCHO中的C原子和O原子同時(shí)與Ti原子成鍵。HCHO在石墨二炔-Sc和石墨二炔-Ti上最穩(wěn)定構(gòu)型的吸附能分別是2.59和2.24 eV，遠(yuǎn)大于在純石墨二炔上的吸附能 0.43 eV。相比之下，我們還研究了HCHO分子在石墨二炔-K、石墨二炔-Ca和石墨二炔-Cr體系上的吸附情況，吸附能分別是1.66、1.60和1.58 eV，均低于在石墨二炔-Sc和石墨二炔-Ti上的吸附能。相比之下，HCHO分子在摻雜Sc和Ti的石墨二炔上的吸附相當(dāng)強(qiáng)。此外，我們還研究了多個(gè)HCHO分子在石墨二炔-Sc和石墨二炔-Ti上的吸附。幾何優(yōu)化說明，在吸附一個(gè)HCHO分子后石墨二炔-Sc·HCHO 和石墨二炔-Ti·HCHO體系對(duì)于額外HCHO分子產(chǎn)生排斥效應(yīng)。從能量角度看，以上結(jié)果表示Sc、Ti吸附原子只能結(jié)合一個(gè)HCHO分子。這個(gè)特點(diǎn)為單分子的傳感提供了前提條件。

圖5 HCHO分子在石墨二炔-Sc和石墨二炔-Ti上吸附的側(cè)視及俯視圖Fig. 5 The side and top views of HCHO adsorption of graphdiyne-Sc and graphdiyne-Ti systems.

為了研究石墨二炔-Sc和石墨二炔-Ti體系對(duì)HCHO分子探測(cè)的靈敏度，我們比較這兩個(gè)體系了吸附HCHO分子前后，兩個(gè)體系的載流子濃度。在第一性原理計(jì)算中，導(dǎo)帶中的電子濃度n為

其中，能量 E的積分下限為導(dǎo)帶底(CBM)，D(E)為態(tài)密度，EF為費(fèi)米能量。價(jià)帶中的空穴濃度 p為

其中，能量 E的積分上限為價(jià)帶頂(VBM)?？傒d流子濃度N = n + p。為了較準(zhǔn)確地估算載流子濃度，各個(gè)體系的能帶均由雜化泛函HSE06的計(jì)算得到，如圖6所示。通常，PBE泛函的計(jì)算結(jié)果會(huì)低估半導(dǎo)體能帶，由雜化泛函HSE06計(jì)算的結(jié)果與PBE相比更精確。PBE泛函和HSE06泛函的計(jì)算結(jié)果見表5。

對(duì)于純凈石墨二炔的2 × 2超胞，能帶表現(xiàn)出本征半導(dǎo)體特性(即n = p)，由PBE和HSE06計(jì)算得到的Γ點(diǎn)直接帶隙分別為Eg= 0.48和0.89 eV。室溫300 K下由PBE和HSE06計(jì)算得到相應(yīng)的載流子濃度為 N = 3.96 × 108，1.29 × 105cm?2(見表5)。吸附HCHO分子后，體系帶隙略微減小到Eg= 0.46/0.88 eV (PBE/HSE06)，室溫下相應(yīng)的載流 子 濃 度 N = 4.71 × 108/1.47 × 105cm?2(PBE/HSE06)比吸附前的純凈石墨二炔高16%/12%，增加的載流子濃度N主要是由于HCHO分子導(dǎo)致的帶隙減小。

在石墨二炔-Sc及石墨二炔-Ti體系中，Sc、Ti吸附原子為石墨二炔貢獻(xiàn)了額外的電子，使導(dǎo)帶移至費(fèi)米面下使石墨二炔轉(zhuǎn)變成導(dǎo)體(圖6b、c)。室溫300 K下，石墨二炔-Sc體系載流子濃度為N =1.21 × 1013/1.30 × 1013cm?2(PBE/HSE06)，石墨二炔-Ti體系載流子濃度為 N=1.06 × 1013/1.23 × 1013cm?2(PBE/HSE06)。我們選取的 Sc、Ti與 2 × 2石墨二炔超胞體系的摻雜濃度c = 3.11 × 1013cm?2，相當(dāng)于每個(gè) Sc原子對(duì)石墨二炔貢獻(xiàn)了 0.42個(gè)電子，而每個(gè)Ti原子對(duì)石墨二炔貢獻(xiàn)了0.4個(gè)電子。電荷分析表明，石墨二炔上的 Sc原子及 Ti原子的凈電荷分別為+0.74e和+0.55e。以上結(jié)果展示了石墨二炔上的Sc、Ti吸附原子對(duì)電荷的貢獻(xiàn)。石墨二炔-Sc及石墨二炔-Ti體系比純石墨二炔體系載流子濃度大得多。

吸附一個(gè)HCHO分子后，石墨二炔-Sc·HCHO與石墨二炔-Ti·HCHO 體系中的導(dǎo)帶都相對(duì)石墨二炔-Sc及石墨二炔-Ti中的向上移動(dòng)(圖6b、c)。這一結(jié)果展示了HCHO分子的電負(fù)性。HCHO分子從石墨二炔中吸引電子，降低載流子濃度。當(dāng)溫度為300 K，石墨二炔-Sc·HCHO體系呈現(xiàn)出金屬性，載流子濃度N = 1.28 × 1013cm?2(HSE06泛函的結(jié)果)，略小于石墨二炔-Sc的(見表5)。對(duì)于石墨二炔-Ti·HCHO體系，導(dǎo)帶甚至移到費(fèi)米面上，體系由金屬性轉(zhuǎn)變?yōu)閚型半導(dǎo)體特性。在室溫300 K下，其載流子濃度N = 6.19 ×106cm?2(HSE06泛函的結(jié)果)，比石墨二炔-Ti體系低幾個(gè)量級(jí)，見圖6c。這表面HCHO分子幾乎全部帶走了Ti原子對(duì)石墨二炔的電子貢獻(xiàn)。所以，石墨二炔-Ti對(duì)于吸附 HCHO分子的電子響應(yīng)比石墨二炔-Sc更敏感且強(qiáng)得多。

表5 由PBE和HSE06兩種泛函計(jì)算得到的石墨二炔相關(guān)體系的載流子濃度N = n + p。Table 5 Carrier concentration N = n + p of graphdiyne systems at the level of PBE and HSE06.

圖6 采用HSE06泛函計(jì)算得到的純凈石墨二炔2 × 2超胞、石墨二炔-Sc、石墨二炔-Ti在吸附HCHO 前后的能帶圖Fig. 6 The energy bands of pristine 2 × 2 graphdiyne, graphdiyne-Sc and graphdiyne-Ti at the level of HSE06.

5 石墨二炔與氨基酸的相互作用

為了研究石墨二炔與氨基酸分子的相互作用，考慮到氨基酸的R基側(cè)鏈分為疏水、親水兩類。此外，氨基酸也按極性、非極性分類。在這里，我們選取甘氨酸(Gly)、谷氨酸(Glu)、苯丙氨酸(Phe)和組氨酸(His)作為每類氨基酸的代表。從氨基酸 R基側(cè)鏈的極性來看，Phe是非極性不帶電具有疏水性，Gly是極性不帶電具有親水性，His是極性帶正電具有親水性，Glu是極性帶負(fù)電具有親水性。從化學(xué)結(jié)構(gòu)來看，Gly及Glu為脂肪族氨基酸，Phe為芳香族氨基酸，His為雜環(huán)族氨基酸。所選氨基酸優(yōu)化后的結(jié)構(gòu)如圖7a所示。

在計(jì)算中，我們?nèi)匀徊捎? × 2的石墨二炔晶胞作為襯底模型。此外，為了與氨基酸在石墨烯上的吸附進(jìn)行比較，我們還將同樣的模擬方法用于7 × 7的石墨烯晶胞(邊長1.72 nm)。對(duì)石墨二炔和石墨烯單層，加上了2.0 nm的真空層。所選擇的2 × 2的石墨二炔晶胞和7 × 7的石墨烯晶胞的尺寸，能夠保證周期性盒子之間氨基酸分子的間距大于 0.9 nm。這樣我們可以認(rèn)為氨基酸分子是近獨(dú)立的。為了測(cè)試極化雙ζ基組的可靠性，我們以石墨烯為測(cè)試對(duì)象，用極化單ζ、極化雙ζ和極化叁 ζ基組進(jìn)行了結(jié)構(gòu)優(yōu)化，得到的碳碳鍵長分別為0.144、0.142和0.142 nm?？梢钥闯觯瑯O化雙ζ和極化叁 ζ基組已經(jīng)足夠合適。為了節(jié)省計(jì)算時(shí)間，在計(jì)算中我們采用的是極化雙ζ基組。為了研究范德華力的作用，我們還使用了vdW-DF52,53泛函與PBE-D2泛函進(jìn)行比較。

圖7 (a)計(jì)算中涉及的四種氨基酸分子，甘氨酸(Gly)、谷氨酸(Glu)、苯丙氨酸(Phe)和組氨酸(His)；(b)室溫下Gly分子在石墨二炔表面的遷移；(c)室溫下Gly分子在石墨烯表面的遷移；(d)分子動(dòng)力學(xué)模擬中的總能量變化；(e) Gly、Glu、His和Phe在石墨二炔表面的最穩(wěn)定吸附位置Fig. 7 (a) The structure of glycine (Gly), glutamic acid (Glu), histidine (His) and phenylalanine (Phe) molecules;(b) The migration of Gly on graphdiyne surface at room temperature; (c) The migration of Gly on graphene surface at room temperature; (d) The change of total energy in molecular dynamics simulations; (e) The most stable adsorption configurations of Gly, Glu, His and Phe on graphdiyne surface.

氨基酸在材料表面的吸附主要是由于范德華力的作用。范德華力分為取向力、誘導(dǎo)力和色散力三類。其中，取向力為極性分子與極性分子間的靜電力。誘導(dǎo)力為極性分子靠近非極性分子后，將非極性分子極化后產(chǎn)生的靜電力。色散力是一種特殊的力，它由偶極矩的量子漲落引起。一般的密度泛函方法(例如 PBE泛函36)不能很好的描述色散力。當(dāng)分子距離較遠(yuǎn)時(shí)色散力占主要地位，若泛函不能描寫量子漲落，計(jì)算出的分子間作用力會(huì)過快地趨于零。為了解決這個(gè)問題，人們發(fā)明了修正的非局域 vdW 泛函(如這里所使用的 vdW-DF)，但這種泛函的計(jì)算量比一般的 GGA泛函大。此外，近年來流行的 DFT-D方法(如本文中使用的PBE-D2)，將一般GGA泛函與經(jīng)驗(yàn)勢(shì)相結(jié)合，也是一種好的方案。下面我們采用PBE、vdW-DF泛函和 PBE-D2泛函計(jì)算氨基酸在石墨二炔表面的吸附能，并進(jìn)行比較，揭示相互作用力的本質(zhì)。

首先，我們以Gly為例，研究Gly分子與石墨二炔的相互作用，并與石墨烯進(jìn)行比較。我們采用了分子動(dòng)力學(xué)(MD)模擬的方法進(jìn)行了初探。初始時(shí)，我們將一個(gè)位置及取向隨機(jī)的 Gly分子置于石墨二炔表面進(jìn)行結(jié)構(gòu)優(yōu)化，然后在恒溫300 K下進(jìn)行了分子動(dòng)力學(xué)模擬。在模擬過程中，我們抽選某些體系勢(shì)能較低的瞬態(tài)進(jìn)行結(jié)構(gòu)優(yōu)化。從優(yōu)化后的幾何結(jié)構(gòu)中找出勢(shì)能最低，也就是 Gly吸附在石墨二炔表面的穩(wěn)定結(jié)構(gòu)。根據(jù)結(jié)果，最穩(wěn)定的石墨二炔吸附Gly復(fù)合體系如圖7e所示。在上述分子動(dòng)力學(xué)模擬中，我們觀察到在300 K溫度下，Gly分子持續(xù)在石墨二炔表面遷移并轉(zhuǎn)動(dòng)。一些結(jié)構(gòu)在展示與圖 7b。在室溫下，Gly分子遷移過程中與石墨二炔表面總是保持一定距離，但不會(huì)離開石墨二炔表面。在模擬中，總能量隨時(shí)間的波動(dòng)并保持在一定范圍內(nèi)(圖7d2)，該結(jié)果反映出體系處于平衡態(tài)。我們也進(jìn)行了800 K下的分子動(dòng)力學(xué)模擬，Gly分子仍然束縛于石墨二炔表面。這一結(jié)果證明 Gly在石墨二炔表面的吸附是穩(wěn)定的。

為了研究Gly在石墨二炔表面的最穩(wěn)定吸附，我們選取MD模擬過程中勢(shì)能較低的幾個(gè)構(gòu)型進(jìn)行優(yōu)化。對(duì)PBE、PBE-D2和vdW-DF泛函，所得到的Gly最大吸附能分別為0.59、0.90和1.10 eV(表6)。計(jì)算結(jié)果表明，PBE泛函給出的吸附能明顯地小于PBE-D2或vdW-DF泛函，這說明色散力在氨基酸吸附中占了很大一部分比例。PBE泛函對(duì)色散力的描寫并不太好，不能表示長程-C6/R6衰減(R為分子間距)。而PBE-D2與vdW-DF泛函結(jié)果接近，說明這里的計(jì)算結(jié)果較為可靠。對(duì)于各個(gè)不同的吸附位，PBE-D2泛函給出的Gly與石墨二炔的吸附能范圍是0.68–0.90 eV。由于Gly在石墨二炔表面的吸附構(gòu)型較為復(fù)雜，模擬中很難遍歷所有結(jié)構(gòu)，故以上吸附能是大致的估計(jì)。在最穩(wěn)定的吸附構(gòu)型(圖7e)中，Gly分子在石墨二炔表面位于大的碳環(huán)附近頂點(diǎn)。在 Gly分子中，由于 O原子與N原子具有強(qiáng)的負(fù)電性，O與N原子附近電子密度較大。當(dāng) Gly分子接近時(shí)，石墨二炔表面將被O或N原子極化，產(chǎn)生誘導(dǎo)力而被吸引。根據(jù)以上得到的Gly吸附構(gòu)型，―COOH與―NH2同時(shí)靠近石墨二炔表面時(shí)具有最大的吸附能(圖7e)。對(duì)于僅―COOH或―NH2一個(gè)靠近石墨二炔面的構(gòu)型，吸附能則較低。

表6 采用PBE、PBE-D2和vdW-DF三種不同的泛函得到的Gly、Glu、His和Phe在石墨二炔表面的吸附能(單位：eV)Table 6 The adsorption energy (in eV) of Gly, Glu, His and Phe on graphdiyne surface at the level of PBE,PBE-D2 and vdW-DF.

為了與石墨烯進(jìn)行比較，研究 Gly在石墨烯表面的吸附。我們采用了與石墨二炔吸附 Gly類似的方法進(jìn)行了模擬。模擬的情況與石墨炔類似，在模擬中總能量隨時(shí)間的波動(dòng)并保持在一定范圍內(nèi)(圖7d1)，結(jié)果反映出體系處于平衡態(tài)。圖7c展示了Gly在石墨烯表面運(yùn)動(dòng)的過程。對(duì)PBE、PBED2和vdW-DF泛函，所得到的Gly最大吸附能分別為0.23、0.53和0.54 eV (表7)。結(jié)果同樣表明，PBE泛函給出的吸附能明顯地小于 PBE-D2或vdW-DF泛函，色散力在氨基酸吸附中占了很大一部分比例。對(duì)PBE泛函，所得到的Gly在石墨烯表面最穩(wěn)定構(gòu)型的吸附能與文獻(xiàn)中的計(jì)算結(jié)果54一致。最穩(wěn)定的吸附構(gòu)型是N、O、C原子所在平面與石墨烯面平行的狀態(tài)。根據(jù)以上結(jié)果，Gly與石墨二炔的吸附明顯強(qiáng)于與石墨烯的吸附。與石墨烯中烯鍵相比，石墨二炔中的炔鍵具有額外的π電子可能是導(dǎo)致與氨基酸分子的強(qiáng)吸引。

表7 采用PBE、PBE-D2和vdW-DF三種不同的泛函得到的Gly、Glu、His和Phe在石墨烯表面的吸附能(單位：eV)Table 7 The adsorption energy (in eV) of Gly, Glu, His and Phe on graphene surface at the level of PBE,PBE-D2 and vdW-DF.

我們用同樣的方法研究Glu、His及Phe在石墨二炔表面的吸附。模擬溫度為300 K。從任意的Glu、His和Phe的構(gòu)型出發(fā)，從模擬中找出Glu、His和Phe勢(shì)能較低的瞬態(tài)構(gòu)型并進(jìn)行結(jié)構(gòu)優(yōu)化。從中找出吸附能最大的構(gòu)型。與Gly類似，在300 K的模擬中，Glu、His及Phe分子在石墨二炔表面持續(xù)遷移并轉(zhuǎn)動(dòng)，與表面保持0.2–0.3 nm的距離，但不離開石墨二炔表面。在模擬中，總能量隨時(shí)間在一定范圍內(nèi)波動(dòng)(圖7d3、d4、d5分別表示Glu、His和Phe的情況)。這種波動(dòng)反映出體系已處于平衡態(tài)。接下來，我們對(duì)分子動(dòng)力學(xué)模擬中能量低的瞬態(tài)結(jié)構(gòu)進(jìn)行優(yōu)化，計(jì)算吸附能，并與石墨烯比較。

對(duì)PBE、PBE-D2和vdW-DF泛函，所得到的Glu在石墨二炔表面最大吸附能分別為0.54、0.90和1.14 eV (表6)。該結(jié)果表明，PBE泛函給出的吸附能明顯地小于PBE-D2或vdW-DF泛函，色散力在氨基酸吸附中占了很大一部分比例。在石墨烯表面，PBE、PBE-D2和vdW-DF泛函所得到的Glu最大吸附能分別為0.43、0.84和0.94 eV (表7)。可以看出，石墨烯對(duì)Glu的吸附明顯弱于石墨二炔。對(duì)于在石墨二炔表面各種不同的位置和取向，PBE-D2泛函所得到的Glu分子與石墨二炔的吸附能分布值為0.82–0.90 eV。當(dāng)Glu接近石墨二炔六元環(huán)時(shí)，吸附能最大，相應(yīng)構(gòu)型如圖7e所示。與Gly類似的是，當(dāng)Glu分子的―COOH 和―NH2基團(tuán)接近石墨二炔表面時(shí)的它具有較大的吸附能。在室溫下，可以迅速發(fā)生在低能構(gòu)型和高能構(gòu)型間的轉(zhuǎn)變，因?yàn)樵觿?dòng)能足夠越過勢(shì)壘。

對(duì)His分子，PBE、PBE-D2和vdW-DF泛函所得到的在石墨二炔表面最大吸附能分別為0.73、1.22和 1.46 eV (表6)。而在石墨烯表面，PBE、PBE-D2和vdW-DF泛函所給出的His最大吸附能分別為0.45、1.00和1.14 eV (表7)?？梢钥闯?，石墨烯對(duì)His的吸附弱于石墨二炔。當(dāng)His分子的雜環(huán)接近石墨二炔的六元環(huán)并且―COOH和―NH2基團(tuán)靠近石墨二炔表面時(shí)，體系具有最大的吸附能，如7e所示。與Gly類似，作用力主要是―COOH和―NH2基團(tuán)產(chǎn)生的吸引力，這自來O原子與N原子強(qiáng)的負(fù)電性。PBE-D2泛函給出的 His分子在石墨二炔各種吸附構(gòu)型的吸附能范圍為1.18–1.22 eV。

對(duì)Phe分子，PBE、PBE-D2和vdW-DF泛函所得到的在石墨二炔表面最大吸附能分別為0.77、1.27和1.53 eV (表 6)。而在石墨烯表面，PBE、PBE-D2和vdW-DF泛函所給出的Phe最大吸附能分別為0.56、1.24和1.37 eV (表7)?？梢钥闯觯?duì)Phe的吸附弱于石墨二炔。PBED2泛函給出的Phe分子在石墨二炔各種吸附構(gòu)型的吸附能范圍為1.25–1.27 eV。Phe分子吸附于石墨二炔表面的最穩(wěn)定結(jié)構(gòu)是 Phe的芳香環(huán)與石墨二炔的六元環(huán)平行并形成AB式堆疊，如圖7e所示。這個(gè)最大吸附能比Gly，Glu及His的最大吸附能都大。這種芳香環(huán)間的相互作用應(yīng)來自π–π相互作用。

圖8 (a)石墨二炔吸附Gly體系的量子電子輸運(yùn)計(jì)算模型，陰影區(qū)為半無限大電極，y、z方向施加周期性邊界條件；(b)石墨二炔吸附氨基酸體系與純石墨二炔的I–Vb曲線Fig. 8 (a) Calculation model of quantum electronic transport in graphdiyne-Gly system. The shadow region denotes semi-infinite electrodes. Periodic boundary conditions are applied on y and z directions; (b) I–Vb curves of graphdiyne-amino-acid systems and pristine graphdiyne.

為了考慮用石墨二炔作為生物傳感器材料的可能性，我們研究了石墨二炔—氨基酸復(fù)合體系的介觀電子輸運(yùn)性質(zhì)。圖8a所示為石墨二炔-Gly體系的電子輸運(yùn)模擬示意圖，電極部分用陰影標(biāo)出，在y、z方向上施加了周期性邊界條件。散射區(qū)為電極原胞的 4倍大小，氨基酸分子位于散射區(qū)中間部分。在模擬中采用每種氨基酸分子最穩(wěn)定的吸附構(gòu)型。純石墨二炔、石墨二炔-Gly、石墨二炔-Glu、石墨二炔-His、石墨二炔-Phe的電流電壓曲線如圖8b所示，都呈現(xiàn)出半導(dǎo)體特性，導(dǎo)通電壓約為0.2 V。當(dāng)電壓在0.2 V以下，對(duì)于無論是否吸附氨基酸電流都趨近于零。當(dāng) Vb大于 0.2 V以后，電流隨電壓增大而增大。在同一Vb下，石墨二炔-氨基酸體系比純石墨二炔的電導(dǎo)率明顯要低。在電壓Vb在1.2 V時(shí)，石墨二炔-Gly，石墨二炔-Glu和純石墨二炔體系可由不同的電流值區(qū)分出來，然而石墨二炔—His與石墨二炔-Phe體系的電流值很接近。根據(jù)以上分析結(jié)果，不同的氨基酸可由特定電壓下的不同的電流值檢測(cè)出。這說明石墨二炔可被用于設(shè)計(jì)和制備生物傳感器。

6 結(jié)論

石墨二炔作為一種新型二維材料，具有獨(dú)特的物理性質(zhì)，可用作小分子傳感器的基礎(chǔ)材料。由于石墨二炔具有一定的化學(xué)惰性，需通過摻雜來改變其分子吸附特性，提高石墨二炔對(duì)分子的響應(yīng)能力。本文首先研究了石墨二炔納米帶的電子輸運(yùn)性質(zhì)。接下來，研究了通過摻雜3d金屬原子來調(diào)節(jié)石墨二炔物理性質(zhì)的方法。在3d金屬中，我們挑選出了在石墨二炔表面吸附能最大的Sc、Ti原子，發(fā)現(xiàn)摻雜Sc、Ti原子使石墨二炔能夠作為氣體分子檢測(cè)相應(yīng)器件。我們先從原子遷移勢(shì)壘的角度，確定了Sc、Ti摻雜石墨二炔在室溫下的穩(wěn)定性。接下來，選擇HCHO分子作為典型代表，從載流子濃度、介觀輸運(yùn)等方面探討了Sc、Ti摻雜石墨二炔對(duì)HCHO分子的響應(yīng)。計(jì)算結(jié)果表明，摻雜 Sc、Ti原子的石墨二炔可高效吸附HCHO。根據(jù)摻雜改性后石墨二炔吸附HCHO分子的電流—電壓曲線的改變，指出基于Sc、Ti摻雜石墨二炔可以設(shè)計(jì)二維氣體分子傳感器。

我們揭示氨基酸分子與石墨二炔的相互作用規(guī)律及物理機(jī)制，發(fā)現(xiàn)色散力在其相互作用中占主要成分。我們研究了氨基酸分子在石墨二炔表面的吸附，分析了氨基酸分子的吸附能和熱運(yùn)動(dòng)形式。在同樣電壓下，石墨二炔吸附不同氨基酸后的電流值不同。該結(jié)果指明了石墨二炔作為氨基酸分子電子學(xué)檢測(cè)器件的基礎(chǔ)。本研究對(duì)今后基于石墨二炔的納米電子學(xué)器件、生物傳感器、納米藥物載體等研究提供了重要參考依據(jù)和理論指導(dǎo)。