<div dir="ltr"><br>Hi everyone,<br><br>I am trying to implement quadrupole damping and anion deformations for classical force field calculations and right now I am just testing what is already working.<br><br>While doing so I realised that as soon as I set the MAX_MULTIPOLE_EXPANSION to Quadrupole in &POISSON &EWALD <br>and define Quadrupole polraizabilities (CPOL),  the conjugate_gradient method for the polarization cycle doesn't seem to work anymore.<br><br>Putting POL_SCF to self_consistent reproduces the forces of the reference, conjugate_gradient however is far of (see datasets attached).<br>However reducing to Dipole expansion and and deleting cpol terms reproduces the respective reference just fine, no matter the algorithm (as it should be, I guess?)<br>What could be the matter here?<br><br>I am using the latest svn checkout.<br><br>Attached the input files, the two sample outputs (forces, one using conjugate_gradient, one using self consistent), and the reference forces.<br><br><br>The only section I am changing at all is<br><br>   &EWALD <br>        ...<br>     &MULTIPOLES T<br>         MAX_MULTIPOLE_EXPANSION QUADRUPOLE<br>         <span style="color: rgb(255, 153, 0);">POL_SCF SELF_CONSISTENT</span><br>         EPS_POL 1.0e-6<br>         MAX_IPOL_ITER 100<br>       &END MULTIPOLES<br>     &END EWALD<br><br><br>Thanks very much for every hint and help! <br>Johannes<br><br></div>