Why does using deuterium mass for hydrogen improve convergence of molecular dynamics simulations of water? I've come across many papers performing ab initio MD (Carr-Parinello or Born-Oppenheimmer) simulations of water that use the deuterium mass for hydrogen. And they state that this helps improve convergence. I haven't seen a derivation for it and I'm not able to come up with an intuitive reason for why this might be. I would greatly appreciate some insights.
Here's an example of a paper that does so (I'm linking this because it doesn't seem to need institutional credentials to access): https://www.spiedigitallibrary.org/journals/Journal-of-Photonics-for-Energy/volume-1/issue-1/016002/iAb-initio-i-modeling-of-water-semiconductor-interfaces-for-photocatalytic/10.1117/1.3625563.full?SSO=1 
Thanks! 
 A: My intuitive reason is based on the need for time-scale separation between the evolution of the electronic degrees of freedom and the nuclear motion. 
The electrons are given fictitious "masses" for the evolution of the quantum degrees of freedom, and are typically maintained in their ground state by a low-temperature thermostat. The nuclei move around at speeds consistent with the physically desired temperature (e.g. ambient temperature). This means that, technically, the system is not at equilibrium. However, energy flow between the two subsystems is kept under control if their natural timescales are sufficiently different: usually called "adiabaticity". 
Giving the lightest nuclei an inflated mass is a simple way of slowing them down. It is used to allow a longer timestep in purely classical MD; but here I think that it also makes it easier to give the electrons a correspondingly higher mass,  without compromising the adiabaticity, and hence use a longer timestep in the quantum part of the calculation (which is the expensive part).
