ROHF RPA TDDFT Excited Energies Calculation of Sc Failed on Mac using PBE0 and aug-cc-pvdz

From NWChem

Viewed 1181 times, With a total of 11 Posts
Jump to: navigation, search

Forum Regular
Threads 56
Posts 268
The error is
...
Entering Davidson iterations
Unrestricted

 Iter   NTrls   NConv    DeltaV     DeltaE      Time   
---- ------ ------ --------- --------- ---------
1 10 0 0.15E+00 0.10+100 104.5
------------------------------------------------------------------------
tddft_diagon: negative excitation energy 0
------------------------------------------------------------------------
------------------------------------------------------------------------
current input line :
31: task tddft energy
------------------------------------------------------------------------
------------------------------------------------------------------------
This type of error is most commonly associated with calculations not reaching convergence criteria,
------------------------------------------------------------------------
and the input is
echo
memory total 14000 stack 6500 heap 500 global 7000 mb
charge 0


geometry
 Sc         0.00000000     0.00000000     0.00000000

end


basis spherical
Sc library aug-cc-pvdz
end

scf
rohf
end

dft
xc pbe0
mult 2
iterations 500
end

tddft
rPA
nroots 10
end

task tddft energy,

but the ROHF CIS calculation was successful with the following
 ----------------------------------------------------------------------------
