Modeling the Effect of Polymer Chain Stiffness on the Behavior of Polymer Nanocomposites

Controlling sodium chloride concentration modulates the supramolecular structure and sol features of wheat starch-acetylated starch binary matrix

The sol performance of wheat starch (WS) matrix incorporating acetylated starch (AS) is crucial for the processing and quality features of wheat products. From a supramolecular structure view, how regulating salt (sodium chloride) concentration modulates the sol features, e.g., pasting, zero-shear viscosity (ZSV) and thixotropy of WS-AS binary matrix was explored. Compared to the salt-free counterpart, the saline matrices exhibited a delayed pasting profile and a decreased viscoelasticity. Thereinto, the sol at 0.02 M NaCl exhibited the smallest ZSV (23,710 Pa·s) and the greatest in-shear recovery ratio (33.7 %). Such variations could be attributed to the weakened coil-helix, nematic-smectic and isotropy-anisotropy transitions from a side-chain liquid-crystalline perspective. Meanwhile, the correlation length (ξ) and radius of gyration (Rg) obtained from small angle X-ray scattering analysis were increased by 5.2 and 9.6 Å respectively, which disclosed a restrained entanglement and an enhanced chain mobility. These results would provide a reference for the design of fluid/semisolid products with optimized qualities.

Attractive crowding effect on passive and active polymer looping kinetics

Loop formation in complex environments is crucially important to many biological processes in life. In the present work, we adopt three-dimensional Langevin dynamics simulations to investigate passive and active polymer looping kinetics in crowded media featuring polymer-crowder attraction. We find polymers undergo a remarkable coil-globule-coil transition, highlighted by a marked change in the Flory scaling exponent of the gyration radius. Meanwhile, looping time as a function of the crowder's volume fraction demonstrates an apparent non-monotonic alteration. A small number of crowders induce a compact structure, which largely facilitates the looping process. While a large number of crowders heavily impede end-to-end diffusion, looping kinetics is greatly inhibited. For a self-propelled chain, we find that the attractive crowding triggers an unusual activity effect on looping kinetics. Once a globular state is formed, activity takes an effort to open the chain from the compact structure, leading to an unexpected activity-induced inhibition of looping. If the chain maintains a coil state, the dominant role of activity is to enhance diffusivity and, thus, speed up looping kinetics. The novel conformational change and looping kinetics of both passive and active polymers in the presence of attractive crowding highlight a rather distinct scenario that has no analogy in a repulsive crowding counterpart. The underlying mechanism enriches our understanding of the crucial role of attractive interactions in modulating polymer structure and dynamics.

The Role of Polymer Chain Stiffness and Guest Nanoparticle Loading in Improving the Glass Transition Temperature of Polymer Nanocomposites

The impact of polymer chain stiffness characterized by the bending modulus (k_θ) on the glass transition temperature (T_g) of pure polymer systems, as well as polymer nanocomposites (PNCs), is investigated using molecular dynamics simulations. At small k_θ values, the pure polymer system and respective PNCs are in an amorphous state, whereas at large k_θ values, both systems are in a semicrystalline state with a glass transition at low temperature. For the pure polymer system, T_g initially increases with k_θ and does not change obviously at large k_θ. However, the T_g of PNCs shows interesting behaviors with the increasing volume fraction of nanoparticles (f_NP) at different k_θ values. T_g tends to increase with f_NP at small k_θ, whereas it becomes suppressed at large k_θ.

Understanding the balance between additives' miscibility and plasticisation effect in polymer composites: a computational study

Plasticisers are small organic molecules routinely added to polymer composites that modify the processability of the compounds by adsorbing on the filler's surface or dispersing into the polymer matrix. Here using a simple yet chemically specific coarse-grained model, we demonstrate that the filler surface coverage and the degree of dispersion into the polymer matrix can be tuned without modifying the chemistry of the plasticisers but only by varying their conformational flexibility. We show that when the adsorption mechanism and clustering into the bulk are entropically driven as in this work, this is a general phenomenon independent on the polymer chemistry and its molecular weight. Our findings suggest a simple practical design rule that requires only minor modifications of the plasticisers' chemistry to achieve maximum adsorption onto the filler surface or dispersion into the polymer matrix.

