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: study of the conformational transitions of heteropolymers in dilute solutions is important for many applications from the chemical industry to biotechnology . directed more towards the former , there has been a significant amount of theoretical work carried out on concentrated copolymer solutions , mixtures and blends using various types of the density formalism . however , these approaches are not valid at infinitely low dilution where the fundamental interactions of the individual macromolecule determine its conformational state . this situation is more relevant for biochemistry . the problem is even harder to address at non  equilibrium conditions typical for biopolymers _ in vivo _ @xcite .    for these reasons we have been working for some time on developing an adequate statistical mechanical technique for studying the equilibrium structure and kinetics across phase transitions in heteropolymers @xcite . the main idea was to extend the gaussian self  consistent ( gsc ) method , originally proposed for the homopolymer ( see e.g. @xcite and references therein ) , to the case of inhomogeneous monomer interactions . this has been achieved in ref . @xcite where we have derived the complete set of non  linear kinetic equations for complex valued equal  time correlation functions of the fourier transforms of the monomer coordinates . there we have analysed in some detail the simplest periodic ( ab ) copolymer , but study of more complicated sequences remained out of our practical computational reach . the relation of the kinetic equations to the equilibrium free energy , as well as the expression for the entropy , were also unknown to us at that stage . thus , the phase diagram have not been elucidated and certain other important issues have not been addressed in that work .    in the following two works @xcite we have achieved further progress and resolved most of these questions , though at the cost of loosing detailed information about individual sequences . namely , we have performed averaging of the gsc equations over the quenched disorder in the monomer amphiphilic strength . this yielded a closed set of kinetic equations for description of random copolymers with a gaussian distribution of the disorder . such an idealised system is very interesting in itself , particularly since there is a hope that its study could shed light on some features of the extremely complex problem of protein folding @xcite .    to render the quenched disorder problem more tractable certain perturbative approximation was necessary @xcite . it became clear in ref . @xcite that this causes some deficiencies in the equilibrium limit of the formalism , and we have alluded as to how these could be alleviated in higher orders of the expansion .    therefore , it seems important to revisit the general formalism of ref . @xcite in order to resolve remaining difficulties , maintaining the rich information about particular sequences and avoiding any further approximations . furthermore , there exists , till now , no simple theoretical procedure capable of giving the equilibrium conformational states of a heteropolymer with arbitrary sequence , that can in a consistent manner also give the full kinetic pathway between these states . the work presented here achieves this . to demonstrate the strength of the extended gsc method we consider a number of interesting examples of heteropolymer sequences . different kinds of interactions , such the hydrodynamic interaction , may be straightforwardly incorporated into the scheme . a rich and nontrivial physical picture emerges as a result of this theoretical progress . traditionally , we proceed from the coarse  grained description of the polymer chain @xcite with the spatial coordinates @xmath0 ascribed to the @xmath1-th monomer . it is assumed that the long timescale evolution of conformational changes is well represented by the phenomenological langevin equation , which upon neglecting the backflow effect of the solvent may be written as , @xmath2 where @xmath3 is the `` bare '' friction constant per monomer . for the discussion of the hydrodynamic interaction we refer the reader to appendix [ app : hydro ] . the thermal fluctuations are incorporated via the gaussian noise which , according to the einstein law , is characterised by the second momentum , @xmath4 where the greek indices denote the spatial components of 3-d vectors .    in the current treatment the solvent is effectively excluded from the consideration @xcite and the resulting monomer interactions are described by the effective free energy functional , @xmath5 where in principle @xmath6 are allowed to have any dependence on the site indices @xmath7 .    the first two terms in eq . ( [ eq : h ] ) describe the connectivity and the stiffness of the chain . their coefficients have the following simple meaning : @xmath8 and @xmath9 with @xmath10 and @xmath11 called the statistical segment length and the persistent length @xcite respectively . apart from these interactions , local along the chain , there are also long  range volume interactions represented by the virial  type expansion @xcite in eq . ( [ eq : h ] ) . the latter reflects the hard core repulsion and weak attraction between monomers , but also the effective interaction mediated by the solvent  monomer couplings @xmath12 . the coefficients in this expansion may be calculated as functions of the temperature and the parameters of molecular interaction . however , for our purposes we do not need to know their explicit form here as we shall keep only a few first terms . appropriate coefficients then may be viewed as independent phenomenological parameters which could be directly related to experimentally measurable quantities .    in refs . @xcite we have discussed that the case of amphiphilic heteropolymers , for which monomers differ only in the monomer  solvent coupling constants , corresponds to the following choice of site  dependent second virial coefficients in eq . ( [ eq : h ] ) , @xmath13 the mean second virial coefficient @xmath14 is associated with the quality of the solvent : positive values correspond to the good solvent ( where effective two  body repulsion of monomers results in the extended coil conformation ) , and the negative values correspond to the bad solvent condition ( where the effective two  body attraction tends to compact the chain ) .    the set of couplings @xmath15 expresses the chemical composition , or using the biological terminology , the primary sequence of a heteropolymer . for simplicity , in the consequent sections we shall consider examples for which these constants can be parametrised in a simple way : @xmath16 , where variables @xmath17 can only take three values : @xmath18 or @xmath19 corresponding to ( a ) hydrophobic , ( b ) hydrophilic and ( c ) `` neutral '' monomers respectively . the parameter @xmath20 is called the degree of amphiphilicity of the chain . note that for more complicated than binary sequences there is another relevant dispersion , @xmath21 , and the combination @xmath22 is a more appropriate variable . in work @xcite we have derived a set ( 18 ) of the gsc kinetic equations for the equal  time correlation functions of the fourier transforms of monomer coordinates @xmath23 . there , at the end of sec . ii , we have mentioned that a polymer with no periodic structure may be described by choosing the number of blocks @xmath24 . the equations for the correlation functions , @xmath25 contain some redundancies and also are intermixed with the diffusive degree of freedom describing the motion of the centre  of  mass . obviously , for a single chain the latter can be easily decoupled from the intra  molecular degrees of freedom by introducing the mean squared distances between pairs of monomers , @xmath26 it turns out that from eq . ( 18 ) in ref . @xcite after certain algebraic manipulations one can obtain a closed set of equations for the quantities @xmath27 . taking into account the additional bending energy contribution in eq . ( [ eq : h ] ) , the gsc equations may be written down in the most general form as follows , @xmath28 where we have introduced the four  point correlation functions , @xmath29 , and the matrix @xmath30 of size @xmath31 with the cofactor @xmath32 , @xmath33 similarly to eq . ( 31 ) in ref . @xcite for the mean energy @xmath34 we have , @xmath35    it is interesting to note that in the case of the homopolymer the right  hand side of the kinetic equations ( [ eq : dkineq ] ) may be rewritten via the instantaneous gradients of the variational free energy by introducing the normal modes ( see eq . ( 5 ) in ref . this establishes connection of the stationary limit of our kinetic approach with the equilibrium theory , in which we recover the gibbs  bogoliubov variational principle . we also note that if first order phase transitions are involved , one has to possess the expression for the free energy in order to determine the phase boundaries by finding the global free energy minimum . the variational free energy , @xmath36 , based on the gaussian _ ansatz _ for equal  time pair correlation functions , contains two terms , @xmath37 . naturally , the mean energy term is given by eq . ( [ eq : e ] ) . the second , entropic contribution , is calculated in appendix [ app : var ] . let us summarise the result of that calculation here . one representation for the entropy that is suitable for numerical analysis is , @xmath38 where we have used the determinant of the @xmath39-dimensional major submatrix of the matrix , @xmath40 and the matrix @xmath41 obviously has the elements equal to @xmath42 . the reason for appearance of the truncated matrix is that we have excluded the zero eigenvalue of @xmath43 related to the translational invariance . here we have also introduced the @xmath39-dimensional orthogonal projector @xmath44 such that , @xmath45 with the matrix elements @xmath46 . this matrix has obviously one zero eigenvalue and @xmath47 degenerate eigenvalues equal to @xmath48 . however , for analytical treatment it is more convenient to obtain a slightly different expression by regularising the zero eigenvalue ( i.e. by imposing the constraint @xmath49 as a `` soft '' condition in eq . ( [ eq : tr ] ) ) , @xmath50 since the matrix @xmath51 is nondegenerate we can easily differentiate eq . ( [ se ] ) so that , @xmath52 where inside the trace @xmath53 over the @xmath39-dimensional subspace projected out by @xmath44 the matrix @xmath43 becomes invertible with the inverse denoted as @xmath54 . this allows us to obtain the combination which appears in the kinetic equations , @xmath55    these preliminaries are sufficient to prove the desired relation . indeed , using eq . ( [ eq : sder ] ) by direct and tedious differentiation of eq . ( [ eq : e ] ) one can show that the kinetic equations ( [ eq : dkineq ] ) may be expressed through the instantaneous gradients @xcite of the variational free energy as , @xmath56}{\partial d_{mm''}(t)}- \frac{\partial { \cal a}[d(t)]}{\partial d_{m'm''}(t ) } \right).\ ] ] this formula , together with eqs . ( [ eq : e],[se ] ) , is the key formal result of the current work . the structure of eq . ( [ eq : kinmain ] ) is sufficiently nontrivial to be guessed from phenomenological arguments and has been derived in a systematic manner proceeding from eq . ( [ eq : dkineq ] ) . we would like to comment here that although , for simplicity , we have presented the explicit formulae above for a ring polymer , our current formalism is general and covariant . in fact , the kinetic equations ( [ eq : kinmain ] ) are valid for any topology of the chain . thus , it is straightforward to consider more complicated topologies such as a star , brush , network , branched chain and so on . for that it is sufficient to modify only the spring and stiffness terms in eq . ( [ eq : h ] ) , and respectively in eq . ( [ eq : e ] ) . for example , to describe an open chain one has simply to suppress the term with @xmath57 in the connectivity contribution , and the terms with @xmath58 in the stiffness contribution to the energy . we undertake a detailed comparison of ring versus open homopolymers in kinetics at the collapse transition of the homopolymer in a separate work @xcite . another interesting possibility is that the general equation ( [ eq : kinmain ] ) is also valid for models with different ways of representing the connectivity and stiffness . for example , the freely rotating model @xcite can be obtained by suppressing the first two terms in eq . ( [ eq : e ] ) and instead keeping fixed the following mean squared distances : @xmath59 and @xmath60 , where @xmath61 and @xmath62 are the bond length and angle . this is easy to prove by adding appropriate `` soft '' constraints to the free energy functional and taking the consequent limit in eq . ( [ eq : kinmain ] ) , so that they become delta  functions in the partition function . finally , let us introduce two main observables : the mean squared radius of gyration and the degree of micro  phase separation @xcite , @xmath63 the second parameter has a meaning of the dimensionless correlation of the matrices of the relative two  body interaction and the mean squared distances . for heteropolymers with two types of monomers it characterises the difference between the mean squared radii of gyration of hydrophilic and hydrophobic species , and for the symmetric composition with their numbers equal we have a simple relation : @xmath64 . the self  consistent kinetic equations ( [ eq : kinmain ] ) have been studied numerically using the explicit formulae ( [ eq : c1],[eq : c2],[eq : c4 ] ) for the effective potentials ( see appendix [ app : der ] ) , and the expression given by the first term in the right  hand side of eq . ( [ eq : dkineq ] ) for the entropic contribution . we used the fifth order adaptive step runge  kutta method @xcite to improve stability of the solution which , for large amphiphilicity parameter , @xmath20 , is characterised by a rather rugged free energy landscape . we note that such kinetic method for finding the equilibrium distributions is more reliable and efficient than the standard methods of free energy minimisation if there are lots of mountains and valleys on its surface . we also refrain from study of the influence of the hydrodynamics in this paper for the analysis and results are complicated enough already . besides , the hydrodynamics in the preaveraged approximation does not affect the equilibrium state , which is recovered by taking the stationary limit in the kinetic equations . we include the volume interactions up to the three  body terms only , i.e. @xmath65 for @xmath66 . as can be seen from eqs . ( [ eq : kinmain],[eq : c2 ] ) the computational time per time step scales with the chain length as @xmath67 here . this performance is intermediate between that of the homopolymer @xmath68 @xcite and that of the random copolymer @xmath69 @xcite . the performance of the formalism in ref . @xcite was @xmath70 , where @xmath71 and @xmath72 are the block length and number of blocks respectively . besides , that formalism relied on the use of complex variables , and the unitary transformation to the real basis was not an easily automatised task for complicated sequences . moreover , the treatment of the diffusive mode was nontrivial and sequence dependent . thus , in every respect , the current scheme is most attractive for study of heteropolymer sequences from the computational point of view . it is natural to work here with combinations @xmath73 and @xmath74 as the units of size and time . we choose @xmath75 , @xmath76 and @xmath77 to fix @xmath78 and @xmath79 equal to unity . in addition , we fix the following interaction parameters : the third virial coefficients @xmath80 and the stiffness @xmath81 as we did in ref . @xcite . now , we turn to discussion of concrete results . we have studied ring chains of the lengths @xmath82 . we shall present here the most complete results for @xmath83 and discuss the @xmath84 dependence only briefly . we have examined many different sequences . however , we present our data in this paper only for three particular choices which we found most typical and illustrative of the heteropolymer behaviour . these are chains of : 30 ( ab ) blocks , 10 ( aaabbb ) blocks and 2 ( abacbbcabccbcaaacbcccbbaacbcca ) blocks , which we call the `` short '' blocks , `` long '' blocks and `` random '' sequences respectively ( see the end of section [ sec : model ] for the monomer notations ) . the phase diagrams in terms of the mean second virial coefficient @xmath85 and the amphiphilicity @xmath20 for the above mentioned three sequences are presented in figs . [ fig:1]-[fig:3 ] . for positive values of @xmath14 and comparatively small values of @xmath20 the conformational state of the chain is akin to a homopolymer extended coil ( see fig .  [ fig : dmneq]a ) . by decreasing @xmath14 to the negative region the chain is caused to undergo a continuous ( second  order  like ) collapse transition ( curve ( c ) ) , that is characterised by a rapid fall of the radius of gyration @xcite ( see fig .  1 of ref . @xcite ) and the change of the fractal dimension . proceeding from the collapsed globule at a fixed negative @xmath14 the increase of the amphiphilicity @xmath20 also causes a continuous transition ( curve ( s ) ) to the micro  phase separated ( mps ) globule @xcite . during this transition the system size , @xmath86 , monotonically increases ( see fig . [ fig:4 ] ) and the mean energy @xmath87 decreases in agreement with our earlier results in refs . the former change is more pronounced for the `` long '' blocks compared to other sequences , for which the connectivity constraints impede formation of structures with a hydrophilic shell and hydrophobic core . the mps order parameter @xmath88 ( see fig . [ fig:5 ] ) increases almost linearly for small @xmath20 , then after the transition asymptotically saturates . the transition from the coil at large values of @xmath20 to the mps globule turns out to be more complicated , and essentially dependent on the sequence . in case of the `` long '' blocks ( fig .  [ fig:2 ] ) the collapse transition to the mps globule becomes discontinuous ( first  order  like ) . thus , inside the boundaries of metastability , designated by spinodals i and i  , there are two competing minima of the free energy : the coil and the mps globule . the former minimum is characterised by a large size , @xmath86 , ( on the right of fig .  [ fig:6 ] ) and a small mps order parameter @xmath88 ( fig . [ fig:7 ] ) , while for the latter minimum the situation is reversed ( see the left hand side of figs .  [ fig:6],[fig:7 ] ) . the depths of the free energy minima become exactly equal on the transition curve i in fig . [ fig:2 ] .    for the `` short '' blocks ( fig . [ fig:1 ] ) , as well as for many random sequences ( such as in fig . [ fig:3 ] ) , the phase diagram is much more complicated . starting from some value of @xmath20 and for intermediate values of @xmath14 there appear additional solutions corresponding to local minima of the free energy . the broad region where this could take place is bounded by the curves i and ii  in figs . [ fig:1],[fig:3 ] . with increasing @xmath20 the number of such solutions grows quickly . significantly , in a region of the phase diagram some of these become the main free energy minimum . as the number of such solutions grows roughly exponentially with the chain length , we do not attempt to draw all their boundaries of ( meta)stability . instead , we shall designate them collectively as _ frustrated _ phases , explaining this terminology below . an important point here is that , as our analysis shows , these solutions become dominant in a narrow region of the phase diagram due to a subtle competition between the mean energy and the entropy . the mps globule is entropically unfavourable there because the overall shrinking force is insufficiently strong . the values of @xmath89 and @xmath88 are intermediate for these solutions and lie between those of the coil and mps globule ( see figs . [ fig:8],[fig:9 ] ) . in this sense , we can call them nonfully compacted and misfolded states . in comparison , the mps globule ( see fig .  [ fig : dmneq]d ) has a more compact size and better optimised volume interactions , what is manifested in a higher value of @xmath88 . most interesting is the local structure of these additional phases . let us discuss the particular example of the 30(ab ) sequence . the formalism of ref . @xcite has used heavily the assumption of certain symmetries for the mean squared distances @xmath90 due to which these variables can take only @xmath91 independent values , where @xmath72 is the number of blocks . these symmetries are : the block translation invariance , @xmath92 where @xmath71 is the block length , and the more complicated inversion symmetry discussed in detail in ref . these symmetries have a simple meaning  a renumbering of monomers does change the average properties over the statistical ensemble  the interactions of a ring chain remain the same . however , the maximal possible number of dynamical variables in the gsc method is much larger and equals to @xmath93 . surprisingly , it turns out that the frustrated phases are characterised by spontaneous breaking of these symmetries , and therefore only the current version of the method that takes into account all degrees of freedom can describe them . thus , the property ( [ eq : sym ] ) is no longer valid . this phenomenon describes formation of local frustrated heterogeneities ( see figs .  [ fig : dmneq]b , c ) in which pieces of the chain form mps clusters @xcite that are prevented from further coalescing by their hydrophilic shells and high entropic barriers . the role of spontaneous symmetry breaking is well recognised in equilibrium statistical mechanics . what is striking here is that the number of distinct spontaneously broken states becomes huge for large system sizes . this diversity and a special foliating structure of various branches leads in the thermodynamic limit to what is known as a spin glass like frozen phase @xcite of random copolymers ( see e.g. refs . @xcite and numerous references therein ) . we shall return to the issue of spontaneous symmetry breaking in kinetics in subsection [ subsec : kin ] .    even though the kinematic symmetries are not present from the outset for arbitrary , or random , sequences , the structure of the phase diagram ( fig . [ fig:3 ] ) and behaviour of main observables ( figs . [ fig:8],[fig:9 ] ) remain very similar . it is the particular structure , number and the shape of boundaries of frustrated phases that are very sensitive to the sequence . the symmetry that may be broken in this case has a subtler , dynamic , meaning and may be expressed in terms of the replica formalism @xcite . in a sense , for a very long periodic chain , blocks may be viewed as identical copies of a smaller block  length chain . thus , it is by no means surprising that the replica symmetry breaking in our case for periodic systems takes such an explicit manifestation in the breaking of the block translational symmetry . this important point was completely missed , nor could it be discovered , in our considerations of ref . therefore we have achieved here a new significant insight into the problem by a simple extension of the gsc method . an interesting feature of our phase diagrams is that the region between spinodals i and ii  , designating where the frustrated phases can exist , expands dramatically with increasing chain length . for example , at @xmath94 , for the 15(ab ) sequence it lies approximately between @xmath95 , while for the 30(ab )  between @xmath96 , and the curve ii  goes nearly vertically downwards for larger @xmath84 . according to the above interpretation , for infinitely long chain somewhere in this broad region and close to its boundaries there are the actual glass transition curves , which should be determined using proper glass order parameters . thus , for short chains the curves i and ii  may be viewed as approximate indicators of the freezing transitions . the former distinguishes between the homopolymer  like and the random coil , while the latter  between the homopolymer  like ( liquid ) and the frozen globules . as for the region of stable mps globule , it gets relatively smaller for larger systems . that is not surprising  the frustrated phases expand and the phase separation involving larger spatial scales requires stronger interactions . having understood the identifications for the spinodals , we now can recognise in figs . [ fig:1],[fig:3 ] the main features of the phase diagram , though rather distorted , of random copolymer model presented in ref . @xcite . finally , let us comment on the phase diagram of the `` long '' blocks ( fig . [ fig:2 ] ) . the micro  phase separation is obviously easier in this case and it dominates for large values of @xmath20 , so that the frustrated phases are suppressed . we found that , qualitatively , in order to form a frustrated phase , in a finite range of @xmath20 values , the number of ( not necessarily identical ) pieces of the chain with competing interactions should be larger than some critical number , in principle , weakly dependent on @xmath84 . here are a few examples conforming to this qualitative criterion : the phase diagram of 10(ab ) behaves roughly as for 10(aaabbb ) , but for 15(ab ) it behaves similarly to 30(ab ) ; while for 15(aabb ) the phase diagram looks like for 15(aaabbb ) at small and moderate values of @xmath20 , it becomes as for 30(ab ) for much larger values of the amphiphilicity . it is reasonable to conjecture therefore that for a large number of ( aaabbb ) blocks , as well as for extremely high values of @xmath20 and just 10 blocks , the frustrated phases may be found again .      here we shall consider the time evolution of the conformational state of the system away from its initial equilibrium after it has been subjected to an instantaneous temperature jump that causes the two  body interaction parameters @xmath14 and @xmath20 to change . we are interested in quenches from the homopolymer coil , where all monomers are equally hydrophilic ( @xmath97 and @xmath98 ) , to the region of parameters corresponding to the mps globular state , so that the ` a ' species became strongly hydrophobic and the ` b ' species remained weakly hydrophilic ( @xmath99 and @xmath100 ) . the temporal behaviour of the mean squared radius of gyration , @xmath101 , the mps parameter , @xmath102 , and the instantaneous free energy , @xmath103 , in kinetics of folding for different sequences is presented in figs .  [ fig:10]-[fig:12 ] . for the homopolymer ( curves ( a ) ) @xmath89 and @xmath36 decrease monotonically to their final equilibrium values , while the mps parameter , @xmath88 , remains identically zero for there is no distinction between different monomer species . these curves agree with the earlier results of ref . @xcite and serve for reference purpose here . now let us discuss the curves ( b ) corresponding to the periodic ( ab ) sequence . in the previous subsection we have mentioned that the current formalism yields results consistent with those of ref . @xcite only beyond the parameter region of the frustrated phases , which are characterised by spontaneous breaking of the block symmetry . in kinetics the situation is somewhat similar , but the consistence with the previous simplified formalism is even more limited . really , the mean squared radius of gyration remains close to that of the `` effective '' homopolymer ( curve ( a ) ) during the first kinetic stage ( see eq . ( 74 ) in ref . @xcite ) . interestingly , this is not quite so for more complicated sequences . the mps parameter , @xmath102 , for the ( ab ) copolymer grows slowly in a way similar to the splitting of the fourier modes @xmath104 for large indices @xmath105 ( see fig .  4 in ref . @xcite ) . however , contrary to the strange conclusion in ref . @xcite that the kinetics of the ( ab ) copolymer proceeds faster than for the homopolymer , we now have a , natural from the point of view of monte carlo simulations @xcite , slowing down of copolymer kinetics . analysis of @xmath27 shows that this effect is entirely due to the spontaneous breaking of the block symmetry in kinetics , something that has not been accounted for in ref . @xcite . indeed , in fig . [ fig:13 ] we exhibit the time dependence of the mean squared distances between two nearest hydrophobic monomers @xmath106 for the ( ab ) copolymer in kinetics after a quench to the mps globule region . for early , as well as for late , times these functions for different @xmath107 are exactly equal to each other . however , there is well defined intermediate period in time where the block symmetry is spectacularly broken ( see fig . [ fig : dmnkin ] ) . we remark that the symmetry breaks and restores in a step  like manner ( e.g. @xmath108 and @xmath109 join together at @xmath110 , earlier than with other functions ) and also that this effect is relatively strong . we would like to make a general comment here on the nature of spontaneous symmetry breaking in kinetics . thus , normally in such situations at equilibrium there exist a thermodynamically unstable symmetric free energy minimum and a disjoint set of symmetry broken minima , which may be transformed to each other by the residual subgroup of symmetry transformations . these states may also be obtained kinetically as infinite time limits of the time evolution starting from any , for example , the symmetric initial state , which happens to be the main free energy minimum before the quench .    however , the formal structure of the gsc kinetic equations ( [ eq : kinmain ] ) is such that they yield a symmetric solution at any moment in time provided one proceeds from the symmetric initial condition . a question arises then  how can one obtain the kinetics that could lead eventually to the multitude of final states with broken symmetry ? the answer is clear in the exact theory  there are fluctuations that can transform between different spontaneously broken states in kinetics . the gsc method presents , though an optimised and improved , but still a mean field type theory , where such fluctuations are not properly included . manifestation of the kinetic spontaneous symmetry breaking takes a different form there . namely , at some moment in time the symmetric solution of the kinetic equations becomes unstable with respect to perturbations ( whether of the initial condition , or of the interaction matrix , @xmath111 ) . thus , for example , one can add an infinitesimal symmetry breaking term @xmath112 to the two  body interaction matrix and consider the limit of vanishing perturbation in the solution . different choices of the form of @xmath112 in the unperturbed limit would yield all possible types of spontaneously broken kinetic evolution , which are , of course , equivalent to each other . numerically , such a regularisation procedure is not even necessary . there is always an intrinsic perturbation due to computer round off and numerical integration errors . thus , if the symmetry is favourable to be kinetically broken somewhere , numerically one obtains one of the spontaneously broken solutions there , rather than the unstable symmetric solution , unless the symmetry conditions have been imposed by hand . in this situation improvement of the numerical precision would have no profound effect  the kinetic process either does not change , or can change only up to an experimentally unobservable residual symmetry transformation . we note also that such a procedure of perturbing the interaction matrix here is analogous to introducing an external magnetic field in the ising model . the spontaneous magnetisation of a macroscopic sample of ferromagnet may be achieved by gradually switching off the external magnetic field . in the absence of the field there would remain domains with long  range order , but varying directions of the spontaneous magnetisation canceling each other . unfortunately , in our case it is not evident how to experimentally implement an analogous gradual switching off of @xmath112 .    now returning to our results , in fig . [ fig:10 ] one can see that the folding kinetics for the ( ab ) copolymer is about 3 times slower than for the homopolymer of the same length , the effect being even stronger for other sequences we have considered . from fig . [ fig:11 ] and the phase diagrams in figs . [ fig:1]-[fig:3 ] it is evident that the considered quench results in the final state of the mps globule for the `` long '' blocks ( curve ( c ) ) , whereas in the frustrated phases for two other sequences ( curves ( b ) and ( d ) ) . for the latter phases the mps parameter is smaller and the radius is larger than for the mps globule , as we already know from the equilibrium considerations . it is interesting to note that the final relaxation to the frustrated phases may be rather unusual . for example , in figs . [ fig:10],[fig:11 ] for the `` random '' sequence @xmath101 increases and @xmath102 decreases during the last kinetic stage , what is the converse to the behaviour at final relaxations in other cases . another unusual observation is that for some sequences the parameter @xmath88 may even become negative during some time in kinetics , something we never observed at equilibrium . this shows that the structure of non  equilibrium conformations can be very complicated . the instantaneous free energy @xmath103 depicted in fig . [ fig:12 ] turns out to be the quantity most sensitive to the conformational structure of the non  equilibrium state . from that figure it strikes that , while the homopolymer folding is a single smooth relaxational process , the folding of heteropolymers proceeds through a multistep acceleration  deceleration process . the flat regions of a staircase  like function correspond to temporary kinetic arrest of the system in transient non  equilibrium ( mostly symmetry broken ) conformations . generally speaking , in the gsc method we deal with the time evolution of a statistical ensemble of various initial conformations . the flat regions appear due to transient trappings of various members of the ensemble in their local shallow energy minima . since such minima are encountered at different moments in time for different members of the ensemble , their influence on the overall time evolution of averaged observables is manifested in a smooth characteristic slowing down .    qualitatively then , summarising the data for various sequences , most of which we have suppressed here , we can say that the number of local minima affects the number of stair steps , whereas the depth of minima determines the lengths of the steps . this interpretation is highly supported by the strong dependence of the stair case structure on the chain length , the sequence and the interaction parameters , for it is known how the frequency and depths of local potential energy minima depend on these factors . on the one hand , it is known that for a long heteropolymer chain the number of local minima grows exponentially with the chain length @xmath84 , and indeed we see that the number of steps in @xmath103 grows just as quickly . on the other hand , increasing of the amphiphilicity , @xmath20 , leads to higher depths of local minima as well as an increase in their number , and really the steps in @xmath103 become longer and more numerous . let us give another example of such connection . for the `` long '' blocks the interactions are less frustrated than , say , for the ( ab ) blocks , thus there are much fewer of local energy minima , but the barrier between the coil and mps globule is a higher . as a result , the kinetic process ( curve ( c ) in fig . [ fig:12 ] ) proceeds through only one rather long kinetically arrested step . it is worthwhile to make a comment here on the notion of kinetic stages introduced in earlier works @xcite . those we associated with typical structures of the conformation and accompanying kinetic laws . we distinguished at least the following kinetic stages at the collapse of the homopolymer : early time necklace formation , middle time coarsening and a number of final relaxational processes . the multistep character of folding observed in this work affects the middle kinetic stage , resulting in its considerable complication and splitting into multiple substages with respective complex kinetic laws that are determined by the sequence . universality of such kinetic laws is doubtful , but probably it can be recovered by averaging over certain classes of sequences with similar folding properties .    apparently , in the gsc method the kinetics is a motion in the space of @xmath93 averaged dynamic variables , @xmath113 , and it is determined by the profile of the free energy . it may be instructive , using figs . [ fig:10 ] and [ fig:12 ] , to present the kinetics via a parametric plot of @xmath36 vs @xmath86 , the latter being the main , though not the only , relevant `` coordinate '' of the system . that would produce a kind of `` bottleneck '' picture that was much discussed by p.  wolynes and others @xcite in relation to the protein folding problem . this indicates that our method produces behaviour that permits interpretation in terms of phenomenological energy landscape models . finally , let us utilise once again the connection between the kinetic evolution of the free energy and the ruggedness of the potential energy landscape . this would allow us to shed some light on the general structure of the energy landscape for complex heteropolymers . thus , a typical example of the `` random '' sequence kinetics ( curve ( d ) in fig .  [ fig:12 ] ) shows first a fast drop , then appearance of short steps that are growing longer with time , until the last step becomes infinite . this translates to the following equivalent energy landscape versus some collective coordinate : first there is a rapid drop away from the unstable coil state , then the surface flattens and small wrinkles appear on it , they are growing larger in the amplitude gradually becoming high mountains and deep valleys , until there is eventually a very deep ravine corresponding to the `` ground '' state of the system separated by a very high barrier from other minima . in our view , the latter picture bears a remarkable resemblance to a typical mean field energy landscape in spin glass systems @xcite recently discussed by g.  parisi @xcite . this observation seems encouraging to us and indicates that the gsc method is capable of describing very complex systems , though it is too detailed and expensive in the complete form for description on a macroscopic scale . nevertheless , by taking the quenched disorder averaging ( similarly to ref . @xcite ) it is possible to achieve an adequate description for heteropolymers that is alternative to the replica formalism . we believe that the underlying connection between both types of approaches has been clarified to some extent by the current work . in the present work we have extended the gaussian self  consistent method to permit study of heteropolymers with complicated primary sequences . this has been achieved by transforming the kinetic gsc equations , earlier derived for the fourier modes , to the form containing only the mean squared distances between pairs of monomers and relating these equations to the instantaneous gradients of the variational free energy calculated based on the gibbs  bogoliubov principle . the revised gsc formalism possesses important computational advantages , but also it is fundamentally superior to its predecessor in that it can describe phenomena of spontaneous symmetry breaking and formation of structural heterogeneities . we then have applied the extended gsc method to some particular amphiphilic heteropolymers . the equilibrium phase diagrams for these have been obtained in a systematic manner . apart from the coil and two globular states we have discovered that in a wide intermediate region of the phase diagram there may be a large number of frustrated partially misfolded states . despite the fact that such states , the number of which grows exponentially with the chain length , are mostly metastable , some of them become the dominant state in their rather narrow domains . the corresponding potential energy profile of the system has strong resemblance to that of a typical spin glass system . thus , we may conclude that the transition to these states in the thermodynamic limit corresponds to a glassy freezing transition . we note that an observation that for sufficiently long sequences of alternating monomers ( `` short '' blocks ) the micro  phase separated state is displaced by the region of glassy frustrated phases has been known for quite some time and well understood in the framework of the density theories @xcite as a result of the density fluctuations which modify the mean  field behaviour . the block translational symmetry breaking in our approach is another manifestation of the glassy phenomena earlier extensively studied for melts of random heteropolymers . in addition , we can see how the destruction of the micro  phase structures occurs for finite  sized systems and how exactly it depends on a particular sequence . the folding kinetics is found to be strongly affected by the presence of these transient frustrated states along the kinetic pathway . this leads to a complicated kinetic process consisting of multiple steps with pronounced slowing down and then acceleration in the folding rate . it is interesting to note that a typical fairly heterogeneous chain with weakly hydrophilic and strongly hydrophobic units folds kinetically first to one of the misfolded states and can undergo a consequent nucleation process to the main thermodynamic state of the micro  phase separated globule . these results for equilibrium and kinetics confirm many predictions of the earlier treatment adopted for random heteropolymers in ref . @xcite based on a quasi  perturbative averaging of the gsc equations over the quenched disorder . the latter approach describes a large set of random heteropolymer sequences with a gaussian distribution of the monomer amphiphilicity . the kinetics of folding seems to be consistent between the previous and current approaches with respect to main observables such as the radius of gyration , energy and the mps parameter . a smooth slow glass  like relaxation in the former approach is produced by averaging over all different sequences with their particular multistep processes of passage through the frustrated states . although here we observe strong dependence of the phase boundaries in the equilibrium phase diagram on the primary sequence , the major conclusion of refs . @xcite that there are three distinct globular states of the homopolymer  like ( liquid ) , the frozen ( glassy ) and the micro  phase separated ( folded ) globules , and an additional state of the random coil , as well their relative location in the diagram , remain valid . we have previously discussed that the equilibrium limit of the theory in the treatment of ref . @xcite had certain deficiencies due to the additional perturbative approximation , which we attempted to remedy phenomenologically . the current fundamental scheme is free of any such problems and provides a reliable test ground for development of further approximations . we have also discovered that heteropolymer sequences of equal length can be roughly divided into distinct classes possessing similar phase diagrams and kinetic folding properties . for instance , periodic polymers with a few long blocks easily undergo micro  phase separation and for them the frustrated states are suppressed at equilibrium and in kinetics . there is a view in the scientific community that complex random sequences also permit more refined classification . the identification of good folding sequences ( see e.g. ref . @xcite ) is believed to be an important prerequisite for unraveling the mysteries of proteins .    in the current paper the gsc method has been applied to a single chain problem . however , since the kinetic equations ( [ eq : kinmain ] ) are covariant it would be relatively straightforward to apply it to solutions of many polymers . formation of nontrivial mesoscopic globules in copolymer solutions at low concentrations has been predicted in the framework of the method @xcite and recently observed experimentally @xcite . generally , the gsc method may be viewed as extension to the realm of kinetics of the variational treatment , which is a better theory than standard mean  field approaches due to a much larger number of variational parameters . we should admit , however , that in the complete form numerical solution of the gsc equations is computationally expensive for systems of a few hundred of monomers . nevertheless , as we have recently shown @xcite , in some approximation the gsc equations can be reduced to those of a simple mean  field theory such as the flory  huggins one , and some extra corrections to the latter . importantly , the formalism in terms of the mean  squared distances is valid for description of the extended coil as well as of the globular states , whereas the density formalism theories are limited to relatively high densities of the system . the weakest point of the method is the gaussian form of the trial hamiltonian , and thus of the correlation functions , a matter that has been discussed at some length in the past @xcite . however , based on the covariant form of the kinetic equations derived here , we hope , in some not too distant future , to take into account the non  gaussian corrections in an analogy to the treatment of the hard  core repulsion for van der waals systems . the main strength of the revised gsc approach lies in the ability to describe a given heteropolymer sequence of finite length , rather than an _ infinite ensemble _ of sequences characterised in a certain probabilistic way . this is an inevitable preliminary step in developing adequate techniques for theoretical modelling of complex biopolymers . understanding of the relation between the chemical composition ( primary sequence ) and the 3-d equilibrium ( tetriary ) structure , as well as of the kinetics of folding , in proteins is one of the great challenges of the modern biotechnological science . we hope that methods like the one we have presented here will take their right place in the collection of new tools for bioinformatics . the authors acknowledge interesting discussions with professor c.a . angell , professor k.  kawasaki , professor g.  parisi , and our colleagues dr a.v . gorelov and a.  moskalenko . the langevin equation with account of the hydrodynamic interaction includes the oseen tensor @xcite , @xmath114\ ,   \phi^{\alpha'}_{n'},\ ] ] where @xmath115 is the right  hand side of eq . ( [ eq : lang ] ) and the noise has the second momentum proportional to the inverse oseen tensor .    in the preaveraged approximation we have , @xmath116 and the analogs of formulae ( a5,a6 ) of ref . @xcite will be , @xmath117 explicit expressions for a few first terms in the virial expansion of the effective potentials , @xmath118 may be found in appendix [ app : der ] . we introduce the trace as the integration over all monomers coordinates subject to the constraint that the reference point is fixed in the centre - of - mass of the chain , @xmath119 then the partition function is obtained as @xmath120 with @xmath121 given by eq . ( [ eq : h ] ) . the delta  function in ( [ eq : tr ] ) removes the trivial divergence in @xmath122 due to the translational invariance . we choose the trial hamiltonian as an arbitrary linear combination , @xmath123 where we have also introduced arbitrary sources @xmath124 . using the delta  function one can exclude , say , @xmath125 and derive , @xmath126=(2\pi)^{3(n-1)/2}({\rm det}\tilde{k}^{(n-1)})^{-3/2 } \exp\left(\frac{1}{2}\sum_{mm'}{}'({\bf j}_m-{\bf j}_0 )   ( \tilde{k}^{-1})_{mm'}({\bf j}_{m'}-{\bf j}_0 ) \right),\ ] ] where the @xmath47 dimensional matrix @xmath127 is @xmath128 we can calculate the averages by , @xmath129}{\partial { \bf j}_n \partial { \bf j}_{n'}}=(\tilde{k}^{-1})_{nn'},\ ] ] where again the latter identity holds only for @xmath130 . from these quantities it is straightforward to obtain the mean squared distances using eq . ( [ eq : ddef ] ) which also holds only for @xmath130 .    finally , we want to express all independent parameters of the matrix @xmath127 via the quantities @xmath90 . for this we compute the following sums applying the above delta  function constraint , @xmath131 substituting @xmath132 from the second formula here to eq . ( [ rtodd ] ) and recalling the relation ( [ d4 ] ) we derive the desired inverse relation ( [ rtodd ] ) for @xmath133 .    from the partition function ( [ z0 ] ) we obtain the `` entropic '' term @xmath134 $ ] that yields precisely eq . ( [ sr ] ) . the gibbs  bogoliubov variational principle is then based on minimising this variational free energy @xmath135 with respect to the @xmath93 independent variational parameters @xmath90 .      for the homopolymer we have the translational invariance along the chain so that @xmath136 , where @xmath137 . then the matrix ( [ rtodd ] ) may be rewritten as , @xmath138 let us apply a unitary transformation to the fourier coordinates generated by the matrix , @xmath139 by a direct evaluation and recalling the relation for the normal modes , @xmath140 we see that the matrix @xmath141 is diagonal with the eigenvalues @xmath142 and @xmath143 . thus , the logarithm of determinant of the @xmath39-dimensional submatrix @xmath144 yields the standard `` entropy '' of the homopolymer ( see eq . ( 3 ) in ref . @xcite ) up to a trivial constant . two first terms in the virial expansion of the derivatives ( [ eq : v ] ) of the mean energy are explicitly given by , @xmath145 where @xmath146 . it is usually assumed that for negative @xmath14 the stability of the system is ensured by the positiveness of @xmath147 . it turns out however that for the whole range of @xmath148 this does not hold and there are additional pathological solutions with singular free energy . the reason for this deficiency of the model ( [ eq : h ] ) is that we have discarded the terms with two coinciding indices in the three  body contribution . this standard trick is not satisfactory therefore , but fortunately it could be remedied by using another prescription for these terms : @xmath149 in ref . @xcite we show that the addition of the latter term , which is subdominant in the large @xmath84 limit anyway , does not change the earlier results for the homopolymer , and even quantitatively improves the agreement with known results for the dense globular state @xcite . importantly , for heteropolymers this term removes spurious solutions and makes the theory well defined . the corresponding contributions to the mean energy and the effective potentials are , @xmath150    m.  l.  huggins , j. chem . * 9 * , 440 ( 1941 ) ; p.  j.  flory , j. chem . phys . * 9 * , 660 ( 1941 ) ; _ polymer blends _ , edited by d.  r.  paul and s.  neuman ( academic press , 1987 ) ; _ multicomponent polymer systems _ , edited by i  .s .  miles and s.  rostami ( longman scientific and technical , singapore , 1992 ) . d.  j.  meier , j. polym . * c26 * , 81 ( 1969 ) ; d.  j.  meier , j. polym . prepr . * 11 * , 400 ( 1970 ) ; * 15 * , 171 ( 1974 ) ; a.  n.  semenov , sov . jetp * 61 * , 733 ( 1985 ) ; t.  ohta and k.  kawasaki , macromolecules * 19 * , 2621 ( 1986 ) ; l.  leibler , macromolecules * 13 * , 1602 ( 1980 ) ; g.  h.  fredrickson and e.  helfand , j. chem . phys . * 87 * , 697 ( 1987 ) ; t.  hashimoto , macromolecules * 20 * , 465 ( 1987 ) . shakhnovich and a.m.  gutin , j. phys . ( france ) * 50 * 1843 , ( 1989 ) ; s.v . panyukov and s.i . kuchanov , sov . jetp * 72 ( 2 ) * , 368 ( 1991 ) ; a.v .  dobrynin and i.ya . erukhimovich , pisma zh . * 53 ( 11 ) * 545 ( 1991 ) , j. phys . * i 5 * 365 ( 1995 ) . g.  wolynes , in _ spin glass ideas in biology _ , edited by d.  stein ( world scientific , singapore , 1991 ) ; _ protein folding _ , edited by t.  e.  creighton ( wiley , new york , 1992 ) ; _ new developments in theoretical studies of proteins _ , edited by r.  elber ( world scientific , singapore , 1994 ) .                  note that one may use the symmetry properties @xmath151 and @xmath152 only after the differentiation of the energy ( [ eq : e ] ) , which is generally speaking expressed through the four  point correlation functions . one should use the standard differentiation rule , @xmath153 . the structure of the kinetic equations is such that the symmetry properties are automatically satisfied on the solution of these equations at each moment in time .        apparently , for a finite length chain , the collapse transition is characterised by a finite width , and therefore its representation in figs . [ fig:1]-[fig:3 ] is sufficiently schematic with a suitable criterion @xcite adopted to draw a definite curve for it . e.g.  timoshenko , yu.a . kuznetsov , and k.a . dawson , physica * a 240 * 432 ( 1997 ) ; proceedings of the international conference on morphology and kinetics of phase separating complex fluids , ed . by f. mallamace ( messina , italy , 1997 ) . a.  kuznetsov , e.  g.  timoshenko , and k.  a.  dawson , in _ proceedings of the international conference on morphology and kinetics of phase separating complex fluids _ , edited by f.  mallamace ( messina , italy , 1997 ) , to be published in il nuovo cimento .
Abstract:

Answer: in this paper we extend the gaussian self  consistent method to permit study of the equilibrium and kinetics of conformational transitions for heteropolymers with any given primary sequence . the kinetic equations earlier derived by us are transformed to a form containing only the mean squared distances between pairs of monomers . these equations are further expressed in terms of instantaneous gradients of the variational free energy . the method allowed us to study exhaustively the stability and conformational structure of some periodic and random aperiodic sequences . a typical phase diagram of a fairly long amphiphilic heteropolymer chain is found to contain phases of the extended coil , the homogeneous globule , the micro phase separated globule , and a large number of frustrated states , which result in conformational phases of the random coil and the frozen globule . we have also found that for a certain class of sequences the frustrated phases are suppressed . the kinetics of folding from the extended coil to the globule proceeds through non  equilibrium states possessing locally compacted , but partially misfolded and frustrated , structure . this results in a rather complicated multistep kinetic process typical of glassy systems .