This code benefits (i) simultaneous inclusion of charge fragments in the CDFT-AIMD calculations to lower simulation wall clock time, (ii) reduction in computational time by efficient construction of the constraint weight function, (iii) solving electronic structure based on orbital transformation 33 method in a three-tiered self-consistent field (SCF) approach, two for energy minimization and one for constraint Lagrangian maximization, and (iv) a hybrid Gaussian and plane-wave basis set to efficiently perform the simulations of low symmetry systems.