Dispersion and Diffusion Mechanism of Nanofillers with Different Geometries in Bottlebrush Polymers: Insights from Molecular Dynamics Simulation

The dispersion and diffusion mechanism of nanofillers in polymer nanocomposites (PNCs) are crucial for understanding the properties of PNCs, which is of great significance for the design of novel materials. Herein, we investigate the dispersion and diffusion behavior of two geometries of nanofillers, namely, spherical nanoparticles (SNPs) and nanorods (NRs), in bottlebrush polymers by utilizing coarse-grained molecular dynamics simulations. With the increase of the interaction strength between the nanofiller and polymer (ε_np), both the SNPs and NRs experience a typical "aggregated phase-dispersed phase-bridged phase" state transition in the bottlebrush polymer matrix. We evaluate the validity of the Stokes-Einstein (SE) equation for predicting the diffusion coefficient of nanofillers in bottlebrush polymers. The results demonstrate that the SE predictions are slightly larger than the simulated values for small SNP sizes because the local viscosity that is felt by small SNPs in the densely grafted bottlebrush polymer does not differ much from the macroscopic viscosity. The relative size of the length of the NRs (L) and the radius of gyration (R_g) of the bottlebrush polymer play a key role in the diffusion of NRs. In addition, we characterize the anisotropic diffusion of NRs to analyze their translational and rotational diffusion. The motion of NRs in the direction perpendicular to the end-to-end vector is more hindered, indicating that there is a strong coupling between the rotation of NRs and the motion of the polymer. The NR motion shows stronger anisotropic diffusion at short time scales because of the steric effects generated by side chains of the bottlebrush polymer. In general, our results provide a fundamental understanding of the dispersion of nanofillers and the microscopic mechanism of nanofiller diffusion in bottlebrush polymers.

Passive and active tracer dynamics in polymer solutions with isotropic-to-nematic phase transition

Macromolecular crowding plays a crucial role in determining the dynamics in a living cell. We adopt Langevin dynamics simulations to investigate the anomalous diffusion dynamics of passive and active particles in a solution of polymer chains with tunable stiffness. The solution's anisotropic feature is modulated by changing both the polymer stiffness and volume fraction, where isotropic-to-nematic phase transition is involved. Our results demonstrate the significant impact of polymer flexibility on the dynamics of both passive and active probes. The distinct diffusion mechanism for an active particle is clarified by the interplay between polymer stiffness, crowdedness and activity. Polymer stiffness leads to a global inhibition effect on passive particle diffusion. The diffusion coefficient exhibits an intriguing non-monotonic variation at increasing polymer stiffness, which is due to the fact that the alignment of polymer chains is beneficial for diffusion along the nematic direction but unfavorable for that in the direction perpendicular to it. In sharp contrast, polymer stiffness plays a dominant role in facilitating active particle diffusion. Self-propulsion of the particle can utilize stiffness-induced elastic interactions more efficiently, which promotes its mobility in both directions. Meanwhile, an active particle might have a stronger ability to take advantage of the polymer alignment, contributing substantially enhanced diffusivity. In addition, the diffusion coefficient of an active particle is subject to a tendency of degeneration against varying volume fraction. This counter-intuitive behavior is due to the contrasting factors that increasing crowdedness induces a lower particle speed but a longer persistent motion time.

Local structuring of diketopyrrolopyrrole (DPP)-based oligomers from molecular dynamics simulations

The microscopic structure of high mobility semiconducting polymers is known to be essential for their performance but it cannot be easily deduced from the available experimental data. A series of short oligomers of diketopyrrolopyrrole (DPP)-based materials that display high charge mobility are studied by molecular dynamics simulations to understand their local structuring at an atomic level. Different analyses are proposed to compare the ability of different oligomers to form large aggregates and their driving force. The simulations show that the tendency for this class of materials to form aggregates is driven by the interaction between DPP fragments, but this is modulated by the other conjugated fragments of the materials which affect the rigidity of the polymer and, ultimately, the size of the aggregates that are formed. The main structural features and the electronic structure of the oligomers are fairly similar above the glass transition temperature and at room temperature.

Glass transition and dynamics of semiflexible polymer brushes

