QE:利用ph.x 计算声子谱input 注意事项(包括单个q点的计算)
material phonon
&INPUTPH
prefix = 'A', !此处的a可以用物质代替比如我算mg的 就写成mg
outdir = '.',
tr2_ph=1.0d-12, !1.0d-12 默认值 自恰的阀值
fildyn = 'A.dyn' !default matdyn 但是我们会自己设置 这是是动力学矩阵写入的文件
fildvscf = 'dv', ! potential variation 写入的文件 随后会用到电声计算中 与fildrho有关
fildrho = 'drho', ! 电子密度的反馈文件,出存在${outdir}/_ph0/${prefix}.${fildrho} 中
ldisp = .true., !ph计算的时候声子的选择可以单独算也可以成格算 默认单独算 但是如果 ldisp=true 则需要设置nq1 nq2 nq3 并且在这个设置为 true之后对下面的lshift_q也会成立
lshift_q = .true.,
nq1 = 4,
nq2 = 4,
nq3 = 4,
第二种如果想单独算某个高对称点的声子 就是单个q点
则输入文件如下
phonons of Si at Gamma 算si的g点处的q&inputphtr2_ph=1.0d-14, 默认阀值是12 现在是14 更精确prefix='silicon',epsil=.true., 此处只实用于非金属并且q=0才可以是true (可能是G点) 而对于默认是false 如果是金属并且q不等于0 就一定要用默认值amass(1)=28.08, 如果不加这行 可以从默认文件中读入outdir='/public/home/lszhang/Softwares/qe-6.5/tempdir/',fildyn='si.dynG',/
0.0 0.0 0.0 !此处是G点
具体的参数设置参考https://www.quantum-espresso.org/Doc/INPUT_PH.html#nq1
ldisp | LOGICAL |
---|---|
Default: | .false. |
If .true. the run calculates phonons for a grid of q-points specified by nq1, nq2, nq3 - for direct calculation of the entire phonon dispersion. |
ldisp | LOGICAL |
---|---|
Default: | .false. |
If .true. the run calculates phonons for a grid of q-points specified by nq1, nq2, nq3 - for direct calculation of the entire phonon dispersion. |
nq1, nq2,nq3 | INTEGER |
---|---|
Default: | 0,0,0 |
Parameters of the Monkhorst-Pack grid (no offset) used when ldisp=.true. Same meaning as for nk1, nk2, nk3 in the input of pw.x. |
epsil | LOGICAL |
---|---|
Default: | .false. |
If .true. in a q=0 calculation for a non metal the macroscopic dielectric constant of the system is computed. Do not set epsil to .true. if you have a metallic system or q/=0: the code will complain and stop. |
QE:利用ph.x 计算声子谱input 注意事项(包括单个q点的计算)相关推荐
- vasp+phonopy计算声子谱及后处理
记录一下最近在计算的三维Na2Te的声子谱以及后处理得到部分声子态密度等信息.计算声子谱主要有两种方法:有限位移法和密度泛函微扰理论(DFPT).vasp version = 5.4.4; phono ...
- Keras之MLPR:利用MLPR算法(1to1+【Input(1)→8(relu)→O(mse)】)实现根据历史航空旅客数量数据集(时间序列数据)预测下月乘客数量问题
Keras之MLPR:利用MLPR算法(1to1+[Input(1)→8(relu)→O(mse)])实现根据历史航空旅客数量数据集(时间序列数据)预测下月乘客数量问题 目录 输出结果 设计思路 实现 ...
- matlab2018中变压器模块,利用MATLAB中Sim+Power+Systems模库时变压器模型的参数计算及其仿真结果比较...
[实例简介] 变压器模型 matlab 仿真 参数计算 第21卷第1期向秋风,等:利用 MATLAB中 Sim Power System模库时变压器模型的参数计算及其仿真结果比较 17 其标幺值:R= ...
- 腾讯云微计算实践:从Serverless说起,谈谈边缘计算的未来
2019独角兽企业重金招聘Python工程师标准>>> 欢迎大家前往云+社区,获取更多腾讯海量技术实践干货哦~ 作者:黄文俊,腾讯云高级产品经理,曾经历过企业级存储.企业级容器平台等 ...
- python 科学计算基础教程电子版-自学Python 编程基础、科学计算及数据分析
自学Python 编程基础.科学计算及数据分析 epub pdf mobi txt 下载 自学Python 编程基础.科学计算及数据分析 epub pdf mobi txt 下载 ☆☆☆☆☆ 李金 著 ...
- 中用BBP公式计算_【真课堂】7年级信息技术:数据计算
教师介绍: 汪小红,一级教师,滨江区科技节scratch项目优秀指导教师,曾获杭州市信息技术优质课二等奖. [教材分析] 数据计算是数据处理过程中的一个重要步骤,数据计算能对数据进行归纳和提炼.学生已 ...
- 输入圆的的半径r 和一个整型数k,当k=1 时,计算圆的面积;当k=2 时,计算圆的周长;当k = 3 时,既要求求圆的周长也要求出圆的面积
题目:输入圆的的半径r 和一个整型数k,当k=1 时,计算圆的面积:当k=2 时,计算圆的周 长:当k = 3 时,既要求求圆的周长也要求出圆的面积. 解析:输入半径r和整型数k,由题目分析的,由于 ...
- 测试方法介绍-计算模型复杂度(GMac)、模型大小(M)、计算速度(FPS)
PRNet-V 计算复杂度为 48.76GMac 参数数量为34.73M (PRNet测试结果)(IEO在12345层) 参数数量为27.57M (PRNet测试结果)(IEO在345层) 计算图片读 ...
- pq分解法潮流计算的matlab,第四节PQ分解法潮流计算
<第四节PQ分解法潮流计算>由会员分享,可在线阅读,更多相关<第四节PQ分解法潮流计算(17页珍藏版)>请在人人文库网上搜索. 1.第四节 PQ分解法潮流计算一 .PQ分解法的 ...
最新文章
- 浪潮集团执行总裁王柏华:这次人工智能产业发展大潮“是真的”
- python 逻辑回归
- Linux环境编程 用户层定时器使用一 timerfd的使用
- windows dos 常用命令
- win7乱码 字符集解决方案
- 控制台怎么查看错误的详细信息_Java 程序该怎么优化?命令篇
- 使用sphinx快速为你python注释生成API文档
- 宝元系统u盘使用说明_教你如何使用U盘安装电脑系统
- 2017.10.12 礼物(zjoi2011) 失败总结
- 1032. 挖掘机技术哪家强(20)-PAT乙级真题
- python newbie——PE No.3
- PL/SQL developer连接远程ORACLE
- Ruby 2.0 发布首个预览版
- revit2019 导出obj_Revit模型转obj模型插件
- python的xlwt模块_Python xlwt模块
- 打印机 打印机驱动 打印机如何与PC通信 什么通信协议 蓝牙打印机项目 蓝牙协议栈
- 基于Rook+Ceph的云原生存储架构剖析
- U盘在手,忘记任何密码都可找回!!
- Linux上安装VIM
- python远程操控超微服务器(IPMI)、电脑开机(WOL)