24小时热门版块排行榜    

查看: 3707  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 时间戳又变了 +11 wuchongjun 2026-08-20 16/800 2026-08-20 19:26 by wuchongjun
[基金申请] 科研孤儿太难了 +9 我4大白菜 2026-08-20 9/450 2026-08-20 18:23 by 2549091707
[基金申请] 只有每年这种时候来逛逛小木虫 +21 yaoyewhu2008 2026-08-20 22/1100 2026-08-20 17:03 by yaoyewhu2008
[基金申请] 今天系统多次维护,明天很可能放榜! +9 zju2000 2026-08-16 10/500 2026-08-20 15:56 by Leogzhya
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +16 Ldrop2023 2026-08-13 20/1000 2026-08-20 11:07 by xskun
[基金申请] 今天基金会出结果吗?20260819 +14 kkkl_v 2026-08-19 15/750 2026-08-20 08:21 by sunzitan
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
请填处理意见