If you have problems during the execution of MRCC, please attach the output with an adequate description of your case as well as the followings:
  • the way mrcc was invoked
  • the way build.mrcc was invoked
  • the output of build.mrcc
  • compiler version (for example: ifort -V, gfortran -v)
  • blas/lapack versions
  • as well as gcc and glibc versions

This information really helps us during troubleshooting :)

RS-DH ADC(2) calculation

  • TeichHyperion
  • Topic Author
  • Offline
  • New Member
  • New Member
More
8 hours 25 minutes ago #1669 by TeichHyperion
RS-DH ADC(2) calculation was created by TeichHyperion
Hi everyone,

I have been trying to calculate some excited states with some double excitation character using RS-DH ADC(2). But the calculation keeps crashing for various reasons. I wanted to ask what the recommended input for that type of calculation is and if there is anything specific that is wrong here? I have been playing around with various settings and keyword combos, but it usually fails at the drpa step. 

Sincerely 
Lukas

calc=RS-PBE-P86
dhexc=SCS-ADC(2)
scfiguess=restart

nsing=3
guido_ct=on
nto=on

occri=on                 
basis=def2-SVP
fmm=on
pssp=on

core=frozen
mem=480GB
unit=angs
charge=0
mult=1

geom=xyz
134



 

Please Log in or Create an account to join the conversation.

More
6 hours 36 minutes ago #1670 by mesterdavid
Replied by mesterdavid on topic RS-DH ADC(2) calculation
Dear Lukas,

Thank you for your interest in the method!

You can find example inputs in the MTEST folder, for instance:
MINP_H2O_aug-cc-pVDZ_RS-PBE-P86_SOS-ADC(2)_direct
The main issue with your input is that you should use the combination calc=TDA and dft=RS-PBE-P86 instead of calc=RS-PBE-P86.

Please also consider the following points:
1. SCF initial guess: If you use scfiguess=restart, please make sure that the appropriate guess is available.
2. SCF acceleration: RS-PBE-P86/SCS-ADC(2) calculations are relatively expensive, so the computational cost of the SCF solution is likely negligible compared to the overall calculation. Therefore, I would recommend disabling the keywords intended to accelerate the SCF procedure, such as occri, fmm, and pssp. Their compatibility with excited-state calculations has not been tested. Alternatively, you could perform the ground-state calculation separately and then restart the excited-state calculation with scfiguess=off, ensuring you retain the necessary files from the ground-state calculation.
3. ADC(2) settings and spin scaling: I would recommend removing the dhexc=SCS-ADC(2) keyword, since the ADC(2) ansatz is used by default. However, the spin-scaling factors must be specified explicitly. Unfortunately, this is not entirely clear. According to the paper DOI: 10.1021/acs.jctc.1c01100, the appropriate factors are scsps=1.24 and scspt=0.64.

calc=TDA
dft=RS-PBE-P86
scsps=1.24
scspt=0.64

Please feel free to contact me if you have any further questions or encounter additional issues.

Best regards,
Dávid
The following user(s) said Thank You: TeichHyperion

Please Log in or Create an account to join the conversation.

Time to create page: 0.037 seconds
Powered by Kunena Forum