Arxiv:Cond-Mat/9804050V1 [Cond-Mat.Soft] 5 Apr 1998 Ino H N-Oedvector End-To-End Distribu- the Probability the of Representation, Tion Spectral and Surface

Home , Polymer physics

arXiv:cond-mat/9804050v1 [cond-mat.soft] 5 Apr 1998 ino h n-oedvector end-to-end distribu- the probability the of representation, tion spectral and surface. connec- fixed its problem curved, by the a only solve characterized on to chain, tivity, necessary polymer is free it a effort, of embark- an Before such on theory. ing group tech- renormalization variational or analysis, niques circum- perturbation these via gener- obtained Under be ally can results for. analytical in- approximate accounted curved only be stances, of to fluctuations need and terfaces segments polymer tween super- motion. II Brownian Type rotational in the lines and (e.g. vortex [3] free of other spaces physics conductors of in connected the mechanics [2]), arise multiply quantum Ref. issues and see the of curved mem- well: types in as the These particles physics on in [1]. of captured cell topology area biomolecule a and of a geometry brane by of wrap played dynamics vesi- role polymers the and the way macromolecules and rodlike the cles, bi- and to in spherical significance related around of problems are interfaces ological curved in embedded neatn)Gusa oye yn nacre man- curved a on lying ifold. polymer Gaussian interacting) ieepesdb h ensur n-oeddistance end-to-end mean-square the by size—expressed eateto hsc srnm,adPlmrSine&Eng & Science Polymer and Astronomy, & Physics of Department ( R ntepeetpprw eie sn path-integral using derive, we paper present the In be- interaction excluded-volume situations, realistic In polymers of characteristics topological and Statistics RWINMTO N OYE TTSISO ETI CURVED CERTAIN ON STATISTICS POLYMER AND MOTION BROWNIAN − R ′ ) 2

n-oeddsac sjs h asinvalue Gaussian the just is distance end-to-end emtia oaiainarea localization geometrical distance peecs,i slna in linear is it case, sphere al n ueial uigtedﬀso qainadtep the and equation G( diﬀusion function the distribution (using numerically and cally nlsdo ufc meddi h ulda pc.We t When space. Euclidean the in embedded D surface a on enclosed rpriso h ufc;a ag cls( scales large at surface; the of properties n h n position end the and sidcda l eghsae ytepeec ftevertex, the of presence case the the by In scales torus. length curved all the at for induced surface) is the of geometry the

( nteabetsaeo niel(non- ideal an space—of ambient the in R iesosadteclne,tecn n h uvdtrsin torus curved the and cone the cylinder, the and dimensions nti ae ehv osdrdaGusa oye hi fl of chain polymer Gaussian a considered have we paper this In − .INTRODUCTION I. R ′

) ( 2 R

. − R ′ ) 2

fteplmrcan u nig r:tecraueo h su the of curvature the are: findings Our chain. polymer the of R G ′ R nvriyo ascuet,Ahrt A01003 MA Amherst, Massachusetts, of University ftecan xlctcosvrfrua r eie o G( for derived are formulas crossover Explicit chain. the of ( | L R R auP Mondescu P. Radu tsotsae ( scales short at ; ′ o h yidradrahsdffrn osatvle a f a (as values constant different reaches and cylinder the for | ; R L ′ fteedt-n vector end-to-end the of ) ; L n the and ) MANIFOLDS niern Center, Engineering L L ∞ → → ∗ Ll 1 n .Muthukumar M. and )teplmris polymer the 0) ), ( sphere der L oecue oueo te ye fitrcin r ac- paper. are this interactions in of for types mentioned, other counted As or volume probability. excluded this no by obeyed equation sion beti h peia ufc,wiei s2–dimensional is surfaces. it other while the surface, in spherical the in object eetto 45 ftepoaiiydsrbto function distribution probability G( the of [4,5] resentation yet. here answered address been we of not problems have theory specific the the from available motion, already Brownian are chain Gaussian results ideal many an [4], of statistics the and problem walk Fg D in 1D) (Fig. nFg .Teplmri nintrinsic an is polymer The 1. Fig. in n at one hi with chain in ie sa xlctpt nerl[5] integral path explicit func- an as distribution given probability tion, following the to according l

( neigDprmn n aeil eerhSineand Science Research Materials and Department ineering h uvdsrae tde nti ae r h S the are paper this in studied surfaces curved The norcluain ehv sdtept-nerlrep- path-integral the used have we calculations our In random- the between equivalence the to due Although, e scnie a consider us Let unlnt) needn ftemetric the of independent length), Kuhn = = R R Fg B,tecircular the 1B), (Fig. | Nl − R R in R ′ R ; ihtefis edlctdat located bead first the with , − t-nerlapoc)teprobability the approach) ath-integral safnto fteoeigage2 angle opening the of function a as L ′ h n-oedvector end-to-end The . D ) ftecn,cnrcino h chain the of contraction cone, the of faGusa oye hi n h diffu- the and chain polymer Gaussian a of ) 2 R

