计算化学公社

标题: 关于加电场的问题 [打印本页]

作者
Author: debai 时间: 2015-11-16 15:04
标题: 关于加电场的问题
在高斯加电场计算时，手册说最小是10的负4次方au的电场，能否加10的负5次方及以下的电场，如果能加格式是什么呢？

作者
Author: sobereva 时间: 2015-11-16 23:47
关键词直接用诸如field=x+300这种格式最少也得加0.0001a.u.
如果要更小，可以用field=read，在分子坐标末尾空一行写外加电场的XYZ分量，此时可以小于0.0001。

作者
Author: 十年磨练 时间: 2015-11-18 09:01

sobereva 发表于 2015-11-16 23:47
关键词直接用诸如field=x+300这种格式最少也得加0.0001a.u.
如果要更小，可以用field=read，在分子坐标末 ...

sob老师，gaussian里面可以加交变电场吗？

作者
Author: sobereva 时间: 2015-11-18 09:13

十年磨练发表于 2015-11-18 09:01
sob老师，gaussian里面可以加交变电场吗？

不能

作者
Author: 十年磨练 时间: 2015-11-18 10:40

sobereva 发表于 2015-11-18 09:13
不能

嗯嗯，谢谢sob老师

作者
Author: debai 时间: 2015-11-18 17:09
谢谢！

作者
Author: wuyiyyyy 时间: 2020-6-8 13:48

sobereva 发表于 2015-11-16 23:47
关键词直接用诸如field=x+300这种格式最少也得加0.0001a.u.
如果要更小，可以用field=read，在分子坐标末 ...

是不是说field=x+1，这个沿着X正方向施加0.0001a.u.的电场？
field=x+100，和field=read这两种写法有没有什么区别呀？我改用第二种写法有时候会报错

作者
Author: sobereva 时间: 2020-6-8 18:53

wuyiyyyy 发表于 2020-6-8 13:48
是不是说field=x+1，这个沿着X正方向施加0.0001a.u.的电场？
field=x+100，和field=read这两种写法有没 ...

此时是沿着负方向施加0.0001a.u.的电场

field=read允许你直接定义外电场矢量，报错显然是你格式写错了

作者
Author: debai 时间: 2023-2-22 18:30

sobereva 发表于 2015-11-16 23:47
关键词直接用诸如field=x+300这种格式最少也得加0.0001a.u.
如果要更小，可以用field=read，在分子坐标末 ...

sob老师，再问一下，按照field=read，在分子坐标末尾空一行写0.00000 0.00000 0.00144或者0.00145，看输出文件时写的是 An electric field of 0.0000 0.0000 0.0014，是程序默认不认14后面的数吗？谢谢！

作者
Author: wzkchem5 时间: 2023-2-22 20:18

debai 发表于 2023-2-22 11:30
sob老师，再问一下，按照field=read，在分子坐标末尾空一行写0.00000 0.00000 0.00144或者0.00145，看输 ...

应该是显示的有效数字问题，而不是读取的有效数字问题。
高斯是用fortran写的，fortran读取一个完整的数的代码比读取一个数的前几位的代码要好写，所以高斯没有理由只读到0.0014。但是输出一个数是相反的，输出小数点后固定几位的代码比输出所有非零数字的代码要好写，所以一个小数在显示的时候很可能被四舍五入，即便是有限小数

作者
Author: debai 时间: 2023-2-23 09:19

wzkchem5 发表于 2023-2-22 20:18
应该是显示的有效数字问题，而不是读取的有效数字问题。
高斯是用fortran写的，fortran读取一个完整的数 ...

感谢！

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