optimized_configuration = OptimizeGeometry(
configuration=molybdenite,
constraints=[FixStrain(False, False, True), SpaceGroupConstraint()],
trajectory_filename='Molybdenite_results.hdf5',
trajectory_object_id='optimize_trajectory',
optimize_cell=True,
restart_strategy=restart_strategy,
)