自旋多重度的定义及常用软件设置方法.docx

本文格式为Word版，下载可任意编辑自旋多重度的定义及常用软件设置方法 dmol中设置磁构型有如下几种方法： 1）可以更改自旋多重度，根据体系的未配对电子数来更改多重度，S=n+1,譬如说体系的未配对电子数是2个，那么自旋 多重度只能取1， 3， 5， 7......，算完后，对比能量 2）在modify里更改上下自旋的搭配来实现不同磁构型，不过只能是共线型的，dmol不能考虑非共线磁布局 Enck: make p1后设置指定的反铁磁排布，然后再探索回对称性，但在探索对称性之前要在options里把formalspin勾上接着把Use formal spin as initial前打√，然后再算 Stou: 在优化自旋的时候，初值确定要取的适合，不然得到的自旋结果很可能是不正确的原子少的体系可能还不要紧，但是多原子体系的话很轻易展现这样的问题 Enck: 在bulid-symmetry-find symmetry里，就会看到options这时候，Use formal spin as initial前打√是务必的，作用是启用前面的那些预设自旋 encke,分外感谢你的扶助，另外额外想问你一个问题，对于自旋的设置，我一向都搞不太明白。

譬如对于我研究的氧化铁体系，我在structure中导入布局，make P1 后，在modify/electronic configuration中设置了每一层原子不同的自旋方向，我的问题是：当我开启electronic configuration对话框，里面不是有个formal charge吗？那个值显示的是3，是不是就是氧化铁中的铁的化合价？然后假设点automatic，就自动变成2了，相应的low spin 和high spin也相应的发生变化，2代表什么？到底应不理应在automatic前打上√呢？对于spin state中改选high还是low我也搞不领会，我通常都是分别计算这两种处境下的能量 看哪个能量低就选哪个，我不知道这样做可行不成行？ 期望你的答复，万分感恩！ Enck: 假设不是要求分外分外高，也就是正常处境下，计算软件本身就会对元素的根本性质举行判断，你选择自动加载电荷即可而对于自旋，你的做法是对的，理论上就理应做不同的设置找出最小值自旋中的high和low只不过是磁性元素在不同的体系中的磁矩大小，其实high和low本身没什么严格界限设置时，你可以手动设定磁性原子的最大可能的自旋（其实high就对应着最大自旋），即孤立原子时的自旋，譬如Fe，价电子为3d6 4s2，其最大自旋为10-6=4uB（假设不知道为什么的话再问）。

而实际上，在化合物中，磁性原子自旋往往低于该值，在磁性材料中Fe一般为2.2uB，所以也可以设置接近测验值的数值，也就是3uB当然，不管设置什么（1-4），你可能会察觉，他们都收敛于一致的基态，其能量变化不会太大，并且最终的磁矩也不会有太大差异那么，既然都能收敛到一致的基态，设定初始为0是否也可以呢，貌似设定为0时的处境不太一样，不会收敛到一致的基态这个时候，谁能量最低就取谁其实，对于一个未知体系，要获得确定的基态，还理应设定（不同的）铁磁、反铁磁等状态来计算，而这几种状态是确定要对比的，即在铁磁、反铁磁、亚铁磁以及顺磁态下对比能量上下当然了，假设结合测验，你可以很轻易通过测验判断或者摈弃某些磁性状态 至于电荷的设定，也可以采取类似的方法处理，看起来很繁杂，事实上往往可以确定某些参数使问题变得简朴本质上，电荷在一个原子周边的大小分布，都没有一个十足的值，说电荷密度，说轨道杂化，说处于成键态，都是一回事，因此也就没有了十足意义上的离子键、共价键或金属键，都是相对成分大小 记得有个重点议论贴特意议论铁磁的设置吧，可以翻出来看看 有spin不为0的磁性原子，且这些原子spin方向一致，这个是铁磁。

有spin不为0的磁性原子，且这些原 假设总的spin近似是0，这个是反铁磁 假设总的spin不为0，那么就是亚铁磁 【共享】自旋多重度的定义及常用软件设置方法 : zzgyb (站内联系TA) 发布: 2022-01-01 定义: 自旋多重度（spin multiplicity）＝2S+1, S=n*1/2,n为单电子数所以，关键是单电子的数目是多少 当有偶数个电子时，例如 O2，共有16个电子，那么单电子数目可能是0，即8个alpha和8个beta电子配对，对应单重态，但是也可能是有9个α电子和7个β电子，那么能成对的是7对，还剩2个α没有配对，于是n=2,对应的是多重度3同理还可以有多重度5，7，9， ...一般而言,是多重度低 的能量低,最稳定,所以,一般来说,偶数电子的体系多重度就是1但是也有例外，假设O2就是一个大家都知道的例子，它的基态是三重态，其单重态反而是激发态 所以，总结一下，就是 电子数目是偶数，未成对电子数目n=0,2,4,6,...自旋多重度是1,3,5,7,... 电子数目是奇数，未成对电子数目n=1,3,5,7,...自旋多重度是2,4,6,8,... 多数处境是多重度低的能量低，有时（更加是有“磁”性的时候，例如顺磁的O2，以及Fe啊什么的），可能会高多重度的能量低。

