RSPCTL WARNING-MAXIMUM NUMBER OF MICROITERATIONS

Find answers or ask questions regarding Dalton calculations.
Please upload an output file showing the problem, if applicable.
(It is not necessary to upload input files, they can be found in the output file.)

Post Reply
ndhieu92
Posts: 31
Joined: 04 Sep 2019, 07:51
First name(s): Hieu
Middle name(s): Dinh
Last name(s): Nguyen
Affiliation: Kunsan Univ
Country: Korea, Republic of

RSPCTL WARNING-MAXIMUM NUMBER OF MICROITERATIONS

Post by ndhieu92 » 22 Oct 2019, 09:05

WIsh you have a good day everyone.

I am trying to calculate for the spin-orbit coupling constant, but I have got a problem. Eventhough the calculation finished without any error, but the result of SOC value between 2 times of calculation is different between each other. Firstly, I checked the input file (mol. and dal. file) and I am sure that it is correct. So I checked the output file and then I always see the warning is: RSPCTL WARNING-MAXIMUM NUMBER OF MICROITERATIONS, 60. So now I am trying to add into dal. file: .MAX MICRO ITERATIONS and increase it into 200. Will this fix my problem? I attached 2 output file and also the new dal.file with my thinking. Please let me know your idea.

Thank you so much
Attachments
SPINORBIT_rsplrso_AMFI_quadratic_autosymmetry_MAXMICROITERATIONS.dal
(441 Bytes) Downloaded 17 times
SPINORBIT_rsplrso_AMFI_quadratic_autosymmetry_SOC_BDELA.out
(404.05 KiB) Downloaded 21 times
SPINORBIT_rsplrso_AMFI_quadratic_autosymmetry_2_SOC_BDELA.out
(450.68 KiB) Downloaded 17 times

hjaaj
Posts: 334
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: RSPCTL WARNING-MAXIMUM NUMBER OF MICROITERATIONS

Post by hjaaj » 23 Oct 2019, 13:55

Your spin-orbit coupling constants are garbage, because your calculation is triplet unstable (complex eigenvalues in the calculation of lowest triplet excitation). The differences between the two calculations are caused by numerical noise because of the non-convergence.

Either your DFT iterations do not converge to the ground state (do you get same total energy if you try with another program?), or B3LYP is a bad model for this case.

By the way, your ".MAX IT" keyword has no effect because the ".MAX IT" in that place specify the max. number of geometry iterations in a geometry optimization. But it does not help to increase the correct max iterations for response, because the convergence is stuck because of the above mentioned problems.

ndhieu92
Posts: 31
Joined: 04 Sep 2019, 07:51
First name(s): Hieu
Middle name(s): Dinh
Last name(s): Nguyen
Affiliation: Kunsan Univ
Country: Korea, Republic of

Re: RSPCTL WARNING-MAXIMUM NUMBER OF MICROITERATIONS

Post by ndhieu92 » 23 Oct 2019, 15:05

hjaaj wrote:
23 Oct 2019, 13:55
Your spin-orbit coupling constants are garbage, because your calculation is triplet unstable (complex eigenvalues in the calculation of lowest triplet excitation). The differences between the two calculations are caused by numerical noise because of the non-convergence.

Either your DFT iterations do not converge to the ground state (do you get same total energy if you try with another program?), or B3LYP is a bad model for this case.

By the way, your ".MAX IT" keyword has no effect because the ".MAX IT" in that place specify the max. number of geometry iterations in a geometry optimization. But it does not help to increase the correct max iterations for response, because the convergence is stuck because of the above mentioned problems.
Thank you Dr. Hans for your reply. I am sorry but may I ask you a question: Do you mean the SCF energy which was shown in the output file is total energy? And I am thinking to try with another program like molpro to confirm it. I will try to change to MCSCF to see how it works.

hjaaj
Posts: 334
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: RSPCTL WARNING-MAXIMUM NUMBER OF MICROITERATIONS

Post by hjaaj » 23 Oct 2019, 15:23

I mean this final energy:

Code: Select all

@    Final DFT energy:          -2824.469146640524  
and these final HOMO-LUMO energies:

Code: Select all

    E(LUMO) :    -0.13437205 au (symmetry 1)
  - E(HOMO) :    -0.14494639 au (symmetry 1)
  ------------------------------------------
    gap     :     0.01057434 au
Note that the HOMO and LUMO orbital energies are very close, so no wonder there is triplet instability.

ndhieu92
Posts: 31
Joined: 04 Sep 2019, 07:51
First name(s): Hieu
Middle name(s): Dinh
Last name(s): Nguyen
Affiliation: Kunsan Univ
Country: Korea, Republic of