The glass transition and dynamics of densely grafted semiflexible polymer brushes are studied by molecular dynamics simulation. The glass transition temperature (T_g) increases with the polymer rigidity. The local glass transition temperature (T_g,local) is estimated from the temperature-dependent dynamics of individual segments including the lateral position fluctuation and lateral mean square displacement. Different from the flexible polymer brush, T_g,local of semiflexible polymer brushes is roughly independent of the segment height. Our simulation reveals that the glass transition is in synchronism with an abrupt change of the chain conformation in semiflexible polymer brushes. When the temperature drops to near T_g, the semiflexible polymer chains elongate, tilt, and become more ordered. Moreover, enhanced segmental dynamics is observed at temperatures just above T_g for the semiflexible polymer brushes.

Molecular Dynamics of Janus Nanodimers Dispersed in Lamellar Phases of a Block Copolymer

We investigate structural and dynamical properties of Janus nanodimers (NDs) dispersed in lamellar phases of a diblock copolymer. By performing molecular dynamics simulations, we show that an accurate tuning of the interactions between NDs and copolymer blocks can lead to a close control of NDs' space distribution and orientation. In particular, NDs are preferentially found within the lamellae if enthalpy-driven forces offset their entropic counterpart. By contrast, when enthalpy-driven forces are not significant, the distribution of NDs, preferentially observed within the inter-lamellar spacing, is mostly driven by excluded-volume effects. Not only does the degree of affinity between host and guest species drive the NDs' distribution in the polymer matrix, but it also determines their space orientation. In turn, these key structural properties influence the long-time dynamics and the ability of NDs to diffuse through the polymer matrix.

Theory of polymer diffusion in polymer-nanoparticle mixtures: effect of nanoparticle concentration and polymer length

The dynamics of polymer-nanoparticle (NP) mixtures, which involves multiple scales and system-specific variables, has posed a long-standing challenge on its theoretical description. In this paper, we construct a microscopic theory for polymer diffusion in mixtures based on a combination of the generalized Langevin equation, mode-coupling approach, and polymer physics ideas. The parameter-free theory has an explicit expression and remains tractable on a pair correlation level with system-specific equilibrium structures as input. Taking a minimal polymer-NP mixture as an example, our theory correctly captures the dependence of polymer diffusion on NP concentration and average interparticle distance. Importantly, the polymer diffusion exhibits a power law decay as the polymer length increases at dense NPs and/or a long chain, which marks the emergence of entanglement-like motion. The work provides a first-principles theoretical foundation to investigate dynamic problems in diverse polymer nanocomposites.

Nanoparticle assembly modulated by polymer chain conformation in composite materials

Mixing nanoparticles into a strategically selected polymer matrix yields nanocomposites with well-controlled microstructures and unique properties and functions. The modulation of nanoparticle assembly by polymer chain conformation can play a dominant role in determining nanocomposite structures, yet such a physical mechanism remains largely unexplored. We hypothesize that highly ordered microdomains of rigid linear polymers provide a template for nanoparticle assembly into open fractal structures. We conducted mesoscopic computer simulations and physical experiments to elucidate how polymer chain conformation regulates the dynamic evolution of nanoparticle structures during the drying processing of polymer nanocomposite films. The evaporation of polymer-nanoparticle mixtures with varying chain stiffnesses was simulated using dissipative particle dynamics. The formation of distinguished nanoparticle assemblies as a result of matrix selection was further corroborated by probing nanoparticle aggregation in different polymer nanocomposite coatings. The results show that polymer conformation not only influences the dispersion states of individual particles (dispersed vs. aggregated), but also modulates the morphologies of large-scale assembly (globular vs. fractal). The emergence of nematically ordered polymer clusters when the chain rigidity is increased creates local solvent-rich "voids" that promote anisotropic particle aggregates, which then percolate into open fractal structures upon solvent evaporation. The nanoparticle dynamics also exhibits an intriguing non-monotonic behavior attributed to the transitions between the coupling and decoupling with polymer dynamics. The nanoparticle assembly morphologies obtained in simulations match well with the electron microscopy images taken in physical experiments.

Performance of Coarse Graining in Estimating Polymer Properties: Comparison with the Atomistic Model

Combining atomistic and coarse-grained (CG) models is a promising approach for quantitative prediction of polymer properties. However, the gaps between the length and time scales of atomistic and CG models still need to be bridged. Here, the scale gaps of the atomistic model of polyethylene melts, the bead–spring Kremer–Grest model, and dissipative particle dynamics with the slip-spring model were investigated. A single set of spatial and temporal scaling factors was determined between the atomistic model and each CG model. The results of the CG models were rescaled using the set of scaling factors and compared with those of the atomistic model. For each polymer property, a threshold value indicating the onset of static or dynamic universality of polymers was obtained. The scaling factors also revealed the computational efficiency of each CG model with respect to the atomistic model. The performance of the CG models of polymers was systematically evaluated in terms of both the accuracy and computational efficiency.

Entropy-induced Separation of Binary Semiflexible Ring Polymer Mixtures in Spherical Confinement

Coarse-grained molecular dynamics simulations are used to investigate the conformations of binary semiflexible ring polymers (SRPs) of two different lengths confined in a hard sphere. Segregated structures of SRPs in binary mixtures are strongly dependent upon the number density of system (ρ), the bending energy of long SRPs (Kb, long), and the chain length ratio of long to short SRPs (α). With a low ρ or a weak Kb, long at a small ratio α, long SRPs are immersed randomly in the matrix of short SRPs. As ρ and bending energy of long SRPs (Kb, long) are increased up to a certain value for a large ratio α, a nearly complete segregation between long and short SRPs is observed, which can be further characterized by the ratio of tangential and radial components of long SRPs velocity. These explicit segregated structures of the two components in spherical confinement are induced by a delicate competition between the entropic excluded volume (depletion) effects and bending contributions.

Monte Carlo simulation on the dynamics of a semi-flexible polymer in the presence of nanoparticles

The dynamics of a semi-flexible polymer chain in the presence of periodically distributed nanoparticles is simulated by using off-lattice Monte Carlo simulations. For repulsive or weak attractive nanoparticles, the dynamics are slowed down monotonically by increasing the chain stiffness k_θ or decreasing the inter-particle distance d. For strong attractive nanoparticles, however, the dynamics show nonmonotonic behaviors with k_θ and d. An interesting result is that a stiff polymer may move faster than a flexible one. The underlying mechanism is that the nanoparticle's attraction is weakened by the chain stiffness. The nonmonotonic behavior of the polymer's dynamics with k_θ is explained by the competition between the weakening effect of the chain stiffness on the nanoparticle's attraction and the intrinsic effect of chain stiffness which reduces the dynamics of the polymer. In addition, the nonmonotonic behavior of the polymer's dynamics with d is explained by the competition between the nanoparticle-exchange motion of the polymer dominated at small d and the desorption-and-adsorption motion at large d. The excluded volume effect of the nanoparticles plays a more important role for stiffer polymers as the attraction of the nanoparticles is weakened by the chain stiffness.

Tailoring the mechanical properties by molecular integration of flexible and stiff polymer networks

Designing a multiple-network structure at the molecular level to tailor the mechanical properties of polymeric materials is of great scientific and technological importance. Through the coarse-grained molecular dynamics simulation, we successfully construct an interpenetrating polymer network (IPN) composed of a flexible polymer network and a stiff polymer network. First, we find that there is an optimal chain stiffness for a single network (SN) to achieve the best stress-strain behavior. Then we turn to study the mechanical behaviors of IPNs. The result shows that the stress-strain behaviors of the IPNs appreciably exceed the sum of that of the corresponding single flexible and stiff network, which highlights the advantage of the IPN structure. By systematically varying the stiffness of the stiff polymer network of the IPNs, optimal stiffness also exists to achieve the best performance. We attribute this to a much larger contribution of the non-bonded interaction energy. Last, the effect of the component concentration ratio is probed. With the increase of the concentration of the flexible network, the stress-strain behavior of the IPNs is gradually enhanced, while an optimized concentration (around 60% molar ration) of the stiff network occurs, which could result from the dominant role of the enthalpy rather than the entropy. In general, our work is expected to provide some guidelines to better tailor the mechanical properties of the IPNs made of a flexible network and a stiff network, by manipulating the stiffness of the stiff polymer network and the component concentration ratio.

Translational and rotational dynamics of an ultra-thin nanorod probe particle in linear polymer melts

Translational and rotational dynamics of a single rigid ultra-thin nanorod probe particle in linear polymer melts are investigated using coarse-grained molecular dynamics (CG-MD) simulations.