Dear Dr Lannuzzi<div><br /></div><div>Thank your for your response!</div><div><br /></div><div>Following you suggestions, I performed a few test calculations (without the ADDED_MOS)</div><div><br /></div><div>1. Starting from the ATOMIC guess with singlet spin multiplicity. The calculation stopped at the beginning of "Performing wfn mixing" and the error message looks like (I pasted only a few lines) </div><div>"PBLAS ERROR 'Array subscript out of bounds: JX = 7, DESCX[N_] = 4'</div><div>from {0,7}, pnum=7, Contxt=0, in routine 'PDCOPY'.</div><div>PBLAS ERROR 'Parameter number 4 had an illegal value'</div><div>from {0,7}, pnum=7, Contxt=0, in routine 'PDCOPY'.</div><div>{0,7}, pnum=7, Contxt=0, killed other procs, exiting with error #-4.</div><div>"</div>
<div><br /></div><div>2. if the calculation started using the quintet wave function, it immediately stopped and the error message is "Number of occupied MOs on restart unit larger than allocated MOs"</div><div><br /></div><div>I also checked the examples in the tests, but I cannot figure out what the problem is. Any insight? Thank you again!</div><div><br /></div><div>Best wishes</div><div><br /></div><div>John</div>
<div class="gmail_quote"><div dir="auto" class="gmail_attr">On Tuesday, October 24, 2023 at 7:10:30 AM UTC-4 Marcella Iannuzzi wrote:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Dear John<div><br></div><div>I think that your input might work if you remove the key ADDED_MOS from SCF.</div><div>Anyway, the input structure for wfn_mix has slightly changed in most recent version of the code. </div><div><br></div><div>Regards</div><div>Marcella</div><div><br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Saturday, October 21, 2023 at 7:52:11 PM UTC+2 John Zhang wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<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></blockquote></div></blockquote></div>
<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/6fe07914-b266-4fa4-a97f-2da42b9b0f0dn%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/6fe07914-b266-4fa4-a97f-2da42b9b0f0dn%40googlegroups.com</a>.<br />