Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: rna is an important biopolymer critical to all living systems  @xcite and may be the crucial entity in prebiotic evolution  @xcite . like for dna , there are four types of nucleotides ( or bases ) a , c , g , and u which , when polymerized can form double helical structures consisting of stacks of stable watson - crick pairs ( a with u or g with c ) . however unlike a long polymer of dna , which is often accompanied by a complementary strand and forms otherwise featureless double helical structures , a polymer of rna usually `` operates '' in the single - strand mode . it bends onto itself and forms elaborate spatial structures in order for bases located on different parts of the backbone to pair with each other , similar conceptually to how the sequence of an amino acid encodes the structure of a protein . understanding the encoding of structure from the primary sequence has been an outstanding problem of theoretical biophysics . most theoretical work in the past decade have been focused on the problem of protein folding , which is very difficult analytically and numerically  @xcite . here , we study the problem of rna folding , specifically the formation of rna _ secondary structures_. for rna , the restriction to secondary structures is meaningful due to a separation of energy scales . it is this restriction that makes the rna folding problem amenable to detailed analytical and numerical studies  @xcite . there exist efficient algorithms to compute the exact partition function of rna secondary structures  @xcite . together with the availability of carefully measured free energy parameters  @xcite describing the formation of various microscopic structures ( e.g. , stacks , loops , hairpins , etc . ) , the probable secondary structures formed by any given rna molecule of up to a few thousand bases can be obtained readily . on the experimental side , rna molecules of @xmath1  @xmath2 bases in length are available . furthermore , the restriction to secondary structures can be physically enforced in a salt solution with monovalent ions , e.g. , @xmath3 , so that controlled experiments are in principle possible  @xcite .    in this work , we are not concerned with the structure formed by a specific sequence . instead , we will study the statistics of secondary structures formed by the ensemble of _ long _ _ random _ rna sequences ( of at least a few thousand bases in length in practice ) . such knowledge may be of value in detecting important structural components in messenger rnas which may otherwise be regarded as random from the structural perspective , in understanding how functional rnas arise from random rna sequences  @xcite , or in characterizing the response of a long single - stranded dna molecule to external pulling forces  @xcite . more significantly from the theoretical point of view , the rna secondary structure problem presents a rare tractable model of a random heteropolymer where concrete progress can be made regarding the thermodynamic properties  @xcite . nevertheless , there are many gaps in our understanding . this paper is a detailed report of our on - going effort in this regard . it provides a self - contained introduction of the random rna problem to statistical physicists as a novel problem of disordered systems , and depicts several approaches we have tried to characterize this system . the manuscript is organized as follows : in sec . [ sec_secondarystructure ] , we provide a detailed introduction to the phenomenology of rna secondary structure formation . we review the key simplifications which form the basis of efficient computing as well as exact solutions in some cases . we also review the properties of the  molten phase " , which is the simplest possible phase of the system assuming sequence disorder is not relevant . in sec . [ sec_randomnesshight ] , we consider the effect of sequence disorder at high temperatures . we show numerical evidence that the random rna sequence is in the molten phase at sufficiently high temperatures , and support this conclusion by solving the two - replica system which can be regarded as a perturbative study on the stability of the molten phase . in sec . [ sec_randomnesslowt ] , we provide a scaling argument , and show why the molten phase should break down at low enough temperatures . this is followed by a detailed numerical study of the low temperature regime . we apply the droplet picture and characterize the statistics of large - scale , low - energy excitations of the secondary structures from the ground state structure . our results support the existence of a very weak ( i.e. , marginal ) glass phase characterized by logarithmic excitation energies . finally , we describe the intermediate temperature regime where the system makes the transition from the glass phase to the molten phase . the solution of the two - replica problem is relegated to the appendices . we present two approaches : in appendix  [ app_intuitive ] , we provide a mapping of the two - replica problem to the denaturation of an effective single rna in @xmath0-dimensional embedding space ; this approach highlights the connection of the rna problem to the self - consistent hartree theory and should be most natural to field theorists . in appendices  [ app_solution ] and  [ app_derivez1 ] , we present the exact solution . it is hoped that the two - replica solution may be helpful in providing the intuition needed to tackle the full @xmath4-replica problem . the secondary structure of an rna describes the configuration of base pairings formed by the polymer . if the pairing of the @xmath5 and @xmath6 bases in a polymer of @xmath7 total bases is denoted by @xmath8 with @xmath9 , then each secondary structure @xmath10 is defined by a list of such pairings , with each position appearing at most once in the list , and with the pairs subject to a certain restriction to be described shortly below . each such structure can be represented by a diagram as shown in fig . [ fig_diagrams ] , where the solid line symbolizes the backbone of the molecule and the dashed lines stand for base pairings . the structure shown can be divided into _ stems _ of consecutive base pairs and _ loops _ which connect or terminate these stems . in naturally occurring rna molecules , the stems typically comprise on the order of five base pairs . they locally form the same double helical structure as dna molecules . however , while the latter typically occur in complementary pairs and bind to each other , rna molecules are mostly single - stranded and hence must fold back onto themselves in order to gain some base pairings .            by a secondary structure , one often considers only the restricted set of base pairings where any two base pairs @xmath8 and @xmath11 in a given secondary structure are either independent , i.e. , @xmath12 , or nested , i.e. , @xmath13 . this excludes the so - called pseudo - knots ( as exemplified by fig .  [ fig_pseudoknot ] ) and makes analytical and numerical studies much more tractable . for an rna molecule , the exclusion of pseudo - knots is a reasonable approximation because the long pseudo - knots are kinetically difficult to form , and even the short ones occur infrequently in natural rna structures  @xcite . the latter is due to their relatively low binding energies for short sequences and the strong electrostatic repulsion of the backbone  because the polymer backbone is highly charged and pseudo - knotted configurations increase the density of the molecule , their formation can be relatively disfavored in low salt solution . similarly , the tertiary structures which involve additional interactions of paired bases are strongly dependent on electrostatic screening and can be  turned off " experimentally by using monovalent salt solution  @xcite . indeed , the pseudo - knots are often deemed part of the tertiary structure of an rna molecule . throughout this study , we will exclude pseudo - knots in our definition of secondary structures . without the pseudo - knots , a secondary structure can alternatively be represented by a diagram of non - crossing arches or by a  mountain " diagram as shown in fig . [ fig_strreps ] .    : in ( a ) the solid line symbolizes the stretched out backbone of the molecule while the dashed arches stand for the base pairs formed . due to the no pseudo - knot constraint two arches never cross . ( b ) shows an equivalent representation as a `` mountain diagram '' . it is a line derived from the arch diagram by going along the backbone from left to right and going one step up for every beginning arch , horizontally for each unbound base , and one step down for each ending arch . such a mountain always stays above the baseline and comes back to the baseline at base @xmath7 . ]      in order to calculate boltzmann factors within an ensemble of secondary structures , we need to assign an energy @xmath14 $ ] to each structure @xmath10 . each secondary structure can be decomposed into elementary pieces such as the stems of base pairs and the connecting loop regions as shown in fig . [ fig_diagrams ] . a common approach is to assume that the contributions from these structural elements to the total energy are independent of each other and additive .    within a stem of base pairs , the largest energy contribution is the _ stacking energy _ between two adjacent base pairs ( g - c , a - u , or g - u ) , and the total energy of the stem is the sum of stacking energies over all adjacent base pairs . since each secondary structure is defined as a _ single _ state in our ensemble , it is necessary to integrate out all other microscopic degrees of freedom of the bases within a given secondary structure and use an effective energy parameter for each base stacking . the most convenient one to use is the gibbs free energy of stacking  @xcite , which contains an enthalpic term due to base stacking , and an entropic term due to the loss of single - stranded degrees of freedom ( as well as the additional conformational change of the backbone and even the surrounding water molecules ) due to base pairing . the magnitudes of these stacking free energies actually depend on the identity of all four bases forming the two base pairs bracketing the stack and are dependent on temperature themselves . while their typical values are on the order of @xmath15 at room temperature , the enthalpic and entropic contributions are each on the order of @xmath16 . thus , upon moderately increasing the temperature from room temperature to about @xmath17 , the stacking free energies become repulsive and the rna molecule denatures . the stacking free energies account for most but not all of the entropic terms for a given secondary structure . there is an additional ( logarithmic ) `` loop energy '' term associated with the entropy loss of each _ closed loop _ of single - stranded rna formed by the secondary structure , as well as the energy necessary to bend the single strand . all of these energy parameters have been measured in great detail  @xcite . when incorporated into an efficient dynamic programming algorithm ( to be described below ) , they can rather successfully _ predict _ the secondary structures of many rna molecules of up to several hundred bases in length  @xcite .    in this paper , we investigate the statistical properties of long , random rna sequences far below the denaturation temperature . we are interested in generic issues such as the existence of a glass phase and various scaling properties . guided by experiences with other disordered systems  @xcite , we believe these generic properties of the system should not depend on the specific choice of the model details . since the full model used in refs  @xcite makes analytical and numerical studies unnecessarily clumsy , we will examine a number of simplified models , while preserving the most essential feature of the system , namely , the pattern of matches and mismatches between different positions of the sequence . as in the realistic model described above , we choose our reference energy to be the unbound state , so that each unbound base in a secondary structure is assigned the energy @xmath18 . we will neglect the logarithmic loop energy , which is important very close to the denaturation transition  @xcite where the average binding energy is close to zero , but not far below the denaturation temperature where most bases are paired . moreover , we will radically simplify the energy rules for base pairing : we neglect the stacking energies and instead associate an interaction energy @xmath19 with every pairing @xmath8 . thus , @xmath20=\sum_{(i , j)\in s}\varepsilon_{i , j}\ ] ] is the total energy of the structure @xmath10 .    within this model , it remains to be decided how to choose the energy parameters @xmath19 s . one possibility is to choose each of the bases @xmath21 randomly from the ` alphabet ' set @xmath22 and then assign @xmath23 with @xmath24 being the match or mismatch energy respectively . here , the value of @xmath25 is actually not essential as long as it is repulsive , since the two bases always have the energetically preferred option to not bind at all . thus the energetics of the system is set by @xmath26 . in our numerical study to be reported in secs . [ sec_randomnesshight ] and  [ sec_randomnesslowt ] , we will primarily use this model ( as long as it is larger than @xmath27 as explained below ) . indeed the alphabet size and the choice of the matching rule can be used as tuning parameters to change the strength of sequence disorder . but in our study , we choose to minimize the number of parameters and tune the effective strength of disorder by changing temperature . ] with @xmath28 . we will refer to this as the `` sequence disorder '' model .    for analytical calculations , it is preferable to treat all the @xmath19 s as _ independent _ identically distributed random variables , i.e. , to assume @xmath29=\prod_{1\le i < j\le n}\rho(\varepsilon_{i , j})\ ] ] for the joint distribution function @xmath30 $ ] of all the @xmath19 s . this choice neglects the correlations between @xmath19 and @xmath31 which are generated through the shared base @xmath32 ; it is an additional approximation on the model ( [ eq_energychoice ] ) . however , we do not anticipate universal quantities to depend on such subtle correlation of the @xmath19 s . this will be tested numerically by comparing the behavior of the model ( [ eq_energychoice ] ) with that of the model defined by eq . ( [ eq_uncorrelated ] ) together with @xmath33 this distribution is chosen to mimic the random sequence model ( [ eq_energychoice ] ) with a 4-letter alphabet , but it does not contain any correlation between the different @xmath19 s . we will refer to this model as defined by eqs . ( [ eq_uncorrelated ] ) and ( [ eq_uncorrelateddisorder ] ) as the `` energy disorder '' model .    in the actual analytical calculations , we will go even one step further and take the @xmath19 to be _ gaussian _ random variables specified by @xmath34 where @xmath35 is the average binding energy and @xmath36 is the variance . in this model ( referred to below as  gaussian disorder " model , ) the parameter @xmath36 provides us with a convenient measure of the disorder strength . again , universal quantities should not depend on the choice of the distribution functions . we will test this directly by performing numerical studies for these gaussian random energies , with @xmath37 chosen to match the first two moments of the distribution eq . ( [ eq_uncorrelateddisorder ] ) . in contrast to prior numerical studies  @xcite , we do not exclude base pairing between neighboring bases @xmath38 , i.e. , we do not set a minimal allowed length for the hairpins . since the results are qualitatively identical to the results of the simpler model presented here , we do not show this data . ] . setting a constraint on the minimal hairpin length would make the analytical study much more cumbersome . however , in the study by pagnani _ et al . _  @xcite , it has been argued that the system will not be frustrated ( and hence will not form a glass ) without this additional constraint . we believe this is an artifact of the 2-letter alphabet used by pagnani _ _ in order to generate the binding energy @xmath19 s via a rule similar to eq . ( [ eq_energychoice ] ) : it is simple to see that for any 2-letter sequence in which the like letters repel and unlike letters attract , one can always find the minimal total binding energy by pairing up neighboring bases of opposite types and removing them from the sequence if no additional constraints such as the minimal hairpin length are enforced . as we will discuss in detail in sec . [ sec_proofglass ] this is not a problem if the alphabet size is larger than 2 . thus , in our study , we use the sequence disorder model with a 4-letter alphabet , or the energy disorder model , without enforcing the minimal hairpin length constraint . while the minimal hairpin length ( of 3 bases ) is known for real rna folding , it should not change the universal properties of long rna sequences . once the energy of each secondary structure is defined , we can study the partition function @xmath39}\ ] ] of the molecule where @xmath40 denotes the set of _ all _ allowed secondary structures of a polymer of @xmath7 bases , and @xmath41 . to calculate this partition function , it is useful to study the _ restricted _ partition function @xmath42 of the substrand from position @xmath43 to position @xmath44 of the rna molecule . given the model ( [ eq_energy ] ) , the restricted partition functions can be split up according to the possible pairings of position @xmath44 . this leads to the recursive equation  @xcite @xmath45 with @xmath46 being the total partition function of the molecule . in terms of the arch diagrams introduced in fig . [ fig_strreps](a ) this can be represented as @xmath47{receq.eps}}\ ] ] where the wavy lines stand for the restricted partition functions . this is easily recognized as a hartree equation . since the restricted partition functions on the right hand side of this equation all correspond to shorter pieces of the rna molecule than the left hand side , this equation allows one to calculate the exact partition function of an rna molecule of length @xmath7 with _ arbitrary _ interactions @xmath19 in @xmath48 time . this is accomplished by starting with the partition functions for single bases and recursively applying eq . ( [ eq_partfuncrec ] ) , and is known as a dynamic programming algorithm  @xcite . this algorithm allows one to compute numerically the partition function involving all secondary structures , for arbitrary rna molecules of up to @xmath49 bases . it also forms the basis of analytical approaches to the problem as we will see shortly . apart from the partition function itself , we will use additional observables in order to characterize the behavior of rna secondary structures . one such quantity of interest is the binding probability @xmath50 , i.e. , the probability that positions @xmath43 and @xmath44 are paired given the @xmath19 s , @xmath51 where @xmath52 is given by the recursion equation  ( [ eq_partfuncrec ] ) and @xmath53 is the partition function of the sequence @xmath54 . the latter can be calculated as the quantity @xmath55 when applying the recursion eq . ( [ eq_partfuncrec ] ) to the duplicated sequence @xmath56 . thus , all @xmath57 such constraint partition functions can be calculated with the same recursion in @xmath48 time . the logarithms @xmath58 of these binding probabilities have a natural interpretation : they can be read as the `` pinching free energies '' , i.e. , as the free energy cost of a pinch between positions @xmath43 and @xmath44 and the unperturbed state . we will make extensive use of this concept of pinched structures in our discussion of the low temperature behavior of rna secondary structures in sec . [ sec_randomnesslowt ] . in our numerical investigations , we will choose as a representative of all the pinching energies for different positions by @xmath59 which is the free energy cost of the largest possible pinch that splits the molecule of length @xmath7 into two pieces of length @xmath60 each . another quantity which describes a secondary structures is its `` size profile '' . as an intrinsic measure of the size of a given secondary structure @xmath10 , we use the `` ladder distance '' @xmath61 between the base at position @xmath62 and the base at position @xmath43 , which is the the number of pairings ( or ladders ) one has to cross to go from a pair involving base @xmath62 to the base @xmath43 ; see fig . [ fig_hdef ] . it can be defined for each secondary structure @xmath10 as the total number of pairings @xmath63 that bracket @xmath43 , i.e. , @xmath64 this quantity can be very easily visualized as the `` height '' at position @xmath43 of the mountain representation of the secondary structure @xmath10 as shown in fig . [ fig_strreps](b ) . a quantity characterizing the full ensemble of secondary structures is the _ thermal average _ @xmath65 of this size profile over all secondary structures with their respective boltzmann factors ; it can be straightforwardly calculated from the probabilities @xmath66 as @xmath67 since we expect all positions in the sequence to behave in a similar way , in our numerics we will summarize the properties of the size profile by the ladder distance from the first to the middle base , i.e. , we will study @xmath68 as a quantity representing the overall `` size '' of an ensemble of secondary structures .     of a secondary rna structure : the size profile measures the extension of the structure if drawn as a planar diagram . as an intrinsic definition of @xmath69 which captures this notion of the size of a secondary structure at position @xmath43 , we use the `` ladder distance '' of base @xmath43 from the end of the molecule , i.e. , the number of base pairs which have to be crossed when connecting position @xmath43 to position @xmath62 along the folded structure as indicated by the dashed line . ]        if sequence disorder does not play an important role , we may describe the rna molecule by replacing all the binding energies @xmath19 by some effective value @xmath70 . as we will see later , this will be an adequate description of our random rna models at high enough temperatures ( but before denaturation . ) for the real rnas , this provides a coarse grained description of repetitive , self - complementary sequences , e.g. , cagcag ... cag , which are involved in a number of diseases  @xcite . we will refer to rna which is well described by this model without sequence disorder as being in the `` molten '' phase . it serves as a starting point for modeling non - specific self - binding of rna molecules , and its properties will form the basis of our study of the random rna at low temperatures .      since in the absence of sequence disorder , the energy of a structure @xmath10 depends only on the number of paired bases @xmath71 of this structure , we can write the partition function in the molten phase as @xmath72\ ] ] the partition functions of the sub - strands @xmath42 become translationally invariant and can be written as @xmath73 where @xmath74 is only a function of the length @xmath7 . the recursion equation  ( [ eq_partfuncrec ] ) then takes the form @xmath75 where @xmath76 upon introducing the @xmath77-transform @xmath78 the convolution can be eliminated and the recursion equation turns into a quadratic equation @xmath79 with the solution @xmath80 performing the inverse @xmath77-transformation in the saddle point approximation yields the expression  @xcite @xmath81 in the limit of large @xmath7 , with the exponent @xmath82 and the non - universal quantities @xmath83 and @xmath84^{1/2}$ ] . this result characterizes the state of the rna where a large number of different secondary structures of equal energy coexist in the thermodynamic ensemble , and the partition function is completely dominated by the configurational entropy of these secondary structures . while the result is derived specifically for the special case @xmath85 , we will argue below that it is applicable also to random @xmath19 s at sufficiently high temperatures , in the sense that for long rna molecules , the partition function is dominated by an exponentially large number of secondary structures all having _ comparable _ energies ( within @xmath86 ) that are smoothly related to each other in configuration space . the latter is what we meant by the `` molten phase '' . the exponent @xmath82 is an example of a scaling exponent characteristic of the molten phase . this and other exponents can be derived in a geometric way by the `` mountain '' representation of secondary structures as illustrated in fig . [ fig_strreps](b ) . each such mountain corresponds to exactly one secondary structure . in the molten phase , the weight of a secondary structure @xmath10 is simply given by @xmath87 . this can be represented in the mountain picture by assigning a weight of @xmath88 to every upward and downward step and a weight of @xmath62 to every horizontal step . since the only constraints on these mountains are ( i ) staying above the baseline , and ( ii ) returning to the baseline at the end , the partition function of an rna of length @xmath7 is then simply that of a random walk of @xmath7 steps , constrained to start from and return to the origin , in the presence of a _ hard wall _ at the origin , with the above weights ( @xmath89 or @xmath62 ) assigned to each allowed step . this partition function is well - known to have the characteristic @xmath90 behavior which we formally derived in the last section  @xcite .    in this framework , it also becomes obvious why imposing a minimal hairpin length does not change the universal behavior of rna at least in this molten phase : if the minimal allowed size of a hairpin is @xmath91 , this enforces a potentially strong _ penalty _ for the formation of a hairpin , since with every hairpin @xmath91 bases are denied the possibility of gaining energy by base pairing . this tends to make branchings less favorable and thus leads to longer stems . however , this additional constraint translates in the mountain representation into the rule that an upwards step may not be followed by a downwards step within the next @xmath91 steps . this is clearly a _ local _ modification of the random walk . thus , it does not change universal quantities although the above mentioned suppression of branchings will require much longer sequences in order to observe the asymptotic universal behavior . for real rna parameters , the crossover length is very long due to this effect . for example , it is several hundred nucleotides for the cag repeat , and even longer for some other repeats . another characteristic exponent describes the scaling of the ladder size @xmath92 with the sequence length @xmath7 . as already mentioned in its definition ( [ eq_hdef ] ) , @xmath92 is equivalent to the average `` height '' of the mid - point of the sequence in the mountain picture . in the molten phase , the random walk analogy immediately yields the result @xmath93 where @xmath94 denotes ensemble average in the molten phase .    as should be clear from the coarse - grained view depicted in fig . [ fig_hdef ] , the ensemble of rna secondary structures in the molten phase can be mapped directly to the ensemble of _ branched polymers_. these branched polymers are _ rooted _ at the bases @xmath95 and @xmath96 of the rna . in this context , @xmath82 is known as the configuration exponent of the rooted branched polymer  @xcite . additionally from the result ( [ eq_hscaling ] ) , we see that the ladder length of the branched polymer scales , a well - known result for the branched polymer in the absence of self - avoidance  @xcite . ] as @xmath97 . because of the very visual analogy of the secondary structures to branched polymer , we refer to the configurational entropy of the secondary structures as the `` branching entropy '' .    finally , the binding probabilities @xmath50 defined in eq . ( [ eq_returnprobability ] ) only depend on the distance @xmath98 in the molten phase , i.e. , @xmath99 . the behavior of this function can be derived explicitly by inserting the result eq . ( [ eq_moltenz ] ) for the partition function into eq . ( [ eq_returnprobability ] ) . alternatively , one just needs to recognize that @xmath100 corresponds in the random walk analogy to the first - return probability of a random walk after @xmath101-steps . in either case , one finds the result @xmath102 i.e. , the return probability decays with the separation @xmath101 of the two bases as a power law with the configuration exponent @xmath82 . for the pinching free energy @xmath103 , we simply set @xmath104 and obtain @xmath105 for large @xmath7 , i.e. , it scales _ logarithmically _ in the molten phase . this logarithmic dependence merely reflects the loss in branching entropy due to the pinching constraint and is a manifestation of the configuration exponent @xmath82 . there are in principle three different scenaria for the behavior of long random rna sequences . ( i ) disorder is irrelevant at any finite temperature , so that the molten phase description presented in sec . [ sec_molten ] applies to long rnas at all temperatures . ( ii ) disorder is relevant at all temperatures , and the molten phase description would be completely inadequate . ( iii ) there is a finite temperature @xmath106 above which the molten description of random rna is correct , while below @xmath106 a qualitatively different description is needed . in accordance with the statistical physics literature , we will refer to the non - molten phase as the glass phase , and @xmath106 as the glass transition . the purpose of the study is to determine which of these three scenaria is actually realized , and to characterize the glass phase if either ( ii ) or ( iii ) is realized .    in this section , we study the high temperature behavior and demonstrate that the molten phase is stable with respect to weak sequence disorder . this ensures that the molten description of rna given in sec . [ sec_molten ] is at least valid at high enough temperatures , thereby ruling out scenario ( ii ) . we will address the question of whether there is a glass phase at low but finite temperatures in sec . [ sec_randomnesslowt ] . before we engage in detailed calculations , we want to convince ourselves with the help of some numerics that weak disorder does not destroy the molten phase . to this end , we study the observables introduced in sec . [ sec_observables ] . we generate a large number of disorder configurations , i.e. , interaction energies @xmath19 using the 3 models introduced in sec . [ sec_intenergies ] : sequence disorder , energy disorder , and gaussian disorder as described by eq . ( [ eq_energychoice ] ) ,  eqs . ( [ eq_uncorrelated ] ) and ( [ eq_uncorrelateddisorder ] ) , and eqs . ( [ eq_uncorrelated ] ) and ( [ eq_gaussiandisorder ] ) respectively , with @xmath107 . then , we calculate the observables @xmath92 and @xmath103 for each disorder configuration at the relatively large temperature of @xmath108 and average the obtained values over many disorder configurations . in order to keep the numerical effort manageable , we average over @xmath109 random sequences for @xmath110 , over @xmath111 sequences for @xmath112 , and over @xmath113 sequences for @xmath114 and @xmath115 . [ fig_moltennumerics ] shows the results ; disorder averaged quantities are denoted by an overline throughout the text . we see that the data for @xmath116 follows a power law with a fitted exponent @xmath117 , with the exponent value decreasing for larger @xmath7 s . this result is consistent with the prediction eq . ( [ eq_hscaling ] ) for the molten phase . also , the pinching free energy follows the predicted logarithmic behavior eq . ( [ eq_moltenpinch ] ) without any noticeable difference between the three choices of disorder . taken together , these results indicate that the three models of disorder belong to the same universality class , i.e. , the molten phase description of the uniformly attracting rna , at high temperatures .     of the sequence at @xmath108 . each plot shows data for three different choices of disorder according to eqs . ( [ eq_energychoice ] ) , ( [ eq_uncorrelateddisorder ] ) , and  ( [ eq_gaussiandisorder ] ) . ( a ) shows the scaling of the average size @xmath116 and the dashed line is the best fit @xmath117 to a power law . ( b ) shows the free energy of the largest pinch as defined in  ( [ eq_largestpinchfreeenergy ] ) . the dashed line is up to an additive constant the logarithmic behavior @xmath118 predicted in eq . ( [ eq_moltenpinch ] ) . the statistical fluctuations are smaller than the size of the symbols in both plots . all plots suggest that the behavior of rna secondary structures at high temperatures is well described by the molten phase picture and independent of the disorder . ] now , we will establish the stability of the molten phase against weak disorder by an analytical argument . we will use gaussian disorder characterized by eqs . ( [ eq_uncorrelated ] ) and  ( [ eq_gaussiandisorder ] ) . as we have shown above , the different microscopic models of binding energies all yield the same scaling behaviors . with the uncorrelated gaussian energies , it is possible to perform the ensemble average of the partition function @xmath119 of @xmath4 rna molecules sharing the same disorder . the disorder - averaged free energy can then in principle be obtained via the `` replica - trick '' @xmath120 , by solving the @xmath4-replica problem  @xcite . the @xmath4-replica partition function can be written down formally as    @xmath121}\\ & = & \sum_{\{s_1\}}\ldots\sum_{\{s_n\ } } \prod_{k=1}^n\exp\left[-\beta\avereps|s_k|\right ] \exp\left[\frac{1}{2}\beta^2\sum_{k=1}^n\sum_{l=1}^n\sum_{(i , j)\in s_k } \sum_{(r , s)\in s_l } \overline{(\varepsilon_{i , j}-\avereps)(\varepsilon_{r , s}-\avereps ) } \right]\\ & = & \sum_{\{s_1\}}\ldots\sum_{\{s_n\}}\prod_{k=1}^n \exp\left[-\beta\avereps|s_k|\right ] \exp\left[\frac{1}{2}\beta^2d\sum_{k=1}^n\sum_{l=1}^n|s_k\cap s_l|\right]\\ & = & \sum_{\{s_1\}}\ldots\sum_{\{s_n\}}\prod_{k=1}^n q^{|s_k| } \prod_{1\le k < l\le n}\widetilde{q}^{|s_k\cap s_l| } \ ] ]    where @xmath122 are the two relevant `` boltzmann factors '' . this effective partition function has a simple physical interpretation : it describes @xmath4 rna molecules subject to a _ homogeneous _ attraction with effective interaction energy @xmath123 between any two bases of the same molecule . as before , this effective attraction is characterized by the factor @xmath124 . in addition , there is an inter - replica attraction characterized by the factor @xmath125 for each bond _ shared _ between any pair of replicas . the inter - replica attraction is induced by the same sequence disorder shared by all replicas . for example , if the base composition in one piece of the strand matches particularly well with another piece , then there is a tendency to pair these pieces together in all replicas . thus , the inter - replica attraction can potentially force the different replicas to `` lock '' together , i.e. , to behave in an correlated way . indeed , the distribution of inter - replica correlations , usually measured in terms of `` overlaps '' , is a common device used to detect the existence of a glass phase in disordered systems  @xcite . the full @xmath4-replica problem is difficult to solve analytically . we will examine this problem in the regime of small @xmath36 , aiming to resolve the relevancy of disorder in a perturbative sense . since the lowest order term of the fully random problem in a perturbation expansion in @xmath36 corresponds to the two - replica ( @xmath126 ) problem , we will focus on the latter in order to study the small-@xmath36 behavior of the full problem . the solution of the two - replica problem will also illustrate explicitly the type of interaction one is dealing with , thereby providing some intuition needed to tackle the full problem . it turns out that the two - replica problem can be solved exactly . here , we outline the saline features of the solution . details of the calculation and analysis are provided in the appendices . we will find that the two - replica system has a phase transition between the molten phase in which the two replicas are uncorrelated and a nontrivial phase in which the two replicas are completely locked together in the thermodynamic limit . the transition occurs at a finite temperature @xmath127 which approaches zero as @xmath128 . thus , the effect of weak disorder is _ irrelevant _ at finite temperatures . let us denote the two - replica partition function @xmath129 for two strands each of length @xmath7 by @xmath130 , where we keep the dependence on @xmath124 implicit . then , @xmath131 the key observation which allows us to solve the two - replica problem is that for each given pair of secondary structures , the bonds shared by two replicas ( hereafter referred to as `` common bonds '' ) form a valid secondary structure by themselves ( see fig .  [ fig_combonds ] . ) thus , we can rearrange the summation over the pairs of secondary structures in the following way : we first sum over all possible secondary structures of the common bonds . for a _ given configuration _ of the common bonds , we then sum over the remaining possibilities of intra - replica base pairings for each replica , with the constraint that no new common bonds are created . note that the common bonds partition the diagram into a number of `` bubbles '' , shown as the shaded regions in fig . [ fig_combonds ] . due to the exclusion of pseudo - knots from the valid secondary structures , only bases belonging to the same bubble can be paired with each other . thus , the two - replica partition function can be written as @xmath132 where the factor @xmath133 is the weight of each common bond , and @xmath134 is the sum of all possible intra - replica pairings of of the @xmath5 bubble of @xmath135 bases in @xmath10 , with the restriction that there are no common bonds . it should be clear that @xmath136 neither depends on the number of stems branching out from the bubble @xmath43 nor the positions of these stems relative to the bases within the bubble . it depends on @xmath10 only through the number of bases @xmath135 in the bubble and is given by a single function @xmath137 independent of @xmath43 . this function can be written explicitly as @xmath138    with eqs . ( [ eq_gintermsofloops ] ) and ( [ eq_defq ] ) , the two - replica problem is reduced to an effective _ single _ homogeneous rna problem , with an effective boltzmann weight @xmath139 for each pairing , and an effective weight @xmath137 for each single stranded loop . as described in appendix  [ app_intuitive ] , this problem becomes formally analogous to that of an rna in the vicinity of the denaturation transition , with @xmath137 being the weight of a single polymer loop fluctuating in @xmath0-dimensional embedding space . the competition between pairing energy and the bubble entropy leads to a _ phase transition _ for the two - replica problem , analogous to the denaturation transition for a single rna . the details of this transition are given in appendix  [ app_solution ] , where the partition function ( [ eq_gintermsofloops ] ) is solved exactly . the exact solution exploits the relation @xmath140 which follows from the definitions ( [ eq_defg ] ) and ( [ eq_defq ] ) , and turns eq .  ( [ eq_gintermsofloops ] ) into a recursive equation for @xmath141 . the solution is of the form @xmath142 for large @xmath7 , with two different forms for @xmath143 and @xmath144 depending on whether @xmath125 is above or below the critical value @xmath145 here @xmath74 is the molten phase partition function , whose large @xmath7 asymptotics is given by eq . ( [ eq_moltenz ] ) and whose values for small @xmath7 can be calculated explicitly from the recursion eq . ( [ eq_g0recursion ] ) . thus , the actual value of @xmath146 can be found for any given @xmath124 .    for @xmath147 , we have @xmath148 and @xmath149 where @xmath150 according to eqs . ( [ eq_gatsingdef ] ) and  ( [ eq_moltenzresult ] ) . in this regime , the two - replica partition function @xmath141 is essentially a product of two single - replica partition functions @xmath151 . compared to eq . ( [ eq_moltenz ] ) , we can identify @xmath143 as @xmath152 , and @xmath144 as a modified version of @xmath153 . since there is no coupling of the two replicas beyond a trivial shift in the free energy per length , @xmath154 , we conclude that the disorder coupling is irrelevant . hence the two - replica system is in the molten phase in this regime .    for @xmath155 , we have @xmath156 and @xmath144 is given as the implicit solution of an equation involving only single - replica partition functions as shown in eqs . ( [ eq_zkdef ] ) and  ( [ eq_zetaresult ] ) . here , the partition function of the two - replica system is found to have the _ same _ form as that of the single - replica system in ( [ eq_moltenz ] ) . this result implies that the two replicas are _ locked _ together via the disorder coupling , and the molten phase is no longer applicable in this regime .    of course , as already explained above , only the weak - disorder limit ( i.e. , @xmath157 ) of the two - replica problem is of relevance to the full random rna problem . in this limit , @xmath158 while @xmath146 is found by evaluating eq . ( [ qc ] ) with @xmath159 . it can be easily verified that @xmath160 as long as @xmath124 is finite . thus in the weak disorder limit , we have @xmath147 , indicating that the molten phase is an appropriate description for the random rna . unfortunately , the two - replica calculation can not be used in itself to deduce whether the molten phase description breaks down at sufficiently strong disorder or low temperature . based on this analysis , we can not conclude whether the type of phase transition obtained for the two - replica problem is present in the full problem . having established the validity of the molten phase description of random rna molecules at weak disorder or high temperatures , we now turn our focus onto the low temperature regime . first , we will give an analytical argument for the existence of a glass phase at low temperatures . then , we will present extensive numerical studies confirming this result and characterizing this glass phase . we will start by showing that the molten phase can not persist for all temperatures down to @xmath161 . to this end , we will _ assume _ that long random rna is in the molten phase for _ all _ temperatures , i.e. , that the partition function for any substrand of large length @xmath162 is given by @xmath163\ ] ] with some effective temperature - dependent prefactor @xmath164 and free energy per length @xmath165 . then , we will show that this assumption leads to a contradiction below some temperature @xmath166 . this contradiction implies that the molten phase description breaks down at some finite @xmath167 . to be specific , we will consider the sequence disorder model ( [ eq_energychoice ] ) in this analysis . the quantity we will focus on is again the free energy @xmath168 of the largest possible pinch . under the assumption that the random sequences are described by the molten phase , it is given by @xmath169 for large @xmath7 and all @xmath170 _ independently _ of the values of the effective prefactor @xmath164 and the free energy per length @xmath165 ( see eq . ( [ eq_moltenpinch ] ) . )    on the other hand , we can study this pinching free energy for each given sequence of bases drawn from the ensemble of random sequences . for each such sequence , we can look for a continuous segment of @xmath171 watson - crick pairs @xmath32@xmath172@xmath173@xmath174 where the bases @xmath175 are within the first half of the molecule and the bases @xmath176 are in the second half ( see fig .  [ fig_goodmatch](a ) . ) for random sequences , the probability of finding such exceptional segments decreases exponentially with the length @xmath101 , with the largest @xmath101 in a sequence of length @xmath7 being typically @xmath177 for exact complementary matches , the proportionality constant is known to be @xmath178  @xcite .     and @xmath179 and the other of which is between positions @xmath180 and @xmath7 . such a piece of length @xmath181 can be found for almost all sequences . ( b ) shows how restricting configurations to those in which the good match forms watson - crick base pairs splits the molecule into two loops which can still form base pairs within the loops independently from each other . ] now , we calculate the pinching free energy @xmath182 by evaluating the two terms separately . the partition function for the unpinched sequence contains _ at least _ all the configurations in which the two complementary segments @xmath183 and @xmath176 are completely paired ( see fig .  [ fig_goodmatch](b ) ) . thus , @xmath184 where @xmath185 is the free energy of the ensemble of structures in which the two complementary segments are paired . the latter is the sum of the free energy of the paired segments and those of the two remaining substrands @xmath186 of length @xmath187 and @xmath188 ( wrapping around the end of the molecule ) of length @xmath189 , i.e. , @xmath190.\ ] ]    the free energy @xmath191 in the presence of the pinch is , by the assumption of the molten phase , the interaction energy of the pinched base pair @xmath192@xmath193 plus the molten free energy of the substrand @xmath194 and the molten free energy of the substrand @xmath195 , i.e. , according to eq . ( [ eq_moltenzeff ] ) @xmath196 up to terms independent of @xmath7 . combining this with eqs . ( [ eq_deltafasdifference ] ) , ( [ eq_unpinchedandpaired ] ) , and  ( [ eq_fpaired ] ) , we get @xmath197 + \ell [ u_{\mathrm{m}}\!+\!2f_0(t)].\ ] ] using the result ( [ eq_lambda ] ) and that @xmath198 and @xmath199 are typically proportional to @xmath7 , we finally obtain @xmath200 \lambda^{-1 } \ln n\ ] ] for large @xmath7 . this is only consistent with eq . ( [ eq_assumedmoltenpinch ] ) if @xmath201.\ ] ]    now , @xmath165 is a free energy and is hence a monotonically decreasing function of the temperature . thus the validity of the inequality ( [ eq_assumedmoltenpinch ] ) depends on the behavior of its right hand side at low temperatures . as @xmath202 , the inequality can only hold if @xmath203 . since the average total energy at @xmath204 is @xmath205 times the average total number of matched pairs of a random sequence , then @xmath206 is simply the fraction of matches and @xmath207 is the fraction of bases not matched ( for @xmath208 . ) clearly , @xmath207 can not be negative , and the inequality ( [ eq_assumedmoltenpinch ] ) must fail at some finite temperature unless @xmath209 . we can make a simple combinatorial argument to show that in most cases the fraction @xmath207 of unbound bases must strictly be positive . to illustrate this , let us generalize the `` alphabet size '' of the sequence disorder model of sec . [ sec_intenergies ] from @xmath210 to an arbitrary even integer @xmath211 . we will still adopt the energy rule ( [ eq_energychoice ] ) where each of the @xmath212 bases can form a `` watson - crick '' pair exclusively with one other base . let us estimate the number of possible sequences for which the fraction of unmatched bases @xmath207 is zero in the limit of long sequence length @xmath7 at @xmath204 . since at @xmath204 , only watson - crick ( w - c ) pairs can be formed , we only need to count the number of sequences for which the fraction of w - c paired bases is @xmath62 . this means that except for a sub - extensive number of bases , all have to be w - c paired to each other . from the mountain picture ( fig . [ fig_strreps ] ) , it is clear that the number of possible secondary structures for such sequences must scale like @xmath213 , since the fraction of horizontal steps is non - extensive so that at each step , there are only the possibilities for the mountain to go up or down . for each of the @xmath179 pairings in one of these @xmath213 structures , there are @xmath212 ways of choosing the bases to satisfy the pairing . so for each structure , there are @xmath214 ways of choosing the sequence that would guarantee the structure . since there are a total of @xmath215 sequences , it is clear that the fraction of sequences with all ( but a sub - extensive number of ) w - c pairs becomes negligible if @xmath216 thus , for @xmath217 , we must have @xmath218 .    for @xmath219 , the left hand side of eq . ( [ eq_frustration ] ) grows faster than its right hand side . this reflects the absence of frustration in this simple two - letter model as already discussed at the end of sec . [ sec_intenergies ] . one way to retain frustration is to introduce additional constraints , e.g. , the minimal hairpin length used in ref .  @xcite . with this constraint , a structure with a sub - extensive number of unmatched bases can only contain a sub - extensive number of hairpins . in the mountain picture , this means that except for a sub - extensive number of steps , there is only one choice to go up or down at every step . this changes eq . ( [ eq_frustration ] ) to @xmath220 . it ensures frustration since @xmath218 for _ all _ @xmath212 . since a minimal length of @xmath221 bases is necessary in the formation of a real hairpin , real rna is certainly frustrated by this argument . the random sequence model which we study in this paper is marginal since @xmath222 and there is no constraint on the minimal hairpin length . in this case , all the prefactors on the two sides of eq . ( [ eq_frustration ] ) ( e.g. , the overcounting of sequences that support more than one structure ) must be taken into account . we will not undertake this effort here , but will verify numerically in sec . [ sec_estimatets ] that @xmath218 also in this case .    in all cases with @xmath218 , it follows that there is some unique temperature @xmath223 below which the consistency condition  ( [ eq_consistencycondition ] ) breaks down , implying the inconsistency of the molten phase assumption in this regime . from this we conclude that there must be a phase transition away from the molten phase at some critical temperature @xmath224 . the precise value of the bound @xmath223 depends on @xmath225 which in turn depends on the stringency of the condition we impose on the rare matching segments . for instance , if we relax the condition of exact complementarity between two segments to allow for matches within each segment , then the constant @xmath225 will be reduced from @xmath226 and the value of @xmath223 will increase . this will be discussed more in sec . [ sec_estimatets ] .      the above argument does not provide any guidance on the properties of the low temperature phase itself . in order to characterize the statistics of secondary structures formed at low temperatures , we re - do the simulations reported in sec . [ sec_hightnumerics ] at @xmath227 in energy units set by @xmath107 . at this temperature , an unbound base pair is penalized with a factor @xmath228 relative to a watson - crick base pair , and a non watson - crick base pair is penalized even more . thus , only the minimal energy structures contribute ( for the sequence lengths under consideration here ) , and we may regard this effectively as at @xmath204 . as in sec . [ sec_hightnumerics ] , we average over @xmath109 realizations of the disorder for @xmath110 , over @xmath111 realizations for @xmath112 and over @xmath113 realizations for @xmath229 . [ fig_lowtnum ] shows the results for the ladder size @xmath230 of the structures for the three models of disorders . the ladder size still scales algebraically with the length of the sequences , with numerically determined exponents ranging from @xmath231 to @xmath232 for the different choices of disorder . they are clearly different from the square root behavior ( dotted line ) expected of the molten phase . thus this result reaffirms our expectation that the secondary structures of a random rna sequence at zero temperature indeed belongs to a phase that is different from the molten phase .     of secondary structures in the low temperature phase with the length @xmath7 of the sequences : the plot shows data for three different choices of disorder according to eqs . ( [ eq_energychoice ] ) , ( [ eq_uncorrelateddisorder ] ) , and  ( [ eq_gaussiandisorder ] ) at @xmath233 . the average system size follows a power law . however , the best fit of the data for sequence disorder at @xmath234 to a power law indicated by the dashed line leads to an exponent of @xmath232 . the corresponding fits for energy and gaussian disorder yield exponents of @xmath232 and @xmath231 , respectively . this is distinctively different from the square root behavior of the molten phase indicated by the dotted line . the comparison of this plot with its counterpart in fig . [ fig_moltennumerics](a ) suggests that the behavior of rna secondary structures at low temperatures is different from the molten phase . ] a key question in characterizing the thermodynamic properties of disordered systems is whether the zero - temperature behavior persists for a range of finite temperatures . if it does , then the system is said to have a finite - temperature glass phase . one way to address this question is to study the overlap between different replicas of the rna molecule as mentioned earlier . if a non - trivial distribution of these overlaps with significant weight on large overlaps persists into finite temperatures , then the finite - temperature glass phase exists . this approach was taken by previous numerical studies  @xcite . unfortunately , the results are inconclusive and even contradictory due to the weakness of the proposed phase transition  only the fourth temperature derivative of the free energy seems to show an appreciable singularity . moreover , due to limitations in the sequence lengths probed , it was difficult to get a good estimate of the asymptotic behavior of the overlap distribution .    in our study , we adopt a different approach based on the droplet theory of fisher and huse  @xcite . in this approach , one studies the `` large scale low energy excitations '' about the ground state . this is usually accomplished by imposing a deformation over a length scale @xmath235 and monitoring the minimal ( free ) energy cost of the deformation . this cost is expected to scale as @xmath236 for large @xmath101 . a positive exponent @xmath237 indicates that the deformation cost _ grows _ with the size . if this is the case , the thermodynamics is dominated by a few low ( free ) energy configurations in the thermodynamic limit , and the statistics of the zero - temperature behavior persists into finite temperatures . on the other hand , if the exponent @xmath237 is negative , then there are a large number of configurations which have low overlap with the ground state but whose energies are similar to the ground state energy in the thermodynamic limit . at any finite temperature @xmath170 , a finite fraction of these configurations ( i.e. , those within @xmath86 of the ground state energy ) will contribute to the thermodynamics of the system . the zero temperature behavior is clearly not stable to thermal fluctuations in this case , and no thermodynamic glass phase can exist at any finite temperature . the analysis of the previous section indicate the existence of a glass phase ; thus we expect to find the excitation energy to _ increase _ with the deformation size . it should be noted that this criterion for glassiness is purely thermodynamical in nature and does not make any statement about _ kinetics_. a system which is not glassy thermodynamically can still exhibit very large _ barriers _ between the many practically degenerate low energy configurations , leading to a _ kinetic glass_. a study of the kinetics of rna , e.g. , in terms of barrier heights , is naturally dependent on the choice of allowed dynamical pathways to transform one rna secondary structure into another one  @xcite . since the latter is a highly non - trivial problem , we will restrict ourselves to thermodynamics and use the droplet picture explained above as our criterion for the existence of a glass phase .      according to the criterion for glassiness just presented , our goal is to determine the value of the exponent @xmath237 for random rnas numerically . to this end , the choice of `` large scale low energy excitations '' needs some careful thoughts . as in every disordered system , there is a very large number of structures which differ from the minimal energy structure only by a few base pairs and which have an energy only slightly higher than the minimum energy structure . these structures are clearly not of interest here . instead we need to find a controlled way of generating droplet excitations of various sizes . we propose to use the pinching method introduced in sec . [ sec_observables ] as a way to generate the deformation , and regard the difference between the minimal energy pinched structure and the ground state structure as the droplet excitation . there are several desirable features about these pinch - induced deformations : first , it gives a convenient way of controlling the size of the deformation . if @xmath238 is a base pair that is bound anyways in the ground state , pinching this base pair does not have any effect and @xmath239 . if we pinch base @xmath43 with some other base @xmath240 , then we force at least a partial deformation of the ground state , for bases in the vicinity of @xmath43 , @xmath241 , and @xmath44 . this is illustrated in fig . [ fig_pinchingstep ] with the deformed region indicated by the shade . as we move the pinch further away from the ground state pairing , we systematically probe the effect of larger and larger deformations ( provided that a pinch only induces _ local _ deformation as we will show ) . second , the minimal energy or the free energy of the secondary structures subject to the pinch constraint is easily calculable numerically by the dynamic programming algorithm as shown in eqs . ( [ eq_returnprobability ] ) and  ( [ eq_pinchfreeenergy ] ) . third , the pinching of the bases in a sense mimics the actual dynamics of the rna molecule at low temperatures . in order for the molecule to transform from one secondary structure to another at a temperature where all matching bases should be paired , the bases have to make local rearrangements of the secondary structures much like the way depicted in fig .  [ fig_pinchingstep ]  @xcite . thus , the pinching energy provides the scale of variation in the local energy landscape for such rearrangements to the pinch free energies . ] . finally , `` pinching '' of a real rna molecule can be realized in the pulling of a long molecule through a pore  @xcite . and @xmath241 ( open circles ) form a base pair in the indicated minimal energy structure ( a ) . thus , forcing these two bases to be bound by pinching does not affect the structure at all . pinching of @xmath8 on the other hand will lead to a local rearrangement ( shaded region ) of the structure as shown in ( b ) . the effect of such a pinch depends on the number of base pairs of the minimal energy structure the pinch is incompatible with . as indicated by the arrow in ( a ) , this number is given by the ladder distance @xmath242 between @xmath43 and @xmath44 ; in this example @xmath243 . ]    a key question to the utility of these pinch deformations is whether the deformation is confined to the local region of the pinch as depicted in fig . [ fig_pinchingstep ] or whether it involves a global rearrangement of the structure . to test this aspect , we numerically study the changes in pinch free energy as a function of the `` size '' of a pinch . here , the definition of the pinch size needs some thought . consider a specific sequence whose minimal energy structure is @xmath244 . if the binding partner of base @xmath43 is base @xmath241 in the minimal energy structure , a natural measure for the size of a pinch @xmath245 with @xmath246 would be the ladder distance @xmath242 between base @xmath43 and base @xmath44 ; see fig . [ fig_pinchingstep ] . from the mountain representation ( fig . [ fig_strreps](b ) ) , it is easy to see that this is just the difference of the respective ladder distances of base @xmath44 and base @xmath43 from base @xmath62 as defined in eq . ( [ eq_defmountainheight ] ) , i.e. , @xmath247 . to find how the excitation energy depend on the pinch size , we just need to follow how the pinching free energy @xmath248 s depend statistically on the size @xmath242 s . to do so , we choose a large number of random sequences , and determine the minimal energy structure @xmath244 for each of these sequences . then , we compute the pinch free energies @xmath248 and the pinch size @xmath242 for _ all _ possible pinches @xmath8 for each sequence . afterwards , we average over all @xmath248 s with the same pinch size @xmath242 over all of the generated random sequences to obtain the function @xmath249 the results obtained at @xmath227 for a large range of sequence sizes from @xmath250 to @xmath115 are shown in fig . [ fig_fofh](a ) . we see that the data for different @xmath7 s fall on top of each other for small @xmath251 s , with @xmath252 this behavior explicitly shows that the pinch deformation is a _ local _ deformation . in particular , we see that for small @xmath251 s , the free energy cost is _ independent _ of the overall length @xmath7 of the molecule .     of minimal energy structure base pairs that are incompatible with the pinch for random sequence disorder at @xmath233 : ( a ) shows the raw data . for small @xmath253 the pinching free is independent of the length @xmath7 of the molecule and obeys a power law @xmath254 ( dashed line . ) this is consistent with the expectation that pinching at small @xmath251 leads to _ local _ rearrangements of the secondary structure . the apparent non - monotonic behavior at large @xmath251 is due to the small number of sequences in which such a value of @xmath251 is realized . ( b ) shows the same data , but the scaling of @xmath253 with @xmath7 is chosen in accordance with fig . [ fig_lowtnum ] while the scaling of the pinching free energy is chosen in accordance with the power law dependence estimated above . ] it is interesting to see at which @xmath251 the entire sequence is involved . one expects @xmath255 since @xmath230 gives the typical scale of the maximum ladder length . to test this , we normalized @xmath251 by @xmath256 and @xmath257 by @xmath258 ( such that the relation ( [ eq_fh ] ) is preserved for small @xmath251 s ) . the result is shown in fig . [ fig_fofh](b ) . we see that the data is approximately collapsed onto a single curve , indicating that pinching is indeed a good way of imposing a controlled deformation from the ground state .      the scaling plot of fig . [ fig_fofh](b ) indicates strongly that the energy associated with the pinch deformation _ increases _ with the size of the deformation , i.e. , @xmath259 . however , the effective exponent involved is small , making the result very susceptible to finite size effects . in order to decide on the glassiness of the system , we want to focus on the energy scales associated with the largest pinch deformations from the ground state . assuming that there is only a single energy scale associated with large pinches , we again study the free energy @xmath103 of the largest pinch as defined in eq . ( [ eq_largestpinchfreeenergy ] ) and average this over the ensemble of sequences . this quantity yields an _ upper bound _ estimate of the energy associated with large scale pinches for each sequence length @xmath7 . we find @xmath260 and @xmath261 to have the same scaling behavior , and thus present only data for the latter . ] . the results are shown in fig . [ fig_pinch](a ) for the three models of disorders . although a weak power law dependence of @xmath261 on @xmath7 can not be excluded , the fitted exponents obtained for the three models are different from each other , ranging from @xmath262 to @xmath263 . this is a strong sign of concern , since the exponents are expected to be independent of details of the models . in fig . [ fig_pinch](b ) , the same data is plotted on a log - linear scale . the data fall reasonably on a straight line for each of the models ( especially for large @xmath7 s ) , suggesting that the pinching free energy may actually increase _ logarithmically _ with the sequence length , similar to what is expected of the behavior in the molten phase ! however in this case , the prefactor of the logarithm depends on the choice of the model and is much larger than the factor @xmath264 expected of the molten phase ; see eq . ( [ eq_moltenpinch ] ) . for example , for the numerical data obtained at @xmath227 , the prefactor is approximately @xmath265 for the sequence disorder model , while the expected slope for the molten phase is @xmath266 at this temperature . having different logarithmic prefactors for the different models is not a concern , since a prefactor is a non universal quantity . thus , our numerical results favor a logarithmically increasing pinch energy , with a prefactor much exceeding @xmath267 at low temperature .    . the exponents are @xmath268 , @xmath269 , and @xmath258 for the sequence , energy , and gaussian disorder respectively . ( a ) shows the same data in a single logarithmic plot together with the best fits to a logarithmic dependence on @xmath7 . the statistical error of the data points is about the size of the symbols for large @xmath7 and smaller than that for @xmath270 . due to the apparent systematic bending of the data in the double logarithmic plot ( b ) we conclude that a logarithmic dependence fits the data better although we can not exclude a power law behavior with a very small power . ] what does this tell us about the possible glass phase of the random rna ? in order to answer this question , we should remind ourselves that rather special deformations are chosen in this study . for our choice of pinch deformations , we observe a logarithmic dependence of the gap between the ground state energy and the energy of the excited configurations on the length of the sequence or deformation . this corresponds to the marginal case of the droplet theory where the exponent @xmath237 vanishes . since the pinching free energies are increasing with length , we can not exclude a glass phase in the case @xmath271 . we can say , though , that the increase of the excitation energy with length is at most a power law with a very small exponent and most probably even less than any power law . therefore , a possible glass phase of rna has to be _ very weak_. if it turns out that the excitation energy is indeed a logarithmic function of length , with a non - vanishing prefactor as @xmath202 as our numerics suggest , then the low - temperature phase would be categorized formally as a _ marginal _ glass phase , analogous to behaviors found in some well - studied model of statistical mechanics  @xcite . in any case , we should note that the actual difference in the excitation energy is only a factor of @xmath210 across two - and - a - half decades in length . thus the glassy effect will be weak for practical purposes . on the other hand , the weak dependence of the excitation energy on length may be the underlying cause of discrepancies in the literature  @xcite regarding the existence of the glass phase for the random rna . now that we have studied in great detail the behavior of random rna in the low and the high temperature phase , we describe its behavior at intermediate temperatures . to this end , we again study the pinch free energies @xmath261 defined in eq . ( [ eq_largestpinchfreeenergy ] ) , but this time over a large range of temperatures . we concentrate on the sequence disorder model eq . ( [ eq_energychoice ] ) with @xmath107 , and study sequences of lengths up to @xmath114 .    from secs . [ sec_scaling ] , [ sec_hightnumerics ] , and  [ sec_weakness ] , we know that the pinch free energy @xmath261 depends logarithmically on the sequence length @xmath7 at both low and high temperatures . indeed , this logarithmic behavior seems to hold for _ all _ temperatures studied . the data for each temperature can easily be fitted to the form @xmath272 the prefactor @xmath273 is found to depend on temperature in a non - monotonic way as shown in fig . [ fig_logfactors ] . the figure contains values of @xmath273 extracted by fitting the data for @xmath234 to the form eq . ( [ eq_logform ] ) . the uncertainty of this fit is on the order of the size of the symbols or smaller . for high temperatures , we find @xmath274 ( dashed line in fig .  [ fig_logfactors ] ) as expected for the molten phase . at low temperatures , it starts from a finite value of the order @xmath62 and decreases linearly with temperature , as @xmath275 ( dotted line in fig . [ fig_logfactors ] ) . if we identify the glass transition temperature @xmath276 as the intersection of the dashed and the dotted lines , we get @xmath277     of the logarithmic dependence of @xmath278 on @xmath7 for random rna sequences generated by the sequence disorder model : at high temperatures , the prefactor indicated by the circle is well - described by the dashed line @xmath279 expected of the molten phase . at low temperatures , it again has a linear temperature dependence and is empirically fitted by the dotted line , @xmath280 . the numerical uncertainty in @xmath273 is of the order of or smaller than the size of the symbols . ] it is interesting to compare this estimate with the lower bound @xmath223 for the glass transition temperature given in sec . [ sec_proofglass ] . according to the consistency condition  ( [ eq_consistencycondition ] ) , this lower bound is defined by @xmath281=\frac{3}{2}k_b t^*\ ] ] with @xmath282 . it is necessary to determine the temperature dependence of the quantity @xmath283 on the left hand side of this equation numerically . to do this , we measure the total free energy of each sequence generated . averaging these free energies over all the sequences of a given length @xmath7 and temperature @xmath170 , and dividing the results by the respective lengths @xmath7 , we obtain an estimate @xmath284 of the free energy per length which approaches the desired @xmath165 for large @xmath7 .    : the data shown is taken at @xmath227 . for small @xmath7 it is dominated by the statistical fluctuations in the number of bases according to eq . ( [ eq_smallnf0 ] ) ( the dashed line ) . at large @xmath7 , it saturates to a positive constant . ] [ fig_ump2f0ofn ] shows how these estimates depend on the sequence length @xmath7 for the lowest temperature @xmath227 studied . instead of the free energy per length itself , the figure shows the fraction of unbound bases @xmath285 . for short sequences these estimates show a clear dependence on the sequence length @xmath7 this can be understood in terms of sequence - to - sequence fluctuations in the maximum number of possible pairings , due to fluctuations in the actual number of each type of bases present in a given sequence , even if all four bases are drawn with equal probability . this effect can be quantified by assuming that there is no frustration for small @xmath7 , i.e. , for any given sequence of the four bases @xmath286 , @xmath287 , @xmath151 , and @xmath288 , a secondary structure with the maximal number of watson - crick base pairs can be formed . if we denote by @xmath289 the number of times that the base @xmath290 appears in the sequence , the maximal number @xmath291 of pairings is given by @xmath292 . the fraction of unbound bases @xmath293 due to this effect can be computed straightforwardly approximating the multinomial distribution of @xmath294 by an appropriate gaussian distribution , with the result @xmath295 we expect this effect to be responsible for the increase in @xmath296 found in fig . [ fig_ump2f0ofn ] . indeed this effect , as indicated by the dashed line in the figure , explains the @xmath7 dependence of @xmath296 well for @xmath297 . however , we also see from the figure a clear saturation effect at large @xmath7 . this saturation reflects the finite fraction of unbound bases which is a frustration effect forced upon the system through the restriction on the type of allowed pairings in a allowed secondary structure . the unbound fraction @xmath298 is _ finite _ asymptotically as expected in sec . [ sec_proofglass ] . note that this value is _ remarkably small _ , as it implies that in the ground state structure of our toy random sequence , more than @xmath299 of the maximally possible base pairs are formed . but this is an artifact of the very simple energy rule used in our toy model . this fraction certainly will become smaller if the realistic energy rules are used , making the system more frustrated hence more glassy .    in order to obtain the temperature dependence of the quantity @xmath283 on the left hand side of eq . ( [ eq_defoftstar ] ) , we will use its value at @xmath114 as an estimate of its asymptotic value . the results are shown in fig . [ fig_ump2f0 ] . the behavior at low temperatures can be described by a linearly decreasing function , shown as the dotted line in fig . [ fig_ump2f0 ] . according to eq . ( [ eq_defoftstar ] ) , the temperature @xmath223 is obtained as the intersection of this curve and @xmath300 , shown as the dashed line in fig . [ fig_ump2f0 ] for @xmath178 . we find @xmath301 which is consistent with the estimate ( [ eq_tgestimate ] ) , but is a rather weak bound . improved bounds on @xmath106 can be made by relaxing the condition of perfect complementarity of the two segments imposed in sec . [ sec_proofglass ] . this leads to larger values of the prefactor @xmath302 in eq . ( [ eq_defoftstar ] ) , hence a smaller slope for the dashed line in fig . [ fig_ump2f0 ] , and a larger value of @xmath223 . while the details of improved bounds will be discussed elsewhere , let us remark here that from fig . [ fig_ump2f0 ] it is clear that no matter what the slope of the dashed line becomes , we will never have @xmath223 larger than the temperature of @xmath303 where the quantity @xmath283 goes below zero . thus , these estimates will always be consistent with the observed glass transition temperature of @xmath304 . moreover , we note that the low temperature behavior @xmath305 as indicated by the dotted line in fig . [ fig_ump2f0 ] , appears to be roughly related to the behavior of @xmath273 ( dotted line in fig . [ fig_logfactors ] ) in the same temperature range , by a single scaling factor of approximately @xmath306 . thus , it is possible that @xmath307\ ] ] if it turns out that @xmath308 for @xmath309 . if this is the case , then it means the procedure we used to estimate the pinch energy in sec . [ sec_proofglass ] is quantitatively correct , implying that the ground state of a random rna sequence indeed consists of the matching of rare segments independently at each length scale . it will be useful to pursue this analysis further using a renormalization group approach similar to what was developed for the denaturation of heterogeneous dnas by tang and chat  @xcite .    : the symbols show numerical estimates of the quantity @xmath283 for different temperatures . the estimates are obtained by averaging the numerically determined free energy of @xmath113 random sequences of length @xmath114 generated by the sequence disorder model eq . ( [ eq_energychoice ] ) with @xmath107 . the low temperature behavior can be described reasonably well by the expression @xmath310 ( dotted line ) . the consistency condition  ( [ eq_consistencycondition ] ) for the molten phase breaks down when this line intersects @xmath311 ( the dashed line ) . this yields @xmath312 as a lower bound for the glass transition temperature of this system . ] in this manuscript , we studied the statistical properties of random rna sequences far below the denaturation transition so that bases predominantly form base pairs . we introduced several toy energy models which allowed us to perform detailed analytical and numerical studies . through a two - replica calculation , we show that sequence disorder is perturbatively irrelevant , i.e. , an rna molecule with weak sequence disorder is in a molten phase where many secondary structures with comparable total energy coexist . a numerical study of the model at high temperature recovers scaling behaviors characteristic of the molten phase . at very low temperatures , a scaling argument based on the extremal statistics of rare matches suggest the existence of a different phase . this is supported by extensive numerical results : forced deformations are introduced by pinching distant monomers along the backbone ; the resulting excitation energies are found to grow very slowly ( i.e. , logarithmically ) with the deformation size . it is likely that the low temperature phase is a marginal glass phase . the intermediate temperature range is also studied numerically . the transition between the low temperature glass phase and the high temperature molten phase is revealed by a change in the coefficient of the logarithmic excitation energy , from being disorder dominated to entropy dominated .    from a theoretical perspective , it would be desirable to find an analytical characterization of the low temperature phase . if the excitation energy indeed diverges only logarithmically , one has the hope that this may be possible , e.g. , via the replica theory , as was done for another well known model of statistical physics  @xcite . it should also be interesting to include the spatial degrees of freedom of the polymer backbone ( via the logarithmic loop energy ) , to see how sequence disorder affects the denaturation transition . another direction is to include sequence _ design _ which biases a specific secondary structure , e.g. , a stem - loop  @xcite . from a numerical point of view , it is necessary to perform simulations with realistic energy parameters to assess the relevant temperature regimes and length scales where the glassy effect takes hold . to make potential contact with biology , one needs to find out whether a molten phase indeed exists between the high temperature denatured phase and the low temperature glass phase for a real random rna molecule , and which phase the molecule is in under normal physiological condition . finally , it will be very important to perform _ kinetic _ studies to explore the dynamical aspects of the glass phase . despite the apparent weakness of the thermodynamic glassiness , the kinetics at biologically relevant temperatures is expected to be very slow for random sequences  @xcite . the authors benefitted from helpful discussions with u.  gerland and d.  moroz . t.h .  acknowledges an earlier collaboration with d.  cule which initiated this study , and is indebted to l .- h . tang for a stimulating discussion during which the simple picture of sec .  [ sec_proofglass ] emerged . this work is supported by the national science foundation through grants no . dmr-9971456 and dbi-9970199 . the authors are grateful to the hospitality of the institute for theoretical physics at uc santa barbara where this work was completed . before we describe the exact solution for the two - replica problem , as defined by the partition function @xmath141 in eq . ( [ eq_gintermsofloops ] ) and the bubble weight @xmath137 in eq . ( [ eq_defq ] ) , we first provide here a heuristic derivation of the qualitative results . this mainly serves to give a flavor of the two - replica problem in the language of theoretical physics .    to this end , we define the quantity @xmath313 to be the partition function over all two replica configurations of a sequence of length @xmath314 under the constraint that base @xmath62 and base @xmath314 form a common bond . it is easy to see that @xmath315 where we set @xmath316 thus , the critical behavior of @xmath317 is identical to the critical behavior of @xmath313 which we will study in the following .    due to the no pseudo knot constraint of the secondary structures , @xmath313 has a very simple structure , @xmath318 as illustrated in fig .  [ fig_pirecurs ] .    : while the first and last of the @xmath7 bases described by @xmath319 always form a base pair this outermost base pair can be followed by a loop with @xmath320 outgoing stems . each of the stems is described by @xmath321 itself while the loop is characterized by its total length that can be split into the different pieces @xmath135 in different ways . ]    to simplify the above equation , it is useful to introduce the @xmath77-transforms @xmath322&=&\sum_{n=1}^\infty \pi(n ) e^{-\mu n}\\ \qhat[y]&=&\sum_{n=1}^\infty \tworq(n ) e^{-yn}\end{aligned}\ ] ] of @xmath321 and @xmath137 . now applying the @xmath77-transform to both sides of eq . ( [ eq_pi.1 ] ) , we obtain @xmath322e^{2\mu}\!\!\!&=&\!\!\ ! ! \int\!\ ! \frac{\mathrm{d}y}{2\pi i } \qhat[y ] \khat[y,\mu ] \big [ 1\ ! + \ ! \khat[y,\mu]\pihat[\mu]e^{2\mu}e^{-y } \label{eq_pi.2}\\\nonumber & & \hspace*{25 mm } + \ ! \left(\khat[y,\mu]\pihat[\mu]e^{2\mu}e^{-y}\right)^2 \!\!\!\!+\ldots \big]\end{aligned}\ ] ] where @xmath323 = \frac{1}{e^{-y+\mu}-1},\ ] ] and the inverse transform @xmath324 e^{y\ell}$ ] was used . ( [ eq_pi.2 ] ) can be simplified greatly to the following form , @xmath325e^{2\mu } = \qhat\cdot \int\frac{\mathrm{d}y}{2\pi i } \frac{\qhat[y]}{\khat^{-1}[y,\mu]-\pihat[\mu]e^{2\mu}e^{-y}}. \label{eq_pi}\ ] ] this is reminiscent of the well - known hartree solution to the @xmath326-theory , or equivalently the self - consistent treatment of the self - interacting polymer problem  @xcite , if we identify @xmath327 as the interaction parameter , @xmath328 $ ] as the  propagator " . the usual form of the hartree equation @xmath325 = \qhat\cdot \int \frac{\mathrm{d}^dk}{(2\pi)^d } \ , \frac{1}{k^2 + \mu -\pihat[\mu ] } \label{eq_pi1}\ ] ] corresponds to the small-@xmath329 , small-@xmath330 limit of eq . ( [ eq_pi ] ) , with @xmath331 playing the role of the square of the  wave number " @xmath332 . note that @xmath333 $ ] plays the role of the density of ( spatial ) states , i.e. @xmath334 = \frac{\mathrm{d}^dk}{(2\pi)^d}$ ] , where @xmath335 denotes the _ dimensionality _ of the  embedding space " .    in the context of rna , de gennes used this approach to describe the denaturation of uniformly attracting rna more than 30 years ago  @xcite . recently , this approach has been extended by moroz and hwa to study the phase diagram of rna structure formation  @xcite . the analysis of a self - consistent equation of the type ( [ eq_pi ] ) is well known  @xcite . the analytical properties of @xmath336 $ ] depend crucially on the form of @xmath333 $ ] . let the singular part of @xmath337 be @xmath338 where @xmath339 is the position of the singularity of @xmath333 $ ] . [ note that @xmath340 by comparing the forms of the hartree equations ( [ eq_pi ] ) and ( [ eq_pi1 ] ) . ] for @xmath341 , there is only one solution for all @xmath342 , with a square root singularity in @xmath336 $ ] at some finite value of @xmath330 . for @xmath343 , there are two possible solutions depending on the value of @xmath327 . the square - root singularity exists for @xmath327 exceeding some critical value depends through @xmath337 on @xmath124 but not on @xmath125 . ] @xmath344 , while for @xmath345 , the square - root singularity disappears and @xmath336 $ ] is governed by the singularity of @xmath337 given in eq . ( [ eq_qsing ] ) . performing the inverse transform and using eq . ( [ eq_gandpi ] ) , we get @xmath346 where @xmath144 is a non - universal parameter given by the location of the singularity , while the exponent @xmath347 characterizes the phase of the system and is given by the singularity of @xmath336 $ ] : we have @xmath348 if @xmath336 $ ] is dominated by the square - root singularity and @xmath349 if @xmath336 $ ] is dominated by @xmath337 . the interpretation of the two phases with @xmath348 and @xmath349 are straightforward : the phase with @xmath350 describes the usual rna secondary structure ( see eq . ( [ eq_moltenz ] ) ) ; here the bubbles described by @xmath137 are irrelevant . in the other phase , the result that @xmath349 indicates that base pairing is not relevant and the system behaves as a single bubble . in the context of the original two - replica problem , the irrelevancy of the bubbles in the @xmath350 phase indicates that the two replicas are locked together , behaving as a single replica in this phase . in the other phase , the attraction of the common bonds is irrelevant , and the two replicas become independent of each other .    as explained in sec . [ sec_replica ] , the purpose of the two - replica calculation is to determine whether the inter - replica attraction , characterized by @xmath327 here , is irrelevant , i.e. , whether the system will not yet be in the @xmath350 phase for a value of @xmath351 . this is only possible if @xmath352 . from the solution of the problem described above , this depends crucially on the singularity of @xmath337 , specifically , on whether @xmath353 . the difficulty in ascertaining the form of @xmath337 lies in the no - common - bond constraint ( i.e. , @xmath354 ) in the definition of @xmath137 ( [ eq_defq ] ) . however , we note that for @xmath355 , the two replica partition function @xmath356 is simply the square of the single replica partition function @xmath74 . thus , @xmath357 according to eq . ( [ eq_moltenz ] ) . since we just convinced ourselves that @xmath347 can take on only two possible values , namely @xmath358 and @xmath359 and since @xmath360 we conclude @xmath361 and moreover @xmath362 . thus , we do expect the phase transition to occur at @xmath363 . however , it is not clear from this calculation if the system at @xmath355 ( or @xmath364 ) is exactly at or strictly below the phase transition point . we leave it to the exact solution of the two replica problem presented in the next two appendices to establish that @xmath364 is indeed strictly below the phase transition point and that therefore disorder is perturbatively irrelevant . we note that in the context of the @xmath326-theory or the self - consistent treatment of the self - interacting polymer , the result @xmath365 implies that the embedding spatial dimension is @xmath366 . thus , the two - replica problem corresponds to the denaturation of a single rna in @xmath0 spatial dimensions . the bubbles @xmath137 s of fig . ( [ fig_combonds ] ) which originate from the branching entropy of the individual rnas play the role of spatial configurational entropy of the single - stranded rna in the denaturation problem . in this appendix , we present the exact solution of the two replica problem . while most of the details are given here , the most laborious part is further relegated to app .  [ app_derivez1 ] . we start by introducing an auxiliary quantity @xmath367 . this is a restricted two - replica partition function , summing over all independent secondary structures of a pair of rnas of length @xmath314 bases in which there are exactly @xmath368 exterior bases of the common bond structure without disrespecting the no pseudo - knot constraint . ] all of which are completely unbound in both replicas . since the exterior bases form one of the bubbles of the common bond structure , the possible binding configurations of these exterior bases are described by @xmath369 . thus , the full partition function of the two replica problem can be calculated from this restricted partition function as @xmath370 now , let us formulate a recursion relation for @xmath371 by adding one additional base @xmath7 to each of the two rnas . we can separate the possible configurations of the new function @xmath372 according to the possibilities that the new base @xmath7 is either not involved in a common bond or forms a common bond with base @xmath373 . this yields the recursion relation @xmath374 for @xmath375 and @xmath376 . the applicable boundary conditions are : @xmath377 , @xmath378 , and @xmath379 for each @xmath380 and @xmath381 .    at this point , it is convenient to introduce the @xmath77-transforms in order to decouple the discrete convolution in eq . ( [ eq_trrecurs ] ) . they are @xmath382 using eq . ( [ eq_trrecurs ] ) and the boundary conditions we get @xmath383 this can be solved for @xmath384 with the result @xmath385 together with the boundary condition @xmath386 , we get @xmath387^{-n}.\ ] ] if we now multiply eq . ( [ eq_gbywandq ] ) by @xmath388 and sum both sides over @xmath7 we get @xmath389 which , upon inserting eq . ( [ eq_hatw ] ) , becomes an implicit equation @xmath390 for the full partition function @xmath391 , provided that we know the function @xmath392 . since @xmath392 does not depend on @xmath125 , we can find its form using the following strategy : if @xmath355 , a common bond does not contribute any additional boltzmann factor . thus , the two replica partition function for this specific value of @xmath125 is just the square of the partition function of a single uniformly attracting rna molecule , i.e. , @xmath393 since we know @xmath74 through the exact expression ( [ eq_singlez ] ) for its @xmath77-transform @xmath394 , we can regard @xmath395 as a known function , even though a closed form expression is not available . from eq . ( [ eq_trimplicit ] ) , we have @xmath396 this is an equation for @xmath392 in terms of the known function @xmath395 . after we solve it for @xmath392 below , we can use eq . ( [ eq_trimplicit ] ) to solve for the only leftover unknown @xmath391 for arbitrary values of @xmath125 .      in the thermodynamic limit , it is sufficient to consider only the singularities of the @xmath77-transform @xmath391 . from the form of @xmath397 in the vicinity of the singularity @xmath398 , the two - replica partition function @xmath317 is readily obtained by the inverse @xmath77-transform , with the result @xmath399 the result immediately yields the free energy per length , @xmath400 . more significantly , the exponent @xmath143 reveals which phase the two - replica system is in : for @xmath355 ( i.e. no disorder ) , the two - replica system is just a product of two independent single - replica systems and we must have @xmath148 as implied by the single - replica partition function @xmath74 in eq . ( [ eq_moltenz ] ) . on the other hand , for @xmath401 , the two replicas are forced to be locked together and behave as a single replica . in this case , we must have @xmath402 . as we will see , @xmath148 and @xmath403 are the only values this exponent can take on for this system ; it indicates whether or not the two replicas are locked , and hence whether or not the effect of disorder is relevant . the singularity @xmath398 of @xmath391 is given implicitly by eq . ( [ eq_trimplicit ] ) , which we now analyze in detail . we start by recalling the solution of the homogeneous single rna problem , eq . ( [ eq_moltenz ] ) . from the relation ( [ eq_g2 ] ) , we have @xmath404 for large @xmath7 , with @xmath82 , @xmath405 , and @xmath406 given in sec . [ sec_moltenpartfunc ] . hence , the @xmath77-transform @xmath407 is defined on the interval @xmath408 . it is a monotonously decreasing function of @xmath77 , terminating with a singularity at @xmath409 which produces the @xmath148 singularity in @xmath356 .    due to ( [ eq_trimplicitqt1 ] ) , the same singularity must occur in @xmath410 at @xmath411 , where @xmath412 is a positive number and does not depend on anything else but @xmath124 . since @xmath413 is a smooth monotonously increasing function which maps the interval @xmath408 into the interval @xmath414 , it follows from eq . ( [ eq_trimplicitqt1 ] ) that @xmath410 is a smooth , monotonously decreasing function which maps the interval @xmath414 into the interval @xmath4150,\gatsing]$ ] . now that we have characterized @xmath410 in detail , we can proceed to study @xmath391 for arbitrary @xmath125 . clearly , according to eq . ( [ eq_trimplicit ] ) @xmath391 has a singularity leading to @xmath148 at @xmath416 , defined implicitly by @xmath417 because @xmath410 has this singularity at @xmath418 . again according to eq . ( [ eq_trimplicit ] ) , we have @xmath419 independent of @xmath125 . this leads to one of the key results @xmath420 if @xmath416 is the only singularity of @xmath391 , it would imply that there is only one phase with @xmath148 , and the free energy per length of the two - replica system is given by @xmath421\ ] ] for all values of @xmath125 . by differentiating this with respect to @xmath125 , we obtain the fraction of common contacts @xmath422 in this phase as a function of @xmath125 . for very large disorder , i.e. , for large @xmath125 , this fraction converges to one . however , since it is the fraction of bonds divided by the total number of bases and every base pair has two bases , it has to be bounded from above by @xmath423 . thus , we conclude that eq . ( [ eq_freeenergy0 ] ) can not be the free energy of the two - replica system for all @xmath125 s . at least for large @xmath125 , there must be another singularity of @xmath391 which will yield a different expression for the free energy to give physically reasonable fraction of common bonds .    in order to find this other singularity , we introduce the inverse function @xmath424 of @xmath391 . from eq . ( [ eq_trimplicit ] ) , it follows that for any @xmath125 and any @xmath4250,\gatsing]$ ] , @xmath426\\ \rightarrow\quad g&= & \widehat{\tworq}[\widehat{\tworz}(g;\widetilde{q})-q^2\widetilde{q } g]\\ \rightarrow\quad \widehat{\tworq}^{-1}(g)&= & \widehat{\tworz}(g;\widetilde{q})-q^2\widetilde{q } g.\end{aligned}\ ] ] since @xmath427 does not depend on @xmath125 , we can eliminate it by evaluating the last equation above at the special value @xmath355 and write @xmath428 where @xmath429 is the inverse of the known function @xmath395 . ( [ eq_zrelation ] ) is now an explicit solution for the inverse of @xmath391 for _ arbitrary _ @xmath125 , and the singularity of @xmath430 , located at @xmath431 and @xmath432 , yields the free energy of the two - replica system , i.e. , @xmath433 , in this phase .    while app . [ app_derivez1 ] derives the position of the dominant singularity present in eq . ( [ eq_zrelation ] ) rigorously , we will resort to some intuitive argument here . since @xmath434 is a monotonously decreasing , convex function , so is its inverse @xmath429 . the latter function ends at the point @xmath435 with some slope @xmath436 in a singularity which produces the @xmath148 behavior indicative of two independently fluctuating uniformly self - attracting replicas . this is shown in fig . [ fig_zhatsketch ] as the solid line . according to eq . ( [ eq_zrelation ] ) , we can obtain the corresponding function @xmath424 for arbitrary @xmath125 by simply adding a linear function @xmath437 to this function . if the slope @xmath438 of this linear function is less than the smallest slope of @xmath429 , i.e. , if @xmath439 , adding this linear function does not qualitatively change anything ( see the short dashed line in fig .  [ fig_zhatsketch ] . ) the only singularity is still the one at @xmath440 , the corresponding value in @xmath77 , i.e. , @xmath441 , is trivially shifted by an amount @xmath442 from @xmath443 as @xmath125 varies . thus , the scaling behavior is characterized by @xmath148 as if @xmath355 ( i.e. , absence of disorders ) , although the free energy @xmath444 is shifted as already derived in eq . ( [ eq_freeenergy0 ] ) .     of the two replica system at various values of the common bond interaction @xmath125 : the solid line shows the free system without any interaction of common bonds . in the presence of an interaction , the inverse function of the partition function can be obtained by adding a linear function to the free system function . if the interaction is not too strong ( @xmath445 , short dashed line ) adding the linear function with a small slope does not change the qualitative form of the partition function . in this case , the two replica system is controlled by the singularity at @xmath440 which is independent of @xmath125 . at stronger interactions ( @xmath446 , long dashed line ) the inverse function develops a minimum . beyond this minimum it is not invertible any more and the two replica system is then dominated by the singularity arising from this minimum . ]    if the final slope @xmath447 of the right hand side of eq . ( [ eq_zrelation ] ) is _ positive _ the situation becomes much different : upon adding the linear function , @xmath424 develops a minimum at some position @xmath448 . thus , the inverse function @xmath391 has to be calculated from the left ( small @xmath449 ) branch of @xmath424 and has a _ square root singularity _ at @xmath450 . this square root singularity implies that the characteristic exponent becomes @xmath403 , consistent with the picture that in this phase the two replicas are locked together and fluctuate as one single effective rna molecule . the position @xmath451 of the minimum is determined by the root of the derivative , i.e. , by @xmath452^{-1}\!\!\!\!\!\!\!\!.\ ] ] the corresponding value @xmath453 determines the location of the square root singularity @xmath450 of @xmath391 , i.e. , the free energy per length of the two replica problem . @xmath450 can be conveniently expressed in terms of the auxiliary quantity @xmath454 defined through @xmath455 as @xmath456 comparing this expression with eq . ( [ eq_freeenergy0 ] ) which is valid for small @xmath125 , we can summarize the complete solution in terms of @xmath457\singlezsing^2,\infty[$ that ful-}\   & \widetilde{q}>\widetilde{q}_{\mathrm{c}}\\ \mbox{fils } \frac{\mathrm{d}}{\mathrm{d}z}\widehat{\tworg}(z;1 ) \!=\!\!-\!1\!/[q^2(\widetilde{q}\!-\!1)]&\\[10pt ] \singlezsing^2=(1 + 2\sqrt{q})^2&\widetilde{q}\le\widetilde{q}_{\mathrm{c } } \end{array}\right.\ ] ] where @xmath458 and @xmath459    in terms of this @xmath460 , the smallest singularity of @xmath461 is located at @xmath462 the free energy per length of the two replica system is given by @xmath463\ ] ] and the fraction of common bonds is @xmath464 this fraction of common bonds turns out to be continuous at the phase transition but it exhibits a jump in its slope at @xmath465 .    in the case @xmath466 , these simplify to eqs . ( [ eq_freeenergy0 ] ) and  ( [ eq_contacts0 ] ) ( or eqs . ( [ eq_moltenzinmaintext ] ) and  ( [ eq_defginmaintext ] ) respectively ) , with @xmath467 independent of @xmath125 . the type of singularity of the laplace transformed partition function is the same as at @xmath355 , resulting in @xmath148 . for @xmath446 , we can not write down a closed form expression for @xmath398 any more . but it is given implicitly in terms of the solution of eq . ( [ eq_zkdef ] ) ; it involves only single replica quantities and can thus be evaluated numerically . moreover , we have seen that the dominant singularity is in this regime a square root singularity implying @xmath403 . in this appendix we give a derivation of the position of the non - trivial singularity in the laplace transform of the partition function @xmath391 . a more intuitive , graphical derivation of this result was given in app . [ app_solution ] . using eq . ( [ eq_zrelation ] ) , we start by calculating @xmath468^{-1}\\ & = & \left[\left.\frac{\mathrm{d}}{\mathrm{d}g } \right|_{g=\widehat{\tworg}(z;\widetilde{q})}\!\!\!\ ! \widehat{\tworz}(g;1)+q^2(\widetilde{q}-1 ) \right]^{-1}\\ & = & \left[\left[\left.\frac{\mathrm{d}}{\mathrm{d}z } \right|_{z=\widehat{\tworz}(\widehat{\tworg}(z;\widetilde{q});1)}\!\!\!\!\widehat{\tworg}(z;1 ) \right]^{-1}\!\!\!\!\!\!\!\!\!+q^2(\widetilde{q}-1)\right]^{-1}\!\!\!\!.\end{aligned}\ ] ] this expression obviously has a singularity at @xmath450 which is defined by @xmath469 since @xmath470 , this is only possible for @xmath471 for smaller values of @xmath125 , there is no other singularity and the free energy per length is given by eq . ( [ eq_freeenergy0 ] ) .    if @xmath472 the additional singularity @xmath450 exists and is  as we will see below  always smaller than the singularity @xmath473 . thus , the free energy per length is given by the singularity @xmath450 in the strong coupling phase , i.e. , for @xmath472 . on the first sight , eq .  ( [ eq_z1def ] ) still looks as if @xmath450 could only be calculated if the full function @xmath391 is known . however , for any @xmath472 we can define @xmath454 to be the unique solution of the equation @xmath474 this quantity depends only on the function @xmath434 . according to eq . ( [ eq_z1def ] ) , @xmath450 and @xmath454 are related by @xmath475 . this implies that @xmath476 . on the other hand , eq . ( [ eq_zrelation ] ) applied to @xmath477 yields @xmath478 which is finally an expression which involves only quantities of the non - interacting system . since @xmath479 is a monotonous function on the interval @xmath480 $ ] , we always have @xmath481 with equality if and only if @xmath467 , i.e. , for @xmath482 . therefore , the free energy per length is indeed given by @xmath483\ ] ] for any @xmath472 . 99 see , e.g. , w. saenger , _ principles of nucleic acid structure _ ( springer - verlag , new york , 1984 ) or b. alberts _ et al . _ , _ molecular biology of the cell _ ( garland publishing , new york , 1994 ) . see , e.g. , r.f . gesteland and j.f . atkins _ eds . _ , _ the rna world : the nature of modern rna suggests a prebiotic rna world _ ( cold spring harbor laboratory press , cold spring harbor , 1993 ) . k.a . dill _ et al . _ , _ protein sci . _ * 4 * , 561 ( 1995 ) . onuchic , z. luthey - schulten , and p.g . wolynes , _ annu . chem . _ * 48 * , 545 ( 1997 ) . t. garel , h. orland , and e. pitard , _ j. phys . i _ * 7 * , 1201 ( 1997 ) . shakhnovich , _ curr . op . biol . _ * 7 * , 29 ( 1997 ) . for a review see p.g . higgs , _ quart . biophys . _ * 33 * , 199253 ( 2000 ) . m. zuker and p. stiegler , _ nucl . res . _ * 9 * , 133148 ( 1981 ) . mccaskill , _ biopolymers _ * 29 * , 11051119 ( 1990 ) . hofacker , w. fontana , p.f . stadler , s. bonhoeffer , m. tacker , and p. schuster , _ monatshefte f. chemie _ * 125 * , 167188 ( 1994 ) . freier , r. kierzek , j.a . jaeger , n. sugimoto , m.h . caruthers , t. neilson , and d.h . turner , _ proc . usa _ * 83 * , 93739377 ( 1986 ) . i. tinoco jr and c. bustamante , _ biol . _ * 293 * 271281 ( 1999 ) and references therein . b. maier , d. bensimon , and v. croquette , _ proc . usa _ * 97 * , 1200212007 ( 2000 ) . de gennes , _ biopolymers _ * 6 * , 715 ( 1968 ) . higgs , _ phys . lett . _ * 76 * , 704707 ( 1996 ) . r. bundschuh and t. hwa , _ phys . lett . _ * 83 * , 14791483 ( 1999 ) . a. pagnani , g. parisi , and f. ricci - tersenghi , _ phys . * 84 * , 20262029 ( 2000 ) . hartmann , _ phys . _ * 86 * , 13821382 ( 2001 ) . a. pagnani , g. parisi , and f. ricci - tersenghi , _ phys . rev . lett . _ * 86 * , 13831383 ( 2001 ) . see , e.g. , p. young _ ed . _ , _ spin glasses and random fields _ ( world scientific , singapore , 1998 ) . d. moroz and t. hwa , in preparation . waterman , _ advan . in math . studies _ * 1 * , 167212 ( academic press , new york , 1978 ) . see , e.g. , m. mitas , _ nucl . acids res . _ * 25 * , 2245 - 2254 ( 1997 ) . see , e.g. , w. feller , _ an introduction to probability theory and its applications _ ( john wiley & sons , new york , 1950 ) . lubensky , j.  isaacson , and s.p . obukhov , _ j. physique _ * 42 * , 15911601 ( 1981 ) . s.f . edwards and p.w . anderson , _ _ * f5 * , 965 ( 1975 ) . see , e.g. , m. mzard , g. parisi , and m.a . virasoro , _ spin glass theory and beyond _ ( world scientific , singapore , 1986 ) . r. arratia , l. gordon , and m. waterman , _ ann . stat . _ * 18 * , 539570 ( 1990 ) . d.a . huse and d.s . j. physique i _ * 1 * , 621625 ( 1991 ) . morgan and p.g . * 105 * , 71527157 ( 1996 ) . morgan and p.g . j. phys . a : math . gen . _ * 32 * , 31533170 ( 1998 ) . c.  flamm , w.  fontana , i.l . hofacker , and p.  schuster , _ rna _ * 6 * , 325338 ( 2000 ) . h.  isambert and e.d . siggia , _ proc . . usa _ * 97 * , 65156520 ( 2000 ) . u. gerland , r. bundschuh , and t. hwa , in preparation . cardy and s.  ostlund , _ phys . b _ * 25 * , 68996909 ( 1982 ) . s.e . korshunov , _ phys . b _ * 48 * , 39693975 ( 1993 ) . t.  hwa and d.  fisher , _ phys . rev . b _ * 49 * , 31363154 ( 1994 ) . tang and h. chat , _ phys . lett . _ * 86 * , 830 ( 2001 ) . h.  isambert , personal communication . j. zinn - justin , _ quantum field theory and critical phenomena _ ( oxford university press , united kingdom , 1989 ) .
Abstract:

Answer: the formation of secondary structures by a random rna sequence is studied as a model system for the sequence - structure problem omnipresent in biopolymers . several toy energy models are introduced to allow detailed analytical and numerical studies . first , a two - replica calculation is performed . by mapping the two - replica problem to the denaturation of a single homogeneous rna in @xmath0-dimensional embedding space , we show that sequence disorder is perturbatively irrelevant , i.e. , an rna molecule with weak sequence disorder is in a _ molten phase _ where many secondary structures with comparable total energy coexist . a numerical study of various models at high temperature reproduces behaviors characteristic of the molten phase . on the other hand , a scaling argument based on the extremal statistics of rare regions can be constructed to show that the low temperature phase is unstable to sequence disorder . we performed a detailed numerical study of the low temperature phase using the droplet theory as a guide , and characterized the statistics of large - scale , low - energy excitations of the secondary structures from the ground state structure . we find the excitation energy to grow very slowly ( i.e. , logarithmically ) with the length scale of the excitation , suggesting the existence of a marginal glass phase . the transition between the low temperature glass phase and the high temperature molten phase is also characterized numerically . it is revealed by a change in the coefficient of the logarithmic excitation energy , from being disorder dominated to entropy dominated .