 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 4532 | 回復(fù): 11 | ||
[求助]
幾何結(jié)構(gòu)優(yōu)化時(shí)能量震蕩，如何處理？ 已有1人參與
|
||
|
現(xiàn)在做的摻雜體系，12.5%，25%，50%，其中12.5%，50%都能收斂，就是25%不能收斂，什么原因呢？ 不太可能是結(jié)構(gòu)不合理吧，畢竟兩邊濃度都成功了。請(qǐng)大家指教。 像這種能量震蕩的該如何調(diào)試呢？ castep文件：幾何優(yōu)化只有一個(gè)yes，3個(gè)no。 +------------+-------------+-------------+-----------------+ <-- min BFGS | Step | lambda | F.delta | enthalpy | <-- min BFGS +------------+-------------+-------------+-----------------+ <-- min BFGS | previous | 0.000000 | 0.002781 | -29834.670729 | <-- min BFGS | trial step | 1.000000 | 0.001926 | -29834.776186 | <-- min BFGS | line step | 3.251631 | -0.003682 | -29837.973701 | <-- min BFGS | quad step | 1.773369 | -0.002702 | -29838.006664 | <-- min BFGS +------------+-------------+-------------+-----------------+ <-- min BFGS BFGS: finished iteration 100 with enthalpy= -2.98380067E+004 eV +-----------+-----------------+-----------------+------------+-----+ <-- BFGS | Parameter | value | tolerance | units | OK? | <-- BFGS +-----------+-----------------+-----------------+------------+-----+ <-- BFGS | dE/ion | 8.339837E-002 | 5.000000E-005 | eV | No | <-- BFGS | |F|max | 1.288961E-001 | 1.000000E-001 | eV/A | No | <-- BFGS | |dR|max | 6.082554E-004 | 5.000000E-003 | A | Yes | <-- BFGS | Smax | 4.065663E+000 | 2.000000E-001 | GPa | No | <-- BFGS +-----------+-----------------+-----------------+------------+-----+ <-- BFGS BFGS : WARNING - Geometry optimization failed to converge after 100 steps ================================================================================ BFGS : Final Configuration: ================================================================================ ------------------------------- Unit Cell ------------------------------- Real Lattice(A) Reciprocal Lattice(1/A) 7.8142248 0.0000000 0.0000000 0.8040702 0.0000000 0.0000000 0.0000000 7.8142248 0.0000000 0.0000000 0.8040702 0.0000000 0.0000000 0.0000000 8.1112610 0.0000000 0.0000000 0.7746250 Lattice parameters(A) Cell Angles a = 7.814225 alpha = 90.000000 b = 7.814225 beta = 90.000000 c = 8.111261 gamma = 90.000000 Current cell volume = 495.290711 A**3 ------------------------------- Cell Contents ------------------------------- xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx x Element Atom Fractional coordinates of atoms x x Number u v w x x----------------------------------------------------------x x O 1 0.000000 0.000000 0.043401 x x O 2 0.000000 0.000000 0.516500 x x O 3 0.000000 0.243773 0.271102 x x O 4 0.000000 -0.248546 0.775364 x x O 5 0.000000 -0.243773 0.271102 x x O 6 0.000000 0.248546 0.775364 x x O 7 -0.243773 0.000000 0.271102 x x O 8 0.248546 0.000000 0.775364 x x O 9 0.243773 0.000000 0.271102 x x O 10 -0.248546 0.000000 0.775364 x x O 11 0.500000 0.000000 0.021770 x x O 12 0.500000 0.000000 0.520766 x x O 13 0.000000 0.500000 0.021770 x x O 14 0.000000 0.500000 0.520766 x x O 15 0.500000 0.250430 0.283155 x x O 16 0.500000 -0.257225 0.761425 x x O 17 -0.500000 -0.250430 0.283155 x x O 18 -0.500000 0.257225 0.761425 x x O 19 -0.250430 0.500000 0.283155 x x O 20 0.257225 0.500000 0.761425 x x O 21 0.250430 -0.500000 0.283155 x x O 22 -0.257225 -0.500000 0.761425 x x O 23 0.500000 0.500000 0.024175 x x O 24 0.500000 0.500000 0.533453 x x Na 1 -0.250144 0.250144 0.498527 x x Na 2 0.250144 -0.250144 0.498527 x x Na 3 -0.250144 -0.250144 0.498527 x x Na 4 0.250144 0.250144 0.498527 x x K 1 -0.248651 -0.248651 0.003625 x x K 2 0.248651 0.248651 0.003625 x x K 3 0.248651 -0.248651 0.003625 x x K 4 -0.248651 0.248651 0.003625 x x Fe 1 0.000000 0.000000 0.255229 x x Fe 2 0.500000 0.500000 0.744987 x x Nb 1 0.000000 0.000000 0.741251 x x Nb 2 0.500000 0.000000 0.252241 x x Nb 3 0.500000 0.000000 0.748703 x x Nb 4 0.000000 0.500000 0.252241 x x Nb 5 0.000000 0.500000 0.748703 x x Nb 6 0.500000 0.500000 0.249257 x xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx BFGS : Final Enthalpy = -2.98380067E+004 eV BFGS : Final <frequency> = 7075.84270 cm-1 BFGS : Final bulk modulus unchanged from initial value ********************************************************** *** There were at least 1 warnings during this run *** *** => please check the whole of this file carefully! *** ********************************************************** ***************************** Symmetrised Forces ***************************** * * * Cartesian components (eV/A) * * -------------------------------------------------------------------------- * * x y z * * * * O 1 0.00000 0.00000 -0.01499 * * O 2 0.00000 0.00000 0.05276 * * O 3 0.00000 0.05652 -0.01386 * * O 4 0.00000 -0.03243 -0.01695 * * O 5 0.00000 -0.05652 -0.01386 * * O 6 0.00000 0.03243 -0.01695 * * O 7 -0.05652 0.00000 -0.01386 * * O 8 0.03243 0.00000 -0.01695 * * O 9 0.05652 0.00000 -0.01386 * * O 10 -0.03243 0.00000 -0.01695 * * O 11 0.00000 0.00000 0.02100 * * O 12 0.00000 0.00000 -0.00609 * * O 13 0.00000 0.00000 0.02100 * * O 14 0.00000 0.00000 -0.00609 * * O 15 0.00000 -0.00585 -0.03748 * * O 16 0.00000 0.08450 -0.01428 * * O 17 0.00000 0.00585 -0.03748 * * O 18 0.00000 -0.08450 -0.01428 * * O 19 0.00585 0.00000 -0.03748 * * O 20 -0.08450 0.00000 -0.01428 * * O 21 -0.00585 0.00000 -0.03748 * * O 22 0.08450 0.00000 -0.01428 * * O 23 0.00000 0.00000 0.00823 * * O 24 0.00000 0.00000 0.03740 * * Na 1 0.06417 -0.06417 0.09154 * * Na 2 -0.06417 0.06417 0.09154 * * Na 3 0.06417 0.06417 0.09154 * * Na 4 -0.06417 -0.06417 0.09154 * * K 1 -0.02515 -0.02515 -0.01297 * * K 2 0.02515 0.02515 -0.01297 * * K 3 0.02515 -0.02515 -0.01297 * * K 4 -0.02515 0.02515 -0.01297 * * Fe 1 0.00000 0.00000 -0.08922 * * Fe 2 0.00000 0.00000 -0.07155 * * Nb 1 0.00000 0.00000 0.04927 * * Nb 2 0.00000 0.00000 0.01781 * * Nb 3 0.00000 0.00000 -0.02764 * * Nb 4 0.00000 0.00000 0.01781 * * Nb 5 0.00000 0.00000 -0.02764 * * Nb 6 0.00000 0.00000 0.03392 * * * ****************************************************************************** *********** Symmetrised Stress Tensor *********** * * * Cartesian components (GPa) * * --------------------------------------------- * * x y z * * * * x -4.065663 0.000000 0.000000 * * y 0.000000 -4.065663 0.000000 * * z 0.000000 0.000000 0.306267 * * * * Pressure: 2.6084 * * * ************************************************* Pseudo atomic calculation performed for O 2s2 2p4 Converged in 23 iterations to a total energy of -428.2520 eV Pseudo atomic calculation performed for Na 2s2 2p6 3s1 Converged in 21 iterations to a total energy of -1298.1443 eV Pseudo atomic calculation performed for K 3s2 3p6 4s1 Converged in 23 iterations to a total energy of -777.4801 eV Pseudo atomic calculation performed for Fe 3d6 4s2 Converged in 32 iterations to a total energy of -849.3473 eV Pseudo atomic calculation performed for Nb 4s2 4p6 4d4 5s1 Converged in 21 iterations to a total energy of -1537.8413 eV Charge spilling parameter for spin component 1 = 0.28% Charge spilling parameter for spin component 2 = 0.28% |
回復(fù)此樓» 猜你喜歡
求調(diào)劑院校信息
已經(jīng)有4人回復(fù)
-
085600材料與化工306
已經(jīng)有4人回復(fù)
-
286求調(diào)劑
已經(jīng)有10人回復(fù)
-
328求調(diào)劑，英語六級(jí)551，有科研經(jīng)歷
已經(jīng)有9人回復(fù)
-
一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑
已經(jīng)有4人回復(fù)
-
286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考！
已經(jīng)有8人回復(fù)
-
資源與環(huán)境 調(diào)劑申請(qǐng)(333分)
已經(jīng)有5人回復(fù)
-
280求調(diào)劑
已經(jīng)有12人回復(fù)
-
269專碩求調(diào)劑
已經(jīng)有5人回復(fù)
-
材料學(xué)碩301分求調(diào)劑
已經(jīng)有7人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 2011年博士點(diǎn)基金聯(lián)合資助課題立項(xiàng)公布
已經(jīng)有3人回復(fù)
默默的小蟲
木蟲 (正式寫手)
- 應(yīng)助: 4 (幼兒園)
- 金幣: 7398.5
- 散金: 50
- 紅花: 1
- 帖子: 769
- 在線: 199.5小時(shí)
- 蟲號(hào): 3819752
- 注冊(cè): 2015-04-19
- 性別: GG
- 專業(yè): 無機(jī)材料化學(xué)
4樓2016-06-02 09:31:01
默默的小蟲
木蟲 (正式寫手)
- 應(yīng)助: 4 (幼兒園)
- 金幣: 7398.5
- 散金: 50
- 紅花: 1
- 帖子: 769
- 在線: 199.5小時(shí)
- 蟲號(hào): 3819752
- 注冊(cè): 2015-04-19
- 性別: GG
- 專業(yè): 無機(jī)材料化學(xué)
8樓2016-06-03 10:37:39
小雨熊
木蟲 (正式寫手)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 1513.5
- 散金: 120
- 紅花: 2
- 帖子: 546
- 在線: 318.9小時(shí)
- 蟲號(hào): 1431950
- 注冊(cè): 2011-10-08
- 專業(yè): 半導(dǎo)體材料
|
這個(gè)可能是因?yàn)槟憬Y(jié)構(gòu)不合理，原子偏離平衡位置太遠(yuǎn)了，所以一直在震蕩，你試一下手動(dòng)調(diào)低精度，或者，如果原子位置不好調(diào)節(jié)，那你試下改變abc的值，然后再優(yōu)化看看 發(fā)自小木蟲Android客戶端 |
9樓2016-06-03 14:55:14
10樓2016-06-05 10:36:23
2樓2016-05-31 22:08:08
3樓2016-06-01 20:18:18
小雨熊
木蟲 (正式寫手)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 1513.5
- 散金: 120
- 紅花: 2
- 帖子: 546
- 在線: 318.9小時(shí)
- 蟲號(hào): 1431950
- 注冊(cè): 2011-10-08
- 專業(yè): 半導(dǎo)體材料
5樓2016-06-02 15:19:17
|
感謝您的回復(fù)： 1.我就是覺得構(gòu)型不應(yīng)該有問題啊，因?yàn)橥瑯拥臉?gòu)型12.5%和50%兩種濃度都可以成功，唯獨(dú)中間濃度25%不行。如圖，就是Fe取代b位鈮。 2.您的意思把截?cái)嗄埽瑂mearing都調(diào)小一點(diǎn)嗎？事實(shí)上我就是用的coarse精度，還需要調(diào)小嗎？k點(diǎn)呢？ 3.設(shè)置： *********************** Exchange-Correlation Parameters *********************** using functional : Local Density Approximation Divergence correction : off DFT+D: Semi-empirical dispersion correction : off ************************* Pseudopotential Parameters ************************** pseudopotential representation : reciprocal space <beta|phi> representation : reciprocal space **************************** Basis Set Parameters ***************************** plane wave basis set cut-off : 280.0000 eV size of standard grid : 1.5000 size of fine gmax : 12.8590 1/A largest prime factor in FFT : 5 finite basis set correction : automatic number of sample energies : 3 sample spacing : 5.0000 eV **************************** Electronic Parameters **************************** number of electrons : 310.0 net charge of system : 0.000 net spin of system : 10.00 number of up spins : 160.0 number of down spins : 150.0 treating system as spin-polarized number of bands : 193 ********************* Electronic Minimization Parameters ********************** Method: Treating system as metallic with density mixing treatment of electrons, and number of SD steps : 1 and number of CG steps : 4 total energy / atom convergence tol. : 0.1000E-04 eV eigen-energy convergence tolerance : 0.1000E-05 eV max force / atom convergence tol. : ignored convergence tolerance window : 3 cycles max. number of SCF cycles : 100 number of fixed-spin iterations : 6 smearing scheme : Gaussian smearing width : 0.1000 eV Fermi energy convergence tolerance : 0.1000E-06 eV periodic dipole correction : NONE ************************** Density Mixing Parameters ************************** density-mixing scheme : Pulay max. length of mixing history : 20 charge density mixing amplitude : 0.5000 spin density mixing amplitude : 2.000 cut-off energy for mixing : 280.0 eV charge density mixing g-vector : 1.500 1/A spin density mixing g-vector : 1.500 1/A *********************** Population Analysis Parameters ************************ Population analysis with cutoff : 3.000 A ********************** Geometry Optimization Parameters *********************** optimization method : BFGS variable cell method : fixed basis quality max. number of steps : 100 estimated bulk modulus : 500.0 GPa estimated <frequency> : 1668. cm-1 geom line minimiser : on with line minimiser tolerance : 0.4000 with spin fully able to relax for all steps total energy convergence tolerance : 0.5000E-04 eV/atom max ionic |force| tolerance : 0.1000 eV/A max ionic |displacement| tolerance : 0.5000E-02 A max |stress component| tolerance : 0.2000 GPa convergence tolerance window : 2 steps backup results every : 5 steps 4.報(bào)錯(cuò)信息： 我試過兩種，一種就直接失敗了，不收斂： *Warning* max. SCF cycles performed but system has not reached the groundstate. Current total energy, E = -29833.17159871 eV Current free energy (E-TS) = -29833.31934436 eV (energies not corrected for finite basis set) NB est. 0K energy (E-0.5TS) = -29833.24547154 eV **************************************************************************** Warning: electronic minimisation did not converge when finding ground state. **************************************************************************** Writing model to Fe_25.check Error in geom_get_forces - electronic_minimisation of current_cell failed Current trace stack: geom_get_forces geom_BFGS geometry_optimise castep Error in geom_get_forces - electronic_minimisation of current_cell failed Current trace stack: geom_get_forces geom_BFGS geometry_optimise castep 還有一種是幾何優(yōu)化4項(xiàng)，3個(gè)no1個(gè)也yes： +-----------+-----------------+-----------------+------------+-----+ <-- BFGS | Parameter | value | tolerance | units | OK? | <-- BFGS +-----------+-----------------+-----------------+------------+-----+ <-- BFGS | dE/ion | 8.339837E-002 | 5.000000E-005 | eV | No | <-- BFGS | |F|max | 1.288961E-001 | 1.000000E-001 | eV/A | No | <-- BFGS | |dR|max | 6.082554E-004 | 5.000000E-003 | A | Yes | <-- BFGS | Smax | 4.065663E+000 | 2.000000E-001 | GPa | No | <-- BFGS +-----------+-----------------+-----------------+------------+-----+ <-- BFGS BFGS : WARNING - Geometry optimization failed to converge after 100 steps BFGS : Final Enthalpy = -2.98380067E+004 eV BFGS : Final <frequency> = 7075.84270 cm-1 BFGS : Final bulk modulus unchanged from initial value ********************************************************** *** There were at least 1 warnings during this run *** *** => please check the whole of this file carefully! *** ********************************************************** |
6樓2016-06-03 09:25:36
7樓2016-06-03 09:26:31
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 anita55 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料求調(diào)劑 +5 | @taotao 2026-03-21 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 0703化學(xué)調(diào)劑 ，六級(jí)已過，有科研經(jīng)歷 +14 | 曦熙兮 2026-03-15 | 14/700 |
|
|
|
[考研] 0805 316求調(diào)劑 +3 | 大雪深藏 2026-03-18 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +5 | 十三加油 2026-03-21 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +3 | 13341 2026-03-20 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | 白QF 2026-03-21 | 3/150 |
|
|
|
[考研] 一志愿西南交大，求調(diào)劑 +5 | 材化逐夢(mèng)人 2026-03-18 | 5/250 |
|
|
|
[考研] 294求調(diào)劑材料與化工專碩 +15 | 陌の森林 2026-03-18 | 15/750 |
|
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4 | 葵梓衛(wèi)隊(duì) 2026-03-18 | 6/300 |
|
|
|
[考研] 288求調(diào)劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 5/250 |
|
|
|
[考研] 353求調(diào)劑 +3 | 拉鉤不許變 2026-03-20 | 3/150 |
|
|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11 | Ymlll 2026-03-18 | 15/750 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +7 | 困于星晨 2026-03-17 | 9/450 |
|
|
|
[考研] 生物學(xué)調(diào)劑招人�。�！ +3 | 山海天嵐 2026-03-17 | 4/200 |
|
|
|
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進(jìn)步了0608 2026-03-16 | 6/300 |
|
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|