密度泛函理論研究鳥氨酸鹽在Mg(0001)表面的吸附機(jī)理

王 群, 張繼巍, 周乃武, 趙永馳

(1. 綿陽師范學(xué)院生命科學(xué)與技術(shù)學(xué)院, 綿陽 621006; 2. 綿陽師范學(xué)院歷史文化與旅游管理學(xué)院,綿陽 621006; 3. 綿陽師范學(xué)院數(shù)學(xué)與計(jì)算機(jī)科學(xué)學(xué)院, 綿陽 621006)

密度泛函理論研究鳥氨酸鹽在Mg(0001)表面的吸附機(jī)理

王 群1, 張繼巍1, 周乃武2, 趙永馳3

(1. 綿陽師范學(xué)院生命科學(xué)與技術(shù)學(xué)院, 綿陽 621006; 2. 綿陽師范學(xué)院歷史文化與旅游管理學(xué)院,綿陽 621006; 3. 綿陽師范學(xué)院數(shù)學(xué)與計(jì)算機(jī)科學(xué)學(xué)院, 綿陽 621006)

鎂及鎂合金具有與人體自然骨的密度相近, 良好的生物相容性等特點(diǎn), 但作為生物醫(yī)用材料植入人體之后, 及易降解和腐蝕. 本論文用密度泛函理論研究了鳥氨酸鹽在Mg(0001)表面的吸附, 通過吸附能, 態(tài)密度, 電荷差分密度等分析發(fā)現(xiàn)鳥氨酸鹽在Mg(0001)表面存在強(qiáng)烈的Mg-O和Mg-N相互作用, 從而鳥氨酸鹽有望在鎂金屬的表面可以形成一層致密的氧化膜, 達(dá)到改善鎂金屬材料的耐腐蝕性, 降低降解速度的目的.

鎂及鎂合金; 鳥氨酸鹽; 耐腐蝕性; 密度泛函理論; 氧化膜

1 引 言

近年來鎂及鎂合金生物醫(yī)用材料的研究開發(fā)因其明顯的優(yōu)點(diǎn)受到了人們的廣泛關(guān)注[1], 鎂和鎂合金的密度與人體的骨皮質(zhì)的密度相似, 都大約在1.75 g/cm3左右[2]; 與人體具有良好的生物相容性, 骨骼組織中儲(chǔ)存著大量的生物鎂, 同時(shí)鎂也是組成骨的主要元素之一, 能夠促進(jìn)鈣的沉積, 引導(dǎo)新骨組織的生長(zhǎng), 如果缺乏鎂會(huì)阻止骨組織生長(zhǎng), 降低造骨細(xì)胞和破骨細(xì)胞的活性, 并導(dǎo)致骨質(zhì)脆弱[3].

但將鎂及其合金作為生物醫(yī)用材料最大的缺陷就是降解速度太快, 易腐蝕, 尤其是在生理環(huán)境中, 含有大量的氯離子, 氧化鎂表面膜會(huì)遭到一定程度的損壞, 導(dǎo)致腐蝕速度變得更加快速[4,5], 因此, 近年來, 研究者們對(duì)改善鎂合金的耐腐蝕性進(jìn)行了大量的研究, 比如用離子注入的方法對(duì)鎂合金表面進(jìn)行改性[6-9], 因其在材料表面內(nèi)側(cè)形成一層新的表面合金層, 從而改變表面狀態(tài)[10], 因此它們的耐腐蝕性得到了很大改善, 但離子注入所得到的改性層很薄, 研究者們也嘗試用別的方法進(jìn)行改性, 比如在鎂合金的溶液中加入有機(jī)添加劑, 可以在其表面生成一層致密的氧化膜, 使其耐腐蝕性大大增加[11-16].

