For our first application of the MiMiC framework, we have coupled the CPMD 27 and GROMACS codes to enable electrostatic embedding QM/MM-MD. GROMACS was chosen based on its high performance, versatility, and the availability of different force fields, as well as its active development community which ensures long-term support. The use of the CPMD code as the main driver enables MD within the CP and BO variants, with efficient orbital extrapolation schemes to improve the computational efficiency for the latter, using a plane wave/pseudopotential implementation of KS-DFT. Density functional perturbation theory calculations can be performed, and electronically excited states can be accessed within the linear response or real-time time-dependent DFT (TDDFT) frameworks as well as the restricted open-shell KS (ROKS) formalism. Nonadiabatic effects 40 can be included either within an Ehrenfest scheme or by virtue of trajectory surface hopping (TSH), while nuclear quantum effects can be taken into account within a path integral (PI) formalism. The code is highly parallelized, offers access to a large number of stateof-the-art exchange-correlation functionals and a computationally efficient treatment of hybrid functionals, and allows for long time scale simulations by virtue of efficient multiple time step (MTS) acceleration algorithms. The choice of CPMD therefore allows for optimal flexibility and high efficiency in the description of the QM subsystem.