Therefore, analyzing the derivative of potential energy with respect to the coupling parameter, dU dλ , is essential for accurate and precise RBFE outcomes from alchemical simulations. By exploring this derivative, researchers can evaluate simulation reliability and identify areas for enhancement in their alchemical methodologies. In this study, we specifically investigated the influence of timestep selection on alchemical simulations by thoroughly examining the reproducibility of dU dλ under various timestep conditions.