 |
 |
 |
| 24小時(shí)熱門版塊排行榜 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2993 | 回復(fù): 7 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
lhtlx95新蟲 (小有名氣)
|
[求助]
用castep計(jì)算單斜晶系彈性常數(shù)的問題 已有1人參與
|
||
|
最近利用castep計(jì)算某一層狀硅酸鹽（單斜晶系）的彈性常數(shù)，計(jì)算過程中遇到一些問題，想請(qǐng)大家指教！ 1、castep中彈性常數(shù)的算法是在給定的一系列的應(yīng)變模式下進(jìn)行結(jié)構(gòu)優(yōu)化，而且是只優(yōu)化原子位置，而不優(yōu)化晶格常數(shù)。然后計(jì)算應(yīng)力，再利用應(yīng)力應(yīng)變進(jìn)行擬合，得到彈性常數(shù)。（不知道我理解的對(duì)不對(duì)，請(qǐng)大家指教）。而對(duì)于單斜晶系，晶角beta不是90度，所以計(jì)算出的彈性常數(shù)應(yīng)該和坐標(biāo)系的選擇有關(guān)系。我測(cè)試選擇了兩種坐標(biāo)系（在orientation convention里修改），一種是C along Z,A in XZ plane（即Y和AC平面垂直），另一種是A along X,B in XY plane（即Z和AB平面垂直），顯然兩種坐標(biāo)系Z方向不同。然后分別進(jìn)行彈性模量的計(jì)算，但我發(fā)現(xiàn)在同一應(yīng)變模式，比如Z方向加應(yīng)變0.003，兩種坐標(biāo)系下加應(yīng)變后，晶格長(zhǎng)度完全相同�。。� C along Z,A in XZ plane坐標(biāo)系下elastic constants 21 cpingxingz 6 11.7_cij__3__1.castep文件中的晶格常數(shù)： Lattice parameters(A) Cell Angles a = 5.175224 alpha = 91.304789 b = 8.993333 beta = 101.075310 c = 9.368014 gamma = 89.792157 A along X,B in XY plane坐標(biāo)系下elastic constants 21 apingxingx 6 11.7_cij__3__1.castep文件中的晶格常數(shù)： Lattice parameters(A) Cell Angles a = 5.175224 alpha = 91.304789 b = 8.993333 beta = 101.075310 c = 9.368014 gamma = 89.792157 但我認(rèn)為在不同坐標(biāo)系下，都給Z方向加0.003的應(yīng)變后，晶格常數(shù)應(yīng)該有所不同！所以我懷疑是不是我坐標(biāo)系設(shè)置沒有對(duì)! 但進(jìn)一步看，不同坐標(biāo)系下，都給Z方向加0.003的應(yīng)變后，計(jì)算得到的應(yīng)力又不相同。 所以很困惑，castep中應(yīng)變是按照什么原理加的？？？？ 2、計(jì)算彈性模量過程中如果截?cái)嗄苓x600eV，收斂的比較好！如下： +-----------+-----------------+-----------------+------------+-----+ <-- BFGS | Parameter | value | tolerance | units | OK? | <-- BFGS +-----------+-----------------+-----------------+------------+-----+ <-- BFGS | dE/ion | 1.452405E-008 | 1.000000E-006 | eV | Yes | <-- BFGS | |F|max | 1.271911E-003 | 2.000000E-003 | eV/A | Yes | <-- BFGS | |dR|max | 3.237645E-005 | 1.000000E-004 | A | Yes | <-- BFGS +-----------+-----------------+-----------------+------------+-----+ <-- BFGS BFGS: Geometry optimization completed successfully. 但如果其他參數(shù)不變，只是把截?cái)嗄芴岣叩?000eV，優(yōu)化的結(jié)果就不收斂，還出現(xiàn)了以下的警告： BFGS: Warning - looks like this system is as converged as possible. Maybe your geometry convergence tolerances are too tight? BFGS: finished iteration 23 with enthalpy= -1.16912502E+004 eV +-----------+-----------------+-----------------+------------+-----+ <-- BFGS | Parameter | value | tolerance | units | OK? | <-- BFGS +-----------+-----------------+-----------------+------------+-----+ <-- BFGS | dE/ion | 8.287700E-008 | 1.000000E-006 | eV | Yes | <-- BFGS | |F|max | 2.499152E-003 | 2.000000E-003 | eV/A | No | <-- BFGS | |dR|max | 4.845477E-004 | 1.000000E-004 | A | No | <-- BFGS +-----------+-----------------+-----------------+------------+-----+ <-- BFGS BFGS: Geometry optimization completed successfully. 請(qǐng)問大家這是為什么？好像彈性模量對(duì)截?cái)嗄芎蚄點(diǎn)的要求比較高，為什么提高截?cái)嗄芊炊皇諗浚?br /> 3、我計(jì)算的彈性模量結(jié)果（先按三斜晶系計(jì)算）： id i j Cij (GPa) 1 1 1 187.87374 +/- 3.438 2 2 2 195.89126 +/- 3.629 3 3 3 62.89919 +/- 3.909 4 4 4 38.06517 +/- 1.711 5 5 5 29.04281 +/- 2.239 6 6 6 77.15969 +/- 0.653 7 1 2 47.91224 +/- 1.692 8 1 3 11.66727 +/- 2.152 9 1 4 -0.98254 +/- 0.958 10 1 5 -32.19545 +/- 1.928 11 1 6 -9.58425 +/- 0.946 12 2 3 6.24273 +/- 1.583 13 2 4 5.91387 +/- 1.170 14 2 5 -13.29736 +/- 2.272 15 2 6 2.85261 +/- 1.039 16 3 4 4.28390 +/- 1.030 17 3 5 5.48421 +/- 2.970 18 3 6 -2.06920 +/- 0.943 19 4 5 -3.76064 +/- 1.397 20 4 6 -10.32344 +/- 0.630 21 5 6 -3.20257 +/- 1.044 后面的偏差和文獻(xiàn)比較起來偏大，不知該如何調(diào)整？ 4、我看到論壇里有蟲友說可以用提腳本文件的方法計(jì)算彈性模量，不知利用腳本文件和直接計(jì)算有什么不同？ |
回復(fù)此樓» 猜你喜歡
數(shù)一英一271專碩（085401）求調(diào)劑，可跨
已經(jīng)有8人回復(fù)
-
295求調(diào)劑
已經(jīng)有4人回復(fù)
-
0856求調(diào)劑
已經(jīng)有5人回復(fù)
-
材料學(xué)碩333求調(diào)劑
已經(jīng)有10人回復(fù)
-
299求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿鄭州大學(xué)，080500學(xué)碩，總分317分求調(diào)劑
已經(jīng)有6人回復(fù)
-
289求調(diào)劑
已經(jīng)有12人回復(fù)
-
070305高分子化學(xué)與物理 304分求調(diào)劑
已經(jīng)有6人回復(fù)
-
0703化學(xué)調(diào)劑，求導(dǎo)師收
已經(jīng)有10人回復(fù)
-
各位老師好，我的一志愿為北京科技大學(xué)085601材料專碩
已經(jīng)有5人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 用CASTEP計(jì)算彈性常數(shù)出錯(cuò)的原因
已經(jīng)有5人回復(fù)
- E-V曲線擬合后怎么找到能量最低點(diǎn)對(duì)應(yīng)的V值
已經(jīng)有4人回復(fù)
- Castep優(yōu)化結(jié)果的晶格常數(shù)最終結(jié)果怎么看？
已經(jīng)有13人回復(fù)
- 求助審稿人提出的有關(guān)晶格優(yōu)化的問題
已經(jīng)有17人回復(fù)
- 關(guān)于六角晶系的彈性常數(shù)計(jì)算
已經(jīng)有16人回復(fù)
- 晶格常數(shù)變化的討論c/a 還是 b/a?
已經(jīng)有9人回復(fù)
- CASTEP 計(jì)算介電常數(shù) 數(shù)值的意義
已經(jīng)有6人回復(fù)
- 為什么采用CASTEP計(jì)算CaTiO3的光學(xué)性質(zhì)其中的靜態(tài)介電常數(shù)與實(shí)驗(yàn)值相差很大
已經(jīng)有14人回復(fù)
- CASTEP計(jì)算彈性常數(shù)C11，C12，C44嚴(yán)重偏大
已經(jīng)有13人回復(fù)
- 計(jì)算不同晶系的彈性常數(shù)
已經(jīng)有260人回復(fù)
- 【求助】castep計(jì)算彈性常數(shù)問題
已經(jīng)有14人回復(fù)
- 【求助】請(qǐng)教手動(dòng)進(jìn)行晶格常數(shù)的優(yōu)化
已經(jīng)有4人回復(fù)
- 【求助】CASTEP計(jì)算彈性常數(shù)得到應(yīng)力應(yīng)變曲線
已經(jīng)有17人回復(fù)
- 【求助】三方晶系彈性常數(shù)計(jì)算求助
已經(jīng)有4人回復(fù)
- 【求助成功】晶格優(yōu)化的結(jié)果很詭異
已經(jīng)有13人回復(fù)
- 【求助】如何從Murnaghan方程的結(jié)果中得到優(yōu)化的晶格參數(shù)？
已經(jīng)有11人回復(fù)
- 【求助】castep中的Optical Permittivity和DC Permittivity
已經(jīng)有3人回復(fù)
- 【求助】CASTEP計(jì)算介電常數(shù)公式中的一個(gè)疑問
已經(jīng)有8人回復(fù)
- 【求助】CASTEP服務(wù)器彈性常數(shù)的計(jì)算
已經(jīng)有18人回復(fù)
【答案】應(yīng)助回帖
|
對(duì)于樓主的第二個(gè)問題：為什么增加截?cái)嗄軙?huì)不收斂，我理解的是因?yàn)樵黾咏財(cái)嗄軐?duì)應(yīng)的計(jì)算精度就會(huì)增加，就相當(dāng)于給結(jié)構(gòu)一個(gè)更小的偏差標(biāo)準(zhǔn)，所以才會(huì)不收斂，我最近在幾何優(yōu)化時(shí)也出現(xiàn)了這個(gè)問題，想問一下你有沒有在不改變精度的情況下解決了這個(gè)問題？ “BFGS: Warning - looks like this system is as converged as possible. Maybe your geometry convergence tolerances are too tight?” |
7樓2015-03-20 17:15:33
lhtlx95
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 4024
- 帖子: 114
- 在線: 341小時(shí)
- 蟲號(hào): 1622309
- 注冊(cè): 2012-02-17
- 性別: GG
- 專業(yè): 固體力學(xué)
2樓2013-11-12 10:13:49
3樓2013-11-13 09:46:41
4樓2013-11-13 09:48:15
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 22408 359分調(diào)劑 +4 | Qshers 2026-03-27 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 312，生物學(xué)求調(diào)劑 +3 | 小譯同學(xué)abc 2026-03-28 | 3/150 |
|
|
|
[考研] 0856求調(diào)劑 +11 | zhn03 2026-03-25 | 12/600 |
|
|
|
[考研] 085600 286分 材料求調(diào)劑 +7 | 麻辣魷魚 2026-03-27 | 8/400 |
|
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +7 | hanamiko 2026-03-27 | 7/350 |
|
|
|
[考研] 291求調(diào)劑 +6 | HanBeiNingZC 2026-03-24 | 6/300 |
|
|
|
[考研] 材料求調(diào)劑 一志愿哈工大總分298分，前三科223分 +5 | dongfang59 2026-03-27 | 5/250 |
|
|
|
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4 | Ripcord上岸 2026-03-27 | 4/200 |
|
|
|
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑 +9 | vv迷 2026-03-22 | 9/450 |
|
|
|
[考研] 0856調(diào)劑 +5 | 求求讓我有書讀?/a> 2026-03-26 | 6/300 |
|
|
|
[考研] 0703化學(xué)338求調(diào)劑！ +6 | Zuhui0306 2026-03-26 | 7/350 |
|
|
|
[考研] 284求調(diào)劑 +11 | junqihahaha 2026-03-26 | 12/600 |
|
|
|
[考研] 317求調(diào)劑 +7 | 蛋黃咸肉粽 2026-03-26 | 7/350 |
|
|
|
[考研] 求調(diào)劑 一志愿 本科 北科大 化學(xué) 343 +6 | 13831862839 2026-03-24 | 7/350 |
|
|
|
[考研] 321求調(diào)劑 +6 | Ymlll 2026-03-24 | 6/300 |
|
|
|
[考研] 329求調(diào)劑 +5 | 1() 2026-03-22 | 5/250 |
|
|
|
[考研] 一志愿天津大學(xué)339材料與化工求調(diào)劑 +3 | 江往賣魚 2026-03-26 | 3/150 |
|
|
|
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研經(jīng)歷 +6 | cccchenso 2026-03-23 | 6/300 |
|
|
|
[考研] 26考研-291分-廈門大學(xué)（085601）-柔性電子學(xué)院材料工程專業(yè)求調(diào)劑 +3 | min3 2026-03-24 | 4/200 |
|
|
|
[考研] 318求調(diào)劑 +3 | plum李子 2026-03-23 | 3/150 |
|