Can Intel MKL make NWCHEM faster in fluorescence calculation?


Forum Vet
Can Intel MKL make NWCHEM faster in fluorescence calculation using TDDFT?


Very Best Regards?

Forum Vet
GAMESS TDDFT gives the excitation energy of 250.4 nm(b3lyp and aug-cc-pvtz), preoptimized using B3LYP and 6-311G++(2d,2p) of uracil in ethanol, and S1-S0 of 321 nm(b3lyp and 6-31G(d,p)) in water through geometry search, whereas Solvent Effects on the Steady-state Absorption and Fluorescence Spectra of Uracil, Thymine and 5-Fluorouracil gives experimental 259 nm in ethanol and 312 nm in water, respectively. The reference gives identical UV absorption peaks in a few solutions, and the calculated here is several nanometers different from the experimental one, within the UV chemical accuracy, and the fluorescence discrepancy here is also within it, both of which can be said “in good agreement with each other”.


NIST gives the UV peak of 258nm at pH3.0 or 7.0, from reference The Ultraviolet Absorption Spectra of Some Pyrimidines. Chemical Structure and the Effect of pH on the Position of Xmar.

Actually, this S1-S0 peak is very weak, as stated in Absorption and Fluorescence Spectra of Uracil in the Gas Phase and in Aqueous Solution: A TD-DFT Quantum Mechanical Study (the same corresponding author), giving calculated S1 and S2 peaks of 303 nm (very weak) and 269.5nm in water, respectively.

Very Best Regards!

Forum Vet
GAMESS calculations B3LYP RHF using accd all give 4.65eV and 5.13eV, and reference Absorption and Fluorescence Spectra of Uracil in the Gas Phase and in Aqueous Solution: A TD-DFT Quantum Mechanical Study gives TDDFT calculated gaseous n/pi 4.80eV and Pi/Pi 5.26eV.

GAMESS cr-eomccsd(T) with aug-cc-pvdz can give two gaseous peaks of 4.96 and 5.62 eVs, the lowest two among all kinds of the first two calculated, but perhaps not the most accurate ones.


According to the reference:Theoretical Study of the Electronic Spectra of Uracil and Thymine
"The present paper reports the results of a similar study on the pyrimidine bases uracil and thymine. The pi-pi* transitions lowest in energy are predicted to be located at 5.0,5.8, 6.5, and 7.0 eV for uracil. ...In addition, we locate four n-pi* transitions and find the first one to be close in energy to the lowest pi-pi* transition."


Very Best Regards!

Forum Vet
The input with a planar symmetry from QA test tce_uracil_creomact gives 5.15 and 5.49 eVs, respectively, uisng 6-31+G(d,p) by NWCHEM7.0.0. Pehaps the rechoice of active orbitals can improve the second one, because I arbitrarily choose them identical to those of the first one from the QA test.


Very Best Regards!

Forum Vet
The second excitation energy is perhaps an hard one for a lot of functionals.

Very Best Regards!


Forum >> NWChem's corner >> Running NWChem