 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2988 | 回復(fù): 7 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
lhtlx95新蟲 (小有名氣)
|
[求助]
用castep計算單斜晶系彈性常數(shù)的問題 已有1人參與
|
|
|
最近利用castep計算某一層狀硅酸鹽（單斜晶系）的彈性常數(shù)，計算過程中遇到一些問題，想請大家指教！ 1、castep中彈性常數(shù)的算法是在給定的一系列的應(yīng)變模式下進(jìn)行結(jié)構(gòu)優(yōu)化，而且是只優(yōu)化原子位置，而不優(yōu)化晶格常數(shù)。然后計算應(yīng)力，再利用應(yīng)力應(yīng)變進(jìn)行擬合，得到彈性常數(shù)。（不知道我理解的對不對，請大家指教）。而對于單斜晶系，晶角beta不是90度，所以計算出的彈性常數(shù)應(yīng)該和坐標(biāo)系的選擇有關(guān)系。我測試選擇了兩種坐標(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)行彈性模量的計算，但我發(fā)現(xiàn)在同一應(yīng)變模式，比如Z方向加應(yīng)變0.003，兩種坐標(biāo)系下加應(yī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è)置沒有對! 但進(jìn)一步看，不同坐標(biāo)系下，都給Z方向加0.003的應(yīng)變后，計算得到的應(yīng)力又不相同。 所以很困惑，castep中應(yīng)變是按照什么原理加的？？？？ 2、計算彈性模量過程中如果截斷能選600eV，收斂的比較好！如下： +-----------+-----------------+-----------------+------------+-----+ <-- 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ù)不變，只是把截斷能提高到1000eV，優(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. 請問大家這是為什么？好像彈性模量對截斷能和K點(diǎn)的要求比較高，為什么提高截斷能反而不收斂？ 3、我計算的彈性模量結(jié)果（先按三斜晶系計算）： 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、我看到論壇里有蟲友說可以用提腳本文件的方法計算彈性模量，不知利用腳本文件和直接計算有什么不同？ |
回復(fù)此樓» 猜你喜歡
303求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿武理085500機(jī)械專業(yè)總分300求調(diào)劑
已經(jīng)有7人回復(fù)
-
材料學(xué)求調(diào)劑
已經(jīng)有5人回復(fù)
-
考研調(diào)劑
已經(jīng)有4人回復(fù)
-
281求調(diào)劑
已經(jīng)有4人回復(fù)
-
0805 316求調(diào)劑
已經(jīng)有6人回復(fù)
-
085601求調(diào)劑總分293英一數(shù)二
已經(jīng)有3人回復(fù)
-
08工學(xué)調(diào)劑
已經(jīng)有17人回復(fù)
-
340求調(diào)劑
已經(jīng)有4人回復(fù)
-
311求調(diào)劑
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 用CASTEP計算彈性常數(shù)出錯的原因
已經(jīng)有5人回復(fù)
- E-V曲線擬合后怎么找到能量最低點(diǎn)對應(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īng)有16人回復(fù)
- 晶格常數(shù)變化的討論c/a 還是 b/a?
已經(jīng)有9人回復(fù)
- CASTEP 計算介電常數(shù) 數(shù)值的意義
已經(jīng)有6人回復(fù)
- 為什么采用CASTEP計算CaTiO3的光學(xué)性質(zhì)其中的靜態(tài)介電常數(shù)與實(shí)驗(yàn)值相差很大
已經(jīng)有14人回復(fù)
- CASTEP計算彈性常數(shù)C11，C12，C44嚴(yán)重偏大
已經(jīng)有13人回復(fù)
- 計算不同晶系的彈性常數(shù)
已經(jīng)有260人回復(fù)
- 【求助】castep計算彈性常數(shù)問題
已經(jīng)有14人回復(fù)
- 【求助】請教手動進(jìn)行晶格常數(shù)的優(yōu)化
已經(jīng)有4人回復(fù)
- 【求助】CASTEP計算彈性常數(shù)得到應(yīng)力應(yīng)變曲線
已經(jīng)有17人回復(fù)
- 【求助】三方晶系彈性常數(shù)計算求助
已經(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計算介電常數(shù)公式中的一個疑問
已經(jīng)有8人回復(fù)
- 【求助】CASTEP服務(wù)器彈性常數(shù)的計算
已經(jīng)有18人回復(fù)
4樓2013-11-13 09:48:15
lhtlx95
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 4024
- 帖子: 114
- 在線: 341小時
- 蟲號: 1622309
- 注冊: 2012-02-17
- 性別: GG
- 專業(yè): 固體力學(xué)
2樓2013-11-12 10:13:49
3樓2013-11-13 09:46:41
lhtlx95
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 4024
- 帖子: 114
- 在線: 341小時
- 蟲號: 1622309
- 注冊: 2012-02-17
- 性別: GG
- 專業(yè): 固體力學(xué)
5樓2013-11-13 15:50:37
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0805 316求調(diào)劑 +6 | 大雪深藏 2026-03-18 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 086003食品工程求調(diào)劑 +5 | 淼淼111 2026-03-24 | 5/250 |
|
|
|
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +6 | 小羊睡著了? 2026-03-23 | 9/450 |
|
|
|
[考研]
 085602 289分求調(diào)劑
+5
|
WWW西西弗斯 2026-03-24 | 5/250 |
|
|
|
[考研] 招08考數(shù)學(xué) +7 | laoshidan 2026-03-20 | 16/800 |
|
|
|
[考研] 081700 調(diào)劑 267分 +9 | 迷人的哈哈 2026-03-23 | 9/450 |
|
|
[考研] 335求調(diào)劑 +4 | yuyu宇 2026-03-23 | 5/250 |
|
|
|
[考研] 一志愿武理材料工程348求調(diào)劑 +6 | ￣^￣゜汗 2026-03-19 | 9/450 |
|
|
|
[考研] 311求調(diào)劑 +3 | 26研0 2026-03-20 | 3/150 |
|
|
|
[考研] 生物學(xué)調(diào)劑 +5 | Surekei 2026-03-21 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +7 | Auroracx 2026-03-22 | 7/350 |
|
|
|
[考研]
一志愿東華大學(xué)控制學(xué)碩320求調(diào)劑
+3
|
Grand777 2026-03-21 | 3/150 |
|
|
|
[考研] 297求調(diào)劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|
|
[考研] 一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3 | 晨昏線與星海 2026-03-18 | 3/150 |
|
|
|
[考研] 一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +4 | 晨昏線與星海 2026-03-19 | 4/200 |
|
|
|
[考研] 290求調(diào)劑 +7 | ^O^乜 2026-03-19 | 7/350 |
|
|
|
[考研] 261求B區(qū)調(diào)劑，科研經(jīng)歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +3 | eation27 2026-03-20 | 3/150 |
|
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊(duì) 2026-03-20 | 3/150 |
|
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|