The quality of 1DHR models depends on various factors: accurate modeling of the moment of inertia of the rotating tops, treatment of the harmonic oscillator that is to be removed in exchange for this hindered rotation (HR), and, most important, a good Fourier fit to an accurate 1D PES. Computing a sufficiently accurate PES can be challenging: methods that compute energies well within chemical accuracy can fail to determine the correct energetical order of conformations on a torsional PES, some examples of which we decribe in the Results section of this study. Because of limited computational resources, highest-level methods that would compute such details reliably cannot be applied throughout to the whole PES. Instead, computational methods at different levels provide different compromises along the uncertainty-efficiency trade-off. For example, force fields can compute energies and gradients within split seconds, with uncertainties of dozens of kJ/mol. Density functional theory (DFT) methods yield energies, geometries and frequencies within the order of seconds or hours; the actual computational time usually scales with the 3rd or 4th power of system size, which makes them still applicable for large systems. Their uncertainties vary between some kJ/mol to outliers or known systematic deviations beyond 20 kJ mol -1 . Gold standard coupled-cluster (CC) methods yields energies of chemical accuracy, 1 kcal mol -1 to 1 kJ mol -1 , often slightly deteriorated but accelerated by some computational tricks (Resolution of the Identity, domain-based local pair natural orbital (DLPNO) ) -yet, they remain applicable to systems with a few dozen atoms.