esrae r h peeS sphere the are surfaces he N ′ iesos(i.1,for 1A, (Fig. dimensions R ed oacntn au nthe in value constant a to tends n h ensur end-to-end mean-square the and R oal ﬂat locally 3 ik fKuhn-length of links ength h emtia aaeesaeshown are parameters geometrical The . 3 ehv acltdanalyti- calculated have we , I FORMULATION II. L D dmninlGusa polymer Gaussian –dimensional sa nrni object intrinsic an as n h mean-square the and cone fcsidcsa induces rfaces R Fg C n the and 1C) (Fig. | R nto of unction R ′ ; − L D l D D and ) − n oa length total and R − 1 R ) the 3); = ′ in 1–dimensional α ′ sdistributed is n h last the and cylin- torus D − 1 DN r =R N 1 2 N − Here φ and E are obtained from the eigenvalue equa- R R D Dr µ µ G( ′; L) = lim d j tion | N 2πlδs r R′ 0= j=1 δs→∞0 Z Y Nδs→=L l ∆φµ(R)= Eµφµ(R), (2.8) N − 2D D 2 exp (rj rj 1) , (2.1) × −2lδs − −  j=1 X with µ some generalized index. In general, the probabil-   ity distribution G(R R′; L) is not translational invariant. | where D is the intrinsic dimension of the polymer, δs The mean-square end-to-end distance can be calculated is the change in the arc-length s along the chain and from r j is the position vector of the j-th bead of the polymer DR DR R R 2 R R 2 d d ′ ( ′) G( ′; L) chain. Note that the distance on the surface is calculated (R R′) = − | . r r 2 − dDR dDR G(R R ; L) as ( j j 1) = guv(uj uj 1) (vj vj 1), where u, v R ′ ′ − − − − − − | are the intrinsic coordinates and the metric tensor guv R (2.9) is discretized in a symmetric form or using the midpoint rule (e.g. see Ref. [ [6,7]]). In RD, the expression (2.1) gives the Gaussian probability distribution: It is to be noted that the calculated mean-square distance represents the distance between the ends of the polymer D 2 measured by an outside observer, located in the ambi- D D 2 G(R R′; L)= exp (R R′) , (2.2) ent space, which is diﬀerent from the geodesic separation | 2πLl −2Ll − measured along the surface.

The probability G(R R′; L) obeys the diﬀusion equation | III. CALCULATIONS AND RESULTS ∂ l ∆ G(R R′; L)=0, (2.3) ∂L − 2D | Due to the known equivalence [4,8] between the dif- subject to the initial condition fusion equation (2.3) and the Schr¨odinger equation for a free particle, we can readily obtain the polymer prob- (D) lim G(R R′; L)= δ (R R′). (2.4) ability distribution G(R, R′; L) for the sphere and the L 0 | − → cylinder cases, from the corresponding probability am- plitudes for a rigid rotor in D dimensions [9] and for a When the polymer lies in a curved manifold, one needs 2–D free particle moving on the surface of a cylinder, to employ the Laplace-Beltrami operator given by (e.g. respectively. see Ref. [ [6]])

1 ij ∆= ∂ig ∂j , (2.5) D 1 RD √g A. The Sphere S − in

ij where g is the inverse of the metric tensor gij of the Let a be the radius of the D 1 dimensional sphere surface and g = det(gij ). We remark that gij is the in- centered at the origin of the coordinate− system. Then RD duced metric from the ambient embedding space ( or R(s) is a vector of magnitude a that measures the posi- R3 ). tion of the bead at s about the center of the sphere. The One can obtain the solution of Eq. (2.3) satisfying the intrinsic dimension of the polymer is D 1. In spherical − initial condition (2.4) by applying a Laplace transform coordinates the Laplace-Beltrami operator (2.5) becomes with respect to L that gives the Green’s function equa- [7, 4.2.3] § tion 1 2 ∆= 2 , (3.1) l (D) a L E ∆ G(R R′; E)= δ (R R′), (2.6) − 2D | − where 2 is the D-dimensional Legendre operator, satis- e L which can be solved by an eigenfunction expansion of G fying the eigenvalue equation followed by an inverse Laplace transform that eventually 2Sm(Ω) = J(J + D 2)Sm(Ω). (3.2) yields the heat kernel: e L J − − J R R R R m G( ′; L)= φµ( )φµ∗ ( ′) exp( EµL). (2.7) The SJ functions are the hyperspherical harmonics or- | − D 1 µ thonormal and complete on the unit S − sphere [10]: X

2 ′ m m ν 1 ν dΩ S (Ω)S ′ (Ω) = δ ′ δ ′ (3.3) cos θ C (cos θ)= [(l + 1) C (cos θ) J J JJ mm l 2(l + ν) l+1 Z M ν ∞ + (2ν + l 1) Cl 1(cos θ)], (3.8) m m (D) − S (Ω)S ∗(Ω′)= δ (Ω Ω′), − J J − m=1 JX=0 X we notice that the orthogonality of the Gegenbauer poly- ν nomials (with C0 (cos θ) = 1) implies that all terms ex- (2J+D 2) (J+D 3)! cept J = 1 are zero. From (3.7) we finally obtain the where M = J−!(D 2)! − and Ω is the solid angle in D dimensions. Using− Eq. (3.2) and the eigenvalue ex- mean-square end-to-end vector for a polymer living on m the surface of a D 1–dimensional sphere as: pansion (2.7), where φµ SJ and Eµ J(J + D 2), − we readily find the probability→ distribution→ function− of a 2 2 Ll polymer chain (or equivalently, of a random walk) that (R R′) =2a 1 exp − − −2a2 starts at R′ and ends at R in N = L/l steps, on the 2 SD 1 sphere: Ll ; Ll a . − ≪ (3.9) ≃ 2a2; Ll a2. M ≫ 1 ∞ m m G(R R′; L)= S (Ω) S ∗(Ω′) (3.4) 2 | aD 1 J J Therefore, a characteristic area 2a called geometrical lo- − J=0 m=1 R R 2 X X calization area emerges in the description of ( ′) , Ll J(J + D 2) − exp − . due to the geometrical restriction of space. × −2a2 D 1 − In contrast to the problem of a stiff polymer chain without any confinement, where the chain expands due The probability found above is normalized properly to to chain stiffness, chain contracts here by geometrical D R3 one ( d R G(R R′; L)=1). In we recover the known confinement. If the area of the confining sphere is small | [11] rotational diffusion result: in comparison with Ll, then the polymer wraps many R 2 times around the sphere so that (R R′) approaches l − 1 ∞ the geometrical localization area. On the other hand if G(R R′; L)= Y (Ω), Y∗ (Ω′) 2 | a2 lm lm a Ll, then the polymer does not sense the curva- l=0 m= l ≫ 2 ture of the sphere so that (R R′) = Ll. The mean- X X− − Ll square end-to-end distance (3.9) normalized at a2 is plot- exp J(J + 1) . (3.5) × −4a2 ted in Fig. 4 as a function of the dimensionless variable Ll R R 2 2 . It isalsoto be noted that ( ′) is independent a − R R of D. Because G( ′; L) is translationally invariant, the end- to-end distance| can be calculated as a conditional ex- pectation value by dropping the integration over R′ in 3 Eq. (2.9). To evaluate the remaining integral, we use the B. The Cylinder in R addition theorem for hyper-spherical harmonics [10]

M Consider a cylinder of radius a and of inﬁnite length, m m 1 (2J + D 2) R3 S (Ω) S ∗(Ω′)= − with the induced metric from (in cylindrical coordi- J J 2 SD D 2 nates, with Oz along the axis of the cylinder): gφφ = a ; m=1 − R X (D 2)/2 gφz = 0; gzz = 1; √g = a. The position vector (s) of C − (cosΨRR′ ), (3.6) × J any bead s has the components (a cos φ, a sin φ, z). The polymer dimension is D = 2. 2πD/2 ν where SD = Γ(D/2) and CJ are the Gegenbauer poly- The Green’s function equation in the Laplace space 2 nomials of argument cos ψRR′ . Writing (R R′) = (2.6) can be solved directly by expanding the G function 2 RR′ − 2a [1 cosΨ ] and inserting (3.4) in (2.9), we get: in the complete′ set of orthogonal functions on the circle − im(φ φ ) e − and Fourier transforming with respecte to z z′, − ∞ which yields 2 2 1 2J + D 2 (R R′) =2a 1 dΩ − − − SD D 2 + J=0 1 ∞ ∞ dk 1 XZ − R R (D 2)/2 G( ′; E)= 2 2 ′ ′ lm lk cosΨRR C − (cosΨRR ) (3.7) | 2πa 2π E + 2 + J m= Z−∞ 4a 4 × X−∞′ ′ Ll J(J + D 2) e im(φ φ ) ik(z z ) exp − . e − e − (3.10) × −2a2 D 1 × − After applying an inverse Laplace transform and eval- Choosing R′ along one of the coordinate axes and apply- uating the integral over k, the spectral expansion of ing the recurrence relation [12, 8.933] G(R R′; L) is found to be § |

3 ′ 1 1 (z−z )2 with the opening angle 2α. In cylindrical coordinates, R R Ll G( ′; L)= e− the position vector R(s) of any point s along the chain | 2πa rπLl is R(s) = (ρ(s)cos φ(s),ρ(s) sin φ(s),ρ(s)cot α) and the ∞ ′ Llm2 im(φ φ ) 2 2 e − e− 4a2 . (3.11) induced metric is: gρρ = 1+cot α; gρφ = 0; gφφ = ρ ; × m= g = ρ√1+cot2 α. As in the previous problem, the X−∞ √ dimension of the polymer is D = 2. An equivalent form of the expression above can be ob- R R tained by using the relation To compute the probability distribution G( ′; L) of the end-to-end vector we start with the Green’s| function + 2 1 ∞ 1 ξ2+i√2λξx equation (2.6) obeyed by the Laplace transform G: exp( λx )= dξ e− 2 (λ> 0), − √2π Z−∞ 2 l 2 ∂ ∂ l ∂ e (3.12) E sin α ρ G(R R′; E) − 4ρ ∂ρ ∂ρ − 4ρ2 ∂φ2 | sin α and the Poisson formula = δ(ρ ρ′) δ(φ φ′)e (3.16) + + ρ − − ∞ ∞ e2iπαm = δ(α n), (m,n Z). (3.13) − ∈ m= n= We solve this equation in a standard manner [13]. First X−∞ X−∞ we eliminate the azimuthal dependence of G by inserting Applying these transformations in Eq. (3.11), we arrive the expansion at the the winding number expansion e + ′ 2 2 ′ 1 ∞ 1 (z−z ) ∞ a (φ−φ +2πn) R R R R Ll Ll G( ′; E)= ˜g (ρ ρ′; E) G( ′; L)= e− e− , (3.14) | 2π m | | πLl m= n= −∞ X−∞ X e exp[im(φ φ′)], (3.17) × − where n (a topological term) is the winding number that counts how many times the chain winds around the cylin- in (3.16), which gives der. A term in the sum represents now the probability R R that a given chain starting at and ending at ′ will 1 ∂ ∂ ν2 wind around the cylinder n times (in L/l steps). This x (1 + ) ˜g (x x′; E) x ∂x ∂x − x2 ν | expansion can be used to investigate the statistics of the 4 winding number (e.g. the mean-square winding number = δ(x x′), (3.18) N 2 ). −xl sin α − To calculate the mean-square end-to-end vector, we 2 2 2 where we introduced the new variable x = kρ, with insert (R R′) = (z z′) +2a [1 cos(φ φ′)] and 2 − − − − 2 4E 2 m the distribution function (3.11) in (2.9). Evaluating the k = l sin2 α and ν = sin2 α . The solutions of the ho- integrals, we get mogeneous equation are the modified Bessel functions I ν (x) and K ν (x). Then, the regularity of ˜gν (x x′; E) Ll | | | | | 2 Ll 2 at x = 0, x , the continuity at x = x′ and the jump (R R′) = +2a 1 e− 4a2 − 2 − in the first derivative→ ∞ impose the solution Ll ; Ll h a2. i ≪ (3.15) ≃ Ll ; Ll a2. 4 2 ˜gν (ρ ρ′; E)= I ν (kρ<)K ν (kρ>). (3.19) ≫ | l sin α | | | | As expected, the polymer remains almost Gaussian (a prefactor of 1/2 instead of 1) along the Oz direction. Here ρ< and ρ> denote the smaller and the larger, re- The influence of the geometry appears through the ra- spectively, between the variables ρ and ρ′; we have also dius a, related to the mean curvature of the cylinder. used the fact that the Wronskian W of the modified 1 Once again, the polymer becomes geometrically localized. Bessel functions is W[Iν , Kν ](x)= . Inserting the so- − x The mean-square end-to-end distance (normalized at a2) lution found back in Eq. (3.17) and performing the in- is represented in Fig. 4 as a function of the dimensionless verse Laplace transform using the formula [14, 5.16.56] Ll § variable a2 . 1 − K [(√a + √b)√E]I [(√a √b)√E] = (3.20) L ν ν − n a+b o 3 1 a b C. The Cone in R e− 2L Iν ( − ) ( a, b> 0), 2L 2L ℜ ℜ

Now let us consider the surface of cone with the ver- we arrive at the exact spectral expansion of G(R R′; L) tex at the origin O, centered about the Oz axis and for a polymer on a cone: |

4 + 1 ∞ 2 Fig. 1). The torus is centered about the Oz axis and G(R R′; L)= I m (2ρρ′/(Ll sin α)) | πLl sin α sin α symmetric with respect to the xOy plane. We introduce m= | | X−∞ the toroidal coordinates η,θ,φ , [15, 2.13] where φ is 2 2 { } § ρ + ρ′ the usual azimuthal angle—in cylindrical coordinates— exp exp[im(φ φ′)], (3.21) × −Ll sin2 α − about Oz: c sinh η cos φ c sinh η sin φ R R R x = ; y = ; which is properly normalized to one ( d G( ′; L)= cosh η cos θ cosh η cos θ 1). For α = π we recover the expression in polar| coordi- − − 2 c sin θ nates of the propagator of a GaussianR polymer chain in z = , (3.24) cosh η cos θ the plane [7,8]. − As the total partition sum gives the area of the cone, where 0 η < , 0 θ < 2π and 0 φ < 2π. The R we calculate the end-to-end distance by fixing ′ and surface of≤ the torus∞ corresponds≤ to a fixed≤ value η of the R 0 integrating only over . Using cylindrical coordinates, coordinate η. Then we have the following geometrical R R 2 2 2 2 with ( ′) = (ρ + ρ′ )/ sin α +2ρρ′[1 cos(φ φ′)] relations: and applying− various formulae [12] involving− definite− in- c 2 2 2 tegrals of Bessel functions, we obtain b = c coth η0; a = ; c = b a . (3.25) sinh η0 − 2 2 2 (R R′) = Ll 1+2 t csc α − √π csc α (1 sin2 α) F ( 1/2, 1; t2 csc2 α) c2 1 1 The induced metric tensor is gθθ = 2 ; − − − − (cosh η0 cos θ) 2 2 − Γ(3/2+1/2csc α) (1 csc α) (1+csc α) c sinh η0 gφφ = 2 ; gθφ = gφθ = 0; g = 2 (sin α) − t (3.22) (cosh η0 cos θ) √ − Γ(1 + csc α) 2 − c sinh η0 2 2 2 . We recall that the length interval on the F [(csc α 1)/2, 1+csc α; t csc α] , (cosh η0 cos θ) 1 1 − 2 2 2 × − − surface is ds = gθθ dθ +gφφ dφ and the volume element 2r 1/2 is d = √g dθ dφ. where t = ρ′/(Ll) and 1F1 is the confluent hypergeo- metric function. The calculation of the propagator of a Gaussian chain on a torus in R3 (which is equivalent to solving the prob- Rather surprisingly, the mean-square distance depends lem of a free quantum particle moving on the same sur- on the position of the chain through the parameter t. To face) cannot be done in closed form, as involves find- interpret this finding, we calculate the limiting values of 2 ing the eigenvalues of a Hill-type equation. Still, using (R R′) for t (polymer far on the conical sur- − → ∞ the path-integral formalism and making certain “sloppy” face or α 0), t 0 (polymer with one end attached approximations, we were able to derive an approximate to the cone→ vertex)→ and α = π (polymer on a 2–D plane 2 form of G(R R′; L) that behaves physically correct when surface). One finds that | computing moments of the distribution of R R′, which − 2 π are the experimentally relevant quantities we are inter- (R R′) = Ll for t = 0; t ; α = , (3.23) − { → ∞ 2 } ested in. This approach will be justified a posteriori by showing that one recovers the usual flat surface case and which is just the Gaussian result in all cases. Physi- the expected asymptotic regimes obtained for spheres cally, when the polymer chain is either fixed at the origin and cylinders. (t = 0) or far away from it (t ), there are no other We start be recalling the path-representation (2.1) → ∞ R R 2 preferred points on the cone and ( ′) has the written in a compact form as planar value (we should also remark that− a cone is geometrically flat). In the crossover region, the singular rN =R N character of the cone vertex becomes dominant and in- G(R R′; L)= [r] exp Sj , (3.26) | r R′ D −  duces the contraction of the chain (possibly by winding 0= j=1 Z X about the origin). This behavior is illustrated in Fig. 2   R R 2 and Fig. 3, where the mean-square distance ( ′) where the integration measure is given explicitly by normalized at Ll is plotted as a function of the t−param- N N 1 eter and the angle α. − 1 2 [r] = lim d rj (3.27) D N πlδsj →∞ j=1 j=1 δsj 0 Y Y δs→=L D. The Torus in R3 j j P and Sj is the short-length polymer action Consider a torus embedded in R3, with the circular 2 cross-section of radius a. The axial circle containing the (rj rj 1) Sj = − − . (3.28) centers of the circular sections has the radius b>a (see lδsj

5 2πN 1 In the previous expressions, δsj = sj sj 1 (with − r − − dφj sinh η0 (3.33) j = 1,N, s0 = 0, sN = L) and j is the position vec- × 0 j=1 tor about the origin of the j-th element of the discretized Z Y N chain. c2 sinh2 η (∆φ )2. exp 0 j In order to evaluate the path-integral, due to the men- × − l σ  j=1 j tioned equivalence between the statistical mechanics of X polymers and quantum mechanics, we will apply the   method used to calculate the propagator of a free particle Observing that the integrals over the angular coordinates moving on a circle and in the presence of a ring-shaped are similar to the path integral for a free particle moving defect [6,8,16]. on a circle (e.g. see Ref. [ [6,7]]), we have the formula

We start with the action element Sj . Using the length N 1/2 2πN 1 (∆ψ )2 − N j interval on the torus, we write Sj in a symmetrized form α α j=1 σ dψj e− j πσ P c2 1 j=1 j Z0 j=1 Sj = Y Y ′ l δsj (cosh η0 cos θj )(cosh η0 cos θj 1) α ∞ α (ψ ψ +2πn)2 − − − = e− σ − . (3.34) 2 2 2 πσ [(∆θj ) + sinh η0 (∆φj ) ], (3.29) r n= × X−∞ where θj = θ(sj ), φj = φ(sj ) are the toroidal coor- { } Here α is an arbitrary positive constant and ψj = ψ(σj ) dinates of the element located at sj along the chain, is the discretized angular variable on the circle, with ψ′ ∆θj = θj θj 1 and ∆φj = φj φj 1. − − − − and ψ denoting the initial and the final positions. Equiva- Next we apply a local rescaling [7,17], which reads lently, recognizing that the angular path integral contains trajectories that wind a different number of times around σj = δsj (cosh η0 cos θj )(cosh η0 cos θj 1) (3.30) − − − the origin, one can use the covering space mapping [6,7]. and which satisfies the global scaling relation: Returning to Eq. (3.33) of the approximate discretized propagator GN , integrating over the angles using formula N (3.34), and replacing σ with its global value (3.31), we σ = σj = L(cosh η0 cos θ)(cosh η0 cos θ′). (3.31) finally obtain—up to a normalization factor—the approx- − − j=1 X imate distribution function of the end-to-end distance for a polymer chain on a curved torus in R3, expressed as a This constraint is nontrivial, as it must be compatible winding number expansion: with the local scaling (3.30). That it gives the correct answer when applied in the path integral and when one 1 ∞ ∞ c2 uses an anisometric discretization (unequal length inter- G(R R′; L) exp | ≈ πLl −Ll n= m= vals δsj ), was discussed by Inomata [17]. Basically, it X−∞ X−∞ amounts to the prescription [18] of ignoring terms of 1 1+ǫ (3.35) O(δsj ) in the path-integral, if ǫ > 0. Then, the dis- × (cosh η0 cos θ)(cosh η0 cos θ′) cretized functional measure becomes: − 2 2− 2 [(θ θ′ +2πn) + (sinh η )(φ φ′ +2πm) ] × − 0 − [r] = (cosh η cos θ)(cosh η cos θ′) D 0 − 0 − N N 1 − Another compact and illuminating form is obtained by 1 2 dθj dφj c sinh η0. (3.32) employing the transformations (3.12) and (3.13), which × σ j=1 j j=1 yield the angular momentum expansion Y Y

We apply now the length rescaling in the action Sj given R R 1 (cosh η0 cos θ)(cosh η0 cos θ′) G( ′; L) 2 − 2 − by (3.29) directly, by dropping boldly any terms coming | ≈ 4π c sinh η0 from a rigorous expansion of ∆θj and ∆φj about the new 2 Θ3[(θ θ′)/2, iσ(θ,θ′)l/(4πc )] (3.36) length variable σj . As mentioned, this will be justiﬁed × − 2 2 Θ [(φ φ′)/2, iσ(θ,θ′)l/(4πc sinh η )], by the ﬁnal results. Combining then with the integration × 3 − 0 measure, we get the approximate discretized propagator

GN where we recall that σ(θ,θ′)= L(cosh η0 cos θ) (cosh η0 2 2 2 b − cos θ′), c = b a and cosh η = .Θ is the theta 1 − − 0 a 3 GN (R R′; σ) (cosh η0 cos θ)(cosh η0 cos θ′) function deﬁned as [12] | ≈ c2 − − N 2πN 1 N c2 − c2 (∆θ )2 j ∞ iπτ n2 2izn dθj exp Θ3(z, τ)= e e ( τ > 0). (3.37) × πlσj 0 − l σj  ℑ j=1 Z j=1 j=1 n= Y Y X X−∞   6 First, we stress that this is not the exact Green func- where N(x, t) is given by tion of the heat equation on the curved torus. Still, as 2π necessary, the expression is symmetric in R, R′, transla- 1 tionally invariant with respect to the azimuthal angle φ, N(x, t)= dθ dθ′ 0 (t cos θ)(t cos θ′) and in the limit η0 a 0 reduces properly to ZZ − − →∞≡ → θ θ′ ix(t cos θ)(t cos θ′) the probability distribution of the end-to-end vector for Θ − , − − . (3.44) × 3 2 4π(t2 1) a polymer on a circle (see the angular part in Eq. (3.14)). − Also, it obeys the initial condition The partition function Z(x, t) of the polymer on the torus (2) G(R R′;0) = δ (R R′), (3.38) | − is: 2 2 3 and its trace has, at least numerically, the correct limit Z(x, t)= a (t 1) 2 N(x, t). (3.45) − [19] when Ll 0: a2 → This reduces properly to the area 4π2ab of the torus when Z(x 0,t)= d2R G(R,L R, 0) (3.39) x 1, t 1 (Θ(z, τ) 1): → | ≫ ≫ ≃ ZZ t π3 t (t √t2 1) 1 4π + x − − + o(x2), Z(x, t) 4π2 ab 1 4π2ab (3.46) ≃ x 60 √t2 1 ≃ − t2 ≃ − r where we introduced the dimensionless variables: Physically, for diﬀerent values of x and t, three regimes are expected: b Ll t = cosh η = ; x = . (3.40) 0 a a2 x a) Ll a2

2 2 2c 2 where in the first case we should recover the solution for (R R′) = [cosh η0 − (cosh η0 cos θ)(cosh η0 cos θ′) a polymer in the plane, in the second we must obtain 2 − − the result for a polymer on a cylinder and in the third, sinh η cos(φ φ′) cos(θ θ′)]. (3.41) − 0 − − − the mean square end-to-end distance should reach a constant value (corresponding to a large winding number Using the Green function (3.36) in the formula for the about both the Oz axis and the axial circle of radius b). mean-square end-to-end distance (2.9) with the proper volume element on the surface of the torus, evaluating These regimes are manifest in Fig. 4, where the results for the torus are plotted with thin lines, for t = the integrals over the azimuthal angles φ and φ′ by ap- 2;4;6;8;10 . plying the definition (3.37) of the theta function and the { } formula: Initially (the first regime in Eq. (3.47)), the polymer is too small to explore the geometry of the surface, and we 2π ′ im(φ φ ) 2 recover (numerically) the case of a polymer in a plane. dφ dφ′ cos(φ φ′) e − =2π (δ 1m + δ+1m), − − This is also the behavior found for the other surfaces. ZZ0 (3.42) As the length of the chain increases, it starts winding along the circle of radius a and then along the axial and expressing all quantities in terms of the dimension- circle, about the Oz axis. When t increases, the behav- less variables from Eq. (3.40), we eventually obtain the ior is similar—for a certain range of x values—to that mean-square end-to-end distance as of a chain on a cylinder: the mean-square end-to-end distance is linear in Ll/2. We can check analytically 2 2 2π (R R′) 2(t 1) that the approximate propagator from Eq. (3.36) pro- 2 2 −2 = − dθ dθ′ duces the correct behavior. Because a Ll b , the a N(x, t) 0 ≪ ≪ ZZ ′ term n = 0 dominates in the theta function in Eq. (3.43) x(t−cos θ)(t−cos θ ) 2 2 2− 2 2 t + (1 t )e− 4(t 1) cos(θ θ′) (Θ3(z, τ) 1) and the exponential coefficient of 1 t − − − (3.43) ≈ − × (t cos θ)2(t cos θ )2 can be expanded, so one obtains explicitly − − ′ θ θ′ ix(t cos θ)(t cos θ′) Θ − , − − , R R 2 2 Ll 2 2 × 3 2 4π(t2 1) ( ′) 2a + (a Ll b ), (3.48) − − ≃ 2 ≪ ≪

7 which is just the limit for the polymer on the cylinder of end-to-end distance of a Gaussian chain on the geomet- radius a from Eq. (3.15). rical parameters of the confining surface. The formulas Eventually, at any given t and large enough x, for spheres and cylinders are simple (Eqs. (3.9,3.15)), de- 2 scribing the crossover between free polymer regime and (R R′) departs from the behavior of a polymer on a cylinder− and approaches a constant value, which is a confined polymer regime. For the case of a cone, the final expression is more complicated, but still exact. For function only of the geometrical parameters a and b. This limiting value can be explicitly calculated. Now we have a torus, there are three asymptotic regimes as outlined 2 2 by Eq. (3.47), and the crossover behavior is obtained nu- a < b Ll, thus once again Θ3(z, τ) 1 but the exponential≪ coefficient of 1 t2 in the formula≈ for distance merically from Eq. (3.43). (3.43) is zero, which gives− the limit: A number of open problems come into attention: the geodesic 2 2 2 2 2 calculation of the mean-square end-to-end dis- (R R′) 2 (a + b ) (a

ACKNOWLEDGMENTS IV. CONCLUSIONS

Due to its relevance to studies of the behavior of poly- Acknowledgment is made to the CUMIRP and Ma- mers confined at interfaces, we have considered here the terials Research Science and Engineering Center at the problem of a linear, Gaussian polymer chain embedded in University of Massachusetts, and the NSF Grant DMR the following surfaces: the SD 1 sphere in D dimensions, − 9625485. The authors are grateful to E. Cattani, J. Dou- cylinder cone curved torus R3 the , the and the in . glas, J. Machta and F. Pedit for useful discussions. We obtained closed formulas for the probability distribution function G(R R′; L) of the end-to-end vector of the chain for the sphere,| the cylinder and the cone and an approximate propagator for the torus. We calculated, analytically in the case of the sphere, the cylinder and the cone, and numerically for the torus, the mean- square end-to-end distance (in the embedding space) of ∗ Also at The National Institute of Materials Physics, Bu- the chain. As such, at least in the limiting regimes pre- cure¸sti, Romania. viously described, the results are also valid for a free [1] B. Alberts, D. Bray, J, Lewis, M. Raff, K. Roberts and quantum particle on a torus if one uses the mapping J. D. Watson, Molecular Biology of the Cell (Garland 2 m Publishing, Inc. New York, 1994). L t t′; l i~ (where t and t′ represent the fi- {nal→ and initial− time→ coordinates} of the particle). [2] L. S. Schulman, Techniques and Applications of Path In- tegration (John Wiley & Sons, London, New York, 1981). Our calculations demonstrate the role played by the [3] G. Blatter, M. V. Feigel’man, V. B. Geshkenbein, A. I. geometry (curvature and shape) of the interface in con- Larkin and V. Vinokur, Rev. Mod. Phys. 66, 1125 (1994). trolling the size of the chain. [4] M. Doi and S. F. Edwards, The Theory of Polymer Dy- namics (Oxford Science Publications, New York, 1986). As the size of the confining surface decreases, the poly- [5] K. F. Freed, Adv. Chem. Phys. 22, 1 (1972). mer size is determined by the parameters of the confin- [6] H. Kleinert, Path Integrals in Quantum Mechanics, ing surface instead of chain length. The crossover be- Statistics and Polymer Physics, 2nd ed. (World Scien- tween the free Gaussian chain limit and the confined limit tific, Singapore, 1995). is controlled by a characteristic geometrical localization [7] D. C. Khandekar, S. V. Lawande and K. V. Bhagwat, 2 area A. A is proportional to a for the cases of spheres Path-Integral Methods and their Applications (World Sci- and cylinders, and there are two characteristic areas pro- entific, Singapore, 1993). 2 2 portional to a and b for the case of curved torus. For [8] F. Wiegel, Introduction to Path-Integral Methods in the cone, the localization area is variable and is deter- Physics and Polymer Science (World Scientific, 1986). mined by the position of the singular point (the vertex) [9] C. Grosche and F. Steiner, Path Integrals on Curved about the chain and by the angle α. Explicit closed for- Manifolds, Z. Phys. C–Particles and Fields 36, 699 mulas are derived for the dependence of the mean-square (1987).

8 [10] A. Erdelyi, W. Magnus, F. Oberhettinger, F. G. Tricomi, York, 1954). Higher Transcendental Functions, Vol. II (McGraw Hill, [15] G. Arfken, Mathematical Methods for Physicists, 2nd ed., New York, 1985). (Academic Press, London, New York, 1970). [11] C. Itzykson, Commun. Math. Phys. 36, 19 (1974). [16] R. C. Ramos Jr., C. C. Bernido and M. V. Carpio- [12] I. S. Gradshteyn and I. M. Ryzhik, Tables of Integrals, Bernido, J. Phys. A: Math. Gen. 27, 8251 (1994). Series and Products, 5th ed., (Academic Press, London, [17] A. Inomata in Path Integrals from meV to MeV, 433, M. New York, 1994). C. Gutzwiller et al. (eds.) (World Scientiﬁc, Singapore, [13] R. Courant and D. Hilbert, Methods of Mathematical 1986). Physics, Vol. I (John Wiley & Sons, London, New York, [18] R. P. Feynman and A. R. Hibbs, Quantum Mechanics 1989). and Path Integrals (McGraw-Hill, new York, 1965). [14] A. Erdelyi, W. Magnus, F. Oberhettinger, F. G. Tricomi, [19] H. P. McKean and I. M. Singer, J. Diﬀ. Geom. 1, 43 Tables of Integral Transforms, Vol. I (McGraw Hill, New (1967).

FIG. 1. The geometrical data that characterize the sphere (A), the cylinder (B), the cone (C) and the torus (D): a is the radius of the sphere, of the cylinder and of the circular cross-section of the torus; b is the radius of the axial circle of the torus; α is the angle of the cone at the vertex.

2 FIG. 2. The mean-square end-to-end distance (R − R′) normalized to Ll for a Gaussian polymer chain embedded on a ′ cone, as a function of t = ρ and α. √Ll

2 FIG. 3. The mean-square end-to-end distance (R − R′) normalized to Ll for a Gaussian polymer chain embedded on a ′ cone, as a function of t = ρ and α: 3–D representation. √Ll

2 2 FIG. 4. The mean-square end-to-end distance (R − R′) normalized to a for a Gaussian polymer chain embedded on Ll t b diﬀerent surfaces, as a function of 2 and = . The symbols joined by thin lines represent the data for the polymer on the 3 a a 3 surface of the curved torus in R : thick solid line, the cylinder in R ; ◦—◦ torus, t = 10; ✷—✷ torus, t = 8; ✸—✸ torus, t = 6; D 1 D △—△ torus, t = 4; ▽—▽ torus, t = 2; thick dotted line, the sphere S − in R .

9 a

A. B.

2α

2a b

C. D.

Fig. 1 (Brownian Motion and Polymer Statistics on Certain Curved Manifolds Mondescu & Muthukumar) Fig. 2 (Brownian Motion and Polymer Statistics on Certain Curved Manifolds − Mondescu/Muthukumar)

1.0

α = {0, π/2} α = π/16 α π 0.8 = /8 α = 3π/16 >/ (L l) 2 α = π/4 α = 5π/16 α = 3π/8 0.6 <(R − R’)

0.4 0.0 0.4 0.8 1.2 1.6 2.0 t = ρ'/( L l)1/2 <(R - R')^2>/(L l)

1 0 0.75 0.5 0.5 1.5 t 1 1 1.5 α 0.5

2 0 Fig.4 Brownian Motion and Polymer Statistics on Certain Curved Manifolds: Mondescu/Muthukumar

L = N l ; t = b/a

200.0 cylinder in R3 torus in R3; t = 10 160.0 t = 8 t = 6 2 t = 4

>/a 120.0 2 t = 2 sphere in RD

80.0 <(R − R’)

40.0

0.0 0.0 1000.0 2000.0 3000.0 4000.0 5000.0 L l/a2