求助Dalton报错QC-SCF not implemented/tested for DFT

Voland

如题，最近我在用Dalton进行磷光光谱计算，这里根据Sob老师的意见，先用Gaussian进行Tn的几何优化，然后再通过Dalton计算振子强度等数据。但是过程中遇到了QC-SCF not implemented/tested for DFT的报错，看起来是QC-SCF目前并不支持DFT方法。那么，在此之前首先是SCF的DIIS不收敛的问题，我看了下手册，结合Sob老师在Gaussian中提到的能够解决SCF不收敛的建议，调整了.SHIFT（目前设置为-4.0），初猜结构也是Gaussian的DFT方法激发态几何优化收敛所得（泛函B3LYP，我注意到了Dalton这里的差别，不过这里不是重点），但是仍然感觉收敛非常缓慢，无意义地增加SCF ITER也不是个好办法（默认60）。请问一下还有什么办法可以加快收敛吗？
实在没辙那我只能用TDHF之类的方法做了，也不知道能不能行。
十分希望收到大家的建议，非常感谢！

hebrewsnabla

轨道读取看起来没什么问题，这一部分已经成功了。

报一个这样的错误Sorry, plot of MOs with Molden is only implemented for spherical GTOs，是我轨道信息读取的不对吗

不是报错。

zhouyulin

hebrewsnabla 发表于 2025-6-29 22:58
已经完全在瞎写了……看dalton手册5.2.4了解一下基本语法

老师，我这边试着改了一下，已经通过dalton得到的磷光计算速率，就是值不太对，以及报一个这样的错误Sorry, plot of MOs with Molden is only implemented for spherical GTOs，是我轨道信息读取的不对吗

zhouyulin

hebrewsnabla 发表于 2025-6-29 22:58
已经完全在瞎写了……看dalton手册5.2.4了解一下基本语法

对不起老师，我之前都是用BDF+MOMAP算的，没有用dalton算过磷光，所以实在是不太清楚，我看dalton里面算磷光是下面这样算的，但是我不知道要怎么把通过MOKIT转化得到的dal文件与dalton示例中计算磷光的dal文件结合起来，我再去尝试一下
ATOMBASIS
HF
test
Atomtypes=2 Angstrom
Charge=1.0  Atoms=1      Basis=aug-cc-pvdz
H          0.0000000000        0.0000000000       -0.5000000000
Charge=9.0  Atoms=1      Basis=aug-cc-pvdz
F          0.0000000000        0.0000000000        0.4500000000

**DALTON INPUT
.RUN RESPONS
**INTEGRALS
.SPIN-ORBIT
**WAVE FUNCTIONS
.DFT
B3LYP
*SCF INPUT
.DOUBLY OCCUPIED
  3 1 1 0
**RESPONS
*QUADRATIC
.PHOSPHORESENCE
.PRINT
3
.ROOTS
  1 1 1 0
**END OF DALTON INPUT

hebrewsnabla

zhouyulin 发表于 2025-6-29 21:15
改成这种形式吗，老师。但是我试了之后还是报错
**DALTON INPUT
.RUN WAVE FUNCTIONS

已经完全在瞎写了……看dalton手册5.2.4了解一下基本语法

zhouyulin

hebrewsnabla 发表于 2025-6-29 18:42
mokit生成的

中间不要插入别的东西，你把**RESPONS插进去，下面的orbital input就无效了。但是可以把H ...

改成这种形式吗，老师。但是我试了之后还是报错
**DALTON INPUT
.RUN WAVE FUNCTIONS
**WAVE FUNCTIONS
.DFT
*ORBITAL INPUT
.MOSTART
FORM18
.PUNCHOUTPUTORBITALS
.RUN RESPONS
**WAVE FUNCTIONS
.DFT
B3LYP
**RESPONS
*QUADRATIC
.PHOSPHORESENCE
.PRINT
3
.ROOTS
  1
**MOLORB (punched by fch2dal of MOKIT Sat Jun 28 17:33:02 2025)
-0.00000031056001 -0.00000212511394  0.00006060471390  0.00000120133812
  0.00000060866362 -0.00000034696842  0.00003191794410  0.00000054888932
  0.00000145203578  0.00000095409512 -0.00000022474850  0.00000040323304
