Can Dalton calculate oscillator strength

 Posts: 290
 Joined: 24 Sep 2014, 08:36
 First name(s): yan
 Last name(s): xiong
 Affiliation: CENTRAL CHINA NORMAL UNIVERSITY
 Country: China
Re: Can Dalton calculate oscillator strength
CAS 1185987443 has been optimized with 6311++G(2d,2p), and the excitation energies in hexane have been calculated with 631+G(d,p), all with GAMESS, finding three weak absoptions together can give the first experimetal peaks in an article published in China using Shimadzu LC一20AD where the middle one, with the largest oscilator strength, has an error of around 0.12 eVs, but the one of 204 nms caused by charge transfer fails.
CAMQTP01 in GAMESS, a functional modified for vertical ionization potential calculation, contributed last year by Dr. Gerasimov, in A.N.Nesmeyanov Institute of Organoelement Compounds of Russian Academy of Sciences, now is being tried.
The optimization using Dalton2018 with b3lyp and 631+G* has been finished, and 631++G(d,p) optimization with that also has been completed.
Very Best Regards!
CAMQTP01 in GAMESS, a functional modified for vertical ionization potential calculation, contributed last year by Dr. Gerasimov, in A.N.Nesmeyanov Institute of Organoelement Compounds of Russian Academy of Sciences, now is being tried.
The optimization using Dalton2018 with b3lyp and 631+G* has been finished, and 631++G(d,p) optimization with that also has been completed.
Very Best Regards!
Last edited by xiongyan21 on 04 Sep 2020, 08:30, edited 12 times in total.

 Posts: 290
 Joined: 24 Sep 2014, 08:36
 First name(s): yan
 Last name(s): xiong
 Affiliation: CENTRAL CHINA NORMAL UNIVERSITY
 Country: China
Re: Can Dalton calculate oscillator strength
During the augccpvdz and b3lyp calculation of the above using 6311++G(d,p) with b3lyp optimized geometry, at the beginnning of the DFT response employing Dalton2018, there is a warning
********* WARNING ORBDIA ***********
DIAGONAL ORBITAL HESSIAN HAS 1 NEGATIVE ELEMENTS
...
Perhaps .INA setting is not very suitable.
The result is completely wrong
...
@ Excited state no: 1 in symmetry 1 ( A )  singlet excitation

@ Excitation energy : 4.54909991E02 au
@ 1.2378731 eV; 9984.1202 cm1; 119.43660 kJ / mol
@ Total energy : 679.67677 au
@ Operator type: XDIPLEN
@ Oscillator strength (LENGTH) : 0.24231228 (Transition moment : 2.8266406 )
@ Operator type: YDIPLEN
@ Oscillator strength (LENGTH) : 2.23750730E02 (Transition moment : 0.85894454 )
@ Operator type: ZDIPLEN
@ Oscillator strength (LENGTH) : 2.78751923E03 (Transition moment : 0.30317391 )
...
Very Best Regards!
********* WARNING ORBDIA ***********
DIAGONAL ORBITAL HESSIAN HAS 1 NEGATIVE ELEMENTS
...
Perhaps .INA setting is not very suitable.
The result is completely wrong
...
@ Excited state no: 1 in symmetry 1 ( A )  singlet excitation

@ Excitation energy : 4.54909991E02 au
@ 1.2378731 eV; 9984.1202 cm1; 119.43660 kJ / mol
@ Total energy : 679.67677 au
@ Operator type: XDIPLEN
@ Oscillator strength (LENGTH) : 0.24231228 (Transition moment : 2.8266406 )
@ Operator type: YDIPLEN
@ Oscillator strength (LENGTH) : 2.23750730E02 (Transition moment : 0.85894454 )
@ Operator type: ZDIPLEN
@ Oscillator strength (LENGTH) : 2.78751923E03 (Transition moment : 0.30317391 )
...
Very Best Regards!

 Posts: 290
 Joined: 24 Sep 2014, 08:36
 First name(s): yan
 Last name(s): xiong
 Affiliation: CENTRAL CHINA NORMAL UNIVERSITY
 Country: China
Re: Can Dalton calculate oscillator strength
An article in China suggested that it has two maximum wavelengths at 204 and 265 nm in 70% hexane, and said the area of the first one is very large, using Shimadzu LC20AD.
CAMQTP01 gives the first one. the fouth, 0.06 eVs deviating from the experimental data, using 631++G(d,p), but it gives the initial two of very small oscillator strengthes and small trasition dipoles, meaning " dark", and the third one divating 0.79 eVs from experimental data.
Very Best Regards!
CAMQTP01 gives the first one. the fouth, 0.06 eVs deviating from the experimental data, using 631++G(d,p), but it gives the initial two of very small oscillator strengthes and small trasition dipoles, meaning " dark", and the third one divating 0.79 eVs from experimental data.
Very Best Regards!
Last edited by xiongyan21 on 04 Sep 2020, 08:25, edited 12 times in total.

 Posts: 290
 Joined: 24 Sep 2014, 08:36
 First name(s): yan
 Last name(s): xiong
 Affiliation: CENTRAL CHINA NORMAL UNIVERSITY
 Country: China
Re: Can Dalton calculate oscillator strength
CAS 864237810, a pesticide, is difficult to optimize, because it requires a lot of Monte Carlo search steps, and the final geometry is obtained using 6311++G(2d,2p) and B3lyp, all with GAMESS. PCM in water, 631++G(d,p) and b3lyp excitation energies can be compared with the expeimental data in an article published in China. According to the calculation, the second absorption can give the data which can be compared with the experimental data because the first one is very weak. The experimental data shows one terrace and one slope nearby, and the calculation gives five absorptions around there where three of them are very weak, but it fails to predict the peak in 200 nms. The functional used for ionization potential fails for the excitation energy. If possible, I will also use Dalton2018 to study whether the two softwares can agree.
Very Best Regards!
Very Best Regards!
Last edited by xiongyan21 on 22 Aug 2020, 01:44, edited 12 times in total.

 Posts: 290
 Joined: 24 Sep 2014, 08:36
 First name(s): yan
 Last name(s): xiong
 Affiliation: CENTRAL CHINA NORMAL UNIVERSITY
 Country: China
