Rh(111)表面NO分子對(duì)多層膜的原子結(jié)構(gòu)?

汪辰超 吳太權(quán) 王新燕 江影

(中國(guó)計(jì)量大學(xué)應(yīng)用物理系,杭州 310018)

Rh(111)表面NO分子對(duì)多層膜的原子結(jié)構(gòu)?

汪辰超 吳太權(quán)?王新燕 江影

(中國(guó)計(jì)量大學(xué)應(yīng)用物理系,杭州 310018)

(2016年9月1日收到;2016年11月3日收到修改稿)

利用第一性原理研究了NO分子對(duì)[(NO)2]分子鏈、分子單層膜,Rh(111)表面上的(NO)2分子單層膜和多層膜的原子結(jié)構(gòu).(NO)2分子單體在虛擬Rh(111)表面自組裝成兩個(gè)穩(wěn)√定的分子鏈,(NO)2分子平行有序排列,氧原子和氮原子都呈現(xiàn)(100)和(111)結(jié)構(gòu).在虛擬Rh(111)-(1×3)上,1.00 ML(molecular layer)覆蓋度時(shí),(NO)2分子自組裝成兩個(gè)穩(wěn)定的分子單層膜(M1和M2),分子膜M1中N-N鍵與襯底的夾角為70?-90?;分子膜M2中N-N鍵平行襯底.在M2/Rh(111)中,(NO)2分子可吸附于頂位、fcc空心位和hcp空心位,通過電荷轉(zhuǎn)移可解釋兩個(gè)空心位的穩(wěn)定性強(qiáng)于頂位.Rh(111)表面(NO)2分子多層膜系統(tǒng)中,(NO)2分子垂直吸附于兩個(gè)空心位,第一層是分子膜M2,N-N鍵平行于襯底,第二層及以上都是分子膜M1,N-N鍵與襯底夾角為70?-90?,分子膜真空層為0.31 nm±0.02 nm.

一氧化氮分子對(duì),單層膜,多層膜,自組裝

1 引 言

分子自組裝(molecular self-assembly)是分子在熱力學(xué)平衡條件下,通過化學(xué)鍵和弱的非共價(jià)鍵的協(xié)同作用自發(fā)地形成有序的結(jié)構(gòu)[1].通過自組裝形成的有序結(jié)構(gòu)是分子間相互作用力和無序性兩種作用機(jī)制相互競(jìng)爭(zhēng),在一定的條件下達(dá)到平衡的結(jié)果.無序的分子在自發(fā)形成有序結(jié)構(gòu)的過程可從外界環(huán)境中引入負(fù)熵,使得自身變得有序.分子自組裝膜是構(gòu)膜分子通過分子間及其與基底材料間物化作用而自發(fā)形成的一種熱力學(xué)穩(wěn)定,排列規(guī)則的單層(或多層)分子膜.在金屬表面處理、生物工程、微電子學(xué)、非線性光學(xué)、分子器件、傳感器件等眾多的領(lǐng)域都有廣泛的應(yīng)用[2-5],是近年來許多領(lǐng)域的科學(xué)家們所研究的熱點(diǎn).

一氧化氮(NO)是大氣污染物中的主要成分之一[6],NO分子很容易形成NO分子對(duì)[(NO)2]分子[7]. Rh在異相催化領(lǐng)域里是眾多催化劑中最有效的催化劑之一,它能催化CO,NO等來控制汽車尾氣,具有重要的工業(yè)價(jià)值[8-14].1986年,Root等[8]利用電子能量損失譜研究了NO/Rh(111)系統(tǒng),結(jié)果表明NO分子吸附于橋位.1998年,Loffreda等[9]利用周期密度泛函理論研究了NO/Rh(111)系統(tǒng),指出fcc空心位和hcp空心位更穩(wěn)定.2006年3月,Wallace等[10]利用偏振調(diào)制紅外反射吸收譜研究了NO/Rh(111)系統(tǒng),結(jié)果表明NO分子可吸附于頂位和兩個(gè)空心位.同年8月,Nakamura等[11]利用紅外反射吸收譜研究了NO/Rh(111)系統(tǒng),其結(jié)果與文獻(xiàn)[10]一致,即NO分子可吸附于頂位和兩個(gè)空心位,但是文獻(xiàn)[8-11]均未涉及(NO)2分子的研究.

