Accurate and efficient machine learning interatomic potentials for finite temperature modeling of molecular crystals