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: the statistical mechanics of two polymers with constraints on their winding angle has been extensively studied in order to understand the behavior of physical polymer systems , like for instance biological macromolecules of dna @xcite or liquid crystals composed of stacks of disk - shaped molecules @xcite , see refs . a detailed review on the subject , together with interesting proposals of how to include in the treatment of topologically entangled polymer link invariants which are more sophisticated than the winding number , can be found in @xcite . up to now , however , despite many efforts , mainly ideal polymer chains or loops winding around each others have been considered , while the repulsive interactions between the monomers have been treated approximatively or exploiting in a clever way scaling arguments integrated by numerical simulations , as for instance in @xcite . here we concentrate ourselves on the case of two directed polymers interacting via a repulsive dirac @xmath0function potential @xcite . we are particularly interested in the average degree of entanglement of the system , which we wish to estimate by computing the square average winding angle of the two polymers . this quantity is also called second moment of the winding angle or simply second moment and is a special example of the topological moments first introduced in ref .  @xcite . to achieve our goals , we develop an approach , which combines quantum mechanical and field theoretical techniques . with respect to previous works , we are able to obtain exact results even if the repulsive interactions are not switched off .    in principle , the average of any observable like the squared winding angle can be derived once the partition function of the system is known , but in our case it turns out that the partition function is simply too complicated to obtain any analytical result . this happens essentially because the full @xmath0function potential is not a central potential , since it mixes both radial and angular variables . for this reason , the usual procedure of going to polar coordinates and then solving the differential equation satisfied by the partition function of the entangled polymers with the method of separation of variables @xcite , does no longer produce simple formulas as in the situations in which only central forces are present .    to avoid these difficulties , one possibility is to approximate the @xmath0function potential with some radial potential , like for instance the hard core potential of ref . however , here we shall adopt a different strategy , based on field theories , which does not require any approximation . this strategy has been developed in @xcite ( see also ref . @xcite for more details ) to cope with ideal closed polymers whose trajectories are concatenated . also such systems are characterized by a non - central potential , which comes out as a consequence of the topological constraints imposed on the trajectories . in the field theoretical formulation of refs . @xcite the computation of the second moment is reduced to the problem of computing some correlation functions of a field theory . a bonus is provided by the fact that this computation requires just a finite number of feynman diagrams to be evaluated . in the present case , due to the presence of the @xmath0function potential , the field theory which we obtain is no longer free as that of refs . nevertheless , we will see that the theory is still linear and thus it can be exactly solved once its propagator is known . luckily , this propagator may be computed exactly using powerful non - perturbative techniques developed in the context of quantum mechanics to deal with hamiltonian containing @xmath0function potentials , see refs . basically , starting from the green function of a particle whose dynamics is governed by a given hamiltonian @xmath1 , these techniques provide an algorithm to construct the green function of a particle corresponding to a _ perturbed _ hamiltonian @xmath2 , where @xmath3 is the @xmath0function potential . one advantage of these methods is that there is a long list of potentials for which the green functions of the unperturbed hamiltonians @xmath4 are known . in this way , it is easy to generalize our treatment including new interactions , which could be relevant in polymer physics , like for instance the coulomb interaction . the price to be paid is that the quantum - mechanical algorithm works when the green functions are expressed as functions of the energy instead of the time . in the polymer analogy , assuming that the ends of the polymers are attached to two planar surfaces perpendicular to the @xmath5axis and located at the positions @xmath6 and @xmath7 , the role of time is played by the distance @xmath8 , while the energy corresponds to the chemical potential conjugated @xmath8 . to recover the original dependence on @xmath8 , one needs to calculate an inverse laplace transform of the field propagator with respect to the energy . in general , this is not a simple task . once the propagator of the linear field theory is known , the correlation functions which enter in the expression of the second moment may be calculated contracting the fields in all possible ways using the wick prescription . at the end , we get in this way an exact formula of the second moment as a function of the _ energy _ , which , we remember , has here the meaning of the chemical potential conjugated to the distance @xmath8 . in the @xmath8 space , due to the problems of computing the inverse laplace transform of the propagator mentioned above , only an approximated expression of the second moment will be given in the limit of large values of @xmath8 and assuming that the strength of the @xmath0function potential is weak enough to allow a perturbative approach . our results allow both a qualitative and quantitative understanding of the way in which the repulsive interactions affect the entanglement of two directed polymers . the corrections introduced by these interactions in the expression of the second moment of ideal polymers have been studied in some interesting limits . first of all , it has been examined the limit of long polymer trajectories , in which we show that repulsive interaction become particularly relevant . moreover , we have investigated also the perturbative regime and the strong coupling limit , which is important to recover the excluded volume interactions . while it is not a problem to take the strong coupling limit within our exact treatment of the repulsive interactions , it turns out that , in this case , the expression of the second moment is particularly complicated from the analytical point of view . for this reason , in the conclusions we will discuss the application of a powerful perturbative method to study field theories at strong coupling due to kleinert @xcite . finally , the consistency of our results with the previous ones has been checked by studying the limit of ideal polymers . the material presented in this paper is divided as follows . in the next section , the problem of computing the second moment of the winding angle of two directed polymers interacting via a @xmath0function potential is briefly illustrated using the path integral approach . a constraint on the winding angle is imposed by coupling the trajectories of the polymers with a suitable external magnetic field , following the strategy of previous works like for instance @xcite . in section [ sec : three ] the second moment is expressed in the form of a finite sum of amplitudes of a linear field theory . these amplitudes may be computed once the propagator of the theory is constructed . in our case , the propagator coincides with the green function of a particle diffusing in a @xmath0function potential . the derivation of this green function in the energy representation using non - perturbative techniques developed in the context of quantum mechanics @xcite is the subject of section [ sec : four ] . the @xmath0function potential is responsible of the appearance of singularities in the propagator at short distances , which have been regulated here with the introduction of a cut - off . this procedure is motivated by the fact that in polymer physics there is no point in considering distances which are smaller than the dimensions of a monomer . a comparison with the more rigorous method of renormalization is made , showing the consistency of the two procedures . the propagator derived in section  ( [ sec : four ] ) has a particularly nice form , in which the contributions coming from the repulsive forces can be separated from the free part of the propagator , which is related to the random walk of ideal polymers . this splitting of the propagator is used in section [ sec : five ] to discuss qualitatively and qualitatively the effects of the @xmath0function interactions on the entanglement of the system . the results of sections [ sec : three ] and [ sec : four ] provide in principle all the ingredients of the second moment . however , the amplitudes of the linear field theory derived in section  [ sec : three ] should still be evaluated . in this task one encounters the typical problems occurring in the evaluation of the analytical expressions of feynman diagrams . in the case of the second moment there are just tree diagrams , but still one has to perform complicated integrations over the spatial coordinates which are transverse to the @xmath5axis . even assuming that polymers are ideal , the analytical evaluation of these integrations requires drastic approximations , see for instance @xcite . to avoid these difficulties , we average the second moment with respect to the positions of the endpoints of the two polymers . this averaged version of the second moment can be computed without any approximation in the energy representation . this is done in section [ sec : six ] . the expression of the averaged second moment in the @xmath9space is provided instead only at the first perturbative order in the strength of the repulsive potential and assuming additionally that the value of @xmath8 is large . we give also an exact formula of the second moment without performing any averaging procedure as a function of @xmath8 . this formula is however explicit only up to the calculation of the inverse laplace transform of the propagator derived in section [ sec : four ] . in section [ sec : seven ] we consider the situation in which the polymers are not interacting in order to allow the comparison with previous results . finally , the discussion of the obtained results and ideas for further developments are presented in the conclusions . our starting point is the action of two directed polymers @xmath10 and @xmath11 : @xmath12 \label{actnorm}\ ] ] where @xmath13 is the potential : @xmath14 the sign of @xmath15 has been chosen in such a way that the interaction associated to the potential @xmath16 is repulsive . the parameters @xmath17 and @xmath8 determine the average length of the trajectories of the polymers . the ends of the polymers are supposed to be fixed on two surfaces perpendicular to the @xmath5axis and located at the heights @xmath6 and @xmath7 . both polymers have a preferred direction along the @xmath18 direction . the vectors @xmath19 and @xmath20 , @xmath21 , measure the polymer displacement along the remaining two directions of the space . the action of eq . ( [ actnorm ] ) resembles that of two quantum particles in the case of imaginary time @xmath18 . to stress these analogies with quantum mechanics , the @xmath5variable will be treated as a pseudo - time and renamed using from now on the letter @xmath22 instead of @xmath18 .    in the system of the center of mass : @xmath23 the action ( [ actnorm ] ) becomes : @xmath24 \label{cmact}\ ] ] the motion of the center of mass , which is a free motion described by the coordinate @xmath25 , will be ignored . we consider the partition function of the above two - polymer system with the addition of a constraint on the entanglement of the trajectories : @xmath26 } \delta(m-\chi ) \label{totparfun}\ ] ] @xmath27 is the so - called winding angle . its expression is given by : @xmath28 where @xmath29 is a vector potential with components : @xmath30 in the above equation we have represented the vector @xmath31 using cartesian coordinates @xmath32 . i.  e. @xmath33 . moreover , from now on , middle latin indices @xmath34 will label the directions which are perpendicular to the @xmath35axis . the definition of the partition function @xmath36 is completed by the boundary conditions at @xmath37 and @xmath38 : @xmath39 the quantity in eq . ( [ wnum ] ) becomes a topological invariant if the polymer trajectories are closed . in the present case , in which the trajectories are open , @xmath27 just counts the angle with which one polymer winds up around the other . thus , the partition function @xmath36 gives the formation probability of polymer paths winding up of an angle @xmath40    exploiting the fourier representation of dirac @xmath0functions @xmath41 eq . ( [ totparfun ] ) can be rewritten as follows : @xmath42 where @xmath43 the lagrangian @xmath44 is that of a particle immersed in the magnetic potential associated to the vector field ( [ vecpot ] ) : @xmath45 the fourier transformed partition function @xmath46 is the grand canonical version of the original partition function @xmath36 , in which the number @xmath47 is allowed to take any possible value . @xmath46 coincides with the propagator @xmath48 , which satisfies the following pseudo - schrdinger equation : @xmath49{\cal g}_\lambda(l;{{\mbox{\bf r}}}_1,{{\mbox{\bf r}}}_0 ) = 0 \label{pseudseq}\ ] ] @xmath50 is the hamiltonian of the system , computed starting from the lagrangian ( [ newaction ] ) : @xmath51 eq . ( [ pseudseq ] ) is completed by the boundary condition at @xmath52 : @xmath53    the average degree of entanglement of the two polymers can be estimated computing the topological moments of the winding angle @xmath54 , @xmath55 @xcite . once the partition function is known , the @xmath54 may be expressed as follows : @xmath56 the quantities @xmath54 depend on the boundary conditions @xmath57 and , of course , on the parameters @xmath17 and @xmath8 . for practical reasons , we will also consider the following _ averaged _ topological moments : @xmath58 as eq . ( [ aderr ] ) shows , the average is performed with respect to the relative positions @xmath57 of the endpoints . this is equivalent to an average over the positions of the endpoints @xmath59 at the instants @xmath37 and @xmath38 , because the coordinates of the center of mass have been factored out from the partition function and thus they do not play any role . the advantage of the averaged topological moments is that , a posteriori , it will be seen that their computation is easier than that of the topological moments given in eq . ( [ sectopmomdef ] ) . from the physical point of view , the averaged topological moments measure the entanglement of two polymers , whose ends at the instants @xmath37 and @xmath38 are free to move .    here we will be interested only in the second moment @xmath60 and in the averaged second moment @xmath61 , i. e. in the case @xmath62 of eqs . ( [ sectopmomdef ] ) and ( [ aderr ] ) . the second moment is in fact enough in order to estimate the formation probability of entanglement with a given winding angle and to determine how the winding angle grows with increasing polymer lengths .    in the following it will be useful to work in the so - called energy representation , i.  .e considering the laplace transformed of the partition function @xmath48 with respect to @xmath8 : @xmath63 the new partition function @xmath64 describes the probability of two entangled polymers of any length subjected to the condition that the relative positions of the polymer end at the initial and final instants @xmath65 and @xmath66 are given by the vectors @xmath67 and @xmath68 . with respect to the formulation in the @xmath9 space , however , the distance @xmath69 is no longer exactly equal to @xmath8 , but is allowed to vary according to a distribution which is governed by the boltzmann - like factor @xmath70 . thus , @xmath71 plays the role of the chemical potential conjugated to the end - to - end distance of the polymer trajectories in the @xmath35direction . it is worth to remember that , roughly speaking , small values of @xmath71 correspond to large values of @xmath8 , while large values of @xmath71 correspond to small values of @xmath8 . starting from eq . ( [ pseudseq ] ) and recalling the boundary conditions ( [ bbdc ] ) , it is easy to check that @xmath64 satisfies the stationary pseudo - schrdinger equation : @xmath72{\cal g}_\lambda(e;{{\mbox{\bf r}}}_1,{{\mbox{\bf r}}}_0)=\delta({{\mbox{\bf r}}}_1-{{\mbox{\bf r}}}_0 ) \label{pseenerep}\ ] ] where @xmath50 is always the hamiltonian of eq . ( [ ham ] ) . in this section we wish to evaluate the expression of the second moment as a function of the energy @xmath71 using a field theoretical formulation of the polymer partition function . the starting point is provided by the formula of the second moment in the @xmath9space suitably rewritten in the following way : @xmath73 for consistency with eq . ( [ sectopmomdef ] ) , the numerator @xmath74 and the denominator @xmath75 appearing in eq . ( [ secmomnumden ] ) must be of the form : @xmath76 and @xmath77    using eq . ( [ laplrep ] ) , it is now straightforward to compute the laplace transform of @xmath74 and @xmath75 : @xmath78 @xmath79 once the functions @xmath80 and @xmath81 are known , one can construct the ratio : @xmath82 which is nothing but the second moment of the winding angle expressed as a function of the chemical potential @xmath71 . we remark that the green function @xmath64 is related to the feynman propagator of the spin@xmath83 aharonov - bohm problem in quantum mechanics . in principle , this green function can be computed exactly starting from eq . ( [ pseenerep ] ) @xcite , but its final expression is too complicated for our purposes . moreover , the method used in @xcite to renormalize the singularities coming from the presence of the @xmath0function potential is valid only in a restricted region of the domain of @xmath84 . this is incompatible with our requirements , because , to derive the second moment , one has to integrate @xmath64 with respect to @xmath84 over the whole real line . for this reason , we prefer here to use a field theoretical representation of this green function . this is achieved by noting that @xmath85 coincides with the inverse matrix element of the operator @xmath86 : @xmath87 and may be expressed in a functional integral form in terms of replica fields : @xmath88 in the above equation @xmath89 are multiplets of replica fields : @xmath90 with action @xmath91\psi \label{erepact}\ ] ] the symbol @xmath92 in eq . ( [ ftrep ] ) denotes summation over the replica index . for example @xmath93 . below it will be used also the convention @xmath94 . the details of the derivation of eq . ( [ ftrep ] ) can be found in previous publications on the subject @xcite and will not be provided here .    in order to proceed , it will be convenient to expand the action ( [ erepact ] ) in powers of @xmath84 : @xmath95 where we have put : @xmath96 \label{szero}\ ] ] @xmath97 \label{sone}\ ] ] @xmath98    at this point we come back to the computation of the quantities @xmath80 and @xmath81 appearing in the expression of the second moment . exploiting the new form of the partition function given by eqs . ( [ ftrep][stwo ] ) , together with the relation @xmath99 and the fact that @xmath100 , it is possible to rewrite eqs . ( [ momno ] ) and ( [ momnt ] ) as follows : @xmath101e^{-s_0(\psi^*,\psi ) } \label{nsdo}\ ] ] @xmath102 the right hand sides of eqs . ( [ nsdo ] ) and ( [ nsdt ] ) represent vacuum expectation values of a field theory governed by the action @xmath103 of eq . ( [ szero ] ) . in the formulation in terms of quantum operators we have : @xmath104 @xmath105 the correlation functions have a subscript @xmath106 to remember that , according to the replica method , they should be computed first assuming that the number of replicas @xmath106 is an arbitrary positive integer and then taking the limit for @xmath106 going to zero .    the above correlators may be evaluated using standard field theoretical methods . one could be tempted to use a perturbative approach assuming that the value of @xmath15 appearing in the action @xmath103 of eqs . ( [ nsdo ] ) and ( [ nsdt ] ) is small , but this is not necessary . as a matter of fact , if it is true that @xmath103 does not describe free fields because of the presence of the @xmath0function potential , it is also true that it is just quadratic in the fields . as a consequence , one is allowed to define a propagator @xmath107 associated with this action . it is easy to check that @xmath107 satisfies the equation : @xmath108g(e;{{\mbox{\bf x}}},{{\mbox{\bf y}}})=\delta({{\mbox{\bf x}}},{{\mbox{\bf y } } } ) \label{fdfjh}\ ] ] using the above propagator , one can evaluate the amplitudes in eqs . ( [ nsdoqop ] ) and ( [ nsdtqop ] ) exactly by contracting the fields in all possible ways according to the wick theorem . after straightforward calculations , one finds : @xmath109 @xmath110 @xmath111 where @xmath112\label{amptwoone}\\ i_2({{\mbox{\bf r}}}_1,{{\mbox{\bf r}}}_0 ) & = & + \frac1{2c^2}\int d^2xd^2y\left [ a_i({{\mbox{\bf x}}})(\nabla_{\mathbf x}^ig(e;{{\mbox{\bf x}}},{{\mbox{\bf r}}}_1))g(e;{{\mbox{\bf y}}},{{\mbox{\bf x } } } ) a_j({{\mbox{\bf y}}})(\nabla_{\mathbf y}^jg(e;{{\mbox{\bf r}}}_0,{{\mbox{\bf y}}}))\right]\label{amptwotwo}\\    i_3({{\mbox{\bf r}}}_1,{{\mbox{\bf r}}}_0)&= & + \frac1{2c^2}\int d^2xd^2y\left [ a_i({{\mbox{\bf x}}})g(e;{{\mbox{\bf r}}}_0,{{\mbox{\bf x}}})(\nabla_{\mathbf x}^i \nabla_{\mathbf y}^jg(e;{{\mbox{\bf x}}},{{\mbox{\bf y}}}))a_j({{\mbox{\bf y}}})g(e;{{\mbox{\bf y}}},{{\mbox{\bf r}}}_1)\right]\label{amptwothree}\\ i_4({{\mbox{\bf r}}}_1,{{\mbox{\bf r}}}_0 ) & = & -\frac1{2c^2}\int d^2xd^2y \left [ a_i({{\mbox{\bf x}}})g(e;{{\mbox{\bf r}}}_0,{{\mbox{\bf x}}})(\nabla_{\mathbf x}^ig(e;{{\mbox{\bf x}}},{{\mbox{\bf y } } } ) ) a_j({{\mbox{\bf y}}})(\nabla_{\mathbf y}^jg(e;{{\mbox{\bf y}}},{{\mbox{\bf r}}}_1))\right]\label{amptwofour}\end{aligned}\ ] ] from the physical point of view , the above equations may be interpreted in the following way . the fields @xmath113 and @xmath114 contain operators which , inside each replica sector , create and annihilate segments of the two polymers , whose relative positions are given by the vector @xmath115 . the two polymer system has been projected in the two - dimensional plane perpendicular to the @xmath35axis . for this reason , there appear only the transverse coordinates @xmath115 . the only remnant of the third dimension is the dependence on the energy @xmath71 . the correlation functions ( [ ampone ] ) and ( [ amptwo ] ) describe the fluctuations of the two polymers immersed in the @xmath0function potential and subjected to the interactions represented by the vector potential ( [ vecpot ] ) . we recall that the origin of the latter interactions is the presence of the constraint on the winding angle in the partition function ( [ totparfun ] ) . to evaluate the correlation functions ( [ ampone ] ) and ( [ amptwo ] ) , one needs to consider only a finite number of feynman diagrams , corresponding to the relevant processes with which the two polymers interact together . the result , after the analytical evaluation of these diagrams , is provided by eqs . ( [ amptwoone][amptwofour ] ) . let us note that in these equations the repulsive interactions due to the @xmath0function potential are hidden in the propagators @xmath107 . the feynman diagrams related to the amplitudes of eqs . ( [ amptwoone])([amptwofour ] ) are all topologically equivalent to the diagram of fig . [ processone ] .    ) ([amptwofour ] ) . the two polymers @xmath10 and @xmath11 start at a distance @xmath116 from each other and interact twice with the the external field @xmath117 . at the end , the relative position of the end points at the instant @xmath38 is given by @xmath68 . the three - vertices appearing in the figure are related to the interaction described by eq . ( [ sone ] ) . , title="fig : " ] +    the amplitude of eq . ( [ amponek ] ) is related instead to the feynman diagram of fig . [ processtwo ] .    ) . the two polymers @xmath10 and @xmath11 start at a distance @xmath116 from each other and interact with the the external field @xmath117 . at the final instant @xmath38 , the relative position of the end points is given by @xmath68 . the four - vertex appearing in this figure is related to the interaction described by eq . ( [ stwo ] ) . , title="fig : " ] +    the vectors @xmath118 and @xmath119 denote the relative positions of the end points of the two polymers at the initial and final instants , as already mentioned . the integration variables @xmath120 and @xmath121 appearing in eqs . ( [ amponek][amptwofour ] ) may be regarded as the vectors which give the relative positions of the trajectories of the two polymers at the instants in which they interact together via the external vector potential @xmath122 of eq . ( [ vecpot ] ) . there is no restriction on the domain of integration of @xmath120 and @xmath121 , so that the components of these relative vectors are allowed to take any value . this implies that the distance between the polymer segments when the interaction with @xmath122 occurs can be arbitrarily large .    now that the correlation functions which are present in the expressions of @xmath123 and @xmath124 given in eqs . ( [ nsdoqop ] ) and ( [ nsdtqop ] ) have been evaluated , see eqs . ( [ amponek][amptwofour ] ) , we may put everything together and give to the second moment of eq . ( [ secmomerep ] ) a more explicit form : @xmath125 in conclusion , the initial problem of computing the second moment of the winding angle @xmath60 has been reduced to the evaluation of a finite number of integrals , which are given in eqs . ( [ amponek][amptwofour ] ) . of course , to make these integrals really explicit , we still need to derive the propagator @xmath107 , which is so far the only missing ingredient . this will be done in the next section . let @xmath126 be the solution of the differential equation : @xmath127 for a given hamiltonian @xmath4 . when @xmath52 , @xmath126 satisfies the boundary condition : @xmath128 in the case of a hamiltonian @xmath50 , obtained by adding to @xmath1 a @xmath0function potential as a perturbation : @xmath129 we consider the analogous differential problem : @xmath130 @xmath131 we wish to compute @xmath132 starting from the green function @xmath126 , which is supposed to be known . it is possible to show that @xmath132 and @xmath126 are related by the integral equation @xcite : @xmath133 we see that in the right hand side of the above equation the presence of the @xmath0function forces us to consider the functions @xmath126 and @xmath132 evaluated at the points @xmath134 and/or @xmath135 . usually , at these points green functions may be not well defined due to the presence of singularities . a concrete procedure to remove these singularities will be indicated later . for the moment , we go further with formal manipulations , assuming that some kind of consistent regularization of the possible divergences has been introduced .    first of all , we perform the integration over @xmath136 in eq . ( [ gkk ] ) : @xmath137 the integral in @xmath138 appearing in the right hand side of eq . ( [ beflap ] ) is a convolution which can be better treated after a laplace transform . thus , we transform both sides of this equation with respect to @xmath8 : @xmath139 where @xmath140 and @xmath141 at this point , it is easy to extract from eq . ( [ lapltragf ] ) the expression of @xmath107 : @xmath142 the above formula may be used in order to solve eq . ( [ fdfjh ] ) . in this case , @xmath4 coincides with the free action : @xmath143 and the function @xmath144 is given by : @xmath145 here @xmath146 denotes the modified bessel function of the second kind of order zero . clearly , we can not apply directly eq . ( [ finald ] ) without introducing a regularization . as a matter of fact , if not treated , the naive denominator in the second term of the right hand side is equal to infinity , i. e. @xmath147 . this is due to the fact that @xmath146 diverges logarithmically in the limit @xmath148 : @xmath149    a natural regularization is suggested by the fact that , in polymer physics , it has no sense to consider lengths which are smaller than the size of the molecules which compose the polymers . thus , it seems reasonable to regulate ultraviolet divergences by introducing a cut - off @xmath150 at short distances . the length @xmath150 is comparable with the molecular size . according to this prescription , by inserting the green function of eq . ( [ gzspec ] ) in eq . ( [ finald ] ) , we obtain : @xmath151 the symbol @xmath152 means that the quantity in the left hand side of an equation has been replaced in the right hand side with its regulated version . the above green function is what we need in order to evaluate explicitly the amplitudes of eqs . ( [ ampzero][amptwo ] ) . the infinities coming from the @xmath0function potential should be treated with some care in order to avoid ambiguities . for this reason , we would like to compare the naive prescription used here to derive eq . ( [ gffin ] ) with the more rigorous procedure of renormalization . it is known in fact that the renormalization of the infinities coming from @xmath0function interactions produces physically sensible results @xcite . the divergences will be regulated introducing a cut - off @xmath153 in the momentum space . as a consequence , it will be convenient to express the free green function of eq . ( [ gzspec ] ) in momentum space . to this purpose , we use the following formula : @xmath154 to evaluate the green function at the singular point @xmath155 we need to compute the following divergent integral : @xmath156 using the above cut - off prescription to eliminate the ultraviolet singularities we get , in the assumption @xmath157 : @xmath158 now , according to the spirit of renormalization , we subtract the infinities from the bare parameters of the theory . in our case , after choosing an arbitrary mass scale @xmath159 , which gives the renormalization point , we renormalize the bare coupling constant @xmath15 . actually , it will be better to call it @xmath160 instead of @xmath15 in order to distinguish it from the effective coupling constant @xmath15 appearing in eq . ( [ gffin ] ) . the subtraction of infinities is performed in such a way that the quantity : @xmath161 becomes finite . we choose a sort of minimal subtraction scheme , in which the renormalized coupling constant @xmath162 is related to the bare coupling constant @xmath160 as follows : @xmath163 applying the last two above equations back to eq .  ( [ finald ] ) , we get as a result : @xmath164 eqs . ( [ gffin ] ) and ( [ gffinren ] ) are reciprocally compatible . in fact , since @xmath150 is very small , because it is the smallest possible length scale in our polymer problem , one can use the following approximation ( see eq . ( [ fdsf ] ) ) in the denominator of the second term of eq . ( [ gffin ] ) : @xmath165 comparing with the analogous denominator in eq . ( [ gffinren ] ) , it is possible to relate @xmath150 with the inverse of the mass @xmath159 : @xmath166 moreover , the effective coupling constant @xmath15 of eq .  ( [ gffin ] ) may be identified with the renormalized coupling constant @xmath162 , which gives the strength of the repulsive interaction ( [ dirpot ] ) at distance scales of order @xmath150 .    before concluding this section , we make a small digression about the translational invariance of the free hamiltonian ( [ freeham ] ) and consequently of the free green function ( [ gzspec ] ) . clearly , this is not the same translational invariance that was already present in the original action ( [ actnorm ] ) due to the translational invariance of the potential ( [ dirpot ] ) . this new invariance is rather related to the fact that the physics of the two polymer system in the absence of any interaction does not change if we modify the relative positions of the polymer ends at @xmath37 and @xmath38 in a symmetric way . an example of such transformations is the translation of both ends of polymer @xmath10 at the initial and final points by a constant vector @xmath167 : @xmath168 as a result of the translations ( [ strazerexone][straelexone ] ) , the relative vector @xmath169 of eq .  ( [ cmcoor ] ) at the instants @xmath37 and @xmath38 changes as follows : @xmath170 clearly the propagator ( [ gzspec ] ) is invariant under the above transformations . this kind of invariance can be explained as follows . as far as the two polymers @xmath10 and @xmath11 do not interact , each of them may be treated as an independent system . if we translate for instance both ends of polymer @xmath10 at @xmath37 and @xmath38 in the symmetrical way shown by eqs . ( [ strazerexone ] ) and ( [ straelexone ] ) , the number of configurations of polymer @xmath10 and consequently the configurational entropy of the whole system do not change , because the transformation is equivalent to a translation of polymer @xmath10 in the space . of course , this invariance disappears as soon as the two polymers start to interact or if they are entangled together . indeed , if one adds to the free hamiltonian ( [ freeham ] ) a @xmath0function potential , the propagator stops to be translational invariant as shown by the green function of eq . ( [ finald ] ) , which does not depend on the difference @xmath171 . in principle we have at this point all the ingredients which are necessary to compute the second moment of eq . ( [ secmomerep ] ) . in eqs . ( [ nsdoqop ] ) and ( [ nsdtqop ] ) , in fact , the quantities @xmath80 and @xmath81 are written as linear combinations of the amplitudes of eqs . ( [ ampzero][amptwo ] ) , which can be explicitly evaluated using the propagator @xmath172 given in eq . ( [ gffin ] ) appearing in the definition of the propagator of eq . ( [ gffin ] ) are replaced by the vectors @xmath173 . this is to avoid confusions with the notation of eqs . ( [ amponek][amptwofour ] ) , where the same pair of vectors @xmath174 has been exploited to denote dummy integration variables . ] and the formulas of eqs . ( [ amponek][amptwofour ] ) . the remaining task is to perform the integrations over the coordinates @xmath115 and @xmath175 in eqs . ( [ amponek][amptwofour ] ) . from the analytical point of view , the evaluation of these integrals poses severe technical problems , which can be solved only with the help of drastic approximations . however , the difficulties become milder if we average the second moment over the endpoints of the polymers as shown in eq . ( [ aderr ] ) . in the energy representation , which we are using , this means that we have to consider the following averaged version of the second moment in eq . ( [ secmomerep ] ) : @xmath176 where @xmath177 accordingly , we need to integrate the quantities @xmath178 and @xmath179 , @xmath180 of eqs . ( [ amponek][amptwofour ] ) with respect to @xmath68 and @xmath67 . putting : @xmath181 we obtain from eqs . ( [ ampone ] ) and ( [ amptwo ] ) the following expressions of @xmath182 and @xmath183 : @xmath184 it will also be convenient to split the propagator @xmath185 of eq . ( [ gffin ] ) into two contributions : @xmath186 where @xmath187 is the free propagator of eq . ( [ gzspec ] ) , which is invariant with respect to the transformations ( [ strainvzer ] ) and ( [ strainvone ] ) , while @xmath188 in the above equation we have isolated in the expression of @xmath189 the factor : @xmath190 it is clear that the origin of the term @xmath189 in the propagator is due to presence of the @xmath0function interaction ( [ dirpot ] ) in the polymer action ( [ actnorm ] ) . in fact , if @xmath191 , this term vanishes identically . thus , using the splitting of the propagator of eq . ( [ propsplit ] ) , it is now possible to separate in eqs . ( [ amponek][amptwofour ] ) the contributions to entanglement given by repulsive forces . it seems natural to expand the quantities @xmath183 , @xmath192 and @xmath193 defined in eqs . ( [ deee ] ) , ( [ koneee ] ) and ( [ ionee ] ) with respect to @xmath194 as follows : @xmath195 where the superscript @xmath196 , with @xmath197 , denotes the order in @xmath194 . there are no higher order terms with @xmath198 , so the above expansions are exact . it is easy to show how @xmath192 and the @xmath193 s depend on the pseudo - energy @xmath71 . after a rescaling of the integration variables @xmath199 and @xmath175 in eqs . ( [ ennee ] ) and ( [ deee ] ) , one finds in fact that : @xmath200 @xmath201 where the factors @xmath202 s and the @xmath203 s are functions of the parameters @xmath150 and @xmath17 , but not of @xmath71 or @xmath15 . in fact , the coupling constant @xmath15 appears only inside the powers of @xmath194 . let us note in eqs . ( [ behaone ] ) and ( [ behatwo ] ) the presence of the overall factor @xmath204 in eqs . ( [ behaone ] ) and ( [ behatwo ] ) . looking at eq . ( [ enneeexp ] ) , it is clear that the whole function @xmath182 is characterized by the leading scaling behavior @xmath205 . in the @xmath9space , after an inverse laplace transform , this behavior corresponds to the following scaling law , which is typical of ideal polymers : @xmath206 . the powers of @xmath194 , appearing in the expressions of @xmath207 and @xmath208 , introduce corrections to this leading behavior that are at most logarithmic in @xmath71 . as a matter of fact , if the condition @xmath209 is satisfied , we have that : @xmath210 naively , the above seems the only logarithmic correction which is possible in the expressions of @xmath182 and @xmath183 when @xmath71 is small . however , that this is not true . in fact , in deriving eqs . ( [ behaone ] ) and ( [ behatwo ] ) , we have not considered the divergences which arise in some of the integrations over the variables @xmath211 and @xmath68 . after regulating these divergences with some prescription , as for instance the ultraviolet cut - off @xmath150 used in eq . ( [ gffin ] ) , we will see in section  [ sec : six ] that the naive rescaling of variables exploited in order to obtain eqs . ( [ behaone ] ) and ( [ behatwo ] ) does no longer work and one should add extra logarithmic factors to these equations . ( [ behaone ] ) and ( [ behatwo ] ) may be also useful to study the case of polymers in confined geometries . as a matter of fact , for large values of @xmath71 , one recovers the limit of small values of @xmath8 , in which the region between the initial and final height is very narrow . looking at eqs . ( [ behaone ] ) and ( [ behatwo ] ) , it is clear that the only interesting corrections when @xmath71 is large come from the powers of @xmath194 . to evaluate these corrections , one should note that the modified bessel function of the second kind @xmath146 goes very fast to zero for large values of @xmath18 . as a consequence , already in the domain of parameters in which @xmath212 , it is possible to make the very interesting approximation @xmath213 unfortunately , it turns out that the values of the energy for which the above equation is satisfied are not physical , as it will be shown below . other useful information on the influence of repulsive forces on the winding angle can be obtained studying the form of the function @xmath189 of eq . ( [ nontrainvpro ] ) . we remember in fact that all the effects of the repulsive forces are concentrated in this component of the propagator . supposing for example that the value of @xmath214 is very large , i.  e. : @xmath215 we have the following approximate expression of @xmath216 : @xmath217 a relation analogous to eq . ( [ appexplx ] ) may be written also for the variable @xmath218 . in practice , ( [ appexplx ] ) means that the repulsive interactions do not play any particularly relevant role in polymer configurations in which the ends of the trajectories at some point are very distant . this is not a surprise . if the trajectories at some height @xmath22 are very far from each other , they will have little or no chance to interact together via the repulsive interactions of eq . ( [ dirpot ] ) , which are of short range . ( [ appexplx ] ) gives the concrete law with which the contributions of the repulsive forces are suppressed in configurations of this kind . in particular , if the distance between the trajectories is much greater than the _ characteristic length scale _ @xmath219 the influence of the repulsive forces ceases to be relevant . of course , even if at some points the trajectories are very distant , polymers will always have a chance to get near enough to be able to interact if they are sufficiently long . as a consequence , we expect that the characteristic length @xmath220 increases with the increasing of the lengths of the trajectories . it is easy to check that this is exactly the case . to show that , let us consider the dependence of @xmath220 on the polymer length . one parameter which determines this length is the distance @xmath8 between the ends of the polymers along the @xmath35axis . indeed , a trajectory connecting the two ends of a polymer must be very long if these ends are located at very distant heights . in the energy representation , large values of @xmath8 correspond to small values of @xmath71 . for example , in the limit @xmath221 , which corresponds to infinite polymer lengths , we have that @xmath222 , confirming our intuitive expectations . another confirmation comes from eq . ( [ randisell ] ) below , where a rough estimation of the behavior of @xmath220 with respect to the distance @xmath8 is given . the dependence on @xmath8 is not the whole story . as a matter of fact , during their random walk in the @xmath35direction from @xmath37 to @xmath38 , polymers are also allowed to wander in the remaining two directions . loosely speaking , the variations in the length of the trajectories associated to the fluctuations in these transverse directions are taken into account by the parameter @xmath17 . smaller values of @xmath17 correspond to longer trajectories and vice - versa , see @xcite . it is now easy to realize from eq . ( [ chalensca ] ) that the characteristic length @xmath220 increases when @xmath17 decreases as expected . taking into account all the above considerations , it is possible to conclude that repulsive forces give relevant contributions to the second moment only in the case of configurations of the system in which the trajectories of the two polymers are not too far from each other . as a matter of fact , in the propagators appearing in the amplitudes of eqs . ( [ amponek][amptwofour ] ) all configurations in which the distance between the trajectories at some height in the @xmath35axis is bigger than a few characteristic lengths @xmath220 are exponentially suppressed according to eq . ( [ appexplx ] ) . one may also add that this suppression becomes milder in the case of long polymers , because we have seen that the characteristic lengths @xmath220 grows with the length of the polymers with a law which has been given in eq . ( [ randisell ] ) .    in the rest of this section we will analyze some interesting limiting cases , in which repulsive interactions become particularly weak or strong . to this purpose , it would be appealing to consider the quantity @xmath223 , where @xmath194 has been given in eq . ( [ lambe ] ) , as an energy dependent effective or running coupling constant of the repulsive interactions . this could be suggested by the expansions of eqs . ( [ dexp][expioe ] ) and by the fact that the quantity @xmath223 has the right dimension to be a coupling constant . indeed , we will see that there are cases in which the strength of @xmath194 really determines the strength of the repulsive forces . however , this is not true in general , as it should be because @xmath194 is just a parameter which has been factored out from the expression of @xmath189 and thus its meaning does not coincide with that of a running coupling constant . keeping that in mind , we start to study the perturbative regime , in which @xmath15 is very small . in the part of the propagator in which there are the contributions of the repulsive forces , i. e. the function @xmath189 , @xmath15 is present only inside @xmath194 . expanding this quantity in powers of @xmath15 , we obtain : @xmath224 we see that , at the leading order in @xmath15 , @xmath194 is proportional to @xmath15 and thus , as it could have been expected , @xmath189 may be treated as a small perturbation with respect to the free propagator @xmath225 . let s now go back to eq . ( [ fapp ] ) . in that equation it turns out that @xmath194 has the same behavior as in the perturbative regime , even if eq  ( [ fapp ] ) has been derived in the hypothesis that @xmath212 , but without supposing that @xmath15 is small . before dwelling upon the physical meaning of this coincidence , let s see what is the significance of the condition @xmath212 . to this purpose , we make the following approximations : @xmath226 as mentioned before , it is quite reasonable to assume that the length @xmath8 is proportional to the inverse of the energy @xmath71 , while the second approximation implies that polymers are very flexible . for example , in polyethylene the kuhn length@xmath227 is of the order of molecular sizes . exploiting eq . ( [ draapp ] ) , it turns out that the condition @xmath212 is equivalent to the condition @xmath228 . this would mean that our system is squeezed in a volume whose height @xmath8 is smaller than the size of a monomer . clearly , this situation is not very physical .    since we have been able to compute the exact form of the propagator @xmath229 , it is not difficult to study also the strong coupling limit @xmath230 . as in the perturbative case , the only affected part of the propagator ( [ gffin ] ) is the factor @xmath194 appearing in @xmath189 . after a trivial calculation , one finds that , in the strong coupling limit , the form of @xmath189 is given by : @xmath231 assuming that polymers are very long , let us study the left hand side of the above equation . this is a ratio of modified bessel functions of the second kind . since @xmath150 is a very small quantity and these functions have a logarithmic singularity if their argument is small , see eq . ( [ fdsf ] ) , it is licit to suppose that @xmath232 unless @xmath233 and/or @xmath234 . on the other side , we know from eq . ( [ cslm ] ) that , if @xmath235 decays exponentially . in other words , in the left hand side of eq . ( [ scnontrainvpro ] ) the denominator will dominate over the numerator whenever the distance between the polymer trajectories is not of the order of a few molecular sizes . thus , if @xmath15 is large , the major contributions to winding angle coming from the repulsive interactions occur when the trajectories are very near to each other . this could be expected from the fact that , in the strong coupling limit , one recovers the excluded volume interactions . finally , let us study the domain of the parameters @xmath71 and @xmath17 in which the condition @xmath236 is verified . we will see that this domain is particularly interesting , because if condition ( [ condd ] ) is verified , the corrections of the repulsive interactions to the entropy dominated behavior of ideal polymers become relevant . it has been already shown that under the assumption made in eq . ( [ condd ] ) , the parameter @xmath194 is approximated as in eq . ( [ logcorr ] ) . even if it is not strictly necessary , we suppose here that @xmath15 has some finite value , while polymers are so long that the following inequality is satisfied : @xmath237 this further assumption is to eliminate the dependence on @xmath15 , which could introduce confusion in the following discussion due to possible interferences of condition ( [ condd ] ) with those of the perturbative and strong coupling regimes . in the @xmath9space , ( [ fjfsfd ] ) corresponds to the inequality @xmath238 . now @xmath189 may be approximated as follows : @xmath239 as promised , the above equation does not contain the parameter @xmath15 . we see from the left hand side of eq . ( [ domhigene ] ) that the function @xmath189 is logarithmically suppressed , due to the presence of @xmath240 in the denominator . this suppression effect is counterbalanced only at short distances by the two modified bessel functions of the second kind appearing in the numerator , which diverge logarithmically whenever @xmath241 and/or @xmath242 . the total result of these opposite effects in the expression of the averaged second moment will be presented in section  [ sec : six ] after performing the explicit computation of the amplitudes of eqs . ( [ amponek][amptwofour ] ) .    to conclude this section , let us give some concrete values of the involved parameters . first of all , let us estimate the values of @xmath8 , for which the two polymer system is in the regime ( [ condd ] ) . using the approximations made in eq . ( [ draapp ] ) , we may conclude that , if the relation ( [ condd ] ) is satisfied , the length @xmath8 needs such that @xmath243 , i.  e. @xmath8 is at least of the order of hundred molecular lengths or more : @xmath244 . moreover , it is possible to give a rough estimation of the maximum distance of the end points , after which the two polymers are too far from each other to allow a relevant contribution to the winding angle due to repulsive interactions . using eq . ( [ cslm ] ) , in fact , it turns out that the repulsive interactions are relevant only in the range of distances : @xmath245 finally , the situation opposite to condition ( [ condd ] ) is not realistic , because it leads to the constraint @xmath246 . this would corresponds to the case of a polymer which is shorter than the size of the molecules composing it . at this point we are ready to compute the quantities @xmath182 and @xmath183 of eqs . ( [ ennee ] ) and ( [ deee ] ) . we start with @xmath183 . using eqs . ( [ deeeexp ] ) , ( [ dexp ] ) and the splitting ( [ propsplit ] ) of the propagator , one has at the zeroth order in @xmath194 : @xmath247 after a shift of variables , the above equation gives : @xmath248 where @xmath249 is the total surface of the system in the two dimensional space , which is transverse to the @xmath35axis . using the identity @xmath250 one finds : @xmath251 this expression of @xmath252 has the following interpretation : we are performing here an average of the second moment with respect to all possible initial and final positions of the endpoints of the polymers and @xmath183 _ counts _ the number of these configurations . the component @xmath252 of @xmath183 depends only on the free propagator @xmath144 , which is translational invariant in the sense discussed after eq . ( [ gzspec ] ) . this invariance explains why the number of configurations grows proportionally to the surface @xmath253 . the reason is that , for each configuration of the polymers , one can obtain other equivalently probable configurations by the symmetric translation of their ends on the surface @xmath253 at the initial and final instants . let us now apply to @xmath252 an inverse laplace transform , in order to go back to the @xmath9space . after a simple calculation we obtain : @xmath254 i. e. @xmath255 does not depend on @xmath8 . the next and last contribution to @xmath183 is given by : @xmath256 exploiting eq . ( [ intbas ] ) to integrate out the variables @xmath119 and @xmath257 , we get : @xmath258 we remark that the above contribution to @xmath183 vanishes in the limit @xmath191 . this could be expected due to the fact that @xmath259 collects all contributions coming from the repulsive interactions . these interactions break explicitly the translational invariance of the free part of the action and , as a consequence , @xmath259 is no longer proportional to the surface @xmath253 as @xmath252 . unfortunately , it is not easy to compute the inverse laplace transform of @xmath259 without making some approximation . to this purpose , we assume that the repulsive interactions are weak , i.  e. @xmath260 , and that the value of @xmath8 is large . in this case , since we are in the domain of small @xmath71 s , it is possible to expand @xmath259 up to the second order in @xmath15 as follows : @xmath261 in order to obtain the above equation we have used both eqs . ( [ fdsf ] ) and ( [ pertapp ] ) . the inverse laplace transform of eq . ( [ deeoneexp ] ) gives : @xmath262l+\frac c{2\pi } v_0 ^ 2l\log l \label{deeonel}\ ] ] where @xmath263 is the euler constant . putting eqs . ( [ deezerointe ] ) and ( [ contdeeone ] ) together , we obtain : @xmath264 this is an exact result . an approximated expression of @xmath265 can be derived instead from eqs . ( [ dzerol ] ) and ( [ deeonel ] ) . now we turn to the derivation of @xmath182 . we start by computing order by order in @xmath194 the contributions to the quantities @xmath192 and @xmath193 of eqs . ( [ expke ] ) and ( [ expioe ] ) respectively . at the zeroth order we have for @xmath192 : @xmath266 after performing an easy integrations over the coordinates @xmath57 , one obtains : @xmath267 the remaining integral with respect to the @xmath115 coordinate is both ultraviolet and infrared divergent and needs to be regulated . we have already seen that the singularities in the ultraviolet domain may consistently be eliminated with the introduction of the small distance cut - off @xmath150 . a large distance cut - off is instead motivated by the fact that the size of a real system is necessarily finite . implicitly , we have already used this kind of infrared regularization in eq . ( [ deezero ] ) , where we have assumed that the total surface @xmath253 of the system in the directions which are transverse to the @xmath35axis is finite . supposing that the shape of @xmath253 is approximately a disk of radius @xmath268 , so that @xmath269 , we may write : @xmath270 substituting eq . ( [ polinfultint ] ) in eq . ( [ kzerot ] ) , one obtains the following expression of @xmath271 : @xmath272 the inverse laplace transform of @xmath271 gives : @xmath273    we have now to compute the quantities @xmath274 , with @xmath180 . the expressions of the @xmath275 s may be obtained from eqs . ( [ ionee ] ) and ( [ amptwoone][amptwofour ] ) , by substituting everywhere the propagator @xmath107 with its free version @xmath144 . it is easy to show that : @xmath276 this vanishing , which is actually a double vanishing , is due to the fact that each of the @xmath275 s contains an integral of a total divergence together with an integral which is zero for symmetry reasons . for some values of @xmath277 , like for instance when @xmath278 , to isolate such integrals it is necessary to perform some integrations by parts . this is allowed because the @xmath274 s are not affected by divergences , contrarily to @xmath192 . as an example , we work out explicitly the case of @xmath279 . the first vanishing integral is the following : @xmath280 this is of course zero due to symmetry reasons . the second vanishing integral in @xmath279 is of the form : @xmath281 after performing the integration over @xmath68 with the help of a shift of variables and of eq . ( [ intbas ] ) , we have , apart from a proportionality factor : @xmath282 since @xmath283 is a divergenceless vector potential , i.  e. @xmath284 , @xmath285 can be rewritten as the integral of a total divergence : @xmath286 clearly , the left hand side of the above equation is zero . this fact can be also checked passing to the fourier representation . exploiting eq . ( [ ftsmbf ] ) and the formula @xmath287 in eq . ( [ deeint ] ) , one obtains : @xmath288 thus @xmath289 because @xmath290 . in an analogous way one shows that also @xmath291 and @xmath292 are identically equal to zero . we are now ready to compute the contributions to @xmath182 , which are linear in @xmath194 . first of all , we treat the term @xmath293 , which is given by : @xmath294 \label{konee}\ ] ] the integrations over @xmath67 and @xmath68 may be easily performed using eq . ( [ intbas ] ) and give as a result a factor which is proportional to @xmath204 . after that , only the following integral in @xmath115 remains to be done : @xmath295 here the ultraviolet divergence , which is present in the left hand side , has been regulated in the usual way with the introduction of the short distances cut - off @xmath150 . infrared divergences are absent . going to polar coordinates , the right hand side of the above equation becomes : @xmath296 putting everything together , one arrives at the final result : @xmath297 if the quantity @xmath298 is small , it is possible to derive the following asymptotic expression of @xmath293 : @xmath299 to go from eq . ( [ konefin ] ) to eq . ( [ koneeer ] ) , we have used the asymptotic formula : @xmath300 which is valid for small values of @xmath298 . we see from eqs . ( [ konefin ] ) and ( [ leaordaasyexp ] ) that the presence of ultraviolet divergences , together with the needed regularization , has modified the naive form of @xmath293 as a function of the pseudo - energy @xmath71 given in eq . ( [ behaone ] ) . the modification consists in the appearance of the factor @xmath301 , which exhibits a square logarithmic singularity in the limit @xmath302 . the inverse laplace transformed of @xmath293 can be derived only making some approximation . as in the case of @xmath259 , we will work in the double limit , in which @xmath15 is very small and @xmath8 is very large . after a few calculations we obtain : @xmath303\left ( \log s+c \right)+\right.\nonumber\\&&\left . \frac14l\log^2(2ca^2)+\frac 12 \log(2ca^2 ) \left [ ( c-1)l - l\log l \right ] \right\ } \label{ddd}\end{aligned}\ ] ]    at this point we have to compute the expressions of the @xmath304 s , @xmath305 . it is possible to show that these contributions vanish identically , i.  e. : @xmath306 the motivations of this vanishing are similar to the motivations for which there are no contributions at the zeroth order : all terms which appear in the quantities @xmath304 contain at least one integral of a total divergence or one integral , which is zero for dimensional reasons . as in the case of the @xmath274 s , there are some values of @xmath277 for which it is necessary to perform an integration by parts in order to isolate these vanishing integrals . once again , this is allowed because the @xmath304 s do not contain divergences . at the next order in @xmath194 , we have the last contribution to @xmath192 : @xmath307 after performing the integrations in @xmath68 and @xmath67 with the help of eq . ( [ intbas ] ) , eq . ( [ lcke ] ) becomes : @xmath308 the integral in @xmath115 is divergent and needs a regularization . going to polar coordinates , we obtain the result : @xmath309 also in this case , we note that the presence of the regularization modifies the dependence of @xmath310 on the pseudo - energy @xmath71 with respect to the naive formula of eq . ( [ behaone ] ) . the correction consists in the factor @xmath311 . in the limit @xmath302 , this factor diverges as powers of @xmath240 .    to conclude the analysis of the contribution to @xmath182 at the second order in @xmath194 , we show that the @xmath312 s are identically equal to zero . as a matter of fact , it is easy to verify that for @xmath313 each @xmath312 contains terms of the kind : @xmath314 these terms vanish identically because of the following identity : @xmath315 substituting eq . ( [ auxfskj ] ) in eq . ( [ blmn ] ) and using the explicit expression of the vector potential @xmath283 of eq . ( [ vecpot ] ) , we get : @xmath316 clearly , the right hand side of the above equation is zero because @xmath317 . if @xmath278 , instead , the vanishing function @xmath318 of eq . ( [ blmn ] ) may be isolated in the expression of @xmath319 only after an integration by parts . finally , at the third order in @xmath194 we have only the quantities @xmath320 s , since @xmath192 has at most quadratic powers of @xmath194 . it is easy to realize that : @xmath321 because the following relations hold : s is separately equal to zero , because these quantities contain terms of the kind given in eq . ( [ blmn ] ) . ] @xmath322 as a consequence of eq . ( [ threezero ] ) , it is clear that there are no contributions to @xmath182 at this order .    using eqs . ( [ kzerfin ] ) , ( [ konefin ] ) and ( [ ktwofin ] ) , we arrive at the final result for @xmath182 : @xmath323 we can now insert in the formula of the second moment of eq . ( [ secmomaveerep ] ) the functions @xmath183 and @xmath182 given in eqs . ( [ deeefin ] ) and ( [ nefin ] ) respectively . the outcome is : @xmath324 } { s+\frac \pi c\lambda(e)e^{-1 } } \label{mefin}\ ] ] in the @xmath9space , the already mentioned difficulties with the computation of the inverse laplace transform of @xmath183 and @xmath182 allow an analytical result only in the double limit of weak coupling constant @xmath15 and of large values of @xmath8 . at the first order in @xmath15 , the expression of @xmath61 reads as follows : @xmath325 where @xmath326 has been given in eq . ( [ ddd ] ) . so far , we have considered the averaged second moment of eq . ( [ secmomerep ] ) , corresponding to the case in which the polymer ends are not fixed . in the energy representation , we have seen that this version of the second moment can be exactly computed . to conclude this section , we would like to show that it is possible to provide also an exact expression of the second moment @xmath327 in the @xmath9 space and with fixed polymer ends up to an inverse laplace transform of the propagator given in eq . ( [ gffin ] ) . the starting point is the exact formula of the second moment @xmath328 of eq . ( [ feesm ] ) . all the ingredients of this formula are defined in eqs . ( [ secmomerep ] ) , ( [ nsdoqop])([nsdtqop ] ) and ( [ ampzero])([amptwofour ] ) . looking at eq . ( [ feesm ] ) , it is clear that : @xmath329 and @xmath330 let us note that the functions @xmath179 are all equal up to integrations by parts , which can shift the differential operators @xmath331 and @xmath332 in eqs . ( [ amponek])([amptwofour ] ) . this fact will be used in order to simplify the expression of the inverse laplace transformed of @xmath80 in the @xmath9space . to compute the inverse laplace transforms of both @xmath80 and @xmath81 , we use the following property of the inverse laplace transform of the product of two functions @xmath333 and @xmath334 : @xmath335 applying eq .  ( [ convilp ] ) to evaluate the inverse laplace transforms of @xmath192 and of the @xmath179 in eqs . ( [ nefull ] ) and ( [ defull ] ) , we obtain after some calculations : @xmath336 the second term in the right hand side of eq . ( [ nlfull ] ) is the contribution given by the functions @xmath179 , @xmath180 , while the first term comes from @xmath178 . remembering the definition ( [ secmomnumden ] ) of the second moment in terms of @xmath74 and @xmath75 , we get : @xmath337^{-1}\left[\frac 2c \int d^2x { { \mbox{\bf a}}}^2({{\mbox{\bf x}}})\int_0^l ds g(l - s;{{\mbox{\bf r}}}_1,{{\mbox{\bf x}}})g(s;{{\mbox{\bf x}}},{{\mbox{\bf r}}}_0)\right.\nonumber\\ & & \!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\ ! \!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\ ! \left.-\frac{2}{c^2}\int d^2x\int d^2 y\int_0^lds g(l - s;{{\mbox{\bf x}}},{{\mbox{\bf r}}}_1)\int_0^s ds ' \partial_x^i\partial_y^jg(s - s';{{\mbox{\bf y}}},{{\mbox{\bf x } } } ) g(s';{{\mbox{\bf r}}}_0,{{\mbox{\bf y } } } ) a_i({{\mbox{\bf x}}})a_j({{\mbox{\bf y}}})\right]\label{formtl}\end{aligned}\ ] ] if we knew how to compute the propagator @xmath132 starting from its laplace transformed ( [ gffin ] ) , we could evaluate explicitly the expression of the second moment in the @xmath9space . unfortunately , it is too complicated to perform the inverse laplace transform of the propagator @xmath107 . due to this technical difficulty , eq . ( [ formtl ] ) is only formal . progress can be made however in the limit @xmath191 , in which the propagator is given by the green function @xmath144 of eq . ( [ gzspec ] ) . this will be done in the next section . in order to allow the comparison with previous results , this section is dedicated to the case of ideal chains in which @xmath191 . first of all , we discuss the formula of the averaged second moment derived in the previous section , eq . ( [ mlfin ] ) . in the limit @xmath191 , ( [ mlfin ] ) becomes : @xmath338 the presence of a geometrical factor like the surface @xmath253 of the system in the expression of @xmath339 has been already related to the translational symmetry of eqs . ( [ strainvzer ] ) and ( [ strainvone ] ) . assuming that this surface has approximately the shape of a disc of radius @xmath268 , we can put @xmath340 as in eq . ( [ polinfultint ] ) . ( [ mlzerofin ] ) predicts that the average degree of entanglement scales as follows with respect to the distance @xmath268 : @xmath341 the meaning of eq . ( [ scalerisr ] ) is the following . we remember that the averaged second moment @xmath339 describes the entanglement of two closed polymers whose ends on the surfaces at @xmath37 and @xmath38 are not fixed . in this way , the polymers are allowed to move freely and it is natural to suppose that , the bigger will be the volume @xmath342 in which the polymers fluctuate , the bigger will be also the average distance between them . thus , if the surface @xmath253 increases its area , the probability of entanglement must decrease . the exact law of this decreasing is given by eq . ( [ scalerisr ] ) .    on the other side , one would expect that the probability of getting entangled is higher for long polymers than for short polymers . ( [ mlzerofin ] ) gives a result which is in agreement with the above expectation , because the second moment @xmath343 scales as follows with respect to the parameters @xmath8 and @xmath17 , which determine the polymer length : @xmath344 in particular , one can show that the total length of a polymer increases proportionally to @xmath8 and it is inversely proportional to the square root of @xmath17 @xcite . accordingly , we see from eq . ( [ scalereslc ] ) that @xmath343 increases proportionally to @xmath8 and inversely proportional to @xmath17 .    at this point we wish to study the second moment @xmath345 of polymers with fixed endpoints . the subscript @xmath346 has been added to the symbol of the second moment to remember that we are working in the limit @xmath191 . since we are dealing with ideal polymers , we have to substitute everywhere in eq . ( [ formtl ] ) the full propagator @xmath132 with the free one . the result of this operation is : @xmath347^{-1}\left[\frac 2c \int d^2x { { \mbox{\bf a}}}^2({{\mbox{\bf x}}})\int_0^l ds g_0(l - s;{{\mbox{\bf r}}}_1,{{\mbox{\bf x}}})g_0(s;{{\mbox{\bf x}}},{{\mbox{\bf r}}}_0)-\right.\nonumber\\ & & \!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\ ! \!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\left.\frac{2}{c^2}\int d^2x\int d^2 y\int_0^lds g_0(l - s;{{\mbox{\bf x}}},{{\mbox{\bf r}}}_1)\int_0^s ds ' \partial_x^i\partial_y^jg_0(s - s';{{\mbox{\bf y}}},{{\mbox{\bf x } } } ) g_0(s';{{\mbox{\bf r}}}_0,{{\mbox{\bf y } } } ) a_i({{\mbox{\bf x}}})a_j({{\mbox{\bf y}}})\right]\label{mlfullvez}\end{aligned}\ ] ] we notice that , as it could be expected , eq . ( [ mlfullvez ] ) coincides with the expression obtained in @xcite for the second moment of one polymer winding up around an infinitely long straight wire lying along the @xmath5axis . luckily , the propagator @xmath348 can be explicitly constructed upon computing the inverse laplace transform of the propagator @xmath349 of eq . ( [ gzspec ] ) : @xmath350 it is easy to check that the second term in the right hand side of eq . ( [ mlfullvez ] ) , which is associated with the contributions coming from the @xmath179 s , does not grow with increasing values of @xmath8 . as a matter of fact , after a rescaling of variables , the numerator of this term gives : @xmath351\nonumber\\ & & \frac 1{t ' } e^{-\frac c{2(t')}\left|{\mbox{\mbld{y}}}'-\frac{\mathbf r_0}l\right|^2 } \label{mlsmfullvez}\end{aligned}\ ] ] in the limit @xmath352 , the quantity in the right hand side of the above equation scales as @xmath353 , where @xmath10 is a constant . moreover , the propagator ( [ freelprop ] ) , which is in the denominator , scales as @xmath354 . thus , the ratio between the right hand side of eq . ( [ mlsmfullvez ] ) and the propagator ( [ freelprop ] ) does not depend on @xmath8 . this completes the proof of our statement .    as a consequence of this statement , as far as the scaling of @xmath345 for large values of @xmath8 is concerned , it is possible to make the following approximation : @xmath355^{-1 } \int d^2x { { \mbox{\bf a}}}^2({{\mbox{\bf x}}})\int_0^l ds g_0(l - s;{{\mbox{\bf r}}}_1,{{\mbox{\bf x}}})g_0(s;{{\mbox{\bf x}}},{{\mbox{\bf r}}}_0 ) \label{star}\ ] ] unfortunately , despite of the fact that we are treating ideal polymers , the integral in @xmath356 appearing in the above equation is still complicated and requires some approximation to be evaluated analytically . we will apply to this purpose the strategy used in ref . @xcite to compute the second moment of three dimensional polymers , adapting it to our two - dimensional case . first of all , let us note that the integral in ( [ star ] ) is ultraviolet divergent . however , the infrared divergences which appeared in the energy representation are absent . this is due to the behavior of the propagator @xmath126 , which is much milder at infinity than the behavior of the green function @xmath144 . to regulate the singularities at small distances , we proceed as usual by introducing the cut - off @xmath150 . after a rescaling of all variables similar to that of eq . ( [ mlsmfullvez ] ) , we get : @xmath357^{-1}\times \nonumber\\ & & \!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\ ! \int_{|x'|\ge\frac{a\sqrt c}{\sqrt l } } \frac{d^2x'}l \frac 1{\mathbf x'^2 }   \int_0 ^ 1 \frac{ds'}{s'(1-s ' ) } \left ( \frac c{2\pi } \right)^2 e^{-\frac 1{2(1-s')}\left|\mathbf x'-\mathbf r_1\sqrt{\frac      cl}\right|^2 } e^{-\frac 1{2s'}\left|\mathbf x'-\mathbf r_0\sqrt{\frac cl}\right|^2 } \label{starst}\end{aligned}\ ] ] to go further , following @xcite , we assume that the relevant contribution to the integral in @xmath358 comes from a narrow region around the singularity in @xmath359 . thus , we may put @xmath360 after making the above crude approximation , we obtain : @xmath361^{-1}\log\left(\sqrt{\frac    lc}a\right ) \int_0 ^ 1ds'\left [ \frac 1{1-s ' } e^{-\frac 1{2(1-s')}\mathbf r_1 ^ 2\frac      cl } e^{-\frac 1{2s'}\mathbf r_0 ^ 2\frac      cl}+ \right.\nonumber\\ & & \left . \frac 1{s ' } e^{-\frac 1{2(1-s')}\mathbf r_1 ^ 2\frac      cl } e^{-\frac 1{2s'}\mathbf r_0 ^ 2\frac      cl}\right ] \label{starstar}\end{aligned}\ ] ] in deriving the above equation we have used the simple relation @xmath362 . let us now study the integral @xmath363 the other integral in @xmath364 appearing in ( [ starstar ] ) can be treated in the same way after the change of variables @xmath365 . it is not to allowed to take in the right hand side of eq . ( [ inttrian ] ) the limit @xmath352 because in this way the integral will not be convergent due to the singularity in @xmath366 . for this reason , we split the domain of integration as follows : @xmath367 where @xmath368 . clearly , the second integral converges after performing the limit @xmath352 in the integrand and gives : @xmath369 the first integral instead diverges logarithmically with growing values of @xmath8 . however , now it is possible to expand the exponential @xmath370 in powers of its argument , because the singularity in @xmath371 lies outside the interval @xmath372 $ ] . keeping only the leading order term with respect to @xmath8 , we get : @xmath373 where ei@xmath374 is the exponential - integral function . when @xmath8 is large , this function may be approximated as follows : ei@xmath375 and , as a consequence : @xmath376 the second integral which we have left in eq . ( [ starstar ] ) gives the same result . putting everything together in the expression of the second moment of eq . ( [ starstar ] ) , we obtain the final result : @xmath377 this is exactly the behavior of the second moment derived in ref .  @xcite . in this article we have studied the entanglement of two directed polymers from a non - perturbative point of view . our formulas of the second moment , a quantity which measures the degree of entanglement of the two polymers , take into account the repulsive forces acting on the segments of the polymers and are exact . the averaged second moment defined in eq . ( [ aderr ] ) , a version of the second moment corresponding to the situation in which the end points of the polymers are free to move , has been computed in eq . ( [ mefin ] ) as a function of the chemical potential @xmath71 conjugated to the distance @xmath8 between the end points in the @xmath35direction . the case of free ends is relevant in the treatment of nematic polymers and polymers in a nematic solvent @xcite . let us note that also the expression of the second moment without any averaging and in the @xmath8 space can be computed . this has been done in eq . ( [ formtl ] ) . however , this equation is explicit only up to the inverse laplace transform of the propagator ( [ gffin ] ) , which is too hard to be obtained in closed form . eq .  ( [ mefin ] ) shows that the averaged second moment is of the form @xmath378 . the overall factor @xmath379 coincides with the scaling power law of two ideal polymers . the correction @xmath333 to this fundamental behavior due to the repulsive interactions is a complicated function of @xmath71 , whose analysis would require numerical methods . nevertheless , it is possible to identify a dominance of the repulsive interactions in the domain of parameters in which the condition @xmath380 is satisfied . this corresponds roughly speaking to the situation in which polymers are very long . in this region , the scaling laws with respect to the energy @xmath71 of the numerator and denominator appearing in the right hand side of eq . ( [ mefin ] ) are corrected by factors which are logarithmic powers of @xmath240 , see for instance eq . ( [ leaordaasyexp ] ) . one advantage of our approach is that it is easy to separate within the expression of the second moment the contributions of purely entropic origin which are typical of free polymers from the contributions coming from the presence of the @xmath0function potential in the polymer action . this is essentially due to the splitting ( [ appexplx ] ) of the propagator @xmath229 appearing in the amplitudes ( [ amponek][amptwofour ] ) . the component @xmath187 of the propagator coincides with the propagator of ideal polymers , while the component @xmath189 takes into account the effects of the interactions . thanks to the splitting ( [ appexplx ] ) , it has been possible to study the way in which the repulsive forces affect the average degree of entanglement of the two polymers . this has been done in section  [ sec : five ] . our results are in agreement with the intuition . the precise law with which the effects of the repulsive forces on the entanglement decrease when the distance between the trajectories increases is given by eq . ( [ appexplx ] ) . in section  [ sec : five ] it has been discussed also the strong coupling limit , which should be taken to recover the limit of excluded volume interactions . in our exact approach , it is not difficult to consider the case in which the coupling constant @xmath15 is large . for instance , the component @xmath189 of the propagator , which is responsible of the effects due to the repulsive interactions , has been given in the strong coupling limit in eq . ( [ scnontrainvpro ] ) . studying the form of this component assuming that polymers are very long , it has been argued that , at strong coupling , the major contributions to the winding angle coming from the repulsive interactions occur when the trajectories are very near to each other . many other qualitative and quantitative characteristics of the behavior of the two polymer system under consideration have been presented in section  [ sec : five ] . the case of ideal polymers , in which @xmath191 , has been discussed at the end of section [ sec : seven ] in order to make comparison with previous works . the scaling of the averaged second moment for large values of @xmath8 obtained in eq . ( [ scalereslc ] ) is in agreement with the results of @xcite , if one takes into account the fact that , after the averaging procedure of eq . ( [ aderr ] ) and the infrared regularization of eqs . ( [ deezero ] ) and ( [ polinfultint ] ) . one is effectively treating a system of polymers confined in a cylinder of finite volume @xmath342 . in section [ sec : seven ] we have evaluated the second moment , always of two ideal polymers , using the approach of ref . the outcome of this calculation , namely the scaling behavior of @xmath381 at the leading order in @xmath8 , is reported in eq . ( [ lllsecmomfreefin ] ) . this result is in agreement with the square logarithmic behavior obtained in @xcite , but not with the logarithmic behavior predicted in @xcite . however , this discrepancy can be expected due to the fact that , in section [ sec : seven ] , we have assumed , following ref . @xcite , that the most relevant contribution to the second moment coming from the integral in eq . ( [ starst ] ) is concentrated in a narrow region near the singularity in @xmath359 . this clashes with the assumptions of ref . @xcite , in which instead it is argued that the main increase in the winding angle does not occur when the polymer trajectories are near , but rather when they are far one from the other . finally , there is also an apparent discrepancy between the linear scaling with respect to @xmath8 of the averaged second moment @xmath339 and the square logarithmic scaling of the second moment @xmath345 . this disagreement is explained by the fact that , in the first case , the ends of the polymers are free to fluctuate , while in the second case they are fixed . it is therefore licit to expect that two polymers with free ends are more likely to entangle than two polymers whose ends are constrained . concluding , we would like to discuss possible further developments of this work , together with some problems which are still left open . first of all , the number of entangling polymers has been limited to two . to go beyond this restriction , one should explore the possibility of replacing the external vector potential @xmath283 of eq . ( [ vecpot ] ) with chern  simons fields . abelian chern - simons field theories have been already successfully applied in order to impose topological constraints to the trajectories of an arbitrary number of closed polymer rings in @xcite . we hope to extend those results also to the case of directed polymers in a forthcoming publication . of course , if the polymer trajectories are open , the constraints among them are no longer of topological nature as in @xcite , so that the application of chern - simons field theory to directed polymers should be considered with some care .    we have also not made any attempt to introduce in the treatment of polymer entanglement more sophisticated constraints than those which can be imposed with the help of the winding angle . this is in effect still an unsolved problem , despite the fact that two powerful and strategies have been proposed for its solution @xcite . in the first approach , pioneered independently by kleinert , kholodenko and one of the authors @xcite . the constraints are expressed via the wilson loop amplitudes of non - abelian chern - simons field theories . some progresses toward a concrete realization of this program in polymer physics have been made in refs . @xcite . in the second approach , developed by nechaev and coworkers , see @xcite and references therein . polymer trajectories are mapped on a complex plane with punctures . the link invariants necessary to impose the constraints are then constructed using the properties of conformal maps . another possible development is the treatment of attractive interactions , in which the strength @xmath15 in eq . ( [ dirpot ] ) takes negative values . in this case , the @xmath0function potential admits a bound state @xcite and the propagator of eq . ( [ gffin ] ) develops a singularity , in which @xmath382 , at the energy corresponding to this bound state . it would be extremely interesting to investigate how these facts affect the polymers entanglement . another issue which deserves attention is that of hairpin turns . hairpins are important in nematic solvents @xcite and can be included with the help of field theories @xcite . we note also that in our formalism it is also possible to study the entanglement of polymers in confined geometries . for example , values of @xmath71 which are near to @xmath383 ( @xmath384 ) correspond roughly speaking to the situation in which polymers fluctuate in a quasi two - dimensional environment , in which the height in the @xmath35direction is of the order of a few molecular sizes .    finally , an open problem , which has not be discussed here because we were mainly interested in the second moment . is the derivation of the full partition function @xmath64 of eq . ( [ laplrep ] ) . as anticipated in the introduction , it is not an easy task to compute @xmath64 because the repulsive potential of eq . ( [ dirpot ] ) is not central . we note however that the expression of @xmath64 coincides with the green function of a spin @xmath385 aharonov - bohm problem in the imaginary time formulation of quantum mechanics . this green function has been already derived in @xcite using sophisticated techniques developed in refs . @xcite , which bypass all the difficulties of dealing with a non - central potential . thus , in principle , the expression of the partition function @xmath64 is known . unfortunately , some of the consistency conditions imposed on the parameters in ref . @xcite seem to be incompatible with the requirements of our physical problem , as noted in section [ sec : three ] . for these reasons , the computation of the full partition function @xmath64 is still a problem which needs further investigations . luckily , the knowledge of the partition function is not necessary if one is interested to study the excluded volume interactions , which arise in the strong coupling limit . in fact , in this case it is possible to apply a powerful method due to kleinert @xcite . this method turns the weak coupling expansion into a strong coupling expansion which is convergent for large values of @xmath15 and is able to accommodate also the anomalous dimensions of quantum field theories . the convergence of this strong coupling expansion is mostly very fast , so that only a few coefficients of the weak coupling expansion must be known , see refs . @xcite for more details . these coefficients can be easily computed starting from the well known partition function of the aharonov - bohm problem without the insertion of the @xmath0function potential @xcite and treating this potential as a small perturbation assuming that the value of @xmath15 is small . the application of kleinert s method in order to complete the brief analysis of the strong coupling limit made in this paper is work in progress . 99 s. a. wasserman and n. r. cozzarelli , _ science _ * 232 * ( 1986 ) , 951 ; s. d. levene , c. donahue , t. c. boles and n. r. cozzarelli , _ biophys . * 69 * ( 1995 ) , 1036 ; f. livolant and y. bouligand , * j. phys . * ( paris ) * 47 * ( 1986 ) , 1813 ; r. podgornik , d. c. rau and v. a. parsegian , _ macromolecules _ * 22 * ( 1989 ) , 1780 . s. chandrasekhar , b. k. sadashiva and k. a. suresh , _ pramana _ * 9 * ( 1977 ) , 471 ; nguyen huu tinh , h. gasparoux and c. destrade , _ mol . _ * 68 * ( 1981 ) , 101 ; t. k. attwood , j. e. lyndon and f. jones , _ liq . crys . _ * 1 * ( 1986 ) , 499 . c. blisle , _ ann . * 17 * ( 1989 ) , 1377 . j. w. pitman and m. yor , _ ann . * 14 * ( 1986 ) , 733 . a. comtet , j. desbois and c. monthus , _ jour . * 73 * ( 1993 ) , 433 . a. grosberg and s. nechaev , _ polymer topology _ , advances in polymer science * 106 * , ( springer verlag , 1993 ) , 1 . b. drossel and m. kardar , in _ theoretical and mathematical models in polymer research _ , a. grosberg ( editor ) , ( academic press , 1998 ) , 187 ; cond - mat/9610019 ; _ phys . rev . _ * e53 * ( 1996 ) , 5861 . f. tanaka , _ prog . * 68 * ( 1982 ) , 148 . d. r. nelson , in _ observation , prediction and simulation of phase transitions in complex fluids _ , m. baus , l. f. rull and j. p. ryckaert ( eds . ) , ( kluwer , the netherlands , 1995 ) . j. d. moroz and r. d. kamien , _ nucl . _ * b506 * [ fs ] ( 1997 ) , 695 . a. l. kholodenko and t. a. vilgis , _ phys . _ * 298 * ( 1998 ) , 251 . f.  ferrari and i.  lazzizzera , phys . b * 444 * ( 1998 ) 167 ; f. ferrari , h. kleinert and i. lazzizzera , _ eur . phys . j. _ b * 18 * ( 2000 ) , 645 . m. otto , _ j. phys . a : math . * 34 * ( 2001 ) , 2539 ; _ topological interactions in systems of mutually interlinked polymer rings _ , cond - mat/0402104 . m. g. brereton and s. shah , _ a _ : math . * 15 * ( 1982 ) , 989 . s. edwards , _ proc . soc . _ * 91 * ( 1967 ) , 513 ; _ a_1 ( 1968 ) , 15 . f. spitzer , _ trans . _ * 87 * ( 1958 ) , 187 . s. prager and h. l. frisch , _ j. chem . * 46 * ( 1967 ) , 1475 . n. saito and y. chen , _ * 59 * ( 1973 ) , 3701 . h. kleinert , _ path integrals in quantum mechanics , statistics , polymer physics , and financial markets _ , ( world scientific publishing , 3nd ed . , singapore , 2003 ) . m. kardar and y .- c . zhang , _ phys . * 58 * ( 1987 ) . 2087 ; m. kardar , g. parisi and y .- c . zhang , _ phys . * 56 * ( 1986 ) , 889 . r. d. kamien , p. le doussal and d. r. nelson , _ phys . rev . _ * a 45 * ( 12 ) ( 1992 ) , 8727 . k. iwata and t. kimura , _ phys . _ * 74 * ( 1981 ) , 2039 . h. bethe and r. peierls , _ proc . ( london ) * 148a * ( 1935 ) , 146 ; e. fermi , _ ricerca scientifica _ * 7 * ( 1936 ) , 13 , [ english translation : e. fermi , _ collected papers _ vol italy 19211938 ( university of chicago press , chicago , 1962 ) ] . f. berezin and l. faddeev , _ sov . * 2 * ( 1961 ) , 372 . s. albeverio , f. geztesy , r. hoeg - krohn and h. holden , _ solvable models in quantum mechanics _ ( springer verlag , berlin , 1988 ) ; s. albeverio , z. brzeniak andl . dabrowski , _ jour . * 27 * ( 1994 ) , 4933 . p. gerbert , _ phys . * 40 * ( 1989 ) , 1346 . c. r. hagen , _ phys . * 64 * ( 1990 ) , 503 ; _ int . _ a*4 * ( 1991 ) , 3119 . c. grosche , _ phys . * 71 * ( 1993 ) , 1 ; _ ann . phys . _ * 3 * ( 1994 ) , 283 . d. k. park , _ jour . * 36 * ( 10 ) ( 1995 ) , 5453 . r. p. feynman and a. r. hibbs , _ quantum mechanics and path integrals _ , ( mcgraw - hill , new york , 1965 ) . l. s. schulman , _ techniques and applications of path integrals _ , ( wiley , new york , 1981 ) . r. jackiw , _ delta  function potential in two- and three - dimensional quantum mechanics _ , in m. a. bg memorial volume , a. ali and p. hoodbhoy ( eds . ) , ( worlds scientific , singapore , 1991 ) . h. kleinert , _ phys . rev . _ * d57 * ( 1998 ) , 2264 . h. kleinert , _ phys . lett . _ * b434 * ( 1998 ) , 74 , [ arxiv : cond - mat/9801167 ] . h. kleinert and v. schulte - frohlinde , _ critical properties of @xmath386theories _ , world scientific , singapore 2001 . f. ferrari , _ phys . lett . _ * a323 * ( 2004 ) , 351 , cond - mat/0401104 . f. ferrari , h. kleinert and i. lazzizzera , int .  j.  mod .  phys . b*14 * ( 2000 ) 3881 [ arxiv : cond - mat/0005300 ] . h. kleinert , _ pfadintegrale in quantenmechanik , statistik un polymerphysik _ , spektrum akademischer verlag , heidelberg , 1st edition 1993 , 850p . , isbn 3 - 86025 - 613 - 0 . a. kholodenko , _ trends in chem . 3 , research trends inc . trivandrum , india ( 1994 ) , pp . s. nechaev and r. voituriez , _ conformal geometry and invariants of 3-strand brownian braids _ , [ arxiv : hep - th/0407008 ] . f. ferrari , _ chern - simons field theories with non - semisimple gauge group of symmetry _ _ jour . math . * 44 * ( 1 ) ( 2003 ) , 138 , hep - th/0210100 . f. ferrari , _ topological field theories with non - semisimple gauge group of symmetry and engineering of topological invariants _ , published in _ trends in field theory research _ , o. kovras ( ed . ) , nova science publishers ( 2004 ) . j. cardy , _ jour . phys . a _ * 16 * ( 1983 ) , l355 .
Abstract:

Answer: in this article we study from a non - perturbative point of view the entanglement of two directed polymers subjected to repulsive interactions given by a dirac @xmath0function potential . an exact formula of the so - called second moment of the winding angle is derived . this result is used to provide a thorough analysis of entanglement phenomena in the classical system of two polymers subjected to repulsive interactions and related problems . no approximation is made in treating the constraint on the winding angle and the repulsive forces . in particular , we investigate how repulsive forces influence the entanglement degree of the two - polymer system . in the limit of ideal polymers , in which the interactions are switched off , we show that our results are in agreement with those of previous works .