[交流] 计算线性三原子分子的振动频率，弯曲振动频率是个虚频？已有9人参与

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

临渊羡鱼，不如退而结网

1楼2017-08-04 11:30:37

回帖支持 ( 显示支持度最高的前 50 名 )

wangf44

禁虫 (小有名气)

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

本帖内容被屏蔽

13楼2017-08-12 22:33:42

普通回帖

gengxin60

2楼2017-08-04 21:30:46

121099

银虫 (小有名气)

应助: 6 (幼儿园)
金币: 327.9
散金: 3
红花: 2
帖子: 76
在线: 76.4小时
虫号: 1299508
注册: 2011-05-18
专业: 理论和计算化学

★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
孙1986: 金币+3 2017-08-06 15:09:08

把角度的自由度打开重新优化一下结构看看

3楼2017-08-05 02:20:06

孙1986

金虫 (正式写手)

应助: 1 (幼儿园)
金币: 1046.9
散金: 1087
红花: 19
帖子: 455
在线: 243.1小时
虫号: 1405621
注册: 2011-09-17
专业: 原子和分子物理

引用回帖:

3楼: Originally posted by 121099 at 2017-08-05 02:20:06
把角度的自由度打开重新优化一下结构看看

好的，我试一下

临渊羡鱼，不如退而结网

4楼2017-08-06 15:09:27

孙1986

金虫 (正式写手)

应助: 1 (幼儿园)
金币: 1046.9
散金: 1087
红花: 19
帖子: 455
在线: 243.1小时
虫号: 1405621
注册: 2011-09-17
专业: 原子和分子物理

引用回帖:

3楼: Originally posted by 121099 at 2017-08-05 02:20:06
把角度的自由度打开重新优化一下结构看看

角度自由度是打开的

临渊羡鱼，不如退而结网

5楼2017-08-06 15:12:44

121099

银虫 (小有名气)

应助: 6 (幼儿园)
金币: 327.9
散金: 3
红花: 2
帖子: 76
在线: 76.4小时
虫号: 1299508
注册: 2011-05-18
专业: 理论和计算化学

★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
孙1986: 金币+3 2017-08-08 11:29:08

引用回帖:

5楼: Originally posted by 孙1986 at 2017-08-06 15:12:44
角度自由度是打开的...

三原子分子结构自由度数目并不多，如果几何优化给出的是全局最小的话（所有自由度都打开），那么就可以检查一下每个频率单点计算看看是否有收敛错的单点结果。

我个人碰到这样的情况的时候，一般选择用笛卡尔坐标系c1对称性优化，来确保优化是全局最小，然后再用优化所得笛卡尔坐标进行频率计算。

如果用的方法是DFT的话，个人经历过很多情况设置grid对几何优化及频率计算影响很大，可以选择更小的grid数值。
如果是WF办法的话，除了以上提到的问题暂时想不出其他的可能原因，你可以把input贴一下。

6楼2017-08-08 03:16:17

孙1986

金虫 (正式写手)

应助: 1 (幼儿园)
金币: 1046.9
散金: 1087
红花: 19
帖子: 455
在线: 243.1小时
虫号: 1405621
注册: 2011-09-17
专业: 原子和分子物理

引用回帖:

6楼: Originally posted by 121099 at 2017-08-08 03:16:17
三原子分子结构自由度数目并不多，如果几何优化给出的是全局最小的话（所有自由度都打开），那么就可以检查一下每个频率单点计算看看是否有收敛错的单点结果。

我个人碰到这样的情况的时候，一般选择用笛卡尔坐 ...

这是输入文件

TIM图片20170808113729.png

临渊羡鱼，不如退而结网

7楼2017-08-08 11:34:26

121099

银虫 (小有名气)

应助: 6 (幼儿园)
金币: 327.9
散金: 3
红花: 2
帖子: 76
在线: 76.4小时
虫号: 1299508
注册: 2011-05-18
专业: 理论和计算化学

★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
孙1986: 金币+4 2017-08-11 16:16:48

引用回帖:

7楼: Originally posted by 孙1986 at 2017-08-08 11:34:26
这是输入文件

TIM图片20170808113729.png
...

用DZ优化了一下结构。
发现如果打开夹角至175以下，则可以找到能量更低的PES极小点，
DZ优化结果为
R1=                1.20567667 ANGSTROM
R2=                1.89890745 ANGSTROM
A1=             123.08527193 DEGREE

分别对应以上的数值。

可能原因是180算是一个非常不稳定的局部极小点，几何优化步长太小跳不出这个极小点（优化步长为1度，非常小），而在频率计算时因为转换成了笛卡尔坐标系进行displancement计算，displancement为0.01，改变幅度较大，在PES的角度上来看harmonic frequency model已经不适用，所以出现虚频。
如果一定要计算角度为180的这个局部极小点，可以考虑多点计算PES，用多项式拟合振动波函数，再计算出harmonic frequency跟其他非简谐数值等。

» 本帖已获得的红花（最新10朵）

孙1986

8楼2017-08-11 04:43:28

孙1986

金虫 (正式写手)

应助: 1 (幼儿园)
金币: 1046.9
散金: 1087
红花: 19
帖子: 455
在线: 243.1小时
虫号: 1405621
注册: 2011-09-17
专业: 原子和分子物理

送红花一朵

引用回帖:

8楼: Originally posted by 121099 at 2017-08-11 04:43:28
用DZ优化了一下结构。
发现如果打开夹角至175以下，则可以找到能量更低的PES极小点，
DZ优化结果为
R1=                1.20567667 ANGSTROM
R2=                1.89890745 ANGSTROM
A1=             ...

这是一个排斥激发态的的结构，可能确实不好处理。但是按照楼主的讲解不太会处理，不过楼主的讲解倒是给了我一点启发，我先把角度稍微改一下，再试试。

临渊羡鱼，不如退而结网

9楼2017-08-11 16:23:46

wangf44

禁虫 (小有名气)

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

本帖内容被屏蔽

» 本帖已获得的红花（最新10朵）

孙1986

10楼2017-08-11 16:51:26