Re: Can Dalton calculate oscillator strength
It should be noted that if the small basis set is used, the excited energies calculated near 200 nms may not be compared with experimental data in terms of eV., and perhaps even using augccpvdz and any functional, TDDFT results may also fail in the VUV and other regions in terms of eV.
Very Best Regards!
Very Best Regards!
Last edited by xiongyan21 on 21 Aug 2020, 11:08, edited 12 times in total.

 Posts: 290
 Joined: 24 Sep 2014, 08:36
 First name(s): yan
 Last name(s): xiong
 Affiliation: CENTRAL CHINA NORMAL UNIVERSITY
 Country: China
Re: Can Dalton calculate oscillator strength
Are there methods in Dalton 2018 to treat double excitations and above in the existence of solvation?
I don‘t know whether Prof. Truhlar's new method, the combination of MCSCF and DFT together, which may always overcome the difficulties of KSDFT in multiconfigurational problems, transition metal chemistry, charge transfer, etc., or MCSCF in GAMESS can include solvation, because I haven't tried this.
Prof. Truhlar's group has demonstrated its abilities in several pubilished articles, e.g,, excitation spectra of retinal in vacuum, which has experimental data.
Very Best regards!
I don‘t know whether Prof. Truhlar's new method, the combination of MCSCF and DFT together, which may always overcome the difficulties of KSDFT in multiconfigurational problems, transition metal chemistry, charge transfer, etc., or MCSCF in GAMESS can include solvation, because I haven't tried this.
Prof. Truhlar's group has demonstrated its abilities in several pubilished articles, e.g,, excitation spectra of retinal in vacuum, which has experimental data.
Very Best regards!
Last edited by xiongyan21 on 18 Aug 2020, 00:45, edited 1 time in total.

 Posts: 290
 Joined: 24 Sep 2014, 08:36
 First name(s): yan
 Last name(s): xiong
 Affiliation: CENTRAL CHINA NORMAL UNIVERSITY
 Country: China
Re: Can Dalton calculate oscillator strength
TDCIS in GAMESS using 631++G(d,p) predicts the oscillator strength of the first and the fourth excitation energies are very small, and the transition dipole of the z direaction and the oscilator strength of the fourth is only the half of the first, and the excitation energy of the first is far from the experimental data, but the fifth, which can be compared with the experiment data(error: 0.27 eVs and augccpvdz: smaller) if the order is correct, and sixth, not shown completely in the experiment, are very strong. TDCIS fails to give the relatively large oscilator strength of the first absorption, the terrance as well as the sople.
It seems the fifth one is of charge transfer character, and perhaps of multiconfiguration nature, thus whether TDCIS is useful for this is questionable., but b3lyp can reproduce the peaks above 200 nms.
I wonder whether Dalton2018 can do TDCISD.
Very Best Regards!
It seems the fifth one is of charge transfer character, and perhaps of multiconfiguration nature, thus whether TDCIS is useful for this is questionable., but b3lyp can reproduce the peaks above 200 nms.
I wonder whether Dalton2018 can do TDCISD.
Very Best Regards!
Last edited by xiongyan21 on 22 Sep 2020, 03:38, edited 20 times in total.

 Posts: 395
 Joined: 27 Jun 2013, 18:44
 First name(s): Hans Jørgen
 Middle name(s): Aagaard
 Last name(s): Jensen
 Affiliation: Universith of Southern Denmark
 Country: Denmark
Re: Can Dalton calculate oscillator strength
As long as you can specify the desired CISD with the RAS specifications, and it fits in memory and does not exceed some internal limits, then you can do a TDCISD.
To do that you ask for .CI in the wavefunction part, and then under **RESPONSE the keywords for excitation calculation and oscillator strengths.
To do that you ask for .CI in the wavefunction part, and then under **RESPONSE the keywords for excitation calculation and oscillator strengths.

 Posts: 290
 Joined: 24 Sep 2014, 08:36
 First name(s): yan
 Last name(s): xiong
 Affiliation: CENTRAL CHINA NORMAL UNIVERSITY
 Country: China
Re: Can Dalton calculate oscillator strength
Dear Prof. Jensen
Thanks a lot for your explanation.
I first try TDDFT in TDA with Dalton2018, know TDCIS should be done also with the specification of .CI, and am going to try TDCISD.
By the way, spinflip TDDFT geometry search with GAMESS using b3lyp and 631G(d,p) of phenyl methyl ketone can give T1S0 vertical gap in ethanol, where the ground state targeted is not spincontaminated at all, within 0.06 eVs deviating from the 00 band summarized by Prof. Wagner et al. previously in Communications to the Editor in JACS, and there are two imaginary frequencies during the following Hessian analysis below 380 cm^1. The imaginary frequencies perhaps should be eliminated. Perhaps I should follow the test using spinflip TDDFT instead of DFT to do Hessian, so the Hessian analysis here is not proper for ZPVE.
I will use Dalton2018 to calculate the phosphorescence lifetime.
I don't know whether Monte carlo minimum search, illustrating Monte Carlo AI, can be applied first to facilitate an triplet geometry optimization.
Very Best regards!
Thanks a lot for your explanation.
I first try TDDFT in TDA with Dalton2018, know TDCIS should be done also with the specification of .CI, and am going to try TDCISD.
By the way, spinflip TDDFT geometry search with GAMESS using b3lyp and 631G(d,p) of phenyl methyl ketone can give T1S0 vertical gap in ethanol, where the ground state targeted is not spincontaminated at all, within 0.06 eVs deviating from the 00 band summarized by Prof. Wagner et al. previously in Communications to the Editor in JACS, and there are two imaginary frequencies during the following Hessian analysis below 380 cm^1. The imaginary frequencies perhaps should be eliminated. Perhaps I should follow the test using spinflip TDDFT instead of DFT to do Hessian, so the Hessian analysis here is not proper for ZPVE.
I will use Dalton2018 to calculate the phosphorescence lifetime.
I don't know whether Monte carlo minimum search, illustrating Monte Carlo AI, can be applied first to facilitate an triplet geometry optimization.
Very Best regards!
Last edited by xiongyan21 on 05 Sep 2020, 15:21, edited 10 times in total.

 Posts: 290
 Joined: 24 Sep 2014, 08:36
 First name(s): yan
 Last name(s): xiong
 Affiliation: CENTRAL CHINA NORMAL UNIVERSITY
 Country: China
