where ρ * is the reweighed density, β is equal to 1/k B T , V t i is the bias in torsion t, D t=1 V t i represents the sum of the concurrent biases acting on the D torsions, for configuration i. Practically, we evaluate ρ, the biased density, from configurations generated in a quasistatic bias regime, as the bulk of the bias is deposited during the early stages of the simulation and the bulk of sampled configurations are visited when bias deposition is negligible. We therefore apply the final bias approximation to obtain a time-independent value of V t i acting on configuration i. To mitigate the noise introduced by exponential reweighting , the density of each point is then averaged over hyperspherical domains of radius 0.1 rad. This step generates a new smoothed set of densities ρ * i , at the cost of a small controllable loss in spatial resolution. The free energy F i , associated with configuration i (thus termed per point), is computed as