where λ serves as a coupling parameter that smoothly transitions between the initial (λ=0) and final states (λ=1). In the second integral, the derivative of the potential energy U with respect to λ is calculated for every snapshot of the simulations, with the averaging brackets denoting the equilibrium sampling of the state corresponding to the λ parameter.