Comment on "Elasticity Model of a Supercoiled DNA Molecule"

We perform simulations to numerically study the writhe distribution of a stiff polymer. We compare with analytic results of Bouchiat and Mezard (PRL 80 1556- (1998); cond-mat/9706050).Comment: 1 page, 1 figure revtex

Non-equilibrium hydrodynamics of a rotating filament

The nonlinear dynamics of an elastic filament that is forced to rotate at its base is studied by hydrodynamic simulation techniques; coupling between stretch, bend, twist elasticity and thermal fluctuations is included. The twirling-overwhirling transition is located and found to be strongly discontinuous. For finite bend and twist persistence length, thermal fluctuations lower the threshold rotational frequency, for infinite persistence length the threshold agrees with previous analytical predictions

Getting DNA twist rigidity from single molecule experiments

We use an elastic rod model with contact to study the extension versus rotation diagrams of single supercoiled DNA molecules. We reproduce quantitatively the supercoiling response of overtwisted DNA and, using experimental data, we get an estimation of the effective supercoiling radius and of the twist rigidity of B-DNA. We find that unlike the bending rigidity, the twist rigidity of DNA seems to vary widely with the nature and concentration of the salt buffer in which it is immerged

Diffusion of a ring polymer in good solution via the Brownian dynamics

Diffusion constants D_{R} and D_{L} of ring and linear polymers of the same molecular weight in a good solvent, respectively, have been evaluated through the Brownian dynamics with hydrodynamic interaction. The ratio $C=D_{R}/D_{L}$, which should be universal in the context of the renormalization group, has been estimated as $C= 1.11 \pm 0.01$ for the large-N limit. It should be consistent with that of synthetic polymers, while it is smaller than that of DNAs such as $C \approx 1.3$. Furthermore, the probability of the ring polymer being a nontrivial knot is found to be very small, while bond crossings may occur at almost all time steps in the present simulation that realizes the good solvent conditions.Comment: 11 pages, 4 figure

Torsional Directed Walks, Entropic Elasticity, and DNA Twist Stiffness

DNA and other biopolymers differ from classical polymers due to their torsional stiffness. This property changes the statistical character of their conformations under tension from a classical random walk to a problem we call the `torsional directed walk'. Motivated by a recent experiment on single lambda-DNA molecules [Strick et al., Science 271 (1996) 1835], we formulate the torsional directed walk problem and solve it analytically in the appropriate force regime. Our technique affords a direct physical determination of the microscopic twist stiffness C and twist-stretch coupling D relevant for DNA functionality. The theory quantitatively fits existing experimental data for relative extension as a function of overtwist over a wide range of applied force; fitting to the experimental data yields the numerical values C=120nm and D=50nm. Future experiments will refine these values. We also predict that the phenomenon of reduction of effective twist stiffness by bend fluctuations should be testable in future single-molecule experiments, and we give its analytic form.Comment: Plain TeX, harvmac, epsf; postscript available at http://dept.physics.upenn.edu/~nelson/index.shtm

Universality in the diffusion of knots

We have evaluated a universal ratio between diffusion constants of the ring polymer with a given knot $K$ and a linear polymer with the same molecular weight in solution through the Brownian dynamics under hydrodynamic interaction. The ratio is found to be constant with respect to the number of monomers, $N$, and hence the estimate at some $N$ should be valid practically over a wide range of $N$ for various polymer models. Interestingly, the ratio is determined by the average crossing number ($N_{AC}$) of an ideal conformation of knotted curve $K$, i.e. that of the ideal knot. The $N_{AC}$ of ideal knots should therefore be fundamental in the dynamics of knots.Comment: 8 pages, 14 figure

Scattering functions of knotted ring polymers

We discuss the scattering function of a Gaussian random polygon with N nodes under a given topological constraint through simulation. We obtain the Kratky plot of a Gaussian polygon of N=200 having a fixed knot for some different knots such as the trivial, trefoil and figure-eight knots. We find that some characteristic properties of the different Kratky plots are consistent with the distinct values of the mean square radius of gyration for Gaussian polygons with the different knots.Comment: 4pages, 3figures, 3table

Mechanical response of plectonemic DNA: an analytical solution

We consider an elastic rod model for twisted DNA in the plectonemic regime. The molecule is treated as an impenetrable tube with an effective, adjustable radius. The model is solved analytically and we derive formulas for the contact pressure, twisting moment and geometrical parameters of the supercoiled region. We apply our model to magnetic tweezer experiments of a DNA molecule subjected to a tensile force and a torque, and extract mechanical and geometrical quantities from the linear part of the experimental response curve. These reconstructed values are derived in a self-contained manner, and are found to be consistent with those available in the literature.Comment: 14 pages, 4 figure

On the Limits of Analogy Between Self-Avoidance and Topology-Driven Swelling of Polymer Loops

The work addresses the analogy between trivial knotting and excluded volume in looped polymer chains of moderate length, $N<N_0$, where the effects of knotting are small. A simple expression for the swelling seen in trivially knotted loops is described and shown to agree with simulation data. Contrast between this expression and the well known expression for excluded volume polymers leads to a graphical mapping of excluded volume to trivial knots, which may be useful for understanding where the analogy between the two physical forms is valid. The work also includes description of a new method for the computational generation of polymer loops via conditional probability. Although computationally intensive, this method generates loops without statistical bias, and thus is preferable to other loop generation routines in the region $N<N_0$.Comment: 10 pages, 5 figures, supplementary tex file and datafil

Abundance of unknots in various models of polymer loops

A veritable zoo of different knots is seen in the ensemble of looped polymer chains, whether created computationally or observed in vitro. At short loop lengths, the spectrum of knots is dominated by the trivial knot (unknot). The fractional abundance of this topological state in the ensemble of all conformations of the loop of $N$ segments follows a decaying exponential form, $\sim \exp (-N/N_0)$, where $N_0$ marks the crossover from a mostly unknotted (ie topologically simple) to a mostly knotted (ie topologically complex) ensemble. In the present work we use computational simulation to look closer into the variation of $N_0$ for a variety of polymer models. Among models examined, $N_0$ is smallest (about 240) for the model with all segments of the same length, it is somewhat larger (305) for Gaussian distributed segments, and can be very large (up to many thousands) when the segment length distribution has a fat power law tail.Comment: 13 pages, 6 color figure