国产成人午夜精品一区二区三区,成 人3d动漫在线观看,老师露出两个奶球让我吃奶头,啦啦啦中文免费视频高清观看

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

熱線:400-152-6858

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

預(yù)存 免費(fèi)試測(cè) 登錄
Document
當(dāng)前位置:文庫(kù)百科 ? 文章詳情
物理VS化學(xué)吸附?RDG,讓你的文章與眾不同
來(lái)源:科學(xué)10分鐘 時(shí)間:2021-01-27 10:50:12 瀏覽:11467次

引言

約化密度梯度(reduced density gradient, RDG)是密度泛函理論中用來(lái)描述電子非均勻性的無(wú)量綱參量[注1],其表達(dá)式為

2010年,楊偉濤教授提出基于RDG的非共價(jià)相互作用(noncovalent interactions, NCI)[1]分析方法,可以用來(lái)歸屬原子間或分子間的相互作用。雖然電子定域化函數(shù)(electron localization function,ELF)和分子中的原子(atoms-in-molecules, AIM)理論等能夠用來(lái)分析化學(xué)鍵,但是對(duì)氫鍵、空間位阻、π-π堆疊等非共價(jià)相互作用的討論有一定的局限性,而NCI分析則可以作為揭示這些非共價(jià)相互作用強(qiáng)有力的工具。因?yàn)榉肿拥碾娮用芏瘸尸F(xiàn)指數(shù)衰減,所以在遠(yuǎn)離分子處RDG呈現(xiàn)很大值;相反,在相互作用的區(qū)域RDG的值則很小,其中非共價(jià)相互作用區(qū)域的電子密度和RDG都比較低。

由于氫鍵和空間位阻區(qū)域的電子密度和RDG很接近而難以進(jìn)行區(qū)分,所以需要再借助通過(guò)電子密度Hessian矩陣的第二大本征值的正負(fù)號(hào)(sign(λ2))來(lái)進(jìn)一步判斷[注2];電子密度大小則可以反映出相互作用的強(qiáng)度,一般認(rèn)為vdW作用的區(qū)域ρ<0.005,因此可以通過(guò)對(duì)RDG等值面進(jìn)行sign(λ2)ρ填色或者做RDG vs sign(λ2)ρ的散點(diǎn)圖來(lái)討論所研究體系的相互作用類型。

在Sobereva的博文中已經(jīng)介紹了RDG相關(guān)的概念和使用Multiwfn對(duì)非周期邊界系進(jìn)行NCI分析的詳細(xì)過(guò)程[2]。但固體和表面以及牽扯過(guò)渡金屬的體系的NCI討論較少。因此,本文將通過(guò)若干實(shí)例的計(jì)算和相應(yīng)的分析,對(duì)這些體系使用Quantum ESPRESSO(QE)計(jì)算RDG的過(guò)程和NCI分析方法進(jìn)行介紹。


準(zhǔn)備工作

本文所采用的計(jì)算程序如下:

1. Quantum ESPRESSO 6.7 [3] 

實(shí)際上QE在5.1.0版本就已經(jīng)支持NCI分析,所以即使正在使用舊版本的用戶也不必?fù)?dān)心無(wú)法完成本文的工作。

2. VESTA 3.1.8 [4] 

個(gè)別新版本有bug,所以本文采用舊版本。

測(cè)試體系為:

1. 來(lái)源于Materials Project [5]的結(jié)構(gòu)

MoS2(mp-2815)和Ni(OH)2(mp-27912)

2. H2O2@PMCS來(lái)源于文獻(xiàn)[6]


計(jì)算參數(shù)

交換相關(guān)泛函

vdW-DF3-opt2

贗勢(shì)

超軟

GBRV

動(dòng)能截?cái)?/span>

40 Ry

密度截?cái)?/span>

480 Ry

k網(wǎng)格

MoS2

9×9×3

Ni(OH)2(Ni U=6.2 eV)

9×9×3

H2O2@PMCS

2×2×1

展寬

Gaussian

0.02 Ry

