康大偉，孫孟樂，尤景漢，姬江濤

(1.河南科技大學(xué)a.物理與工程學(xué)院;b.車輛與動力工程學(xué)院，河南 洛陽471003;2.洛陽理工學(xué)院 數(shù)理部，河南 洛陽471023)

0 前言

隨著器件小型化和電路高集成度的發(fā)展，傳統(tǒng)的集成電路將走到其物理極限。為了克服這一難題，一門新的多領(lǐng)域交叉學(xué)科——分子電子學(xué)應(yīng)運而生。分子電子學(xué)是利用組裝和合成工藝直接從分子尺度上“從下到上”的設(shè)計電路。脫氧核糖核酸(DNA)作為一種生物大分子，由于其獨有的自組裝和自識別特性，特別適合分子電子學(xué)的要求，DNA分子中的電荷輸運性質(zhì)已經(jīng)成為生命科學(xué)、化學(xué)、物理學(xué)等多個交叉學(xué)科的研究熱點。隨著納米技術(shù)的發(fā)展，近年來有很多實驗直接將單個或一束DNA分子連接到金屬電極間，來測量它的電導(dǎo)。這些實驗可以更直接的揭示DNA的導(dǎo)電性質(zhì)，也為DNA作為分子器件的可能性進(jìn)行了直接的實驗驗證［1］。這類實驗的結(jié)果比較復(fù)雜，揭示DNA分子可能是絕緣體、半導(dǎo)體、導(dǎo)體甚至超導(dǎo)體［2－7］。導(dǎo)電性的多樣結(jié)果一是源于DNA本身的性質(zhì)，DNA的多種序列和分子形態(tài);二是來源于外界環(huán)境的影響:溫度，溶液中的離子，電極與DNA分子的接觸［8－11］。在直接測量實驗中，金屬電極與DNA分子界面處的良好接觸是實驗可重復(fù)性的一個重要因素，因此有必要研究界面的影響。DNA分子與電極間的連接是實驗結(jié)果可重復(fù)性的一個關(guān)鍵因素，降低界面處的接觸電阻是提高器件導(dǎo)電特性的一個有效途徑。最近的實驗中［12］，DNA分子器件中出現(xiàn)了負(fù)微分電導(dǎo)現(xiàn)象，實驗結(jié)果發(fā)現(xiàn)通過DNA分子的電流隨著電壓的升高反而出現(xiàn)了下降的趨勢。在以前的研究中，往往假設(shè)電勢全部降落在DNA分子與電極的界面處，而在DNA分子上電勢降落為零。如果考慮DNA分子上有電勢降落，負(fù)微分電導(dǎo)的出現(xiàn)是否與界面處的電勢降落有關(guān)需要進(jìn)一步的研究。本文擬采用格林函數(shù)方法和Landauer-Büttiker公式來研究不同的界面電勢降落情況下DNA分子器件中的電荷輸運性質(zhì)。

1 模型和方法

對于金屬電極連接的DNA分子器件(見圖1)，選取1個螺旋周期(包含11個基對)的Poly(G)-Poly(C)DNA分子。整個體系的哈密頓可以用下式來表示

其中Hmol描述了Poly(G)-Poly(C)DNA分子的π軌道。

圖1 DNA分子器件示意圖

金屬電極的哈密頓由下式給出

其中，ε0和 τ分別是電極的在位能和相鄰格點的躍遷積分是電極中在格點(i，j)位置上產(chǎn)生或消滅一個自旋為σ的電子的算符。

Hcoup描述了電極與DNA分子間的耦合，這個耦合只發(fā)生在金屬電極和DNA分子的接觸界面上

其中假定界面耦合是自旋非相關(guān)的。

通過體系的電流就可以通過Landauer-Buttiker公式給出，

其中，h是普朗克常數(shù);T(E)是電子從左電極到右電極的透射系數(shù)

是電極中的費米分布函數(shù)，其中μl(r)是左(右)電極的化學(xué)勢，kB是玻耳茲曼常數(shù)，T是溫度。透射系數(shù)T(E)可以通過Fisher-Lee關(guān)系從體系的晶格格林函數(shù)得出，

其中I是單位矩陣算符。

本文中參數(shù)的選取如下:對于Poly(G)-Poly(C)DNA分子，將其最低占據(jù)能級(LUMO)設(shè)為在位能:εG=1.14 eV，εC=－1.06 eV。一條鏈上相鄰堿基間的躍遷積分設(shè)為tG(C)=0.45 eV，鏈間的互補堿基間的躍遷積分設(shè)為t⊥=0.60 eV，這些都是在參數(shù)選取的合理區(qū)間。金屬電極參數(shù)的選取為:τ=1.50 eV。電極的費米面選在DNA帶隙的中間。界面耦合設(shè)為自旋非相關(guān)的，取值為0.40 eV。本文選取的溫度為室溫T=300 K。

