计算化学公社

标题: 用高斯计算电离出一个质子的有机分子的紫外光谱与实际光谱对不上 [打印本页]

作者
Author: bychenpku 时间: 2018-9-7 10:53
标题: 用高斯计算电离出一个质子的有机分子的紫外光谱与实际光谱对不上
用不同泛函（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

作者
Author: granvia 时间: 2018-9-7 12:23
试试CAM-B3LYP

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

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

作者
Author: bychenpku 时间: 2018-9-8 20:05

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

好的，谢谢sob老师

欢迎光临计算化学公社 (http://bbs.keinsci.com/)