Towards hybrid quantum mechanical/molecular mechanical simulations of Li and Na intercalation in graphite - force field development and DFTB parametrisation

In this work a previously established QM/MM simulation protocol for the treatment of solid-state interfaces was extended towards the treatment of layered bulk materials enabling for instance investigation of metal intercalation in graphitic carbon materials. In order to study the intercalation of Li in graphite, new density functional tight binding (DFTB) parameters for Li have been created. Molecular dynamics (MD) simulations at constant temperatures (273.15, 298.15 and 323.15 K) have been carried out to assess the performance of the presented DFTB MD simulation approach. The intercalation of variable lithium and sodium content was investigated via z-distribution functions and analysis of the diffusivity in the direction parallel to the graphene plane. Both the calculated diffusion coefficients and the activation energy in case of lithium are in good agreement with experimental data. The comparison of the QM/MM MD simulation results provide detailed insights into the structural and dynamical properties of intercalated metal ions. The QM/MM simulation protocol, previously developed for solid-state interfaces, has been extended for the investigation of graphitic carbon materials. To study the intercalation of Li in graphite, new DFTB parameters for Li have been formulated.