Re: RSPCTL WARNING-MAXIMUM NUMBER OF MICROITERATIONS

Post by ndhieu92 » 23 Oct 2019, 16:10

hjaaj wrote:
23 Oct 2019, 15:23
I mean this final energy:

Code: Select all

@    Final DFT energy:          -2824.469146640524  
and these final HOMO-LUMO energies:

Code: Select all

    E(LUMO) :    -0.13437205 au (symmetry 1)
  - E(HOMO) :    -0.14494639 au (symmetry 1)
  ------------------------------------------
    gap     :     0.01057434 au
Note that the HOMO and LUMO orbital energies are very close, so no wonder there is triplet instability.
Maybe I will try to recalculate with gaussian 09 software with same basis set B3LYP/6-31G* in dal.file to confirm the DFT energy and HOMO-LUMO energies.
I will keep updating the result.
Thank you so much Dr. Hans

xiongyan21
Posts: 184
Joined: 24 Sep 2014, 08:36
First name(s): yan
Last name(s): xiong
Affiliation: CENTRAL CHINA NORMAL UNIVERSITY
Country: China

Re: RSPCTL WARNING-MAXIMUM NUMBER OF MICROITERATIONS

Post by xiongyan21 » 02 Nov 2019, 10:55

Dear Prof. Jensen
Could you please give me an rough idea of how large the energy gap between LUMO and HOMO of the triplet state of a molecule would be, because in the uhf spin-flip TDDFT geometry search of thiabendazole with PCM, also can be very small. [/b] The advantage of the spin-flip triplet geometry, MECIs searches lies in it does not need an inital geometry guess.
I will let it supass Cs symmetry beginning and retry.
Last edited by xiongyan21 on 15 Nov 2019, 04:47, edited 4 times in total.

xiongyan21
Posts: 184
Joined: 24 Sep 2014, 08:36
First name(s): yan
Last name(s): xiong
Affiliation: CENTRAL CHINA NORMAL UNIVERSITY
Country: China

Re: RSPCTL WARNING-MAXIMUM NUMBER OF MICROITERATIONS

Post by xiongyan21 » 04 Nov 2019, 02:50

After changing a functional and using C1 symmetry from the very beginning, the apparently successful UHF spin-flip TDDFT search of the triplet state of thiabendazole with PCM using GAMESS gives perhaps satisfactory HOMO and LUMO gaps, respectively. Of course, spin-flip TDDFT may perfer to some functionals. I use wxMacMolplt to view the orbital energies.
Last edited by xiongyan21 on 17 Nov 2019, 06:09, edited 18 times in total.

xiongyan21
Posts: 184
Joined: 24 Sep 2014, 08:36
First name(s): yan
Last name(s): xiong
Affiliation: CENTRAL CHINA NORMAL UNIVERSITY
Country: China

Re: RSPCTL WARNING-MAXIMUM NUMBER OF MICROITERATIONS

Post by xiongyan21 » 04 Nov 2019, 07:52

I have done the virbational analysis of the triplet state using GAMESS and found there is no imaginary frequence, but it has a warning meaning it is not a true stationary point although the convergence criteria are met, thus HESSIAN analysis is not very rigorious. Perhaps with NWCHEM, there is no such warning.
Last edited by xiongyan21 on 24 Nov 2019, 13:24, edited 7 times in total.

xiongyan21
Posts: 184
Joined: 24 Sep 2014, 08:36
First name(s): yan
Last name(s): xiong
Affiliation: CENTRAL CHINA NORMAL UNIVERSITY
Country: China

Re: RSPCTL WARNING-MAXIMUM NUMBER OF MICROITERATIONS

Post by xiongyan21 » 05 Nov 2019, 01:06

It seems this can be understood according to what Dr. Schmidt from Ames Lab told me before that HOMO an LUMO are not very suitable concepts.
Last edited by xiongyan21 on 17 Nov 2019, 06:10, edited 2 times in total.

xiongyan21
Posts: 184
Joined: 24 Sep 2014, 08:36
First name(s): yan
Last name(s): xiong
Affiliation: CENTRAL CHINA NORMAL UNIVERSITY
Country: China

Re: RSPCTL WARNING-MAXIMUM NUMBER OF MICROITERATIONS

Post by xiongyan21 » 05 Nov 2019, 03:34

Dear Dr. Jensen
Yes, you are right.
According to the article of Prof. Anna Krylov in University of Southern California, the traditional single reference method for excited states will fail when the subspaces of the occupied and virtual become nearly degenerate, which is also true for ground states. According to an published article, SF-TDDFT can only deal with one electron excitation and doubly excited states, GAMESS can give energie gaps of the singlet to triplet, triplet to triplet, and triplet to mixed.
Last edited by xiongyan21 on 28 Nov 2019, 15:36, edited 11 times in total.

