vasp输入文件快捷转化为QE输入文件(pwscf.in)

2 篇文章 0 订阅
1 篇文章 0 订阅
文章讲述了如何将vasp的输入文件转换为QuantumESPRESSO的in文件进行DFT计算,重点在于控制参数、系统参数、电子参数和离子参数的设置,并强调了结构优化的任务。转换工具vaspkit简化了这一过程,但用户仍需确保伪势路径和文件名的准确性。
摘要由CSDN通过智能技术生成

经验分享

vasp软件的输入文件包括INCAR,KPOINTS,POSCAR和POTCAR,其中INCAR为计算的设置文件,POSCAR为计算模型的结构信息,KPOINTS为网格密度设置,POTCAR为赝势信息。而在Quantum ESPRESSO软件包中,对于不同的计算任务则需要写入不同的in文件,如以pw.x计算结构优化和自洽时,所需in文件的内容则等同于包括了INCAR、POSCAR和KPOINTS,同时也需要指定好相关赝势的路径和文件名称。

由于pw.x所需in文件内容较为复杂,手动编写麻烦,这里借助vaspkit,将已有vasp输入文件快速转化为in文件

在转化之前准备好POSCAR、KPOINTS文件

借助于vaspkit的416功能

图片

图片

如果没有准备KPOINTS则会报错

图片

生成的pwscf.In文件为模版内容,如有计算需要,则自行修改。

其中 pseudo_dir =和 ATOMIC_SPECIES位置在进行QE计算之前需要完全确定好。

 &CONTROL calculation = 'relax' title = 'pwscf' restart_mode = 'from_scratch' outdir = './' prefix = 'pwscf' wf_collect = .true. pseudo_dir = './' etot_conv_thr = 1.0D-6 forc_conv_thr = 1.0D-6 tstress = .true. tprnfor = .true. / &SYSTEM ibrav = 0 nat=          20 ntyp=           4 ecutwfc = 50 ecutrho = 400 occupations = 'smearing' degauss = 0.001 smearing = 'gaussian' / &ELECTRONS electron_maxstep = 250 conv_thr = 1.0D-10 diagonalization = 'cg' /  &IONS ion_dynamics = 'bfgs' pot_extrapolation = 'second_order' wfc_extrapolation = 'second_order' / &CELL cell_dynamics = 'bfgs' press_conv_thr = 0.1 / ATOMIC_SPECIES P  ATOMIC_MASS    P_PSEUdo S  ATOMIC_MASS    S_PSEUdo Cu ATOMIC_MASS   Cu_PSEUdo In ATOMIC_MASS   In_PSEUdo  K_POINTS { automatic }   5   5   2   0.000000   0.000000   0.000000  CELL_PARAMETERS (angstrom)    3.0713142534    5.3143103819    0.0000000000   -3.0713142534    5.3143103819    0.0000000000   -2.0331424593    0.0000000000   13.7129044246  ATOMIC_POSITIONS (crystal)  P    0.6120167299    0.8952915069    0.3337950929  P    0.5574443271    0.9513622495    0.1694800230  P    0.1047084931    0.3879832701    0.8337950929  P    0.0486377505    0.4425556729    0.6694800230  S    0.8701965644    0.9490888864    0.1305656946  S    0.9325699557    0.5448957766    0.3774602855  S    0.6580069895    0.1846603003    0.3789131289  S    0.2921182450    0.9100464945    0.3786027656  S    0.5295430290    0.6541199142    0.1292835356  S    0.2344427778    0.2906374998    0.1312045491  S    0.0509111136    0.1298034356    0.6305656946  S    0.4551042234    0.0674300443    0.8774602855  S    0.8153396997    0.3419930105    0.8789131289  S    0.0899535055    0.7078817550    0.8786027656  S    0.3458800858    0.4704569710    0.6292835356  S    0.7093625002    0.7655572222    0.6312045491 Cu    0.2902955674    0.5508846117    0.3647423517 Cu    0.4491153883    0.7097044326    0.8647423517 In    0.9123898369    0.2620767557    0.2361130171 In    0.7379232443    0.0876101631    0.7361130171

下面借用某中文智能语言模型对该in文件的解释:

  1. 这段代码是一个用于进行密度泛函理论(DFT)计算的输入文件,使用Quantum ESPRESSO软件包进行处理。以下是代码中各个部分的功能解释:

  2. &CONTROL: 控制参数部分,指定了一些计算的控制参数,如计算模式(calculation)、标题(title)、重启模式(restart_mode)、输出目录(outdir)、前缀(prefix)等。其中,计算模式被设置为'relax',表示将执行结构优化计算;重启模式被设置为'from_scratch',表示从原始输入文件重新开始计算;输出目录和前缀指定了输出文件和日志文件的位置和名称。

  3. &SYSTEM: 系统参数部分,指定了模拟系统的参数,如晶格参数(ibrav)、原子数量(nat)、原子类型(ntyp)、能量截断(ecutwfc和ecutrho)、占据数(occupations)、自洽迭代收敛阈值(conv_thr)等。其中,ibrav设置为0表示使用默认的晶格类型;nat设置为20表示模拟系统中包含20个原子;ecutwfc和ecutrho分别表示平面波基组和密度泛函的能量截断。

  4. &ELECTRONS: 电子参数部分,指定了电子计算的参数,如电子最大步数(electron_maxstep)、自洽迭代收敛阈值(conv_thr)、对角化方法(diagonalization)等。

  5. &IONS: 离子参数部分,指定了离子计算的参数,如离子动力学方法(ion_dynamics)、势能外推方法(pot_extrapolation)、波函数外推方法(wfc_extrapolation)等。

  6. &CELL: 细胞参数部分,指定了晶胞计算的参数,如晶胞动力学方法(cell_dynamics)、压力收敛阈值(press_conv_thr)等。

  7. ATOMIC_SPECIES: 原子种类部分,定义了模拟系统中的四种原子类型,包括元素名称、原子质量和对应的pseudopotential文件。

  8. K_POINTS: k点部分,定义了布里渊区的k点采样。

  9. CELL_PARAMETERS: 晶胞参数部分,定义了晶格参数。

  10. ATOMIC_POSITIONS: 原子位置部分,定义了模拟系统中每种原子的位置坐标。

  11. 这段代码主要的功能是进行密度泛函理论计算,并对一个包含20个原子的固体进行结构优化和性质计算。

愿有所成

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值