Isostaticity and the solidification of semiflexible polymer melts

How Segmental Dynamics and Mesh Confinement Determine the Selective Diffusivity of Molecules in Cross-Linked Dense Polymer Networks

The diffusion of molecules ("penetrants") of variable size, shape, and chemistry through dense cross-linked polymer networks is a fundamental scientific problem broadly relevant in materials, polymer, physical, and biological chemistry. Relevant applications include separation membranes, barrier materials, drug delivery, and nanofiltration. A major open question is the relationship between transport, thermodynamic state, and penetrant and polymer chemical structure. Here we combine experiment, simulation, and theory to unravel these competing effects on penetrant transport in rubbery and supercooled polymer permanent networks over a wide range of cross-link densities, size ratios, and temperatures. The crucial importance of the coupling of local penetrant hopping to polymer structural relaxation and the secondary importance of mesh confinement effects are established. Network cross-links strongly slow down nm-scale polymer relaxation, which greatly retards the activated penetrant diffusion. The demonstrated good agreement between experiment, simulation, and theory provides strong support for the size ratio (penetrant diameter to the polymer Kuhn length) as a key variable and the usefulness of coarse-grained simulation and theoretical models that average over Angstrom scale structure. The developed theory provides an understanding of the physical processes underlying the behaviors observed in experiment and simulation and suggests new strategies for enhancing selective polymer membrane design.

Local and Global Order in Dense Packings of Semi-Flexible Polymers of Hard Spheres

The local and global order in dense packings of linear, semi-flexible polymers of tangent hard spheres are studied by employing extensive Monte Carlo simulations at increasing volume fractions. The chain stiffness is controlled by a tunable harmonic potential for the bending angle, whose intensity dictates the rigidity of the polymer backbone as a function of the bending constant and equilibrium angle. The studied angles range between acute and obtuse ones, reaching the limit of rod-like polymers. We analyze how the packing density and chain stiffness affect the chains' ability to self-organize at the local and global levels. The former corresponds to crystallinity, as quantified by the Characteristic Crystallographic Element (CCE) norm descriptor, while the latter is computed through the scalar orientational order parameter. In all cases, we identify the critical volume fraction for the phase transition and gauge the established crystal morphologies, developing a complete phase diagram as a function of packing density and equilibrium bending angle. A plethora of structures are obtained, ranging between random hexagonal closed packed morphologies of mixed character and almost perfect face centered cubic (FCC) and hexagonal close-packed (HCP) crystals at the level of monomers, and nematic mesophases, with prolate and oblate mesogens at the level of chains. For rod-like chains, a delay is observed between the establishment of the long-range nematic order and crystallization as a function of the packing density, while for right-angle chains, both transitions are synchronized. A comparison is also provided against the analogous packings of monomeric and fully flexible chains of hard spheres.

Crystallization of Flexible Chains of Tangent Hard Spheres under Full Confinement

We present results from extensive Monte Carlo simulations on the crystallization of athermal polymers under full confinement. Polymers are represented as freely jointed chains of tangent hard spheres of uniform size. Confinement is applied through the presence of flat, parallel, and impenetrable walls in all dimensions. We analyze crystallization as the summation of two contributions: one that occurs in the bulk volume of the system (bulk crystallization), and one on the wall surfaces (surface crystallization). Depending on volume fraction initially amorphous (disordered) hard-sphere chain packings transit to the stable crystal phase. The established ordered morphologies consist primarily of hexagonal close-packed (HCP) crystals in the bulk volume and of triangular (TRI) crystals on the surface. As in the case of athermal packings in the bulk (without confinement), a structural competition is observed between the 5-fold local symmetry and the formation of close-packed crystallites. Effectively, the full confinement inside a cube favors the growth of the HCP crystal, as the FCC one is quite incompatible with the imposed spatial constraints. Consequently, we observe the formation of noncompact ordered motifs which grow from the surface to the inner volume of the simulation cell. We further compare the 2D and 3D crystals formed by monomeric hard spheres under the same simulation conditions. Significant differences are observed at low densities that tend to diminish as concentration increases.

Densest versus jammed packings of bent-core trimers

We identify putatively maximally dense packings of tangent-sphere trimers with fixed bond angles (θ=θ_{0}), and contrast them to the disordered jammed states they form under quasistatic and dynamic athermal compression. Incommensurability of θ_{0} with three-dimensional (3D) close packing does not by itself inhibit formation of dense 3D crystals; all θ_{0} allow formation of crystals with ϕ_{max}(θ_{0})>0.97ϕ_{cp}. Trimers are always able to arrange into periodic structures composed of close-packed bilayers or trilayers of triangular-lattice planes, separated by "gap layers" that accommodate the incommensurability. All systems have ϕ_{J} significantly below the monomeric value, indicating that trimers' quenched bond-length and bond-angle constraints always act to promote jamming. ϕ_{J} varies strongly with θ_{0}; straight (θ_{0}=0) trimers minimize ϕ_{J} while closed (θ_{0}=120^{∘}) trimers maximize it. Marginally jammed states of trimers with lower ϕ_{J}(θ_{0}) exhibit quantifiably greater disorder, and the lower ϕ_{J} for small θ_{0} is apparently caused by trimers' decreasing effective configurational freedom as they approach linearity.

Effects of spherical confinement and backbone stiffness on flexible polymer jamming

We use molecular simulations to study jamming of a crumpled bead-spring model polymer in a finite container and compare to jamming of repulsive spheres. After proper constraint counting, the onset of rigidity is seen to occur isostatically as in the case of repulsive spheres. Despite this commonality, the presence of the curved container wall and polymer backbone bonds introduce new mechanical properties. Notably, these include additional bands in the vibrational density of states that reflect the material structure as well as oscillations in local contact number and density near the wall but with lower amplitude for polymers. Polymers have fewer boundary contacts, and this low-density surface layer strongly reduces the global bulk modulus. We further show that bulk-modulus dependence on backbone stiffness can be described by a model of stiffnesses in series and discuss potential experimental and biological applications.

Different bending models predict different dynamics of sedimenting elastic trumbbells

The main goal of this paper is to examine theoretically and numerically the impact of a chosen bending model on the dynamics of elastic filaments settling in a viscous fluid under gravity at low-Reynolds-number. We use the bead-spring approximation of a filament and the Rotne-Prager mobility matrix to describe hydrodynamic interactions between the beads. We analyze the dynamics of trumbbells, for which bending angles are typically larger than for thin and long filaments. Each trumbbell is made of three beads connected by springs and it exhibits a bending resistance, described by the harmonic or - alternatively - by the 'cosine' (also called the Kratky-Porod) bending models, both often used in the literature. Using the harmonic bending potential, and coupling it to the spring potential by the Young's modulus, we find simple benchmark solutions: stable stationary configurations of a single elastic trumbbell and attraction of two elastic trumbbells towards a periodic long-lasting orbit. As the most significant result of this paper, we show that for very elastic trumbbells at the same initial conditions, the Kratky-Porod bending potential can lead to qualitatively and quantitatively different spurious dynamics, with artificially large bending angles and unrealistic shapes. We point out that for the bead models of an elastic filament, the range of applicability of the Kratky-Porod model might not go beyond bending angles smaller than π/2 for touching beads and beyond an even much lower value for beads well-separated from each other. The existence of stable stationary configurations of elastic trumbbells and a family of periodic oscillations of two elastic trumbbells are very important findings on their own.