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

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

熱線：400-152-6858

測試狗科研服務(wù)

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

Document

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

ASE原子模擬環(huán)境之計算氮分子在銅表面的吸附能

來源：時間：2022-12-16 16:32:16 瀏覽：3111次

本節(jié)簡要(但不完整)地概述ASE可以做什么。我們將計算氮分子在銅表面的吸附能。這是通過計算單獨平板和單獨分子的總能量來實現(xiàn)的。然后將吸附劑加入板中并使其松弛，計算出該復(fù)合體系的總能量。吸附能為孤立能減去復(fù)合體系能量之和。這是重新平衡后的系統(tǒng)圖片:

請看下面的源碼：N2Cu.py



from ase import Atoms
from ase.calculators.emt import EMT
from ase.constraints import FixAtoms
from ase.optimize import QuasiNewton
from ase.build import fcc111, add_adsorbate
h = 1.85
d = 1.10
slab = fcc111('Cu', size=(4, 4, 2), vacuum=10.0)
slab.calc = EMT()
e_slab = slab.get_potential_energy()
molecule = Atoms('2N', positions=[(0., 0., 0.), (0., 0., d)])
molecule.calc = EMT()
e_N2 = molecule.get_potential_energy()
add_adsorbate(slab, molecule, h, 'ontop')
constraint = FixAtoms(mask=[a.symbol != 'N' for a in slab])
slab.set_constraint(constraint)
dyn = QuasiNewton(slab, trajectory='N2Cu.traj')
dyn.run(fmax=0.05)
print('Adsorption energy:', e_slab + e_N2 - slab.get_potential_energy())

假設(shè)在你的python和ASE的環(huán)境變量設(shè)置正確，運行源碼文件：



python N2Cu.py

下面我們對代碼進行講解。

Atoms

Atoms對象是原子的集合。下面是通過直接指定兩個氮原子的位置來定義N2分子:



>>> from ase import Atoms
>>> d = 1.10
>>> molecule = Atoms('2N', positions=[(0., 0., 0.), (0., 0., d)])

你還可以構(gòu)建晶體，例如使用晶體模塊來構(gòu)建回與常見晶體結(jié)構(gòu)對應(yīng)的Atoms對象。讓我們做一個Cu(111)表面:

>>> from ase.build import fcc111

>>> slab = fcc111('Cu', size=(4,4,2), vacuum=10.0)

Calculators

ASE調(diào)用計算軟件可以使用ase.calculators，包括emt、Asap、Dacapo、GPAW、Abinit、Vasp等等。有關(guān)完整列表，請參見ASE主頁。

這里我們使用有效介質(zhì)理論(EMT)計算器，因為它非?？?，因此對于入門非常有用。我們可以在之前創(chuàng)建的Atoms對象上附加一個計算器：



>>> from ase.calculators.emt import EMT
>>> slab.calc = EMT()
>>> molecule.calc = EMT()

使用Atoms類中的get_potential_energy()方法計算系統(tǒng)的總能量:

>>> e_slab = slab.get_potential_energy()

>>> e_N2 = molecule.get_potential_energy()

結(jié)構(gòu)優(yōu)化

利用QuasiNewton優(yōu)化器對吸附在銅表面的N2分子結(jié)構(gòu)進行優(yōu)化。首先將吸附劑添加到銅板上，例如在頂部位置:

>>> h = 1.85

>>> add_adsorbate(slab, molecule, h, 'ontop')

為了加快弛豫速度，我們利用constraints模塊中的FixAtoms使銅原子固定。只有N2分子被允許弛豫到平衡結(jié)構(gòu):

>>> from ase.constraints import FixAtoms

>>> constraint = FixAtoms(mask=[a.symbol != 'N' for a in slab])

>>> slab.set_constraint(constraint)

現(xiàn)在將QuasiNewton優(yōu)化器應(yīng)用到系統(tǒng)并保存軌跡文件。運行優(yōu)化器直到體系達(dá)到收斂標(biāo)準(zhǔn)，即所有原子上的力應(yīng)小于某個fmax:

>>> from ase.optimize import QuasiNewton

>>> dyn = QuasiNewton(slab, trajectory='N2Cu.traj')

>>> dyn.run(fmax=0.05)

輸入-輸出

將原子位置寫入文件的過程是使用write()函數(shù)完成的:

>>> from ase.io import write

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

這將寫入一個xyz格式的文件。其他格式有：

格式	描述
xyz	簡單的xyz格式
cube	高斯cube文件
pdb	蛋白質(zhì)數(shù)據(jù)庫文件
traj	ASE自有軌跡格式
py	Python腳本

從文件讀取這樣做：



>>> from ase.io import read
>>> slab_from_file = read('slab.xyz')

如果文件包含多個構(gòu)型，write()函數(shù)的默認(rèn)寫入最后一個構(gòu)型。但是，我們可以通過以下方式加載特定的構(gòu)型:

>>> read('slab.traj')      # 最后一個構(gòu)型

>>> read('slab.traj', -1)  # 同上

>>> read('slab.traj', 0)   # 第一個構(gòu)型

可視化

可視化原子最簡單的方法是view()函數(shù)：



>>> from ase.visualize import view
>>> view(slab)

這將彈出一個ase.gui窗口。其他的可視化軟件可以通過指定可選關(guān)鍵字viewer=…-選擇使用'ase.gui', ' gopenmol ', ' vmd '，或' rasmol '中的一個。(注意，這些替代查看器不是ASE的一部分，需要由用戶單獨安裝。)VMD查看器可以采取一個可選的data參數(shù)來顯示三維數(shù)據(jù):

>>> view(slab, viewer='VMD', data=array)

分子動力學(xué)

我們看看一個使用VelocityVerlet進行分子動力學(xué)計算應(yīng)用到氮氣分子上的例子。我們首先創(chuàng)建一個VelocityVerlet對象，給出它的分子和牛頓定律積分的時間步長。然后，我們通過調(diào)用它的run()方法來進行分子動力學(xué)計算，并給出要采取的步驟數(shù):

>>> from ase.md.verlet import VelocityVerlet

>>> from ase import units

>>> dyn = VelocityVerlet(molecule, dt=1.0 * units.fs)

>>> for i in range(10):

...     pot = molecule.get_potential_energy()

...     kin = molecule.get_kinetic_energy()

...     print('%2d: %.5f eV, %.5f eV, %.5f eV' % (i, pot + kin, pot, kin))

...     dyn.run(steps=20)

上一篇：Python中的數(shù)字型數(shù)組（Numpy）教程

下一篇：ASE原子環(huán)境模擬之原子構(gòu)建和計算單元

相關(guān)文章

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

2020-05-03

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

2021-01-22

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

2021-06-19

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

2019-10-25

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

2019-10-25

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

2020-08-24

項目推薦/Project

第一性原理-自旋密度

第一性原理-鋰硫路徑

第一性原理-能帶

熱門文章/popular

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

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

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

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

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

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

微信掃碼分享文章

意見反饋

有獎舉報

商務(wù)合作

...

更多

公眾號

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

小程序

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

舉報有獎

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

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

意見類型

測試服務(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)建/加入團隊，解鎖定制化權(quán)益

01

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

02

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

03

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

04

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

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

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

Document