arxiv: v1 [math.pr] 10 Nov 2017

Transcription

1 arxiv:7.43v math.pr] Nov 7 Convergence of the Freely Rotating Chain to the Kratky-Porod Model of Semi-flexible Polymers Philip Kilanowski, Peter March, Marko Šamara June 8, 8 Abstract The freely rotating chain is one of the classic discrete models of a polymer in dilute solution. It consists of a broken line of N straight segments of fixed length such that the angle between adjacent segments is constant and the N torsional angles are independent, identically distributed, uniform random variables. We provide a rigorous proof of a folklore result in the chemicahysics literature stating that under an appropriate scaling, as N, the freely rotating chain converges to a random curve defined by the property that its derivative with respect to arclength is a brownian motion on the unit sphere. This is the Kratky-Porod model of semiflexible polymers. We also investigate limits of the model when a stiffness parameter, called the persistence length, tends to zero or infinity. The main idea is to introduce orthogonal frames adapted to the polymer and to express conformational changes in the polymer in terms of stochastic equations for the rotation of these frames. Introduction The simplest caricature of a polymer in dilute solution is a sequence of random vectors R, R,..., R N in R 3 joined together by line segments Q n = R n R n, n N. The freely jointed chain, or bead/rod model, introduced by Kramers in 946 ], assumes the segments, or bonds, Q n are independent, identically distributed, random variables uniformly distributed on a sphere of radius a. A closely related variant, introduced by Rouse in 953 ], called the bead/spring model, assumes the segments Q n are i.i.d gaussian random vectors. The physical interpretation in Kramers model is that statistically significant sections of the polymer behave as if they were inextensible rods that can rotate freely with respect to one another. In Rouse s model, the physical interpretation Dominican House of Studies, Washington DC Rutgers University, Piscataway, NJ Arizona State University, Tempe AZ This research was ported by an Intergovernmental Personnel Act Mobility Assignment NSF DMS 64959

2 is that statistically significant sections of the polymer behave as if they were linear springs connecting N + beads in which the thermal forces of the solvent acting on the beads are equilibrated by the elastic restoring forces of the springs. These models provide a reasonable description of ideal or flexible polymers, in which correlations between monomers decay rapidly with distance along the polymer. That they provide equivalent descriptions of the large scale properties of flexible polymers can be seen from the fact that the bead/rod and bead/spring models are N step random walks that converge to three dimensional brownian motion, as N, under a suitable scaling. However, these are poor models for stiff or semi-flexible polymers, in which the monomer-monomer correlation length is comparable to the length of the polymer itself. To address this issue, Kirkwood and Riseman 3], 4] introduced in 948 a variant of the Kramers chain, called the freely rotating chain, in which geometric constraints are imposed on the segments so that they are correlated with one another and no longer independent. In this model, the bonds have a common length a; the angle between adjacent bonds, or bond angle, has a fixed value θ; and the bond Q n is chosen independently and uniformly at random from the cone of axis Q n, aperture θ and side length a, (cf. Figure ). While the bonds Q n are no longer independent, they do form a Markov chain, namely, a geodesic random walk on the sphere of radius a of step size aθ. We recommend standard references to the classicaolymer physics literature such as 5], 6], 7] and 8] for more in depth discussion of the basic physicahenomena and the standard mathematical models of polymers. Figure : Construction of the discrete polymer from its bonds Q n At the same time, a continuum version of the freely rotating chain was proposed independently by Kratky and Porod 9]. Their model is defined by a

3 3 Hamiltonian for a curve R(s) in terms of its derivative with respect to arclength u(s) = R(s) s, L H(R) = λ u(s) s ds. Here λ is a physical constant, and L is the arclength of the polymer. By definition, the tangent vector u(s) has unit length, but it is often pointed out explicitly that u(s), because in the physical formalism, this fact greatly complicates the implied functional integration over paths that is necessary to obtain the partition function of the model. For this reason, the Kratky-Porod model is considered to be somewhat difficult to work with. A formal calculation relates the Kratky-Porod model to the scaling limit of the freely rotating model. Since Q k Q k = a ( cos θ), it follows: first, that the difference quotients (Q k Q k )/a approximate a derivative with respect to arclength, provided a = L/N; and second, that the sum of squares of these difference quotients is bounded if θ = O(/ N); say, θ = κ/ N. This formal calculation is the basis of a folklore result in the chemical physics literature that the Kratky-Porod model is the limit of the freely jointed model under the indicated scaling. Now, it is difficult to make rigorous sense of the Hamiltonian formalism of Kratky and Porod, and we do not try to do so. Rather, we show that under this scaling, the freely jointed chain converges to a random curve whose derivative with respect to arclength can be identified explicitly as a spherical brownian motion. This is easy to guess, referring back to the probabilistic description of Q n as a geodesic random walk on the sphere of radius a. Evidently, the normalized segments Q n /a lie on the unit sphere and the steps Q n Q n have uniformly distributed directions and are taken along great circles of arclength O(/ N), which is the correct scaling for convergence to spherical brownian motion. The first result of this article makes these arguments precise and rigorously defines the Kratky-Porod model. Note that throughout this work, we use pinned boundary conditions at s =. Thus, the initial bead is always R =, and the initial bond is always in the vertical direction in the laboratory frame, namely, Q = ae 3 where e 3 = (,, ). Theorem. Let R N n, n N, be the freely rotating chain with bond length a = L/N and bond angle θ = κ/ N. Let R N (s), s L be the piecewise linear curve obtained by linear interpolation of the beads R N (nl/n) = R N n. Let Q s (λ), s L be brownian motion on the unit sphere S, generated by λ S, starting from the north pole Q (λ) = e 3 = (,, ) and let = L/κ. Then the processes R N converge weakly, as N, to the Kratky-Porod model defined by

4 4 R s = Q u ( )du, s L. Note that the unit tangent process Q s ( ) is not standard brownian motion on S since it is generated by S rather than the standard S. The reason for the non-standard normalization is to focus attention on the parameter = L/κ, called the persistence length. Since EQ s ( ) Q t ( )] = e t s / () as a standard calculation shows (see Section ), is a correlation length that describes the exponential decay of tangent-tangent correlations along the polymer. Values of that are large compared to the polymer length L are indicative of stiff, rod-like polymers, whereas values that are small compared to L are indicative of flexible, random coiolymers. Thus, the dimensionless ratio L = κ characterizes the intrinsic stiffness of the polymer and gives a physical interpretation of the bond angle parameter κ. The persistence length is the key parameter describing a qualitative transition in the Kratky-Porod model which is recorded in the second main result of this article. Theorem. The Kratky-Porod model with fixed length L exhibits the hard rod/random coil transition, in the sense that: A. (Hard Rod) As, (R s ; s L) (se 3 ; s L) in probability, and ( (R s se 3 ); s L) converges weakly to a Gaussian process (W s ; s L) with components Ws 3, and for i =,, W i s = β i udu, s L where the β i s are independent standard one dimensional brownian motions starting at zero. B. (Random Coil) As, ( 3/ R s ; s L) converges weakly to (W s ; s L), where W is a standard 3-dimensional Brownian motion. The significance of Theorem is that it shows how the Kratky-Porod model interpolates between two extreme polymer types. It confirms, as expected, that the freely rotating chain model for semi-flexible polymers agrees with the freely jointed chain and the bead/spring chain in the large N and small regime. It also confirms the expected hard rod limit as, but the precise form of the gaussian fluctuations on the scale O( ) seems to be new. This regime is known as the weakly bending limit in the polymer physics literature, e.g., 6].

