Molecular dynamics (MD) is a very powerful tool for modeling of molecular systems at normal and elevated temperatures T ≥ 300 K. However, the simulation of nuclear quantum effects (NQEs) is still a challenging topic in MD. At normal temperatures, the path-integral molecular dynamics (PIMD) approach is very effective for modelling equilibrium properties of molecular systems. This technique takes advantage of the exact mapping between the quantum-mechanical system and its representation as a set of N interconnected replicas (also termed as beads) with the N -times elevated temperature. Taking enough beads of this so-called ring polymer, when their respective temperature (T × N ) overcomes the Debye temperature T D = hν/k B of the highest possible vibrational mode with frequency ν, the representation approaches the exact limit. The minimally required number of beads is thus estimated as N min = hν kBT . Therefore at low temperatures, this number is getting very large thus diminishing the applicability of PIMD.