- 博客(31)
- 收藏
- 关注
RSS订阅原创 QE:symmetry tolerance
当你用相同对称性结构却得到不同的动力学矩阵数，多半是symmetry tolerance的问题，QE默认的symmetry tolerance是1e-5,这可能导致他有的对称性找不到。这个值不能在输入文件中改。改完，用完记得改回去。在这里改，改完重新编译。
2024-06-06 21:37:25
117
原创 ABINIT-OPT计算输入输出浅析
重要参数： scalecart1，scalecart2 ， scalecart3体积扩大收缩相关。a.SrB3C3.files-所有要用到的文件目录及tmp目录。这个文件我是通过DFT-CP软件包自动准备的，也可以自己准备。给的是每个原子的WK位置，有几个不同的WK位置给几行。ngfft：这个网格要测试收敛，且是12的倍数。每个参数都可以在这个网址查到具体意思。材料相关，原子类型，个数，能带数。化学压强计算相关参数。
2024-05-01 14:28:05
518
原创 VASP OUTCAR中的各种能量
enthalpy is TOTEN：焓值H=U+PV，同时等于吉布斯自由能：G=U+PV。free energy TOTEN:亥姆霍兹自由能A。energy without entropy： 内能U。这是正常的定义，在VASP的计算中，T=0。吉布斯自由能：G=U-TS+PV。亥姆霍茨自由能：A=U-TS。吉布斯自由能：G=U+PV。亥姆霍茨自由能：A=U。
2024-04-01 15:08:12
1869
1
原创 画conveshull的技术整理
x轴，y轴导入origin之后，手动算一下，cah4的值和cah2加cah6合成cah4的值，哪个低。所以，CaH4是m=1,n=1,横坐标是0。比如这个图，x轴：cah2是0，H2是1。cah4的焓值减cah2+h2/每原子。y轴 ：cah2是0，H2是0。横坐标的值是n/(m+n)
2024-03-10 19:41:02
452
原创 固定金属相空间群及原子位置
有的时候我们在制作晶体结构的时候经常出现文献中值给相的名称，例如fcc,bcc.hcp但不给具体的空间群及不等价原子的WK位置的情况，尤其是老文献。下面是个人整理的不同名称代表的晶体空间群。如有出入，欢迎指正。
2023-11-11 15:54:08
1616
原创 使用EPW计算孤立ME方程
1.自洽声子计算scf.in &control calculation = 'scf' restart_mode = 'from_scratch', prefix = 'pb', pseudo_dir = '../', outdir = './' / &system ibrav = 2, celldm(1) = 9.27, nat = 1, n
2022-05-27 18:02:31
3350
原创 QE动力学矩阵文件的主要内容及单位
matdyn0不可约q点的位置matdyn1PHonon/PH/io_dyn_mat_old.f90 写文件部分代码Mould/cell_base.f90Mould/ions_base.f90 用到的变量及单位基本形式Dynamical matrix fileElectron-phonon coefficients for hydride 2 4 0 1.8897269 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000
2022-05-15 16:18:43
843
4
原创 手动计算高对称路径的方法两则
0：获得高对称点详见文章“QE计算声子谱”第二节1：直接写脚本import numpy as np#每两个高对称点之间插入多少个点number = 30#scf.out里面的倒格失vec_x=np.array([0.000000, 0.322290, 0.322290])vec_y=np.array([0.322290, 0.000000, 0.322290])vec_z=np.array([0.322290, 0.322290, 0.000000])#高对称点的个数num
2022-05-13 20:50:10
705
原创 QE画声子谱
1 q2r.x建立q2r.in&INPUTfildyn = 'matdyn' lshift_q = .true. zasr = 'crystal' flfrc = 'ifc.dat'/q2r.x -in q2r.in > q2r.out2.计算高对称点2.1方法1：MS1.改元胞2.计算高对称点2.2方法2：网站AFLOWhttp://www.aflowlib.org/aflow-online/2.2.1.输入POSCAR2.2.2.勾选
2022-05-09 21:38:16
2440
2
原创 SCTK计算lambda
在SCTK的后处理中，calculation = “lambda_mu_k”可以计算出lambda和mu的值PATH/sctk.x < sctk.in > lambda.outlambda.out文件##### Integrate gapeq ##### mu : 0.36011E+00 lambda : 0.18219E+01 ##### Write mu.frmsf and lambda_frmsf #####
2022-03-08 16:04:15
595
原创 用DOCKER安装及使用SSCHA
# 添加docker用户组，一般已存在，不需要执行sudo groupadd docker# 将登陆用户加入到docker用户组中sudo gpasswd -a $USER docker# 更新用户组newgrp docker# 测试docker命令是否可以使用sudo正常使用docker version#增加权限sudo chmod a+rw /var/run/docker.sock......
2022-02-23 21:52:32
901
原创 SCTK的使用——MgB2
1. 自洽计算——算电子密度export OMP_NUM_THREADS=1 #只允许单线程$ mpiexec -np 29 PATH/pw.x -nk 29 -in scf.in > scf.out #PATH可执行文件路径重要的参量：calculation = “scf”用pw.x进行自洽计算注意：我们可以先计算声子和电声相互作用也可以先计算库伦作用。（大概是说第二步和第三步的计算顺序无所谓？）2. 计算声子和电声相互作用2.1 计算声子频率和变形势（deformation po
2022-01-13 16:10:29
1418
原创 使用QE中的ld1.x制作超软赝势
使用QE中的ld1.x制作超软赝势输入文件inputrlderiv计算对数导数的半径config电子结构当“lsd”=1时，向上自旋和向下自旋状态可能出现两次，各自占位为:3p4 3p2 = 4up，2down。否则，假定采用洪德规则|| | |...
2021-08-17 20:52:51
1172
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人