 |
 |
 |
| 24小時(shí)熱門版塊排行榜 |
- >論壇更新日志 (352)
- >考研 (35)
- >蟲友互識(shí) (21)
- >導(dǎo)師招生 (11)
- >碩博家園 (6)
- >考博 (6)
- >文獻(xiàn)求助 (5)
- >文獻(xiàn)求助完結(jié)子版 (2)
- >教師之家 (2)
- >無(wú)機(jī)/物化 (2)
- >復(fù)合材料 (1)
- >招聘信息布告欄 (1)
- >論文道賀祈福 (1)
- >基金申請(qǐng) (1)
- >論文投稿 (1)
- >休閑灌水 (1)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 4198 | 回復(fù): 20 | |||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||||
ChemiAndy木蟲 (正式寫手)
|
[交流]
【原創(chuàng)】EMASS是CPMD最重要的參數(shù) 已有14人參與
|
||||
|
我們組今天舉辦了CPMD Workshop. Radu Iftime from UdeMon給了Lecture. 本人負(fù)責(zé)的是tutorial部分。不過我也是Beginner。Radu的Lecture給了我很大啟發(fā)。其中印象最深的是Emass是CPMD最重要的參數(shù)。 比較CPMD和BOMD，CPMD最重要的區(qū)別在于波函數(shù)(殼層電子)作為動(dòng)力學(xué)變量在盒子內(nèi)傳播。它有虛擬質(zhì)量fictitious mass，它有速度，自然也有動(dòng)能。那么這個(gè)動(dòng)能，Ekinetic，屬于一個(gè)虛擬量，憑空捏造出來(lái)的，目的是通過波函數(shù)的自然傳播，免去優(yōu)化波函數(shù)的步驟。（BOMD每一步都要優(yōu)化波函數(shù)）。這就產(chǎn)生一個(gè)問題，離子(原子核+內(nèi)層電子贗勢(shì)=離子)自然要傳遞能量給波函數(shù)這個(gè)運(yùn)動(dòng)員。畢竟在一起運(yùn)動(dòng)，碰撞嘛。那么能量在離子和電子之間的交換就是人為的，artifact。在實(shí)際體系中，這種傳遞是不存在的。Born-Openheimer說的就是原子核運(yùn)動(dòng)與電子運(yùn)動(dòng)的分離嘛。現(xiàn)在你耦合在一起，能量傳來(lái)傳去的，這不真實(shí)。而且，電子虛擬質(zhì)量越大，越不真實(shí)。這是CPMD最大的局限。本來(lái)嘛，想省事都是要付出代價(jià)的。 Radu組里面研究了這個(gè)虛擬質(zhì)量的影響。對(duì)于水，如果EMASS虛擬質(zhì)量選擇不當(dāng)，比如用780，那么這個(gè)g(r)和紅外光譜都要和真實(shí)值差好多，就是說對(duì)結(jié)構(gòu)和動(dòng)力學(xué)結(jié)果造成了嚴(yán)重的偏差。而如果EMASS用270, 就好很多，基本和真實(shí)值一致。 那是不是說選擇虛擬質(zhì)量越小越好呢？是的，那樣的話，coupling gap就會(huì)很大。那就能很好地實(shí)現(xiàn)離子和電子運(yùn)動(dòng)的分離，也就越接近于真實(shí)。但是，如果選擇較小的EMASS，也是要付出代價(jià)的。代價(jià)就是，。。。。(回帖可見) CPMD默認(rèn)的EMASS是400,很多人從來(lái)不改，我也是，因?yàn)槲也⒉恢酪趺锤�，改大點(diǎn)好，還是改小點(diǎn)好�，F(xiàn)在，Radu的研究就是告訴你，在可能的情況下越小越好。對(duì)于電子-離子耦合不嚴(yán)重的體系（比如，所有離子的質(zhì)量都很大，至少不包含H原子）你可以使用較大的EMASS，但你一定要警惕，如果你計(jì)算的是動(dòng)力學(xué)量，比如振動(dòng)光譜，或擴(kuò)散系數(shù)，你可能得到較大的偏差。 那么除了以上兩個(gè)參數(shù)，還有什么參數(shù)是重要的？當(dāng)然是Cutoff。這個(gè)相當(dāng)于Gauss里面的基組大小。直接決定了計(jì)算的效率。如果你看文獻(xiàn)，90%都用Cutoff 70+BLYP。但是如果你真做過驗(yàn)證，即計(jì)算Plot Cutoff~Energy的變化，去找那個(gè)Plateau，你就會(huì)發(fā)現(xiàn)，小于120Ry，根本就不可能。Radu的驗(yàn)證推薦，至少要用100Ry，才不會(huì)對(duì)動(dòng)力學(xué)結(jié)果產(chǎn)生嚴(yán)重偏差。 70Ry， 結(jié)構(gòu)量g(r)可能是對(duì)的，但振動(dòng)和擴(kuò)散就未必了。 以上討論全部基于MT贗勢(shì)，即Martin-Trulier Pseudopotential，cutoff一般要>70Ry 用超軟贗勢(shì)的話，一般贗勢(shì)選取25Ry。 超硬贗勢(shì)的話，一般Cutoff選擇140Ry。 謝謝valenhou001 指正。 [ Last edited by ChemiAndy on 2012-2-26 at 13:48 ] |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 模型 - 模擬 - 統(tǒng)計(jì) - 熱力學(xué) | 我學(xué)習(xí)計(jì)算的一些帖子 |
» 猜你喜歡
一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷
已經(jīng)有8人回復(fù)
-
求調(diào)劑
已經(jīng)有5人回復(fù)
-
求化學(xué)調(diào)劑
已經(jīng)有5人回復(fù)
-
085600，材料與化工321分求調(diào)劑
已經(jīng)有10人回復(fù)
-
0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分
已經(jīng)有5人回復(fù)
-
一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑
已經(jīng)有3人回復(fù)
-
085600，專業(yè)課化工原理，320分求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料與化工272求調(diào)劑
已經(jīng)有16人回復(fù)
-
290求調(diào)劑
已經(jīng)有3人回復(fù)
-
327求調(diào)劑
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- TiO2晶胞參數(shù)優(yōu)化問題
已經(jīng)有4人回復(fù)
- TiO2 晶胞參數(shù)優(yōu)化問題
已經(jīng)有3人回復(fù)
- 【美語(yǔ)怎么說 How to say it 】001-山寨
已經(jīng)有8人回復(fù)
- NPT系綜動(dòng)力學(xué)
已經(jīng)有19人回復(fù)
- 求廢舊塑料的主要參數(shù)指標(biāo)
已經(jīng)有15人回復(fù)
- 【求助】高嶺土模擬廢水濁度的測(cè)定
已經(jīng)有5人回復(fù)
- 【求助】2個(gè)cpmd問題，一個(gè)是時(shí)間步長(zhǎng)，一個(gè)是ficitious electronic mass！
已經(jīng)有11人回復(fù)
- 【求助】CPMD計(jì)算的問題
已經(jīng)有3人回復(fù)
|
[ Last edited by ChemiAndy on 2011-4-15 at 14:43 ] |
3樓2011-04-15 14:35:06
xk6891
至尊木蟲 (著名寫手)
- 模擬EPI: 1
- 應(yīng)助: 0 (幼兒園)
- 金幣: 13735.8
- 散金: 379
- 紅花: 9
- 帖子: 1514
- 在線: 546小時(shí)
- 蟲號(hào): 890342
- 注冊(cè): 2009-11-01
- 專業(yè): 理論和計(jì)算化學(xué)
2樓2011-04-15 11:36:45
ChemiAndy
木蟲 (正式寫手)
- 模擬EPI: 40
- 應(yīng)助: 52 (初中生)
- 金幣: 2573.3
- 散金: 18
- 紅花: 104
- 帖子: 590
- 在線: 553.8小時(shí)
- 蟲號(hào): 836907
- 注冊(cè): 2009-08-30
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
4樓2011-04-15 14:42:52
xk6891
至尊木蟲 (著名寫手)
- 模擬EPI: 1
- 應(yīng)助: 0 (幼兒園)
- 金幣: 13735.8
- 散金: 379
- 紅花: 9
- 帖子: 1514
- 在線: 546小時(shí)
- 蟲號(hào): 890342
- 注冊(cè): 2009-11-01
- 專業(yè): 理論和計(jì)算化學(xué)
就是就是 5樓2011-04-15 15:09:41
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 300求調(diào)劑，材料科學(xué)英一數(shù)二 +10 | leaflight 2026-03-24 | 10/500 |
|
|---|---|---|---|---|
|
|
[考研]
 一志愿武漢理工，總分321，英一數(shù)二，求老師收留。
