Arturo Narros

Fluids of semiflexible ring polymers: effective potentials and clustering

We present a computational investigation of the structural properties of a fluid of semiflexible ring polymers. Stiffness is introduced by implementing intramolecular barriers. Because of these barriers, shrinkage of the rings is energetically unfavourable, and the ring size can exhibit a non-monotonic density dependence. At high concentrations the rings can swell and adopt open configurations that facilitate interpenetration and clustering. We obtain effective potentials between the centers-of-mass of the rings at infinite dilution, and explore their validity over the whole range of concentrations. Except for the limit of small rings, the effective fluid of ultrasoft particles provides a good description of the real system over a considerable range of densities, even above the overlap concentration. In particular the clustering behaviour predicted by the effective description is observed in the real system for a certain range of molecular masses. However, the effective description is incomplete. Inspection of the clusters of real rings reveals that these can arrange in a complex disordered phase formed by long columns of oblate rings, which are penetrated by bundles of elongated prolate rings. These anisotropic features of the real system are not captured by the standard effective approach, which only considers macromolecular centers-of-mass. This suggests the need to include the relative orientation between rings in the effective potentials.

Effects of Knots on Ring Polymers in Solvents of Varying Quality

We employ extensive computer simulations to investigate the conformations and the interactions of ring polymers under conditions of worsening solvent quality, in comparison with those for linear polymers. We determine the dependence of the Θ-temperature on knotedness by considering ring polymers of different topologies. We establish a clear decrease of the former upon changing the topology of the polymer from linear to an unknotted ring and a further decrease of the same upon introducing trefoil- or 5-fold knots but we find no difference in the Θ-point between the two knotted molecules. Our results are based on two independent methods: one considering the scaling of the gyration radius with molecular weight and one based on the dependence of the effective interaction on solvent quality. In addition, we calculate several shape-parameters of the polymers to characterize linear, unknotted, and knotted topologies in good solvents and in the proximity of the Θ-point. The shape parameters of the knotted molecules show an interesting crossover at a degree of polymerization that depends on the degree of knottedness of the molecule.

Influence of topology on effective potentials: coarse-graining ring polymers

We employ computer simulations and integral equation theory techniques to perform coarse-graining of self-avoiding ring polymers with different knotedness and to derive effective interaction potentials between the centers of mass of these macromolecular entities. Different microscopic models for the monomer-monomer interactions and bonding are employed, bringing about an insensitivity of the effective interactions on the microscopic details and a convergence to a universal form for sufficiently long molecules. The pair effective interactions are shown to be accurate up to within the semidilute regime with additional, many-body forces becoming increasingly important as the polymer concentration grows. The dramatic effects of topological constraints in the form of interaction potentials are brought forward and critically discussed.

Hydrogen motions in the α-relaxation regime of poly(vinyl ethylene): A molecular dynamics simulation and neutron scattering study

The hydrogen motion in poly(vinyl ethylene) (1,2-polybutadiene) in the alpha-relaxation regime has been studied by combining neutron spin echo (NSE) measurements on a fully protonated sample and fully atomistic molecular dynamics simulations. The almost perfect agreement between experiment and simulation results validates the simulated cell. A crossover from Gaussian to non-Gaussian behavior is observed for the intermediate scattering function obtained from both NSE measurements and simulations. This crossover takes place at unusually low Q values, well below the first maximum of the static structure factor. Such anomalous deviation from Gaussian behavior can be explained by the intrinsic dynamic heterogeneity arising from the differences in the dynamics of the different protons in this system. Side group hydrogens show a markedly higher mobility than main chain protons. Taking advantage of the simulations we have investigated the dynamic features of all different types of hydrogens in the sample. Considering each kind of proton in an isolated way, deviations from Gaussian behavior are also found. These can be rationalized in the framework of a simple picture based on the existence of a distribution of discrete jumps underlying the atomic motions in the alpha process.

Cluster glasses of semiflexible ring polymers

We present computer simulations of concentrated solutions of unknotted nonconcatenated semiflexible ring polymers. Unlike in their flexible counterparts, shrinking involves a strong energetic penalty, favoring interpenetration and clustering of the rings. We investigate the slow dynamics of the centers-of-mass of the rings in the amorphous cluster phase, consisting of disordered columns of oblate rings penetrated by bundles of prolate ones. Scattering functions reveal a striking decoupling of self- and collective motions. Correlations between centers-of-mass exhibit slow relaxation, as expected for an incipient glass transition, indicating the dynamic arrest of the cluster positions. However, self-correlations decay at much shorter time scales. This feature is a manifestation of the fast, continuous exchange and diffusion of the individual rings over the matrix of clusters. Our results reveal a novel scenario of glass formation in a simple monodisperse system, characterized by self-collective decoupling, soft caging, and mild dynamic heterogeneity.

