Author Topic: CCSD on triplet excited state  (Read 982 times)

garix

  • Newbie
  • *
  • Posts: 6
  • Karma: +0/-0
CCSD on triplet excited state
« on: October 20, 2020, 04:26:47 PM »
Dear all,

I would like to run a CCSD calculation for the triplet excited state, however, the output file provides me with the following error, just after the optimization of the GS state amplitude:

  ----------------------------------------------------------------------------

    Check for potential start vectors for CCSD eigenvalue problem:

     type  irrep  state  model     file
    ------------------------------------------------------------------------
 Illegal state input for CCSD     excitation energies:
 irrep, column:                     1                     1
 multiplicity :                     3
 charge diff. to GS:                     0

========================
 internal module stack:
------------------------
    ccsdf12
    ccsdrspim
========================

 Illegal state input in ccsdrspim!
 ccsdf12 ended abnormally


I ran the same calculation for the singlet and I did not go through any problem.
The part of the control file regarding the excited state calculation is the following:

$ricc2
  ccsd
$excitations
   irrep=a multiplicity=3 nexc=1


Thanks for your help.

Gaetano



Marek Sierka

  • Developers
  • Jr. Member
  • *
  • Posts: 18
  • Karma: +0/-0
    • Sierkalab
Re: CCSD on triplet excited state
« Reply #1 on: October 21, 2020, 09:46:07 AM »
Do you run open shell or closed shell calculation?

garix

  • Newbie
  • *
  • Posts: 6
  • Karma: +0/-0
Re: CCSD on triplet excited state
« Reply #2 on: October 21, 2020, 03:18:55 PM »
I tried both RHF and UHF reference. The former provided me with that error. The latter stopped just after the calculation of the CCS excitation energies, without showing any error:

======== LINEAR CC RESPONSE SOLVER ========

    threshold for convergence: 1.00E-06
    maximum number of simultaneous Jacobi matrix transformations:  1

     summary of start vectors generation:
     -------------------------------------------
     type of solution vectors      : RE0
     symmetry                      : a
     number of vectors requested   :     1
     number of vectors generated   :     1
     -------------------------------------------


    Iter  #Vectors  time (min)    max. residual


Marek Sierka

  • Developers
  • Jr. Member
  • *
  • Posts: 18
  • Karma: +0/-0
    • Sierkalab
Re: CCSD on triplet excited state
« Reply #3 on: October 21, 2020, 04:56:25 PM »
I guess you cannot calculate triplet excited states with RHF reference.  For the UHF reference, did you get any error message in the log file?

garix

  • Newbie
  • *
  • Posts: 6
  • Karma: +0/-0
Re: CCSD on triplet excited state
« Reply #4 on: October 21, 2020, 05:07:22 PM »
No, no error. I tried to rerun the UHF reference calculation increasing the memory. So far, the only error is ccsdf12 step seems to have been in serious trouble but the calculation keeps running. Let's see.

Thank you

Arnim

  • Developers
  • Sr. Member
  • *
  • Posts: 239
  • Karma: +0/-0
Re: CCSD on triplet excited state
« Reply #5 on: October 22, 2020, 12:05:40 PM »
Hello,

at this time (version 7.5) the implementation is currently restricted to closed–shell singlet excited states.

And indeed tripet excitations can (in principle) be calculated from RHF wave functions.

Best
Arnim

garix

  • Newbie
  • *
  • Posts: 6
  • Karma: +0/-0
Re: CCSD on triplet excited state
« Reply #6 on: October 30, 2020, 11:57:40 AM »
Thank you very much for your help guys.

Best,

Gaetano