+11
|
nnnnnnn5 2026-03-25 | 11/550 |
|
|
[基金申請(qǐng)] 面上5B能上會(huì)嗎？ +4 | redcom 2026-03-29 | 4/200 |
|
|
|
[考研] 289求調(diào)劑 +5 | BrightLL 2026-03-29 | 5/250 |
|
|
|
[考研] 材料學(xué)碩333求調(diào)劑 +11 | 北道巷 2026-03-24 | 11/550 |
|
|
|
[考研] 0856求調(diào)劑 +13 | zhn03 2026-03-25 | 14/700 |
|
|
|
[考研] 085701求調(diào)劑初試286分 +4 | secret0328 2026-03-28 | 4/200 |
|
|
|
[考研] 320分，材料與化工專業(yè)，求調(diào)劑 +9 | 一定上岸aaa 2026-03-27 | 13/650 |
|
|
|
[考研] 311求調(diào)劑 +9 | lin0039 2026-03-26 | 9/450 |
|
|
|
[考研] 274求調(diào)劑 +17 | 顧九笙要謙虛 2026-03-24 | 23/1150 |
|
|
|
[考研] 材料調(diào)劑 +8 | 匹克i 2026-03-23 | 8/400 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
|
[考研] 機(jī)械學(xué)碩310分，數(shù)一英一，一志愿211本科雙非找調(diào)劑信息 +3 | @357 2026-03-25 | 3/150 |
|
|
|
[考研] 289求調(diào)劑 +17 | 碩星赴 2026-03-23 | 17/850 |
|
|
|
[考研] 0854人工智能方向招收調(diào)劑 +4 | 章小魚567 2026-03-24 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|
|
|
[考研] 材料調(diào)劑 +3 | iwinso 2026-03-23 | 3/150 |
|
|
|
[考研] 一志愿吉林大學(xué)材料與化工303分求調(diào)劑 +4 | 為學(xué)666 2026-03-24 | 4/200 |
|
|
|
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +7 | 小羊睡著了? 2026-03-23 | 10/500 |
|