- Posts: 5
- Thank you received: 0
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:
This information really helps us during troubleshooting
- 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
Switching to another root during CCSD calculation of N2 triplet
- Jos Suijker
- Topic Author
- Offline
- New Member
Less
More
3 years 1 month ago #1163
by Jos Suijker
Switching to another root during CCSD calculation of N2 triplet was created by Jos Suijker
Dear Forum,
We calculate the energy of the lowest triplet of N2 by MCSCF followed by the CCSD calculation. Using the patches posted on the MRCC website convergence of my calculations require much less iterations then in the un-patched version and we get convergence also at longer bondlength's. In the un-patched version the calculation does not converge at longer N-N bond length's.
It seems however that the calculation switches to another root during its iterations. for example at bond length R1=1.7775 The calculation is almost converged in iteration 23 to an energy that seems right to me (comparing with the energy found at a little smaller bond length) at iteration 24 however the solution makes a jump and converges to this root. (output is attached). The same effect is also visible in the un-patched version although there is no convergence in that case
Is there a setting that might prevent this jump?
Thanks in advance for any help.
Jos Suijker
Technical University of Eindhoven
The Netherlands
We calculate the energy of the lowest triplet of N2 by MCSCF followed by the CCSD calculation. Using the patches posted on the MRCC website convergence of my calculations require much less iterations then in the un-patched version and we get convergence also at longer bondlength's. In the un-patched version the calculation does not converge at longer N-N bond length's.
It seems however that the calculation switches to another root during its iterations. for example at bond length R1=1.7775 The calculation is almost converged in iteration 23 to an energy that seems right to me (comparing with the energy found at a little smaller bond length) at iteration 24 however the solution makes a jump and converges to this root. (output is attached). The same effect is also visible in the un-patched version although there is no convergence in that case
Is there a setting that might prevent this jump?
Thanks in advance for any help.
Jos Suijker
Technical University of Eindhoven
The Netherlands
Attachments:
Please Log in or Create an account to join the conversation.
- kallay
- Offline
- Administrator
- Mihaly Kallay
3 years 1 month ago #1164
by kallay
Best regards,
Mihaly Kallay
Replied by kallay on topic Switching to another root during CCSD calculation of N2 triplet
Dear Jos,
You should try root following, that is, set diag=follow in the input.
You should try root following, that is, set diag=follow in the input.
Best regards,
Mihaly Kallay
The following user(s) said Thank You: Jos Suijker
Please Log in or Create an account to join the conversation.
- Jos Suijker
- Topic Author
- Offline
- New Member
Less
More
- Posts: 5
- Thank you received: 0
3 years 1 month ago #1165
by Jos Suijker
Replied by Jos Suijker on topic Switching to another root during CCSD calculation of N2 triplet
Dear Prof. Kallay
Thank you very much, this solution worked perfect. The calculation converged easily last night. We should have found this keyword ourselves in the manual it is there.
Regards,
Jos Suijker
Thank you very much, this solution worked perfect. The calculation converged easily last night. We should have found this keyword ourselves in the manual it is there.
Regards,
Jos Suijker
Please Log in or Create an account to join the conversation.
Time to create page: 0.038 seconds