鎂合金氧化膜的制備過程中, 實(shí)驗(yàn)已經(jīng)證實(shí)可以在電解質(zhì)溶液中添加各種有機(jī)添加劑改善膜的致密性、孔洞分布的均勻性等, 但具體的作用機(jī)理還不是相當(dāng)明確, 計(jì)算機(jī)模擬的方法為我們提供了一條好的路徑, 目前用計(jì)算機(jī)模擬手段進(jìn)行機(jī)理研究的報(bào)道已經(jīng)很多, 并取得了一些有價(jià)值的研究成果[17-19]. 我們?cè)O(shè)想是否可以在堿性電解液溶液中添加一種有機(jī)添加劑鳥氨酸鹽, 是否可以改變這種性質(zhì)呢? 鳥氨酸是一種堿性氨基酸, 存在于生物體多種組織和細(xì)胞中, 具有生物活性, 在生物代謝中起重要作用, 如果可行, 保護(hù)膜為鳥氨酸的鎂合金生物醫(yī)用材料就具有能夠改善其耐腐蝕性, 降低降解率, 在人體內(nèi)無毒無害的優(yōu)點(diǎn). 現(xiàn)目前為止, 實(shí)驗(yàn)研究還未見報(bào)道, 本論文用DFT詳細(xì)研究了鳥氨酸鹽在Mg(0001)表面的吸附機(jī)制, 發(fā)現(xiàn)它們之間有很強(qiáng)的相互作用, 說明鳥氨酸鹽有望在鎂合金表面形成一層致密的氧化膜, 本理論研究可以為實(shí)驗(yàn)研究者提供有價(jià)值的理論參考.

2 模型的建立和計(jì)算方法

2.1 模型的建立

鳥氨酸陰離子(Orn anion)模型: 由于在堿性溶液中鳥氨酸鹽與鎂起反應(yīng)的主要是鳥氨酸陰離子, 因此我們構(gòu)建其陰離子模型, 其化學(xué)式是-OOCCH(NH2)CH2CH2CH2NH2, 有一個(gè)α-NH2, 一個(gè)δ-NH2和α-COO-, 如圖1a.

Mg(0001)模型: 由于Mg晶體的密排面為(0001)面, 此面的能量相對(duì)最低, 因此在構(gòu)筑晶胞時(shí), 我們構(gòu)建了一個(gè)包含72個(gè)鎂原子的Mg(0001)晶胞, 表面大小為9.628 ?×16.676 ?. 賦予Mg周期性, 并在其上方加上厚度為25 ?的真空層. 如圖1b, 因?yàn)橹唤⒘艘粋€(gè)Mg的表面, 為了簡(jiǎn)化, 在以后的模型建立中我們就直接用Mg代替Mg(0001)表面.

圖1 (a) 鳥氨酸陰離子的三維結(jié)構(gòu); (b) Mg(0001)的三維結(jié)構(gòu) (碳原子：灰色; 氧化子：紅色; 氮原子：藍(lán)色; 氫原子：白色; 鎂原子：綠色)Fig. 1 The three-dimensional structure models of Orn anion (a) and Mg(0001) (b) . Color Codes: Carbon, gray; Oxygen, red; Nitrogen, blue; Hydrogen, white; Magnesium, green

Mg-Orn-O-standing: 結(jié)構(gòu)優(yōu)化之前, Orn中α-COO-站在Mg(0001)表面, 如圖2a, 在結(jié)構(gòu)優(yōu)化之后, 如圖3a.

Mg-Orn-N-O-standing: 結(jié)構(gòu)優(yōu)化之前, Orn中α-COO-和α-NH2站Mg(0001)表面, 如圖2b, 在結(jié)構(gòu)優(yōu)化之后, 如圖3b.

Mg-Orn-N-lying: 結(jié)構(gòu)優(yōu)化之前, Orn中α-NH2和δ-NH2平躺在Mg(0001)表面, 如圖2c, 在結(jié)構(gòu)優(yōu)化之后, 如圖3c.

Mg-Orn-N-standing: 結(jié)構(gòu)優(yōu)化之前, Orn中δ-NH2站在Mg(0001)表面, 如圖2d, 在結(jié)構(gòu)優(yōu)化之后, 如圖3d.

圖2 結(jié)構(gòu)優(yōu)化前: (a) Mg-Orn-O-standing模型; (b) Mg-Orn-N-O-standing模型; (c) Mg-Orn-N-lying模型; (d) Mg-Orn-N-standing模型Fig. 2 Before the geometry optimization: (a) Mg-Orn-O-standing model; (b) Mg-Orn-N-O-standing model; (c) Mg-Orn-N-lying model; (d) Mg-Orn-N-standing model

圖3 結(jié)構(gòu)優(yōu)化后: (a) Mg-Orn-O-standing模型; (b) Mg-Orn-N-O- standing模型; (c) Mg-Orn-N-lying 模型; (d) Mg-Orn-N-standing 模型Fig. 3 After the geometry optimization: (a) Mg-Orn-O- standing model; (b) Mg-Orn-N-O-standing model; (c) Mg-Orn-N-lying model; (d) Mg-Orn-N-standing model

