Melt viscosities of lattice polymers using a Kramers potential treatment

Kramers relaxation times $\tau_{K}$ and relaxation times $\tau_{R}$ and $\tau_{G}$ for the end-to-end distances and for center of mass diffusion are calculated for dense systems of athermal lattice chains. $\tau_{K}$ is defined from the response of the radius of gyration to a Kramers potential which approximately describes the effect of a stationary shear flow. It is shown that within an intermediate range of chain lengths N the relaxation times $\tau_{R}$ and $\tau_{K}$ exhibit the same scaling with N, suggesting that N-dependent melt-viscosities for non-entangled chains can be obtained from the Kramers equilibrium concept.Comment: submitted to: Journal of Chemical Physic

Polyethylene under tensile load: strain energy storage and breaking of linear and knotted alkanes probed by first-principles molecular dynamics calculations

The mechanical resistance of a polyethylene strand subject to tension and the way its properties are affected by the presence of a knot is studied using first-principles molecular dynamics calculations. The distribution of strain energy for the knotted chains has a well-defined shape that is very different from the one found in the linear case. The presence of a knot significantly weakens the chain in which it is tied. Chain rupture invariably occurs just outside the entrance to the knot, as is the case for a macroscopic rope.Comment: 8 pages, 11 figures, to appear on J. Chem. Phy

Scaled Particle Theory for Hard Sphere Pairs. I. Mathematical Structure

We develop an extension of the original Reiss-Frisch-Lebowitz scaled particle theory that can serve as a predictive method for the hard sphere pair correlation function g(r). The reversible cavity creation work is analyzed both for a single spherical cavity of arbitrary size, as well as for a pair of identical such spherical cavities with variable center-to-center separation. These quantities lead directly to prediction of g(r). Smooth connection conditions have been identified between the small-cavity situation where the work can be exactly and completely expressed in terms of g(r), and the large-cavity regime where macroscopic properties become relevant. Closure conditions emerge which produce a nonlinear integral equation that must be satisfied by the pair correlation function. This integral equation has a structure which straightforwardly generates a solution that is a power series in density. The results of this series replicate the exact second and third virial coefficients for the hard sphere system via the contact value of the pair correlation function. The predicted fourth virial coefficient is approximately 0.6 percent lower than the known exact value. Detailed numerical analysis of the nonlinear integral equation has been deferred to the sequel (following paper

Lagrangian statistics in forced two-dimensional turbulence

We report on simulations of two-dimensional turbulence in the inverse energy cascade regime. Focusing on the statistics of Lagrangian tracer particles, scaling behavior of the probability density functions of velocity fluctuations is investigated. The results are compared to the three-dimensional case. In particular an analysis in terms of compensated cumulants reveals the transition from a strong non-Gaussian behavior with large tails to Gaussianity. The reported computation of correlation functions for the acceleration components sheds light on the underlying dynamics of the tracer particles.Comment: 8 figures, 1 tabl

Emergent states in dense systems of active rods: from swarming to turbulence

Dense suspensions of self-propelled rod-like particles exhibit a fascinating variety of non-equilibrium phenomena. By means of computer simulations of a minimal model for rigid self-propelled colloidal rods with variable shape we explore the generic diagram of emerging states over a large range of rod densities and aspect ratios. The dynamics is studied using a simple numerical scheme for the overdamped noiseless frictional dynamics of a many-body system in which steric forces are dominant over hydrodynamic ones. The different emergent states are identified by various characteristic correlation functions and suitable order parameter fields. At low density and aspect ratio, a disordered phase with no coherent motion precedes a highly-cooperative swarming state at large aspect ratio. Conversely, at high densities weakly anisometric particles show a distinct jamming transition whereas slender particles form dynamic laning patterns. In between there is a large window corresponding to strongly vortical, turbulent flow. The different dynamical states should be verifiable in systems of swimming bacteria and artificial rod-like micro-swimmers.Comment: 14 pages, 8 figure

Stiff polymer in monomer ensemble

We make use of the previously developed formalism for a monomer ensemble and include angular dependence of the segments of the polymer chains thus described. In particular we show how to deal with stiffness when the polymer chain is confined to certain regions. We investigate the stiffness from the perspectives of a differential equation, integral equations, or recursive relations for both continuum and lattice models. Exact analytical solutions are presented for two cases, whereas numerical results are shown for a third case.Comment: 10 pages, including 6 figure

Universal decay of scalar turbulence

The asymptotic decay of passive scalar fields is solved analytically for the Kraichnan model, where the velocity has a short correlation time. At long times, two universality classes are found, both characterized by a distribution of the scalar -- generally non-Gaussian -- with global self-similar evolution in time. Analogous behavior is found numerically with a more realistic flow resulting from an inverse energy cascade.Comment: 4 pages, 3 Postscript figures, submitted to PR