Eudaimonia 发表于 2022-2-4 10:16 老师 您好 请问晶胞长度3 3 11 k点5 5 2 是怎么得到的呀 |
卡开发发 发表于 2022-2-17 14:11 好的，谢谢老师！ |
不会扣篮的后卫 发表于 2022-2-17 13:57 2、应该也有，你可以按照化学式来检索，然后再mag order那一行里面有显示磁序。 4、建模的时候修改xyz和你的cell就行，或者也可以通过先转换为其他格式用可视化程序修改好再转换回来，可视化程序可选的有ase的gui，或者你也可以用M$处理好导出cif。 |
卡开发发 发表于 2022-2-16 13:29 谢谢老师的建议！ 关于第2条，我是想请问老师，Material Projects会不会有类似的AFM结构，可以用来参考磁矩，然后应用到我的输入文件? 如果有的话，请问老师该如何搜索？ 关于第4条，请问老师，用超晶格算的话，在CP2K中的输入关键词是什么？ |
不会扣篮的后卫 发表于 2022-2-16 12:43 1、设置的应该没问题，但操作上太繁琐，事实上所有spin up的Fe和spin down的Fe分别额外设置这两种即可，其余原子要么不是d价层要么是d满层，基本上不会带局域磁矩，因此不进行初猜也没关系，自洽迭代后应该能得到正确结果，MP网站上给出的结果带很少的磁矩是因为布居分析的时候划分到上面的。 2、我没明白你所谓的找到更多磁矩是什么意思。 3、最好不要去修改元素的价电子数目吧，而且那个也不是改价电子数目，而是进行constraint DFT来约束SCF中原子的布居电荷。 4、那个Decrease in energy我检查下来你体系的能量是一直上升的，仔细一看其实SCF也没收敛，那个多k点的事情是没办法避免，要么就直接用超晶格算，否则根本不可能得到合理的结果。另外你的DFT+U的参数也没正确设置上去，即便用了多k点因为（半）局域交换-关联泛函自作用的问题也会导致低估磁矩。你体系若不牵扯到要优化块体（虽然我还是很建议那么做一下），那直接去算大一些的表面再通过工具或者手工换算到小的表面也没什么。 另注：以后发邮件截图还是慎重一些，避免带个人信息，免得造成不必要的麻烦。 |
卡开发发 发表于 2022-2-14 15:21 好的，谢谢老师的建议！ 因为我老板也是刚接触CP2K，可能他也需要时间来把理论与CP2K完美结合。 我会转告他来自于论坛上老师，同学讲解的CP2K新知识。 |
| 参与人数Participants 1 | eV +2 | 收起 理由Reason |
|---|---|---|
|
| + 2 | 欢迎讨论 |
不会扣篮的后卫 发表于 2022-2-14 13:01 后面那个警告 WARNING in s_square_methods.F:99 :: Found non-uniform occupation 但是，需要注意的是RELAX_MULTIP_1 WARNING in qs_mo_occupation.F:357 :: Multiplicity changed from 1 to 103 |
卡开发发 发表于 2022-2-13 16:35 谢谢老师的建议。因为最终会造大表面，所以小结构就先不考虑k点问题了。 老师，我根据您的建议，采取了6个chemcial unit CuFeO2的结构，用CP2K跑了两个新的输入文件， 一个自旋多重度参数设置为是RELAX_MULTIP 1，想要让磁矩弛豫，由于在计算过程中发现了报错，我终止这个运行过程， 出现了以下类似的报错： *** WARNING in qs_mo_occupation.F:357 :: Multiplicity changed from 141 to *** *** 139 ----- 报错1 另一个自选多重度参数设置为是MULTIPLICITY 1，成功跑出输出文件，但出现一个 *** WARNING in s_square_methods.F:99 :: Found non-uniform occupation ***-----报错2 ------------------------------------------------------------------------------------------------------------------------------------- 我把FIXED_MAGNETIC_MOMENT的这个参数设置加入3个chemcial unit CuFeO2结构的输入文件中， （MULTIPLICITY 16），因为把3个Fe 都当成是+5的磁矩，MULTIPLICITY =2* （15/2) + 1 =16, 想测试是否可以把 *** WARNING in s_square_methods.F:99 :: Found non-uniform occupation ***避免, 没有FIXED_MAGNETIC_MOMENT 15的这个参数设置的话， *** WARNING in s_square_methods.F:99 :: Found non-uniform occupation ***又出现了。 但是有FIXED_MAGNETIC_MOMENT 15的这个参数设置的话， *** WARNING in s_square_methods.F:99 :: Found non-uniform occupation ***就不会出现。 请问老师,报错1，报错2为什么会出现呢？ （附上输入和输出文件) |
292.54 KB, 下载次数 Times of downloads: 6
本帖最后由 卡开发发 于 2022-2-13 16:44 编辑 不会扣篮的后卫 发表于 2022-2-13 14:02 1、61个电子占据64轨道，那说明有电子数占据是<1的，这种引起的原因是使用了展宽，也合理。 2、这个体系建议用RELAX_MULTIP替代MULTIPLICITY让磁矩弛豫，一般最终的磁矩可能不是正好是5。 3、AFM的结构应该不会是3 chemical unit，这样应该不符合化学计量数，应该是某个偶数倍，spin up和spin down的电子数目应当是一样的才是符合AFM的。需要小心的是如果你拿到的结构是实验结构，这个结构可能是不区分磁性的，因此实际考虑磁性区分后的结构可能比这个结构大2N倍。我不确定你做的是表面还是块体，如果是表面也得注意这个问题，最好先找到AFM下考虑区分磁矩的单胞结构再来做表面。 4、小的表面也应当考虑多k点，但是我不确定现在的版本DFT+U能否同时设置k点（n年没用过这程序，某些设计一直认为比较离谱），如果不能使用多k点你还得造大的表面。 5、SCF没收敛，虽然程序继续跑了下去  ，建议从上面的问题先下手。 |
卡开发发 发表于 2022-2-5 11:42 好的，谢谢老师！ |
不会扣篮的后卫 发表于 2022-2-5 11:10 占据数目指的就是电子在轨道上占据了多少，按道理是一个轨道占据1电子（不考虑自旋极化是2），由于使用了展宽，在Fermi能级附近的轨道占据数可能是分数，Fermi能级以上没有电子填充的轨道就是空轨道（或者叫空带）。 一般大体系计算做Fock矩阵对角化都是做子空间的，只需要求解比较低的几个占据态即可。对于无磁性的绝缘体，电子全部配对，需要的展宽也不大，这种情况对空轨道需求比较少；反之则需要稍多一些的空轨道。如果要增加空轨道，可以用ADDED_MOS（https://manual.cp2k.org/trunk/CP ... /SCF.html#ADDED_MOS）。 |
卡开发发 发表于 2022-1-29 16:50 请问老师，空带数目和占据数指的是什么？在跑3个chemical unit的CuFeO2的输入结构文件时，遇到了All alpha MOs are occupied的报错。 感觉和之前的老师提到的“自旋极化之后alpha或beta上可能会排更多的电子数”现象类似。 （附上3个chemical unit的CuFeO2输入和输出文件） |
4.23 KB, 下载次数 Times of downloads: 2
20.6 KB, 下载次数 Times of downloads: 1
All alpha MOs are occupied的报错
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2024-11-24 16:26 , Processed in 0.181191 second(s), 27 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites