[CP2K-user] [CP2K:16659] Ignored Convergence Criteria For UKS Triplet Calculation

Jürg Hutter hutter at chem.uzh.ch
Thu Mar 3 12:43:39 UTC 2022


Hi

your SCF didn't converge. The CG stepsize was zero and so the energy change also was zero.
There are many problems in your input that might have caused this behavior.
I have attached an updated version (but not tested!) that I would have used as a stating
point. Please fine tune all the changes to your needs.

regards

Juerg Hutter

________________________________________
From: cp2k at googlegroups.com <cp2k at googlegroups.com> on behalf of Josh W <jwheels91 at gmail.com>
Sent: Wednesday, March 2, 2022 10:33 PM
To: cp2k
Subject: [CP2K:16657] Ignored Convergence Criteria For UKS Triplet Calculation

Hi all,

I'm trying to run a single point calculation of a triplet organometallic complex. Under my DFT section I use the following keywords:

        &DFT
                BASIS_SET_FILE_NAME BASIS_MOLOPT
                POTENTIAL_FILE_NAME GTH_POTENTIALS
                UKS TRUE
                MULTIPLICITY 3
                ....
                ....
        &END

Although I have set EPS_DEFAULT to 1E-7, the change in energy during my SCF cycle goes down to about 1E-14 but the calculation continues regardless. If I change my multiplicity to 1, then the calculation has no problems (except that now I have a singlet). Normaly I would just take the energy from the last SCF cycle, but I plan on performing geometry optimizations in the future and would like to have CP2K to move onto the next step once the iteration has reached convergence.

I have attatched my input and output, I'm using ver 7.1.


On another note, I came across this in the manual:

Section COUPLING<https://manual.cp2k.org/cp2k-7_1-branch/index.html#CP2K_INPUT/FORCE_EVAL/MIXED/COUPLING.html>

  *   Coupling between two force_eval: E=(E1+E2 - sqrt((E1-E2)**2+4*H12**2))/2

  *   Section path: CP2K_INPUT<https://manual.cp2k.org/cp2k-7_1-branch/CP2K_INPUT.html> / FORCE_EVAL<https://manual.cp2k.org/cp2k-7_1-branch/CP2K_INPUT/FORCE_EVAL.html> / MIXED<https://manual.cp2k.org/cp2k-7_1-branch/CP2K_INPUT/FORCE_EVAL/MIXED.html> / COUPLING<https://manual.cp2k.org/cp2k-7_1-branch/CP2K_INPUT/FORCE_EVAL/MIXED/COUPLING.html>

Is there a reference for the above equation? Is H12 defined anywhere?

Thanks!

--
You received this message because you are subscribed to the Google Groups "cp2k" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+unsubscribe at googlegroups.com<mailto:cp2k+unsubscribe at googlegroups.com>.
To view this discussion on the web visit https://groups.google.com/d/msgid/cp2k/5c17a847-c688-4f98-b4a0-a2a554b03488n%40googlegroups.com<https://groups.google.com/d/msgid/cp2k/5c17a847-c688-4f98-b4a0-a2a554b03488n%40googlegroups.com?utm_medium=email&utm_source=footer>.

-- 
You received this message because you are subscribed to the Google Groups "cp2k" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+unsubscribe at googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/cp2k/ZR0P278MB0759A55FA081491085861D529F049%40ZR0P278MB0759.CHEP278.PROD.OUTLOOK.COM.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: input.inp
Type: application/octet-stream
Size: 3629 bytes
Desc: input.inp
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20220303/5f568ab5/attachment.obj>


More information about the CP2K-user mailing list