2.2 計(jì)算方法

Materials Studio(Accelrys, San Diego, CA)中的基于第一性原理DFT的Dmol3進(jìn)行計(jì)算研究. 該模塊中用數(shù)值基組表示物理波函數(shù)[20,21]. 采用等同于Gaussian中6-31G**基組的DNP雙重?cái)?shù)值基組(DNP double numerical basis set), 用DFT半核贗勢(shì)( DFT semicore pseudopotentials)處理內(nèi)層電子[22], 用Perdew-Burke-Ernzerhof (PBE) 廣義梯度近似(GGA)計(jì)算交換關(guān)聯(lián)能[23], 根據(jù)Monkhorst-Pack方案, 布里淵區(qū)域(Brillouin zone)內(nèi)的K點(diǎn)網(wǎng)格設(shè)置為2×2×1[24], 費(fèi)米拖尾效應(yīng)(Fermi smearing)值設(shè)置為0.005 Ha(1Ha = 27.2114 eV), 全局軌道截?cái)?global orbital cutoff)值設(shè)置為4.8 ?. 體系的幾何構(gòu)型優(yōu)化和能量計(jì)算收斂標(biāo)準(zhǔn)為：(a)自洽循環(huán)數(shù)量級(jí)達(dá)到1.0×10-6Ha/atom及以下;(b)能量值數(shù)量級(jí)達(dá)到1.0×10-5Ha /atom及以下;(c)最大應(yīng)力≤0.002 Ha /?;(d)最大位移≤0.005 ?.

2.3 吸附能

吸附能(adsorption energy,Eads)是反映Orn和Mg(0001)表面之間相互作用強(qiáng)度的參數(shù), 可以根據(jù)方程(1)計(jì)算得到

Eads=EOrn+Mg-(EOrn+EMg)

(1)

其中EOrn+Mg,EOrn,EMg分別代表吸附體系的總能Etotal, Orn的能量和Mg(0001)表面的能量. 如果吸附能Eads值為負(fù)值, Orn就可以自發(fā)的吸附在Mg(0001)表面, 而且負(fù)值越大, 吸附越穩(wěn)定. 表1給出了所有吸附體系的吸附能.

表1 鳥氨酸在Mg(0001)面的吸附能

3 結(jié)果和討論

3.1 吸附能的分析

Mg-Orn-O-standing: 結(jié)構(gòu)優(yōu)化之前, Orn中羧基的兩個(gè)氧原子站在Mg(0001)表面, 如圖2a, 結(jié)構(gòu)優(yōu)化之后, 羧基中的兩個(gè)氧原子仍然站在Mg(0001)表面(3a). Orn與Mg(0001)之間的吸附能很大(-4.72 eV), 這表明這兩者之間有很強(qiáng)的作用力. 而且羧基中氧原子和Mg之間的距離分別是2.118 ?, 2.953 ?, 2.178 ?, 這進(jìn)一步說明Orn中α-COO-氧原子和Mg形成強(qiáng)烈的Mg-O相互作用.

Mg-Orn-N-O-standing: 結(jié)構(gòu)優(yōu)化之前, Orn中α-COO-和α-NH2站Mg(0001)表面, 如圖2a, 結(jié)構(gòu)優(yōu)化之后, 結(jié)構(gòu)沒有太大的變化, 只是Mg把羧基拉的更靠近其表面(圖3b). 這樣由于形成了Mg-N(2.358 ?), Mg-O(2.198 ?, 2.413 ?)之間的強(qiáng)烈的相互作用, 它們之間的吸附能變得更大了(-5.04 eV).

Mg-Orn-N-lying: 結(jié)構(gòu)優(yōu)化之前, Orn中α-NH2和δ-NH2平躺在Mg(0001)表面(圖2c), 結(jié)構(gòu)優(yōu)化之后, 結(jié)構(gòu)變化很大, 如圖3c, 此時(shí)結(jié)構(gòu)幾乎發(fā)生了翻轉(zhuǎn), α-COO-和δ-NH2平躺在金屬的表面, 它們之間形成了很強(qiáng)的Mg-N(2.323 ?), Mg-O(2.151 ?, 2.206 ?, 2.966 ?), 由于相互作用的鍵長(zhǎng)比第二個(gè)模型的鍵長(zhǎng)更短, 因此吸附能進(jìn)一步增大, 變?yōu)?5.32 eV.

