| 查看: 1129 | 回复: 1 | |||||
beefly专家顾问 (职业作家)
地沟油冶炼专家
|
[交流]
Gaussian计算超极化率的符号错误 已有1人参与
|
|
今天偶然读到下面的文章 曾薇，丁培江，赵可清，Gaussian程序计算的一阶超极化率的符号问题，四川师范大学学报（自然科学版），33(2)，228，2010。可下载 http://www.paper.edu.cn/index.ph ... 395(2010)02-0228-03 文章说Gaussian程序计算的一阶超极化率差一个负号。一开始感觉不可思议，因为感觉Gaussian不大可能犯这种低级错误，而且即使犯了，也会很快被众多用户发现并得到及时改正，不可能保留这么多年（如果从已知可以计算超极化率的G92算起，那么也有20年了）。 但是经过简单的计算后，发现作者是对的。这里使用不同于文献中的算法，直接从定义出发，即，一阶超极化率（beta）等于能量对电场的三阶导数，再乘以-1。可以用5点公式计算数值三阶导数。五个点的能量和电场位移量已经在文章的表2给出，三阶导数公式的因子可以在这里查到（第一个表） http://en.wikipedia.org/wiki/Finite_difference_coefficient MATLAB流程如下： e=[-15.2698590257 -15.2732037424 -15.2768727762 -15.2808584389 -15.2851530988]; c=[-.5 1 0 -1 .5]; beta=-c*e/(0.005^3) 计算得到的beta=-61.2796与表1中FF结果非常接近。有一点小差异是因为从公式(2)计算的一阶超极化率包含了高阶项的贡献。可见Gaussian程序计算的一阶超极化率确实差一个负号。也许有人会认为Gaussian程序使用了不同的定义，但是迄今为止，文献中所有计算能量对电场n阶导数的公式都有一个负号。既然偶极矩（1阶导数）和极化率（2阶导数）公式都有负号，对一阶超极化率（3阶导数）何必用不同的定义呢？实际上，Gaussian的超极化率文献用的就是带负号的公式。 其实不光文章中指出的CPHF和ENONLY方法（正确的叫法应该是解析方法和纯数值方法），利用偶极矩数值求导的半解析方法一样会得到错误的超极化率符号。也就是说，Gaussian程序超极化率符号的错误与算法无关。这个错误在Gaussian09 C.01中依然存在。 经过检查，还发现GAMESS(US)程序计算一阶超极化率的符号是正确的。 [ Last edited by beefly on 2012-11-1 at 03:14 ] |
回复此楼» 猜你喜欢
散金求好运
已经有315人回复
-
博士生招生, 材料模拟方向，University of South Carolina
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有273人回复
-
Postdoc position, materials simulations/machine learning, U. of South Carolina
已经有0人回复
-
祈福青基
已经有34人回复
-
瑞典-博后-固态电池-粘结剂-锂离子电池方向
已经有53人回复
-
牛来！米来！面来！
已经有8人回复
-
广东省面上有什么动态么，啥时候公布呀？
已经有1人回复
-
2027申博咨询
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 介电函数曲线与电子跃迁
已经有7人回复
- 关于键级和AIM分析
已经有4人回复
- 我理论计算出了ECD，但不知道如何画分子轨道图
已经有6人回复
- 关于使用Multiwfn计算ELF的Basin积分的问题
已经有6人回复
- 用cbs-qb3计算I时出错
已经有9人回复
- 使用Multiwfn基于完全态求和(SOS)方法计算极化率和超极化率
已经有8人回复
- 使用Multiwfn分析Gaussian09的极化率、第一超极化率的输出
已经有12人回复
- 关于LANL2DZ加f极化函数的问题
已经有5人回复
- 高斯算频率的问题
已经有5人回复
- 高斯能不能计算分子极化率
已经有12人回复
- 高斯软件文件存储chk问题
已经有17人回复
- Gaussian中输出文件哪一个是可以用于下一步计算的文件啊
已经有5人回复
- 【小卒个人文集】用NLO-Calculator来计算二阶非线性光学响应系数张量及其各向异性
已经有32人回复
- 关于化合物超极化率的计算方法
已经有6人回复
- 请问在Gaussian计算中怎样考虑弱相互作用？
已经有23人回复
- 【求助】高斯计算td的问题
已经有9人回复
- 【求助】拉曼 极化率
已经有7人回复
- 【求助】审稿意见求助
已经有9人回复
- 【求助】Gaussian09作业中断，可否接着断点继续计算？
已经有11人回复
- 【求助】优化计算后偶极距选择
已经有12人回复
- 【求助】如何在以下高斯03输出文件中得到二阶超极化率
已经有15人回复
lihb734
铁杆木虫 (职业作家)
站在计算化学入门的门槛上
- QC强帖: 2
- 应助: 132 (高中生)
- 金币: 7495.4
- 散金: 9378
- 红花: 74
- 帖子: 3107
- 在线: 1177.3小时
- 虫号: 448403
- 注册: 2007-11-01
- 性别: GG
- 专业: 理论和计算化学
2楼2012-11-01 09:23:37
30