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, 1982b,a)
where S2 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 S2f and τf are the amplitude and timescale of the faster of the two motions whereas S2s and τs are those of the slower motion. S2f and S2s are related by the formula S2 = S2f⋅S2s.
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 2019-03-08.