2007年,Nakai等[12]利用近邊X射線吸收精細(xì)結(jié)構(gòu)和從頭計(jì)算方法研究了(NO)2/Rh(111)系統(tǒng),其吸附位的結(jié)論與文獻(xiàn)[10,11]一致,并指出Rh(111)表面上分子第一層膜中NO分子不成對(duì),NO分子垂直吸附于表面的頂位和兩個(gè)空心位,分子第二層膜中NO分子開始成對(duì),N-N鍵和(NO)2分子平面都是傾斜于表面,分子多層膜中(NO)2分子的N-N鍵垂直于表面.文獻(xiàn)[12]中給出了(NO)2分子的N-N鍵約為0.230 nm,N-N-O角約為96.0°(記為∠NNO=90°),但是未分析(NO)2分子間的相互作用,也未考慮(NO)2分子在吸附過程中結(jié)構(gòu)的變化情況.2008年,Jansen和Popa[13]利用貝葉斯統(tǒng)計(jì)和密度泛函理論研究了不同覆蓋度的NO/Rh(111)系統(tǒng),其吸附位的結(jié)論與文獻(xiàn)[10-12]一致,但是文獻(xiàn)[13]強(qiáng)調(diào)吸附位與覆蓋度有關(guān),在覆蓋度為0.5 ML(molecular layer)時(shí)NO分子主要吸附在頂位.

2012年,Wu等[14]利用Castep模塊研究了(NO)2/Rh(111)系統(tǒng),給出了(NO)2單體、分子鏈、單層膜和吸附系統(tǒng)的結(jié)構(gòu)參數(shù),其中分子穩(wěn)定的吸附位為兩個(gè)空心位和頂位,該結(jié)論與文獻(xiàn)[10-13]一致.文獻(xiàn)[14]中討論了(NO)2分子間的相互作用,并強(qiáng)調(diào)研究吸附系統(tǒng)時(shí)需考慮分子單層膜的整體吸附,同時(shí)指出了(NO)2分子的結(jié)構(gòu)參數(shù)在單體、分子鏈、單層膜和吸附系統(tǒng)中可取不同的數(shù)值,如當(dāng)∠NNO=90°時(shí),(NO)2分子的N-N鍵取值范圍為0.19-0.31 nm.文獻(xiàn)中計(jì)算(NO)2分子鏈時(shí)未考慮Rh(111)結(jié)構(gòu),解離Rh(111)面之前未優(yōu)化Rh晶體,也未給出(NO)2分子單層膜在Rh(111)表面的自組裝過程,還未分析(NO)2分子多層膜的原子結(jié)構(gòu)與相互作用.

1995年,Brown等[15]利用反射吸收紅外光譜研究了低溫下的(NO)2/Ag(111)系統(tǒng),指出Ag(111)表面上(NO)2分子單層膜中N-N鍵平行于表面,而(NO)2分子多層膜中N-N鍵垂直于表面.而文獻(xiàn)[12]中(NO)2/Rh(111)系統(tǒng)第二層(NO)2分子膜與其他分子層不同.因此本文在文獻(xiàn)[14]的基礎(chǔ)上,以(NO)2分子為基本結(jié)構(gòu)單元,研究Rh(111)表面上(NO)2分子多層膜的結(jié)構(gòu)參數(shù)及多層膜內(nèi)的相互作用.

2 模型參數(shù)

CASTEP(ambridge sequential total energy package)是一個(gè)基于密度泛函方法的第一性原理的程序.本文中原子結(jié)構(gòu)的計(jì)算基于密度泛函理論.交換關(guān)聯(lián)泛函采用的是最流行的Perdew-Burke-Ernzerhof廣義梯度近似[16].核心電子和價(jià)電子之間的作用采用超軟贗勢(shì)描述.