Mg-Orn-N-standing: 結(jié)構(gòu)優(yōu)化之前, Orn中δ-NH2站在Mg(0001)表面(圖2d), 結(jié)構(gòu)優(yōu)化之后, 和結(jié)構(gòu)3c類似, 金屬表面的Mg原子把α-COO-也拉到了表面與它發(fā)生強(qiáng)烈的相互作用, 如圖3d. 此時(shí)的吸附能比結(jié)構(gòu)3c稍大(-5.35 eV), 這是因?yàn)樗鼈冎gMg-O之間的距離比前三個(gè)模型都短(2.063 ?, 2.059 ?), 而且它們之間也形成了強(qiáng)烈的Mg-N(2.318 ?)作用力.

總結(jié)：從這四個(gè)優(yōu)化之后的模型我們可以看到, 當(dāng)鳥氨酸中氨基更接近金屬表面時(shí), 金屬表面的Mg原子更傾向于把鳥氨酸中α-COO-拉的更靠近其表面, 這也說明了α-COO-中的氧原子更易和金屬表面的Mg原子形成很強(qiáng)的Mg-O作用力, 但同時(shí)鳥氨酸中氨基和金屬表面的Mg原子形成Mg-N作用力時(shí), 會(huì)使鳥氨酸和Mg之間的相互作用力變得更強(qiáng), 從這點(diǎn)來看, Mg-N之間的相互作用力也不容忽視.

3.2 布局電荷分析

布局電荷分析有助于幫助我們了解鳥氨酸陰離子在吸附到Mg(0001)面之后電子的轉(zhuǎn)移情況. 電子之間的轉(zhuǎn)移情況如表2, 從表中我們可以看到, 與鳥氨酸相互作用的Mg原子失去了電子, 鳥氨酸中O和N原子得到了電子, 也就是說Mg原子上電子轉(zhuǎn)移到了N和O原子上, 兩者之間形成了Mg-N, Mg-O相互作用.

表2 原子的布局電荷

3.3 差分電荷密度分析

電子密度圖也可以直觀地展示生物分子和生物材料之間的相互作用, 電子密度差揭示了在吸附過程中電子密度的改變, 如下圖公式(2)所示.在電子密度等值面圖中紅色代表該區(qū)域原子得到電子, 藍(lán)色代表該區(qū)域原子失去電子.

Δρ=ρOrn+Mg-(ρOrn+ρMg)

(2)

鳥氨酸在Mg(0001)面上的吸附, 差分電荷密度圖如4a-d所示, 從圖中可以看到, 鳥氨酸的羧基和氨基都是紅色, 表示他們?cè)诤蚆g(0001)表面相互作用時(shí)得到了電子, 電子從Mg(0001)表面轉(zhuǎn)移到了鳥氨酸中氨基的N原子和羧基的氧原子上, 這個(gè)結(jié)果與電荷布局分析的結(jié)果也是一致的.而且, 我們可以發(fā)現(xiàn)在鳥氨酸羧基氧原子上紅色區(qū)域更大, 說明比羧基中的O原子更易得到電子, 這也可以證實(shí)當(dāng)鳥氨酸與Mg相互作用時(shí),Mg-O之間的相互作用是主要的作用.

3.4 態(tài)密度的分析