圖2 DNA分子器件在不同的電勢降落下的電壓電流曲線

2 結(jié)果與討論

在金屬電極連接的DNA分子器件電荷輸運性質(zhì)的直接實驗測量中，分子與電極間良好的界面接觸是實驗數(shù)據(jù)可重復(fù)性測量的必要條件。當(dāng)DNA分子與金屬電極形成良好的接觸后，在金屬電極施加外加偏壓，這時候左右兩個電極處在高低兩個不同電勢下。在以往的研究中，往往假設(shè)電勢對稱的降落在DNA分子與金屬電極的界面處，具體情況為左右兩個界面處各降落V/2。這種情況下系統(tǒng)的輸運性質(zhì)如圖2中所示。當(dāng)電壓比較小時，體系的電流基本為零，此時金屬電極的費米面處于DNA分子禁帶內(nèi)，因此沒有分子軌道進(jìn)入導(dǎo)電窗口。當(dāng)電壓大于0.7 V時，電流開始顯著增大，對應(yīng)著DNA分子的導(dǎo)電能級開始進(jìn)入導(dǎo)電窗口。隨著電壓的進(jìn)一步增大，更多的導(dǎo)電分子能級進(jìn)入導(dǎo)電窗口，這導(dǎo)致電流進(jìn)一步增大。因此，此時的電壓電流曲線呈現(xiàn)出非線性的單調(diào)遞增的趨勢。由電流隨著電壓的上升趨勢可以看出:在某個導(dǎo)電能級進(jìn)入導(dǎo)電窗口時，電流上升趨勢很快;當(dāng)電極費米面處在兩個導(dǎo)電能級間時，電流隨著電壓的上升趨勢較緩。體系通過的電流在電壓為4 V時達(dá)到了50 μA，此結(jié)果與最近的實驗［12］符合。

在實際的實驗中，如果分子與電極的界面處電勢降落小于電極上施加的電壓，則有部分電壓會降落在DNA分子上。電勢降落不只是發(fā)生在界面處的情況下，圖2的虛線和點劃線分別給出了DNA分子中的電勢降落達(dá)到總電勢降的1/3和1/2兩種情況下的電壓電流曲線。由圖2可以看出:DNA分子上的電勢降落壓制了通過其中的電流。在低電壓下，電流將小于無電勢降落的情況。而且隨著電壓的繼續(xù)升高，當(dāng)電壓大于3 V時電流反而出現(xiàn)了下降的情況，即出現(xiàn)了負(fù)微分電導(dǎo)現(xiàn)象。當(dāng)DNA分子上的總電勢降落所占的比重更大時，即電勢有一半降在分子上時，電流進(jìn)一步降低，而且在電壓大于2.6 V時就開始出現(xiàn)負(fù)微分電導(dǎo)現(xiàn)象。所以在DNA分子上的電勢降落越大，體系的電流越小。而在相同的電勢降落下，隨著電壓的升高，都會出現(xiàn)負(fù)微分電導(dǎo)現(xiàn)象。對于DNA分子中電壓降落Vdrop=V/3這種情況，當(dāng)電壓達(dá)到3 V左右時，電流開始隨著電壓的升高而降低。而當(dāng)Vdrop=V/2時，當(dāng)電壓達(dá)到2.5 V左右時，電流就會隨著電壓的升高而降低?？梢姡S著電勢降落在DNA分子上的比例的升高，其負(fù)微分電導(dǎo)出現(xiàn)的電壓值越低。

為了具體分析產(chǎn)生負(fù)微分電導(dǎo)的原因，畫出了在不同的電勢降落下，LUMO能級(最低未占據(jù)能級)對應(yīng)的波函數(shù)的形狀。如圖3a所示，其中縱坐標(biāo)為波函數(shù)的數(shù)值，橫坐標(biāo)為DNA分子堿基對應(yīng)的坐標(biāo)。當(dāng)DNA分子中不存在電壓降落時，LUMO能級對應(yīng)的波函數(shù)是擴展在整個分子中的，而擴展的分子軌道可以提供較大的電子透射幾率，此時通過體系的電流較大。當(dāng) DNA分子中存在電勢降落時，如圖3b所示，LUMO能級對應(yīng)的波函數(shù)局域在DNA分子的右端，局域軌道阻礙了電子的隧穿，從而降低了通過體系的電流。負(fù)微分電導(dǎo)的出現(xiàn)可以解釋為:隨著偏壓的增大，DNA分子上的電勢降落持續(xù)增大，因而各個導(dǎo)電的分子軌道傾向于局域化，這個效應(yīng)是減少電流的趨勢。另一方面，電壓增大導(dǎo)致更多的導(dǎo)電通道進(jìn)入導(dǎo)電窗口，這是電流增加的趨勢。當(dāng)兩方面的因素相互達(dá)到平衡時，電流隨著電壓的升高開始出現(xiàn)拐點，當(dāng)電壓繼續(xù)增大時，分子軌道局域化的效應(yīng)開始壓制更多導(dǎo)電通道進(jìn)入導(dǎo)電窗口的效應(yīng)，因此，電流開始隨著電壓的升高而下降。