在本文中,(NO)2分子單體的結(jié)構(gòu)見圖1.該原胞為正方體,其長(zhǎng)寬高標(biāo)記為a,b,c,其中N-O間距標(biāo)記為L(zhǎng)NO,N-N間距標(biāo)記為L(zhǎng)NN,O-O間距標(biāo)記為L(zhǎng)OO,夾角∠NNO標(biāo)記為α.考慮到(NO)2分子單體的平面結(jié)構(gòu)尺寸(0.120 nm×0.251 nm),為了忽略分子間的相互作用,分子單體N-N鍵平行于原胞底面的對(duì)角線.此時(shí)原胞邊長(zhǎng)a,b,c均取值1.400 nm.

圖1 (網(wǎng)刊彩色)(NO)2分子單體的結(jié)構(gòu)示意圖Fig.1. (color online)The sketch map of(NO)2monomer.

圖2 Rh(111)-(1×)結(jié)構(gòu)示意圖 1,頂位;2,fcc空心位;3,橋位;4,hcp空心位Fig.2.The sketch map of Rh(111)：1,top site;2,fcc hollow site;3,bridge site;4,hcp hollow site.

當(dāng)(NO)2分子單體吸附于Rh(111)表面時(shí),通?；瘜W(xué)吸附能ΔEchem表示為

其中E(NO)2/Rh(111)為(NO)2/Rh(111)吸附系統(tǒng)的能量,ERh(111)為Rh(111)襯底的能量,E(NO)2monomer為(NO)2分子單體的能量.

當(dāng)(NO)2分子單體組成團(tuán)簇(如分子鏈和單層膜)時(shí),結(jié)合能ΔEb表示為

其中E(NO)2molecular chain和E(NO)2monolayer分別是(NO)2分子鏈和單層膜的能量.

因?yàn)槲接赗h(111)表面的(NO)2分子單體容易組成了單層膜,則化學(xué)吸附能附能應(yīng)改寫為

吸附系統(tǒng)中Rh(111)表面上分子多層膜間的相互作用可表示為

其中E(n+1)monolayer/Rh(111)是n+1層單層膜在Rh(111)表面吸附的能量,而Enmonolayer/Rh(111)是n層單層膜在Rh(111)表面吸附的能量.

由于分子單層膜的能量低于分子單體,對(duì)于分子自組裝吸附系統(tǒng),必須考慮單層膜的整體吸附,而不能只簡(jiǎn)單考慮單體與襯底的作用.要考慮單層膜的整體吸附,必須引入虛擬表面.虛擬表面和第一性原理已被廣泛應(yīng)用于確定金屬表面分子自組裝膜的原子結(jié)構(gòu)[14,17-21],由這些文獻(xiàn)可知分子膜的原子結(jié)構(gòu)主要由相鄰分子決定,吸附襯底對(duì)其的影響較弱.

在計(jì)算過程中,為了獲得滿意的結(jié)果我們采用了截止能量為380.0 eV的贗勢(shì)波.精細(xì)精度的收斂判據(jù)如下：Hellmann-Feynaman力小于0.001 eV/nm;能量變化小于0.000005 eV;位移變化小于0.00005 nm.對(duì)結(jié)構(gòu)模型結(jié)構(gòu)優(yōu)化過程中k點(diǎn)取為精細(xì),具體參數(shù)列于表1.所有的計(jì)算都采用Materials Studio 8.0環(huán)境下[22]的CASTEP模塊[23,24].

表1 (NO)2單體、Rh晶體、Rh(111)面、分子鏈、單層膜和吸附系統(tǒng)的結(jié)構(gòu)參數(shù)Table 1.The structural parameters of(NO)2monomer,Rh crystal,Rh(111)surface,molecular chains,monolayers and(NO)2/Rh(111)system.

3 結(jié)果與討論

3.1 虛擬Rh(111)表面(NO)2分子鏈的原子結(jié)構(gòu)

