亚洲欧美日本视频在线观看-国产精品粉嫩一区二区-欧美高清在线不卡视频-日韩乱码在线高清观看

<th id="2egsc"></th>

快捷下單入口關(guān)于合作招聘新人手冊(cè) 會(huì)員中心

熱線：400-152-6858

測(cè)試狗科研服務(wù)

預(yù)存免費(fèi)試測(cè) 登錄

Document

當(dāng)前位置：文庫(kù)百科 ? 文章詳情

ASE原子模擬環(huán)境入門(mén)之操作原子

來(lái)源：時(shí)間：2022-12-16 16:36:30 瀏覽：3024次

銀吸附在鎳板上

我們建立一個(gè)由四個(gè)Ni原子組成的單層平板和一個(gè)Ag吸附原子構(gòu)成的模型。

定義原子板


>>> from math import sqrt
>>> from ase import Atoms
>>> a = 3.55
>>> atoms = Atoms('Ni4',
...               cell=[sqrt(2) * a, sqrt(2) * a, 1.0, 90, 90, 120],
...               pbc=(1, 1, 0),
...               scaled_positions=[(0, 0, 0),
...                                 (0.5, 0, 0),
...                                 (0, 0.5, 0),
...                                 (0.5, 0.5, 0)])
>>> atoms.center(vacuum=5.0, axis=2)

看看晶胞和原子的位置:


>>> atoms.cell
Cell([[5.020458146424487, 0.0, 0.0], [-2.5102290732122423, 4.347844293440141, 0.0], [0.0, 0.0, 10.0]])
>>> atoms.positions
array([[ 0.        ,  0.        ,  5.        ],
        [ 2.51022907,  0.        ,  5.        ],
        [-1.25511454,  2.17392215,  5.        ],
        [ 1.25511454,  2.17392215,  5.        ]])
>>> atoms[0]
Atom('Ni', [0.0, 0.0, 5.0], index=0)

將這個(gè)結(jié)構(gòu)寫(xiě)入一個(gè)文件，并通過(guò)調(diào)用ase.gui來(lái)繪制整個(gè)系統(tǒng):

>>> from ase.visualize import view

>>> atoms.write('slab.xyz')

>>> view(atoms)

在查看器(ase gui)中，可以在所有三個(gè)方向上重復(fù)單元格(使用Repeat ? View window)。從命令行，使用ase gui -r 3,3,2 slab.xyz。

我們現(xiàn)在在一個(gè)高度為h=1.9?的三原子中心位點(diǎn)添加一個(gè)吸附原子 :


>>>h = 1.9
>>> relative = (1 / 6, 1 / 6, 0.5)
>>> absolute = np.dot(relative, atoms.cell) + (0, 0, h)
>>> atoms.append('Ag')
>>> atoms.positions[-1] = absolute

現(xiàn)在的結(jié)構(gòu)是這樣的:


view(atoms)

構(gòu)造界面

現(xiàn)在，我們用Ni(111)和水做一個(gè)界面。首先我們需要一層水。在這個(gè)腳本W(wǎng)L.py中構(gòu)建了一層水，并保存在文件WL.traj中。

WL.py


import numpy as np
from ase import Atoms
p = np.array(
    [[0.27802511, -0.07732213, 13.46649107],
     [0.91833251, -1.02565868, 13.41456626],
     [0.91865997, 0.87076761, 13.41228287],
     [1.85572027, 2.37336781, 13.56440907],
     [3.13987926, 2.3633134, 13.4327577],
     [1.77566079, 2.37150862, 14.66528237],
     [4.52240322, 2.35264513, 13.37435864],
     [5.16892729, 1.40357034, 13.42661052],
     [5.15567324, 3.30068395, 13.4305779],
     [6.10183518, -0.0738656, 13.27945071],
     [7.3856151, -0.07438536, 13.40814585],
     [6.01881192, -0.08627583, 12.1789428]])
c = np.array([[8.490373, 0., 0.],
              [0., 4.901919, 0.],
              [0., 0., 26.93236]])
W = Atoms('4(OH2)', positions=p, cell=c, pbc=[1, 1, 0])
W.write('WL.traj')

現(xiàn)在運(yùn)行WL.py腳本，從traj文件中讀取原子對(duì)象:


from ase.io import read
>>> W = read('WL.traj')

使用視圖的看看結(jié)構(gòu)

我們看看這個(gè)單位晶胞


W.cell
Cell([8.490373, 4.901919, 26.93236])

我們需要一個(gè)鎳(111)平板，它與水盡可能接近。一個(gè)2 × 4正交的fcc111 超晶胞應(yīng)該足夠好了。


>>> from ase.build import fcc111
>>> slab = fcc111('Ni', size=[2, 4, 3], a=3.55, orthogonal=True)


slab.cell
Cell([5.020458146424487, 8.695688586880282, 0.0])

看看這兩個(gè)單元晶胞，我們可以看到，如果我們把其中一個(gè)單元晶胞在平面上旋轉(zhuǎn)90度，它們的匹配度相差大約2%。讓我們旋轉(zhuǎn)單元晶胞:

W.cell = [W.cell[1, 1], W.cell[0, 0], 0.0]

用rotate()旋轉(zhuǎn)分子：

>>> W.rotate(90, 'z', center=(0, 0, 0))

再將原子放回到晶胞里

>>> W.wrap()

wrap()方法僅在啟用周期性邊界條件時(shí)有效。我們?cè)贜i(111)和水之間有2%的晶格失配，所以我們?cè)谄矫嫔峡s放水的晶胞以匹配平板的晶胞。參數(shù)scale_atoms=True表示原子位置應(yīng)該隨單元格縮放。默認(rèn)值是scale_atoms=False，表示單元格更改時(shí)笛卡爾坐標(biāo)保持不變。


