× 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 :)

"Invalid spatial symmetry!" problem in CCSD(T)

More
7 years 3 months ago #205 by mradon
Dear All:

I am trying to perform CCSD(T) calculations for an open-shell system based on UHF reference under the D2h symmetry. In my MINP file, I used the following keywords to define spin multiplicity, occupation numbers, and spatial symmetry of the required electronic state: mult, symm, occ. I am pretty sure that I used these keywords properly because the SCF converged to the desired electronic state and the energy is the same as computed with other program using the same basis set and geometry.

However, when going to the CC part, the progam gives the following error message and quits:

CC(2)(3) calculation
(...)
Number of alpha electrons: 28
Number of beta electrons: 24
Spin multiplicity: 2
z-component of spin: 2.0
(...)
Invalid spatial symmetry!

I do not understand why the spatial symmetry is invalid for the CC because it is correctly specified and the same input was OK for the SCF calculations. Moreover, the numbers of alpha and beta electrons printed above is wrong and inconsistent with the ones used during the SCF. Hence, the multiplicity is wrong: 2 should be 6 (five unpaired electrons), the z-component of spin is 2.0, should be 2.5 (BTW, the spin multiplicity 2 and z-component of spin 2.0 are inconsistent with each other).

My question is then: how can I pass the information about spin multiplicity and spatial symmetry to the CC program appropriately? Apparently, the keywords: mult, symm, and occ which work for the SCF part are not respected by the CC part. Should I use some additional keywords to define my wave function properly?

Thank you in advance for any suggestions/ comments!

With best regards,
Mariusz

--
Dr Mariusz Radon, Ph.D.
Coordination Chemistry Group
Faculty of Chemistry
Jagiellonian University
ul. Ingardena 3, 30-060 Krakow, Poland
www2.chemia.uj.edu.pl/~mradon

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

More
7 years 3 months ago #206 by kallay
Dear Mariusz,
In principle, the reference determinant is the HF determinant, so I don't understand what happens here. If you post the input file, I will take a look at it.
Anyway, you can set the reference determinant using the refdet keyword, that should work even in this case.

Best regards,
Mihaly Kallay
The following user(s) said Thank You: mradon

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

More
7 years 3 months ago #207 by mradon
Thank you! Your solution with refdet works fine.

Best regards,
Mariusz RadoĊ„

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

Time to create page: 0.018 seconds
Powered by Kunena Forum