Re: Can Dalton calculate oscillator strength
Using the triplet geometry found by SpinFlip TDDFT in ethanol with GAMESS and with the following input for Dalton 2018 and 631G
**DALTON INPUT
.DIRECT
.PARALLEL
.RUN RESPONSE
**INTEGRALS
.SPINORBIT
**WAVE FUNCTIONS
.DFT
B3LYP
*SCF INPUT
**RESPONS
*QUADRATIC
.ECPHOS
.PHOSPHORESCENCE
.PRINT
10
.ROOTS
3
**END OF INPUT,
the spinorbit gap got in vacuum is within 0.04 eVs deviating from that got in ethanol with GAMESS SPINFLIP TDDFT, and the phosphorescence lifetime is around 0.3 ms(in terms of effective charge spinorbit integrals) deviating from the experimental 8 ms in ethanol, i.e., the shortlived emission. It is estimated that in solvents, the longlived phosphorescence is 7576 kcal and the lifetime is from 0.35 to 2.25 seconds in solvents or gel; Dalton2018 gives one agreeing with the above energy and the other of about 1 second to 1.3 seconds in vacuum.
Prof. Wagner suggested the shortlived one was an mixture of the two triplet states, and the longerlived one, which was calculated (not directly from experiments) by Prof. Lamola，who published the experimental data in another reference, in ethanol, might well involve the nonequilibrium loss an alpha proton for other ketones, but Prof. Lamola gave another explanantion first.
We can see spinflip TDDFT method from the group headed by Prof.A.I.Krylov, once an M.Sc in Moscow State University, a PhD in The Hebrew University of Jerusalem, a post PhD in Prof. M. HeadGordon's group, and now a US university distinguished professor having winned Dirac medal,Theoretical Chemistry Award, and Bessel Research Award, can facilitate the phosphorecence emission wavelength estimation through the vertical T1S0 gaps directly provided after one run. They have developed and are developing other methods, eg. spinflip coupled cluter methods(not in GAMESS), by spinflip ansatz used for the singlettriplet gaps of diradicals, triradicals, bond breaking situations, etc. as well as for the search of conical intersections, and are applying the tools in combustion, bioimaging, solar energy, etc., collaborating with experimental laboratories.
Thanks a lot for the suggestion by Prof. Ruud of using the triplet structure found by other software to calculate its phosphorescence lifetime.
Very Best Regards!
**DALTON INPUT
.DIRECT
.PARALLEL
.RUN RESPONSE
**INTEGRALS
.SPINORBIT
**WAVE FUNCTIONS
.DFT
B3LYP
*SCF INPUT
**RESPONS
*QUADRATIC
.ECPHOS
.PHOSPHORESCENCE
10
.ROOTS
3
**END OF INPUT,
the spinorbit gap got in vacuum is within 0.04 eVs deviating from that got in ethanol with GAMESS SPINFLIP TDDFT, and the phosphorescence lifetime is around 0.3 ms(in terms of effective charge spinorbit integrals) deviating from the experimental 8 ms in ethanol, i.e., the shortlived emission. It is estimated that in solvents, the longlived phosphorescence is 7576 kcal and the lifetime is from 0.35 to 2.25 seconds in solvents or gel; Dalton2018 gives one agreeing with the above energy and the other of about 1 second to 1.3 seconds in vacuum.
Prof. Wagner suggested the shortlived one was an mixture of the two triplet states, and the longerlived one, which was calculated (not directly from experiments) by Prof. Lamola，who published the experimental data in another reference, in ethanol, might well involve the nonequilibrium loss an alpha proton for other ketones, but Prof. Lamola gave another explanantion first.
We can see spinflip TDDFT method from the group headed by Prof.A.I.Krylov, once an M.Sc in Moscow State University, a PhD in The Hebrew University of Jerusalem, a post PhD in Prof. M. HeadGordon's group, and now a US university distinguished professor having winned Dirac medal,Theoretical Chemistry Award, and Bessel Research Award, can facilitate the phosphorecence emission wavelength estimation through the vertical T1S0 gaps directly provided after one run. They have developed and are developing other methods, eg. spinflip coupled cluter methods(not in GAMESS), by spinflip ansatz used for the singlettriplet gaps of diradicals, triradicals, bond breaking situations, etc. as well as for the search of conical intersections, and are applying the tools in combustion, bioimaging, solar energy, etc., collaborating with experimental laboratories.
Thanks a lot for the suggestion by Prof. Ruud of using the triplet structure found by other software to calculate its phosphorescence lifetime.
Very Best Regards!
Last edited by xiongyan21 on 10 Sep 2020, 12:26, edited 8 times in total.

 Posts: 290
 Joined: 24 Sep 2014, 08:36
 First name(s): yan
 Last name(s): xiong
 Affiliation: CENTRAL CHINA NORMAL UNIVERSITY
 Country: China
