计算同位素如何设置

AutMaple

高阁 发表于 2025-4-25 20:46
(1)我用的phonopy算的声子谱，之前算想算β再算错差的很多，目前是通过pdos对Fe积分算。偏差不算离谱，但 ...

（7）北京高压的刘锦老师吗？你们组只有你一个做计算有点麻烦。。。
你去求求刘锦老师，说想去地化所的刘耘老师那里问一问怎么做这些计算。刘耘老师那里，我们有专门适合VASP计算β的脚本，直接处理vasprun.xml文件，基于公式（1）。关于高压的问题，你可以顺便问问地化所的张飞武老师。这些东西在网上交流起来不太方便，但是线下就方便得多。
（6）一系列的实验晶格参数和实验频率，和这些实验晶格参数用软件计算出的压力，两个压力能得到一个依赖关系，这个方法比较糙。细致一点的方法可以问问张飞武老师。
（5）没有晶格参数，就最好用别人推荐的U值
（4）好吧哈哈哈
（3）声学支？不是光学支吗？使用公式（1）计算的时候，带入的都是光学支的频率。原胞的声学支是3，2*2*2就是8个原胞，也就是24个声学支，也没问题呀。
（2）我看不太懂matlab的代码，我只能读懂python的，而且我想看的是你基于公式（1）给出的计算脚本
（1）先算β，再用β算<F>?是不是本末倒置了？

高阁

AutMaple 发表于 2025-4-25 23:16
（7）北京高压的刘锦老师吗？你们组只有你一个做计算有点麻烦。。。
你去求求刘锦老师，说想去地化所的 ...

(7)是高科的刘老师；
(6)感谢提醒；
(5)我看网上说U值依赖软件和赝势(不知道是否依赖体系、温压)，找到的文献也都用的QE。我只找到一篇文献推荐了U，使用后电子步很难收敛，之前有人提醒我这说明结构可能不合理。请问应该怎么使用推荐的U呢？
(3)不是声学支，是声子谱的条数。原胞里有n个原子就有3n条声子谱，用j标记。222超胞有8个波矢k，其它波矢k是由周期边界条件引入，非物理。巨正则系综配分函数是对量子态i连乘，i应该既包括k又包括j，我觉得应该把8组连乘起来开8次方以归一化，但是不太清楚怎么在高对称路径上把物理k找到。您说的加权这个是怎么找到物理k并确定权重呢？
(2)这个脚本就是根据(1)算的，可以简述一下：同位素声子谱计算用的就是您讲的方法。phonopy给的band.dat是从左到右再从右到左……这样作图，脚本里先给它顺序调整一下再排成一列。轻重元素都这样就排成了两列，我默认两个声子谱排列方式是对应的，再严格照(1)计算，频率单位取了hz。Γ点声学支会导致分母上3个0，就乘了1。之前尝试用(6)算，phonopy对pdos归一化后为3，需要除掉。请问用(1)算是否也要除以3呢？除以3后似乎很接近正确值。
(1)用公式(1)得到的是β，但是实验测到的是F，大家作图也是F，所以我想算个F好做比较。目前感受是算β再得F准一点，因为用了温度，用(6)算F不准，因为没加温度，应改用准简谐近似声子谱。
此外用(1)还是有些问题，我之前说的是[size=1em]https://doi.org/10.1016/j.gca.2020.11.025，它说β几乎与1/T^2成正比。(1)中β是0K的，用它算2000K下F应该还是有偏差。

感谢指导！

AutMaple

这个应该是u=hv/kbT吧，你没有带上温度吗？

高阁

AutMaple 发表于 2025-4-26 10:03
这个应该是u=hv/kbT吧，你没有带上温度吗？

有温度，是图里的分母上的2。不好意思，太久了忘记了。请问(1)频率单位是hz吗？用hz应该是这样，结果不需要除以3，是359 N/m.之前不知道单位改来改去乱了。(1)是无量纲量，我不太确定应该用什么单位。