計(jì)算參數(shù)選擇的依據(jù):從文獻(xiàn)[7]測(cè)試結(jié)果來(lái)看,vdW-DF3-opt2在S22×5和S66×8的表現(xiàn)都還不錯(cuò),原文中所使用的是GBRV超軟贗勢(shì)[8],該贗勢(shì)本身也有很好的可靠性和可移植性[9,10]。建議密度截?cái)啵╡cutrho)選取充分以保證計(jì)算結(jié)果可靠以及避免RDG計(jì)算過(guò)程產(chǎn)生噪聲。


計(jì)算過(guò)程和數(shù)據(jù)處理

A1 使用PW模塊

進(jìn)行結(jié)構(gòu)優(yōu)化,并進(jìn)行自洽計(jì)算產(chǎn)生波函數(shù)以及電子密度

A2 通過(guò)PW/tools/pwo2xsf將out文件轉(zhuǎn)換為xsf格式[注3],如

A3 使用PP模塊進(jìn)行后處理,通過(guò)plot_num=19產(chǎn)生RDG格點(diǎn),以及 output_format=5來(lái)產(chǎn)生xsf文件

A4 類似過(guò)程A3,使用PP模塊通過(guò)plot_num=20產(chǎn)生sign(λ2)ρ的xsf文件,但建議提前將文件夾拷貝一份或者通過(guò)fileout更改xsf輸出格點(diǎn)文件的名稱避免上一步的結(jié)果被覆蓋。

B1 打開(kāi)VESTA程序,將A2中out文件產(chǎn)生的xsf文件拖入VESTA窗體

B2菜單欄中Edit->Edit Data->Volumetric Data,分別點(diǎn)擊Isosurface和Surface coloring旁的Import分別導(dǎo)入RDG網(wǎng)格和sign(λ2)ρ的xsf文件。

B3 菜單欄中Objects->Properties->Isosurfaces,將Isosurfaces中的Isosurface level改為0.5,以及將Surface coloring中的Saturation level的Max改為0.02,Min改為-0.04,這樣能夠和文獻(xiàn)[1]較好符合。

B4 調(diào)整其他顯示部分,比如菜單中去除勾選Objects-> Volumetric Data->Show Sections,以及在Objects->Properties->Atoms對(duì)原子顏色進(jìn)行調(diào)整,在Edit->Bonds對(duì)顯示成鍵的判斷進(jìn)行調(diào)整,使得作圖更加美觀,但具體調(diào)整的細(xì)節(jié)這里不進(jìn)行贅述。如果對(duì)操作過(guò)程有問(wèn)題也歡迎大家討論。


案例

MoS2是層狀結(jié)構(gòu)材料,MoS2層間出現(xiàn)了RDG等值面圍成的區(qū)域,且層間S……S連線的區(qū)域呈現(xiàn)淡綠色(ρ<0.005),證明了層間存在相互作用,且通過(guò)vdW作用相結(jié)合并主要是以層間S……S的相互作用為主導(dǎo);而層內(nèi)的三個(gè)Mo原子之間顯示為紅色,表現(xiàn)出位阻效應(yīng)。

Ni(OH)2與MoS2情況類似,但Ni(OH)2層間作用以O(shè)-H……O-H的氫鍵作用為主導(dǎo);層內(nèi)的Ni原子之間也呈現(xiàn)出位阻效應(yīng)。在Ni-O鍵中形成的藍(lán)色環(huán)形區(qū)域是因?yàn)镹i-O之間相互作用比較強(qiáng),而pp.x程序在處理RDG的時(shí)候?qū)ⅵ?gt;0.05的網(wǎng)格直接設(shè)置為100[注4],導(dǎo)致共價(jià)鍵或離子鍵在圖中不被顯示。

H2O2@PMCS這一類結(jié)構(gòu)常見(jiàn)于單原子催化過(guò)程中,也是我們計(jì)算較多的一類結(jié)構(gòu)。文獻(xiàn)[6]給出的吸附能為-0.45 eV,僅從吸附能來(lái)看可能是物理吸附。盡管本文采用的方法與文獻(xiàn)不同,但給出了-0.47 eV的結(jié)果與文獻(xiàn)[6]十分接近。

