GaAs超胞自洽能Siesta计算输入文件正确性问询

整体可行性判断 #

你的输入文件可以完成受扰动GaAs超胞的自洽能计算，核心的SCF循环、基组、交换关联泛函、晶格/原子坐标等关键设置均已包含。下面针对细节给出修正建议和补充说明：

关键设置检查与修正 #

1. 移除冗余动力学设置 #

你设置了MD.TypeOfRun CG和MD.NumCGsteps 50，但MD.Steps 0会让Siesta不执行任何结构优化或动力学步骤，这部分属于冗余设置，建议删除以下两行以避免混淆：

MD.TypeOfRun            CG          # Type of dynamics:
MD.NumCGsteps  50

2. SCF收敛判据合理性说明 #

你同时设置了自由能（SCF.FreeE.Tolerance 1.d-6 eV）和密度矩阵（SCF.DM.Tolerance 1.d-4 eV）的收敛判据，Siesta会自动采用更严格的标准。对于自洽能计算，这个组合是合理的；若后续需要高精度结构优化，可考虑将密度矩阵判据收紧至1.d-5 eV。

3. 基组与计算精度补充 #

PAO.SplitNorm 0.26是合理的默认值，若后续收敛测试中结果不稳定，可尝试在0.2-0.3范围内调整该参数。DZP基组足以支撑GaAs的自洽能计算，后续收敛测试可验证是否需要升级至TZDP基组。

4. k网格设置建议 #

当前4×4×4的Monkhorst-Pack网格是16原子超胞的合适起点，后续收敛测试可逐步提升k点密度（如6×6×6），直到自洽能变化小于1e-6 eV。

5. 其他细节优化 #

• LatticeConstant 1.00 Ang设置正确，因你已用埃单位给出完整晶格矢量，该参数仅作为单位缩放因子。
• SaveHS .true.和SaveRho .true.会保存哈密顿量、重叠矩阵和电荷密度，若无需后续调用这些文件，可改为.false.节省磁盘空间。
• ElectronicTemperature 300. K有助于SCF收敛，对半导体GaAs体系无负面影响。

修正后的输入文件（移除冗余项） #

SystemName      Bulk GaAs
SystemLabel     GaAs

PAO.BasisSize           DZP
PAO.SplitNorm           0.26

#  >>>> KPOINTS <<<<

%block kgrid_Monkhorst_Pack
   4  0  0  0.0
   0  4  0  0.0
   0  0  4  0.0
%endblock kgrid_Monkhorst_Pack

xc.functional      GGA    # Default vaxc.authors         
xc.authors         PBE    # Default value

# >>>> self-consistent-field stuff <<<<

ElectronicTemperature   300. K
SCF.Mix                 Hamiltonian
SCF.Mix.First          .true.
SCF.Mix.First.Force    .false.
SCF.Mixer.Method        Pulay
SCF.Mixer.Weight        0.05 # 0.1
SCF.Mixer.History       6

MaxSCFIteration         1000
SCF.FreeE.Converge     .true.
SCF.FreeE.Tolerance     1.d-6 eV
SCF.DM.Converge        .true.
SCF.DM.Tolerance        1.d-4 eV
SaveHS                 .true.
SaveRho                .true.

Mesh.Cutoff             300. Ry     # Reduce eggbox effect

# >>>> Eigenvalue problem: order-N or diagonalization <<<<

SolutionMethod          diagon      # OrderN or Diagon

#
# Definition of Atomic Species
#
NumberOfSpecies       2
%block ChemicalSpeciesLabel
  1  31  Ga
  2  33  As
%endblock ChemicalSpeciesLabel

#
# Atoms
#
NumberOfAtoms         16

AtomicCoordinatesFormat   Ang # Ang|Bohr|Fractional
%block AtomicCoordinatesAndAtomicSpecies
  15.9714956   9.2251079   6.5243845   1
   9.9805742   1.1624815   3.2515290   1
   9.9998782   3.4597186   0.0189113   1
   4.0158845   4.6226940   3.2679597   1
  11.9955181   9.2120818   6.5270022   1
  10.0131978   8.0744755   3.2672024   1
   9.9843919   5.7806334   6.5138485   1
   7.9993333   4.6273260   3.2797882   1
   1.9906870   1.1415551   0.8025453   2
   4.0015415   2.3195738   4.0641599   2
   4.0091927   4.5998543   0.8183037   2
   5.9923684   5.7572029   4.0649693   2
   6.0124279   1.1630422   0.8234672   2
   7.9795257   2.2963281   4.0748852   2
   8.0148156   4.6141221   0.8180230   2
   9.9968920   5.7599471   4.0969952   2
%endblock AtomicCoordinatesAndAtomicSpecies

LatticeConstant      1.00 Ang
%block LatticeVectors
       7.9950000   0.0000000   0.0000000
       3.9975000   6.9238730   0.0000000
       3.9975000   2.3079580   6.5278900
%endblock LatticeVectors

内容的提问来源于stack exchange，提问作者user29730484