用高斯计算电离出一个质子的有机分子的紫外光谱与实际光谱对不上

bychenpku

用不同泛函（B3LYP,PEB0,M06-2X）计算了3,5二羟基苯甲酸（DHBA）在不同PH下的紫外光谱。酸性条件下，DHBA不电离，这时候用B3LYP算的光谱和实验光谱拟合的最好，其他两个泛函都出现了蓝移。中性条件下，DHBA的羧基电离出质子，实验光谱是240nm和296nm，这时候这几个泛函算出来的紫外光谱全都不同程度的蓝移了，B3LYP蓝移的最少，但是还是差了有20nm了。溶剂模型是SMD隐式溶剂模型，溶剂是水。是因为氢键的原因还是泛函选的不好？
下面是Input
%chk=DHBZopt.chk
%mem=1000MB
%nprocs=10
# opt freq B3LYP/def2tzvp scrf=smd

DHBZopt

-1 1
C                 -1.30644111   -1.26235777   -0.03555408
C                  0.09255937   -1.18062969   -0.04188818
C                  0.72148020    0.07126302   -0.00804772
C                 -0.04859945    1.24142766    0.03212702
C                 -1.44759994    1.15969961    0.03846030
C                 -2.07652077   -0.09219310    0.00461977
H                  0.68053228   -2.07407637   -0.07256211
H                  0.43159555    2.19727558    0.05796329
H                 -3.14468869   -0.15459433    0.00945678
O                 -2.23339551    2.35374519    0.07945397
H                 -2.41492893    2.64825773   -0.81603918
O                 -1.94819706   -2.53979934   -0.07008444
H                 -2.08611522   -2.85297362    0.82685505
C                  2.25884339    0.16107403   -0.01500620
O                  2.82358862    1.28522260    0.01538083
O                  3.04463897   -1.03297159   -0.05599837


%chk=DHBZtd.chk
%mem=1000MB
%nprocs=20
# b3lyp/def2tzvp td(nstates=50) scrf=smd

DHBZtd

-1 1
C                 -1.29503700    1.19735900   -0.00002600
C                  0.09535400    1.20690400   -0.00002300
C                  0.79367700   -0.00000400   -0.00000300
C                  0.09535600   -1.20690400    0.00002100
C                 -1.29504300   -1.19735900    0.00001900
C                 -2.00413600   -0.00000300   -0.00000300
H                  0.62432200    2.14898500   -0.00005600
H                  0.62431900   -2.14898700    0.00005400
H                 -3.08878900    0.00000500    0.00000500
O                 -1.93962500   -2.41345600    0.00004600
H                 -2.89687900   -2.27509400    0.00001500
O                 -1.93963500    2.41345400   -0.00002700
H                 -2.89688700    2.27508100   -0.00008700
C                  2.31850600    0.00000100    0.00000000
O                  2.90100100   -1.11829700   -0.00014000
O                  2.90099000    1.11830500    0.00014200

granvia

试试CAM-B3LYP

sobereva

20nm的误差对于TDDFT来说已经可以接受了
阴离子，记得带弥散函数
还有其它角度可以试图改进结果，诸如特定态溶剂模型、考虑显式水、用双杂化泛函做TDDFT、用w调控的泛函等

PS：用CAM-B3LYP只会蓝移更厉害

bychenpku

sobereva 发表于 2018-9-8 03:45
20nm的误差对于TDDFT来说已经可以接受了
阴离子，记得带弥散函数
还有其它角度可以试图改进结果，诸如特 ...

好的，谢谢sob老师