24小时热门版块排行榜    

查看: 3070  |  回复: 19

fatpig8832

铁杆木虫 (著名写手)


[交流] 【求助】请问分子轨道系数的正负号全部调转表示发生了什么?

用molpro算激发态,中间要换一下轨道(rotate命令),我用一个C-C轨道替换前面的一个C-H轨道。

换完之后发现只是将要换的C-H轨道的系数正号变负号,负号变正号,系数的绝对值则几乎不变...这...是换了还是没换啊...

但是这样一搞,最后得到的几个激发态的数据跟我想要的数据很接近,所以我也不知道这样换正确与否...如果不这样搞,数据则差得比较远...

所以我想请教一下,正负号调转到底是啥意思?
回复此楼

» 猜你喜欢

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

» 本主题相关价值贴推荐,对您同样有帮助:

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

love254246267

木虫 (著名写手)



fatpig8832(金币+1):谢谢参与
期待高手解答,我也想知道
2楼2011-03-07 14:00:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mchen10

金虫 (正式写手)



fatpig8832(金币+1):谢谢参与
这是交坏两个轨道时用了rotation matrix, 旋转90度

    0    -1                 O1                -O2
[               ]   *   [          ]  =  [            ]
     1    0                 O2                 O1

至于为什么这么做,我也不知道
猜测是为了保持正交性
3楼2011-03-08 09:38:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mchen10

金虫 (正式写手)


突然想起来,不过也是猜的呵

矩阵交换两行之后 det(A)变成原来的负数 所以其中一行应取负以保持det(A)不变
4楼2011-03-08 09:52:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fatpig8832

铁杆木虫 (著名写手)


引用回帖:
Originally posted by mchen10 at 2011-03-08 09:38:51:
这是交坏两个轨道时用了rotation matrix, 旋转90度

    0    -1                 O1                -O2
[               ]   *   [          ]  =  [            ]
     1    0                 O2         ...

如果照这样解释的话,那应该是先用 C-C 轨道 替换掉 C-H 轨道的位置后,再将 C-C 轨道的所有正负号改变才对...

但事实上,C-C轨道 和 C-H轨道都呆在原来的位置上,且只有C-H轨道的所有正负号发生了改变...
5楼2011-03-08 10:01:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mchen10

金虫 (正式写手)


引用回帖:
Originally posted by fatpig8832 at 2011-03-08 10:01:27:
如果照这样解释的话,那应该是先用 C-C 轨道 替换掉 C-H 轨道的位置后,再将 C-C 轨道的所有正负号改变才对...

但事实上,C-C轨道 和 C-H轨道都呆在原来的位置上,且只有C-H轨道的所有正负号发生了改变...

所以说嘛。。我都是猜的

可以更精确地描述一下“原来的位置”吗?有点含糊

还有。。交换轨道只对occ closed之间的交换有意思

而最终打印出来的轨道顺序是按能量高低排的,所以你会看到顺序不变 (如果两个轨道都是occ的话)

这样回答合理吗?
6楼2011-03-08 10:16:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fatpig8832

铁杆木虫 (著名写手)


引用回帖:
Originally posted by mchen10 at 2011-03-08 10:16:24:
所以说嘛。。我都是猜的

可以更精确地描述一下“原来的位置”吗?有点含糊

还有。。交换轨道只对occ closed之间的交换有意思

