| 查看: 622 | 回复: 1 | |||
[交流]
【求助】晶格优化崩溃?
|
|
各位朋友。本人驰豫原子位置的时候,坐标一直变大?偏离实验数据跟文献上的计算数据很大,麻烦大家给我看看,我的输入文件有什么不妥的地方,万分感激。。。 以下是我的输入文件内容: &control calculation='vc-relax', restart_mode='from_scratch', prefix='Ba' tstress = .true. tprnfor = .true. pseudo_dir = './', outdir='./' nstep = 105 , etot_conv_thr = 1.0E-6 , forc_conv_thr = 1.0D-5 , dt = 50 , / &system ibrav = 0, nat = 5, ntyp = 3, ecutwfc = 40.0 , ecutrho = 400.0 , nosym = .false. , nbnd = 80, occupations = 'smearing' , degauss = 0.05, smearing = 'marzari-vanderbilt' , nspin = 1 , / &electrons mixing_beta = 0.7 conv_thr = 1.0d-8 / &IONS / &CELL cell_dynamics = 'bfgs' press = 1000 / 。。。 。。。 。。。 K_POINTS automatic 8 8 8 0 0 0 |
» 猜你喜欢
面上提前没消息,有中的吗
已经有14人回复
这种情况还有戏吗
已经有12人回复
微信指数没变化,科研之友没阅读
已经有19人回复
准备明年的基金了
已经有7人回复
面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇!
已经有24人回复
2026年国自然面上资助率
已经有21人回复
基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊
已经有10人回复
2027年申博
已经有3人回复
一个有机合成实验室都需要哪些设备?
已经有7人回复
HXDI做水性聚氨酯乳液,是不是特别容易出渣
已经有3人回复
» 本主题相关价值贴推荐,对您同样有帮助:
关于优化晶格常数
已经有4人回复
晶格常数优化
已经有7人回复
【求助】请高人指点一下掺氮TiO2的制备
已经有19人回复
【求助】CASTEP计算弹性常数得到应力应变曲线
已经有17人回复
【求助】vasp优化晶格参数a
已经有13人回复
【求助】简单问题----关于vasp晶格优化
已经有7人回复
【求助成功】晶格优化的结果很诡异
已经有13人回复
【求助成功】晶格优化中,murn.f的使用方法
已经有7人回复
【求助】侯老师的vasp入门的晶格优化部分没看懂,问一下
已经有25人回复
【原创】晶格优化中ISIF取值的探究(六方MgO为例)
已经有9人回复
» 抢金币啦!回帖就可以得到:
2027年博士生(硕博连读/直博)吉林大学-苏州国家实验室光电半导体智能设计与应用团队
+1/474
山东征女友,坐标济南
+1/184
坐标广州,征女友
+2/138
南京大学医工交叉方向招聘脂质纳米材料与RNA递送方向博士后1-2名
+1/85
★★★中山大学(深圳)院士团队招收27届推免硕士生(低维材料与器件方向)
+1/79
采用本质安全型氢气流量控制的化学合成 - Alicat ISMC本安防爆质量流量控制器
+2/78
Postdoctoral fellow in Medicinal Chemistry, McGill University
+1/71
2026年University of South Carolina继续招收博士生
+1/64
连续氢化反应中流量控制器、背压控制器的应用 | Alicat 气液压力测控全局解决方案
+2/64
中科院东莞材料所依托中科大招收2027届推免生
+1/35
西湖大学拓扑光学、非厄米光学、太赫兹方向博士后招聘(长期有效)
+5/30
招收博士后(神经退行性疾病、认知障碍疾病)
+1/27
西湖大学物理学、光学、电子信息方向博士生招生(长期有效)
+5/25
诚征女友( 西安 )
+1/18
西湖大学拓扑光学、非厄米光学、太赫兹方向博士后招聘(长期有效)
+5/15
北理工集成电路杰青团队 | 诚招科助理
+1/10
(Scopus期刊)Biofunctional Materials《生物功能性材料》招募青年编委及征稿
+1/4
国自然发榜在即:此刻的你,是什么心情?
+1/2
有没有管理类博士招生
+1/1
这些Scopus收录的英文国际期刊正在征稿(快速审稿)2026年底前免APC!
+1/1
输出
|
以下是输出: Forces acting on atoms (Ry/au): atom 1 type 1 force = 0.00000000 0.00000000 0.00000000 atom 2 type 2 force = -0.00076425 0.00076425 0.00000000 atom 3 type 2 force = 0.00076425 -0.00076425 0.00000000 atom 4 type 3 force = -0.00079442 -0.00079442 0.43102097 atom 5 type 3 force = 0.00079442 0.00079442 -0.43102097 Total force = 0.000000 Total SCF correction = 0.000000 entering subroutine stress ... total stress (Ry/bohr**3) (kbar) P= 989.12 0.00670743 0.00000026 -0.00001451 986.70 0.04 -2.13 0.00000026 0.00670743 -0.00001451 0.04 986.70 -2.13 -0.00001451 -0.00001451 0.00675680 -2.13 -2.13 993.96 number of scf cycles = 11 number of bfgs steps = 10 enthalpy old = -666.0630713336 Ry enthalpy new = -666.0638799665 Ry CASE: enthalpy_new < enthalpy_old new trust radius = 0.0033684758 bohr new conv_thr = 0.0000000010 Ry new unit-cell volume = 379.17935 a.u.^3 ( 56.18858 Ang^3 ) CELL_PARAMETERS (alat= 3.96250000) 1.627610836 -0.000000404 0.000086123 -0.000000404 1.627610836 0.000086123 0.813717998 0.813717998 2.300657444 |
2楼2010-12-11 19:27:24










回复此楼