In this paper, We will combine the direct NVE scheme with both conventional DFT-based AIMD and CNEO-MD. The CNEO-MD framework was recently developed by our group to incorporate nuclear quantum effects, especially zero-point effects due to quantum delocalization, in practical molecular simulations. Its working equations are m d⟨x⟩ dt = ⟨p⟩ (10) and