Root 1 a 0.008076665 a.u. 0.2198 eV
<S2> = 0.9169
----------------------------------------------------------------------------
...
----------------------------------------------------------------------------
Root 2 a 0.008079542 a.u. 0.2199 eV
<S2> = 0.9168
----------------------------------------------------------------------------
...GAMESS 2016R1 UHF TDDFT gave two 0.048eV degenerate states , whereas DALTON 2016.2 ROHF TDDFT (modified according to Prof. Jensen's proposal on DALTON forum) gave
@ Excited state no: 1 in symmetry 1 ( Ag )
----------------------------------------------

@ Excitation energy : 2.90740354E-03 au
@ 7.91144749E-02 eV; 638.10131 cm-1; 7.6333868 kJ / mol

@ Total energy : -760.42849 au
...
@ Excited state no: 2 in symmetry 1 ( Ag )
----------------------------------------------

@ Excitation energy : 4.31554186E-02 au
@ 1.1743187 eV; 9471.5195 cm-1; 113.30453 kJ / mol

@ Total energy : -760.38825 au

with two excited states (b3g and b2g) degenerate succeeding.
Edited On 8:18:09 PM PST - Mon, Nov 27th 2017 by Xiongyan21

Forum Regular
Threads 56
Posts 268
How can I compare the results and how can I make the above NWCHEM calculation successful?

Very Best Regards!
Edited On 9:07:05 AM PST - Sat, Nov 25th 2017 by Xiongyan21

  • Niri Forum:Admin, Forum:Mod, NWChemDeveloper, bureaucrat, sysop
    Profile
    Send PM
Forum Regular
Threads 4
Posts 207
Hi, Did you use the same exchange-correlation functional in the GAMESS and Dalton calculations as you did for the NWChem one ?

Best,
-Niri Govind

niri.govind@pnnl.gov

Forum Regular
Threads 56
Posts 268
Dear Dr. Govind
 I am sorry not to mention one thing, i.e., since NWCHEM is significantly slower than either
GAMESS(default grids, etc.) or DALTON on this background, the previous calculations by
NWCHEM used PBE0 with aug-cc-pvdz, while the other calculations all used PBE0 with
aug-cc-pvtz.
Now NWCHEM CIS calculation uses PBE0 with aug-cc-pvtz gives:
Root 1 a 0.008804949 a.u. 0.2396 eV
 <S2> =   0.9125

...
 ----------------------------------------------------------------------------
Root 2 a 0.008826526 a.u. 0.2402 eV
<S2> = 0.9127
...


Very Best Regards!
Edited On 3:15:45 PM PST - Wed, Nov 29th 2017 by Xiongyan21

Forum Regular
Threads 56
Posts 268
RPA NWCHEM6.6 calculation is successful with PBE0 and aug-cc-pvtz which gives
Root 1 a 0.000702294 a.u. 0.0191 eV
 <S2> =   0.9599
...
Root 2 a 0.000810231 a.u. 0.0220 eV
<S2> = 0.9603
... ,
the closest to that of the NIST data.
Edited On 8:00:20 PM PST - Wed, Jan 31st 2018 by Xiongyan21

Forum Regular
Threads 56
Posts 268
There is something happening about this calculation with NWCHEM6.8 on macOS High Sierra with GCC 7.2.0, mpich 3.2.1 and Xcode 9.2.
First, the calculation using aug-cc-pvtz could not converge
...
Entering Davidson iterations
 Unrestricted

 Iter   NTrls   NConv    DeltaV     DeltaE      Time   
---- ------ ------ --------- --------- ---------
1 10 0 0.15E+00 0.10+100 328.3
------------------------------------------------------------------------
tddft_diagon: negative excitation energy 0
------------------------------------------------------------------------.
Second, the one with aug-cc-pvdz was successful for one time
...


 Root   1         a              0.000479490 a.u.                0.0130 eV 
<S2> = 0.9633
----------------------------------------------------------------------------,
...
----------------------------------------------------------------------------
Root 2 a 0.000721485 a.u. 0.0196 eV
<S2> = 0.9615
----------------------------------------------------------------------------
...
----------------------------------------------------------------------------
Root 3 a 0.001720248 a.u. 0.0468 eV
<S2> = 0.8747
----------------------------------------------------------------------------
...


 Root   4         a              0.010212053 a.u.                0.2779 eV 
<S2> = 0.7555
----------------------------------------------------------------------------
...
but failed when the computer was restarted using the same input
...Entering Davidson iterations
 Unrestricted

 Iter   NTrls   NConv    DeltaV     DeltaE      Time   
---- ------ ------ --------- --------- ---------
1 10 0 0.15E+00 0.10+100 151.3
2 30 0 0.79E-01 0.35E-01 205.2
------------------------------------------------------------------------
tddft_diagon: negative excitation energy 0
------------------------------------------------------------------------
------------------------------------------------------------------------
current input line :
29: task tddft energy
------------------------------------------------------------------------
...

Finally, the repeated CIS TDDFT calculations had a small difference with aug-cc-pvdz:
1. ...


 Root   1         a              0.008091830 a.u.                0.2202 eV 
<S2> = 0.9172
----------------------------------------------------------------------------
...
----------------------------------------------------------------------------
Root 2 a 0.008125684 a.u. 0.2211 eV
<S2> = 0.9172
----------------------------------------------------------------------------
...
2. ...


 Root   1         a              0.008077333 a.u.                0.2198 eV 
<S2> = 0.9181
----------------------------------------------------------------------------
...


 Root   2         a              0.008115493 a.u.                0.2208 eV 
<S2> = 0.9181
----------------------------------------------------------------------------
...

GAMESS gave 0.048eV, 0.048eV and 0.284eV,which could be reproduced from different versions.

Do these mean the rohf calculation using NWCHEM was very unstable in this case and a patch is required?

Thanks a lot.

Very Best regards!
Edited On 7:37:44 PM PST - Wed, Jan 31st 2018 by Xiongyan21

Forum Regular
Threads 56
Posts 268
The GAMESS Tamm-Dancoff approximation calculation gave two degenerate states of 0.238 eVs, which agree with the NWCHEM CIS (TDA) calculation results.

It is remarkabe that this calculation has a great difference between the TDDFT results with and without TDA.

In fact, it is told by GAMESS aug-cc-pvdz is an illegal minimal basis function requested for Sc in this calculation.

Today, the RPA TDDFT calculation also fails with NWCHEM6.8.

By the way, it seems that geometry 0.0 0.0 0.0 gets all 0.0eVs results using NWCHEM6.8 for this calculation, which is not the case for GAMESS.
Edited On 7:16:13 PM PST - Wed, Jan 31st 2018 by Xiongyan21

Forum Regular
Threads 56
Posts 268
There is something really unusual for this calcualtion. For the first time, NWCHME6.8 compiled on Ubuntu17.10 was successful and has somewhat different results from those obtained by that compiled on macOS High Sierra, but the repeated ones all failed.
What is the reason here?
Edited On 8:25:16 PM PST - Sun, Jan 21st 2018 by Xiongyan21

Forum Regular
Threads 56
Posts 268
I started a new calculation using the similar input to that of tddft_he, which can add 2nd order and 3rd order derivatives for many of the functionals, but it still failed at the second root calculation, and the first one is much more larger.
Edited On 7:39:35 PM PST - Wed, Jan 31st 2018 by Xiongyan21

  • Niri Forum:Admin, Forum:Mod, NWChemDeveloper, bureaucrat, sysop
    Profile
    Send PM
Forum Regular
Threads 4
Posts 207
Hi Xiongyan,

I am looking into this and will get back to you.

Best,
-Niri Govind

Forum Regular
Threads 56
Posts 268
Dear Dr. Niri
  Have you found what has caused the porblem here?

  Thanks a lot!

  Very Best Regards?
Edited On 3:14:50 AM PDT - Fri, Apr 20th 2018 by Xiongyan21

Forum Regular
Threads 56
Posts 268
A possible conclusion should be drawn as soon as raising the question: the negative excitation energy is irrational and wrong here, becuase it could not be below the reference state here.
Edited On 3:15:00 PM PDT - Wed, Mar 13th 2019 by Xiongyan21


Forum >> NWChem's corner >> Running NWChem



Who's here now Members 0 Guests 1 Bots/Crawler 0


AWC's: 2.5.10 MediaWiki - Stand Alone Forum Extension
Forum theme style by: AWC