基于能帶理論的態(tài)密度(densityofstates,DOS)分析可以幫助我們進(jìn)一步揭示Orn與Mg之間相互作用的本質(zhì).如果相鄰原子的局域態(tài)密度在同一個(gè)能量上同時(shí)出現(xiàn)了尖峰, 則我們將其稱之為雜化峰(hybridizedpeak), 它可以直觀地向我們展示相鄰原子之間的作用強(qiáng)弱[25]. 從圖5a-5d中的DOS圖譜可以看出Mg與Orn之間都可以形成相互作用力. 四個(gè)模型中模型三和模型四的相互作用力最強(qiáng), 那么我們仔細(xì)分析這兩個(gè)模型相互作用的相鄰原子的PDOS圖, 如圖模型三的PDOSMg24-N80在-10eV, -8.42eV, -7.26eV, -6.12eV, -4.92eV位置有峰重疊 (圖6a),Mg40-O76在-8.93eV, -7.18eV, -6.13eV位置有峰重疊 (圖7a),Mg48-O81在-10.78eV, -10.07eV, -8.09eV, -7.23eV, -6.14eV, -3.19eV, -1.40eV, 0.32eV位置有峰重疊(圖8a),說明上面的兩兩原子之間在這些位置有相互作用的可能性.另外, 我們?cè)俜治瞿Ｐ退牡腜DOS圖, 從圖中我們?nèi)匀豢梢钥吹?Mg46-N80在-9.79eV, -8.39eV, -6.78eV, -4.64eV位置有峰重疊 (圖6b),Mg62-O76在-8.65eV,-7.34eV, -6.82eV, -5.43eV, 0.03eV位置有峰重疊 (圖7b),Mg70-O81在-9.25eV, -8.60eV, -7.33eV, -2.27eV, -1.57eV, 0.08eV位置有峰重疊 (圖8b),說明上面的兩兩原子之間在這些位置有相互作用的可能性, 這些PDOS結(jié)果的分析也和上述的分析結(jié)果一致.

圖4 差分電子密度圖. (a) Mg-Orn-O-standing模型; (b) Mg-Orn-N-O-standing模型; (c) Mg-Orn-N-lying模型; (d) Mg-Orn-N-standing模型. 紅色代表該區(qū)域原子得到電子, 藍(lán)色代表該區(qū)域原子失去電子Fig. 4 A three-dimensional isosurface plot of the deformation electron density with an isovalue of 0.05 e/? 3. (a) Mg-Orn-O-standing model; (b) Mg-Orn- N-O-standing model; (c) Mg-Orn-N-lying model; (d) Mg-Orn-N-standing model. The charge accumulation and charge depletion are represented by red and blue, respectively

圖5 (a) Orn, Mg, and Mg-Orn-O-standing的態(tài)密度圖; (b) Orn, Mg, and Mg-Orn- N-O-standing的態(tài)密度圖; (c) Orn, Mg, and Mg-Orn-N-lying的態(tài)密度圖; (d) Orn, Mg, and Mg-Orn- N-standing的態(tài)密度圖Fig. 5 DOSs for the (a) Orn, Mg, and Mg-Orn-O-standing; (b) Orn, Mg, and Mg-Orn- N-O-standing; (c) Orn, Mg, and Mg-Orn-N-lying; (d) Orn, Mg, and Mg-Orn-N- standing

圖6 (a) Mg-Orn-N-lying模型中Mg24,N80的分態(tài)密度圖; (b) Mg-Orn-N-standing中模型Mg46,N80的分態(tài)密度圖Fig. 6 The PDOSs for the (a) Mg24, N80 atoms in the Mg-Orn-N-lying model; (b) Mg46, N80 atoms in the Mg-Orn-N-standing model

圖7 (a) Mg-Orn-N-lying模型中Mg40,O76的分態(tài)密度圖; (b) Mg-Orn-N-standing中模型Mg62,O76的分態(tài)密度圖Fig. 7 The PDOSs for the (a) Mg40, O76 atoms in the Mg-Orn-N-lying model; (b) Mg62, O76 atoms in the Mg-Orn-N-standing model

圖8 (a) Mg-Orn-N-lying模型中Mg48,O81的分態(tài)密度圖; (b) Mg-Orn-N-standing中模型Mg70,O81的分態(tài)密度圖Fig. 8 The PDOSs for the (a) Mg48, O81 atoms in the Mg-Orn-N-lying model; (b) Mg70, O81 atoms in the Mg-Orn-N-standing model

4 結(jié) 論

本論文首次用DFT的方法研究了有機(jī)添加劑鳥氨酸鹽在堿性溶液中與Mg(0001)表面相互作用的機(jī)理, 發(fā)現(xiàn)鳥氨酸鹽可以與Mg(0001)表面發(fā)生強(qiáng)烈的相互作用, 主要是存在Mg-O和Mg-N相互作用; 且Mg-O強(qiáng)于Mg-N相互作用, 是主要的相互作用. 由于這些強(qiáng)烈的靜電相互作用導(dǎo)致了鳥氨酸鹽可以在Mg及其合金表面生成一層致密的有機(jī)氧化膜, 從而達(dá)到改善腐蝕性和降低降解速率的目的.

[1]LiN,ZhengY.Novelmagnesiumalloysdevelopedforbiomedicalapplication:Areview[J].J.Mater.Sci.Technol., 2013, 29 (6): 489.

