Molecular dynamics (MD) simulations are well positioned to elucidate the aspects of electrospray ionization (ESI) and high-energy collision dissociation (HCD), as well as give insight into processes that involve neutral species that cannot be observed experimentally in ESI, HCD, and collision-induced dissociation (CID). Here, we utilize temperature dissociation molecular dynamics (TDMD) to model the HCD/CID of lithium formate clusters carrying a single positive charge. These simulations successfully reproduce the experimental ESI HCD spectra of lithium formate solutions and also support the existence of magic number clusters (MNCs) that have been observed. The simulations also provide strong evidence that the main fragmentation channel of such clusters involves neutral (LiHCOO)(2) dimers.