如果電壓降落在DNA分子上的比例相對于降落在界面處的比例很小時，波函數(shù)的局域效應(yīng)較弱，而隨著電壓升高更多的導(dǎo)電通道進(jìn)入導(dǎo)電窗口，這個效應(yīng)占主要，所以電流隨著電壓升高會持續(xù)上升。在早期的DNA分子器件導(dǎo)電實驗中［13－15］，是通過物理吸附將DNA分子與電極連接，DNA分子與金屬電極間的連接不是化學(xué)連接，所以接觸電阻較大。此時，電勢將主要降落在界面處。如果界面處接觸很好，電壓將被迫降落到DNA分子上，當(dāng)電壓降的比例升高到一定程度，就會出現(xiàn)負(fù)微分電導(dǎo)的現(xiàn)象。

圖3 LUMO能級的波函數(shù)

3 結(jié)論

研究了DNA分子器件中電勢降落分布對電荷輸運性質(zhì)的影響。當(dāng)電勢降落全部在界面時，電流隨著電壓升高是持續(xù)上升的;當(dāng)電勢開始降落在DNA分子上時，電流相比DNA分子上沒有電勢降落的情況會有所降低，而且開始出現(xiàn)負(fù)微分電導(dǎo)的現(xiàn)象;當(dāng)電勢降落在DNA分子上的比例進(jìn)一步提高時，負(fù)微分電導(dǎo)出現(xiàn)時對應(yīng)的電壓值更小，而且下降的趨勢更明顯。分析了產(chǎn)生負(fù)微分電導(dǎo)的原因，發(fā)現(xiàn)電勢降落引起的分子軌道的局域是主要因素，電勢降落越大，分子局域現(xiàn)象越明顯，這就進(jìn)一步促使了負(fù)微分電導(dǎo)的出現(xiàn)。

［1］Endres R G，Cox D L，Singh R R P.Colloquium:The Quest for High-conductance DNA［J］.Rev Mod Phys，2004，76:195 －214.

［2］Fink H W，Schonenberger C.Electrical Conduction Through DNA Molecules［J］.Nature，1999，398:407 －410.

［3］Porath D，Bezryadin A，Vries S de，et al.Direct Measurement of Electrical Transport Through DNA Molecules［J］.Nature，2000，403:635 －638.

［4］Roy S，Vedala H，Roy A D，et al.Direct Electrical Measurements on Single-molecule Genomic DNA Using Single-walled Carbon Nanotubes［J］.Nano Lett，2008(8):26 － 30.

［5］Song B，Elstner M，Cuniberti G.Anomalous Conductance Response of DNA Wires Under Stretching［J］.Nano Lett，2008(8):3217－3220.

［6］Xu B Q，Zhang P M，Li X L，et al.Direct Conductance Measurement of Single DNA Molecules in Aqueous Solution［J］.Nano Lett，2004(4):1105 －1108.

［7］Yoo K H，Ha D H，Lee J O，et al.Electrical Conduction Through Poly(dA)-Poly(dT)and Poly(dG)-Poly(dC)DNA Molecules［J］.Phys Rev Lett，2001，87(19):198102.

［8］Cuniberti G，Craco L，Porath D，et al.Backbone-induced Semiconducting Behavior in Short DNA Wires［J］.Phys Rev B，2002，65:241314.

［9］Guo A M，Xiong S J.Effects of Contact and Efficient Charge Transport in G4-DNA Molecules［J］.Phys Rev B，2009，80:035115.

［10］Maciá E，Triozon F，Roche S.Contact-dependent Effects and Tunneling Currents in DNA Molecules［J］.Phys Rev B，2005，71:113106.

［11］Wang X F，Chakraborty T.Charge Transfer Via a Two-Strand Superexchange Bridge in DNA［J］.Phys Rev Lett，2006，97:106602.

［12］Ning K，Artur E，Elke S.Observation of Negative Differential Resistance in DNA Molecular Junctions［J］.Appl Phys Lett，2010，96:023701.

［13］Cohen H，Nogues C，Naaman R，et al.Direct Measurement of Electrical Transport Through Single DNA Molecules of Complex Sequence［J］.Proc Nat Acad Sci，2005，102:11589.

［14］Ullien D，Cohen H，Porath D.The Effect of the Number of Parallel DNA Molecules on Electric Charge Transport Through‘Standing DNA’［J］.Nanotechnology，2007，18:424015.

［15］Nogues C，Cohen S R，Daube S，et al.Sequence Dependence of Charge Transport Properties of DNA［J］.J Phys Chem B，2006，110:8910.1a21a1b