-0.00000029556427 -0.00000009210741 -0.00000029050104 -0.00000005080526

hebrewsnabla

zhouyulin 发表于 2025-6-29 14:22
还是不行，老师。运行之后显示有报错
@ MPI MASTER, node no.:     0
@ Reason: Illegal "*" option un ...

mokit生成的

1. 
   
2. **WAVE FUNCTIONS
   
3. .HF
   
4. *ORBITAL INPUT
   
5. .MOSTART
   
6. FORM18
   
7. .PUNCHOUTPUTORBITALS
   

复制代码

中间不要插入别的东西，你把**RESPONS插进去，下面的orbital input就无效了。但是可以把HF改成DFT。
**RESPONS 相关的关键词写后面。

zhouyulin

hebrewsnabla 发表于 2025-6-29 13:32
试试

还是不行，老师。运行之后显示有报错
@ MPI MASTER, node no.:     0
@ Reason: Illegal "*" option under **RESPONS
以及
WARNING for ERROR  (DALTON program stopped with exit code 100)
- THEREFORE scratch directory /public1/home/scfa8196/dalton/dalton/install/tmp/DALTON_scratch_scfa8196/T1new_127316 is not deleted
  by this script so you may restart the calculation, if possible
  ( /public1/home/scfa8196/dalton/dalton/install/tmp/DALTON_scratch_scfa8196/T1new_127316 may be deleted by your queueing system )
- REMEMBER to delete /public1/home/scfa8196/dalton/dalton/install/tmp/DALTON_scratch_scfa8196/T1new_127316 yourself !!!

Output from 'grep -n ERROR' :
31379:  --- SEVERE ERROR, PROGRAM WILL BE ABORTED ---
----------------------------------------------------------
2 informational messages have been issued by DALTON,
output from 'grep -n INFO'  (max 10 lines):
30200:    INFORMATION: Switched to ".FCKTRA" integral transformation, because MPI parallel
30643:@ INFO: Sorry, plot of MOs with Molden is only implemented for spherical GTOs

hebrewsnabla

zhouyulin 发表于 2025-6-28 21:43
改成这样吗，老师，这是改完的dal文件(T1new.dall)

试试

zhouyulin

hebrewsnabla 发表于 2025-6-28 21:01
不要设置 .DOUBLY OCCUPIED

另外，当前的计算是c1的，所以roots什么的也是一个数，不是四个数。

改成这样吗，老师，这是改完的dal文件(T1new.dal[url=]l[/url])

hebrewsnabla

zhouyulin 发表于 2025-6-28 19:19
老师，这是我的输入输出文件

不要设置 .DOUBLY OCCUPIED

另外，当前的计算是c1的，所以roots什么的也是一个数，不是四个数。

zhouyulin

hebrewsnabla 发表于 2025-6-28 18:40
在dalton里面要做b3lyp的话，生成fch的高斯任务也要用相同的泛函。如果确实已经用了相同的泛函，上传一下 ...

老师，这是我的输入输出文件

hebrewsnabla

zhouyulin 发表于 2025-6-28 17:58
老师，我已经试着通过MOKIT得到相关的mol文件和dal文件，并且跑一遍之后正常结束。但是当我想把dal文件改 ...

在dalton里面要做b3lyp的话，生成fch的高斯任务也要用相同的泛函。如果确实已经用了相同的泛函，上传一下dalton输出文件。

zhouyulin

hebrewsnabla 发表于 2025-6-27 22:55
可以。

通过MOKIT得到相关的mol文件和dal文件，并且跑一遍之后正常结束，并且也看到了SCF收敛提示，但是我不知道我要怎么改dal文件，从而使它能够通过dalton计算磷光

zhouyulin

hebrewsnabla 发表于 2025-6-27 22:55
可以。

老师，我已经试着通过MOKIT得到相关的mol文件和dal文件，并且跑一遍之后正常结束。但是当我想把dal文件改成计算磷光的文件，即T1new(MOKIT).dal，运行时又出现了QC-SCF not implemented/tested for DFT的报错，是我的dal文件哪里有问题吗