The variety of torsion fitting procedures used in prior studies shows that a standard procedure for carrying out this step is currently lacking. This is in part because each approach involves making a different set of compromises when fitting the approximate potential, and it is challenging to assess the impact of an approach on the accuracy of simulated properties. This and other challenges motivates the creation of automated benchmarking tools that are the subject of current research but beyond the scope of this paper.