5 5 The main idea in the proofs of these results is a representation Q n = Z n Q of the n th bond of the freely rotating chain as a random rotation of the initial bond. The Markovian structure of the bond process Q n leads to a simple representation of the rotation process Z n in terms of i.i.d. rotations. It develops the fact that Z n is the solution of a simple stochastic difference equation in the Lie group SO(3) that corresponds, in the scaling limit, to the Stratonovich SDE, Z s = Z s B s () Here B(s) is a brownian motion in the Lie algebra of 3-by-3 anti-symmetric matrices given by B(s) = β s βs (3) βs βs where the β i s are independent standard one dimensional brownian motions starting at. The weak convergence of (Bk N, ZN k, RN k ) to (B(s), Z(s), R(s)) is guaranteed by a theorem of Kurtz and Protter ] and the identification of Q s ( ) = Z(s)e 3 as a spherical brownian motion is standard exercise in stochastic calculus. At this point it is useful to recall that each rotation in SO(3) can be identified with an orthonormal frame in R 3, namely, the frame obtained by rotation of a fixed reference frame. Under this identification the SDE for Z above can be understood as an equation for a moving frame. This SDE is not a stochastic version of the Frenet-Serret frame equations, as one might have thought at the outset, but rather it is a stochastic version of the parallel transport or Bishop frame equations ]. (Experts in stochastic differential geometry will recognize the stochastic Bishop frame Z(s) as a brownian motion on the orthonormal frame bundle of S and Z(s)e 3 as the bundle projection onto the base). Even in the deterministic case, the Bishop frame equations are more general than the Frenet-Serret frame equations in the sense that the former require less regularity of the underlying curve than the latter. This is crucial in the present case, since the Bishop frame parameters that define the Lie algebra process B(s) are a pair of independent one-dimensional brownian motions whereas neither the curvature nor the torsion of R(s) exists in any obvious sense. It is possible to prove the convergence of the freely rotating chain by adapting the results of Pinsky on isotropic transport processes ], 3]. However, we prefer to develop an approach based on stochastic Bishop frames because it is self-contained and allows us easily to derive the asymptotic behavior of the Kratky-Porod model as the persistence length tends to zero or infinity.

6 6 Convergence Theorems. Construction of the Model Consider a polymer chain of N segments of length a, with each pair of consecutive segments meeting at the same planar angle. Denote by θ the angle between two consecutive bonds (see figure ). Denote the vectors comprising the segments of the polymer by Q, Q,..., Q N, and let the position of each bead on the polymer, or bond between segments, be denoted by R, R,..., R N. Without loss of generality, we can specify the position and orientation of the polymer in R 3 by choosing our coordinates appropriately. Fix one end of the polymer (the zeroth bead or fixed end) at the origin, so that R =. It then follows that n R n = Q i. (4) i= Next, orient the polymer so that the first segment always points straight up, in the positive z-direction. That is, R = Q = ae 3. (5) Pinning down of one end of the polymer specifies boundary conditions for the model. Since the bond angle θ is fixed, the bond Q n can be expressed as a rotation of Q n through angle θ about an axis chosen uniformly at random from the unit vectors in the plane perpendicular to Q n. We denote that rotation by G n and note that in the freely rotating chain, the G n s are assumed to be mutually independent. Let Z n = G n G n... G and note that Z n = G n Z n. Since Q n = G n Q n, it follows that Q n = Z n Q. Thinking of Z n as an orthogonal change of basis, we see that the conjugate rotation H n = Z n G nz n (6) is a rotation of Q through an angle θ about an axis chosen uniformly at random from the unit vectors in the plane perpendicular to Q = ae 3, and that the H n s are mutually independent (see figure ). Thus, Z n = G n Z n (7) = Z n Zn nz n (8) = Z n H n (9) = H H... H n. () In the sequel, we work with an explicit model of the freely rotating chain by defining H n to be the rotation through an angle θ about the axis (cos φ n, sin φ n, )

