Single-simulation determination of phase boundaries: A dynamic Clausius–Clapeyron integration method