Atomic motions in the αβ-merging region of 1,4-polybutadiene: A molecular dynamics simulation study

We present fully atomistic molecular dynamics simulations on 1,4-polybutadiene in a wide temperature range from 200 to 280 K, i.e., in the region where the alpha- and beta-relaxations merge and above. A big computational effort has been performed-especially for the lowest temperatures investigated-to extend the simulation runs to very long times (up to 1 mus for 200 K). The simulated sample has been carefully validated by using previous neutron scattering data on the real sample with similar microstructure. Inspecting the trajectories of the different hydrogens in real space, we have observed a heterogeneous dynamical behavior (each kind of hydrogen moves in a different way) with signatures of combined hopping and diffusive motions in the whole range investigated. The application of a previously proposed model [Colmenero et al., Europhys. Lett. 71, 262 (2005)] is successful and a characterization of the local motions and diffusion is possible. The comparison of our results to those reported in the literature provides a consistent scenario for polybutadiene dynamics and puts into a context the different experimental observations. We also discuss the impact of the hopping processes on the observation and interpretation of experimentally accessible magnitudes and the origin of the deviations from Gaussian behavior in this system.

Hydrogen motions and the α-relaxation in glass-forming polymers: Molecular dynamics simulation and quasi-elastic neutron scattering results

The combination of molecular dynamics simulations and neutron scattering measurements on three different glass-forming polymers (polyisoprene, poly(vinyl ethylene) and polybutadiene) has allowed to establish the existence of a crossover from Gaussian to non-Gaussian behavior for the incoherent scattering function in the α-relaxation regime. The deviation from Gaussian behavior observed can be reproduced assuming the existence of a distribution of discrete jump lengths underlying the sublinear diffusion of the atomic motions during the structural relaxation.

Architecture-Induced Size Asymmetry and Effective Interactions of Ring Polymers: Simulation and Theory

We investigate, by means of Monte Carlo simulations, the role of ring architecture and topology on the relative sizes of two interacting polymers as a function of the distance between their centers-of-mass. As a general rule, polymers swell as they approach each other, irrespectively of their topologies. For each mutual separation, two identical linear polymers adopt the same average size. However, unknotted rings at close separations adopt different sizes, with the small one being “nested” within the large one over long time intervals, exchanging their roles in the course of the simulation. For two rings of different architectures and identical polymerization degree, the knotted one is always smaller, penetrating the unknotted one. On the basis of these observations, we propose a phenomenological theory for the effective interactions between rings, modeling them as unequal-sized penetrable spheres. This simple approximation provides a good description of the simulation results. In particular, it rationalizes the non-Gaussian shape and the short-distance plateau observed in the effective potential between unknotted ring polymers and pairs of unequal-sized unknotted/knotted ones. Our results demonstrate the crucial role of the architecture on both the effective interactions and the molecular size for strongly interpenetrating polymers.

Effective interactions of knotted ring polymers.

In the present article, we review recent computational investigations on the properties of ring polymers in solution. In particular, we focus on effective interactions obtained by means of coarse-graining techniques. We discuss the relative importance of the self-avoidance and the topological contributions in the qualitative features of the effective potential. We extend our previous results on identical rings and determine the effective potential between dissimilar ring polymers of distinct topology and size. The results obtained reveal the dramatic effects of the specific topology on the effective interactions, and hence in the structural correlations, of polymeric systems.

Self‐Atomic Motions in Glass‐Forming Polymers: Neutron Scattering and Molecular Dynamics Simulations Results

By a combined effort of neutron scattering and molecular dynamics simulation, we have recently shown the existence of a crossover from Gaussian to non‐Gaussian character of the self‐correlation function of hydrogen atoms in the α‐relaxation regime of a glass forming polymer, polyisoprene (PI). Following these previous results, here we present new data displaying different features of the above mentioned crossover. In particular, neutron scattering and molecular dynamics simulations corresponding to another glass‐forming polymer, poly(vinyl ethylene) (PVE), and simulation results about the effect of temperature and density on the crossover in PI. We found that PVE displays a similar behavior to PI, showing the generality of the crossover picture. Moreover, we have also found that for PI the main features of the crossover do not depend on temperature and density. These results can be rationalized in the framework of the anomalous jump diffusion model which was previously introduced by us for self atomic mot...