Here, wide-angle X-ray scattering, molecular dynamics (MD) simulations, and integral equation theory are used to study the structure of poly(diethylsiloxane) (PDES), poly(ethylmethylsiloxane) (PEMS), and poly(dimethylsiloxane) (PDMS) melts. The structure functions of PDES, PEMS, and PDMS are similar, but systematic trends in the intermolecular packing are observed. The local intramolecular structure is extracted from the experimental structure functions. The bond distances and bond angles obtained, including the large Si-O-Si angle, are in good agreement with the explicit atom (EA) and united atom (UA) potentials used in the simulations and theory and from other sources. Very good agreement is found between the MD simulations using the EA potentials and the experimental scattering results. Good agreement is also found between the polymer reference interaction site model (PRISM theory) and the UA MD simulations. The intermolecular structure is examined experimentally using an appropriately weighted radial distribution function and with theory and simulation using intermolecular site/site pair correlation functions. Finally, experiment, simulation, and theory show systematic increases in the chain/chain packing distances in the siloxanes as the number of sites in the pendant side chains is increased.
We have enhanced our parallel molecular dynamics (MD) simulation software LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator, lammps.sandia.gov) to include many new features for accelerated simulation including articulated rigid body dynamics via coupling to the Rensselaer Polytechnic Institute code POEMS (Parallelizable Open-source Efficient Multibody Software). We use new features of the LAMMPS software package to investigate rhodopsin photoisomerization, and water model surface tension and capillary waves at the vapor-liquid interface. Finally, we motivate the recipes of MD for practitioners and researchers in numerical analysis and computational mechanics.
Using molecular dynamics simulations, a constitutive model for the chemical aging of polymer networks was developed. This model incorporates the effects on the stress from the chemical crosslinks and the physical entanglements. The independent network hypothesis has been modified to account for the stress transfer between networks due to crosslinking and scission in strained states. This model was implemented in the finite element code Adagio and validated through comparison with experiment. Stress relaxation data was used to deduce crosslinking history and the resulting history was used to predict permanent set. The permanent set predictions agree quantitatively with experiment.
The Independent Network Model (INM) has proven to be a useful tool for understanding the development of permanent set in strained elastomers. Our previous work showed the applicability of the INM to our simulations of polymer systems crosslinking in strained states. This study looks at the INM applied to theoretical models incorporating entanglement effects, including Flory's constrained junction model and more recent tube models. The effect of entanglements has been treated as a separate network formed at gelation, with additional curing treated as traditional phantom contributions. Theoretical predictions are compared with large-scale molecular dynamics simulations.
Constitutive models for chemically reacting networks are formulated based on a generalization of the independent network hypothesis. These models account for the coupling between chemical reaction and strain histories, and have been tested by comparison with microscopic molecular dynamics simulations. An essential feature of these models is the introduction of stress transfer functions that describe the interdependence between crosslinks formed and broken at various strains. Efforts are underway to implement these constitutive models into the finite element code Adagio. Preliminary results are shown that illustrate the effects of changing crosslinking and scission rates and history.