小木虫 --- 500万硕博科研人员喜爱的学术科研平台
&&查看话题
计算声子谱时老是出错。求大神解答
如题，计算声子谱就会出现下面的问题
我是在优化后的晶胞基础上计算的。
错误提示如下：
*ERROR* in PHONON calculation:&&Prerequisite GEOMETRY OPTIMISATION& && & ?
?& &&&has not been performed or is not converged to sufficient accuracy.& && & ?
?& && && && && && && && && && && && && && && && && && && && && && && && && && &?
?& && &&&Maximum force |F|max& &&&= 0.1276& &&&eV/A& && && && && && && && && & ?
?& && &&&Parameter GEOM_FORCE_TOL = 0.5000E-01 eV/A& && && && && && && && && & ?
?& && && && && && && && && && && && && && && && && && && && && && && && && && &?
?& &This could be because of any of the following& && && && && && && && && && &?
?& &&&(a) No geometry optimisation was performed& && && && && && && && && && & ?
?& &&&(b) A geometry optimisation failed to converge& && && && && && && && && &?
?& &&&(c) An error was made copying the optimised geometry to this .cell file&&?
?& &&&(d) Cutoff, k-point sampling or other convergence parameter differ& && & ?
?& && && &from the geometry optimisation run to this one.& && && && && && && & ?
?& &Please check your input parameters and if necessary run or re-run a& && &&&?
?& &geometry optimisation calculation.& && && && && && && && && && && && && &&&?
?& && && && && && && && && && && && && && && && && && && && && && && && && && &?
?& &Alternatively increase value of GEOM_FORCE_TOL to 0.1276& &&&or greater& & ?
?& &to allow this run to proceed at risk of inaccurate phonon frequencies.& &&&?
?& && && && && && && && && && && && && && && && && && && && && && && && && && &?
?& &Typically a force tolerance convergence of 0.025 eV/A is necessary to& && &?
?& &obtain well converged phonon frequencies althout a smaller value can be& & ?
?& &required in the presence of weak forces, e.g. in molecular crystals& && &&&?
声子谱.jpg
我已经优化过了，还是这样的，
phonopy是单独的一个计算声子谱的软件还是要通过其他软件使用的一个模块啊？
http://phonopy.sourceforge.net/这是主页，你看看就知道了
是的，我确定
好的，我看看
浏览器进程
打开微信扫一扫
随时随地聊科研小木虫 --- 500万硕博科研人员喜爱的学术科研平台
&&查看话题
请教声子谱的计算
老板让自己摸索用vasp算声子谱，不知用vasp算声子谱怎么样，难不难学？想向算过的同学请教一下！希望同学们多多发表意见，我也好好学习一下:D
你好，我最近在算声子的时候一直被一个问题所烦恼，就是关于R3c的对称性它的超胞该怎么设，刚开始我们设的2*2*2，然后放上去用VASP跑，一跑就掉下来了，后来我们发现这样不对，要算六边形的区域，但是这样我就不知道要怎么设超胞了，请问，您有做过这样的对称性，该怎么设超胞，或者有没有相关的资料，在此谢过！！
一跑就掉下来，有可能是内存不够。如果硬件不够，建议使用Castep, Abinit, or Quantum Espresso.
一跑就掉下来应该是内存不足的问题，要算六边形的区域是什么意思？声子谱只能体系整体的啊
嗯，检查了下是因为内存不够好像，但是我们的hexagonal下如何设超胞呢，好像我们之前设的2*2*2是不对的，谢谢赐教！
昨天发现是因为内存不够好像，不过要解决的主要问题是我们的hexagonal下如何设超胞呢，这个因为是初学，怎么弄都没怎么搞懂，我们之前设的2*2*2是不对的，hexagonal下好像不能这么设，cubic是可以的，谢谢赐教！
你可以尝试一下只在C轴上超胞，a,b不动
具体我也不懂
我们主要动的就是a、b，C轴我们反而无所谓，我们现在如果a*2,b*2，那ab面就是一个斜的大“田”字形，但是这样是不对的，我们只要中间的六边形，所以不知道如何设超胞，关键是ab面，谢谢！
浏览器进程
打开微信扫一扫
随时随地聊科研小木虫 --- 500万硕博科研人员喜爱的学术科研平台
&&查看话题
【求助】abinit计算声子频率的问题
急求一个abinit计算声子频率的详细资料！本人对abinit软件一窍不通！所以力求详尽，能立竿见影，对本人能有直接迅速的帮助！谢谢！:rol:;)
浏览器进程
打开微信扫一扫
随时随地聊科研() update: 现在这个教程已经严重过时，已经不适合最新的 phonopy 版本了。
1. Phonopy 简介
Phonopy 是一个由 python 实现的的晶体声子分析程序。它是目前提供了 VASP 的 Wien2k 的接口用来计算原子受力。它的主要功能有：
计算声子色散谱；
计算声子态密度，包括分立态密度；
声子热力学性质，包括自由能，热容量，焓；
Phonopy 通过力常数的方法计算声子谱。力常数由计算原子在超晶胞中被移动后的受力得到(Parlinsk-Li-Kawasoe 方法)。
同样类型的程序还有 phon, fropho, phonon. 其中 phonon 是商业软件，卖的很贵，fropho 和 phonopy 的代码其实都是来自于 phon, fropho 是为了代替 phon 而开发的，目的是为了使用 phon 更方便，phon 和 fropho 主要都是由 fortran 开发的，而现在 fropho 已经停止开发，由 python 开发的 phonopy 代替了 fropho, phonopy 在使用上更为方便，在计算量上更为减少。因为 phon, fropho 和 phonon 在移动原子位置时都是一次只移动一个原子的一个方向，而 phonopy 则可以一次移动一个原子的多个方向，所以和其它程序相比, phonopy 最多可以减少 2/3 的计算量。
2. Phonopy 的安装
phonopy 主要由 python 开发，其中的找晶体空间群功能是由 C 语言(spglib)实现。在安装 phonopy 时需要安装以下它所依赖的软件包：
Python 的开发包
numpy , numpy-dev
matplotlib
python-lxml
python-yaml
如果是使用的 Ubuntu/Debian 系统只需：
sudo aptitude install python-dev python-numpy
python-matplotlib python-tk python-lxml python-yaml
如果是 OpenSUSE 则可以去 http://software.opensuse.org/search 搜索下载安装相关软件包。
去 https://sourceforge.net/projects/phonopy/ 下载 phonopy 的最新软件包.
解压后进入 软件包的 c 目录，
python setup.py install
编译安装所需要的 python 模块。
将 PATH/phonopy-0.6.1 $PATH 里：
比如： export PATH=~/opt/phonopy-0.6.1:$PATH
3. Phonopy 的简单使用
下面我以 BC5 为例说明如何和 VASP 做接口计算声子谱：
a, 建立超胞
准备 POSCAR 文件
单胞和原胞都可以
准备 INPHON 文件
ATOM_NAME = B C
NDIM = 2 2 1
LSUPER = .TRUE.
简单说明：
第一行为原子符号名 B C ，注意要和 POSCAR 中的顺序相同。
第二行为要建的超胞的维数。
第三行是指要建立移动原子的 POSCAR.
然后输入命令：
phonopy.py
就会输出 DISP, POSCAR-* (001, 002, 003, …)
在本例出会输出 12 个 POSCAR-*, 依次为 POSCAR-001,…POSCAR-012
b, 计算力常数
将这些 POSCAR-* 做 SCF, 将输出文件 vasprun.xml 保存下来。
我的计算脚本为：
for a in 01 02 03 04 05 06 07 08 09 10 11 12
cp POSCAR-0$a POSCAR
/opt/intel/impi/3.1/bin64/mpiexec -n 8 /share/apps/vasp/bin/vaspmpi_mkl10_O1 &&out.vasp 2&&err.vasp &
cp vasprun.xml vasprun.xml-0$a
我的 INCAR 为：
SYSTEM = BC5
PREC = Accurate
ENCUT = 520
EDIFF = 1e-5
ISMEAR = 0
SIGMA = 0.05
LCHARG = FALSE
LWAVE = FALSE
这样计算结束后会得到 vasprun.xml-001, …, vasprun.xml-012
然后，执行命令：
phonopy.py -f vasprun.xml-*
就会得到文件 FORCES
c, 得到声子谱
有了文件 FORCES 后就可以计算声子谱了，这时将 INPHON 改为以下内容：
ATOM_NAME = B C
NDIM = 2 2 1
NPOINTS = 51
QI = 0.0 0.0 0.5 0.0 0.5 0.5 0.0 0.5 0.0 0.0 0.0 0.0 0.5 0.5 0.0 0.5 0.0 0.0 0.5 0.0 0.5
QF = 0.0 0.5 0.5 0.0 0.5 0.0 0.0 0.0 0.0 0.5 0.5 0.0 0.5 0.0 0.0 0.5 0.0 0.5 0.5 0.5 0.5
简单说明：
ND 是指有几条路径；
NPOINTS 是指在每两个高对称性点间插入多少个点，包括边界；
QI 和 QF 给出路径的起始和结束点。
其中我用的 Q 点顺序为：
0.0 0.0 0.5 Z
0.0 0.5 0.5 T
0.0 0.5 0.0 Y
0.0 0.0 0.0 G
0.5 0.5 0.0 S
0.5 0.0 0.0 X
0.5 0.0 0.5 U
0.5 0.5 0.5 R
phonopy.py -p
即可画出声子谱， 如果执行 phonopy.py -p -s 即可将声子谱曲线保存为 PDF 文件。
其中 声子谱的数据是以 yaml 格式保存在 band.yaml 里的。
如果你想用 Origin, Qtiplot, Xmgrace 等工具画图的话，可以用我写的 python 脚本，(http://www.zhuli.name/file/readband.py) 将 band.yaml 转换一下，输出文件有两个 band.dat, qpoint.dat.
其中 band.dat 是声子谱数据， qpoint.dat 是高对称性点的数据。
以上只是 phonopy 的其本用法，如果要计算声子态度度和热力学性质的话可以参考
http://phonopy.sourceforge.net 上的文档，还是比较简单的。
我的计算结果:
Filed under:
Thank you! I am a student of College of Science , SCAU .
Leave a Reply114网址导航}