在文獻(xiàn)[14]的基礎(chǔ)上對(duì)(NO)2分子單體進(jìn)行結(jié)構(gòu)優(yōu)化,獲得分子單體內(nèi)鍵長(zhǎng)LNO為0.120 nm,LNN為0.199 nm,LOO為0.251 nm,鍵角α為103°,分子單體能量為-1408.805 eV.以優(yōu)化后的分子單體為基構(gòu)建兩種最可能的分子鏈(C1和C2),相鄰(NO)2分子中O-O間距標(biāo)記為D,如圖3所示.改變分子間距D,從1.000 nm減小到0.220 nm,對(duì)這些分子鏈進(jìn)行結(jié)構(gòu)優(yōu)化,結(jié)果列于表2,“-”表示該分子鏈結(jié)構(gòu)不穩(wěn)定,對(duì)應(yīng)的結(jié)合能曲線繪于圖4.

圖3 (網(wǎng)刊彩色)(NO)2分子鏈C1和C2結(jié)構(gòu)示意圖Fig.3. (color online)The sketch map of(NO)2molecular chains C1 and C2.

表2 隨分子間距D從1.000 nm減小到0.220 nm時(shí)(NO)2分子鏈的能量Table 2.Energy of(NO)2molecular chains as distanceDfrom 1.000 nm to 0.220 nm.

比較表2和圖4中兩種分子鏈的結(jié)合能,可知當(dāng)分子間距D取值在0.240-0.380 nm區(qū)間時(shí)分子鏈C1結(jié)構(gòu)穩(wěn)定,而不論分子間距D取值如何,分子鏈C2結(jié)構(gòu)都不穩(wěn)定.由文獻(xiàn)[12,14]可知(NO)2分子為平面分子,呈梯形狀,氮原子呈弱正電性,氧原子呈弱負(fù)電性,兩個(gè)氮原子和兩個(gè)氧原子間都是弱排斥,氧原子和氮原子間是弱吸引,當(dāng)弱吸引和弱排斥平衡時(shí),分子單體結(jié)構(gòu)穩(wěn)定.在兩種分子鏈中逐漸減小分子間距D,當(dāng)間距較大,相鄰(NO)2分子作用力很弱,結(jié)合能約為零,分子鏈能量約等于分子單體的能量.在分子鏈C1中,當(dāng)分子間距D≤0.380 nm時(shí),相鄰分子內(nèi)的排斥開始小于吸引,分子鏈C1開始趨于穩(wěn)定;當(dāng)分子間距D≤0.220 nm(小于分子單體內(nèi)間距LOO)時(shí),(NO)2分子的結(jié)構(gòu)被破壞,分子鏈C1結(jié)構(gòu)不穩(wěn)定.在分子鏈C2中,當(dāng)分子間距D≤0.280 nm時(shí),相鄰分子內(nèi)的氧原子的排斥使得分子鏈C2結(jié)構(gòu)不穩(wěn)定.在分子鏈C1中當(dāng)分子間距取值為0.280 nm±0.040 nm時(shí),分子單體內(nèi)LNO為0.120 nm,LNN為0.193 nm,LOO為0.278 nm,鍵角α為110°.該數(shù)據(jù)和文獻(xiàn)[14]一致.

圖4 分子間距D從1.000 nm減小到0.220 nm時(shí)(NO)2分子鏈的結(jié)合能ΔEb1Fig.4.Energy ΔEb1of(NO)2molecular chain as distanceDfrom 1.000 nm to 0.220 nm.

表3 (NO)2分子鏈C3-a和C3-b的能量與結(jié)構(gòu)參數(shù)Table 3.Energy and structural parameters of(NO)2molecular chains C3-a and C3-b.

圖5 (網(wǎng)刊彩色)(NO)2分子鏈C3-a和C3-b的結(jié)構(gòu)示意圖Fig.5.(color online)The sketch maps of(NO)2molecular chains C3-a and C3-b.

