In theory, dH/dL (or dV/dL in thermodynamic integration) originated from the difference in enthalpies of two distinct states lambda and lambda+dL. I wonder how the Gromacs module calculates the value of dH/dL for each lambda point by performing only one simulation at lambda.