Re: Can Dalton calculate oscillator strength
Spinflip TDDFT with GAMESS using 631G(d,p) of gaseous pyrazine give T1S0 gap around 0.12 eVs deiviating from Prof. Nguyen's TDDFT theoretical work with 631+G(d,p) published in JCP before, and about 0.18 eVs deviating from the experimental data indexed there. This triplet geometry is very easy to be found.
The triplet structure perhaps has an imaginary frequncy during latter UHF DFT analysis which perhaps could be eliminated.
I will use Dalton2018 to calculate the phosphorescence.
Very Best Regards!
The triplet structure perhaps has an imaginary frequncy during latter UHF DFT analysis which perhaps could be eliminated.
I will use Dalton2018 to calculate the phosphorescence.
Very Best Regards!
Last edited by xiongyan21 on 23 Sep 2020, 00:43, edited 18 times in total.

 Posts: 290
 Joined: 24 Sep 2014, 08:36
 First name(s): yan
 Last name(s): xiong
 Affiliation: CENTRAL CHINA NORMAL UNIVERSITY
 Country: China
Re: Can Dalton calculate oscillator strength
According to Prof. Nguyen: "A mean absolute error ~MAE! of 0.08 eV is obtained by
B3LYP/6311G(d,p) for 37 experimental energies of organic chromophores. Interestingly, the MAE increases to
0.14 eV after including zeropoint energy corrections."
The increase of convergence tightness and the use of 631+G(d,p) makes the later seminumerical Hessian analysis of pyrazine triplet geometry found by spinflip TDDFT have no imaginary frequency, in the existence of spinflip TDDFT, and the T1S0 vertical gap within 0.21 eVs deviating from the experimental value; the adiabatic theoretical/ experimental discrypency of T1S0 gap with ZPVE within 0.06 eVs. ZPVE is very samll for this molecule, and the ground state geometry found with C1 symmetry is nonplanar.
Monte Carlo UHF minimum searches have been successful for the triplet phenyl methyl ketone in vacuum and in ethanol, further spinflip triplet geometry optimization and its Hessian anaysis will be done.
I will use UHF DFT to find the triplet structure of pyrazine beginning with the one got by spinflip TDDFT.
Very Best Regards！
B3LYP/6311G(d,p) for 37 experimental energies of organic chromophores. Interestingly, the MAE increases to
0.14 eV after including zeropoint energy corrections."
The increase of convergence tightness and the use of 631+G(d,p) makes the later seminumerical Hessian analysis of pyrazine triplet geometry found by spinflip TDDFT have no imaginary frequency, in the existence of spinflip TDDFT, and the T1S0 vertical gap within 0.21 eVs deviating from the experimental value; the adiabatic theoretical/ experimental discrypency of T1S0 gap with ZPVE within 0.06 eVs. ZPVE is very samll for this molecule, and the ground state geometry found with C1 symmetry is nonplanar.
Monte Carlo UHF minimum searches have been successful for the triplet phenyl methyl ketone in vacuum and in ethanol, further spinflip triplet geometry optimization and its Hessian anaysis will be done.
I will use UHF DFT to find the triplet structure of pyrazine beginning with the one got by spinflip TDDFT.
Very Best Regards！
Last edited by xiongyan21 on 29 Aug 2020, 10:00, edited 19 times in total.

 Posts: 290
 Joined: 24 Sep 2014, 08:36
 First name(s): yan
 Last name(s): xiong
 Affiliation: CENTRAL CHINA NORMAL UNIVERSITY
 Country: China
Re: Can Dalton calculate oscillator strength
The phosphorescence emission calculated with Dalton2018 wit 631++G** for pyrazine in vacuum using geometry found by spinflip TDDFT with GAMESS is within 0.1 eVs deviating from the experimental value and the phosphorescence liftime less than 0.060.08 seconds. Neverthless, in the terminal, there is such information
1976: *** INFO, negative eigenvalues in reduced matrix. Symmetry 1, triplet F Dimension of reduced matrix = 24
1991: *** INFO, negative eigenvalues in reduced matrix. Symmetry 1, triplet F Dimension of reduced matrix = 30
1998: *** INFO, negative eigenvalues in reduced matrix. Symmetry 1, triplet F Dimension of reduced matrix = 30
2061: *** INFO, negative eigenvalues in reduced matrix. Symmetry 1, triplet F Dimension of reduced matrix = 6
2068: *** INFO, negative eigenvalues in reduced matrix. Symmetry 1, triplet F Dimension of reduced matrix = 6
2083: *** INFO, negative eigenvalues in reduced matrix. Symmetry 1, triplet F Dimension of reduced matrix = 12
2090: *** INFO, negative eigenvalues in reduced matrix. Symmetry 1, triplet F Dimension of reduced matrix = 12
2105: *** INFO, negative eigenvalues in reduced matrix. Symmetry 1, triplet F Dimension of reduced matrix = 18
2112: *** INFO, negative eigenvalues in reduced matrix. Symmetry 1, triplet F Dimension of reduced matrix = 18
Does it matter? Perhaps only spinflip methods can allow negative roots, and of course here the negative eigenvalues have been ignored.
631+G* gives the discrepancy within 0.08 eVs.
Very Best Regards！
1976: *** INFO, negative eigenvalues in reduced matrix. Symmetry 1, triplet F Dimension of reduced matrix = 24
1991: *** INFO, negative eigenvalues in reduced matrix. Symmetry 1, triplet F Dimension of reduced matrix = 30
1998: *** INFO, negative eigenvalues in reduced matrix. Symmetry 1, triplet F Dimension of reduced matrix = 30
2061: *** INFO, negative eigenvalues in reduced matrix. Symmetry 1, triplet F Dimension of reduced matrix = 6
2068: *** INFO, negative eigenvalues in reduced matrix. Symmetry 1, triplet F Dimension of reduced matrix = 6
2083: *** INFO, negative eigenvalues in reduced matrix. Symmetry 1, triplet F Dimension of reduced matrix = 12
2090: *** INFO, negative eigenvalues in reduced matrix. Symmetry 1, triplet F Dimension of reduced matrix = 12
2105: *** INFO, negative eigenvalues in reduced matrix. Symmetry 1, triplet F Dimension of reduced matrix = 18
2112: *** INFO, negative eigenvalues in reduced matrix. Symmetry 1, triplet F Dimension of reduced matrix = 18
Does it matter? Perhaps only spinflip methods can allow negative roots, and of course here the negative eigenvalues have been ignored.
631+G* gives the discrepancy within 0.08 eVs.
Very Best Regards！
Last edited by xiongyan21 on 27 Aug 2020, 15:57, edited 4 times in total.

 Posts: 290
 Joined: 24 Sep 2014, 08:36
 First name(s): yan
 Last name(s): xiong
 Affiliation: CENTRAL CHINA NORMAL UNIVERSITY
 Country: China