Gaussian中自选多重度设置 在用gaussian时,假设分子的基态不是单重态，可以通过在title section里指定spin来计算其他组态如三重态的能量 譬如以下的计算指定氧原子为单重态 -------------------------- # hf/6-31G(d) Oxygen atom 0 1 O 0.0000 0.0000 0.0000 -------------------------- 而以下的输入那么指定氧原子为三重态 -------------------------- # hf/6-31G(d) Oxygen atom 0 3 O 0.0000 0.0000 0.0000 -------------------------- VASP中自选多重度设置 在VASP中，操纵自旋的参数是INCAR中的ISPIN, ISPIN=1表示举行非自旋极化的计算，而ISPIN=2表示举行自旋极化的计算为了进一步得到我们所需要的组态，可以通过NUPDOWN指定自旋向上和自旋向下的电子数差，此方法与直接指定自旋多重度是一样的。

VASP里，C的参考态是S2P2，即C的基态，ISPIN=2理应是这个态 值得留神的是，在举行固定自旋多重态的计算中，假设初始的电荷密度的自旋极矩和INCAR中设置的NUPDOWN不同的话，收敛会很慢假设是从读入的初始波函数开头计算，那么没有类似问题 以下是vasp手册上的原话： There is a word of caution required: If NUPDOWN is set in the INCAR file the initial moment for the charge density should be the same. Otherwise convergence can slow down. When starting from atomic charge density (ICHARG=2), VASP will try to do this automatically by setting MAGMOM to NUPDOWN/NIONS. The user can of course overwrite this default by specifying a different MAGMOM (which should still result in the correct total moment). If one starts from the wavefunctions, the initial moment will be always correct, because VASP will ``push'' the required number of electrons from the down to the up component. If starting from a chargedensity supplied in the CHGCAR file (ICHARG=1), the initial moment is usually incorrect! Hyperchem中自选多重度设置 在Setup设置中举行设置 Charge和Spin Multiplicity在Ab initio和Semi-empirical计算的对话框中展现。

Charge：指定额外的净剩电荷额外电荷定义当前的分子系统是一个电中性系统，正电系统（阳离子），还是一个负电系统（阴离子） Spin Multiplicity：自旋多重度闭壳分子的自旋多重度为1（单重态）一个自由基，有一个未成对电子，自旋多重度为2（双重态）有两个未成对电子的系统，自旋多重度一般为3（三重态）然而在某些处境下，例如两个自由基，两个未成对电子也可能产生单态 State: 描述系统中价电子状态包括指定分子处于第一激发单态(Next lowest)或者Lowest电子态 Lowest 给定自旋多重度的最低电子态它不确定是基态 Next lowest 给定自旋多重度（单重、双重、三重或四重态）的第一电子激发态 在HyperChem 6.0中，semi-empirical方法通常需要给定多重度的最低能态（lowest）或者次最低能态（next lowest）的计算由于偶数个电子的分子没有未成对电子，是闭壳层单态，所以只有最低三重态是有效的，而次最低三重态是无效的例如，苯有偶数个电子，并且基态是闭壳层单重态我们可以计算基态（最低单重态），第一激发单重态（次最低单重态），或者第一激发三重态（最低三重态）。

也就是说，或者HOMO被两个电子占据，或者一个电子在HOMO，另一个电子在LUMO，产生了激发单重态或者三重态对双重态和四重态，只有给定多重度的最低态（lowest）可用UHF选项仅允许给定多重度的最低态（lowest）可用例如，可以用UHF选项研究苯的最低三重激发态，但是不能用来计算单重激发态这是由于HyperChem中的UHF选项不允许任意的轨道占据，也不允许CI计算对于RHF选项，可以计算CI波函数这个计算由一系列计算得到的RHF轨道开头，或者从最低单态（或双重态），或者从half electron的单态和三态（或双态和四重态）轨道从这些初始态可以得到参考态，从这些参考态，产生一系列微观态通过这些微观态适当的线性组合，可以获得一系列更切实的态 MS中自选多重度设置 Spin unrestricted: When checked, indicates that the calculation will be performed using different orbitals for different spins. This is known as a 'spin-unrestricted' or 'spin-polarized' calculation. If unchecked, the calculation uses the same orbitals fo。