太大了小峓子受不了春雨,久久精品a亚洲国产v高清不卡,最近2019手机中文字幕7,久久99国产精品久久99果冻传媒

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

熱線：400-152-6858

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

Document

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

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

來源：時(shí)間：2022-12-16 16:35:08 瀏覽：3474次

許多外部計(jì)算軟件可以與ASE一起使用，包括 GPAW, Abinit, Vasp, Siesta, Quantum ESPRESSO, Asap, LAMMPS 等等, 詳細(xì)請(qǐng)查閱支持計(jì)算軟件完整的清單.

使用ASE調(diào)用計(jì)算軟件

這里將介紹如何在ASE中調(diào)用外部計(jì)算軟件的方法。我們使用Vasp作為例子，其他計(jì)算軟件也可以以類似的方式使用。有關(guān)詳細(xì)信息，請(qǐng)參閱支持計(jì)算軟件的完整清單。

重要提示：ASE不提供VASP的代碼或許可，這些必須從其他地方獲得。ASE只使用VASP創(chuàng)建一個(gè)接口，由此將ASE提供的工具與VASP一起使用。

配置環(huán)境變量

首先，需要告訴ASE哪里能找到VASP的可執(zhí)行文件，以及在哪里找到贗勢(shì)。你需要定義兩個(gè)環(huán)境變量:


$ export ASE_VASP_COMMAND="mpirun -np XX $HOME/vasp/bin/vasp_std"
$ export VASP_PP_PATH=$HOME/vasp/mypps

第一個(gè)環(huán)境變量ASE_VASP_COMMAND是執(zhí)行VASP的默認(rèn)方式，應(yīng)該以與正常執(zhí)行VASP運(yùn)行相同的方式定義它。這里我們假設(shè)VASP可執(zhí)行文件vasp_std位于$HOME/VASP/bin中。注意，如果希望并行執(zhí)行VASP，這個(gè)調(diào)用還應(yīng)該包括MPI可執(zhí)行文件mpirun或mpiexec，在本例中為mpirun。第二個(gè)變量VASP_PP_PATH是到VASP贗勢(shì)庫(kù)的路徑。

vdw_kernel.bindat是一個(gè)額外(可選)變量文件，在進(jìn)行范德華計(jì)算時(shí)需要，其中luse_vdw=True。


$ export ASE_VASP_VDW=$HOME/<path-to-vdw_kernel.bindat-folder>

注意，這是針對(duì)文件夾，而不是文件本身。

運(yùn)行測(cè)試

ASE知道了如何執(zhí)行VASP，我們可以嘗試設(shè)置一個(gè)簡(jiǎn)單的計(jì)算。首先我們建立一個(gè)原子對(duì)象。首先我們建立原子對(duì)象。


from ase.build import molecule
atoms = molecule('N2')
atoms.center(vacuum=5)

為了執(zhí)行VASP DFT計(jì)算，我們現(xiàn)在設(shè)置了一個(gè)計(jì)算器對(duì)象。


from ase.calculators.vasp import Vasp
calc = Vasp(xc='pbe',  # 選擇交換關(guān)聯(lián)函數(shù)
            encut=400, # 平面波截?cái)?/span>
            kpts=(1, 1, 1)) # k點(diǎn)
atoms.calc = calc
en = atoms.get_potential_energy()  # 這個(gè)調(diào)用將開啟計(jì)算
print('Potential energy: {:.2f} eV'.format(en))

這將產(chǎn)生以下輸出:


Potential energy: -16.59 eV

ASE與VASP接口的流程是，ASE處理寫輸入文件，這是運(yùn)行所必需的，然后執(zhí)行ASE_VASP_COMMAND，即執(zhí)行VASP。一旦VASP運(yùn)行完成，ASE然后讀取所有相關(guān)文件，在本例中是OUTCAR、vasprun .xml和CONTCAR，并將屬性存儲(chǔ)在計(jì)算器對(duì)象中。有關(guān)VASP功能的更多信息，請(qǐng)參閱VASP。對(duì)于其他計(jì)算軟件請(qǐng)參考相關(guān)頁(yè)面。

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

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

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

相關(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電位的基本理論、測(cè)試方法和應(yīng)用

2020-08-24

項(xiàng)目推薦/Project

第一性原理-聲子譜

第一性原理-形成能

舉報(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

上傳附件：