| 查看: 3939 | 回复: 8 | |||
justinzhouke新虫 (初入文坛)
|
[交流]
PWscf(QE)的结构优化问题 已有3人参与
|
|
我在用vc-relax做结构优化的时候，还没有达到应力收敛条件的时候就退出了 The maximum number of steps has been reached. 这个怎么处理啊。 求大神帮忙。 ............ Forces acting on atoms (Ry/au): atom 1 type 1 force = 0.00000000 0.00000000 0.00000000 atom 2 type 2 force = 0.00000000 0.00000000 0.00000000 Total force = 0.000000 Total SCF correction = 0.000000 entering subroutine stress ... total stress (Ry/bohr**3) (kbar) P= -0.02 -0.00000013 -0.00000000 0.00000000 -0.02 -0.00 0.00 0.00000000 -0.00000013 -0.00000000 0.00 -0.02 -0.00 0.00000000 -0.00000000 -0.00000013 0.00 -0.00 -0.02 The maximum number of steps has been reached. End of BFGS Geometry Optimization Begin final coordinates new unit-cell volume = 337.13304 a.u.^3 ( 49.95796 Ang^3 ) CELL_PARAMETERS (alat= 11.69699086) 0.000000000 0.472261950 0.472261950 0.472261950 -0.000000000 0.472261950 0.472261950 0.472261950 -0.000000000 ATOMIC_POSITIONS (crystal) In -0.000000000 -0.000000000 0.000000000 As 0.250000000 0.250000000 0.250000000 End final coordinates ............ 输入文件如下： &control calculation = 'vc-relax', restart_mode= 'from_scratch', pseudo_dir = './', outdir = './tmp', prefix='inas', wf_collect = .true. tstress = .true. tprnfor = .true. etot_conv_thr = 1.0D-4 forc_conv_thr = 1.0D-3 nstep = 20000 / &system ibrav= 0, celldm(1)=11.6969908603; nat= 2, ntyp= 2, ecutwfc =48, / &electrons electron_maxstep=2500 mixing_beta = 0.7 conv_thr = 5.0d-12 / &IONS ion_dynamics = 'bfgs' / &CELL cell_dynamics = 'bfgs' press=0.0 press_conv_thr = 5.0D-4 / ATOMIC_SPECIES In 111.818 In.pz-bhs.UPF As 74.920 As.pz-bhs.UPF ATOMIC_POSITIONS crystal In 0.00 0.00 0.00 As 0.25 0.25 0.25 K_POINTS automatic 9 9 9 0 0 0 CELL_PARAMETERS alat 0.000000000 0.5 0.5 0.5 -0.000000000 0.5 0.5 0.5 0.000000000 |
回复此楼» 猜你喜欢
上海工程技术大学【激光智能制造】课题组招收硕士
已经有6人回复
-
带资进组求博导收留
已经有11人回复
-
自荐读博
已经有5人回复
-
求个博导看看
已经有16人回复
-
上海工程技术大学张培磊教授团队招收博士生
已经有4人回复
-
求助院士们，这个如何合成呀
已经有4人回复
-
临港实验室与上科大联培博士招生1名
已经有9人回复
-
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？
已经有6人回复
-
最近几年招的学生写论文不引自己组发的文章
已经有11人回复
-
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- pwscf结构优化报错问题求助
已经有0人回复
- 第一次运用QE计算貌似出现问题，求解答
已经有12人回复
- QE（pwscf）计算计算声子谱问题
已经有1人回复
- Error in routine read_namelists (1): reading namelist ions 求助
已经有3人回复
- quantum-espresso（原先叫pwscf）安装问题
已经有8人回复
- QE_GPU加速的使用范围
已经有1人回复
- QE PWscf pw-gpu 计算速度 终止条件 收敛速度
已经有0人回复
- Quantum-ESPRESSO 5.0.3，Windows 32-bit，串行版
已经有55人回复
- PWscf终端运行报错，菜鸟求助
已经有4人回复
- 有关QE(Pwscf） PHonon编译 计算的一些基本问题
已经有10人回复
- QE(Pwscf) 非共线磁性计算
已经有0人回复
- pwscf-5.0得到的声子谱很乱是怎么回事？
已经有8人回复
- Problem - nscf - c_bands: eigenvalues not converged
已经有20人回复
- pwscf中average.x的使用
已经有0人回复
- 紧急求助：运行pwscf出现对称性操作警告
已经有43人回复
- 求专家解惑：pwscf并行计算cpu分配问题
已经有3人回复
- pwscf 安装测试出问题。。求救
已经有2人回复
- pwscf画能带图求助
已经有6人回复
- pwscf计算integrated charge出错及不收敛
已经有13人回复
liqizuiyang
木虫 (著名写手)
- 应助: 59 (初中生)
- 金币: 6859.5
- 散金: 261
- 红花: 94
- 帖子: 2093
- 在线: 708.9小时
- 虫号: 2366686
- 注册: 2013-03-21
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
0.5巴的压强……这个值低了，难怪不收敛。
2楼2015-08-28 09:10:56
justinzhouke
新虫 (初入文坛)
- 应助: 1 (幼儿园)
- 金币: 3.5
- 红花: 1
- 帖子: 24
- 在线: 6.5小时
- 虫号: 3513863
- 注册: 2014-11-02
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
|
你好， 我把press_conv_thr 改为1.0D-4和1.0D-5后，还是出现 The maximum number of steps has been reached. 其中1.0D-5的部分输出结果为： Forces acting on atoms (Ry/au): atom 1 type 1 force = 0.00000000 0.00000000 0.00000000 atom 2 type 2 force = 0.00000000 0.00000000 0.00000000 Total force = 0.000000 Total SCF correction = 0.000000 entering subroutine stress ... total stress (Ry/bohr**3) (kbar) P= -0.02 -0.00000013 -0.00000000 0.00000000 -0.02 -0.00 0.00 0.00000000 -0.00000013 -0.00000000 0.00 -0.02 -0.00 0.00000000 -0.00000000 -0.00000013 0.00 -0.00 -0.02 The maximum number of steps has been reached. |
3楼2015-08-28 10:47:04
liqizuiyang
木虫 (著名写手)
- 应助: 59 (初中生)
- 金币: 6859.5
- 散金: 261
- 红花: 94
- 帖子: 2093
- 在线: 708.9小时
- 虫号: 2366686
- 注册: 2013-03-21
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2015-08-28 11:18:43
justinzhouke
新虫 (初入文坛)
- 应助: 1 (幼儿园)
- 金币: 3.5
- 红花: 1
- 帖子: 24
- 在线: 6.5小时
- 虫号: 3513863
- 注册: 2014-11-02
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
|
我改为0.1后，优化后的结果为： Total force = 0.000000 Total SCF correction = 0.000000 entering subroutine stress ... total stress (Ry/bohr**3) (kbar) P= -0.05 -0.00000031 -0.00000000 0.00000000 -0.05 -0.00 0.00 -0.00000000 -0.00000031 0.00000000 -0.00 -0.05 0.00 0.00000000 -0.00000000 -0.00000031 0.00 -0.00 -0.05 bfgs converged in 6 scf cycles and 5 bfgs steps (criteria: energy < 0.10E-03, force < 0.10E-02, cell < 0.10E+00) End of BFGS Geometry Optimization Final enthalpy = -16.9554247603 Ry Begin final coordinates new unit-cell volume = 337.14584 a.u.^3 ( 49.95986 Ang^3 ) CELL_PARAMETERS (alat= 11.69699086) 0.000000000 0.472267929 0.472267929 0.472267929 -0.000000000 0.472267929 0.472267929 0.472267929 -0.000000000 ATOMIC_POSITIONS (crystal) In 0.000000000 0.000000000 -0.000000000 As 0.250000000 0.250000000 0.250000000 End final coordinates 没有出现 The maximum number of steps has been reached.可以没有达到应力为0. 当我改为0.01后还是出现 没有出现 The maximum number of steps has been reached. |
5楼2015-08-28 11:51:53
liqizuiyang
木虫 (著名写手)
- 应助: 59 (初中生)
- 金币: 6859.5
- 散金: 261
- 红花: 94
- 帖子: 2093
- 在线: 708.9小时
- 虫号: 2366686
- 注册: 2013-03-21
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
6楼2015-08-28 13:44:54
justinzhouke
新虫 (初入文坛)
- 应助: 1 (幼儿园)
- 金币: 3.5
- 红花: 1
- 帖子: 24
- 在线: 6.5小时
- 虫号: 3513863
- 注册: 2014-11-02
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
7楼2015-08-28 15:03:46
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
楼主我问个问题，像力 Forces acting on atoms (Ry/au): atom 1 type 1 force = 0.00000000 0.00000000 0.00000000 atom 2 type 2 force = 0.00000000 0.00000000 0.00000000 这怎么解释，为啥是三个值，还有怎么检查各个原子上的力，下面是我计算relax结果的力这个具体是什么意思啊 3个值 Forces acting on atoms (au): Ti -0.602690E-03 0.299794E-02 -0.103070E+00 Ti 0.114534E-02 -0.280239E-02 -0.235120E-01 Ti 0.154898E-02 0.768104E-02 0.230297E-01 Ti -0.305820E-02 -0.297500E-02 0.888280E-01 |
8楼2015-08-29 15:00:35
9楼2015-09-08 07:15:39