计算化学公社

标题: 求助：异质结构电荷密度的计算及后处理所对应的输入文件怎么写呢？ [打印本页]

作者
Author: ZhangBigPower 时间: 2023-8-30 15:36
标题: 求助：异质结构电荷密度的计算及后处理所对应的输入文件怎么写呢？
本帖最后由 ZhangBigPower 于 2023-8-30 16:15 编辑

学习过程中看到知乎的一个不详细的教程，设置了石墨4个C原子吸附一个O原子的示例。
他提供了自洽计算输入文件和一个pp后处理的输入文件。
最后需要导出：总体系、碳原子、氧原子三者的电荷密度xsf文件，教程只能用来得到总体系的xsf文件，请问怎么改编输入文件再分别得到两个分项的电荷密度xsf文件呢？
教程提供输入文件见下：
&CONTROL
  pseudo_dir = './'
  calculation = 'scf'
  prefix = 'grap-o'/&SYSTEM
  nat= 5
  ntyp= 2
  ibrav= 0
  ecutwfc= 20.0
/
&ELECTRONS
  mixing_beta = 0.7
  conv_thr =  1.0d-8
/
&IONS
/
&CELL
/
ATOMIC_SPECIES
C 12.011  C.pbe-n-kjpaw_psl.1.0.0.UPF
O 15.999  O.pbe-n-kjpaw_psl.1.0.0.UPF
CELL_PARAMETERS angstrom
   2.46400000       0.00000000       0.00000000
   0.00000000       4.26780000       0.00000000
   0.00000000       0.00000000    12.00000000
ATOMIC_POSITIONS angstrom
C       0.00332640       0.08834346       4.92792000
C       1.23453792       2.27618845       5.26944000
C       1.23520320       0.82829462       5.01600000
C       0.00280896       3.01511534       4.98420000
O       1.23103904       1.62522092       6.50928000
K_POINTS automatic
7 5 1  0 0 0
pp后处理输入文件
&inputpp
prefix = 'grap-o'
outdir = '.'
filplot = 'grap-o'
plot_num= 0
/
&plot
nfile = 1
filepp(1) = 'grap-o'
weight(1) = 1.0
iflag = 3
output_format = 5
fileout = 'grap-o.xsf'
/

作者
Author: ZhangBigPower 时间: 2024-1-18 16:09
已解决

欢迎光临计算化学公社 (http://bbs.keinsci.com/)