為了進(jìn)一步確定H2O2與基底相互作用的類型,我們進(jìn)行了NCI分析,結(jié)果表明過(guò)氧化氫與Zn-N4的相互作用卻是是非共價(jià)相互作用,其中H2O2中一個(gè)O-H與Zn-N4中的N原子形成了氫鍵;而另一個(gè)O-H的O與Zn原子產(chǎn)生了靜電作用,這種靜電作用從顏色上看要更強(qiáng)于氫鍵作用,因此可以認(rèn)為H2O2在PMCS的吸附作用是由這種靜電作用所主導(dǎo)的物理吸附。


Promolecule近似

Promolecule近似就是將球?qū)ΨQ的自由原子密度直接進(jìn)行疊加,即

Promolecule近似下的密度通??梢宰鳛樽郧?chǎng)的初始猜測(cè),也能夠用于大體系的定性計(jì)算。在QE中也能實(shí)現(xiàn)Promolecule近似的計(jì)算,這是因?yàn)镼E中UPF格式的贗勢(shì)文件中已經(jīng)存有自由原子的電子密度。在自洽場(chǎng)計(jì)算過(guò)程中,我們只進(jìn)行初始猜測(cè)而不做自洽過(guò)程,因此將electron_maxstep設(shè)置為0,且startingwfc設(shè)置為’atomic’,就可以得到Promolecule近似下的密度。

由于Promolecule近似做NCI分析也能給出定性可靠的結(jié)果[1],因此在處理更大的體系時(shí)可以使用這種方法。這里仍然以H2O2@PMCS體系為例,展示Promolecule近似下NCI分析的結(jié)果,下圖表明Promolecule近似(等值面取0.035)與自洽計(jì)算后的差距并不大。


散點(diǎn)圖

由于非共價(jià)相互作用區(qū)域電子密度和RDG都比較低,且λ2的符號(hào)為負(fù),通常在散點(diǎn)圖中呈現(xiàn)為一個(gè)或多個(gè)尖峰(spike)。如H2O2@PMCS的散點(diǎn)圖中,左下角橫坐標(biāo)位于-0.04~0.02的spike就對(duì)應(yīng)于H2O2中OH與N的氫鍵。

而promolecule近似下的散點(diǎn)圖與上圖的結(jié)果較為接近

那么如何繪制這樣的散點(diǎn)圖呢?我在這里提供了一個(gè)自己寫的腳本nci_scatter.py能夠直接繪制散點(diǎn)圖。需要預(yù)先安裝python 3.x,調(diào)用到的python組件有numpy和matplotlib,只要通過(guò)pip install就能夠進(jìn)行安裝。只要執(zhí)行nci_scatter.py腳本按照每一步的提示用鼠標(biāo)拖入對(duì)應(yīng)的xsf文件就能夠直接畫出上面的散點(diǎn)圖。

“nci_scatter.py”

鏈接:

https://pan.baidu.com/s/1rtUFsfvlFPtLiqROc2qNHw

提取碼:csgo

請(qǐng)大家請(qǐng)自行下載~


總結(jié)

Quantum ESPRESSO是一款功能十分強(qiáng)大且開(kāi)源免費(fèi)的第一性原理計(jì)算程序,NCI分析也是一種強(qiáng)大的非共價(jià)相互作用分析方法。QE已經(jīng)內(nèi)置了RDG的計(jì)算,但卻很少有人寫QE做NCI分析的教程,本文以固體和表面中的一些常見(jiàn)體系為案例,對(duì)實(shí)現(xiàn)過(guò)程和分析進(jìn)行了一些敘述討論,更多的理論基礎(chǔ)和分析細(xì)節(jié)還是建議大家多研讀文獻(xiàn)。希望本文能夠?yàn)榇蠹业目蒲泄ぷ髦μ聿剩瑫r(shí)也歡迎大家共同參與討論互相進(jìn)步。


注釋

[1] RDG是無(wú)量綱參量,這是因?yàn)?/span>

[2] 這些概念來(lái)源于AIM理論,電子密度的Hessian矩陣是一個(gè)3×3的矩陣,對(duì)角化后得到三個(gè)本征值λ1>λ2>λ3。若λ2<0意味著λ3也小于0,對(duì)應(yīng)于鍵臨界點(diǎn)(bond critical point,BCP);若λ2>0意味著λ1也大于0,對(duì)應(yīng)于環(huán)臨界點(diǎn)(ring critical point, RCP)。

[3] pwo2xsf的作用是將pw計(jì)算得到的out文件轉(zhuǎn)化為xsf文件,可選選項(xiàng)有

-ic 初始結(jié)構(gòu)

-lc 最終結(jié)構(gòu)

-oc 最優(yōu)結(jié)構(gòu)

-a 將軌跡轉(zhuǎn)化為AXSF

[4] 嚴(yán)謹(jǐn)起見(jiàn),我們檢查了QE的pp.x程序的源代碼,發(fā)現(xiàn)在QE 6.x版本分析磁性體系的sign(λ2)ρ和ELF時(shí),電子密度存在bug(這是因?yàn)樽孕龢O化的計(jì)算中,QE 6.x版本的電子密度不是按照spin up和spin down進(jìn)行存儲(chǔ)的,而是按照total和spin進(jìn)行存儲(chǔ)的,但sign(λ2)ρ和ELF這兩段代碼忘記進(jìn)行修改)。這點(diǎn)我們已經(jīng)向QE官方提交了這個(gè)bug并得到了確認(rèn),所以先建議大家不要分析磁性體系,后續(xù)版本會(huì)進(jìn)行改進(jìn)。本文中的案例Ni(OH)2的結(jié)果展示是我們使用修改后的程序,因此結(jié)果可靠。


參考資料

[1] Johnson E R, Keinan S, Mori-Sánchez P, et al. Revealing noncovalent interactions[J]. Journal of the American Chemical Society, 2010, 132(18): 6498-6506.

DOI: 10.1021/ja100936w

[2] http://sobereva.com/68

[3] https://github.com/QEF/q-e/releases

[4] http://www.jp-minerals.org/vesta/en/download.html

[5] https://www.materialsproject.org

[6] Xu B, Wang H, Wang W, et al. A Single‐Atom Nanozyme for Wound Disinfection Applications[J]. Angewandte Chemie, 2019, 131(15): 4965-4970.

DOI: 10.1002/anie.201813994

[7] Chakraborty D, Berland K, Thonhauser T. Next-Generation Nonlocal van der Waals Density Functional[J]. Journal of Chemical Theory and Computation, 2020, 16(9): 5893-5911.

DOI: 10.1021/acs.jctc.0c00471

[8] Garrity K F, Bennett J W, Rabe K M, et al. Pseudopotentials for high-throughput DFT calculations[J]. Computational Materials Science, 2014, 81: 446-452.

DOI: 10.1016/j.commatsci.2013.08.053

[9] Marsman M, Marzari N, Nitzsche U, et al. Reproducibility in density functional theory calculations of solids[J]. Science, 2016, 351: 6280.

DOI: 10.1126/science.aad3000

[10] https://molmod.ugent.be/deltacodesdft


學(xué)術(shù)交流微信群


測(cè)試狗同步輻射咨詢預(yù)約群 / 同步輻射GIWAXS/GISAXS交流群 / OPV鈣鈦礦制備與表征交流群 / 測(cè)試GO-XRD精修交流群 / 催化理論與實(shí)驗(yàn)交流群 / 分子模擬應(yīng)用交流群 / 電催化交流群 / 測(cè)試狗推文創(chuàng)作群/CSG學(xué)術(shù)干貨分享群 / CSG能源材料交流
加微信群方式:添加GO小妹微信15680750806,
請(qǐng)備注“昵稱—學(xué)校/機(jī)構(gòu)”


學(xué)術(shù)交流QQ群


測(cè)試狗EBSD學(xué)習(xí)交流QQ群:597369268
機(jī)器學(xué)習(xí)+催化計(jì)算QQ群:966186610
Lammps分子模擬實(shí)踐交流群:1137524036
材料測(cè)試&模擬計(jì)算QQ群:776761012
科研字幕組(科學(xué)10分鐘)QQ群:703443384
測(cè)試GO軟件資源共享QQ群:567425397
評(píng)論 / 文明上網(wǎng)理性發(fā)言
12條評(píng)論
全部評(píng)論 / 我的評(píng)論
最熱 /  最新
全部 3小時(shí)前 四川
文字是人類用符號(hào)記錄表達(dá)信息以傳之久遠(yuǎn)的方式和工具。現(xiàn)代文字大多是記錄語(yǔ)言的工具。人類往往先有口頭的語(yǔ)言后產(chǎn)生書面文字,很多小語(yǔ)種,有語(yǔ)言但沒(méi)有文字。文字的不同體現(xiàn)了國(guó)家和民族的書面表達(dá)的方式和思維不同。文字使人類進(jìn)入有歷史記錄的文明社會(huì)。
點(diǎn)贊12
回復(fù)
全部
查看更多評(píng)論
相關(guān)文章

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

2020-05-03

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

2021-06-19

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

2021-01-22

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

2019-10-25

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

2019-10-25

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

2020-08-24

項(xiàng)目推薦/Project
化學(xué)吸附儀(TPD/TPR/TPO)

化學(xué)吸附儀(TPD/TPR/TPO)

腐蝕性氣體吸附分析儀(BET)

腐蝕性氣體吸附分析儀(BET)

全自動(dòng)高溫高壓氣體吸附儀

全自動(dòng)高溫高壓氣體吸附儀

第一性原理-吸附能

第一性原理-吸附能

熱門文章/popular

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

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

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

【科研干貨】電化學(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

下單享高額積分,積分當(dāng)錢花

200
200元無(wú)門檻優(yōu)惠券
立即激活 立即下單
已使用
已作廢

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

如下單過(guò)程中有任何疑問(wèn)或需要幫助,請(qǐng)隨時(shí)咨詢專屬顧問(wèn)~
9
9折無(wú)門檻優(yōu)惠券

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

立即領(lǐng)取 立即下單
已使用
已作廢

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

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

雙雙開(kāi)盲盒(100%中獎(jiǎ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
關(guān)于我們 新手幫助 測(cè)試干貨 商務(wù)合作 基金查詢 相關(guān)資質(zhì) 模擬計(jì)算 現(xiàn)場(chǎng)測(cè)試 服務(wù)項(xiàng)目 科研繪圖 同步輻射 電池行業(yè)

聯(lián)系方式/contact

400-152-6858

工作時(shí)間/work time

09:00-18:00

測(cè)試狗公眾號(hào)

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

測(cè)試狗小程序

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

蜀公網(wǎng)安備51010602000648號(hào)

蜀ICP備17005822號(hào)-1

成都世紀(jì)美揚(yáng)科技有限公司

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

业余 自由 性别 视频 视频| 蜜臀av人妻久久无码精品麻豆| 成人做爰高潮a片免费视频| 被多个强壮的黑人灌满浆| 紫川动漫在线观看免费完整版| 天天躁日日躁狠狠躁免费麻豆| 狠狠色丁香婷婷久久综合麻豆| 扒开她粉嫩的小缝的a片| 国产黄a三级三级三级看三级| 日本三级韩国三级韩级| 精品久久久久久亚洲偷窥一区| 久久久久亚洲av无码专区体验| 玩弄高耸白嫩的乳峰a片| 久久热精品视频| 超碰视频在线| 东北浪妇交换俱乐部| 国产乱码精品一区二区三区久久| 美艳人妻办公室抽搐呻吟| 蜜芽va亚洲va欧美va天堂| 白嫩少妇激情无码| 久久久欧美国产精品人妻噜噜| 大学生第一次破苞疼哭了| 亚洲色偷精品一区二区三区| 小雪第一次尝到又大又粗| 教官脔到她哭h粗话h好爽视频| 美女被c到爽哭视频网站| 亚洲av无码国产综合专区| 丰满人妻中伦妇伦精品app| 伸进去摸老妇的毛| 少妇高潮a片无套内谢麻豆传| 玩弄懦弱美人双xing大nai| 久久综合亚洲色hezyo国产| 三个人弄了我一晚上| 爱爱爱爱看视频| 国内精品卡一卡二卡三| 全国最大成人网站| 日本无码成人片在线观看波多| 国精产品一区二区三区公司| 丝瓜视频在线观看| 久久日产一线二线三线| 调教超级yin荡玩物学生男男|