skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Modeling of Entangled Polymer Diffusion in Melts and Nanocomposites: A Review
This review concerns modeling studies of the fundamental problem of entangled (reptational) homopolymer diffusion in melts and nanocomposite materials in comparison to experiments. In polymer melts, the developed united atom and multibead spring models predict an exponent of the molecular weight dependence to the polymer diffusion very similar to experiments and the tube reptation model. There are rather unexplored parameters that can influence polymer diffusion such as polymer semiflexibility or polydispersity, leading to a different exponent. Models with soft potentials or slip-springs can estimate accurately the tube model predictions in polymer melts enabling us to reach larger length scales and simulate well entangled polymers. However, in polymer nanocomposites, reptational polymer diffusion is more complicated due to nanoparticle fillers size, loading, geometry and polymer-nanoparticle interactions.  more » « less
Award ID(s):
1706014
PAR ID:
10178856
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Polymers
Volume:
11
Issue:
5
ISSN:
2073-4360
Page Range / eLocation ID:
876
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The dynamics of polymer melts at the crossover between unentangled and entangled regimes is formalized here through an extension of the Cooperative Dynamics Generalized Langevin Equation (CDGLE) (J. Chem. Phys. 1999, 110, 7574), by including the constraint to the dynamics due to entanglements through an effective intermonomer potential that confines the motion of the chains. As one polymer chain in a melt interpenetrates with a N other chains, with N the degree of chain polymerization, their dynamics is coupled through their potential of mean-force, leading to chains’ cooperative motion and center-of-mass subdiffusive dynamics. When increasing the degree of polymerization, the extended CDGLE approach describes the dynamical behavior of unentangled to weakly entangled systems undergoing cooperative dynamics. By direct comparison of the CDGLE with data of Neutron Spin Echo (NSE) experiments on polyethylene melts, we find that the cooperative dynamics in entangled systems are confined in the region delimited by entanglements. We extend the CDGLE to describe linear dynamical mechanical measurements and use it to calculate shear relaxation for the polyethylene samples investigated by NSE. The effects of cooperative dynamics, local flexibility, and entanglements in the shear relaxation are discussed. It is noteworthy that the theoretical approach describes with accuracy the crossover from unentangled to entangled-global dynamics for polyethylene melts of increasing chain length, covering the regimes of unentangled and weakly entangled (up to 12 entanglements) dynamics in one approach. 
    more » « less
  2. Self-assembled metal nanoparticle-polymer nanocomposite particles as nanoreactors are a promising approach for performing liquid phase reactions using water as a bulk solvent. In this work, we demonstrate rapid, scalable self-assembly of metal nanoparticle catalyst-polymer nanocomposite particles via Flash NanoPrecipitation. The catalyst loading and size of the nanocomposite particles can be tuned independently. Using nanocomposite particles as nanoreactors and the reduction of 4-nitrophenol as a model reaction, we study the fundamental interplay of reaction and diffusion. The induction time is affected by the sequence of reagent addition, time between additions, and reagent concentration. Combined, our experiments indicate the induction time is most influenced by diffusion of sodium borohydride. Following the induction time, scaling analysis and effective diffusivity measured using NMR indicate that the observed reaction rate are reaction- rather than diffusion-limited. Furthermore, the intrinsic kinetics are comparable to ligand-free gold nanoparticles. This result indicates that the polymer microenvironment does not de-activate or block the catalyst active sites. 
    more » « less
  3. null (Ed.)
    Diffusional isotope fractionation has been widely used to explain lithium (Li) isotope variations in minerals and rocks. Isotopic mass dependence of Li diffusion can be empirically expressed as , where is the diffusivity of a Li isotope. The knowledge about temperature and compositional dependence of the factor which is essential for understanding diffusion profiles and mechanisms remains unclear. Based on the potential energy and interatomic forces generated by deep neural networks trained with ab initio data, we performed deep potential molecular dynamics (DPMD) simulations of several Li pseudo-isotopes (with mass = 2, 7, 21, 42 g/mol) in albite, hydrous albite, and model basalt melts to evaluate the factor. Our calculated diffusivities for 7Li in albite and model basalt melts at 1800 K compare well with experimental results. We found that in albite melt decreases from at 4000 K to at 1800 K. The presence of water appears to slightly weaken the temperature dependence of , with decreasing from to in hydrous albite melt. The calculated in model basalt melt takes much smaller values, decreasing from at 4000 K to at 1800 K. Our prediction of in albite and hydrous albite melts is in good agreement with experimental data. More importantly, our results suggest that Li isotope diffusion in silicate melts is strongly dependent on melt composition. The temperature and compositional effects on can be qualitatively explained in terms of ionic porosity and the coupled relationship between Li diffusion and the mobility of the silicate melt network. Two types of diffusion experiments are suggested to test our predicted temperature and compositional dependence of . This study shows that DPMD is a promising tool to simulate the diffusion of elements and isotopes in silicate melts. 
    more » « less
  4. Both entangled and unentangled polymer melts exhibit stress overshoots when subject to shearing flow. The size of the overshoot depends on the applied shear rate and is related to relaxation mechanisms such as reptation, chain stretch, and convective constraint release. Previous experimental work shows that melts subjected to interrupted shear flows exhibit a smaller overshoot when sheared after partial relaxation. This has been shown to be consistent with predictions by constitutive models. Here, we report molecular dynamics simulations of interrupted shear of polymer melts where the shear flow after the relaxation stage is orthogonal to the originally applied flow. We observe that, for a given relaxation time, the size of the stress overshoot under orthogonal interrupted shear is larger than observed during parallel interrupted shear, which is not captured by constitutive models. Differences in maxima are also observed for overshoots in the first normal stress and chain end-to-end distance. We also show that measurements of the average number of entanglements per chain and average orientation at different scales along the chain are affected by the change in shear direction, leading to nonmonotonic relaxation of the off-diagonal components of orientation and an appearance of a “double peak” in the average number of entanglements during the transient. We propose that such complex behavior of entanglements is responsible for the increase in the overshoots of stress components and that models of the dynamics of entanglements might be improved upon by considering a tensorial measurement of entanglements that can be coupled to orientation. 
    more » « less
  5. Eighteen successful diffusion couple experiments in 8-component SiO2–TiO2–Al2O3–FeO–MgO–CaO–Na2O–K2O basaltic melts were conducted at 1260°C and 0.5 GPa and at 1500°C and 1.0 GPa. These experiments are combined with previous data at 1350°C and 1.0 GPa (Guo and Zhang, 2018) to study the temperature dependence of multicomponent diffusion in basaltic melts. Effective binary diffusion coefficients of components with monotonic diffusion profiles were extracted and show a strong dependence on their counter-diffusing component even though the average (or interface) compositions are the same. The diffusion matrix at 1260°C was obtained by simultaneously fitting diffusion profiles of all diffusion couple experiments as well as appropriate data from the literature. All features of concentration profiles in both diffusion couples and mineral dissolution are well reproduced by this new diffusion matrix. At 1500°C, only diffusion couple experiments are used to obtain the diffusion matrix. Eigenvectors of the diffusion matrix are used to discuss the diffusion (exchange) mechanism, and eigenvalues characterize the diffusion rate. Diffusion mechanisms at both 1260 and 1500°C are inferred from eigenvectors of diffusion matrices and compared with those at 1350°C reported in Guo and Zhang (2018). There is indication that diffusion eigenvectors in basaltic melts do not depend much on temperature, but complexity is present for some eigenvectors. The two slowest eigenvectors involve the exchange of SiO2 and/or Al2O3 with nonalkalis. The third slowest eigenvector is due to the exchange of divalent oxides with other oxides. The fastest eigenvector is due to the exchange of Na2O with other oxide components. Some eigenvalues differ from each other by less than 1/3, and their eigenvectors are less well defined. We define small difference in eigenvalues as near degeneracy. In strict mathematical degeneracy, eigenvectors are not uniquely defined because any linear combination of two eigenvectors is also an eigenvector. In the case of near degeneracy, more constraints either in terms of higher data quality or more experiments are needed to resolve the eigenvectors. We made a trial effort to generate a set of average eigenvectors, which are assumed to be constant as temperature varies. The corresponding eigenvalues are roughly Arrhenian. Thus, the temperature-dependent diffusion matrix can be roughly predicted. The method is applied to predict experimental diffusion profiles in basaltic melts during olivine and anorthite dissolution at ~1400°C with preliminary success. We further applied our diffusion matrix to investigate multicomponent diffusion during magma mixing in the Bushveld Complex and found such diffusion may result in an increased likelihood of sulfide and Fe-Ti oxide mineralization. 
    more » « less