 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 6186 | 回復(fù): 17 | ||||||||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||||||||
lhtlx95新蟲 (小有名氣)
|
[求助]
用castep進(jìn)行結(jié)構(gòu)優(yōu)化的遇到的一些問題
|
|||||||
|
各位蟲友： 本人最近在做有關(guān)層狀硅酸鹽的第一性原理計算，應(yīng)用的軟件是MS中的castep，由于接觸第一性原理時間不長，在結(jié)構(gòu)優(yōu)化過程中遇到了一些問題，希望得到同行前輩的指點，小弟先謝過了！ 1、根據(jù)我看到的相關(guān)資料，在利用castep進(jìn)行計算時，需要對參數(shù)（主要是截斷能和K點）進(jìn)行收斂性測試。在論壇中一個帖子上，有同行提到“將體系總能量和晶格常數(shù)組為幾何優(yōu)化的目標(biāo)物理量，但能量截斷值大于某一能量值或k點取樣達(dá)到一定密度時，總能量和晶格常數(shù)將趨于一穩(wěn)定值”，請問各位蟲友，收斂性測試具體如何進(jìn)行？比如，通過什么具體的方法得到晶格常數(shù)對K點采樣的穩(wěn)定值。 2、在一篇相關(guān)的文獻(xiàn)中看到：“The pseudopotentials we used for Si, Al, and O were generated with the core electrons treated using the same Perdew -Wang 91 GGA as for the valence electrons in contrast to the common practice of generating pseudopotentials using the LDA then using the GGA approximation for valence electrons�！� 此論文也用的是castep方法，請問各位蟲友，在castep中如何選擇各元素的贗勢，是不是如下圖，在potentials選項卡中選擇，那么同一元素不同的贗勢有什么區(qū)別，以及如何選擇芯電子和價電子的贗勢？ 3、我所要計算的模型葉臘石的分子式為[Al4Si8O20(OH)4]，請問此結(jié)構(gòu)在做結(jié)構(gòu)優(yōu)化時是否需要考慮自旋極化，如果需要考慮的話，是否需要用LDA+U的方法以及initial spin如何設(shè)定？或者能夠推薦相關(guān)的資料或者書籍，謝謝。 4、結(jié)構(gòu)優(yōu)化后的結(jié)果和實驗結(jié)果對比的話，根據(jù)大家的經(jīng)驗，一般模擬的晶格參數(shù)與實驗結(jié)果相差多少時就認(rèn)為模擬的結(jié)果是合理的？ potentials.jpg |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 第一性原理和電化學(xué) | 量化計算有意義的討論 | 刺客 | good resources |
| 分子模擬 | xuexijisuan |
» 猜你喜歡
考研調(diào)劑
已經(jīng)有3人回復(fù)
-
317分 一志愿南理工材料工程 本科湖工大 求調(diào)劑
已經(jīng)有9人回復(fù)
-
297求調(diào)劑
已經(jīng)有10人回復(fù)
-
一志愿武漢理工，總分321，英一數(shù)二，求老師收留。
已經(jīng)有10人回復(fù)
-
343求調(diào)劑
已經(jīng)有5人回復(fù)
-
085601材料工程找調(diào)劑
已經(jīng)有7人回復(fù)
-
求調(diào)劑，一志愿 南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分
已經(jīng)有6人回復(fù)
-
一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷
已經(jīng)有5人回復(fù)
-
340求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿北化求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 固體物理和MS學(xué)習(xí)
已經(jīng)有28人回復(fù)
wangliangzhe
木蟲 (正式寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 2463.6
- 散金: 5
- 帖子: 598
- 在線: 252.9小時
- 蟲號: 302656
- 注冊: 2006-12-02
- 專業(yè): 區(qū)域環(huán)境化學(xué)
7樓2012-11-09 11:28:52
不在寂寞
鐵桿木蟲 (小有名氣)
- 1ST強(qiáng)帖: 1
- 應(yīng)助: 42 (小學(xué)生)
- 金幣: 7140
- 散金: 1
- 紅花: 9
- 帖子: 147
- 在線: 278.5小時
- 蟲號: 1273748
- 注冊: 2011-04-22
- 性別: GG
- 專業(yè): 無機(jī)非金屬類光電信息與功
【答案】應(yīng)助回帖
★
感謝參與，應(yīng)助指數(shù) +1
WDD880227: 金幣+1, 多謝交流 2012-11-07 19:03:47
感謝參與，應(yīng)助指數(shù) +1
WDD880227: 金幣+1, 多謝交流 2012-11-07 19:03:47
|
我也是剛剛接觸，說下我的理解 1、收斂測試是根據(jù)一系列幾何優(yōu)化，做出能量曲線，固定K點，測試不同截斷能下的能量，然后再固定截斷能，測試不同K點的能量。這樣做出曲線，找能量低的點或相對穩(wěn)定的點，不一定是能量最低的參數(shù)。 2、贗勢的選取不是固定的，你可以對同一體系進(jìn)行幾個函數(shù)的測試，然后選取與實驗值相對接近的，也可根據(jù)文獻(xiàn)中的做，看看別人都用什么，差別不是很大，主要是因為近似方法不同。另外，同一個體系用的是相同的函數(shù)，不是對每個原子進(jìn)行函數(shù)設(shè)置。 3、自旋極化是針對有磁性的體系，如果體系沒有磁性不必選擇。+U是對結(jié)果進(jìn)行修正的，當(dāng)理論結(jié)果與實驗結(jié)果相差較大時才考慮+U，而且U值也不是隨便加的。 4、晶格盡量接進(jìn)吧，這個沒什么標(biāo)準(zhǔn)，相差越小，說服力越強(qiáng)。 只是個人愚見  |
» 本帖已獲得的紅花（最新10朵）
2樓2012-11-07 13:27:41
wangliangzhe
木蟲 (正式寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 2463.6
- 散金: 5
- 帖子: 598
- 在線: 252.9小時
- 蟲號: 302656
- 注冊: 2006-12-02
- 專業(yè): 區(qū)域環(huán)境化學(xué)
3樓2012-11-07 14:39:28
lhtlx95
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 4024
- 帖子: 114
- 在線: 341小時
- 蟲號: 1622309
- 注冊: 2012-02-17
- 性別: GG
- 專業(yè): 固體力學(xué)
送鮮花一朵 |
謝謝您的指點，可我還有一些困惑： 1、有些資料上說收斂性測試是改變截斷能或K點對系統(tǒng)進(jìn)行單點能計算，你說是對系統(tǒng)進(jìn)行結(jié)構(gòu)優(yōu)化，到底哪種方法是正確的？而且，我做過比如固定截斷能，改變K點對同一系統(tǒng)進(jìn)行結(jié)構(gòu)優(yōu)化，如果總能趨于穩(wěn)定，但晶格常數(shù)卻變化很大，那該如何處理？就是說通過什么途徑的到晶格常數(shù)的穩(wěn)定值。 2、上圖給出的選項卡中，應(yīng)該就是對每個元素可以進(jìn)行贗勢的選擇，關(guān)鍵是什么時候需要選擇，選擇的標(biāo)準(zhǔn)是什么？ 3、您說自旋極化是針對有磁性的體系，請問如何判斷體系是否有磁性（不好意思，本人原本不是學(xué)化學(xué)的，可能問的問題比較低級），或者我給出的這個體系[Al4Si8O20(OH)4]是否具有磁性？ 4、晶格參數(shù)和實驗越接近就越好的話，除了調(diào)節(jié)截斷能、K點以及贗勢的種類，還可以通過調(diào)節(jié)哪些參數(shù)提高計算的準(zhǔn)確度，以及選擇這些參數(shù)有什么標(biāo)準(zhǔn)？ 請各位蟲友多多指教�。�！ |
4樓2012-11-08 09:22:42
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 317分 一志愿南理工材料工程 本科湖工大 求調(diào)劑 +9 | 芋泥小鈴鐺 2026-03-28 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 材料學(xué)碩333求調(diào)劑 +11 | 北道巷 2026-03-24 | 11/550 |
|
|
|
[考研] 300求調(diào)劑，材料科學(xué)英一數(shù)二 +8 | leaflight 2026-03-24 | 8/400 |
|
|
|
[考研] 085701求調(diào)劑初試286分 +4 | secret0328 2026-03-28 | 4/200 |
|
|
|
[考研] 343求調(diào)劑 +5 | 愛羈絆 2026-03-28 | 5/250 |
|
|
|
[考研] 083000學(xué)碩274求調(diào)劑 +8 | Li李魚 2026-03-26 | 8/400 |
|
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-27 | 6/300 |
|
|
|
[考研] 311求調(diào)劑 +3 | 希望上岸阿小楊 2026-03-23 | 3/150 |
|
|
|
[考研] 張芳銘-中國農(nóng)業(yè)大學(xué)-環(huán)境工程專碩-298 +4 | 手機(jī)用戶 2026-03-26 | 4/200 |
|
|
|
[考研] 材料與化工085600，總分304，本科有兩篇sci參與，求調(diào)劑 +10 | 幸運的醬醬 2026-03-22 | 12/600 |
|
|
|
[考研] 考研化學(xué)308分求調(diào)劑 +10 | 你好明天你好 2026-03-23 | 12/600 |
|
|
|
[考研] 一志愿211，335分，0856，求調(diào)劑院校和導(dǎo)師 +4 | 傾____蕭 2026-03-27 | 5/250 |
|
|
|
[考研] 求調(diào)劑，一志愿 南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +4 | @taotao 2026-03-26 | 5/250 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +4 | 李李不服輸 2026-03-26 | 4/200 |
|
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-25 | 3/150 |
|
|
|
[考研] 網(wǎng)絡(luò)空間安全0839招調(diào)劑 +4 | w320357296 2026-03-25 | 6/300 |
|
|
|
[考研] 材料調(diào)劑 +3 | iwinso 2026-03-23 | 3/150 |
|
|
|
[考研] B區(qū)考研調(diào)劑 +4 | yqdszhdap－ 2026-03-22 | 5/250 |
|
|
|
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3 | 開開心心沒煩惱 2026-03-23 | 4/200 |
|