|
|
本帖最后由 高阁 于 2025-2-18 22:38 编辑
背景
我尝试用dynaphopy算Mg7O8Fe在2000K，50Gpa下的非简谐声子谱。目前仅在Fe上加了磁矩，在Fe的d轨道上加了U值，计算按照如下步骤进行：
1.线性响应计算U值。初始磁矩设为4，根据Calculate U for LSDA+U - VASP Wiki中方法算U。从这一步开始每一步加压。
2.结构优化；
3.phonopy用有限位移法计算0K时声子谱与FORCE_CONSTANTS。此处单点能计算时根据上一步更新磁矩，扩胞为1*2*2，共64个原子，晶格参数为7.96978、7.92590、7.92590；
4.用NVP系综训练MFLL。此处POTIM=0.7，训练了50000步。从这一步开始每一步加温；
5.用MFLL计算MD，共200000步，POTIM=0.7；
6.Dynaphopy计算非谐性声子谱。
问题
1.请问如何确定哪些物质需要设置磁矩以及初始磁矩的值？目前我只给3d以及更高轨道的物质设置了磁矩，Fe的3d轨道有6个电子，4个轨道未充满，请问设置为4是否正确？
2.请问如何确定哪些物质需要设置U值？目前我只给Fe的d轨道加了U值，请问是否正确？
3.请问训练MFLL时IVDW应该如何设置？Part 2: Machine learning force fields4.2中给出IVDW=10，但是和LDA泛函冲突，目前我使用了PBE泛函，请问是否可行？如果用LDA泛函应该如何设置？
4.请问训练MFLL时能否使用单胞，等MD时再扩胞？上述扩胞大小是否合理？
5.请问用NVP系综训练MFLL时应该如何设置LANGEVIN_GAMMA_L和LANGEVIN_GAMMA？根据线性响应和U值以及涨落、耗散的关系，我猜测LANGEVIN_GAMMA_L为涨落，应该为LANGEVIN_GAMMA（耗散）的一半左右，LANGEVIN_GAMMA具体值应该与原子质量成正比，且应该与U值大小差不多。目前我将O的LANGEVIN_GAMMA取为U值的初值，其余元素按质量取初值，LANGEVIN_GAMMA_L取为U值作为初始值再调整。请问这样做是否合理？
6.请问用MFLL算MD后应该如何将轨迹可视化？目前使用ovito和VMD都会报错。
7.请问在训练MFLL和使用MFLL进行MD计算时有少量电子步未收敛（估计<3%）是否需要重新计算？
8.请问dynaphopy算声子谱时-psm 3和 -psm 1算得结果不同，是否说明参数设置还存在问题？
感谢指导！
|
|
发表于 Post on 2025-2-18 16:06:52
|
只看该作者 Only view this author
|只看大图 View big picture
|正序浏览
|阅读模式 Reading model
收藏 Add to favorites