[2] Rittel D, Wang Z G. Thermo-mechanical aspects of adiabatic shear failure of AM50 and Ti6Al4V alloys [J].Mech.Mater., 2008, 40 (8): 629.

[3] Zhang J, Zong Y, Fu P H,etal. Application and development prospect of magnesium alloys as biomedical materials [J].CJTER., 2009, 13 (29): 5747 (in Chinese) [張佳, 宗陽, 付彭懷, 等. 鎂合金在生物醫(yī)用材料領(lǐng)域的應(yīng)用及發(fā)展前景 [J]. 中國(guó)組織工程研究與臨床康復(fù), 2009, 13 (29): 5747]

[4] Gao J C, Hu D, Song C J. Medical magnesium alloy degradation and its impact on human physiology [J].J.Funct.Mater., 2012, 43 (19): 2577 (in Chinese) [高家誠(chéng), 胡德, 宋長(zhǎng)江. 醫(yī)用鎂合金降解及其對(duì)人體的影響 [J]. 功能材料, 2012, 43 (19): 2577]

[5] Song G L, Song S Z. Corrosion behavior of pure magnesium in a simulated body fluid [J].ActaPhys.Chim.Sin., 2006, 22 (10): 1222 (in Chinese) [宋光鈴, 宋詩哲. 鎂在人體模擬液中的腐蝕行為 [J]. 物理化學(xué)學(xué)報(bào), 2006, 22 (10): 1222]

[6] Lei M K, Li P, Yang H G,etal. Wear and corrosion resistance of Al ion implanted AZ31 magnesium alloy [J].Surf.Coat.Technol., 2007, 201 (9): 5182.

[7] Liu C, Xin Y, Tian X,etal. Corrosion resistance of titanium ion implanted AZ91 magnesium alloy [J].J.Vac.Sci.Technol. A, 2007, 25 (2): 334.

[8] Wang X, Zeng X, Wu G,etal. The effect of Y-ion implantation on the oxidation of AZ31 magnesium alloy [J].Mater.Lett., 2007, 61 (4): 968.

[9] Wang X, Zeng X, Wu G,etal. Effects of tantalum ion implantation on the corrosion behavior of AZ31 magnesium alloys [J].J.Alloy.Compd., 2007, 437 (1): 87.

[10] Huo H W, Li Y. Corrosion and production of magnesium alloys [J].Mater.Rev., 2001, 15 (7): 25 (in Chinese) [霍宏偉, 李瑛. 鎂合金的腐蝕與防護(hù) [J]. 材料導(dǎo)報(bào), 2001, 15 (7): 25]

[11] Cui X, Li Q, Li Y,etal. Microstructure and corrosion resistance of phytic acid conversion coatings for magnesium alloy [J].Appl.Surf.Sci., 2008, 255 (5): 2098.

[12] Dai Y, Li Q, Gao H,etal. Effects of five additives on electrochemical corrosion behaviours of AZ91D magnesium alloy in sodium chloride solution [J].Surf.Eng., 2011, 27 (7): 536.

[13] Guo X, An M, Yang P,etal. Effects of benzotriazole on anodized film formed on AZ31B magnesium alloy in environmental-friendly electrolyte [J].J.Alloy.Compd., 2009, 482 (1): 487.

[14] Liu Y, Wei Z L, Yang F W,etal. Environmental friendly anodizing of AZ91D magnesium alloy in alkaline borate-benzoate electrolyte [J].J.Alloy.Compd., 2011, 509 (22): 6440.

[15] Liu Y, Yang F W, Wei Z L,etal. Anodizing of AZ91D magnesium alloy using environmental friendly alkaline borate-biphthalate electrolyte [J].Trans.NonferrousMet.Soc.China, 2012, 22 (7): 1778.

[16] Liu Y, Yang F W, Fan J,etal. Effect of organic additives on the properties of anodic coatings on magnesium alloys [J].Corros.Sci.Prot.Technol., 2012, 24 (3): 253 (in Chinese) [劉妍, 楊富巍, 樊娟, 等. 有機(jī)添加劑對(duì)鎂合金陽極氧化膜性能的影響 [J]. 腐蝕科學(xué)與防護(hù)技術(shù), 2012, 24 (3): 253]

[17] Zhang F L, Xu B H, Xu J H,etal. Density functional theory study of HCOOH adsorption on Pd-Fe(1 1 1)( nPd: nFe=1 : 1) surfaces [J].J.At.Mol.Phys., 2014, 31 (1): 161 (in Chinese) [張福蘭, 徐伯華, 徐建華, 等. HCOOH在Pd-Fe(1 1 1)( nPd: nFe=1 : 1) 表面吸附的密度泛函理論研究 [J]. 原子與分子物理學(xué)報(bào), 2014, 31 (1): 161]

[18] Huang W Y, Cheng C, Tao T. Density functional theory study of methyl/hydroxyl adsorption on Ir(1 l 1) surface [J].J.At.Mol.Phys., 2013, 30 (3): 451 (in Chinese) [黃武英, 程春, 陶濤. 甲基和羥基在Ir(1 l 1)表面吸附的密度泛函理論研究 [J]. 原子與分子物理學(xué)報(bào), 2013, 30 (3): 451]

[19] Zhang J P, Zhu X L, Zhang Y Y,etal. Investigations of the atomic structure and electronic states for Ni atatom adsorption on Al(1 1 1)surface using first-principles [J].J.At.Mol.Phys., 2013, 31 (3): 463 (in Chinese) [張金平, 朱曉玲, 張洋洋, 等. 第一性原理計(jì)算Ni原子吸附在A1(11 1)表面原子結(jié)構(gòu)和電子態(tài) [J]. 原子與分子物理學(xué)報(bào), 2014, 31 (3): 463]

[20] Delley B. All electron numerical method for solving the local density functional for polyatomic molecules [J].J.Chem.Phys., 1990, 92 (1): 508.

[21] Delley B. From molecules to solids with the DMol3 approach [J].J.Chem.Phys., 2000, 113 (18): 7756.

[22] Delley B. Hardness conserving semilocal pseudopotentials [J].Phys.Rev. B, 2002, 66 (15): 155125.

[23] Perdew J P, Burke K, Ernzerhof M. Generalized gradient approximation made simple [J].Phys.Rev.Lett., 1996, 77 (18): 3865.

[24] Monkhorst H J, Pack J D. Special points for Brillouin-zone integrations [J].Phys.Rev. B, 1976, 13 (12): 5188.

[25] Kusama H, Orita H, Sugihara H. DFT investigation of the TiO2band shift by nitrogen-containing heterocycle adsorption and implications on dye-sensitized solar cell performance [J].Sol.Energ.Mater.Sol.Cells, 2008, 92 (1): 84.

Density functional theory study of the adsorption mechanism of ornithine salt on the Mg(0001) surface

WANG Qun1, ZHANG Ji-Wei1, ZHOU Nai-Wu2, ZHAO Yong-Chi3

(1. Collage of Life Science and Biotechnology, MianYang Normal University, Mianyang 621006, China;2. Collage of Historical Culture and Tourism Management, MianYang Normal University, Mianyang 621006, China;3. Collage of Mathematics and Compute Science, MianYang Normal University, Mianyang 621006, China)

Magnesium (Mg) and its alloys have the characteristics of similar density to human natural bones and the excellent biocompatibility. However, it is easy to degrade and corrode after implanting into the human body as biomedical materials. In this paper, we investigated the adsorption of ornithine (Orn) salt on the Mg(0001) surface using the density functional theory (DFT). It is found that there are strong Mg-O and Mg-N interactions between Orn and Mg(0001) surfaces through analyzing the adsorption energies, density of states (DOS), charge difference density etc. Therefore, Orn salt is expected to form a layer dense oxidation film on the Mg and its alloys surfaces, and then to reach the aims of the larger corrosion resistance and the lower degradation speed.

Magnesium and its alloy; Ornithine salt; Corrosion resistance; Density functional theory; Oxidation film

綿陽師范學(xué)院校級(jí)資助項(xiàng)目 (2013B16,2014A06)；四川省教育廳基金資助項(xiàng)目(14ZB0266)；綿陽師范學(xué)院校級(jí)教改基金資助項(xiàng)目(mnu-JY1408)

王群(1981—),女, 湖北荊州人, 研究生, 講師, 從事化學(xué), 計(jì)算化學(xué)及計(jì)算材料學(xué)方面的研究. E-mail: carlcili@163.com

103969/j.issn.1000-0364.2015.08.028

O647.3

A

1000-0364(2015)08-0686-09

投稿日期: 2015-01-14