7 7 Figure : Commutative diagram for H n and G n where φ,..., φ N are i.i.d. random variables uniformly distributed on, π). Because of the product structure of Z n, the freely rotating chain satisfies a simple system of discrete stochastic equations R n = n Z k (ae 3 ) () k= Z n Z n = Z n (H n I) () R =, Z = I. (3) To analyze the system in the limit N, we observe that if N is large, then θ = κ/ N is small, hence H n is close to the identity. Such rotations can be expressed as matrix exponentials of the form ( ) κ H n = exp b n (4) N where b n = sin φ n cos φ n sin φ n cos φ n This observation plays a key role in the following result. (5) Theorem 3. Let a = L N, θ = κ and = L. For s L, define N κ Z N s = Z N Ns/L], and BN s = Ns/L] k= (H n I).

8 8 Then, the process ((Zs N, Bs N ); s L) converges weakly to the process ((Z s, B s ); s L), where βs B s = βs (6) βs βs and β and β are independent standard one dimensional Brownian motions starting at zero. Furthermore, Z s satisfies the Stratonovich equation with initial condition Z = I. Z s = Z s B s (7) Proof : This theorem is proved in several steps. First, we show that the nonzero entries of Bs N converge to Brownian motions whose variance at time s is s/. Second, we use a theorem of Kurtz and ( Protter ] to show that ) the pair (Bs N, Zs N ; s L) converges weakly to B s, Z s ; s L, where Z s satisfies the Stratonovich equation Z s = Z s B s. For the first step, note that Ns L ] Bs N = n= ( ) κ b n + κ N N b n + r n (8) N where r n eκ/ κ 3 6N. 3/ Note that the non-zero entries of b n are mutually independent random variables with mean and variance κ. Therefore, by Donsker s theorem, the non-zero N entries in the first term of (8) converge to independent Brownian motions whose variance at time s is sκ L = s. Next, by the Law of Large Numbers, the second term in (8) converges to s D, where D = diag( /, /, ). Finally, the third term in (8) is O(N 3/ ) and hence converges ( to zero. process (B N Ns L ] ) is tight and converges weakly to B s ; s L Therefore, ) the. Note that B N s and Z N s are piecewise constant and have jump discontinuities only at

9 9 integer multiples of L nl N. Now, if s = N, then Z N s = Z N s Z N s (9) = Z N n Z N n () = Z N n (H N n I) () = Z n s B N s. () Hence, Z N s satisfies the Stratonovich equation Zs N = Zs B N s N. (3) But by (8), we have B N, B N ] s = Cs for a constant C independent of N. Thus, the sequence (Bs N, Zs N ; s L) is tight in the Skorohod topology. Consequently, by the theorem of Kurtz and Protter ], this sequence converges weakly to (, Z s ; s L) where Z s = Z s B s (4) which is equivalent to equation (7) since B s is continuous. Theorem 4. Let a = L N, θ = κ and = L, and define the process N κ Q s = Z s e 3, s L. Then, Q s ( ) is a spherical Brownian motion generated by S. Consequently, the freely jointed chain (RNs/L] N ; s L) converges weakly to the process R s = Z u e 3 du, s L. Proof : The strategy is to use Lévy s characterization theorem, by which it suffices to show that for any C function f : S R, the process f(q(s)) f(q()) S f(q(u))du is a martingale. Let f : S R be a smooth function and let F : R 3 R be a smooth extension of f such that for some ε > and all ξ S, F (rξ) = f(ξ), ε < r < + ε. (5) Then, the radial derivative of F vanishes, and so the expression of the Laplacian in R 3 in polar coordinates simplifies to R 3F (ξ) = S F (ξ) = S f(ξ) (6)

10 for all ξ S. Let X s = Z s e, Y s = Z s e and Q s = Z s e 3. Then the matrix equation Z s = Z s B s is equivalent to the system X s = Q s β s (7) Y s = Q s β s (8) Q s = ( Xs β s + Y s β s). (9) Applying Itô s formula to Q s thought of as a process in R 3 yields f(q s ) f(e 3 ) = F (Q s ) F (e 3 ) (3) = F (Q s ) Q s (3) = F (Q u ) dq u + Hess F (Q u ) : d Q, Q u (3) The first term above is a martingale. To evaluate the integrand of the second term, note that hence d Q, Q u = (X u X u )du (33) Hess F (Q u ) : d Q, Q u = X u F (Q u ) + Yu F (Q u ) (34) where the directional derivatives are computed ω by ω. Since the unit vector Q u points almost surely in the radial direction, we have and hence Q u F (Q u ) = r F (Q u) = a.s. (35) ( Hess F (Q u ) : d Q, Q u = ) + Q F (Q u ) (36) u Xu + Yu = R 3 F (Q u ) a.s. (37) = S f(q u ). (38)

11 Here we have used the fact that ( R 3F (x) = u u for any orthogonal frame (u, v, w) based at x R 3. ) + vu + wu F (x) (39) 3 Asymptotic Behavior in the Persistence Length One advantage of the representation of the tangent process Q s = dr s ds = Z se 3 in terms of the SO(3)-valued process Z with Z s = Z s B s is that we can analyze the behavior of the Kratky-Porod model when the persistence length tends to or infinity. We first need a lemma about the correlation between vectors Q s and Q t at different times (i.e. arc lengths) s and t. Lemma 5. Let Q s = Z s e 3 be the process as above. Then EQ s Q t ] = e t s. Proof : Suppose s < t. Since Q s Q t is rotationally invariant, it follows that for any ξ S, and rotation O SO(3) with Oξ = e 3, E ξ Q Q s ] = E ξ OQ OQ s ] (4) = E e3 e 3 Q s ] (4) By the Markov property, E e3 Q s Q t ] = E e3 E Qs Q Q t s ]] = E e3 Q Q t s ]. Therefore, to prove the lemma, it suffices to show EQ s e 3 ] = e s, for all s, L]. To this end, define f : R 3 R by f(p ) = P e 3. Alternatively, f has the form f(p ) = r cos φ in polar coordinates of the point P = (r, θ, φ). Then R 3f(P ) = r f(p ) + r r f(p ) + r Sf(P ) (4) = r cos φ + r Sf(P ). (43) On the other hand, R 3f, as f is a linear function of cartesian coordinates of its argument. Thus, S f(p ) = r cos φ = f(p ). (44) Now, applying Itô s formula for Q s, we find that the process

12 M s = f(q s ( )) f(q ( )) S f(q u ( ))du (45) = f(q s ( )) + f(q u ( ))du (46) is a martingale with mean. Taking the expectation in (46), we conclude that the non-random function x s = Ef(Q s ( ))] is the solution of the initial value problem ẋ s = x s, x = ; (47) hence, EQ s e 3 ] = x s = e s. Lemma 6. E R t ] = t Proof : E R t ] t = E = = ) ( e t. t v t v Q u du = t t ] Q v dv (48) E Q u Q v ] du dv (49) e (v u) du dv (5) ) ( e t. (5) Theorem 7. The Kratky-Porod model with fixed length L exhibits the hard rod/random coil transition, in the sense that, A. (Hard Rod) As, (R s ; s L) (se 3 ; s L) in probability, and ( (R s se 3 ); s L) converges weakly to a Gaussian process (W s ; s L) with components Ws 3, and for i =,, W i s = β i udu, s L where the β i s are independent, standard, one dimensional brownian motions starting at zero. B. (Random Coil) As, ( 3/ R s ; s L) converges weakly to (W s ; s L), where W is a standard 3-dimensional Brownian motion.

13 3 Proof of A. Since it is enough to show that R s se 3 = (Z u I)e 3 du, (5) (Z u I) du in probability as. To this end, for a matrix-valued process M = m ij ] 3 i,j= we define a norm M pathwise by the formula M = 3 i,j= s L m ij s (53) and we wilrove that the random variable (Z u I) du converges to in probability. Now, it is not hard to see that every integrable matrix-valued process M satisfies M u du L M. (54) Thus, it suffices to show Z I in probability as. For ε >, we have { P { Z I ε} = P Z u B u ε } ] ε E Z u B u (55) and so it suffices to show that the expected value on the right hand side of (55) is bounded. Note that the matrix Z u B u has entries of the form z ij u β k u, for all i, j 3 and k =,, or perhaps a sum of two such integrals. Assuming C = max E s ] i,j,k zu ij βu k <, (56) then E s L Z u B u 36C, (57) as one can see by squaring out the terms of the sum in the definition of the norm M and estimating the cross term by the Cauchy-Schwarz inequality. Thus, to prove that the expected value on the right hand side of (55) is bounded, it suffices to show that

14 4 E s L z u β u ] < (58) (where we have dropped the erscripts for ease of notation). The conversion of Stratonovich to Îto integrals gives us: Observe that s L z u dβ u + z, β s s L Also, by Jensen s inequality, E s L z u β u = z u dβ u + s L z u dβ u + z, β s. (59) z u dβ u s L z, β s + 4 z, β s. (6) s L z u dβ u ] E s L z u dβ u ] /. (6) Therefore, to prove (58), it suffices to show ] E z u dβ u < (6) s L and that for some constant K, almost surely, s L z, β s K. (63) By Doob s inequality, and using the fact that z u almost surely, E s L ] ] L z u dβ u 4E zu du (64) 4L (65) which shows (6). To prove (63), we note that z, β s can be written as the sum of at most two terms of the form l / p z u βu, i β, with z being an entry

15 5 of Z. Using d β i, β j u = δ ij du, it follows that z, β s is equal to only one term of the form l / p z u du. Thus, z, β s s L s L max zu ij du (66) i,j L (67) which completes the proof of the first part of A. Recalling the definition of the Gaussian proces W s, s L, let us observe that (R s se 3 ) W s = (Z u I)e 3 du = B u e 3 du (68) ( (Z u I) B u ) e 3 du. (69) Let us show that the norm of the expression in (69) converges to in probability, as. However, by (54), it is enough to show that (Z I) B (7) in probability. Now, iterating the equation Z s = l / p Z s B s twice, it follows that for all s, L], ( s Z s I = I + ) u Z v B v B u (7) which yields = B s + u Z v B v B u (7) (Z s I) B s = u Therefore, for all ε >, P { ( ) } Z I B ε 3 ε E Z v B v B u. (73) x ij s i,j= s L (74)

16 6 where x ij s ] = X s = u we observe that each entry x ij s of the form Z v B v B u. As in the proof of the first part of A, of the matrix X s is the sum of at most two terms u z v β k v β l u (75) where z v is an entry of Z v. To complete the proof, it suffices to show that u ] E z v βv k βu l <. (76) s L This is obvious if x ij is equal to just one term of the form (75). If, however, x ij is a sum of two such terms, the sufficency of (76) becomes clear after expanding E x ij s ] and using the Cauchy-Schwarz inequality. Now, to see (76), we convert from Stratonovich to Itô integrals and obtain u z v β k v β l u = = ( u z v dβv k + ) z, βk s βu l (77) z v dβv k dβu l + z v dβv k, β l u + z, β k u dβ l u + 4 In fact, z, β, β k s =, since z, β is of finite variation. Also, s L z v dβv k, β l = δ kl s s L s z, β, β k s. (78) almost surely. Furthermore, by Doob s inequality and (67), z v dv (79) L (8) E s L ] ] L z, β k u dβu l 4 E z, β k u du (8) 4L3. (8) Therefore, squaring (78), and using (67), (8) and (8), we see that to prove (76) for Stratonovich integrals, it suffices to prove the boundedness of the corresponding quantity for Îto integrals, that is, u ] E z v dβv k dβu l <. (83) s L

17 7 But, by Doob s inequality, u ] L u ] E z v dβv k dβ l 4 E z v dβ k v du (84) s L u ] 4L E z v dβv k (85) u L ] L 6L E z u du (86) This completes the proof of part A. 6L. (87) Proof of B. We would like to prove a functional central limit theorem for (R s ; s L) by viewing R s = (R s, R s, R 3 s) as a vector of additive functionals of the underlying spherical Brownian motion: R i s( ) = = Q i u( ) du (88) f i (Q u ( )) du (89) where f i : S R is the i th coordinate function of Q, i.e. f i (Q) = Q i, for i =,, 3. Let g i : S R be a solution of the Poisson equation S g i = f i, and let σ be the normalized Lebesgue measure on the unit sphere. By symmetry, f i (Q)σ(dQ) =, and since f i is smooth, there exists a smooth solution S g i, unique up to an additive constant. By a change of variable, the process (Q s ( ); s ) is equal in distribution to (Q s (); s ), as well as (R s ( ); s ) = D ( R s (); s ). (9) Using Itô s formula and S g i = f i, we find that is a martingale. So, R i s( ) = g i (Q s ( )) g i (Q ( )) f i (Q u ( )) du (9) f i (Q u ( )) du = g i (Q s ( )) g i (Q ( )) + M i s( ) (9) for some martingale M i s( ), i 3. In particular, for any t,

18 8 R i t() = g i (Q t ()) g i (Q ()) + M i t (). (93) Comparing (93) for t = s/ with (9), and using the fact that Q s ( ) D = Q s () along with (9), it follows that Ms(l i p ) = D M i s (). This yields Rs(l i p ) = ] g i (Q s ()) gi (Q ()) + M i s (). (94) In fact, M i s(), which we denote for short by M i s, i 3, is a system of continuous martingales with predictable quadratic covariation M i, M j t = t where Γ is the carré du champ operator: Γ(g i, g j )(Q s ()) ds (95) Γ(ϕ, ψ) = ] S (ϕψ)(q) ϕ(q) S ψ(q) ψ(q) S ϕ(q). (96) Since g i is bounded, the functional central limit theorem for the stochastic as is equivalent to the process / R i s( ) = R i s () = R i s corresponding result for the martingale M i s. According to the functional central limit theorem for martingales (e.g. 4], 5]), this amounts to showing that lim λ λ M i, M j λt = c ij t (97) for some non-negative definite matrix C = c ij ]. By the ergodic theorem for Q, t c ij = lim t t M i, M j t = lim Γ(g i, g j )(Q s ) ds (98) t t = Γ(g i, g j )(Q)σ(dQ) (99) S almost surely and in L. Rather than calculating Γ(g i, g j ) explicitly, we can identify C by a symmetry argument. We will show that OCO T = C for any rotation O SO(3), which implies that C is a constant multiple of the identity matrix. To see this, let O SO(3) and define a process R as R s = OR s = = OQ u () du () Q u () du ()

19 9 where Q = OQ is a spherical brownian motion generated by S, starting at Oe 3. Because R is a Kratky-Porod model, R differs from a martingale, say, M by a bounded term, due to equation (9). Hence, On the other hand, lim t t E R t RT t ] = lim t t E M T t M t ] () = lim t t E M, M t (3) = lim t t E M, M t = C. (4) C = lim t t E R t RT t ] = lim t t OERRT ]O T (5) So, C = ci for some constant c. Now, using Lemma 6, = lim t t OE M, M to T (6) = OCO T. (7) Now, and so c = /3. This completes the proof of Theorem 7. 4 Conclusion lim t t E R t ] = lim t t E R t ] =. (8) lim t t E R t ] = tr C = 3c (9) The Kratky-Porod model for semi-flexible polymers in dilute solution is a continuous model that proves to be the continuum limit of the discrete freely rotating chain model, under an appropriate scaling. While this result had been known for over sixty years in the field of chemicahysics, our Theorem 4 provides the first known rigorous, probabilistic proof of the correspondence between the two models. This is shown by means of an intermediate result (Theorem 3) that relates the tangent vector along the polymer to a sequence of matrix-valued processes, defined by a stochastic differential equation. This sequence is shown to converge by a theorem of Kurtz and Protter, and the limiting vector-valued process is found to be a spherical Brownian motion that matches the continuous Kratky-Porod model.

