<div dir="ltr">Hi Noam,<div><br></div><div>The original files I attached earlier contained no constraints (fixed atoms or intramolecular).</div><div><br></div><div>What I have noticed is that for AD_LANGEVIN thermostat with:</div><div><br></div><div>1) REGION GLOBAL: Does not work in QM, MM  or QM/MM system (with or without constraints). For pure QM, no constraints (H2O-32.inp) I get a cholesky decomposition error which is gone once I change the thermostat settings.</div><div><br></div><div>2) REGION MASSIVE: Fully functional.</div><div><br></div><div>Since we need to use a global thermostat for fixed atoms, so effectively AD_LANGEVIN doesn't seems to work with constraints and/or REGION GLOBAL thermostat setting.</div><div><br></div><div>Best,</div><div>Mayank</div><div> <br><br>On Wednesday, May 6, 2020 at 6:00:13 PM UTC+2, Noam Bernstein wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div><div dir="auto">It’s quite possible that AD_LANGEVIN doesn’t handle constrained atoms correctly. I certainly didn’t have them in mind when I coded it up. Can you confirm that it works OK with REGION GLOBAL if there are no fixed atoms?</div></div><div dir="auto"><br></div><div dir="auto">    Noam</div><div dir="auto"><br></div>
</blockquote></div></div>