7.37. Excited States via DLPNO-STEOM-CCSD

7.37.1. PNO dressing (experimental keyword)

In STEOM-DLPNO-CCSD method most of the steps are done using the powerful PNO approximation except the last one corresponding to the STEOM-CCSD calculation itself. In the canonical version of the dressing (default in ORCA) all the amplitudes previously computed at the PNO level are recanonicalized which increases the calculation cost.PNO dressing uses some of the PNO intermediates to reduce the dressing time. It has however no effect on the storage size.

Using this PNO dressing will reduce the calculation timing but at a reduce accuracy for the excitation energies. We thus recommend this option only for specific large systems (>2500 basis functions).

The keyword \(\textbf{DirectDressSTEOM}\) set to true enable this option.

7.37.2. Keywords from STEOM-CCSD

Most of the keywords from STEOM-CCSD can be used within the DLPNO version, except for core excitation. More information on the Excited States via STEOM-CCSD section.

7.37.3. Tips and Tricks

As written in the typical calculation section of the manual we would recommend this input for standard organic molecules.

As a general guideline we will discuss some of the keywords used there:

! STEOM-DLPNO-CCSD def2-TZVP def2-TZVP/C def2/J TightSCF

%mdci
  NRoots 6
  DoRootWise true
  OThresh 0.005
  VThresh 0.005
  TCutPNOSingles 1e-11
  NDAV 400
  DoStoreSTEOM true
  DoSimpleDens false
  AddL2Term True
  DTol 1e-5
end

* xyz 0 1
  C     0.016227   -0.000000    0.000000
  O     1.236847    0.000000   -0.000000
  H    -0.576537    0.951580   -0.000000
  H    -0.576537   -0.951580   -0.000000
*

  • TIGHTSCF is a must for any CCSD calculation.

  • We will recommend using TIGHTPNO for all molecules as it is not a lot more expensive and helps achieving a better convergence.

  • The OTHRESH, VTHRESH and TCutPNOSingles keywords help with converging the calculations, increasing the percentage active of each root. In contrary to standard STEOM-CCSD, we would acknowledge that the roots are converged when the percentage active character is at least 96%. Of course you have to check that the amplitude and orbitals associated with the excitation are correct (and what you expect). Tightening the 3 keywords mentioned will increase this percentage active character. The most important being the TCutPNOSingles one. Be careful, the computational cost increases exponentially when tightened. In more case \(10^{-12}\) or \(5\times 10^{-12}\) should be enough. For OTHRESH and VTHRESH you should not go below \(10^{-3}\), as the benefits are not so obvious.A trick to achieve a better convergence is to play with the number of roots. It is often not necessary to compute a lot of roots if you are only interested in the first 3 for example. If some high energy roots have some low percentage active character, removing them can help for the convergence of other roots.

  • As shown in the related paper[298], the STEOM density with the L2term is now the default for computing the excited state properties.

  • The choice of the basis set can also speed up the calculation without a significant loss of accuracy. For most organic molecules def2-TZVP(-f) is enough. Trying def2-SVP is also a good idea for preliminary tests.

  • The STEOM-CCSD excitation energies are very dependent from the starting geometry. Geometries optimised with several DFT functional can yield significant differences for the excitation energies (about 0.1 eV).

  • It has happened that some roots were missing with the keyword \(\textbf{Dorootwise}\) set to true. Turning it off solved the issue but this keyword should be on by default.