比較表3中LOO和分子間距D,它們幾乎相等,即分子鏈C3-a和C3-b中,氧原子和氮原子均有序排列,呈現(xiàn)(100)和(111)結(jié)構(gòu).分子鏈C3-a中分子單體結(jié)構(gòu)特征明顯,而分子鏈C3-b中分子單體的每個(gè)NO分子都是等價(jià)的.和分子單體結(jié)構(gòu)相比較,分子鏈C3-a中,LNO保持不變,LNN微減小3%,LOO增大了11%,鍵角α增大了8%,而分子鏈C3-b中,LNO保持不變,LNN增大了19%,LOO增大了12%,鍵角α減小了3%.這意味著(NO)2分子結(jié)構(gòu)在不同的團(tuán)簇中可取不同的值,且比較(NO)2分子單體的能量和NO單分子的能量可知(NO)2分子是否存在[14].由結(jié)合能可知,分子鏈?zhǔn)且粋€(gè)自組裝系統(tǒng).

3.2 虛擬Rh(111)表面(NO)2分子單層膜的原子結(jié)構(gòu)格參數(shù)中的

穩(wěn)定的分子鏈C3中相鄰分子單體間距(0.275 nm)對(duì)應(yīng)Rh晶格參數(shù)中1.由于相鄰分子鏈最可能的作用距離為0.300-0.500 nm,則Rh晶為合適的相鄰分子間距,即分子單層膜的襯底結(jié)構(gòu)為覆蓋度為1.00 ML的令(NO)2分子單體中N-N鍵與表面的夾角為γ.把分子鏈C3-a和C3-b吸附于虛擬Rh(111)表面構(gòu)建出分子單層膜,改變夾角γ從0°每隔10°增大到90°,對(duì)這些分子膜進(jìn)行結(jié)構(gòu)優(yōu)化,結(jié)果可知分子鏈C3-a和C3-b分別組成一個(gè)穩(wěn)定的分子膜(M1和M2).由于當(dāng)γ取值在70°-90°之間時(shí),分子膜的能量相等,且結(jié)構(gòu)參數(shù)一致,故把這些分子膜合起來標(biāo)記為M1.分子膜M2中,NO分子幾乎垂直虛擬表面,而N-N鍵平行虛擬襯底.分子膜M1和M2的結(jié)構(gòu)示意圖見圖6.

圖6 (網(wǎng)刊彩色)(NO)2分子單層膜M1和M2的結(jié)構(gòu)示意圖Fig.6.(color online)The sketch maps of(NO)2monolayers M1 and M2.

表4 (NO)2分子單層膜的能量與結(jié)構(gòu)參數(shù)Table 4.Energy and structural parameters of(NO)2monolayers.

分子膜M1由許多旋轉(zhuǎn)的分子鏈C3-a組成,因其能量略大于分子鏈C3-a的能量(約0.01 eV),則在分子膜M1中相鄰分子鏈呈弱吸引作用.分子膜M2由分子鏈C3-b組成,因其能量等于分子鏈C3-b的能量,則在分子膜M2中相鄰分子鏈之間的相互作用約為零.

3.3 Rh(111)表面(NO)2分子單層膜的原子結(jié)構(gòu)

比較三個(gè)結(jié)構(gòu)的能量,可知fcc空心位的穩(wěn)定性稍強(qiáng)于hcp空心位,空心位的穩(wěn)定性強(qiáng)于頂位.已知氮原子有五個(gè)核外電子,其中兩個(gè)電子與氧原子形成共價(jià)鍵,剩下的三個(gè)電子對(duì)應(yīng)三個(gè)襯底原子時(shí)比對(duì)應(yīng)一個(gè)襯底原子合適.當(dāng)NO分子吸附于hcp空心位時(shí),第二層的Rh原子對(duì)N原子有弱排斥作用.因此(NO)2分子膜M2在Rh(111)表面吸附時(shí),fcc空心位稍強(qiáng)于hcp空心位,他們都強(qiáng)于頂位,而橋位不穩(wěn)定.根據(jù)CASTEP計(jì)算結(jié)果,(NO)2分子單體中NO的電子結(jié)構(gòu)為N+0.16=O-0.16,N原子轉(zhuǎn)移了部分電子給O原子.而分子吸附于頂位時(shí),其電子結(jié)構(gòu)為Rh+0.14-N0=O-0.14,即N原子轉(zhuǎn)移給O原子的部分電子等于N原子從Rh原子獲得的部分電子.當(dāng)分子吸附于兩個(gè)空心位時(shí),其電子結(jié)構(gòu)為Rh+0.34-N-0.18=O-0.16,即N原子從Rh原子獲得的部分電子除了轉(zhuǎn)移給O原子的部分電子還剩下部分電子.根據(jù)電荷轉(zhuǎn)移情況可知空心位吸附時(shí)電荷相互作用強(qiáng)于頂位吸附時(shí),即空心位吸附比頂位更穩(wěn)定,該結(jié)論和文獻(xiàn)[9-14]的結(jié)論一致.

圖7 (網(wǎng)刊彩色)

3.4 Rh(111)表面(NO)2分子多層膜的原子結(jié)構(gòu)

在三個(gè)穩(wěn)定分子膜吸附結(jié)構(gòu)的基礎(chǔ)上,把分子膜M1和M2分別吸附于的頂位,橋位和兩個(gè)空心位構(gòu)建出分子雙層膜吸附系統(tǒng),分子膜之間的真空層H暫取0.30 nm.對(duì)這些系統(tǒng)進(jìn)行結(jié)構(gòu)優(yōu)化,結(jié)果列于表6,分子膜M2在M2/Rh(111)吸附不穩(wěn)定.

表5 M2/Rh(111)吸附系統(tǒng)的能量及結(jié)構(gòu)參數(shù)Table 5.The energy and structural parameters of M2/Rh(111).

表6 M2/Rh(111)上吸附分子膜M1的能量及結(jié)構(gòu)參數(shù)Table 6.The energy and structural parameters of M1 on M2/Rh(111).

從表6可知,兩個(gè)分子膜之間的相對(duì)平移對(duì)分子間相互作用影響很小.在(NO)2分子膜M2吸附于Rh(111)的兩個(gè)空心位的基礎(chǔ)上吸附分子膜M1,則分子膜M1和M2之間的相互作用約0.02 eV,真空層為(0.31±0.02)nm.而在分子膜M2吸附于Rh(111)的頂位的基礎(chǔ)上吸附分子膜M1,則分子膜M1和M2之間的相互作用小于0.01 eV.因此可判斷要吸附分子多層膜,第一層NO分子最佳的吸附位是兩個(gè)空心位,且第二層分子膜中N-N鍵與襯底的夾角為70°-90°.

為了進(jìn)一步研究(NO)2分子多層膜在Rh(111)表面的吸附結(jié)構(gòu),在分子膜M1吸附于M2/Rh(111)的fcc空心位的基礎(chǔ)上,再加一分子膜M1,對(duì)該分子多層膜吸附系統(tǒng)進(jìn)行結(jié)構(gòu)優(yōu)化,結(jié)果見表7,對(duì)應(yīng)的吸附結(jié)構(gòu)見圖8.由表7可知分子多層膜間的相互作用約為0.010 eV.文獻(xiàn)[12,15]所述分子多層膜中N-N垂直襯底,本文的計(jì)算結(jié)果是分子多層膜中N-N鍵與襯底夾角取值在70°-90°之間.

圖8 (網(wǎng)刊彩色)(NO)2分子多層膜在Rh(111)表面的吸附結(jié)構(gòu)Fig.8.(color online)The sketch map of(NO)2multilayer on Rh(111)surface.

表7 (M1+M2)/Rh(111)上吸附分子膜M1的能量及結(jié)構(gòu)參數(shù)Table 7.The energy and structural parameters of M1 on(M1+M2)/Rh(111).

4 結(jié) 論

