The relaxation equations are themselves dependent on the calculation of the spectral density values *J*(*ω*).
Within model-free analysis these are modelled by the original model-free formula (Lipari and Szabo, 1982a,b)

where *S*^{2} is the square of the Lipari and Szabo generalised order parameter and *τ*_{e} is the effective correlation time.
The order parameter reflects the amplitude of the motion and the correlation time in an indication of the time scale of that motion.
The theory was extended by Clore et al. (1990) by the modelling of two independent internal motions using the equation

where *S*^{2}_{f} and *τ*_{f} are the amplitude and timescale of the faster of the two motions whereas *S*^{2}_{s} and *τ*_{s} are those of the slower motion.
*S*^{2}_{f} and *S*^{2}_{s} are related by the formula
*S*^{2} = *S*^{2}_{f}⋅*S*^{2}_{s}.

If these forms of the model-free spectral density functions are unfamiliar, that is because these are the numerically stabilised forms presented in d'Auvergne and Gooley (2008b).
The original model-free spectral density functions presented in Lipari and Szabo (1982a) and Clore et al. (1990) are not the most numerically stable form of these equations.
An important problem encountered in optimisation is round-off error in which machine precision influences the result of mathematical operations.
The double reciprocal
*τ*^{-1} = *τ*_{m}^{-1} + *τ*_{e}^{-1} used in the equations are operations which are particularly susceptible to round-off error, especially when
*τ*_{e} *τ*_{m}.
By incorporating these reciprocals into the model-free spectral density functions and then simplifying the equations this source of round-off error can be eliminated, giving relax an edge over other model-free optimisation software.

The relax user manual (PDF), created 2016-10-28.