E. Bodon,A. Del Popolo,L. Luksan,E. Spedicato
Abstract:Numerical experiments are performed in order to study the performance of ABS codes in solving non-linear systems of equations.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: **Whether random RNA sequences will experience glass transition at a finite temperature**.
Specifically, the author studied the secondary structures formed by random RNA sequences and explored the behavior of these structures at different temperatures. The main problem can be decomposed into the following points:
1. **Stability of the high - temperature phase**:
- First, the author proved the stability of the high - temperature phase (i.e., the molten phase) under weak - disorder conditions through two - replica calculation. They showed that at high temperatures, the behavior of random RNA is similar to that of the uniform - energy model, and the influence of disorder on the system is negligible.
2. **Low - temperature behavior and the existence of glass transition**:
- Next, the author considered the energy contributions of rare complementary subsequence regions, indicating that the high - - temperature phase cannot last until zero temperature. By analyzing these rare regions, they proved that below a certain critical temperature \( T_c \), the system will enter a low - temperature phase different from the molten phase, thus supporting the existence of a finite - temperature glass transition.
### Summary of specific content
#### Introduction
RNA is an important biopolymer, and its secondary structure is formed by base pairing. Unlike DNA, RNA is usually single - stranded and can fold on itself to form complex secondary structures. Understanding how these structures are encoded from the primary sequence is an important problem in theoretical biophysics. This paper studied the secondary structures formed by random RNA sequences and explored the behavior of these structures at different temperatures.
#### Model description
- **Definition of secondary structure**: The secondary structure of RNA is formed by a series of base pairings, and each base can participate in at most one pairing.
- **Energy model**: Each pairing \((i, j)\) has an interaction energy \( \epsilon_{i,j} \), and the total energy \( E[S] \) is the sum of all pairing energies.
- **Simplified model**: To study statistical properties, the author used a simplified energy model, where \( \epsilon_{i,j} \) are independent Gaussian random variables with a mean of \( \epsilon \) and a variance of \( D \).
#### High - temperature phase
- At high temperatures, the influence of disorder in the system is small, and it can be described by a uniform - energy model. Through the z - transformation, the author obtained the partition function \( Z_0(N) \) of the molten phase and showed the free - energy behavior of the system at high temperatures.
#### Numerical results
- The author verified through numerical simulation that at high temperatures, the behavior of random RNA is consistent with that of the molten phase; while at low temperatures, the free - energy of the system shows different behavior, implying the existence of a low - temperature phase.
#### Theoretical analysis
- **Two - replica calculation**: Through two - replica calculation, the author proved that under weak - disorder conditions, the high - temperature phase is stable.
- **Inconsistency of the low - temperature phase**: By considering rare complementary subsequence regions, the author proved that below a certain critical temperature \( T_c \), the molten - phase assumption no longer holds, thus supporting the existence of a finite - - temperature glass transition.
#### Conclusion
Through theoretical analysis and numerical simulation, the author proved that random RNA sequences will undergo glass transition at a finite temperature. This finding helps to understand the behavior of RNA molecules at different temperatures and provides a new perspective for future theoretical research.
### Key formulas
1. Recurrence relation of the partition function:
\[
Z_{i,j} = Z_{i,j - 1}+\sum_{k = i}^{j - 1}Z_{i,k - 1}\cdot e^{-\beta\epsilon_{k,j}}\cdot Z_{k + 1,j - 1}
\]
2. Partition function of the molten phase:
\[
Z_0(N)\sim N^{-\theta_0}\exp[-\beta N f_0(q)]
\]
where \( q = e^{-\beta\epsilon_0} \), \( \theta_0=\frac{3}{2} \), \( f_0(q)=-k_B T \l