xiongyan21
Posts: 184
Joined: 24 Sep 2014, 08:36
First name(s): yan
Last name(s): xiong
Affiliation: CENTRAL CHINA NORMAL UNIVERSITY
Country: China

Re: RSPCTL WARNING-MAXIMUM NUMBER OF MICROITERATIONS

Post by xiongyan21 » 08 Nov 2019, 02:29

The article Determination of thiabendazole in foods and waters by solid-phase transmitted room-temperature phosphorescence said the enhanced 472 nm is its naive phosphorescence.
Last edited by xiongyan21 on 17 Nov 2019, 06:23, edited 8 times in total.

xiongyan21
Posts: 184
Joined: 24 Sep 2014, 08:36
First name(s): yan
Last name(s): xiong
Affiliation: CENTRAL CHINA NORMAL UNIVERSITY
Country: China

Re: RSPCTL WARNING-MAXIMUM NUMBER OF MICROITERATIONS

Post by xiongyan21 » 08 Nov 2019, 03:53

6-311++G(2d,2p) can make thiabendazole produce a small imaginary frequency with b3lyp geometry around 55 cm-1 using Cs symmetry and HF frequency analysis, perhaps caused by rotation, etc., while 6-31+(d,p) and b3lyp cannot when using GAMESS. Although it may really does not matter for the excitation energy calculation, it is perhaps a problem for the THz calculation for it.
Somebody recommended the decrease of the symmetry, e.g., using C1, to eliminate it. I have just tried thiabendazole, which also can be considered as to be planar even in C1, using gamess with 6-311++G(2d,2p) and b3lyp, and the approximate 55 cm-1 imaginary frequency has been removed during HF vibrational anaysis with the warning.
For the molecule thiol, if mp2 is used for its geometry optimization, 6-311++(3df, 2p ) is required to eliminate it with seminumerical mp2 gradients.
Last edited by xiongyan21 on 07 Dec 2019, 04:55, edited 28 times in total.

xiongyan21
Posts: 184
Joined: 24 Sep 2014, 08:36
First name(s): yan
Last name(s): xiong
Affiliation: CENTRAL CHINA NORMAL UNIVERSITY
Country: China

Re: RSPCTL WARNING-MAXIMUM NUMBER OF MICROITERATIONS

Post by xiongyan21 » 10 Nov 2019, 10:15

Thiabendazole SF-TDDFT optimized geometry give the triplet to singlet emission energy having a large discrepency ( shorter about 76 nm )from the experimental enhanced data in Determination of thiabendazole in foods and waters by solid-phase transmitted room-temperature phosphorescence.
It is said in the paper the sample and enhancer solutions were deposited on a slip of paper (Whatman No. 4), previously humidified with pH 4.0 acetic/acetate buffer solution, after which the paper was dried and placed between 2 quartz plates. The diffuse transmitted phosphorescence intensity was measured directly at excitation of 303 nm and emission of 472 nm.

Monte Carlo search perhaps is needed for carbendazole first for the ground state due to its rotation,facilitating the geometry optimization at least.
Last edited by xiongyan21 on 23 Nov 2019, 06:22, edited 7 times in total.

xiongyan21
Posts: 184
Joined: 24 Sep 2014, 08:36
First name(s): yan
Last name(s): xiong
Affiliation: CENTRAL CHINA NORMAL UNIVERSITY
Country: China

Re: RSPCTL WARNING-MAXIMUM NUMBER OF MICROITERATIONS

Post by xiongyan21 » 12 Nov 2019, 15:28

The article Highly sensitive and selective room-temperature phosphorescence determination of thiabendazole by the supramolecular interaction of thiabendazole/β-cyclodextrin/triton X-100 gives λex/λem = 298/481 nm resulting from a 1 ∶ 1 ∶ 1 β-cyclodextrin (β-CD)/thiabendazole (TBZ)/triton X-100 (TX-100) supramolecular ternary inclusion complex.

Flow injection analysis with on-line nylon powder extraction for room-temperature phosphorescence determination of thiabendazole
gives maximum λem = 505 nm.

hjaaj
Posts: 334
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: RSPCTL WARNING-MAXIMUM NUMBER OF MICROITERATIONS

Post by hjaaj » 12 Nov 2019, 16:15

Dear Yan Xiong

Sorry, but your comments are not really related to the original posting in this thread by Hieu Nguyen. Furthermore, you discuss items (as UHF spin-flip and SF-TDDFT) which are not implemented in dalton, and therefore the dalton forum is not the right place for such a discussion.

