24小时热门版块排行榜    

Znn3bq.jpeg
北京石油化工学院2026年研究生招生接收调剂公告
查看: 3106  |  回复: 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
回复此楼

» 猜你喜欢

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

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

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

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫


★ ★
gkf高: 金币+2 2013-12-25 21:48:53
引用回帖:
8楼: Originally posted by 科学岛 at 2013-12-24 21:15:49
1.那非Ryd的自旋布居数不归一化是不是也可以?
2.对于开壳层自旋污染<S2>是从能量就算的输出文件中得到,还是 从NBO计算的文件中得到?另外<S2>算符有什么特殊的含义吗?...

1 只要偏差不太大,也可以不做归一化。
2 这和NBO模块根本无关。直接取量化程序自身输出的。这就是角动量平方的期望值,对于判断自旋污染之类的问题有用,所有量化书里都会介绍
9楼2013-12-24 22:47:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫


【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与,应助指数 +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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

科学岛

铁虫 (初入文坛)

引用回帖:
2楼: Originally posted by sobereva at 2013-12-24 12:11:47
体系有n个自旋平行的未配对电子,分子电子自旋磁矩就为n*μ_B。而自旋相反的电子的自旋磁矩会抵消,所以闭壳层体系无电子自旋磁矩。通过布居分析(如Mulliken、NPA等),得到每个原子,或每个原子轨道上的自旋电子数 ...

谢谢!
一看就知道您对这方面很了解。我刚接触这方面的内容,按我对您回复的理解,以下几点是否正确:
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

科学岛

铁虫 (初入文坛)

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
  • 2013-12-24 16:42:43, 527.6 K
  • 附件 2 : reopt_reopt_Neutral_Geom_6_.gjf
  • 2013-12-24 16:43:03, 4.76 K
4楼2013-12-24 16:43:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 266调剂 +5 daya sun 2026-04-07 6/300 2026-04-08 00:08 by daya sun
[考研] 材料307分求大佬组收留 +6 Hll胡 2026-04-07 6/300 2026-04-07 23:06 by jp9609
[考研] 材料调剂 +13 汉123456 2026-04-07 14/700 2026-04-07 22:53 by 来看流星雨10
[考研] 专硕085403,291分,有两篇专利,一国一奖 +3 哈吉咪哈吉咪 2026-04-07 3/150 2026-04-07 18:21 by 蓝云思雨
[考研] 306求调剂 +3 15287505595 2026-04-03 3/150 2026-04-07 18:08 by 蓝云思雨
[考研] 301求调剂 +18 121. 2026-04-04 18/900 2026-04-07 17:49 by 蓝云思雨
[考研] 333求调剂 +6 合乘杨习夕 2026-04-06 6/300 2026-04-07 09:44 by 猪会飞
[考研] 334分控制工程求调剂 +4 姜尚真sadasd 2026-04-03 4/200 2026-04-07 09:26 by 蓝云思雨
[考研] 复试调剂 +14 呼呼?~+123456 2026-04-05 14/700 2026-04-06 22:50 by chenzhimin
[硕博家园] 0856材料化工求调剂,一志愿211,初试成绩349 +4 江淮北月 2026-04-05 4/200 2026-04-06 22:44 by chenzhimin
[考研] 一志愿211,化学学硕,310分,本科重点双非,求调剂 +13 努力奋斗112 2026-04-04 13/650 2026-04-06 07:13 by jj987
[考研] 308求调剂 +4 maverick^_^ 2026-04-03 4/200 2026-04-05 19:08 by 蓝云思雨
[考研] 一志愿北京交通大学材料工程总分358求调剂 +4 cs0106 2026-04-04 4/200 2026-04-05 18:46 by imissbao
[考研] 348求调剂 +6 wukira 2026-04-04 6/300 2026-04-05 18:11 by 猪会飞
[考研] 一志愿上海海洋大学083200食品学硕,求调剂,接受其他专业083200 +4 what张 2026-04-04 5/250 2026-04-05 14:07 by chw1980_0
[考研] 085600调剂 +9 东照照照 2026-04-04 9/450 2026-04-05 13:44 by ujn_zhuj
[考研] 怎么删帖子啊 +3 缝曦1000 2026-04-04 3/150 2026-04-04 14:20 by 土木硕士招生
[考研] 求调剂不挑专业 +3 xrh030412 2026-04-01 3/150 2026-04-03 14:40 by 氮气气气
[考研] 一志愿北京科技大学085601材料工程英一数二初试总分335求调剂 +8 双马尾痞老板2 2026-04-02 9/450 2026-04-02 14:45 by 5896
[考研] 275学硕081000服从调剂到其他专业,保不住本专业了 +7 一只小小水牛 2026-04-02 8/400 2026-04-02 14:23 by alice-2022
信息提示
请填处理意见