single molecule biomechanical experiments allow to measure the forces generated by biomolecules and their response to applied forces  @xcite . for `` passive '' biopolymers , that do not produce active movement against load , the resulting force - extension curves provide a probe of the nonlinear elasticity of the chains and its relationship to configurational changes .
such studies were carried out for single stranded and double stranded dna  @xcite , for the muscle protein titin  @xcite , the extracellular matrix protein tenascin  @xcite , the polysaccharides dextran  @xcite and xanthan  @xcite , as well as the synthetic polymer poly(ethylene - glycol )  @xcite .
theoretical modeling of the elastic behavior of these polymers is often hampered by lack of detailed knowledge concerning the configurations involved and the relevant interactions .
a direct and detailed confrontation between experiment and theory is however possible for the case of homopolypeptides capable of forming an @xmath4-helix . in this case
the molecular configurations are well understood .
furthermore , one can relate the elastic force law to the known zimm - bragg parameters that characterize the helix - coil transition of the free polypeptides .
as we shall discuss , the extension force law of a long polypeptide just above the melting temperature of the helix , @xmath0 , exhibits two plateaus .
both are associated with a coexistence of helical and coil domains .
the low - tension plateau is due to helix formation induced by the chain extension and the accompanying loss of configurational entropy .
the second , high - tension plateau is due to the extension induced melting of the helices when the end - to - end length exceeds the length of the fully helical chain . for @xmath5
only the second plateau survives but the force law for weak extensions exhibits a `` quasi plateau '' _
i.e. , _ a regime with a weak slope that is not associated with a coexistence .
this last regime is due to the facile alignment of the long persistent regions in the helix . in the following
we analyze the force laws , their dependence on temperature and the corresponding population distribution of helical and coil segments . while we are unaware of single molecule experiments on this system , similar behavior was observed experimentally on fibers of collagen during the 1950s  @xcite .
two different treatments of this problem were recently advanced  @xcite . in the present article we present a unified analysis of the problem , tracing the physical origins of the disagreements between the force laws obtained in the two earlier publications and analyzing their range of validity .
our treatment is distinctive in that it is based on free energy argument instead of the customary transfer matrix approach  @xcite .
this free energy argument allows to recover the results of the transfer matrix method but is physically transparent and relatively simple mathematically . within this approach
we consider the extension of a very long homopolypeptide capable of forming @xmath4-helices .
the discussion is limited to quasi - static extension assuming that the rate of equilibration of the configurations of the chain is much faster than the rate of extension .
two models are explored . in one , the helical domains are viewed as rigid rods while in the second , they are modeled as semiflexible chains . for each model we present a rigorous derivation of the corresponding force law allowing for the polydispersity of the helical and coil domains and the associated mixing entropy ( see fig .
[ multi ] ) .
in addition , we explore two useful approximations allowing for simplified calculations of the force law . in one approximation ,
the mixing entropy of the helices and domains is neglected and the plateau is associated with coexistence of one helical domain with one coil domain ( see fig .  [ phases ] ) . within this `` diblock '' or `` @xmath6 '' approximation
the coexistence regimes involve a first order phase transition . in spite of this erroneous conclusion
, this approximation correctly identifies the important length and force scales in the problem .
the second , mean field approximation , neglects the direct coupling of the distribution of domains sizes and the applied tension . as we shall see , this approximation is actually exact when the helical domains are modeled as long semiflexible chains .
the problem and the analysis are of interest from a number of perspectives .
homopolypeptides are the simplest representatives of helicogenic molecules . in this case
it is not necessary to allow for the sequence heterogeneity which affects the elasticity of heteropolypeptides . nor is it necessary to consider the long range interaction that occur in multiple helices . as a result , the number of parameters invoked in the theory is smaller and the mathematical analysis is simpler .
the force laws obtained are determined by known zimm - bragg parameters with no adjustable parameters . accordingly
this system enables a direct comparison between experiment and theory . in this juncture
it is important to note the availability of synthetic methods for producing long polypeptides with a well defined architecture  @xcite .
furthermore , the formalism described provides a convenient basis for the analysis of more complex helicogenic polymers such as dna and collagen . from a polymer science point of view
the elasticity of homopolypeptides is of interest because it reflects the effects of internal degrees of freedom arising from intrachain self - assembly .
the study of these systems enables the exploration of the distinctive non - linear force laws that are the signatures of such intrachain self - assembly .
finally , this elasticity plays a role in the stabilization of the helical state of grafted polypeptides  @xcite . in turn
this is of importance for the function of fusogenic polypeptides and the design of surfaces that favor spreading and growth of cells  @xcite .
the paper is organized as follows . in section
ii we introduce the free energy argument _ via _ an analysis of the helix - coil transition in free , undeformed homopolypeptides .
a semiquantitative analysis of the effect of extension is presented in section iii .
two approaches are used .
one focuses on the intersections of the free energy curves of a pure coil and a pure helix as a function of the end - to - end distance , @xmath2 .
as we shall see , each intersection is a rough diagnostic for a transition between the two configurations and for the occurrence of a plateau in the force law .
the second , more quantitative approach is the diblock or @xmath6 approximation . within this approximation the homopolypeptide
is assumed to consist of two domains , a helix and a coil , and the mixing entropy associated with the polydispersity of the helical and coil domains is set to zero , @xmath6 .
finally , in section iv we present a rigorous analysis of the force laws allowing for the polydispersity of the coexisting helices and coil domains .
as stated earlier , we implement this analysis for two models . in one
the helical domains are viewed as rigid rods and in the second as semiflexible chains described by the worm like chain ( wlc ) model . in this section
we also examine the validity of the mean field approximation introduced in  @xcite . within this approximation ,
the explicit coupling between the tension and the distribution of domain sizes is ignored . as we shall see
, it recovers the exact result when the helices are semiflexible _
i.e. , _ when the persistence length of the helices is finite and the helical domains are large .
the analysis of the rigid helices model recovers the results obtained by tamashiro and pincus  @xcite using a transfer matrix formalism .
a comparison between the two models and an outline of some open problems are presented in the discussion .
it is helpful to first summarize the main features of the helix - coil transition of free polypeptides in the absence of imposed extension .
we will consider only homopolypeptides , consisting of a single type of amino acid residues , in order to avoid complications introduced by sequence of different residues found in heteropolypeptides .
our discussion focuses on the case of long chains where the degree of polymerization of the polymer @xmath7 is large , @xmath8 . to this end
we first discuss the zimm - brag parameters and then introduce the free energy analysis of the transition . as we shall see , this is equivalent to the transfer matrix formalism  @xcite that is traditionally used to analyze the transition  @xcite .
our discussion focuses on helicogenic polypeptides that are capable of assuming two configurations : a random coil and an @xmath4-helix . in the coil state
the monomers are free to rotate thus leading to a small kuhn length with typical value of @xmath9 corresponding to the length of few monomers . for simplicity we will identify the kuhn length with the span of a single monomer , @xmath10 .
the helical configuration involves hydrogen bonds between residues @xmath11 and @xmath12 .
this constrains the orientation of the three intermediate monomers . as a result
the persistence length @xmath13 of the helical configuration is large , of order @xmath14 . at the same time
the projected length of a monomer along the axis of the helix decreases to @xmath15 .
the length of the polypeptide in a helical configuration is thus smaller than the span of a fully extended coil by a factor of @xmath16  @xcite .
the homopolypeptides are modeled as a two state system .
each monomer can exist either in a helical or a coil state .
choosing the coil state as a reference state , the excess free energy @xmath17 of a monomer within a helical domain reflects two contributions .
first is the change in enthalpy upon formation of an intrachain @xmath18-bond , @xmath19 .
the second is associated with the loss of configurational entropy , @xmath20 , in the helical state .
altogether @xmath21 the enthalpy term favors the helical state while the entropy term favors the coil state .
thus , at low temperatures the free energy @xmath17 is negative and the helical configuration is preferred while at high temperatures @xmath17 is positive and the coil state is dominant .
the two terms are comparable at @xmath22 when @xmath23 .
@xmath0 provides an approximate value for the melting temperature of the @xmath4-helix .
an additional parameter is needed to describe the helix - coil transition .
the number of @xmath18-bonds in a helical domain consisting of @xmath24 monomers is @xmath25 .
it is thus necessary to allow for the special state of the terminal monomers of the helical domain .
these two monomers lose their configurational entropy with no gain due to the formation of @xmath18-bonds .
each of the terminal bonds is thus assigned an excess free energy of @xmath26 with respect to the coil state . in the following
we will use @xmath27 measured with respect to the helical state , that is @xmath28 as we shall discuss later , this definition is consistent with the known experimental results . clearly
, this definition is consistent with the traditional one when the coil state is used as a reference .
@xmath27 plays the role of an interfacial free energy associated with the boundary between helix and coil domains .
this term makes the helix - coil transition cooperative since it favors the formation of large domains .
it is customary to describe the helix - coil transition in terms of the zimm - bragg parameters  @xcite : @xmath29 where @xmath30 is the boltzmann constant .
the @xmath31 dependent @xmath32 represents the boltzmann factor associated with adding one monomer to a helical domain .
@xmath33 at the transition temperature @xmath0 while @xmath34 for @xmath5 and @xmath35 for @xmath36 .
@xmath37 corresponds to the boltzmann factor associated with the creation of a helical domain .
the definition of @xmath27 as @xmath38 suggests that @xmath37 is independent of @xmath31 .
the definition of @xmath27 as @xmath39 does not predict that @xmath37 is independent of @xmath31 but it yields an identical estimate for the numerical value of @xmath37  @xcite .
the experimentally measured @xmath37 , as obtained from plots of the helical content _ vs. _
@xmath31 , are in the range of @xmath40 , depending on the residue , and exhibit a weak @xmath31 dependence . as stated earlier , in the following we will use @xmath41 . as we shall discuss , a plot of the helical content _ vs. _
@xmath32 exhibits a sigmoid behavior and the width of the transition scales with @xmath42 .
thus , the transition becomes sharper , more cooperative , as @xmath37 decreases .
having defined the relevant molecular parameters , @xmath17 and @xmath43 or equivalently @xmath32 and @xmath44 , we are in a position to analyze the helix - coil transition as it occurs in free chains upon changing the temperature @xmath31 .
the minimal description of the transition requires the specification of two properties , the total number of monomers in a helical configuration , @xmath45 and the number of helical domains , @xmath46 . while @xmath47 and @xmath48 are sufficient for a discussion of the helix - coil transition of free chains , the analysis of the effect of extension requires additional information .
in particular , the probability distribution of helical domains , @xmath49 , and of coil domains , @xmath50 , comprising of @xmath24 monomers . in order to determine these parameters we minimize the appropriate phenomenological free energy .
this method recovers the results obtained by the transfer matrix approach and is easily generalized to allow for the effect of applied tension .
the free energy allows for two principal contributions : ( i ) the free energy of the monomers .
all @xmath51 helical monomers , including the terminal monomers of the domains are assigned an excess free energy of @xmath17 .
the @xmath52 terminal monomers are assigned an additional free energy of @xmath53 . as a result
their excess free energy with respect to the coil state is @xmath54 .
( ii ) a mixing entropy term associated with the different possible placements of @xmath51 helical monomers and @xmath52 terminal monomers constituting the domain boundaries .
since the helical and coil domains alternate , there is no entropy due to their mixing .
rather , the mixing entropy arises because of the polydispersity of the domains _ i.e. _ , helical and coil domains of different size @xmath24 are considered distinguishable . since there are @xmath46 helical and @xmath46 coil domains the mixing entropy is @xmath55 the free energy per monomer in the unperturbed chain is thus @xmath56 the lagrange multipliers , @xmath57 and @xmath58 , ensure the normalization @xmath59 and @xmath60 .
the two remaining lagrange multipliers , @xmath61 and @xmath62 , impose the average sizes @xmath63 and @xmath64 of the helical and coil domains .
this description corresponds to a simplified version of the transfer matrix method involving a two by two matrix . in both cases
there is no constraint ensuring that helical domains incorporate at least three monomers  @xcite .
since we focus on the case of @xmath65 and @xmath8 the weight of small domains is negligible and this description yields the correct results .
the probabilities @xmath49 and @xmath50 are determined by minimization of @xmath66 with respect to @xmath49 and @xmath50 ( appendix a ) leading to    @xmath67    note that @xmath49 and @xmath50 decay exponentially with decay constants @xmath68 and @xmath69 $ ] .
thus , @xmath49 and @xmath50 are not sharply peaked at the mean values @xmath63 and @xmath70 .
in particular , there is a significant population of small helices even when the average size of the helical domain @xmath63 is large . in this limit , @xmath71 .
similarly , when @xmath70 is large @xmath72 .    upon substituting @xmath49 and @xmath50 in @xmath73 as given by ( [ eq_free2 ] ) , we obtain @xmath66 as a function of @xmath47 and @xmath48 only  @xcite :    @xmath74    the equilibrium conditions @xmath75 and @xmath76 then lead to    @xmath77    in turn , these yield ( appendix b ) the expression for the equilibrium values of @xmath78 and @xmath48    @xmath79^{-1/2}}{s+1 + \sqrt{(s-1)^{2}+4\,\sigma s}}.\end{aligned}\ ] ]    these expressions are identical to the ones obtained _ via _ the transfer matrix method  @xcite .    the behavior of @xmath47 and @xmath48 as a function of @xmath32 and @xmath37 is plotted in fig .  [ theta0 ] . when @xmath80 we recover the boltzmann distribution for independent two - level systems : @xmath81 .
in marked contrast when @xmath82 , the fraction of helical monomers is steplike : @xmath83 for @xmath34 and @xmath84 for @xmath35 .
the chain is either a pure helix or a pure coil . for intermediate values of @xmath37 ,
the plot of @xmath47 _ vs. _ @xmath32 yields a sigmoid curve with a transition region of width @xmath85 centered around @xmath86 .
the number of domains is maximal in the vicinity of @xmath33 and it decreases with @xmath37 .
the parameter @xmath87 serves as a measure of the temperature , @xmath31 . for low temperatures , @xmath88 and
the polymer is mainly in a helical state . at high temperatures , @xmath89 and
the polymer is mostly in a coil state .
the crossover occurs at @xmath90 for which @xmath86 .
near this crossover , @xmath91 and the width of the crossover regime is @xmath92 .
the most notable feature of the force laws of homopolypeptides are the plateaus . in this section
we present qualitative and semiquantitative analysis of the plateaus .
each plateau in the force law corresponds to a coexistence between two states , a helix and a coil .
the situation is somewhat analogous to a liquid - gas transition which leads to a plateau in the pressure - volume diagram . in our case
, the tension @xmath1 plays the role of the pressure and the end - to - end distance @xmath2 is analogous to the volume .
this analogy is however of limited validity since the homopolypeptides are one - dimensional systems with short range interactions and thus incapable of undergoing a first order phase transition . in such systems
the boundary free energy is independent of the size of the domain . as a result , a first order phase transition involving a coexistence of two domains is prevented by the mixing entropy  @xcite .
accordingly , the average domain sizes @xmath64 and @xmath93 , as obtained in the previous section , are independent of the size of the chain , @xmath7 .
nevertheless , the main features of the system can be recovered if one ignores the role of the mixing entropy . as expected , within this rough approximation the plateaus
are erroneously associated with first order phase transitions . in spite of this deficiency
, this approximation recovers the correct length and force scales but the force laws obtained are wrong in detail .
before we present this approximation we investigate the crossing of the free energy curves associated with the `` pure '' helix and coil states .
this is a rough diagnostic for first order phase transitions . in our case
it provides an insight concerning the physical origin of the plateaus as well as a reasonable estimate for some of the length and force scales involved .
the occurrence of the plateaus is signaled by the crossing of the free energy curves of the pure helix and the pure coil as a function of @xmath2 .
each crossing point signals a coexistence between the two `` phases '' and thus a plateau in the force law .    with the unperturbed coil as a reference state ,
the free energy of the coil is simply its elastic free energy . in this section
, we will utilize the fixed @xmath2 ensemble and we thus denote the elastic free energy by @xmath94 . to allow for the finite extensibility of the chain we use the freely jointed chain ( fjc ) model ( appendix c ) for @xmath94 and
the free energy per monomer is @xmath95 where @xmath96 $ ] , @xmath97 is the reduced end - to - end distance of the polymer and @xmath98 is the reduced force . @xmath99 and @xmath100 are related _ via _ @xmath101 where @xmath102 is the langevin function .
the free energy of the helix reflects two contributions .
first is the excess free energy of a helix , @xmath103 . in the limit of an infinite chain ,
when the contribution of the terminal monomers is negligible , the free energy per monomer is @xmath104 . at this stage
we assume that the persistence length of the helix is infinite and that it may be viewed as a perfectly rigid and unextendable rod .
within this view , extending the helix beyond its natural length is associated with an infinite free energy penalty .
this penalty is important because the helix is shorter than the fully extended coil .
as noted earlier , the contribution of each helical monomer to the overall length of the helix is @xmath105 .
altogether , @xmath104 so long as @xmath106 while in the range @xmath107 it is @xmath108 .
three scenarios are now possible ( fig .
[ crossing ] ) . at low temperatures , @xmath109 and @xmath34 .
accordingly , @xmath110 while @xmath111 . in this regime
the helix is more stable than the coil .
when @xmath107 the free energy of the helix diverges and the coil state becomes preferable .
this curve crossing corresponds to a plateau involving the breakdown of the helix due to the applied tension . above the melting temperature , for @xmath112 and @xmath113
the unstretched polypeptide is in a coil state . upon stretching @xmath114 increases and eventually intersects with @xmath115
this intersection signals a plateau in the force law due to the formation of helical domains . at @xmath116
the free energies @xmath114 and @xmath115 cross again because of the infinite free energy of a helix with @xmath107 .
this second intersection signals a second plateau involving the melting of the helix and the formation of coil domains .
finally , for sufficiently high temperatures , @xmath117 , no intersections occur and the polypeptide remains in a coil state for all @xmath99 values .
the threshold temperature for this regime , @xmath118 , is defined by the equality @xmath119 where @xmath120 is the free energy of a stretched coil with @xmath116 .
the intersection of the free energy curves provides a rough indicator for the occurrence of plateaus in the force law .
however , it does not allow to calculate an explicit force - extension diagram .
a simple approach yielding the main features of the force law is the @xmath6 or diblock approximation  @xcite . within this approximation the force law
is calculated for a diblock copolymer with one helical domain and one coil domain .
the size of the annealed domains varies with the extension and is determined by the equilibrium conditions . as is implied by its name , in this approximation @xmath121
is neglected . as expected ,
setting @xmath6 leads to a first order phase transition associated with a coexistence of two `` phases '' .
consequently , some features of the force law are wrong .
in particular , this approximation yields a perfect plateau with @xmath122 while in reality @xmath1 increases weakly with @xmath2 . these features notwithstanding , this simple approximation does recover the correct length and force scales . within this discussion
we will also briefly examine the effect of the persistence length of the helical domains .
in particular we compare the behavior of a helix endowed with an infinite persistence length to that of a helix with a finite one .    within the @xmath6 or diblock approximation
the free energy @xmath123 of the chain is    @xmath124    here , the first term allows for the excess free energy of the helical monomers . since @xmath6 the helical and coil phases separate into two domains .
the `` interfacial '' free energy of the single helical domain gives rise to the second , constant term .
these two terms are supplemented by an elastic free energy @xmath125 that depends on @xmath47 , @xmath126 and the end - to - end distance , @xmath2 .
different versions of the @xmath6 are possible , depending on the choice of @xmath125 . in the following we focus on the simplest version where @xmath125 is given by the gaussian elasticity of an ideal random coil in the fixed @xmath2 ensemble .
two additional assumptions are invoked : ( i ) the helical domain is perfectly rigid and unextendable and ( ii ) the helical domain is perfectly aligned with the applied force .
these assumptions specify the elastic free energy per monomer @xmath127 this simple version of the diblock approximation permits an explicit analytical solution for the force law .
the corresponding equilibrium condition is @xmath128 within this approximation the equilibrium condition @xmath129 is an identity yielding no additional information .
the solution of this quadratic equation in @xmath130 is @xmath131 .
this defines a critical @xmath17 @xmath132 corresponding to a critical @xmath133 or @xmath118 . while @xmath134 , the equilibrium condition can not be satisfied , that is , the polypeptide remains in a coil state irrespective of the imposed tension .
a critical point occurs for @xmath135 .
when @xmath136 two solutions exist , each defining a plateau in the force _ vs. _ extension curve .
the solution @xmath137 corresponds to the lower plateau , obtained for @xmath138 , where @xmath139 this first regime disappears for @xmath140 or @xmath5 .
the condition @xmath84 defines a lower boundary of the plateau , @xmath141 @xmath142 increases as @xmath17 increases until @xmath143 for @xmath135 .
for @xmath144 the fraction of helices is @xmath84 and the chain is in a purely coil form .
the upper boundary of this plateau occurs at @xmath145 when the chain is fully helical , @xmath83 . at intermediate extensions , @xmath146 the fraction of helical bonds increases as @xmath147 and @xmath148 occurs at the midpoint of the plateau @xmath149 .
the equilibrium free energy of the chain for @xmath144 is @xmath150 while for @xmath151 it is @xmath152 the second plateau , is described by the solution , @xmath153 , leading to @xmath154 the lower boundary of this plateau occurs at @xmath155 when @xmath83 .
the upper boundary , when @xmath84 , occurs at @xmath156 for @xmath157 the fraction of helical monomers decreases as @xmath158 and @xmath148 occurs at the midpoint of the plateau @xmath159 .
notice that the width of the two plateaus shrinks as @xmath160 or @xmath161 .
for @xmath162 the fraction of monomers in a helical state is again @xmath84 and the corresponding free energy is @xmath150 .
the equilibrium free energy of the chain in the range @xmath163 is @xmath164 note that only the lower plateau is obtained if @xmath165 is assumed@xcite .
the @xmath6 approximation allowed us to obtain rough force laws corresponding to the three scenarios identified by the curve crossing argument .
thus , for @xmath134 or @xmath117 , the polypeptide remains in a coil configuration for the whole range of extension .
the `` two plateau scenario '' occurs for @xmath136 corresponding to @xmath166 or @xmath167 .
in this case we can distinguish five regimes ( fig .
[ phases ] ) : ( i ) for low extensions the chain is in a coil state , @xmath84 , and the elastic behavior is gaussian leading to @xmath168 .
( ii ) in the range @xmath146 helix and coil domains coexist according to the lever rule @xmath169 . this coexistence is associated with a plateau with a constant tension @xmath170 ( iii ) at @xmath171 the force curve exhibits a discontinuity associated with a step like increase to a second plateau .
the magnitude of the step , @xmath172 is @xmath173 ( iv ) in the range @xmath157 helix and coil domains coexist according to the lever rule @xmath174 .
this coexistence is associated with a second plateau with a constant tension @xmath175 ( v ) finally , for @xmath176 the chain is in a strongly stretched coil state obeying @xmath168 .    when @xmath177 or @xmath178 the unperturbed chain is in a helical state .
consequently , the lower plateau disappears .
it is replaced by a `` quasi plateau '' which is not associated with a helix - coil coexistence .
rather , it is due to the facile orientation of the helical chain which is modeled as a rigid rod . at @xmath116
a jump occurs to the upper plateau .
the analysis presented above utilized the gaussian elastic free energy in the fixed @xmath2 ensemble to describe the coil .
while this allows us to obtain explicit expressions for all the characteristics of the force curve , it also introduces errors .
clearly , this choice of @xmath125 does not capture the behavior of the coil in the high extension regime , where finite extensibility begins to play a role . as a result ,
the values of @xmath179 are overestimated .
this has two consequences .
first , within this approximation the span of the two plateaus are equal @xmath180 while for more realistic choice of @xmath125 the higher plateau is narrower . a second more important
pathology is the disappearance of regime ( v ) in the @xmath181 case .    to allow for the finite extensibility it is necessary to use the freely jointed chain model ( fjc ) for the coil segment  @xcite .
it is important to note that so long as @xmath182 and @xmath183 the discrepancies between the two schemes are rather small ( fig .  [ force_db ] ) .
another ingredient of our simplified analysis is the modeling of the helices as rigid and unextendable rods .
an alternative description of the elastic properties of helices models them as worm like chains _ i.e. , _ semiflexible chains capable of undergoing undulations ( see appendix d )  @xcite .
as we shall see later , this semiflexible approximation of the helices has a qualitative effect on the results of the rigorous analysis of the system .
however , in the context of the @xmath6 approximation for @xmath183 , it mainly affects the transition between the two plateaus ( fig .
[ force_db ] ) . when the helix is modeled as a rigid rod the transition occurs as a jump that takes place at @xmath116 . on the other hand , when the wlc model is used the jump , while abrupt , takes place over a finite interval _
@xmath184 is replaced by @xmath185 . for the @xmath140 ( fig .
[ force_db2 ] ) the transition at the upper boundary of the `` quasi plateau '' is more gradual when the elasticity of the semiflexible helix is allowed for .
finally , note that the diblock approximation is recovered from the more rigorous treatments presented below when the @xmath186 limit is taken before the @xmath187 limit .
when @xmath188 limit is approached for a finite chain it leads to a single domain .
the analysis presented in the previous section fails to allow for a number of important ingredients .
the most important is the mixing free energy , @xmath189 , associated with the polydispersity of the domains . in this section
we present an analysis that allows for the role of this mixing entropy . because of @xmath121 the helix - coil transitions no longer take place as first order phase transitions . in turn
this has two important consequences : ( i ) the `` perfect '' plateaus , with constant forces @xmath190 and @xmath191 are replaced by `` imperfect '' plateaus with a finite slope .
( ii ) the crossovers between the plateaus and the neighboring regimes become smoother .
the analysis is based on a generalization of the approach described in section ii that is , all chain characteristics , @xmath49 , @xmath50 , @xmath78 and @xmath48 are obtained from the minimization of the appropriate free energy _
i.e. _ , the equilibrium conditions . as in section iii
, the free energy is of the form    @xmath192    where , @xmath66 is the free energy of the unperturbed free chain and @xmath125 , accounts for the elasticity of the chain .
however , in the following we utilize different expressions for these two terms .
since we aim to obtain @xmath49 and @xmath50 it is necessary to utilize @xmath66 as given by equation ( [ eq_free2 ] ) .
the handling of the elastic free energy is significantly different . in section iii
we used the gaussian elasticity to describe the deformation behavior of the coils and assumed that the helices contribute to the elastic behavior as a geometric constraint . while this rough description is sufficient at the level of the @xmath6 approximation , a more refined treatment is needed for the detailed analysis to be undertaken below .
in particular , it is necessary to allow for two additional features : ( i ) the finite extensibility of the chain and ( ii ) the gradual alignment of the helical domains with the applied tension . the freely jointed chain ( fjc ) model allows us to incorporate the finite extensibility of the coil domains . with regard to the helices
we will explore two different descriptions that impose finite extensibility on the helices .
first we will model the helical domains as rigid and unextendable rods and allow for their length dependent alignment by using the fjc model _
i.e. , _ the helical domains are considered as freely jointed rigid rods forming a chain .
this amounts to the assumption of an infinite persistence length for each helical domain .
this picture is reasonable so long as the size of the helical domains is small compared to the persistence length , @xmath13 . in the second approach ,
each of the helical domains is viewed as worm like chain ( wlc ) that is , semiflexible chain characterized by a finite persistence length and capable of undergoing undulations .
it is assumed that the domains are large in comparison to @xmath13 .
an important difference between the fjc and the wlc models is that in the later the persistence length depends on the applied tension .
our approach utilizes the analysis of marko and siggia of the extension elasticity of the wlc  @xcite .
the utilization of these models makes the constant @xmath1 ensemble mathematically convenient .
this choice of ensemble is also preferable because it corresponds better to the experimental situation encountered in single molecule measurements . in each case , we first obtain the equilibrium values of @xmath49 , @xmath50 , @xmath47 and @xmath48 of a chain subject to a given tension , @xmath1 , by minimizing @xmath193 with respect to @xmath49 , @xmath50 , @xmath47 and @xmath48 .
this allows us to determine the equilibrium @xmath193 thus enabling us to calculate the @xmath3 diagram by using @xmath194      when the helical domains are viewed as rigid rods and the fjc is invoked , the reduced end - to - end distance @xmath97 of the polypeptide is @xmath195    where @xmath98 is the reduced tension .
the first part of this expression corresponds to the projection of the coil monomers along the direction of the tension .
the second part reflects to the contribution of the helical domains to @xmath99 .
because the individual monomers within the coil domains are aligned independently , the domain size does not play a role .
the extension of each coil domain is extensive with respect to the number of monomers , @xmath24 .
consequently , this contribution may be expressed in terms of the first moment of the probability distribution @xmath196 . as a result
the full probability @xmath50 is irrelevant . in marked contrast , the monomers within the helical domains are not aligned individually but as part of the rod - like domains .
the non - linearity introduced by the @xmath197 term in the expression for @xmath99 couples @xmath1 and @xmath49 . as a consequence , @xmath49 will be explicitly modified by the applied tension while @xmath50 depends on @xmath1 only implicitly , _ via _ @xmath47 and @xmath48 .
eliminating @xmath50 from ( [ r_rigid ] ) , we obtain @xmath198 this expression for @xmath199 enables us to obtain the corresponding elastic free energy @xmath200 ( appendix c ) @xmath201 where @xmath96 $ ] .
altogether the free energy per monomer is thus @xmath202 since the @xmath1 dependence of @xmath123 is due only to @xmath125 , ( [ rdef ] ) yields @xmath2 as given by ( [ rdef2 ] ) . as before ,
@xmath49 and @xmath50 are determined by minimizing @xmath193 with respect to these probabilities subject to two constraints : ( i ) the normalization of the probabilities @xmath49 and @xmath50 and ( ii ) the average sizes @xmath64 and @xmath93 of the coil and helical domains . upon substituting the resulting @xmath49 and @xmath50 in @xmath123 we obtain @xmath193 as a function of @xmath47 and @xmath48 .
@xmath78 and @xmath48 are then determined by the equilibrium conditions @xmath203 .
the details of the calculation involved are described in appendix e. @xmath47 and @xmath48 for a given tension @xmath1 are    @xmath204    where    @xmath205
the probabilities @xmath50 and @xmath49 are    @xmath206    @xmath49 , as @xmath50 , is essentially an exponentially decaying function .
however , the characteristic decay constants are now dependent on @xmath1 .
the dependence is explicit for @xmath49 while for @xmath50 the @xmath1 dependence is implicit , _ via _ @xmath47 and @xmath48 .    in fig .
[ theta2 ] , the equilibrium @xmath47 and @xmath48 , are plotted as functions of @xmath98 for different values of @xmath37 .
clearly , the results of the @xmath6 or diblock approximation are approached in the limit of strong cooperativity , @xmath188 .
in particular , @xmath48 approaches zero while @xmath47 exhibits a step like behavior . thus , @xmath207 outside the plateaus range , @xmath208 and @xmath209 , while within the plateaus range , @xmath210 we find @xmath211 and @xmath212 .
consequently , the mean size of the helical domains , @xmath93 , becomes very large compared to @xmath13 . in turn
, this suggests that the validity of the model of helical domains as rigid rods is questionable . clearly , the rigid helices model is meaningful only while the persistence length is much larger than @xmath213 .
once the equilibrium characteristics of the chain for a given @xmath1 are found , it is possible to obtain the force - extension diagram . the reduced end - to - end distance , @xmath99 , as a function of the reduced tension , @xmath100 , is ( appendix e )    @xmath214 .\ ] ]    the corresponding force - extension diagram is plotted in fig
.  [ force ] .
again , the results of the @xmath6 or diblock approximation are approached in the @xmath188 limit .
note the `` jump '' in the force at @xmath215  @xcite . at this point , between the two plateaus ,
the polypeptide is in a `` purely '' helical state , @xmath211 and @xmath212 .
this jump is another indication of difficulties due to the approximation of the helices as rigid rods . as we shall see in the next section
, this unphysical feature disappears when the finite persistence length of the helices is allowed for .
it is of interest to compare the exact solution of the model to the mean field approximation ( mf )  @xcite .
the mf approximation involves the removal of the explicit coupling between @xmath50 , @xmath49 and @xmath1 . as noted earlier , this decoupling is exact for  @xmath50 . for @xmath49
the decoupling is achieved by replacing @xmath216 in ( [ rdef2 ] ) by @xmath217 . in effect , all helical domains are assigned an equal size of @xmath93 and the effect of their polydispersity on @xmath125 is ignored . the corresponding free energy is    @xmath218    the equilibrium conditions @xmath219 @xmath220 lead to    @xmath221 , \\ & & \frac{y^{2}}{(\theta -y)(1-\theta -y ) } = \\ \nonumber & & \hskip 0.4 cm   = \frac{\sigma y}{\theta \gamma x } \sinh \left
( \frac{\theta \gamma x}{y}\right ) \exp \left [ -\frac{\theta \gamma x}{y}{\cal l}\left ( \frac{\theta \gamma x}{y}\right ) \right ] .\end{aligned}\ ] ]    the overall features of the mf solution are rather similar to the exact solution ( fig .  [ force ] ) .
the two primary differences are : ( i ) the mf solution indicates that the transition between the plateaus involves a second order phase transition  @xcite in marked distinction to the exact solution .
( ii ) because of ( i ) the plateaus do not exhibit an inflection point and the characteristic value of the slope , at @xmath116 , scales as @xmath222@xcite .
the assumption of perfectly rigid helical sequences implies an infinite persistence length .
this assumption is reasonable so long as the size of the helical domains is small compared to the known values of their persistence length , @xmath13 : a typical value is @xmath223 corresponding to @xmath224 monomers .
as we have seen , the rigid rod approximation for an infinite chain leads to helical domains that are larger . in turn
, this leads to an unphysical jump , located at @xmath225 , in the force - extension curve .
two approaches allow to circumvent this difficulty .
one approach treats each helical monomer as an extendable spring _
i.e. _ , it introduces an elastic penalty associated with the deviation of the length of the helical monomer from its equilibrium length @xmath226 .
this approach was utilized in the analysis of the extension behavior of dna  @xcite and of poly-(ethylene glycol )  @xcite . in the following
we adopt a different approach that proved successful in fitting the force curves of dna . within this approach
the helical domains are considered as worm like chains ( wlc ) instead of rigid rods ( appendix d ) .    in the wlc model
the helical domains are viewed as semiflexible chains undergoing undulations .
the angular correlations along the chain decay because of the thermal undulations . within this model
the decay length of the angular correlations , @xmath227 depends on the tension , @xmath228 for weakly stretched chains @xmath229 is identical to the persistence length , @xmath13 , of the unperturbed chain , @xmath230 while for stronger tensions @xmath231 and the correlations decay faster . a long helical domain of length @xmath232
may thus be considered as a fjc of @xmath233 `` rigid '' segments of length @xmath234 and the reduced end - to - end distance is @xmath235 where @xmath102 is the langevin function and @xmath236  @xcite .
as before , the expression for @xmath99 reflects two contributions .
the first term corresponds to the alignment of the coil monomers with the applied tension using the fjc model .
the second term describes the alignment of the semiflexible helical domains . in the following
we assume that the helical domains are larger than @xmath237 that is , @xmath238 . the error introduced by this assumption diminishes as @xmath1 increases and @xmath237 decreases .
within this approximation each @xmath234 segment in a helical domain is aligned independently .
this is in marked contrast to the rigid helices model where all the monomers within a helical domain align as a unit . as a result we may utilize @xmath239 in order to obtain @xmath240 .
consequently , @xmath99 is independent of the probabilities @xmath50 and @xmath49 .
expression ( [ dist_flex ] ) for @xmath199 enables us to obtain ( appendix d ) the corresponding elastic free energy @xmath200 @xmath241 in marked distinction to the `` rigid helices '' scenario , this elastic free energy is independent of the probabilities @xmath49 and @xmath50 .
accordingly , @xmath49 and @xmath50 are not modified explicitly by the extension and obey ( [ proba_h ] ) and ( [ proba_c ] ) .
in other words , the mean field approximation provides an exact solution of the problem .
altogether , the free energy of a polypeptide , using the wlc model for the helical domains , is @xmath242{\cal l}(\alpha ) .
\nonumber\end{aligned}\ ] ] this expression is similar to @xmath123 for a free , undeformed chain ( [ eq_free1 ] ) .
the two free energies differ in two respects : ( i ) ( [ eq_free_flex ] ) contains an additive term @xmath243 which is independent of @xmath47 and @xmath48 .
( ii ) the force independent @xmath32 in ( [ eq_free1 ] ) is replaced by a force dependent @xmath244 @xmath245 .
\label{stilde}\ ] ] consequently , the equilibrium conditions @xmath203 lead to expressions of the same form as those found for the free chains but with @xmath246 replacing @xmath32 @xmath247 ^{-1/2}}{\tilde{s}+1+\sqrt{(\tilde{s}-1)^{2}+4\,\sigma \tilde{s}}}.\end{aligned}\ ] ] substitution of the equilibrium values of @xmath47 and @xmath48 for a given @xmath100 in ( [ dist_flex ] ) yields the corresponding equilibrium end - to - end distance .
plots of @xmath47 _ vs. _ @xmath248 and @xmath248 _ vs. _ @xmath249 , as obtained for the `` semiflexible helices '' scenario for different @xmath37 , are depicted in fig .
[ thetawlc ] .
the replacement of the `` rigid helices '' by semiflexible ones leads to the disappearance of the `` jump '' in the force law .
another difference with respect to fig .
[ theta2]a is that @xmath83 is never attained .
again , we recover the appropriate diblock approximation in the limit of @xmath188 .
all the force - extension curves in fig .
[ thetawlc]b cross at three different points .
the external points correspond to a reduced force @xmath250 satisfying @xmath251 . in turn , this condition is equivalent to the equilibrium condition @xmath252 for @xmath123 corresponding to the appropriate diblock approximation in which the coil is modeled as fjc and the helix as a wlc  @xcite . in particular , @xmath250 correspond to points in the plateaus where @xmath148 , @xmath253 /2 $ ] and @xmath254 .
the inner crossing point corresponds to @xmath255 for which @xmath244 is maximal .
@xmath255 satisfies the condition @xmath256 $ ] _ i.e. _ , to an equality between the projected length of the coil and the helical monomers along the force . from ( [ dist_flex ] ) and ( [ stilde ] ) ,
it is apparent that this condition is independent of @xmath37 .
the plateaus depicted in fig .  [ thetawlc]b exhibit a small but finite slope .
its characteristic value , at the inflection points @xmath250 , scales with @xmath257  @xcite .
our analysis of the extension behavior of polypeptides utilized two models .
the two differ in the description of the helical domains . in one model , the helical domains are viewed as perfectly rigid and unextendable rods .
the `` rigid helices '' model is reasonable so long as the length of the helical domains is small compared to the persistence length . in the second approach
the helical domains are treated as long semiflexible polymers within the worm like chain ( wlc ) model .
this approach is preferable when the length of the helical domains is larger than the force dependent decay length @xmath258 . for each of the models
, we presented a rigorous solution and two approximations . within the diblock approximation one
neglects the mixing entropy associated with the polydispersity of the domains . in the mean field approximation
there is no direct coupling between the distribution of sizes of the domains and the applied tension . in all cases considered ,
our analysis is based on the minimization of a phenomenological free energy .
this approach is equivalent to the transfer matrix method but it is physically transparent and mathematically simple .
the general features of the extension behavior of the polypeptides are independent of the model and the approximation scheme .
three different scenarios are possible for the extension of long polypeptides : ( a ) for @xmath117 the polypeptide remains in a coil configuration for all extensions .
( b ) in the range @xmath167 the chain is initially in a coil state but it undergoes a helix - coil transition upon extension .
the helical domains eventually melt upon stronger extension .
consequently , the force law exhibits two plateaus associated with a helix - coil coexistence .
( c ) for @xmath5 the unstretched polypeptide is a helix which melts upon stretching thus leading to a force law with a single plateau .    in the second scenario , @xmath167 , it is possible to distinguish between five regimes : ( i ) for weak extensions the polypeptide behaves as a random coil exhibiting a linear response .
( ii ) stronger extensions , and the associated loss of configurational entropy , favor the formation of helical domains .
the associated coexistence of helical and coil domains gives rise to a lower plateau .
( iii ) a steep crossover corresponding to the stretching of a stiff , mainly helical , polymer leads to ( iv ) a second , higher plateau due to the break - up of the helices when the imposed end - to - end distance is larger then the length of a fully helical chain .
this plateau is also associated with helix - coil coexistence .
( v ) when the helical content is reduced to zero the chain exhibits a strong extension behavior of a random coil , with deviations from the linear response behavior due to finite extensibility effects .    in the third scenario , @xmath5
, the regimes ( i ) and ( ii ) disappear and regime ( iii ) extends down to @xmath259 .
the small slope of the force law in this regime ( fig .
[ f_s1.1 ] ) is reminiscent of a plateau .
however , this quasi - plateau is not associated with a helix - coil coexistence .
rather , it is due to the easy alignment of the persistent segments in the helix .    the two models , `` rigid '' and `` semiflexible '' helices , differ with regard to the precise details of the force diagram .
thus , region ( iii ) occurs as a `` jump '' within the `` rigid helices '' model but is a steep yet gradual increase in the `` semi - flexible helices '' model .
the slope of the plateaus within the rigid helices model scales as @xmath260 ^{2}$ ] @xcite while in the `` semi - flexible helices '' model it scales as @xmath257 .
overall , the crossovers between the various regimes are sharper in the `` rigid helices '' model .
the diblock approximation , with the appropriate elastic term , allows to recover the main features , length and force scales , of each model .
however , within this @xmath6 approximation the transition takes place as a first order phase transition . as a result
it yields perfect plateaus , @xmath261 or @xmath262 , and sharp crossovers .
the mean field approximation , where the polydispersity is decoupled from the tension is exact for the `` semiflexible helices '' model . in the case
of the `` rigid helices '' model this approximation , while it captures the overall features of the force law also gives rise to a number of artifacts  @xcite .
neither of the models considered captures the physics completely .
the `` rigid helices '' model fails to describe the behavior of helical domain that are long compared to the persistence length . within the `` semiflexible helices ''
model all helical domains are assumed to be long in comparison to @xmath237 .
accordingly , this model is likely to mishandle short helical domains that do behave as rigid rods .
the importance of this deficiency diminishes however as @xmath1 increases and @xmath237 decreases .
a full description of the system requires a model combining the positive features of the two models described above .
while our analysis focused on the case of homopolypeptides forming @xmath263 -helices , our approach can be extended to more complicated systems . in the case of heteropolypeptides
it is necessary to allow for the residue dependence of @xmath17 . for double stranded dna
it is necessary to allow for long range interactions as given by the stockmayer formula for @xmath44  @xcite .
we would like to thank m. n. tamashiro and p. pincus as well as i. rouzina and v. a. bloomfield for providing us with their results prior to publication .
the work of ab was supported by the marie curie fellowship hpmf - ct-1999 - 00328 .
the minimization of the free energy ( [ eq_free2 ] ) with respect to the probabilities @xmath49 and @xmath50 yields    @xmath264    these are solved by @xmath265 and @xmath266 where    @xmath267    the parameters @xmath268 , @xmath269 , @xmath270 and @xmath271 are determined using the constraints    @xmath272    yielding @xmath273 , @xmath274 , @xmath275 and @xmath276 thus leading to the expressions ( [ proba_h ] ) for @xmath49 and ( [ proba_c ] ) for @xmath50 .
the equilibrium conditions @xmath277 with @xmath66 given by ( [ eq_free1 ] ) lead to    @xmath278    from ( [ appendixa2 ] ) , we obtain    @xmath279\ ] ]    or    @xmath280    where @xmath281 . the second condition ( [ appendixa1 ] ) yields    @xmath282    combining these two equations allows us to eliminate @xmath48 and obtain an equation for @xmath78 as a function of @xmath32 and @xmath44 ,    @xmath283    leading to    @xmath284
the freely jointed chain model may be considered as the macromolecular analog of the langevin theory of paramagnetism  @xcite .
a flexible chain comprising @xmath7 monomers is the counterpart of @xmath7 non interacting magnetic dipoles .
the monomer length @xmath285 is the analog of the magnetic moment @xmath286 , the applied tension @xmath287 plays the role of the magnetic field @xmath288 and the induced magnetic moment , @xmath289 is the analog of the end to end distance @xmath290 each monomer is assigned an orientational energy @xmath291 where @xmath292 is the angle between the force and the segment .
the reduced end - to - end distance @xmath249 of the polymer , at a temperature @xmath31 , is the mean value of @xmath293 :    @xmath294    where @xmath98 is the reduced force and @xmath295 is the langevin function .
the corresponding elastic free energy in the fixed @xmath1 ensemble , as given by @xmath296 , leads to the following free energy per monomer : @xmath297 where @xmath298 $ ] , is the integrated langevin function .
the elastic free energy in the fixed @xmath2 ensemble @xmath299 is the legendre transform of @xmath300 where @xmath97 and @xmath100 are related by @xmath301 . for weak extensions
@xmath94 reduces to the familiar gaussian form @xmath302 while for strong extensions it diverges as @xmath303 .
within the wormlike chain model ( wlc ) the polymer is viewed as a bendable rod  @xcite . in the following
we consider the case of a chain of constant contour length @xmath304 and a single bending modulus @xmath305 where @xmath13 is the persistence length of the unperturbed chain .
the chain trajectory is described by the position of a point on the chain , @xmath306 at the contour length @xmath307  the unit tangent vector @xmath308 specifies the local orientation of the chain .
a straight rod corresponds to @xmath309 or @xmath310 accordingly , the bending energy is a quadratic function of @xmath311 . since @xmath312 ,
the only quadratic form is @xmath313 and the conformational distribution of the chain , its partition function , up to a contour length @xmath32 , is @xmath314\ ] ] where @xmath315 is a function of @xmath316 and @xmath317 _ i.e , _
this describes a gaussian process with the constraint @xmath319 . as
such it is analogous to rotational diffusion and obeys @xmath320 where @xmath321 is the rotational operator that is related to the angular momentum operator in quantum mechanics , @xmath322 , as @xmath323 .
the correlations along the unperturbed chain decay as @xmath324 when the chain is subjected to an external field @xmath325 the expression for @xmath326 becomes @xmath327\ ] ] the effect of tension @xmath287 is specified by @xmath328 thus leading to @xmath329 it is possible to find eigenstate satisfying @xmath330 . for a long chain
the free energy is extensive and the eigenfunction expansion of @xmath331 involves terms of the form @xmath332 where @xmath333 is the free energy per unit length due to the bending fluctuations .
an analytical solution of equation ( [ eq_psi ] ) is not available .
marko and siggia  @xcite argued that the wlc free energy is determined by the lowest @xmath334 in the spectrum and utilized a variational approach in order to determine @xmath334 .
the force law is then recovered from @xmath335 . to this end
they chose the trial function @xmath336 $ ] where @xmath4 is a variational parameter and @xmath337 is the angle between @xmath338 and @xmath287 .
as we shall see later , @xmath339 where @xmath340 is a force dependent decay length that replaces @xmath13 in characterizing the decay of @xmath341 .
assuming @xmath342 @xmath343 this leads to  @xcite @xmath344 where @xmath295 is the langevin function and @xmath98 is the reduced force .
the argument @xmath345 that minimizes @xmath334 is specified by @xmath346 which leads to : @xmath347{\cal l}^{\prime } ( \alpha ) .\ ] ] where @xmath348 . here
the end - to - end distance is simply deduced from : @xmath350 this end - to - end distance is identical to the one obtained from the fjc model if the reduced force @xmath351  @xcite is replaced by @xmath352 . for small forces , when @xmath353 and @xmath354 ( [ eq_alpha ] ) leads to @xmath355 or @xmath356 . for high
tension @xmath4 increases more slowly . in this limit @xmath357 while @xmath358 and ( [ eq_alpha ] ) yields @xmath359 or @xmath360 .
the difference between the variational approximation and the exact solution of the wlc model is less than @xmath361 over the whole range of forces . the associated elastic free energy , @xmath362 is @xmath363 since @xmath364 and where @xmath334 corresponds to the explicit expression in ( [ eq_g ] ) with @xmath365 according to ( [ eq_alpha ] ) . in eq .
( [ eq_fel ] ) this elastic free energy is used for a helical domain of length @xmath366 with @xmath24 monomers .
in this appendix we determine the probabilities @xmath367 and @xmath368 of a coil and a helical segment with @xmath24 monomers as well as the fractions of helical monomers @xmath78 and of helices @xmath48 .
the probability @xmath49 is modified by the stretching and is more difficult to obtain .
to do so we consider the two remaining equilibrium conditions @xmath374 and @xmath375 after substitution of the equilibrium form of @xmath50 as given by ( [ prob_pc ] )                                            to obtain the end - to - end distance @xmath2 as a function of the tension @xmath1 we use @xmath393 which leads to @xmath394 the second term may be tackled exactly using the definition @xmath395 : @xmath396                                                                to be precise , this approach is equivalent to a transfer matrix description utilizing a second order transfer matrix . in turn , this description is sufficient for the analysis of the helix - coil transition in an infinite chain .
the final expression for the mixing entropy for the free chain can be calculated directly  @xcite .
@xmath406 where @xmath407 is the product of two factors : @xmath408 the number of possible placements of @xmath46 boundaries to separate @xmath409 helical bonds into @xmath46 helical domains and @xmath410 the number of ways to place @xmath46 boundaries to separate @xmath411 coil bonds into @xmath46 coil domains . using the sterling formula
@xmath412 we obtain @xmath413 $ ] where @xmath414 stands for @xmath415 .        in the fixed @xmath1 ensemble ,
the free energy is @xmath416 where @xmath96 $ ] ( see appendix c ) and the boundaries of the plateau are specified by the equilibrium condition @xmath417 while the end - to - end distance is : @xmath418 $ ] .    in this case
the @xmath6 approximation leads to ( see appendix d and section ivb ) @xmath419 { \cal l}(\alpha)$ ] .
the domain boundaries are specified by @xmath420{\cal l } ( \alpha _ { \pm } ) /4n_{h}$ ] with @xmath421 and the end - to - end distance is : @xmath422 $ ] .                to characterize the slope we expand @xmath244 around @xmath250 where @xmath251 . in the immediate vicinity of @xmath250
we have @xmath426 where @xmath427 is a constant .
accordingly , @xmath428 $ ] and @xmath429\right\}$ ] leading to @xmath430 .
the variational eigenstate @xmath431 $ ] where @xmath289 is the normalization factor . with @xmath432 , the normalization condition @xmath433 yields @xmath434 . since there is no @xmath435 dependence @xmath436 where @xmath437 is the appropriate unit vector .
consequently ( [ eq_psi ] ) leads to @xmath438