关于做PDOS图所需的空轨道计算问题

Papillon

大家好，我想得到带有空轨道信息的.molden文件来做PDOS图，计算的是49个水分子液滴中的Fe-S体系。之前已经用OT的方法跑完MD，想基于其RESTART.wfn文件用对角化算单点能得到.molden文件，加入ADDED_MOS 100 100之后却提示非整数占据。


部分输入文件如下：
————————————————————————————————————————————————————
&SCF
      MAX_SCF 200
      EPS_SCF 1.0E-05 #Convergence threshold of density matrix of inner SCF
     SCF_GUESS RESTART #Use wavefunction from WFN_RESTART_FILE_NAME file as initial guess
      &DIAGONALIZATION
        ALGORITHM STANDARD #Algorithm for diagonalization
      &END DIAGONALIZATION
      &MIXING #How to mix old and new density matrices
        METHOD BROYDEN_MIXING #PULAY_MIXING is also a good alternative
        ALPHA 0.4 #Default. Mixing 40% of new density matrix with the old one
        NBROYDEN 8 #Default is 4. Number of previous steps stored for the actual mixing scheme
      &END MIXING
      ADDED_MOS   100   100 #Number of virtual MOs to solve for alpha and beta spins
      &PRINT
        &RESTART #Note: Use "&RESTART OFF" can prevent generating .wfn file
          BACKUP_COPIES 0 #Maximum number of backup copies of wfn file. 0 means never
        &END RESTART
      &END PRINT

    &END SCF


报错如下：
————————————————————————————————
*** WARNING in qs_scf.F:576 :: SCF run NOT converged ***




  Integrated absolute spin density  :                               5.2260420845


*** WARNING in s_square_methods.F:99 :: Found non-uniform occupation ***




*** WARNING in s_square_methods.F:101 :: Found non-uniform occupation ***


  Ideal and single determinant S**2 :                    8.750000       8.891584




请问各位老师这是什么原因呢？谢谢大家！

sobereva

先确保依然用OT做单点任务在读取wfn文件的情况下SCF一轮就能收敛，免得单点任务计算的关键设置和动力学过程不一样，或者有什么其它疏漏

Papillon

sobereva 发表于 2024-8-15 00:45
先确保依然用OT做单点任务在读取wfn文件的情况下SCF一轮就能收敛，免得单点任务计算的关键设置和动力学过程 ...

好的，谢谢sob老师，我去试一下