3 Degree of agreement with the experiment should never be simply taken as the sole criterion for selecting the functional. One should consider all possible sources of discrepancy with experimental observation (e.g. conformation, solvation effect, protonation state, vibronic effect, double-excitation, and so on)
Yes, I am aware of that but right now referring solely to the excitation which is probably the most relevant at the moment of discussing whether O3LYP is fine?
]]>2 If you can afford, certainly aug-cc-pVTZ is much better choice than aug-cc-VDZ. aug-cc-pVDZ is just the lowest acceptable basis set for this purpose.
3 Degree of agreement with the experiment should never be simply taken as the sole criterion for selecting the functional. One should consider all possible sources of discrepancy with experimental observation (e.g. conformation, solvation effect, protonation state, vibronic effect, double-excitation, and so on).
]]>To calculate electron affinity, diffuse functions are absolutely needed.
That's why I switched from the regular def2-TZVP to ma-def2-TZVP. Apparently, on another forum, people are arguing that Karlsruhe's basis sets are developed in a way that they don’t require diffuse functions as much as, for example, Pople’s basis sets. As many options exist, there are as many opinions... However, to avoid receiving negative comments from reviewers, I decided to stick with an augmented basis set.
aug-cc-pVDZ is not expensive. Dunning's basis sets are not preferential choice for DFT calculation of thermodynamic data, but for electron excitation studies, it is usable.
Okay, but shouldn’t I actually use aug-cc-pVTZ? I read somewhere that "For DFT calculations, the aug-cc-pVnZ basis set family may not be the best choice (non-ideal contraction for DFT, overly large and often poor SCF energies)". My system is also relatively small, with no more than 30 heavy atoms.
O3LYP is a rarely used functional, I do not recommend it. Generally, for calculating electronic absorption spectrum, DFT functionals with HF composition lower than 20% are not recommended, because excitation energies are usually notably underestimated. O3LYP only contains 11.61% HF composition.
I am aware that O3LYP is uncommon, but only it and TPSSh yield reasonable excitation energies compared to other functionals—no more than 10 nm greater than experimental results. I tested a bunch of them.
]]>aug-cc-pVDZ is not expensive. Dunning's basis sets are not preferential choice for DFT calculation of thermodynamic data, but for electron excitation studies, it is usable.
O3LYP is a rarely used functional, I do not recommend it. Generally, for calculating electronic absorption spectrum, DFT functionals with HF composition lower than 20% are not recommended, because excitation energies are usually notably underestimated. O3LYP only contains 11.61% HF composition.
]]>In the meantime, I checked some simpler functional - O3LYP instead of TPSSh and as for now looks fine. Given the absorption spectrum are coherent with the experimental one, I suppose the swap is validated?
Best,
]]>Best,
Tian
]]>! LibXC(TPSSh) ma-def2-TZVP AutoAux TightSCF D4 SMD(H2O)
Thanks
------------------------
DAVIDSON-DIAGONALIZATION
------------------------
Dimension of the eigenvalue problem ... 34500
Number of roots to be determined ... 40
Maximum size of the expansion space ... 400
Maximum number of iterations ... 100
Convergence tolerance for the residual ... 2.500e-07
Convergence tolerance for the energies ... 2.500e-07
Orthogonality tolerance ... 1.000e-14
Level Shift ... 0.000e+00
Constructing the preconditioner ... o.k.
Building the initial guess ... o.k.
Number of trial vectors determined ... 400
****Iteration 0****
Time for iteration : TOTAL=69.2 TRAFO=1.2 RIJ=6.8 COSX=27.6 XC=21.4
Size of expansion space: 120
Lowest Energy : -94740.194830773849
Maximum Energy change : 94740.194830773849 (vector 0)
Maximum residual norm : 1023356330656.517578125000