双腿间已经湿成一片,亚洲啪啪综合AV一区,厨房的春潮A片,日本一道久久高清国产

<li id="7ox1o"><center id="7ox1o"><em id="7ox1o"></em></center></li>

<label id="7ox1o"><tt id="7ox1o"><big id="7ox1o"></big></tt></label>

<menu id="7ox1o"><form id="7ox1o"><thead id="7ox1o"></thead></form></menu>

<menu id="7ox1o"></menu>

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

熱線：400-152-6858

測試狗科研服務(wù)

預(yù)存免費試測登錄

Document

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

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

來源：時間：2022-12-16 16:36:30 瀏覽：2534次

銀吸附在鎳板上

我們建立一個由四個Ni原子組成的單層平板和一個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)

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

>>> from ase.visualize import view

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

>>> view(atoms)

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

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


>>>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)和水做一個界面。首先我們需要一層水。在這個腳本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)在運行WL.py腳本，從traj文件中讀取原子對象:


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

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

我們看看這個單位晶胞


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

我們需要一個鎳(111)平板，它與水盡可能接近。一個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])

看看這兩個單元晶胞，我們可以看到，如果我們把其中一個單元晶胞在平面上旋轉(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()方法僅在啟用周期性邊界條件時有效。我們在Ni(111)和水之間有2%的晶格失配，所以我們在平面上縮放水的晶胞以匹配平板的晶胞。參數(shù)scale_atoms=True表示原子位置應(yīng)該隨單元格縮放。默認(rèn)值是scale_atoms=False，表示單元格更改時笛卡爾坐標(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')

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

計算狗·模擬計算

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

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

相關(guān)文章

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

2019-10-25

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

2021-01-22

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

2021-06-19

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

2019-10-25

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

2020-05-03

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

2020-08-24

項目推薦/Project

第一性原理-聲子譜

第一性原理-吸附能

熱門文章/popular

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

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

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

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

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

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

微信掃碼分享文章

意見反饋

有獎舉報

商務(wù)合作

...

更多

公眾號

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

小程序

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

舉報有獎

TEL: 191-3608-6524

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

舉報內(nèi)容：

200

上傳附件：

文件格式不正確，請重新上傳文件格式不正確，請重新上傳文件格式不正確，請重新上傳

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

聯(lián)系方式

姓名

電話

提交意見

意見反饋

Suggestions

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

意見類型

測試服務(wù)

網(wǎng)站功能

財務(wù)報賬

其他類型

意見內(nèi)容

200

聯(lián)系方式

姓名

電話

提交意見

收起

01

專屬信用額度，先測后付0元下單

02

下單享高額積分，萬千好禮免費兌

￥200

200元無門檻優(yōu)惠券

立即激活立即下單

已使用

已作廢

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

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

如下單過程中有任何疑問或需要幫助，請隨時咨詢專屬顧問~

9折

9折無門檻優(yōu)惠券

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

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

已使用

已作廢

免費測+驚喜盲盒+高額福利，多重福利大放送

新人免費測

雙雙開盲盒（100%中獎）

邀請人專享

受邀人專享

立即邀請

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

01

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

02

萬元大額信用額度，享先測后付

03

團(tuán)隊成員統(tǒng)一開票報銷；

04

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

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

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

Document 大悟县| 永丰县| 宜君县| 东乡族自治县| 铁岭县| 大港区| 赞皇县| 双桥区| 柳河县| 南川市| 岢岚县| 封丘县| 拉萨市| 虎林市| 渑池县| 沾益县| 岳普湖县| 常熟市| 五华县| 正安县| 大兴区| 迁安市| 庄河市| 无为县| 博罗县| 阳山县| 延津县| 中方县| 舞阳县| 长乐市| 铁岭县| 托里县| 玉林市| 大庆市| 泾阳县| 咸丰县| 井陉县| 原阳县| 石渠县| 汾阳市| 禹州市|

<td id="9xs0w"><i id="9xs0w"><small id="9xs0w"></small></i></td>

<strike id="9xs0w"><form id="9xs0w"><kbd id="9xs0w"></kbd></form></strike>

<td id="9xs0w"><table id="9xs0w"></table></td>

<sup id="9xs0w"></sup>