vverawei 发表于 2022-3-1 14:16 您好，我也遇到了和您相同的问题，可以加个微信交流一下吗，15735717121 |
卡开发发 发表于 2022-3-1 12:07 好哒好哒，谢谢老师的耐心解答，十分感谢！ |
vverawei 发表于 2022-3-1 09:24 1、不+U对于一些体系有可能会估计出来低自旋态或无自旋态（NM）更稳定。这种稳定性一般看自洽场能量可以判断，用杂化泛函也是为了避免U值的影响来讨论这件事。不见得所有的过渡金属的化合物都有磁性，比如TiO2。当然，如果实验上能确定磁序是什么，如果能直接判断肯定更好。 2、不同磁性就意味着可能是基态也可能不是基态，对于电子结构相关的问题影响很大，比如能带图和DOS图，热电因为也牵涉到，所以也会受到影响。 3、是要这么做，如果实验上能确定结构的甚至可以简单拿单点能比较，因为确实个别磁序会出现算不出结果的情况，比如SCF收敛不了或者收敛不到这个磁序上。 |
本帖最后由 vverawei 于 2022-3-1 09:59 编辑 卡开发发 发表于 2022-3-1 08:49 老师不好意思再打扰您一下，还有一些疑问麻烦您多指教。1.您刚刚的意思是对于闭壳层体系如3d10来说，加不加U值对其磁性按理说影响不大（我是都不+U来比较的），如果我想进一步探究到底是哪种磁性（顺磁还是反铁磁）的话是需要用杂化泛函来验证对吗？那如果我想知道她是否具有磁性该如何判断呢？是不是对于过渡金属来说都是具有磁性的？ 2.如果对磁性的判断不准确，即不确定是顺磁还是反铁磁设置ISPIN=1或者是2对于其能带图、dos图以及热电图会有很大影响吗？ 3.对于Cu3C6S6这样的导电配位聚合物，它是二价铜即3d9的体系，一定有自旋（设ISPIN=2)对吗？我把它分为两种情况设置magmom，即3*1（铁磁）和2*1 1*-1（铁磁）这两种，分别进行结构优化比较 能量高低来判断采用哪种，这种方法对吗？ |
卡开发发 发表于 2022-3-1 08:49 好滴，谢谢老师您的指导 |
本帖最后由 卡开发发 于 2022-3-1 08:51 编辑 vverawei 发表于 2022-3-1 08:28 建议用杂化泛函的原因就在这：其实很难保证当前的U值做出来的结果是NM磁序稳定但更大的U却是FM或AFM磁序稳定，而且极端情况下NM可能和FM或AFM磁序对应的最佳U值并不相同（虽然这种情况比较少见），从而有给出这两种磁序稳定性估计出来的顺序与实际相反。尽管杂化泛函也存在精确交换的比例任意性的问题，但毕竟普适性还是比DFT+U的U值更好一些，当然杂化泛函计算量确实很大，计算条件允许的情况建议杂化泛函检验下便于说服审稿人。 |
卡开发发 发表于 2022-2-28 22:51 好的老师，谢谢您。我还是有几个疑问，烦请老师多多指教，想请问一下这样的情况我是不是可以不用开自旋了？说明它是顺磁性对吗？在没有使用杂化泛函的情况下我不断增大U值发现它的带隙在变大，这样的话还需要再用杂化泛函去再做一次吗？ |
vverawei 发表于 2022-2-28 21:58 d10满层应该DFT+U按道理影响不会太大，但是慎重起见体系不大还是可以杂化泛函验证下。 |
卡开发发 发表于 2021-5-1 09:08 老师您好，请问在一个一价铜的导电聚合物体系中，亚铜离子为3d10，没有单电子存在，这个时候我用LDAUTYPE=2分别设置ispin=1或者是2来做结构优化和SCF计算，后比较能量发现不开自旋的能量更低，那么此时我可以不开自旋但是+U的方法来继续算吗？这样是不是就说明该导电配位聚合物是顺磁性的？ |
卡开发发 发表于 2021-5-1 09:08 十分感谢！ |
XWH 发表于 2021-4-30 13:41 DFT+U原则上说并不是一定要针对有自旋极化的系统，比如非自旋极化的系统需要通过Hubbard U修正的例子也不少，比如TiO2。CASTEP界面上限定了DFT+U必须开启自旋极化，但我不确定是否能够不通过图形界面而是简单修改*.param文件使得DFT+U不去使用自旋极化就能够进行。如果实在是无法这样操作，你可以开启自旋极化并且对每个+U的原子磁矩猜测为0，如果体系真的没有磁性，这样计算的结果也不会有磁矩。 CASTEP某些方面的功能确实太弱了。 |
| 参与人数Participants 1 | eV +3 | 收起 理由Reason |
|---|---|---|
|
| + 3 |
jiangning198511 发表于 2021-4-30 17:03 是这意思，感谢！ |
XWH 发表于 2021-4-30 16:45 你的意思是用第一性原理得到的结果去拟合出势函数，然后再跑分子动力学吗？这个我没做过，但好多文献有提到，拟合参数一般就是通过第一性原理的结果或实验值来完成 |
jiangning198511 发表于 2021-4-30 16:39 好的，我再看看，可以再问个题外话吗？没有人能根据MS算某一个氧化物相关性质参数，然后去写势函数，用到lammps里面这种的吧？ |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-20 00:52 , Processed in 0.192741 second(s), 30 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites