SCTK 屏蔽库仑相互作用

Question

使用SCTK软件计算屏蔽库伦相互作用报错
我修改了输入文件如下，但还是报错了
步骤 1： scf 计算，文件如下

&CONTROL
   calculation='scf', pseudo_dir='./', outdir='./tmp', verbosity='high',
  tprnfor=.true., tstress=.true., forc_conv_thr=1.0d-5, nstep=100, prefix='CeYH6'
/
&SYSTEM
  ibrav= 0, nat= 8, ntyp= 3,
  occupations = 'smearing', smearing = 'gauss', degauss = 1.0d-2,
  ecutwfc = 50, ecutrho = 500,
  la2F=.true.
/
&ELECTRONS
  conv_thr = 1.0d-8
  mixing_beta = 0.7d0
/
&IONS
ion_dynamics='bfgs'
/
&CELL
  press_conv_thr=0.1
  press=500
/
ATOMIC_SPECIES
  Ce 140.116 Ce.pbe-spdn-rrkjus_psl.1.0.0.UPF
  Y 88.90585 Y.pbe-spn-rrkjus_psl.1.0.0.UPF
  H 1.00794 H.pbe-rrkjus_psl.1.0.0.UPF
CELL_PARAMETERS (angstrom)
   3.7671000000  0.0000000000  0.0000000000
   0.0000000000  3.7671000000  0.0000000000
   0.0000000000  0.0000000000  3.7671000000
ATOMIC_POSITIONS (crystal)
  Ce  0.0000000000  0.0000000000  0.0000000000
   H  0.2518800000  0.5000000000  0.0000000000
   H  0.7481200000  0.5000000000  0.0000000000
   H  0.5000000000  0.0000000000  0.2518800000
   H  0.0000000000  0.2518800000  0.5000000000
   Y  0.5000000000  0.5000000000  0.5000000000
   H  0.0000000000  0.7481200000  0.5000000000
   H  0.5000000000  0.0000000000  0.7481200000
K_POINTS {automatic}
  6 6 6 0 0 0

步骤 2 bands.in 文件 使用 pw.x 计算后用脚本 处理 输出文件如下：

  tprnfor=.true., tstress=.true., forc_conv_thr=1.0d-5, nstep=100, prefix='CeYH6'
/
&SYSTEM
  ibrav= 0, nat= 8, ntyp= 3,
  occupations = 'smearing', smearing = 'gauss', degauss = 1.0d-2,
  ecutwfc = 50, ecutrho = 500,
  la2F=.true.
/
&ELECTRONS
  conv_thr = 1.0d-8
  mixing_beta = 0.7d0
/
&IONS
ion_dynamics='bfgs'
/
&CELL
  press_conv_thr=0.1
  press=500
/
ATOMIC_SPECIES
  Ce 140.116 Ce.pbe-spdn-rrkjus_psl.1.0.0.UPF
  Y 88.90585 Y.pbe-spn-rrkjus_psl.1.0.0.UPF
  H 1.00794 H.pbe-rrkjus_psl.1.0.0.UPF
CELL_PARAMETERS (angstrom)
   3.7671000000  0.0000000000  0.0000000000
   0.0000000000  3.7671000000  0.0000000000
   0.0000000000  0.0000000000  3.7671000000
ATOMIC_POSITIONS (crystal)
  Ce  0.0000000000  0.0000000000  0.0000000000
   H  0.2518800000  0.5000000000  0.0000000000
   H  0.7481200000  0.5000000000  0.0000000000
   H  0.5000000000  0.0000000000  0.2518800000
   H  0.0000000000  0.2518800000  0.5000000000
   Y  0.5000000000  0.5000000000  0.5000000000
   H  0.0000000000  0.7481200000  0.5000000000
   H  0.5000000000  0.0000000000  0.7481200000
K_POINTS {automatic}
  6 6 6 0 0 0
