The GEMS method enables molecular dynamics simulations of large heterogeneous systems at ab initio quality.