We present a computationally efficient general first-principles based method for spin-lattice simulations for solids and clusters. The method is based on a coupling of atomistic spin dynamics and molecular dynamics simulations, expressed through a spin-lattice Hamiltonian, where the bilinear magnetic term is expanded up to second order in displacement. The effect of first order spin-lattice coupling on the magnon and phonon dispersion in bcc Fe is reported as an example, and we observe good agreement with previous simulations. In addition, we also illustrate the coupled spin-lattice dynamics method on a more conceptual level, by exploring dissipation-free spin and lattice motion of small magnetic clusters (a dimer, trimer and quadmer). The here discussed method opens the door for a quantitative description and understanding of the microscopic origin of many fundamental phenomena of contemporary interest, such as ultrafast demagnetization, magnetocalorics, and spincaloritronics. arXiv:1804.03119v2 [cond-mat.mtrl-sci] 9 Jul 2018