利用第一性原理研究了虛擬Rh(111)表面的(NO)2分子鏈、分子單層膜,Rh(111)表面上的(NO)2分子單層膜和多層膜的原子結(jié)構(gòu).(NO)2分子單體在虛擬Rh(111)表面自組裝成兩個(gè)穩(wěn)定的分子鏈,一分子鏈結(jié)合能為0.309 eV,(NO)2分子平面平行有序排列,氧原子和氮原子都呈現(xiàn)(100)結(jié)構(gòu);另一分子鏈結(jié)合能為0.266 eV,(NO)2分子平面旋轉(zhuǎn)后平行有序排列,氧原子和氮原子都呈現(xiàn)(111)結(jié)構(gòu).以這兩個(gè)分子鏈為基礎(chǔ),在覆蓋度為1.00 ML的虛擬上,(NO)2分子自組裝成兩個(gè)穩(wěn)定的分子單層膜(M1和M2),一分子膜M1結(jié)合能為0.318 eV,N-N鍵與襯底的夾角為70°-90°;另一分子膜M2結(jié)合能為0.266 eV,N-N鍵平行襯底,N-O鍵與襯底夾角97°.

在M2/Rh(111)吸附系統(tǒng),(NO)2分子可吸附于頂位、fcc空心位和hcp空心位,通過電荷轉(zhuǎn)移可解釋兩個(gè)空心位的穩(wěn)定性強(qiáng)于頂位.Rh(111)表面(NO)2分子多層膜系統(tǒng)中,(NO)2分子垂直吸附于兩個(gè)空心位,第一層分子膜中N-N鍵平行襯底,第二層及其他層分子膜中N-N鍵與襯底夾角為70°-90°,分子膜間相互作用為0.01 eV,分子層間真空層為(0.31±0.02)nm.

[1]Whitesides G M,Mathias J P,Seto C T 1991Science254 1312

[2]Hickman J J,Ofer D,Laibinis P E,Whitesides G M,Wrighton M S 1991Science252 688

[3]Fujita M,Ibukuro F,Hagihara H,Ogura K 1994Nature367 720

[4]Wang W,Huang L,Zhang Y,Li C M,Zhang H Q,Gu N,Peng L,Zhao L X,Shen H Y,Chen T S,Hao L P 2002Acta Phys.Sin.51 63(in Chinese)[王偉,黃嵐,張宇,李昌敏,張海黔,顧寧,彭力,趙麗新,沈浩瀛,陳堂生,郝麗萍2002物理學(xué)報(bào)51 63]

[5]Hu H L,Zhang K,Wang Z X,Kong T,Hu Y,Wang X P 2007Acta Phys.Sin.56 1674(in Chinese)[胡海龍,張琨,王振興,孔濤,胡穎,王曉平2007物理學(xué)報(bào)56 1674]

[6]Palmer R M J,Ferrige A G,Moncada 1987Nature327 524

[7]Orville-Thomas W J 1954J.Chem.Phys.22 1267

[8]Root T W,Fisher G B,Schmidt L D 1986J.Chem.Phys.85 4679

[9]Loffreda D,Simon D,Sautet P 1998Chem.Phys.Lett.291 15

[10]Wallace W T,Cai Y,Chen M S,Goodman D W 2006J.Phys.Chem.B110 6245

[11]Nakamura I,Kobayashi Y,Hamada H,Fujitani T 2006Surf.Sci.600 3235

[12]Nakai I,Kondoh H,Shimada T,Yokota R,Katayama T,Ohta T 2007J.Chem.Phys.127 024701

[13]Jansen A P J,Popa C 2008Phys.Rev.B78 085404

[14]Wu T Q,Zhu P,Jiao Z W 2012Appl.Surf.Sci.263 502

[15]Brown W A,Gardner P,King D A 1995J.Phys.Chem.99 7065

[16]Perdew J P,Burke K,Ernzerhof M 1996Phys.Rev.Lett.78 3865

[17]Wu T Q,Wang X Y,Jiao Z W,Luo H L,Zhu P 2013Acta Phys.Sin.62 186301(in Chinese)[吳太權(quán),王新燕,焦志偉,羅宏雷,朱萍2013物理學(xué)報(bào)62 186301]

[18]Wu T Q,Wang X Y,Jiao Z W,Luo H L,Zhu P 2014Vacuum101 399

[19]Wu T Q,Wang X Y,Zhou H,Luo H L,Jiao Z W,Zhu P 2014Appl.Surf.Sci.290 425

