<p>Dear Prof. Lannuzzi</p>
<p>The quintet state has 4 singly occupied alpha spin orbitals. I tried to flip two of these alpha spin orbitals to beta spin. therefore, I set the "RESULT_MO_INDEX" as negtive to represent the unoccupied MO space. I actually also tries the opposite way: setting ORIG_MO_INDEX as the virtual MO and RESULT_MO_INDEX as HOMO and HOMO-1. This run was "successful", however, I am not sure I understand the results. The output file shows that "Performing wfn mixing and Writing new MOs to file" at the end of the calculation. In this specific calculation, the spin multiplicity is 5 and I used the ATOMIC guess, the program seems to perform wfn mixing after SCF for the initial quintet state (although not converge), How do i know the wfn mixing is correctly performed? Does the saved wfn file contains the spin-fliped singlet state? Could you please provide some suggestions? Thank you! Attached please find the input and output.</p>
<p> </p>
<p>Best wishes</p>
<p> </p>
<p>John</p><br /><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Thursday, October 26, 2023 at 4:54:26 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 Peng<div><br></div><div>But why have you set ORIG_IS_VIRTUAL as FALSE and given negative indexes to the MOs?</div><div><br></div><div>Regards</div><div>Marcella<br><div><br></div></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Tuesday, October 24, 2023 at 6:37:10 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">Dear Prof. Lannuzzi<div><br></div><div>The previous quintet calculation did not have the "ADDED_MOS". Using ATOMIC guess, I have MO_CUBES in the PRINT section of DFT, but got the same error. Please see attached outputs (I did not include the generated cube files). I am using 2023.2. Maybe my compilation have some issues? or should I try other version of cp2k? Thanks!</div><div><br></div><div>Best wishes</div><div><br></div><div>Peng</div><div><br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Tuesday, October 24, 2023 at 11:45:52 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"><div><br></div><div>Dear John, </div><div><br></div><div>The restart does not work if the previous calculation had ADDED_MOS, but the ATOMIC guess should work. </div><div>You have to keep the MO_CUBES in the PRINT section of DFT<br></div><div>It works for me. </div><div><br></div><div>Regards</div><div>Marcella</div><div><br></div><br><br><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Tuesday, October 24, 2023 at 4:51:09 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">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></blockquote></div></blockquote></div></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/5c1ea429-64f7-4985-ab69-ef5a34c00210n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/5c1ea429-64f7-4985-ab69-ef5a34c00210n%40googlegroups.com</a>.<br />