<div dir="ltr"><p class="MsoNormal" style="font-size:12pt;color:rgb(0,0,0)"><span style="font-size:10pt"><font face="arial, sans-serif">Dear all, </font></span></p><p class="MsoNormal" style="font-size:12pt;color:rgb(0,0,0)"><span style="font-size:10pt"><font face="arial, sans-serif"> </font></span></p><p class="MsoNormal" style="font-size:12pt;color:rgb(0,0,0)"><font face="arial, sans-serif"><span style="font-size:10pt">I'm running a single point calculation on a positively charged rubrene (</span><span style="font-size:10.5pt;border-width:1pt;border-style:none;border-color:windowtext;padding:0cm">C</span>_{<span style="font-size:10pt;border-width:1pt;border-style:none;border-color:windowtext;padding:0cm">42</span>}<span style="font-size:10.5pt;border-width:1pt;border-style:none;border-color:windowtext;padding:0cm">H</span>_{<span style="font-size:10pt;border-width:1pt;border-style:none;border-color:windowtext;padding:0cm">28</span>}<span style="font-size:10pt">) dimer (spin polarized calculation) using PBE+50%HFX, restarting from a converged PBE wavefunction. Unfortunately the SCF does not converge after 350 steps. The "convergence" value jumps to huge values at SCF step 2 ("convergence" at SCF step1 = <b>3.3</b>, "convergence" at SCF step2 = <b>12281</b>) and at longer times (around SCF step 150 onwards), it shows oscillatory behaviour. The behaviour is the same if I use RESTART ATOMIC. </span></font></p><p class="MsoNormal" style="font-size:12pt;color:rgb(0,0,0)"><span style="font-size:10pt"><font face="arial, sans-serif"> </font></span></p><p class="MsoNormal" style="font-size:12pt;color:rgb(0,0,0)"><span style="font-size:10pt"><font face="arial, sans-serif">I note:</font></span></p><ul type="disc" style="color:rgb(0,0,0);font-size:medium"><li class="MsoNormal" style="margin:0cm 0cm 0.0001pt;font-size:12pt;line-height:12.75pt"><span style="font-size:10pt"><font face="arial, sans-serif">The SCF also doesn't converge for the neutral dimer</font></span></li><li class="MsoNormal" style="margin:0cm 0cm 0.0001pt;font-size:12pt;line-height:12.75pt"><span style="font-size:10pt"><font face="arial, sans-serif">The SCF <b>does converges </b>for both the neutral and cationic <b>single rubrene molecule</b></font></span></li><li class="MsoNormal" style="margin:0cm 0cm 0.0001pt;font-size:12pt;line-height:12.75pt"><span style="font-size:10pt"><font face="arial, sans-serif">For the dimer, I have tried a number of settings: different mixing schemes (BROYDEN, DIRECT_P_MIXING), diagonalization and OT method, larger CUTOFF, smaller EPS_SCHWARZ</font></span></li><li class="MsoNormal" style="margin:0cm 0cm 0.0001pt;font-size:12pt;line-height:12.75pt"><span style="font-size:10pt"><font face="arial, sans-serif">I also can't converge the SCF using PBE0, where I observe similar behaviour </font></span></li></ul><p class="MsoNormal" style="font-size:12pt;color:rgb(0,0,0)"><span style="font-size:10pt"><font face="arial, sans-serif">I attach the relevant files. Any help would be much appreciated. </font></span></p><p class="MsoNormal" style="font-size:12pt;color:rgb(0,0,0)"><span style="font-size:10pt"><font face="arial, sans-serif"> </font></span></p><p class="MsoNormal" style="font-size:12pt;color:rgb(0,0,0)"><span style="font-size:10pt"><font face="arial, sans-serif">Best,</font></span></p><p class="MsoNormal" style="font-size:12pt;color:rgb(0,0,0)"><span style="font-size:10pt"><font face="arial, sans-serif"> </font></span></p><p class="MsoNormal" style="font-size:12pt;color:rgb(0,0,0)"><span style="font-size:10pt"><font face="arial, sans-serif">Jan</font></span></p><p class="MsoNormal" style="font-size:12pt;color:rgb(0,0,0)"><span style="font-size:10pt"><font face="arial, sans-serif"><br></font></span></p><p class="MsoNormal" style="font-size: 12pt;"><span style="font-size: 10pt;"><font face="arial, sans-serif" style="" color="#000000"><b><u>EDIT:</u></b></font></span></p><p class="MsoNormal" style=""><font color="#000000"><font face="arial, sans-serif"><span style="font-size: 13.3333px;">I thought I should also add, when I try to use Fermi-Dirac smearing, I get the following errors:  "</span></font><span style="font-family: "Andale Mono"; font-size: 12px; font-variant-ligatures: no-common-ligatures;">WARNING in qs_mo_occupation.F:118 :: Fermi-Dirac smearing includes </span><span style="font-variant-ligatures: no-common-ligatures; font-family: "Andale Mono"; font-size: 12px;">the first MO" </span><span style="font-family: arial, sans-serif; font-size: 13.3333px;">and <span style="background-color: rgb(255, 255, 255);">"</span></span><span style="font-variant-ligatures: no-common-ligatures; font-family: "Andale Mono"; font-size: 12px; background-color: rgb(255, 255, 255);">WARNING in qs_mo_occupation.F:122 :: Fermi-Dirac smearing includes </span></font><span style="color: rgb(0, 0, 0); font-variant-ligatures: no-common-ligatures; font-family: "Andale Mono"; font-size: 12px;">the last MO => Add more MOs for proper smearing.</span><span class="Apple-converted-space" style="color: rgb(0, 0, 0); font-variant-ligatures: no-common-ligatures; font-family: "Andale Mono"; font-size: 12px;"> "</span></p><p class="MsoNormal" style=""><br></p><p class="MsoNormal" style=""><span style="color: rgb(0, 0, 0); font-family: arial, sans-serif; font-size: 13.3333px;">This occurs even when I use a tiny Electronic Temperature (0.1 K) and 100 added MOs for each spin channel (files attached).</span><br></p><p class="MsoNormal" style="font-size:12pt;color:rgb(0,0,0)"><span style="font-size:10pt"><font face="arial, sans-serif"> </font></span></p></div>