20 Moreover, the model is proven to converge to the expected results in two extreme cases: for long persistence lengths (large values of ), the adjacent segments are highly correlated, and so the polymer approximates a straight rod; while for short persistence lengths (small values of ), the adjacent segments are nearly independent, and the polymer approaches the Rouse or freely jointed chain model, in which the tangent vector approaches a three-dimensional Brownian motion. These limiting cases are shown in Theorem 7, which provides the additional result that in the former case, the deviation from the straight rod is itself a Brownian motion in the plane. This result places the Kratky-Porod model on a continuum with the straight rod and the Rouse model as the stiffness parameter, or persistence length, varies. With these results firmly in place, we can expand the Kratky-Porod model into other realms. One such case, in which a constant external force is applied to the polymer, introduces other limiting regimes in which the shape and orientation of the polymer depend both on the persistence length and relative strength and direction of the force. Another case is the Rotational Isomeric State approximation model, in which conformation of the polymer is determined by finite and fixed number of possible states in each step of the conformation. These models, therefore, unite the fields of probability theory and chemical physics, using the methods and tools of the former field to demonstrate a property in the latter. At long last, the relationship between the freely rotating and Kratky-Porod polymer models has been proven. References ] H. A. Kramers, The Behavior of Macromolecules in Inhomogeneous Flow, J. Chem. Phys., 4 (946), pp ] P. E. Rouse, A Theory of the Linear Viscoelastic Properties of Dilute Solutions of Coiling Polymers, J. Chem. Phys., (953), pp ] J. G. Kirkwood and J. Riseman, The Intrinsic Viscosities and Diffusion Constants of Flexible Macromolecules in Solution, J. Chem. Phys., 6 (948), pp ] J. G. Kirkwood and J. Riseman, The Rotatory Diffusion Constants of Flexible Molecules, J. Chem. Phys., (949), pp ] J. P. Flory, Statistical Mechanics of Chain Molecules, Carl Hanser Verlag (989) 6] M. Doi and S. F. Edwards, The Theory of Polymer Dynamics, (986) Oxford University Press Inc., New York

21 7] M. V. Volkenstein, Configurational Statistics of Polymeric Chains, (963) John Wiley & Sons, Ltd. 8] H. Yamakawa, Modern Theory of Polymer Solutions, (97) Harper & Row, New York 9] O. Kratky and G. Porod, Diffuse Small-angle Scattering of X-rays in Colloid Systems, J. Colloid Sci. 4 (949), pp ] Th.G. Kurtz and Ph.E. Protter, Weak convergence of stochastic integrals and differential equations, Probabilistic Models for Nonlinear Partial Differential Equations (Montecatini Terme, 995), Lecture Notes in Math. 67, pp. -4, Springer, Berlin. ] R. L. Bishop, There is More than One Way to Frame a Curve, The American Mathematical Monthly, Vol. 8, No.3 (Mar., 975), pp ] M. A. Pinsky, Isotropic Transport Process on a Riemannian Manifold, Trans. Amer. Math. Soc. 8 (976), pp ] M. A. Pinsky, Homogenization in Stochastic Differential Geometry, Publ. RIMS, Kyoto Univ., 7 (98) pp ] S. N. Ethier and Th. G. Kurtz Markov Processes: Characterization and Convergence (986) J. Wiley & Sons, New York 5] W. Whitt Proofs of the Martingale FCLT, Probability Surveys Probability Surveys Vol. 4 (7) pp. 683