Re: Can Dalton calculate oscillator strength
With GAMESS the UHF optimization of pyrazine in vacuum with the increased convergence tightness, 631+G(d,p) and pbe0 gives an imaginary frequency about 172 cm^1, and when the basis set is changed into 6311++G(2d,2p), that becomes about 260 cm^1. The adiabatic T1S0 gap after ZPVE(~0.21 eVs) is within 0.26 eVs smaller than the experimentl value.
Spinflip has been expanded by Prof. M. S. Gordon, once a PhD student of Nobel Prize Winner John A Pople, now a distinguished US university professor and Ames Lab senior researcher having winned many metals, prizes, and prestigies for achievements, who are leading his group continue to increase famous GAMESS's quantum chemical capabilities rapidly, and his student to spinflip ORMASCI, which still lacks dynamic correlation but gives results almost identical to those obtained using other multireference methods (for molecules in his recent article) and has just been implemented in the newest GAMESS. FCI is the strength of Prof. Gordon's Group and Ames Lab. GAMESS now can also do spinflip MRMP2 using another kind of FCI for energy.
The phosphorescence of this structure with Dalton2018 gives a descrepancy of about 0.23 eVs using 631+G* when compare with experiments, and the lifetime around 0.06 to 0.08 seconds, accompanying the above warnings in the terminal.
If an imaginary frequency exists, the vibrational analysis will exclude it, making ZPVE not that reliable.
Very Best Regards!
Spinflip has been expanded by Prof. M. S. Gordon, once a PhD student of Nobel Prize Winner John A Pople, now a distinguished US university professor and Ames Lab senior researcher having winned many metals, prizes, and prestigies for achievements, who are leading his group continue to increase famous GAMESS's quantum chemical capabilities rapidly, and his student to spinflip ORMASCI, which still lacks dynamic correlation but gives results almost identical to those obtained using other multireference methods (for molecules in his recent article) and has just been implemented in the newest GAMESS. FCI is the strength of Prof. Gordon's Group and Ames Lab. GAMESS now can also do spinflip MRMP2 using another kind of FCI for energy.
The phosphorescence of this structure with Dalton2018 gives a descrepancy of about 0.23 eVs using 631+G* when compare with experiments, and the lifetime around 0.06 to 0.08 seconds, accompanying the above warnings in the terminal.
If an imaginary frequency exists, the vibrational analysis will exclude it, making ZPVE not that reliable.
Very Best Regards!
Last edited by xiongyan21 on 23 Sep 2020, 15:29, edited 27 times in total.

 Posts: 290
 Joined: 24 Sep 2014, 08:36
 First name(s): yan
 Last name(s): xiong
 Affiliation: CENTRAL CHINA NORMAL UNIVERSITY
 Country: China
Re: Can Dalton calculate oscillator strength
All the Dalton2018 calculations on Ubuntu 18.04 may need to be carefully examined, because although they are finished and the results can be compared with experimental data, it may contain serious errors.
Very Best Regards!
Very Best Regards!
Last edited by xiongyan21 on 29 Aug 2020, 09:31, edited 1 time in total.

 Posts: 395
 Joined: 27 Jun 2013, 18:44
 First name(s): Hans Jørgen
 Middle name(s): Aagaard
 Last name(s): Jensen
 Affiliation: Universith of Southern Denmark
 Country: Denmark
Re: Can Dalton calculate oscillator strength
Negative eigenvalues in reduced matrices for triplet properties (as spinorbit operators) are a sign that your reference state is triplet unstable, i.e. that an unrestricted DFT would give a lower energy (unrestricted DFT is not implemented in Dalton).
I think you should be able to get reasonable results by invoking the .TDA approximation, which is often better than the full model for triplet properties because of the triplet instability issue.
I think you should be able to get reasonable results by invoking the .TDA approximation, which is often better than the full model for triplet properties because of the triplet instability issue.

 Posts: 290
 Joined: 24 Sep 2014, 08:36
 First name(s): yan
 Last name(s): xiong
 Affiliation: CENTRAL CHINA NORMAL UNIVERSITY
 Country: China
