24小时热门版块排行榜    

查看: 6980  |  回复: 37
【奖励】 本帖被评价22次,作者zeoliters增加金币 20.2
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

zeoliters

木虫 (正式写手)


[资源] 【原创】-(专栏文章)化学反应体系中Broken symmetry Calculation的简要介绍

1. 什么是Broken symmetry Calculation?

“Broken symmetry”这一名称来源于“分子的对称性与波函数的对称性的差别”。一般的计算,可以根据处理电子的自旋方式不同而分成两种:close shell和open shell(spin polarization)。对于close shell,自旋向上的电子和自旋向下的电子是一样多且占据同一轨道,所以始终保持着同分子一样的对称性,而不存在Broken symmetry Calculation。 对于open shell,由于分别处理自旋向上和自旋向下的电子,就可以通过降低波函数的对称性使得波函数的对称性同分子的对称性是不一样的,即为Broken symmetry Calculation. (当然也可以保持其对称性一样,但是这就不是Broken symmetry Calculation.)

2. 为什么要使用Broken symmetry Calculation?

当描述分子的电子态的时候,对称性保持的计算已经不能够得到准确的能量了,所以引入了Broken symmetry Calculation使得计算更加准确。现以双自由基为例来介绍。
当计算自由基时,有两种可能的状态:单重态双自由基(singlet diradical,=1)和三重态的双自由基(triplet diradical, =2)。这样计算出来的结果有两种可能:
(1)三重态的能量比较低:这种情况下,检查是否在2左右。如果没有偏离太远,(个人理解不超过10%,yzcluster指出在10%以内已经是很大的了,对于分析体系的磁学性质,特别是局域磁矩时,有时是不可容忍的。因为本人没做过局域磁矩,对局域磁矩没有概念),这个三重态比较稳定。
(2)单重态的能量比较低:这种情况下,检查是否在1左右,如果是(个人理解偏差不超过10%),这说明此能量真是单重态的能量,并且此单重态比较稳定。但是经常出现大于1,比如说1.2,1.3等等。这说明此时得到的能量并不是单重态双自由基的能量。从数学的角度来说,之所以出现不等于1,就是因为波函数里面混进了三重态的波函数,波函数不再是自旋平方算符的本征函数。此时的波函数不能够用单个 slater行列式来表示,而是要通过用多个slater行列式来表示,解决的办法就是降低波函数的对称性,即使用Broken symmetry Calculation. 通过使用Broken symmetry Calculation.得到能量,再通过投影算符得到真正的单重态的能量。此能量才是真正的单重态的能量。

(未完,待续,待续内容:如何用Gaussian和VASP进行Broken symmetry Calculation)

[ Last edited by cadick on 2009-12-10 at 23:38 ]
回复此楼

» 收录本帖的淘帖专辑推荐

第一性原理 高斯计算总结 broken symmetry

» 猜你喜欢

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

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

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

yzcluster

金虫 (著名写手)


★★★★★ 五星级,优秀推荐

★ ★ ★ ★ ★
fegg7502(金币+5,VIP+0):thank you very much! 6-17 16:49
鼓励一下。但是有几个小问题
1.自旋污染在10%以内已经是很大的了,对于分析体系的磁学性质,特别是局域磁矩时,有时是不可容忍的。10%的范围有点大。
2.刚看到LZ的标题,我的第一反应是帖子描述的内容是凝聚态理论中的对称性破缺问题。在凝聚态物理的相变理论中,由于对称破缺,会造体系成序参量的出现,由高温无序相向低温有有序相转变。其实,对称性问题在量子场论中有更为广泛的概念。LZ标题可否再具体一些。
3.文中提到“对于open shell,自旋向上的电子和自旋向下的电子是一样多且占据同一轨道,所以始终保持着同分子一样的对称性,而不存在Broken symmetry Calculation。 对于close shell,由于分别处理自旋向上和自旋向下的电子,就可以通过降低波函数的对称性使得波函数的对称性同分子的对称性是不一样的”,是否说反了?希望LZ在好好CHECK一下。
4.个别句子可以再精确些,比如“波函数不再是自旋的本征函数”,最好改为“波函数不再是自旋平方算符的本征函数”,等等。
以上几点还望LZ考虑一下。当然,我的意见仅供参考。

[ Last edited by yzcluster on 2009-6-17 at 15:52 ]
6楼2009-06-17 15:50:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 38 个回答

zeoliters

木虫 (正式写手)


非常感谢提供建议!
7楼2009-06-17 15:53:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zeoliters

木虫 (正式写手)


非常感谢你的建议!第2,3,4非常中肯,第1个问题,我也根据你的建议进行了修改!
引用回帖:
Originally posted by yzcluster at 2009-6-17 08:50:
鼓励一下。但是有几个小问题
1.自旋污染在10%以内已经是很大的了,对于分析体系的磁学性质,特别是局域磁矩时,有时是不可容忍的。10%的范围有点大。
2.刚看到LZ的标题,我的第一反应是帖子描述的内容是凝聚态理 ...

[ Last edited by zeoliters on 2009-6-17 at 09:08 ]
8楼2009-06-17 15:58:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yzcluster

金虫 (著名写手)


★ ★ ★ ★
fegg7502(金币+4,VIP+0):谢谢,,那以后作者写出来之后先关闭回复,并置顶,请大家阅读如何? 6-17 16:51
对于专栏文章这个板块很看好。我觉得有两个改进的地方:
1.类似于PAPER的审稿制度,每个专栏文章出来之前最好能找专业相近的版主CHECK一下,以保证稿子的质量。
2.加大鼓励回帖的虫友,特别是能提出有效问题的虫友,这样才能都有所提高,同时可以使稿子的质量更高。不要回帖都是一堆表情,这是没有意义的,我想在小木虫我们不需要这些社会上浮躁的东西。

如果可能,尽量把这个东西做成一个特色精品,实际的解决科研中的某个问题,比如基本概念,基本物理OR化学图像,或是基本的操作等等。显然,问题也不用太大,但要实际解决。使之成为小木虫的特色。
最后,也佩服一下LZ吃螃蟹的精神。在这一点上,其他人(包括我)都落后了。
9楼2009-06-17 16:18:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 07:28 by ZPa0EcMwuECS
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 07:28 by ZPa0EcMwuECS
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-30 06:55 by ZPa0EcMwuECS
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:58 by ZPa0EcMwuECS
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 01:12 by ZPa0EcMwuECS
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-30 01:03 by ZPa0EcMwuECS
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-30 00:52 by ZPa0EcMwuECS
[基金申请] 国社科申报系统有变化 +3 kynobel 2026-08-23 3/150 2026-08-29 21:57 by 余韵清
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 为什么到现在没收到通知? +4 tannykie 2026-08-29 4/200 2026-08-29 17:14 by lmz0216
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见