 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2987 | 回復(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)系。我測試選擇了兩種坐標(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)變后，晶格長度完全相同�。�！ 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ù)此樓» 猜你喜歡
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)價(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ù)
3樓2013-11-13 09:46:41
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
4樓2013-11-13 09:48:15
lhtlx95
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 4024
- 帖子: 114
- 在線: 341小時(shí)
- 蟲號(hào): 1622309
- 注冊(cè): 2012-02-17
- 性別: GG
- 專業(yè): 固體力學(xué)
5樓2013-11-13 15:50:37
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 300分，材料，求調(diào)劑，英一數(shù)二 +5 | 超贊的 2026-03-24 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 085600材料與化工調(diào)劑 +8 | A-哆啦Z夢 2026-03-23 | 13/650 |
|
|
|
[考研] 07化學(xué)280分求調(diào)劑 +6 | 722865 2026-03-23 | 6/300 |
|
|
|
[考研] 材料專碩331求調(diào)劑 +4 | 鮮當(dāng)牛 2026-03-24 | 4/200 |
|
|
|
[考博] 申博26年 +4 | 八6八68 2026-03-19 | 4/200 |
|
|
|
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +4 | 曼111 2026-03-24 | 5/250 |
|
|
|
[考研] 一志愿國科過程所081700，274求調(diào)劑 +3 | 三水研0水立方 2026-03-23 | 3/150 |
|
|
[考研] 361求調(diào)劑 +3 | Glack 2026-03-22 | 3/150 |
|
|
|
[考研] 生物學(xué)一志愿985，分?jǐn)?shù)349求調(diào)劑 +6 | zxts12 2026-03-21 | 9/450 |
|
|
|
[考研] 336化工調(diào)劑 +4 | 王大坦1 2026-03-23 | 5/250 |
|
|
|
[考研] 263求調(diào)劑 +6 | yqdszhdap－ 2026-03-22 | 9/450 |
|
|
|
[考研] 276求調(diào)劑 +3 | YNRYG 2026-03-21 | 4/200 |
|
|
|
[考研] 323求調(diào)劑 +6 | 洼小桶 2026-03-18 | 6/300 |
|
|
|
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|
|
|
[考研] 303求調(diào)劑 +5 | 安憶靈 2026-03-22 | 6/300 |
|
|
|
[考研] 297求調(diào)劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | .m.. 2026-03-21 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +3 | 白QF 2026-03-21 | 3/150 |
|
|
|
[考研] 330求調(diào)劑 +4 | 小材化本科 2026-03-18 | 4/200 |
|
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|