K_POINTS crystal
432
   0.000000000000000    0.000000000000000    0.000000000000000   1.0
   0.000000000000000    0.000000000000000    0.166666666666666   1.0
   0.000000000000000    0.000000000000000    0.333333333333333   1.0
   0.000000000000000    0.000000000000000   -0.500000000000000   1.0
   0.000000000000000    0.000000000000000   -0.333333333333334   1.0
   0.000000000000000    0.000000000000000   -0.166666666666667   1.0
   0.000000000000000    0.166666666666666    0.000000000000000   1.0
   0.000000000000000    0.166666666666666    0.166666666666666   1.0
   0.000000000000000    0.166666666666666    0.333333333333333   1.0
   0.000000000000000    0.166666666666666   -0.500000000000000   1.0
   0.000000000000000    0.166666666666666   -0.333333333333334   1.0
                                                                                                                

接着计算sctk.in 文件如下：

&CONTROL
  prefix = 'CeYH6'
  outdir = './tmp'
  calculation = 'kel'
/
&KEL
  start_q = 1
  last_q = 24
  nci = 5
  laddxc = 0
  ecutfock = 80
  nq1 = 6
  nq2 = 6
  nq3 = 6
/

产生报错如下：

  Number of symmetries :      24
     Lowest band which contains FS :      12
     Highest band which contains FS :      15
     [opt_tetra]  Optimized tetrahedron method is used.

     #####  Read from STDIN  #####

                      q grid : 6  6  6
         # of Chebyshev int. : 5
                 The first q : 1
                  The last q : 24
                      laddxc : 0
                         lsf : 0
       Cutoff kinetic energy :  0.80000E+02


     Dynamical matrices for ( 6, 6, 6)  uniform grid of q-points
     With a half shift
     (  11 q-points):
       N         xq(1)         xq(2)         xq(3)
       1   0.083333333   0.083333333   0.083333333
       2   0.083333333   0.083333333   0.250000000
       3   0.083333333   0.083333333   0.416666667
       4   0.083333333   0.250000000   0.250000000
       5   0.083333333   0.250000000   0.416666667
       6   0.083333333   0.416666667   0.250000000
       7   0.083333333   0.416666667   0.416666667
       8   0.250000000   0.250000000   0.250000000
       9   0.250000000   0.250000000   0.416666667
      10   0.250000000   0.416666667   0.416666667
      11   0.416666667   0.416666667   0.416666667
     Because shifted q grid is used, q2r will not work !

     EXX grid:     4385 G-vectors     FFT dimensions: (  24,  24,  24)
     Real-space augmentation: initializing EXX grid
     Message from routine real_space_q:
     integrated real space q too different from target q
     Largest difference:   0.475585 for atom type # 1

     #####  Read gvectors.dat & evc.dat  #####

       Total RAM for input WFC :   0.27E+00 GB
ik =     5, ii =     1
In    0.0000000000   0.1471033632   0.1471033632, Required    0.0000000000   0.0000000000  -0.2942067263

 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     Error in routine get_wfcg (5):
     Invalid k-point (order)
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

     stopping ...
"kel.out" 95L, 3526B                                                                      

Answer 1

阿里嘎多学长整理AIGC生成，因移动端显示问题导致当前答案未能完全显示，请使用PC端查看更加详细的解答过程

问题解答

您使用 SCTK 软件计算屏蔽库仑相互作用，但报错了。您已经修改了输入文件，但是问题仍然存在。

可能原因

1. 输入文件格式错误：请检查您的输入文件是否符合 SCTK 的格式要求。
2. 伪势文件路径错误：在 &CONTROL 部分，您设置了 pseudo_dir 为 ./，请确保该目录下存在正确的伪势文件。

解决方案

1. 检查输入文件格式：使用 SCTK 的文档或示例文件来检查输入文件的格式是否正确。
2. 检查伪势文件路径：确保 pseudo_dir 指定的目录下存在正确的伪势文件。
3. 重新运行 SCF 计算：如果以上步骤都没有解决问题，请重新运行 SCF 计算，检查输出文件是否正确。

核心代码

无需代码，问题主要是输入文件格式或伪势文件路径错误。