而最终打印出来的轨道顺序是按能量高低排的,所以你会看到顺序不变 (如果 ...

是这样子D...

前面我指定了
closed,18;occ,26  (举个例子,大概是这样)

此分子有41个电子,所以占掉了21个轨道,其中第 21.1 号轨道是单电子

用orbprint命令打印出来,可以看到类似

(轨道,能量,轨道系数)
...
20.1     -0.9     ...
21.1     -0.8     ...
----------------------------我是分割线------------------------
22.1      0.5     ...
23.1      0.6     ...
24.1      0.8     ...
25.1      0.7     ...
26.1      0.6     ...
----------------------------我是分割线------------------------
27.1     -0.1     ...
28.1      0.1     ...
29.1      0.2     ...
...


那么 22.1~26.1 这5个轨道 是属于occ,但又没电子的轨道。且能量也不是按高低顺序排列的...21.1及之前,27.1及之后都是按顺序排列的...

现在要算激发态,假设 22.1 是个 C-H 轨道, 29.1 是个 C-C 轨道,我用
rotate,22.1,29.1,0
将这两轨道交换

结果就是 22.1 的轨道系数那里乘了(-1),其他所有东西几乎没发生改变(当然,小数点后5,6位还是有点变动的,忽略不计)。包括各个轨道的能量及顺序,均无变化,29.1的能量及轨道系数与rotate前一致...
7楼2011-03-08 11:29:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lihb734

铁杆木虫 (职业作家)



fatpig8832(金币+1):谢谢参与
引用回帖:
Originally posted by boylc789 at 2011-03-08 19:14:43:

请您不要经常纯表情回帖。
9楼2011-03-08 20:57:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mchen10

金虫 (正式写手)


旋转0度等于不旋转,是不是试着旋转90度
27能量为负 是不是应该换到occ里面

不好意思刚刚我说错了 交换生效可以是 occ 和 close之间,也可以是virtual 和 occ/close之间
10楼2011-03-09 01:17:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fatpig8832

铁杆木虫 (著名写手)


引用回帖:
Originally posted by mchen10 at 2011-03-09 01:17:38:
旋转0度等于不旋转,是不是试着旋转90度
27能量为负 是不是应该换到occ里面

不好意思刚刚我说错了 交换生效可以是 occ 和 close之间,也可以是virtual 和 occ/close之间

...这个...虽然程序里写着是0度,但实际上是90度...两种写法是等价的,在output file也可以看到它写了90度...

27.1 虽然能量很小,但它也是一个 C-H 轨道,所以换进去没什么用... 我也试过,换进去和没换结果毫无差别...

我想知道的就是,一个轨道的轨道系数全部乘了(-1)是否就是指转了90度?
如果是的话,那用来rotate的另一个轨道到底有啥用呢...迷惘中...
11楼2011-03-09 09:50:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mchen10

金虫 (正式写手)


引用回帖:
Originally posted by fatpig8832 at 2011-03-09 09:50:28:
...这个...虽然程序里写着是0度,但实际上是90度...两种写法是等价的,在output file也可以看到它写了90度...

27.1 虽然能量很小,但它也是一个 C-H 轨道,所以换进去没什么用... 我也试过,换进去和没换结果 ...

我也很想知道 等待牛人出现
12楼2011-03-09 10:37:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kekexili_08

木虫 (正式写手)



fatpig8832(金币+1):谢谢参与
一个轨道的轨道系数全部乘了(-1)

那就是什么都没发生

波函数本身是几率问题, 乘以常数, 还是自身。
13楼2011-03-09 11:25:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fatpig8832

铁杆木虫 (著名写手)


引用回帖:
Originally posted by kekexili_08 at 2011-03-09 11:25:19:
一个轨道的轨道系数全部乘了(-1)

那就是什么都没发生

波函数本身是几率问题, 乘以常数, 还是自身。

如果是这样的话,那rotate两个轨道其实相当于没有rotate了

但为什么算出来的结果会不一样呢?
14楼2011-03-09 12:03:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hexinchem

铁杆木虫 (著名写手)



fatpig8832(金币+1):谢谢参与
加油~~~~
15楼2011-03-09 12:10:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kekexili_08

木虫 (正式写手)


引用回帖:
Originally posted by fatpig8832 at 2011-03-09 12:03:01:
如果是这样的话,那rotate两个轨道其实相当于没有rotate了

但为什么算出来的结果会不一样呢?

rotate 可不是乘以常数那么简单, 直观想一下, rotate后, 电子云变了。
16楼2011-03-09 15:05:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luqing6879

木虫 (著名写手)



fatpig8832(金币+1): 谢谢参与
楼主你的问题解决了么?我最近也遇到相似的问题。
17楼2014-01-28 14:31:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fatpig8832

铁杆木虫 (著名写手)


引用回帖:
17楼: Originally posted by luqing6879 at 2014-01-28 14:31:44
楼主你的问题解决了么?我最近也遇到相似的问题。

没解决,因为是个小问题,所以解不解决无所谓...文章都已经发了好久了...

不知你遇到了什么问题?
18楼2014-01-28 16:14:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luqing6879

木虫 (著名写手)


引用回帖:
18楼: Originally posted by fatpig8832 at 2014-01-28 16:14:54
没解决,因为是个小问题,所以解不解决无所谓...文章都已经发了好久了...

不知你遇到了什么问题?...

跟你类似的问题。但是我想我可能解决了。不知道你还能不能找出你的输入文件。我猜你是不是只是写了rotate, orb1,orb2,0  这样的语句,而没有加start语句?

细看了下手册,要求交换轨道必须要先读进某个初猜。 所以要先start,record; 从某个记录上先把轨道读进来。
19楼2014-01-29 01:26:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fatpig8832

铁杆木虫 (著名写手)


引用回帖:
19楼: Originally posted by luqing6879 at 2014-01-29 01:26:52
跟你类似的问题。但是我想我可能解决了。不知道你还能不能找出你的输入文件。我猜你是不是只是写了rotate, orb1,orb2,0  这样的语句,而没有加start语句?

细看了下手册,要求交换轨道必须要先读进某个初猜。 所 ...

...不一定要加start的,加了可以从某个record读,不加则是默认从上一步的计算读...
20楼2014-01-29 22:19:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
boylc7898楼
2011-03-08 19:14   回复  
fatpig8832(金币+1):谢谢参与
相关版块跳转 我要订阅楼主 fatpig8832 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见