<p>Hi</p><p>I want to obtain an open-shell singlet state wave function for a C2 molecule using the WFN_MIX function. I first performed a quintet state calculation (converged smoothly) and tried to use this quintet state wave function as the initial guess in the following open-shell singlet state calculation. In the WFN_MIX section, I tried to flip the spin of the original HOMO and HOMO-1 orbitals (Also see the attached input)</p>
<p>    &WFN_MIX</p><p>      &UPDATE</p>
<p>        RESULT_MO_INDEX    -1</p>
<p>        RESULT_SPIN_INDEX  BETA</p>
<p>        RESULT_SCALE       0.00000</p>
<p><br /></p>
<p>        ORIG_MO_INDEX      2</p>
<p>        ORIG_SPIN_INDEX    ALPHA</p>
<p>        ORIG_SCALE         1.0000</p>
<p><br /></p>
<p>        ORIG_IS_VIRTUAL    .FALSE.</p>
<p>      &END UPDATE</p>
<p><br /></p>
<p>      &UPDATE</p>
<p>        RESULT_MO_INDEX    -2</p>
<p>        RESULT_SPIN_INDEX  BETA</p>
<p>        RESULT_SCALE       0.00000</p>
<p><br /></p>
<p>        ORIG_MO_INDEX      1</p>
<p>        ORIG_SPIN_INDEX    ALPHA</p>
<p>        ORIG_SCALE         1.0000</p>
<p><br /></p>
<p>        ORIG_IS_VIRTUAL    .FALSE.</p>
<p>      &END UPDATE</p>
<p>    &END WFN_MIX</p>
<p><br /></p>
<p>but the calculation failed with an error message "If ORIG_IS_VIRTUAL is set, the array unoccupied_orbs must be associated!"</p>
<p>I think my usage of WFN_MIX is not correct. Could someone please help on this issue? Thanks!<br /></p><p><br /></p><p>Best </p><p><br /></p><p>John</p><p><br /></p>

<p></p>

-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/6700672a-8114-43ea-ae38-7b399192f275n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/6700672a-8114-43ea-ae38-7b399192f275n%40googlegroups.com</a>.<br />