Dynaphopy 计算高温声子谱

李江山

你好，请问我得到了quasiparticles_data.yaml文件如何画出声子频率，因为我到了画图那一步报错了，只生成了这个文件，可以手动作图吗

Iridium-191

sobereva 发表于 2023-5-9 23:51
CUTOFF足够大是消除虚频的极其关键要素，绝不仅仅是看优化收敛精度

谐振近似得到的频率若有明显虚频总有 ...

老师好，那对于比如铁电体或者STO这种量子顺电体来说，或者其他高温相，室温相在0K下都是非稳态的相，这种是不是只能通过有限温度下的计算来解决了？谢谢老师

小付小氟

将拟合方式改为fftw,通过指令dynaphopy input OUTCAR -i -psm 3,采用快速傅里叶变换拟合。我想问用端口之间输出了图片。输入什么样的指令导出图片对应的数值呢？使用这种方法拟合，是对应的力常数文件发生改变吗

Shana

小付小氟 发表于 2024-4-8 16:08
将拟合方式改为fftw,通过指令dynaphopy input OUTCAR -i -psm 3,采用快速傅里叶变换拟合。我想问用端口之间 ...

能够直接输出力常数的

九点一定起

您好，我想问下那个input文件中的晶胞与原胞的转换矩阵是怎么得到的呀，如何得到这个矩阵呢

Shana

九点一定起 发表于 2024-6-12 11:03
您好，我想问下那个input文件中的晶胞与原胞的转换矩阵是怎么得到的呀，如何得到这个矩阵呢

你可以将晶胞转换成原胞，然后用原胞扩胞算声子和AIMD。这样转换矩阵就是对角矩阵

九点一定起

Shana 发表于 2024-6-12 15:14
你可以将晶胞转换成原胞，然后用原胞扩胞算声子和AIMD。这样转换矩阵就是对角矩阵

您好 我这个声子和AIMD都算完了 但是这个矩阵具体是怎么来的呢 在哪里看呢

九点一定起

九点一定起 发表于 2024-6-18 11:11
您好 我这个声子和AIMD都算完了 但是这个矩阵具体是怎么来的呢 在哪里看呢

就是dynaphopy 中input中需要输入的转换矩阵 这些我都算完了 不知道这个怎么设置 是怎么算出来的还是怎么得到的呢 麻烦您了

Shana

九点一定起 发表于 2024-6-18 12:37
就是dynaphopy 中input中需要输入的转换矩阵 这些我都算完了 不知道这个怎么设置 是怎么算出来的还是怎么 ...

有个PRIMITIVE MATRIX和SUPERCELL MATRIX，比如你将晶胞先约化成原胞，那么PRIMITIVE MATRIX的对角元就是1，1，1.然后2*2*2扩胞，SUPERCELL MATRIX的对角元就是2，2，2
PRIMITIVE MATRIX
1 0 0
0 1 0
0 0 1
SUPERCELL MATRIX
2 0 0
0 2 0
0 0 2

Shana

九点一定起 发表于 2024-6-18 12:37
就是dynaphopy 中input中需要输入的转换矩阵 这些我都算完了 不知道这个怎么设置 是怎么算出来的还是怎么 ...

phonopy是可以输出PRIMITIVE MATRIX的，可以参考输出界面

17873772769

Shana老师，您好，在dynaphopy的参数设置中写了使用fftw方法时，psm可设置为2、3、4，我在使用dynaphopy input OUTCAR -i -psm 3时，会存在报错“ModuleNotFoundError: No module named 'pyfftw'”，但我使用2时可以正常运行，这结果应该是一样的么（使用后确实无虚频了，且看起来结果合理）

Shana

17873772769 发表于 2024-10-21 17:07
Shana老师，您好，在dynaphopy的参数设置中写了使用fftw方法时，psm可设置为2、3、4，我在使用dynaphopy in ...

从我的使用感受上来看，dynaphopy的几个算法处理的结果是不大一样的，无虚频的话，消除的很成功