计算化学公社

标题: DMol3几何优化，delta E参数是no，是否收敛 [打印本页]

作者
Author: Channing 时间: 2018-2-5 20:09
标题: DMol3几何优化，delta E参数是no，是否收敛
在DMol3里面进行优化，任务完成后给的提示是successfully-completed，但是从outmol文件里看最后的参数delta E No；|F|max Yes；|dR|max  Yes。

这种情况算收敛吗？可以提高精度进行下一步计算吗？

谢谢各位解答！

我用的是medium精度，具体参数如下
# Task parameters
Calculate                   optimize
Opt_energy_convergence       2.0000e-005
Opt_gradient_convergence    4.0000e-003 A
Opt_displacement_convergence  5.0000e-003 A
Opt_iterations             500
Opt_max_displacement       0.3000 A
Lattice_optimization       2 0.0150 A
Symmetry                   on
Max_memory                   2048
File_usage                   smart
Scf_density_convergence    1.000000e-005
Scf_charge_mixing          2.000000e-001
Scf_spin_mixing             5.000000e-001
Scf_diis                   6 pulay
Scf_iterations             1000
# Electronic parameters
Spin_polarization          unrestricted
Charge                      0
Basis                      dnp
Pseudopotential             ecp
Functional                   pbe
Aux_density                hexadecapole
Dftd                         G06
Integration_grid             medium
Occupation                   thermal 0.0100
Cutoff_Global                4.1000 angstrom
# Kpoint definition file (intervals/offset):
Kpoints                      file    2 2 2 0.0000 0.0000 0.0000
Co3O4.kpoints
# Calculated properties

作者
Author: 卡开发发 时间: 2018-2-5 22:49
算收敛。第一步的能量差deltaE=E1-E0，E0是估计出来的所以第一步的deltaE没有参考价值。力收敛了就是收敛了，但是否对应于极小点需要计算频率确定。

作者
Author: Channing 时间: 2018-2-7 16:11

卡开发发发表于 2018-2-5 22:49
算收敛。第一步的能量差deltaE=E1-E0，E0是估计出来的所以第一步的deltaE没有参考价值。力收敛了就是收敛了 ...

非常感谢您的回复！
还有几个问题，希望您有时间能解答：
请问如何进行 “是否对应于极小点需要计算频率决定” ？计算出频率后在Vibrational Analysis中怎样进行呢？

作者
Author: 卡开发发 时间: 2018-2-7 17:16

Channing 发表于 2018-2-7 16:11
非常感谢您的回复！
还有几个问题，希望您有时间能解答：
请问如何进行 “是否对应于极小点需要计算 ...

你计算完频率之后可以在Vibrational Analysis点calculate，之后频率会列出来。

作者
Author: Channing 时间: 2018-2-7 17:46

卡开发发发表于 2018-2-7 17:16
你计算完频率之后可以在Vibrational Analysis点calculate，之后频率会列出来。

感谢感谢！我知道可以得到一个frequency-intensity，还可以得到光谱图。
可是“极小点”怎么找？ “对应于极小点”是什么对应什么呀？

作者
Author: sobereva 时间: 2018-2-7 17:49

Channing 发表于 2018-2-7 17:46
感谢感谢！我知道可以得到一个frequency-intensity，还可以得到光谱图。
可是“极小点”怎么找？ ...

极小点结构靠几何优化得到

作者
Author: Channing 时间: 2018-2-7 18:03

sobereva 发表于 2018-2-7 17:49
极小点结构靠几何优化得到

感谢回答！

极小点是指几何优化过程中能量最低的点的结构吗？还是指优化结束时的结果？

我用DMol3进行优化，没有得到Energy.xcd和Convergence.xcd，以前用castep是有这两个文件的。所以如果极小点是能量最低的点的结构，这个极小点结构要怎么找呀？

作者
Author: sobereva 时间: 2018-2-7 18:07

Channing 发表于 2018-2-7 18:03
感谢回答！

极小点是指几何优化过程中能量最低的点的结构吗？还是指优化结束时的结果？

用最后一个点。优化过程中能量最低点应当就是优化过程最后一个点，如果不是，说明优化算法存在问题
Dmol3我不怎么用

作者
Author: Channing 时间: 2018-2-7 18:08

sobereva 发表于 2018-2-7 18:07
用最后一个点。优化过程中能量最低点应当就是优化过程最后一个点，如果不是，说明优化算法存在问题
Dm ...

get！十分感谢！

作者
Author: zorow 时间: 2018-2-7 18:46
请问楼主的软件是哪个版本，报价是多少？

作者
Author: 卡开发发 时间: 2018-2-7 19:02

Channing 发表于 2018-2-7 18:08
get！十分感谢！

一般看构型优化的最后一步，原则上最后一步能量是最低的，且要求最后一步达到构型优化的收敛标准。

你的情况应该是这样：你的结构一上来就能就能满足收敛标准，这样的话当然只有一步，有时候就不会跳出Energy.xcd和Convergence.xcd。如果你非要得到Energy.xcd和Convergence.xcd，你可以在DMol3 Analysis中去分析Energy Evolution，不过这样做的意义不大。

作者
Author: Channing 时间: 2018-2-7 20:24

zorow 发表于 2018-2-7 18:46
请问楼主的软件是哪个版本，报价是多少？

2017R2，报价不清楚，学校买的…

作者
Author: Channing 时间: 2018-2-7 20:33

卡开发发发表于 2018-2-7 19:02
一般看构型优化的最后一步，原则上最后一步能量是最低的，且要求最后一步达到构型优化的收敛标准。

你 ...

哇懂了懂了，非常感谢大家！

作者
Author: zorow 时间: 2018-2-7 20:33

Channing 发表于 2018-2-7 20:24
2017R2，报价不清楚，学校买的…

好的，感谢您的回复...

欢迎光临计算化学公社 (http://bbs.keinsci.com/)