Re: Can Dalton calculate oscillator strength
Dear Prof. Jessen
Thanks a lot for your elaboration, which makes me know one triplet structure obtained from one operator may not be stable when using others.
For the first excited state of thiabendazole searched in water, both Hessian analysis methods exhibit no imaginary frequency.
As for pyrazine, Prof. Nguyen said all the triplet states of all the molecules studied in his article except for planar stilbene , have no imaginary frequency, but spinflip TDDFT nonplanar structure found by GAMESS does have when the frequency analysis is done with UHF DFT.
Prof. Gordon said in his recent article CIS, RPA and TDDFT notoriously fail to describe the correct topology of conical intersections. Perhaps because of electron correlation, it is very difficult to find the true minimum triplet geometry using UHF DFT and using spinflip TDDFT in terms of later UHF Hessian analysis. I will carefully examine the bond lengths and angles of the triplet geometry found by spinflip TDDFT , compare them with those in several previously published articles, and try the new methods. I have already repeated the test of search the quintet state of TMM in the article using ROHF with spinflip ORMAS.
He also comments on if a PES （1A2 state of a molecule) is very flat, optimization including dynamic correlation to more quantitatively ascertain its shape is required, which I think is also beneficial to THz and low frequency RAMAN analyses.
Very Best Regard!
Thanks a lot for your elaboration, which makes me know one triplet structure obtained from one operator may not be stable when using others.
For the first excited state of thiabendazole searched in water, both Hessian analysis methods exhibit no imaginary frequency.
As for pyrazine, Prof. Nguyen said all the triplet states of all the molecules studied in his article except for planar stilbene , have no imaginary frequency, but spinflip TDDFT nonplanar structure found by GAMESS does have when the frequency analysis is done with UHF DFT.
Prof. Gordon said in his recent article CIS, RPA and TDDFT notoriously fail to describe the correct topology of conical intersections. Perhaps because of electron correlation, it is very difficult to find the true minimum triplet geometry using UHF DFT and using spinflip TDDFT in terms of later UHF Hessian analysis. I will carefully examine the bond lengths and angles of the triplet geometry found by spinflip TDDFT , compare them with those in several previously published articles, and try the new methods. I have already repeated the test of search the quintet state of TMM in the article using ROHF with spinflip ORMAS.
He also comments on if a PES （1A2 state of a molecule) is very flat, optimization including dynamic correlation to more quantitatively ascertain its shape is required, which I think is also beneficial to THz and low frequency RAMAN analyses.
Very Best Regard!
Last edited by xiongyan21 on 19 Sep 2020, 06:36, edited 1 time in total.

 Posts: 290
 Joined: 24 Sep 2014, 08:36
 First name(s): yan
 Last name(s): xiong
 Affiliation: CENTRAL CHINA NORMAL UNIVERSITY
 Country: China
Re: Can Dalton calculate oscillator strength
For the spinflip found triplet geometry with 631+G(d,p) the lengths of NC, CC, CH are 0.056, 0.045, 0.01 Bohrs deviating from the empirically corrected values of 3B3u reported in J. Phys. Chem, and the bond angles are about 2, 1, 1 degrees different, the discrepencies from those of MCSCF STO3G with GAMESS also reported in that reference are 0.117, 0.076, 0.019 Bohrs; 0.63, 0.32, 1.41 degrees.
UHF DFT triplet structure gives all the three bond angles are around 119120 degrees, with the angle deviances of 1.24, 0.61, 0.56 degrees, and the bond length discrepencies are 0.037,0.045, 0.005 Bohrs when compared with the empirically corrected data.
No experimental data of these has been found.
It seems basis set does not make a significant difference, and the imaginary frequency should be eliminated to get a proper match, although the reference didn't do Hessian analysis.
According to an article in Theo. Chem. Acc. in 2011 for spinflip TDDFT in another package, for the 1A3 of COH2, the angle of HCH is 9 degrees smaller then the experimental value when using b3lyp, and the deviances of bond length of CH, CO are 0.038, 0.038, respectively.
Very Best Regards!
UHF DFT triplet structure gives all the three bond angles are around 119120 degrees, with the angle deviances of 1.24, 0.61, 0.56 degrees, and the bond length discrepencies are 0.037,0.045, 0.005 Bohrs when compared with the empirically corrected data.
No experimental data of these has been found.
It seems basis set does not make a significant difference, and the imaginary frequency should be eliminated to get a proper match, although the reference didn't do Hessian analysis.
According to an article in Theo. Chem. Acc. in 2011 for spinflip TDDFT in another package, for the 1A3 of COH2, the angle of HCH is 9 degrees smaller then the experimental value when using b3lyp, and the deviances of bond length of CH, CO are 0.038, 0.038, respectively.
Very Best Regards!

 Posts: 290
 Joined: 24 Sep 2014, 08:36
 First name(s): yan
 Last name(s): xiong
 Affiliation: CENTRAL CHINA NORMAL UNIVERSITY
 Country: China
