 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2994 | 回復(fù): 7 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
lhtlx95新蟲 (小有名氣)
|
[求助]
用castep計(jì)算單斜晶系彈性常數(shù)的問題 已有1人參與
|
|
|
最近利用castep計(jì)算某一層狀硅酸鹽（單斜晶系）的彈性常數(shù)，計(jì)算過程中遇到一些問題，想請大家指教！ 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ù)。（不知道我理解的對不對，請大家指教）。而對于單斜晶系，晶角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è)置沒有對! 但進(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. 請問大家這是為什么？好像彈性模量對截?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ù)此樓» 猜你喜歡
材料科學(xué)與工程求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿大連理工大學(xué)材料求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿北化085600材料專碩275|有文章專利｜求調(diào)劑
已經(jīng)有9人回復(fù)
-
材料求調(diào)劑 一志愿哈工大總分298分，前三科223分
已經(jīng)有8人回復(fù)
-
274求調(diào)劑
已經(jīng)有3人回復(fù)
-
合肥區(qū)域性重點(diǎn)一本招收調(diào)劑
已經(jīng)有4人回復(fù)
-
085600材料與化工調(diào)劑
已經(jīng)有10人回復(fù)
-
22408 359分調(diào)劑
已經(jīng)有7人回復(fù)
-
抱歉
已經(jīng)有4人回復(fù)
-
085602化工求調(diào)劑（331分）
已經(jīng)有8人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- 用CASTEP計(jì)算彈性常數(shù)出錯(cuò)的原因
已經(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ì)算
已經(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ù)
- 【求助】請教手動(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ù)
lhtlx95
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 4024
- 帖子: 114
- 在線: 341小時(shí)
- 蟲號(hào): 1622309
- 注冊: 2012-02-17
- 性別: GG
- 專業(yè): 固體力學(xué)
5樓2013-11-13 15:50:37
lhtlx95
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 4024
- 帖子: 114
- 在線: 341小時(shí)
- 蟲號(hào): 1622309
- 注冊: 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ā)表 | |
|---|---|---|---|---|
|
[考研] 279求調(diào)劑 +12 | j的立方 2026-03-29 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 一志愿鄭州大學(xué)，080500學(xué)碩，總分317分求調(diào)劑 +10 | 舉個(gè)栗子oi 2026-03-24 | 11/550 |
|
|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分 +12 | 辰沐5211314 2026-03-26 | 12/600 |
|
|
|
[考研] 322求調(diào)劑：一志愿湖南大學(xué) 材料與化工（085600），已過六級(jí)。 +9 | XX小鄧 2026-03-29 | 9/450 |
|
|
|
[考研] 286求調(diào)劑 +5 | Faune 2026-03-30 | 5/250 |
|
|
|
[考研] 化學(xué)0703 調(diào)劑 306分 一志愿211 +8 | 26要上岸 2026-03-28 | 8/400 |
|
|
|
[考研] 0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分 +7 | Shern—- 2026-03-28 | 7/350 |
|
|
|
[考研] 材料專碩調(diào)劑 +11 | 椰椰。 2026-03-29 | 11/550 |
|
|
|
[考研] 303求調(diào)劑 +7 | DLkz1314. 2026-03-30 | 7/350 |
|
|
|
[考研] 一志愿211，335分，0856，求調(diào)劑院校和導(dǎo)師 +7 | 傾____蕭 2026-03-27 | 8/400 |
|
|
|
[考研] 290求調(diào)劑 +3 | dfffsar 2026-03-29 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +4 | QiMing7 2026-03-25 | 5/250 |
|
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-27 | 7/350 |
|
|
|
[考研] 化學(xué)調(diào)劑 +4 | 愛吃番茄的旭 2026-03-24 | 5/250 |
|
|
|
[考研] 315調(diào)劑 +4 | 0860求調(diào)劑 2026-03-26 | 5/250 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
|
[考研] 環(huán)境專碩324分求調(diào)劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|
|
|
[考研] 求b區(qū)院校調(diào)劑 +4 | 周56 2026-03-24 | 5/250 |
|
|
|
[考研] 302求調(diào)劑 +4 | 錦衣衛(wèi)藤椒 2026-03-25 | 4/200 |
|
|
|
[考研] 344求調(diào)劑 +3 | desto 2026-03-24 | 3/150 |
|