>>> W.set_cell(slab.cell, scale_atoms=True)
>>> zmin = W.positions[:, 2].min()
>>> zmax = slab.positions[:, 2].max()
>>> W.positions += (0, 0, zmax - zmin + 1.5)

最后我們將水復(fù)制到平板上:

>>> interface = slab + W

>>> interface.center(vacuum=6, axis=2)

>>> interface.write('NiH2O.traj')

兩個(gè)atoms對(duì)象的求和，會(huì)將兩個(gè)atoms對(duì)象的原子位置都包括進(jìn)來(lái)，而使用第一個(gè)atoms對(duì)象的超晶胞。

計(jì)算狗·模擬計(jì)算

上一篇：ASE原子環(huán)境模擬之調(diào)用外部計(jì)算軟件的方法

下一篇：ASE原子模擬環(huán)境入門(mén)之結(jié)構(gòu)優(yōu)化: H2O

相關(guān)文章

基礎(chǔ)理論丨一文了解XPS（概念、定性定量分析、分析方法、譜線結(jié)構(gòu)）

2020-05-03

晶體結(jié)構(gòu)可視化軟件 VESTA使用教程（下篇）

2021-01-22

手把手教你用ChemDraw 畫(huà)化學(xué)結(jié)構(gòu)式：基礎(chǔ)篇

2021-06-19

【科研干貨】電化學(xué)表征：循環(huán)伏安法詳解（上）

2019-10-25

【科研干貨】電化學(xué)表征：循環(huán)伏安法詳解（下）

2019-10-25

Zeta電位的基本理論、測(cè)試方法和應(yīng)用

2020-08-24

項(xiàng)目推薦/Project

第一性原理-聲子譜

第一性原理-吸附能

熱門(mén)文章/popular

基礎(chǔ)理論丨一文了解XPS（概念、定性定量分析、分析方法、譜線結(jié)構(gòu)）

晶體結(jié)構(gòu)可視化軟件 VESTA使用教程（下篇）

手把手教你用ChemDraw 畫(huà)化學(xué)結(jié)構(gòu)式：基礎(chǔ)篇

【科研干貨】電化學(xué)表征：循環(huán)伏安法詳解（上）

電化學(xué)實(shí)驗(yàn)基礎(chǔ)之電化學(xué)工作站篇（二）三電極和兩電極體系的搭建和測(cè)試

【科研干貨】電化學(xué)表征：循環(huán)伏安法詳解（下）

微信掃碼分享文章

意見(jiàn)反饋

有獎(jiǎng)舉報(bào)

商務(wù)合作

...

更多

公眾號(hào)

關(guān)注我們了解更多

小程序

隨時(shí)預(yù)約掌握進(jìn)度

舉報(bào)有獎(jiǎng)

TEL: 191-3608-6524

如：在網(wǎng)絡(luò)上惡意使用“測(cè)試狗”等相關(guān)關(guān)鍵詞誤導(dǎo)用戶點(diǎn)擊、惡意盜用測(cè)試狗商標(biāo)、冒稱官方工作人員等情形，請(qǐng)您向我們舉報(bào)，經(jīng)查實(shí)后，我們將給予您獎(jiǎng)勵(lì)。

舉報(bào)內(nèi)容：

200

上傳附件：

文件格式不正確，請(qǐng)重新上傳文件格式不正確，請(qǐng)重新上傳文件格式不正確，請(qǐng)重新上傳

文件格式：jpg、jpeg、png、gif、tif、doc、docx、ppt、pptx、xls、xlsx、pdf、zip、rar

聯(lián)系方式

姓名

電話

提交意見(jiàn)

意見(jiàn)反饋

Suggestions

您可以在此留下您寶貴的意見(jiàn)，您的意見(jiàn)或問(wèn)題反饋將會(huì)成為我們不斷改進(jìn)的動(dòng)力。

意見(jiàn)類型

測(cè)試服務(wù)

網(wǎng)站功能

財(cái)務(wù)報(bào)賬

其他類型

意見(jiàn)內(nèi)容

200

聯(lián)系方式

姓名

電話

提交意見(jiàn)

收起

01

專屬信用額度，先測(cè)后付0元下單

02

下單享高額積分，萬(wàn)千好禮免費(fèi)兌

￥200

200元無(wú)門(mén)檻優(yōu)惠券

立即激活立即下單

已使用

已作廢

了解更多會(huì)員權(quán)益

全流程在線可視化，便捷高效觸手可及

如下單過(guò)程中有任何疑問(wèn)或需要幫助，請(qǐng)隨時(shí)咨詢專屬顧問(wèn)~

9折

9折無(wú)門(mén)檻優(yōu)惠券

支付一筆訂單后可領(lǐng)取

立即領(lǐng)取立即下單

已使用

已作廢

免費(fèi)測(cè)+驚喜盲盒+高額福利，多重福利大放送

新人免費(fèi)測(cè)

雙雙開(kāi)盲盒（100%中獎(jiǎng)）

邀請(qǐng)人專享

受邀人專享

立即邀請(qǐng)

創(chuàng)建/加入團(tuán)隊(duì)，解鎖定制化權(quán)益

01

1500元團(tuán)隊(duì)專屬優(yōu)惠券

02

萬(wàn)元大額信用額度，享先測(cè)后付

03

團(tuán)隊(duì)成員統(tǒng)一開(kāi)票報(bào)銷；

04

專業(yè)工程師課題專屬服務(wù)

領(lǐng)取成功，請(qǐng)下單

請(qǐng)您支付一筆訂單后才可以領(lǐng)取優(yōu)惠券

Document

<th id="6w0ik"></th>