<div dir="ltr"><div>Dear Daniel,</div><div><br></div><div>Thank you very much for your 
multiple suggestions, and taking the time to also run them. I had no 
idea it is possible (or indeed advisable) to use the amplitudes of a 
different kind of coupled cluster method, and I don't think the nomp2 
option for OQVCCD is documented in the Molpro manual.</div><div> 
Anyways, I will gather the geometries with QVCCD convergence issues 
(less than 1000 out of ~50k, some have HF convergence problems and other
 weirder failures) and try the zero-vector guess first, since it saves 
the time required by the BDCD calculation.</div><div><br></div><div>Best wishes,</div>Tibor<div class="gmail_extra"><br><div class="gmail_quote">2017-05-23 12:25 GMT+02:00 Daniel Kats <span dir="ltr"><<a href="mailto:dnkats@gmail.com" target="_blank">dnkats@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div><div><div><div>Dear Tibor,<br><br></div>There are few things one can try:<br><br></div>1. The best solution would be to go from a single-reference geometry to the multi-reference geometries in small steps, always restarting from the previous calculation, e.g.,<br>file,2,molecule.wfu<br>.....<br>hf<br>{oqvccd(t),shiftp=0.8,shifts=0<wbr>.6;maxit,100;start,5200.2;save<wbr>,5200.2}<br><br></div>2. If it's not possible, one can try to generate a more reasonable starting guess than MP2. <br>First option is to simply use a zero-vector instead of MP2-amplitudes:<br>{oqvccd(t),shiftp=0.8,shifts=0<wbr>.6,nomp2=1;maxit,100}<br><br></div>3. Another option is to use e.g. distinguishable cluster amplitudes as a starting guess for OQVCCD:<br>{bdcd,nomp2=1,shiftp=0.8,shift<wbr>s=0.8;maxit,100;save,5200.2}<br>{oqvccd(t);maxit,100;start,520<wbr>0.2}<br></div><br>You can find outputs of these last two options attached.<br><br></div>Best wishes,<br></div>Daniel <br><div><div><div><div><div><div><div><div><div><div><div><br><div class="gmail_quote"><div dir="ltr">On Tue, May 23, 2017 at 9:01 AM Tibor Győri <<a href="mailto:tiborgyri@gmail.com" target="_blank">tiborgyri@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Yes, I have tried {oqvccd(t);shift,1,1} and {oqvccd(t),shifts=0.2,shiftp=0<wbr>.2,thrdis=1.0;diis,1,1,15,1;ma<wbr>xit,100;}</div><div>Since I only tried these at a single geometry, it is possible that they would help at the other problematic geometries, but they are no silver bullet.<br></div></div><div class="gmail_extra"><br><div class="gmail_quote">2017-05-22 7:42 GMT+02:00 Peterson, Kirk <span dir="ltr"><<a href="mailto:kipeters@wsu.edu" target="_blank">kipeters@wsu.edu</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Tibor,<br>
<br>
so you tried shifting both the singles and pairs via something like:   {oqvccd(t);shift,1,1}   ?<br>
<br>
best,<br>
<br>
-Kirk<br>
<br>
> On May 21, 2017, at 5:54 AM, Tibor Győri <<a href="mailto:tiborgyri@gmail.com" target="_blank">tiborgyri@gmail.com</a>> wrote:<br>
><br>
> Dear Molpro users,<br>
><br>
> Recently I have run into some issues while trying to use OQVCCD(T) for the purpose of calculating energies for PES fitting. The problem is, that to fit a global PES some calculations have to be run at pretty extreme geometries, and I have seen many examples, where ordinary CCSD and CCSD-F12 converges but OQVCCD does not.<br>
><br>
> The awesome multiple bond breaking characteristics of OQVCCD(T) would be really advantageous for fitting a global PES, but these convergence issues are hindering my efforts. I have tried to increase the maximum number of iterations allowed to 100, but a number of calculations still do not converge.<br>
> I also tried to use level shifts and increased DIIS size at one geometry, but it did not work any better.<br>
><br>
> Is there anything I could do to increase the chance of OQVCCD converging in difficult cases? Since these input files are automatically generated, anything that requires manual per-input tuning is pretty much out of the question.<br>
> I have attached an abridged log file of a calculation that had issues.<br>
><br>
> Best wishes,<br>
> Tibor Győri<br>
> MSc student<br>
> University of Szeged<br>
><br>
> <oqvccd.log>__________________<wbr>_____________________________<br>
> Molpro-user mailing list<br>
> <a href="mailto:Molpro-user@molpro.net" target="_blank">Molpro-user@molpro.net</a><br>
> <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__www.molpro.net_mailman_listinfo_molpro-2Duser&d=DwICAg&c=C3yme8gMkxg_ihJNXS06ZyWk4EJm8LdrrvxQb-Je7sw&r=mqaXWsZ0ghQR-pXvPLjY_g&m=OtWbWWEZBfJvvNsgh34ZzXKGiVsLFXfovJUhURi27qM&s=M6yVxvlkDsAB_09PjJGiboroHxeAxsy3-OuOSeigYfk&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.<wbr>com/v2/url?u=http-3A__www.molp<wbr>ro.net_mailman_listinfo_molpro<wbr>-2Duser&d=DwICAg&c=C3yme8gMkxg<wbr>_ihJNXS06ZyWk4EJm8LdrrvxQb-Je7<wbr>sw&r=mqaXWsZ0ghQR-pXvPLjY_g&m=<wbr>OtWbWWEZBfJvvNsgh34ZzXKGiVsLFX<wbr>fovJUhURi27qM&s=M6yVxvlkDsAB_0<wbr>9PjJGiboroHxeAxsy3-OuOSeigYfk&<wbr>e=</a><br>
<br>
</blockquote></div><br></div>
______________________________<wbr>_________________<br>
Molpro-user mailing list<br>
<a href="mailto:Molpro-user@molpro.net" target="_blank">Molpro-user@molpro.net</a><br>
<a href="http://www.molpro.net/mailman/listinfo/molpro-user" rel="noreferrer" target="_blank">http://www.molpro.net/mailman/<wbr>listinfo/molpro-user</a></blockquote></div></div></div></div></div></div></div></div></div></div></div></div></div><span class="m_1316724306294329320HOEnZb"><font color="#888888"><div dir="ltr">-- <br></div><div data-smartmail="gmail_signature"><div dir="ltr"><span class="m_1316724306294329320m_8619129767938860098inbox-lG">Daniel Kats<br>Institute</span> <span class="m_1316724306294329320m_8619129767938860098inbox-lG">for</span> <span class="m_1316724306294329320m_8619129767938860098inbox-lG">Theoretical</span> Chemistry<br>
University of Stuttgart<br>
Pfaffenwaldring 55<br>
D-70569 Stuttgart</div></div>
</font></span></blockquote></div><br></div></div>