使用Amesp 2.1(dev)进行TDDFT-ris计算激发态结果诡异问题

Suchent

如题,本人使用Amesp 2.1版本在TD(ris)-wB97x-D3/def2-TZVP(-f)级别下计算三个分子(Benzophenone,Cinnamic acid,Para-aminobenzoic acid)的前5个单重激发态的时候发现结果非常诡异, 其能量和正常的紫外-可见激发相比明显过于小了(s1<1eV), 甚至有一个是负值. 同样的分子使用ORCA在相同级别下(普通TDDFT)计算则结果正常(s1在3-4eV之间). 我想请问一下这是TDDFT-ris近似的问题,或者是Amesp自身在计算TDDFT-ris时候的bug,又或者我关键词写错了?三个分子的Amesp输入/输出文件(.aip,.aop)+同级别下ORCA输出文件(.out)已经提供

JohnCase

有orca的mkl文件不？我给你续算一下ris试试

zjxitcc

补充一下2L所说，你现在应该是有Amesp的.mo文件，运行

1. m2a xxx.mo

复制代码

即可生成xxx.amo文件。随后运行

1. amo2fch xxx.amo

复制代码

即可生成xxx.fch文件，内含基态SCF轨道，也可以提供fch文件给2L做TDDFT-ris计算。所有量化软件都是互通的，没必要多次重复算基态。m2a是Amesp自带的小程序，amo2fch是MOKIT的小程序。

Suchent

JohnCase 发表于 2025-7-23 18:30
有orca的mkl文件不？我给你续算一下ris试试

久等了,mkl文件已经提供,AMESP计算的.mo文件我转成molden了也在里面

JohnCase

Suchent 发表于 2025-7-23 22:48
久等了,mkl文件已经提供,AMESP计算的.mo文件我转成molden了也在里面

根据orca的波函数，我对mol0用gpu4pyscf的TDDFT-ris计算前10个state，结果是这样的
eV        oscillator strength
3.74131 0.00103
5.12597 0.01713
5.19882 0.02539
5.42537 0.39693
5.53196 0.05450
6.40542 0.04584
6.50319 0.04692
6.50801 0.17981
6.64748 0.36612
6.72457 0.01717

UV-vis图谱看着跟orca TDDFT结果差不多的。




我上传了gpu4pyscf的计算结果。

你也可以用cpu版的代码试试，详细说明在这里
http://bbs.keinsci.com/thread-50702-1-1.html

嫌安装麻烦的话，我最近会推出在线计算TDDFT-ris的网站，请期待。