[20]Wu T Q,Cao D,Wang X Y,Jiao Z W,Jiang Z T,Chen M G,Luo H L,Zhu P 2015Appl.Surf.Sci.339 1

[21]Wu T Q,Cao D,Wang X Y,Jiao Z W,Chen M G,Luo H L,Zhu P 2015Appl.Surf.Sci.330 158

[22]Guo Z H,Yan X H,Xiao Y 2010Phys.Lett.A374 1534

[23]Florence A J,Bardin J,Johnston B,Shankland N,Griffi n T A N,Shankland K 2009Z.Kristallogr.Suppl.30 215

[24]Segall M D,Lindan P J D,Probert M J,Pickard C J,Hasnip P J,Clark S J,Payne M C 2002J.Phys.Condens.Matter14 2717

PACS:63.20.dk,81.16.Dn DOI:10.7498/aps.66.026301

Structure of NO dimer multilayer on Rh(111)?

Wang Chen-ChaoWu Tai-Quan?Wang Xin-Yan Jiang Ying
(Department of Applied Physics,China Jiliang University,Hangzhou 310018,China)

1 September 2016;revised manuscript

3 November 2016)

Molecular self-assembly is the spontaneous organization of molecules under thermodynamic equilibrium conditions into well-defined arrangements via cooperative effects between chemical bonds and weak noncovalent interactions.Molecules undergo self-association without external instruction to form hierarchical structures.Molecular self-assembly is ubiquitous in nature and has recently emerged as a new strategy in chemical biosynthesis,polymer science and engineering.NO monomer is apt to be absorbed on the surfaces of some metals such as Ir(111),Ni(111),Pd(111),Pt(111),Rh(111)and Au(111),and the interactions of NO monomer with the metal surfaces have been extensively studied.When NO monomer is weakly adsorbed on the noble-metal surface,it cannot be reduced completely but forms a stable structure,which is named NO dimer.The first-principle technique is employed to determine the structures of NO dimer((NO)2)molecular chains and monolayers on virtual Rh(111),as well as(NO)2monolayer and multilayer on Rh(111).First,(NO)2monomers are assembled into two stable molecular chains on the virtual Rh(111)surface,whose bind energies are 0.309 and 0.266 eV,respectively.The molecular chains are self-assembly systems,in which(NO)2monomers are parallel and ordered,and the O atoms and N atoms are shown to be of(100)and(111)structures,respectively.Then,the two molecular chains are assembled into two stable monolayers(denoted as M1 and M2)on the virtualand the coverage is 1.00 ML.In the M1 monolayer,the angle between the N—N bond of(NO)2monomer and the substrate is in a range of 70?-90?,and in the M2 monolayer,the N—N bond is parallel to the substrate.

In the adsorption system of M2/Rh(111),(NO)2molecules can be adsorbed on the top as well as the hcp and fcc hollow sites.When(NO)2molecules are adsorbed on the top site,the adsorption system is best described by the electron structure Rh+0.14—N0=O-0.14,and when(NO)2molecules are absorbed on the two hollow sites,the adsorption system is described by the electron structure Rh+0.34—N-0.18=O-0.16.Therefore,(NO)2molecules are more apt to be adsorbed on the two hollow sites than on the top site.In the adsorption systems of M1+M2/Rh(111)and M1+(M1+M2)/Rh(111),(NO)2molecules are adsorbed vertically on the two hollow sites,the N—N bond is parallel to the substrate in the first monolayer,and the angle between the N—N bond and the substrate is in a range of 70?-90?in the second and third monolayers.The interaction between the neighbor monolayers is about 0.01 eV,and the thickness of the vacuum layer is 0.31 nm±0.02 nm.

NO dimer,monolayer,multilayer,self-assembly

:63.20.dk,81.16.Dn

10.7498/aps.66.026301

?浙江省自然科學(xué)基金(批準(zhǔn)號(hào)：LY13E080007)資助的課題.

?通信作者.E-mail：buckyballling@hotmail.com

*Project Supported by the National Natural Science Foundation of Zhejiang Province,China(Grant No.LY13E080007).

?Corresponding author.E-mail：buckyballling@hotmail.com