Study of the structure and dynamics of poly(vinyl pyrrolidone) by molecular dynamics simulations validated by quasielastic neutron scattering and x-ray diffraction experiments

Disentangling Self-Atomic Motions in Polyisobutylene by Molecular Dynamics Simulations

We present fully atomistic molecular dynamics simulations on polyisobutylene (PIB) in a wide temperature range above the glass transition. The cell is validated by direct comparison of magnitudes computed from the simulation and measured by neutron scattering on protonated samples reported in previous works. Once the reliability of the simulation is assured, we exploit the information in the atomic trajectories to characterize the dynamics of the different kinds of atoms in PIB. All of them, including main-chain carbons, show a crossover from Gaussian to non-Gaussian behavior in the intermediate scattering function that can be described in terms of the anomalous jump diffusion model. The full characterization of the methyl-group hydrogen motions requires accounting for rotational motions. We show that the usually assumed statistically independence of rotational and segmental motions fails in this case. We apply the rotational rate distribution model to correlation functions calculated for the relative positions of methyl-group hydrogens with respect to the carbon atom at which they are linked. The contributions to the vibrational density of states are also discussed. We conclude that methyl-group rotations are coupled with the main-chain dynamics. Finally, we revise in the light of the simulations the hypothesis and conclusions made in previously reported neutron scattering investigations on protonated samples trying to address the origin of the dielectric β-process.

Insight into the Structure and Dynamics of Polymers by Neutron Scattering Combined with Atomistic Molecular Dynamics Simulations

Combining neutron scattering and fully atomistic molecular dynamics simulations allows unraveling structural and dynamical features of polymer melts at different length scales, mainly in the intermolecular and monomeric range. Here we present the methodology developed by us and the results of its application during the last years in a variety of polymers. This methodology is based on two pillars: (i) both techniques cover approximately the same length and time scales and (ii) the classical van Hove formalism allows easily calculating the magnitudes measured by neutron scattering from the simulated atomic trajectories. By direct comparison with experimental results, the simulated cell is validated. Thereafter, the information of the simulations can be exploited, calculating magnitudes that are experimentally inaccessible or extending the parameters range beyond the experimental capabilities. We show how detailed microscopic insight on structural features and dynamical processes of various kinds has been gained in polymeric systems with different degrees of complexity, and how intriguing questions as the collective behavior at intermediate length scales have been faced.

Investigation of the dynamics of aqueous proline solutions using neutron scattering and molecular dynamics simulations

We applied quasielastic neutron scattering (QENS) techniques to samples with two different contrasts (deuterated solute/hydrogenated solvent and the opposite label) to selectively study the component dynamics of proline/water solutions. Results on diluted and concentrated solutions (31 and 6 water molecules/proline molecule, respectively) were analyzed in terms of the susceptibility and considering a recently proposed model for water dynamics [Arbe et al., Phys. Rev. Lett., 2016, 117, 185501] which includes vibrations and the convolution of localized motions and diffusion. We found that proline molecules not only reduce the average diffusion coefficient of water but also extend the time/frequency range of the crossover region ('cage') between the vibrations and purely diffusive behavior. For the high proline concentration we also found experimental evidence of water heterogeneous dynamics and a distribution of diffusion coefficients. Complementary molecular dynamics simulations show that water molecules start to perform rotational diffusion when they escape the cage regime but before the purely diffusive behavior is established. The rotational diffusion regime is also retarded by the presence of proline molecules. On the other hand, a strong coupling between proline and water diffusive dynamics which persists with decreasing temperature is directly observed using QENS. Not only are the temperature dependences of the diffusion coefficients of both components the same, but their absolute values also approach each other with increasing proline concentration. We compared our results with those reported using other techniques, in particular using dielectric spectroscopy (DS). A simple approach based on molecular hydrodynamics and a molecular treatment of DS allows rationalizing the a priori puzzling inconsistency between QENS and dielectric results regarding the dynamic coupling of the two components. The interpretation proposed is based on general grounds and therefore should be applicable to other biomolecular solutions.

Molecular ring rotation in poly(vinylferrocene)

We investigate the ring rotation dynamics in poly(vinylferrocene) (PVFc) using incoherent neutron spectroscopy. PVFc contains ferrocene units laterally attached to a polymer backbone, allowing for one cyclopentadienyl ring of the organometallic sandwich structure of ferrocene to undergo rotational jump diffusion. The barrier of rotation is found to be broadly distributed, but the dynamics can be well described using a rotation rate distribution model which is well known from the description of methyl group rotation in glassy polymers. As necessary information for the analysis of quasielastic scattering data, we measure the static structure factor of the polymer using polarized neutron diffraction. Neutron time-of-flight and backscattering data are then combined and consistently modeled over the large temperature range from 80 K to 350 K yielding an Arrhenius behavior of the jump rate distribution. The mean value of potential barrier distribution is found to be 〈E_A〉 = 9.61(2) kJ mol^-1 with a root mean square width of σ_E = 3.12(1) kJ mol^-1, being the result of superposition of constant intramolecular and heterogeneous intermolecular rotational barriers.

Dynamics of tetrahydrofuran as minority component in a mixture with poly(2-(dimethylamino)ethyl methacrylate): A neutron scattering and dielectric spectroscopy investigation

We have investigated a mixture of poly(2-(dimethylamino)ethyl methacrylate) (PDMAEMA) and tetrahydrofuran (THF) (70 wt. % PDMAEMA/30 wt. % THF) by combining dielectric spectroscopy and quasielastic neutron scattering (QENS) on a labelled sample, focusing on the dynamics of the THF molecules. Two independent processes have been identified. The "fast" one has been qualified as due to an internal motion of the THF ring leading to hydrogen displacements of about 3 Å with rather broadly distributed activation energies. The "slow" process is characterized by an Arrhenius-like temperature dependence of the characteristic time which persists over more than 9 orders of magnitude in time. The QENS results evidence the confined nature of this process, determining a size of about 8 Å for the volume within which THF hydrogens' motions are restricted. In a complementary way, we have also investigated the structural features of the sample. This study suggests that THF molecules are well dispersed among side-groups nano-domains in the polymer matrix, ruling out a significant presence of clusters of solvent. Such a good dispersion, together with a rich mobility of the local environment, would prevent cooperativity effects to develop for the structural relaxation of solvent molecules, frustrating thereby the emergence of Vogel-Fulcher-like behavior, at least in the whole temperature interval investigated.