to show statistically significant deviations, suggesting a gradual erosion of accuracy as the timestep increases. These significant p-values at 4 fs, and increasingly at 2.5 fs and 3.3333 fs, highlight the increased risk of inaccuracies in calculated thermodynamic properties when employing larger timesteps, reinforcing the observation of increasing dU dλ magnitude with increasing timestep. The t-test results provide strong statistical support for the conclusion that smaller timesteps (0.5 fs, 1 fs, and 2 fs) offer more consistent and reliable estimates of dU dλ compared to the larger timesteps, especially 4 fs, which should likely be avoided for accurate thermodynamic property calculations.