xiongyan21
Posts: 184
Joined: 24 Sep 2014, 08:36
First name(s): yan
Last name(s): xiong
Affiliation: CENTRAL CHINA NORMAL UNIVERSITY
Country: China

Re: RSPCTL WARNING-MAXIMUM NUMBER OF MICROITERATIONS

Post by xiongyan21 » 12 Nov 2019, 16:36

Dear Dr. Jensen
Sorry for this.

Thanks, anyway.

Very Best Regards!

xiongyan21
Posts: 184
Joined: 24 Sep 2014, 08:36
First name(s): yan
Last name(s): xiong
Affiliation: CENTRAL CHINA NORMAL UNIVERSITY
Country: China

Re: RSPCTL WARNING-MAXIMUM NUMBER OF MICROITERATIONS

Post by xiongyan21 » 16 Nov 2019, 15:35

Since the triplet found in water is far from these published experimental data using BHHLYP, perhaps it is not useful for the phosphorescence of this molecule.

Please note GAMESS has incorrect analytic gradient for meta-GGA with SF-TDDFT for now.

GAMESS gives a deviance of 0.367 eV from 478nm using b3lyp and an sperical Dunning basis set in ethanol, and with sperical 6-31G+(d,p), it becomes 0.33 eV, questionable with the default grid.
Photophysical Properties of Benzimidazole and Thiabendazole and Their Homologues. Effect of Substituents and Solvent on the Nature of the Transition gives 453,478, and 500nm in ethanol at 77K with the lifetime of 157 ms.
The article Calculation of the Phosphorescence of Porphyrins by the Density Functional Method got 1.11 eV of the gap of T1 and S0 using 3-21G and B3LYP by DALTON。 The authors think 3-21G underestimated the energy of the T1 state in comparison with the experimental value( delta E=1.58 eV). Actually, what is calculated in the reference is porphin.
I have used UHF SF-TDDFT to calculate the T1 to So gap of porphyrin at nonplanar initial geometry with 3-21G and BHHLYP in toluene with GAMESS, got the transition energy of about 1.36 eV.
Theoretical DFT study of phosphorescence from porphyrins gives the gap to 3B2u of 1.46 eV using 6-31G(d,p) optimized ground state geometry with b3lyp.
Porphyrin Free Base Phosphorescence gives the heavy atom enhanced phosphorescence of 795 nm by time resolved spectra in 50% EPA and 50% ethyl iodide at 77K for proto-porphyrin.
Last edited by xiongyan21 on 09 Dec 2019, 16:03, edited 33 times in total.

xiongyan21
Posts: 184
Joined: 24 Sep 2014, 08:36
First name(s): yan
Last name(s): xiong
Affiliation: CENTRAL CHINA NORMAL UNIVERSITY
Country: China

Re: RSPCTL WARNING-MAXIMUM NUMBER OF MICROITERATIONS

Post by xiongyan21 » 18 Nov 2019, 11:31

It is noted in the recently released GAMESS source code, including 2019.R1, the test, i.e., UHF SF-TDDFT geometry search of the triplet state of H2CO, not converged for GCC greater than 4.4.
Thus, perhaps, fix is wanted.
Last edited by xiongyan21 on 21 Nov 2019, 17:13, edited 1 time in total.

User avatar
magnus
Posts: 478
Joined: 27 Jun 2013, 16:32
First name(s): Jógvan Magnus
Middle name(s): Haugaard
Last name(s): Olsen
Affiliation: Hylleraas Centre, UiT The Arctic University of Norway
Country: Norway

Re: RSPCTL WARNING-MAXIMUM NUMBER OF MICROITERATIONS

Post by magnus » 18 Nov 2019, 12:23

The Dalton forum is also not the right forum to discuss GAMESS calculations.

xiongyan21
Posts: 184
Joined: 24 Sep 2014, 08:36
First name(s): yan
Last name(s): xiong
Affiliation: CENTRAL CHINA NORMAL UNIVERSITY
Country: China

Re: RSPCTL WARNING-MAXIMUM NUMBER OF MICROITERATIONS

Post by xiongyan21 » 06 Dec 2019, 05:08

UHF-SF-TDDFT using accd and b3lyp gives the deviance of 0.346eV from 478 nm in ethanol with the increased grid, still questionable and requiring further increases, above the default, for thiabendazole
Since it is particularly difficult for the quantum chemical caculation of the signlet-triplet gap with the geometry search of the triplet of thiabendazole, i.e., quantitative mismatich between the theoretical and the experimental results, if allowed, FCI, MCSCF, etc., or explicit solvation model might be carried out for the related search.

Post Reply

Who is online

Users browsing this forum: No registered users and 5 guests