Re: Can Dalton calculate oscillator strength
As for pyridine, spinflip TDDFT, although using C1 symmetry a plananr geometry is found, gives vertical T1S0 gap discrepencies of about 0.12, 0.05 and 0.06 eVs using PBE0, B3LYP and B3LYPV1R(matching electron gas formula in some packages), respectively, when compared with the experimental data of 3B2 indexed in the article of Prof. Bousquet et al in France and Japan. Spinflip TDDFT perhaps improves the gap even using 631G(d.p). But their results are very close to the experimental data indexed by Prof. Nguyen of Air Force Research Laboratory in US, where the geometry is nonplanar and the experimental result matches adiabatic T1S0 gap calculated by TDDFT. The spinflip TDDFT vertical T1S0 gap is 0.3 eVs different from that experimental data. Of course, there is an imaginary frequency of the planar triplet geometry found, making ZPVE unreliable.
For the S1S0 gap of pyrazine, B3lyp with C1 symmetry can give adiabatic one deviating 0.2 eVs, respectively, from the experimental values of 1B1 and 2B1 indexed in the article of Prof. Klerier, in Los Alamos National Laboratory, et al. in US, agreeing very well with those from HF calculations using symmetry there. ZPVE makes the gap unacceptable. The descrepencies of angle CNC and lengths of CN, CC are about 6.9 degrees, 0.038 and 0.048 angstroms when compared with experimental data. Since the reference emphasizes localized excitation, delocalized orbitals are no longer used, and one localization method proposed by distinguished Professor Emerituz Ludenberg in Iowa State University and Ames Lab, etc., gives very similar results but changes the frequencies in THz region. I am not sure whether the localization is suitable here.
Very Best Regards!
For the S1S0 gap of pyrazine, B3lyp with C1 symmetry can give adiabatic one deviating 0.2 eVs, respectively, from the experimental values of 1B1 and 2B1 indexed in the article of Prof. Klerier, in Los Alamos National Laboratory, et al. in US, agreeing very well with those from HF calculations using symmetry there. ZPVE makes the gap unacceptable. The descrepencies of angle CNC and lengths of CN, CC are about 6.9 degrees, 0.038 and 0.048 angstroms when compared with experimental data. Since the reference emphasizes localized excitation, delocalized orbitals are no longer used, and one localization method proposed by distinguished Professor Emerituz Ludenberg in Iowa State University and Ames Lab, etc., gives very similar results but changes the frequencies in THz region. I am not sure whether the localization is suitable here.
Very Best Regards!
Last edited by xiongyan21 on 13 Sep 2020, 02:36, edited 3 times in total.

 Posts: 395
 Joined: 27 Jun 2013, 18:44
 First name(s): Hans Jørgen
 Middle name(s): Aagaard
 Last name(s): Jensen
 Affiliation: Universith of Southern Denmark
 Country: Denmark
Re: Can Dalton calculate oscillator strength
Yes. If you e.g. search for "dft spinorbit dalton" on google, you get some informative links.

 Posts: 290
 Joined: 24 Sep 2014, 08:36
 First name(s): yan
 Last name(s): xiong
 Affiliation: CENTRAL CHINA NORMAL UNIVERSITY
 Country: China
Re: Can Dalton calculate oscillator strength
Dear Prof. Jensen
I have understood the article, Influence of triplet instabilities in TDDFT of Prof.Tozer, an RSC metal winner, et al. in Durham University in United Kindom.
By the way, Prof. Tozer is one of the authors of the diagnostic of excited states natures.
GAMESS UHF or ROHF TDDFT can only handle with the identical multiplicities of the ground and excited states。
Can Dalton2018 do this?
Spinflip TDDFT pbe0 T1S0 vertical gap of furan, almost the same as that in the article of Prof. Bousquet, et al., using TDDFT and 6311++G(2d,2p), is about 0.40 eVs deviating from the experimental value in the article of Prof.Flicker, et al. , in California Institue of Technology in US. Adiabatic pbe0 T1(spinflip TDDFT) S0(UHF DFT) gap is very small, thus failing. Spinflip TDDFT, TDDFT and DFT with pbe0 fail to deal with this. Spinflip CIS gives a larger descrepency of 0.53 eVs when compared with experimental data, and there is no imaginary frequency of the optimized structure during the subsequent Hessian analysis in the existence of spinflip CIS, and the spinsquared operators are 0.2844 and 2.200 where perhaps a little spin contaminations makes the gap underestimted.
Very Best Regards!
I have understood the article, Influence of triplet instabilities in TDDFT of Prof.Tozer, an RSC metal winner, et al. in Durham University in United Kindom.
By the way, Prof. Tozer is one of the authors of the diagnostic of excited states natures.
GAMESS UHF or ROHF TDDFT can only handle with the identical multiplicities of the ground and excited states。
Can Dalton2018 do this?
Spinflip TDDFT pbe0 T1S0 vertical gap of furan, almost the same as that in the article of Prof. Bousquet, et al., using TDDFT and 6311++G(2d,2p), is about 0.40 eVs deviating from the experimental value in the article of Prof.Flicker, et al. , in California Institue of Technology in US. Adiabatic pbe0 T1(spinflip TDDFT) S0(UHF DFT) gap is very small, thus failing. Spinflip TDDFT, TDDFT and DFT with pbe0 fail to deal with this. Spinflip CIS gives a larger descrepency of 0.53 eVs when compared with experimental data, and there is no imaginary frequency of the optimized structure during the subsequent Hessian analysis in the existence of spinflip CIS, and the spinsquared operators are 0.2844 and 2.200 where perhaps a little spin contaminations makes the gap underestimted.
Very Best Regards!
Last edited by xiongyan21 on 22 Sep 2020, 03:49, edited 14 times in total.

 Posts: 395
 Joined: 27 Jun 2013, 18:44
 First name(s): Hans Jørgen
 Middle name(s): Aagaard
 Last name(s): Jensen
 Affiliation: Universith of Southern Denmark
 Country: Denmark
Re: Can Dalton calculate oscillator strength
The answer is yes, using quadratic response. II can see on this forum that you used this yourself in 2016 ...

 Posts: 290
 Joined: 24 Sep 2014, 08:36
 First name(s): yan
 Last name(s): xiong
 Affiliation: CENTRAL CHINA NORMAL UNIVERSITY
 Country: China
