| 查看: 3058 | 回复: 12 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
如何得到原子中某个轨道对分子总磁矩的贡献？ 已有1人参与
|
|||
|
在Gaussian 中用pop=full,即可得到磁矩的数据。 Do spin unrestricted (or restricted open-shell) calculations (e.g.UB3LYP/6-311+G*) to get spin moments for each orbital printed. Be careful about spin-multiplicity though. Use pop=full or pop=regular in the route card for population analysis as well. Spin moments will be printed below Gross population analysis section. You can calculate total magnetic moment from that. 具体网址：http://www.ccl.net/chemistry/res ... .002-dir/index.html 但我想知道某个具体轨道对磁矩的贡献，知道怎么分析吗？例如下图中Au的5d轨道对总磁矩的贡献怎么分析？谢谢！  QQ截图20131224102322.jpg |
回复此楼» 猜你喜欢
PDF在线压缩软件
已经有5人回复
-
材料/化学相关专业2026级学术/专业型硕士研究生
已经有19人回复
-
物理化学论文润色/翻译怎么收费?
已经有175人回复
-
求催化作用基础第三版 （甄开吉) 课件
已经有3人回复
-
26年材料/化学相关专业博士研究生招生 温州大学化学与材料工程学院(26年秋季入学)
已经有0人回复
-
求助电催化剂制备，怎么在碳布表面生长PdCo合金
已经有0人回复
-
求助电催化剂制备，怎么在碳布表面生长PdCo合金
已经有0人回复
-
304求调剂
已经有3人回复
-
上海华东理工大学机械与动力工程学院先进材料与制造研究组招聘公告
已经有0人回复
-
304求调剂
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 为什么某些双原子分子不能把全部电子填入分子轨道
已经有3人回复
- Castep如何计算原子磁矩沿x，y，z轴的分量？
已经有4人回复
- 如何计算分子轨道中各元素对应原子轨道的组成比例
已经有4人回复
- 原子轨道合成分子轨道的问题
已经有3人回复
- 原子磁矩的计算问题
已经有11人回复
- 如何计算原子某轨道的能量？
已经有12人回复
- vasp可以固定单个原子的磁矩吗？
已经有5人回复
- 有哪位高人知道单个原子磁矩在铁磁相和反铁磁相的大小有变化么？
已经有5人回复
- 【求助】关于分子对称性、群论和原子轨道
已经有7人回复
- 【求助】我想计算一个分子中C原子的p轨道键长，用什么软件计算比较好。
已经有6人回复
- 【求助】如何得到能带在每个原子和轨道上的投影？
已经有10人回复
- 【求助】如何判断分子轨道是由各个原子的哪个轨道形成的
已经有11人回复
- 【转帖】MS中原子对前线轨道贡献的计算实例
已经有29人回复
sobereva
至尊木虫 (著名写手)
本人已永久离开小木虫
- QC强帖: 17
- 应助: 226 (大学生)
- 贵宾: 5.4
- 金币: 2180.5
- 红花: 217
- 帖子: 1241
- 在线: 366.8小时
- 虫号: 281271
- 注册: 2006-09-23
- 专业: 宇宙学
★ ★ ★
zhou2009: 金币+3 2013-12-24 18:55:54
zhou2009: 金币+3 2013-12-24 18:55:54
|
1 不是这么简单相乘。比如某轨道上alpha电子布居数是na，beta是nb，那么自旋布居就是na-nb，(na-nb)*μ_B才是这个轨道的电子贡献的磁矩。后面两句话你的理解正确。 2 高斯输出的Gross orbital populations是基于Mulliken方法得到的基函数的布居数，这种计算方式很不推荐使用。尤其是带了弥散函数的情况更不能用（比如CCL上那个人拿6-311+G*举例就很不合适），NPA得到的轨道布居数有意义得多。Mulliken方法的缺陷在《谈谈轨道成份的计算方法》（http://hi.baidu.com/sobereva/item/d52876160110ca8d89a956c5）、《原子电荷计算方法的对比》（http://www.whxb.pku.edu.cn/EN/10.3866/PKU.WHXB2012281）上都有讨论 3 问题不在于程序，而在于原理 4 我的文章里没有专门讨论磁矩的 5 自旋污染几乎是众所周知的问题 |
6楼2013-12-24 18:05:46
sobereva
至尊木虫 (著名写手)
本人已永久离开小木虫
- QC强帖: 17
- 应助: 226 (大学生)
- 贵宾: 5.4
- 金币: 2180.5
- 红花: 217
- 帖子: 1241
- 在线: 366.8小时
- 虫号: 281271
- 注册: 2006-09-23
- 专业: 宇宙学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
zhou2009: 金币+5 2013-12-24 13:58:11
科学岛: 金币+100, ★★★★★最佳答案 2013-12-25 08:26:01
感谢参与，应助指数 +1
zhou2009: 金币+5 2013-12-24 13:58:11
科学岛: 金币+100, ★★★★★最佳答案 2013-12-25 08:26:01
|
体系有n个自旋平行的未配对电子，分子电子自旋磁矩就为n*μ_B。而自旋相反的电子的自旋磁矩会抵消，所以闭壳层体系无电子自旋磁矩。通过布居分析（如Mulliken、NPA等），得到每个原子，或每个原子轨道上的自旋电子数，就能讨论它们各自对分子磁矩的贡献。比如某d原子轨道上自旋电子数为0.4，就说这个d原子轨道上电子自旋磁矩是0.4μ_B。也可以根据电子自旋密度讨论，某点自旋密度为x，就可以说这个点的电子对分子磁矩的贡献是x*μ_B。也因此，可以将空间划分为不同子区域，根据自旋密度在相应区域积分值讨论相应区域内的自旋磁矩。μ_B=e*h_bar/(2m_e)这个量被定义为玻尔磁子。 高斯布居分析只是给出每个基函数的布居数，怎么和原子轨道对应上需要基于你对所用基组的了解。pop=NPA给出的NPA分析结果则可以直接给出每个原子轨道的布居数。由于你用了赝势，所以NPA分析输出的哪个是5d需要自行判断。 另外，分子磁矩是电子自旋磁矩、电子轨道运动产生的磁矩和原子核自旋磁矩共同贡献的，通过电子布居分析讨论的只是电子自旋磁矩部分，一般总磁矩也是由此主导。 |
2楼2013-12-24 12:11:47
|
谢谢！ 一看就知道您对这方面很了解。我刚接触这方面的内容，按我对您回复的理解，以下几点是否正确： 1.NPA得到各原子的各电子层上的电子布居数n，相应的磁矩就是n*μ_B了？某原子的磁矩就是各轨道磁矩之和，各原子磁矩之和就是团簇的总磁矩？ 2.http://www.ccl.net/chemistry/res ... .002-dir/index.html 中说是根据Gross orbital populations之后的数据来计算磁矩，现在主流的是采用“NPA得到各原子的各电子层上的电子布居数n，相应的磁矩就是n*μ_B”还是"Gross orbital populations之后的数据来计算磁矩"思路。 3.网上查询了一下，也有用MS计算磁矩的， 4.您自己有没有这方面的已经发表的文章，到时可以引用引用。 5.另外，您有没有听说过自旋污染？ |
3楼2013-12-24 15:42:06
|
1.附件中分别是NBO计算的输入和输出文件，是不是总的“ NATURAL POPULATIONS: Natural atomic orbital occupancies” 减去Alpha spin orbitals 的“ NATURAL POPULATIONS: Natural atomic orbital occupancies”就是相应的beta 自旋轨道的，之后Alpha，beta spin orbitals 的电子布局数之差就是相应原子轨道的未配对电子？ 2.另外，总磁矩不需要考虑激发态吧， 例如下面的就不需要考虑了吧， 3 Au 1 S Ryd( 7S) 0.00182 2.01677 4 Au 1 S Ryd( 8S) 0.00003 8.33001 5 Au 1 S Ryd( 9S) 0.00000 70.15820 2.若不对，能否在输出文件中指出，万分谢谢！ |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : reopt_reopt_Neutral_Geom_6_.out
- 附件 2 : reopt_reopt_Neutral_Geom_6_.gjf
2013-12-24 16:42:43, 527.6 K
2013-12-24 16:43:03, 4.76 K
4楼2013-12-24 16:43:06
50