|
|
谢谢老师的回复!Multiwfn果然很方便!太感谢啦!
我还想请老师帮忙评估一下,我目前计算分子晶体的精度是否可以接受?还需要进一步提高吗?
我把 vc-relax.in 文件,scf.in 文件,以及ph.in 文件列在下面了,请老师帮忙评估一下哈,万分感谢!
vc-relax.in文件
&CONTROL
title='CBP-A-vc-relax',
calculation='vc-relax',
pseudo_dir='./',
outdir='./tmp',
verbosity='high',
tprnfor=.true.,
tstress=.true.,
forc_conv_thr=1.0d-4,
nstep=300,
/
&SYSTEM
ibrav= 0, nat= 124, ntyp= 3,
occupations = 'fixed',
ecutwfc = 60,
ecutrho = 600,
input_dft= 'pz',
vdw_corr = 'grimme-d2',
/
&ELECTRONS
conv_thr = 1.0d-8
mixing_beta = 0.5d0
/
&IONS
/
&CELL
cell_dynamics = 'bfgs'
press_conv_thr=0.1
/
ATOMIC_SPECIES
C 12.0107 C_ONCV_PBE_sr.upf
H 1.00794 H_ONCV_PBE_sr.upf
N 14.00674 N_ONCV_PBE_sr.upf
scf.in文件
&CONTROL
title='CBP-A-scf',
calculation='scf',
prefix = 'CBP-A-Raman'
pseudo_dir='./',
outdir='./tmp',
verbosity='high',
/
&SYSTEM
ibrav= 0, nat= 124, ntyp= 3,
occupations = 'fixed',
ecutwfc = 80,
ecutrho = 800,
input_dft='pz',
vdw_corr = 'grimme-d2',
/
&ELECTRONS
conv_thr = 1.0d-10
mixing_beta = 0.5d0
/
&IONS
/
&CELL
/
ATOMIC_SPECIES
C 12.0107 C_ONCV_PBE_sr.upf
H 1.00794 H_ONCV_PBE_sr.upf
N 14.00674 N_ONCV_PBE_sr.upf
ph.in 文件
&INPUTPH
tr2_ph = 1.0d-14
prefix = 'CBP-A-Raman'
outdir = './tmp'
fildyn = 'CBP-A-Raman.dyn'
ldisp = .false.
epsil = .true.
trans = .true.
lraman= .true.
qplot = .false.
nq1 = 1, nq2 = 1, nq3 = 1
/
0.0 0.0 0.0
|
|