  The quantum thermal bath (QTB) has been presented as analternative to path-integral-based methods to introduce nuclear quantumeffects in molecular dynamics simulations. The method has proved to beefficient, yielding accurate results for various systems. However, the QTBmethod is prone to zero-point energy leakage (ZPEL) in highly anharmonicsystems. This is a well-known problem in methods based on classicaltrajectories where part of the energy of the high-frequency modes istransferred to the low-frequency modes leading to a wrong energydistribution. In some cases, the ZPEL can have dramatic consequences onthe properties of the system. Thus, we investigate the ZPEL by testing theQTB method on selected systems with increasing complexity in order to studythe conditions and the parameters that influence the leakage. We also analyze the consequences of the ZPEL on the structuraland vibrational properties of the system. We find that the leakage is particularly dependent on the damping coefficient and thatincreasing its value can reduce and, in some cases, completely remove the ZPEL. When using sufficiently high values for thedamping coefficient, the expected energy distribution among the vibrational modes is ensured. In this case, the QTB methodgives very encouraging results. In particular, the structural properties are well-reproduced. The dynamical properties should beregarded with caution although valuable information can still be extracted from the vibrational spectrum, even for large values ofthe damping term. 
Zero-Point Energy Leakage in Quantum Thermal Bath Molecular Dynamics   Simulations