| 查看: 1017 | 回复: 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 ] |
回复此楼» 猜你喜欢
武汉纺织大学全国重点实验室杨应奎教授团队拟招17-19位硕博生！
已经有0人回复
-
上海工程技术大学生态与功能纺织品开发研究团队招收学硕专硕研究生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有134人回复
-
青岛大学化学化工学院2026年招收调剂硕士研究生
已经有0人回复
-
新能源材料方向考研调剂招生（长江师范学院）
已经有0人回复
-
0854电子信息272求调剂！
已经有0人回复
-
化学-化工-材料-食品-生物-环境2026年研究生调剂（可以不考数学）
已经有0人回复
-
化学-化工-材料-食品-生物-环境-农业 2026年研究生调剂（可以不考数学）
已经有1人回复
-
070300一志愿211，312分求调剂院校
已经有4人回复
-
26年储能、电池方向博士申请
已经有0人回复
-
070300化学求调剂
已经有2人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 介电函数曲线与电子跃迁
已经有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
