2、ThiRDoRDer计算三阶力常数 首先需要对高精度优化后的结构扩胞: 与phonopy不同的是参数-d表示考虑哪些近邻原子的受力来计算力常数矩阵。d为正数时表示截断半径(单位nm) 为负数时表示所考虑近邻原子的个数。 一般情况d取<=-10 或 >=0.6。 这会影响到扩胞生成的POSCAR-*的个数与计算精准度。 这里使用有限位...
〇需要准备的软件:VASP+Thirdorder+ShengBTE ①需要准备INCAR\KPOINT\POSCAR\POTCAR四个输入文件,新建一个目录名为shengbte,放进去,这四个文件应该尽量和用phonopy计算二阶力常数矩阵的时候一致。 INCAR可以写为: SYSTEM=HELLO #随便写 PREC=High ISTART=0 ICHARG=2 ISPIN=1 NCORE=4 NELM= 60 NELMIN=4 ISYM=...
2、ThiRDoRDer计算三阶力常数 首先需要对高精度优化后的结构扩胞: 与phonopy不同的是参数-d表示考虑哪些近邻原子的受力来计算力常数矩阵。d为正数时表示截断半径(单位nm) 为负数时表示所考虑近邻原子的个数。 一般情况d取<=-10 或 >=0.6。 这会影响到扩胞生成的POSCAR-*的个数与计算精准度。 这里使用有限位...
这样就能得到二阶力常数文件FORCE_CONSTANTS,以及声子谱的图band.pdf数据a.dat。对于声子谱,我们要保证没有虚频,这样才能保证晶格的稳定性。 2、ThiRDoRDer计算三阶力常数 首先需要对高精度优化后的结构扩胞: 与phonopy不同的是参数-d表示考虑哪些近邻原子的受力来计算力常数矩阵。d为正数时表示截断半径(单位nm) 为...
2、ThiRDoRDer计算三阶力常数 首先需要对高精度优化后的结构扩胞: 与phonopy不同的是参数-d表示考虑哪些近邻原子的受力来计算力常数矩阵。d为正数时表示截断半径(单位nm) 为负数时表示所考虑近邻原子的个数。 一般情况d取<=-10 或 >=0.6。 这会影响到扩胞生成的POSCAR-*的个数与计算精准度。
#扩胞命令模版 thirdorder_vasp.py sow x y z -c (x,y,z是三个方向的扩胞倍数,c是考虑多少个临近原子的受力来计算力常数矩阵) thirdorder_vasp.py sow 2 2 2 -12 根据显示结果这里有72个DFT任务,相对来说较少,但是如果手动一个任务一个任务提交也是非常麻烦且耗时,这里我们通过一个脚本来整理文件,并一...
天玑算-科研服务丨第一性原理 热电教程(5):VASP+Thirdorder计算三阶力常数矩阵 834 -- 16:21 App 天玑算 | 第一性原理:VASPKIT功能介绍11 Density ofStates(DFT+VASP) 1.4万 15 19:31 App 天玑算-科研服务丨第一性原理 热电教程(1):热电材料基本概述浏览...
#扩胞命令模版 thirdorder_vasp.py sow x y z -c (x,y,z是三个方向的扩胞倍数,c是考虑多少个临近原子的受力来计算力常数矩阵) thirdorder_vasp.py sow 2 2 2 -12 根据显示结果这里有72个DFT任务,相对来说较少,但是如果手动一个任务一个任务提交也是非常麻烦且耗时,这里我们通过一个脚本来整理文件,并一...