Freezing transition of the random bond RNA model: statistical properties of the pairing weights

Abstract

To characterize the pairing-specificity of RNA secondary structures as a function of temperature, we analyse the statistics of the pairing weights as follows : for each base $(i)$ of the sequence of length N, we consider the $(N-1)$ pairing weights $w_i(j)$ with the other bases $(j \neq i)$ of the sequence. We numerically compute the probability distributions $P_1(w)$ of the maximal weight, the probability distribution $Π(Y_2)$ of the parameter $Y_2(i)= \sum_j w_i^2(j)$, as well as the average values of the moments $Y_k(i)= \sum_j w_i^k(j)$. We find that there are two important temperatures $T_c<T_{gap}$. For $T>T_{gap}$, the distribution $P_1(w)$ vanishes at some value $w_0(T)<1$, and accordingly the moments $\bar{Y_k(i)}$ decay exponentially in $k$. For $T<T_{gap}$, the distributions $P_1(w)$ and $Π(Y_2)$ present the characteristic Derrida-Flyvbjerg singularities at $w,Y_2=1/n$ for $n=1,2..$. In particular, there exists a temperature-dependent exponent $μ(T)$ that governs these singularities and the decay of the moments $ \bar{Y_k(i)} \sim 1/k^{μ(T)}$. The exponent $μ(T)$ grows from $μ(T=0)=0$ up to $μ(T_{gap}) \sim 2$. The study of spatial properties indicates that the critical temperature $T_c$ where the roughness exponent changes from the low temperature value $ζ\sim 0.67$ to the high temperature value $ζ\sim 0.5$ corresponds to the exponent $μ(T_c)=1$. For $T<T_c$, there exists frozen pairs of all sizes, whereas for $T_c< T <T_{gap}$, there exists frozen pairs, but only up to some characteristic length diverging as $ξ(T) \sim 1/(T_c-T)^ν$ with $ν\simeq 2$. The similarities and differences with the weight statistics in Lévy sums and in Derrida&#39;s Random Energy Model are discussed.