标题: atomicrex数值导数与解析导数能量梯度相差太大的报错该如何解决? [打印本页] 作者Author: zhaoyanx4 时间: 2022-8-13 04:34 标题: atomicrex数值导数与解析导数能量梯度相差太大的报错该如何解决? 最近在尝试用atomicrex拟合势参,但程序运行时总报如下错误:
Function pair (r=10) Numerical derivative: 1680.358481 (err: 960.2007266) Analytical derivative: 0 Deviation: -inf
ERROR: Numerical derivative deviates too much from the analytical energy gradient computed by function "pair"
翻译过来是说数值导数与我在输入文件里提供的一阶导数值相差太大。尝试修改文件的过程中发现报错中的numerical derivative值主要受cutoff影响,其次势参初值会造成微弱影响,但根据我使用的势函数直接计算出的值来看,报错中numerical derivative对应的r值离cutoff很远,且该r值下势函数一阶导值也不为0。被这个问题困扰了好几天,请问有大佬知道这个问题该如何解决吗?附上一部分势函数输入文件
(, 下载次数 Times of downloads: 5)