Re: Can Dalton calculate oscillator strength
I am trying Spinflip ORMASCI, and arbitrarily send the first partition of active space 29 electrons, 36 active electrons in total among which there are 19 alpha and 17 beta electrons based on spinflip TDDFT.
One of the original ORMARS programmer is Dr. M. Schmidt in Ames Lab and Prof. Gordon's group, already retired, whom has been deemed as a longtime tutor and friend by many PhD graduates.
I also will look into the forum and try that using Dalton2018.
Very Best Regards!
One of the original ORMARS programmer is Dr. M. Schmidt in Ames Lab and Prof. Gordon's group, already retired, whom has been deemed as a longtime tutor and friend by many PhD graduates.
I also will look into the forum and try that using Dalton2018.
Very Best Regards!
Last edited by xiongyan21 on 03 Oct 2020, 02:17, edited 1 time in total.

 Posts: 290
 Joined: 24 Sep 2014, 08:36
 First name(s): yan
 Last name(s): xiong
 Affiliation: CENTRAL CHINA NORMAL UNIVERSITY
 Country: China
Re: Can Dalton calculate oscillator strength
Spinflip OMARS geometry search of the triplet structure of furan is OK, where a planar structure is found. The active orbitals are 115 and the active electrons are 36 in 3 partitions with nstate being equal to 3, where the final energy is significant different from the ROHF one. FCI should be applied to the ground state to obtain the discrepancy, but the spinflip ORMAS search uses analytical gradient!
Actually, Prof. Flicker said in his article that the lowest inealstic tranition occurs in the region 3.3 eV to 4.9 eV, with a maxmum intensity at 3.99 eV. The other singlettriplet gap is at 5.22 eV where the electronimpact energyloss spectra peak is higher.
The adibatic T1(sformas)S0(ormas, nstate=4, mult=1, the other same as spinflip) also is vanishly small. I will examine what is wrong here.
The vertical gaps T1(sformasCI)S0(ormasCI, the first of nstate=420, mult=1) is around 0.93 eVs deviating from 3.99 eVs. The inputs should be examined to verify whether the results can match the experiments.
Prof. Bousquet, et al using SACCI give the gap of 4.02 eV. and they believe the experimental value is close to 4.1 eV based on the experiments of Prof. Flicker, et al.
Very Best Regards!
Actually, Prof. Flicker said in his article that the lowest inealstic tranition occurs in the region 3.3 eV to 4.9 eV, with a maxmum intensity at 3.99 eV. The other singlettriplet gap is at 5.22 eV where the electronimpact energyloss spectra peak is higher.
The adibatic T1(sformas)S0(ormas, nstate=4, mult=1, the other same as spinflip) also is vanishly small. I will examine what is wrong here.
The vertical gaps T1(sformasCI)S0(ormasCI, the first of nstate=420, mult=1) is around 0.93 eVs deviating from 3.99 eVs. The inputs should be examined to verify whether the results can match the experiments.
Prof. Bousquet, et al using SACCI give the gap of 4.02 eV. and they believe the experimental value is close to 4.1 eV based on the experiments of Prof. Flicker, et al.
Very Best Regards!
Last edited by xiongyan21 on 23 Sep 2020, 05:33, edited 7 times in total.

 Posts: 290
 Joined: 24 Sep 2014, 08:36
 First name(s): yan
 Last name(s): xiong
 Affiliation: CENTRAL CHINA NORMAL UNIVERSITY
 Country: China
Re: Can Dalton calculate oscillator strength
Using the geometry found by spinflip ORMAS, the vertical T1(MRPT2) S0(MCQDPT, by spinflip MRMP2) can be 0.09 eV deviating from 3.99 eV, so is the gap of T1(MCQDPT) S0(MCQDPT, by spinflip MRMP2) when the active orbitals and active electrons are all two. The triplet state energy calculation includes another kind of FCI created by Ames Lab, a full second order CI with the exact Hessian and a mcscf orbital optimization in advance, where the active orbitals and active electrons are all two. Perhap spinflip MRMP2 only allows some of the methods. MCQDPT is only for CASSCF wavefunctions.
The spinflip ORMAS optimization uses delocalized internal coordinates to accelerate its convergence. The basis set here is 631+G(d,p).
Using the geometry found by spinflip TDDFT PBE0, also a planar one, the vertical T1(MRPT2) S0(MCQDPT, by spinflip MRMP2) can be about 0.29 eV deviating from 3.99 eV.
For furan, the quantum chemical calculation of the T1S0 gap is challenging, but in the reference, professors perhaps use the triplet geometry from experiments or through their advanced chemical attainments. It seems perhaps we can say the vertical gap calculated by MRMP2 and spinflip MRMP2 based on the triplet geometry found by spinflip ormas can agree with the experimental data, i.e., 3.99 eV.
All the above use GAMESS, and natural MOs are used. I will use Dalton2018 to calculate its phosphorescence.
Perhaps this also applies to hemexazole, there being no reported phosphorescence and fluorescence.
Very Best Regards!
The spinflip ORMAS optimization uses delocalized internal coordinates to accelerate its convergence. The basis set here is 631+G(d,p).
Using the geometry found by spinflip TDDFT PBE0, also a planar one, the vertical T1(MRPT2) S0(MCQDPT, by spinflip MRMP2) can be about 0.29 eV deviating from 3.99 eV.
For furan, the quantum chemical calculation of the T1S0 gap is challenging, but in the reference, professors perhaps use the triplet geometry from experiments or through their advanced chemical attainments. It seems perhaps we can say the vertical gap calculated by MRMP2 and spinflip MRMP2 based on the triplet geometry found by spinflip ormas can agree with the experimental data, i.e., 3.99 eV.
All the above use GAMESS, and natural MOs are used. I will use Dalton2018 to calculate its phosphorescence.
Perhaps this also applies to hemexazole, there being no reported phosphorescence and fluorescence.
Very Best Regards!
Last edited by xiongyan21 on 10 Oct 2020, 04:04, edited 4 times in total.
Who is online
Users browsing this forum: No registered users and 5 guests