Polyurea-Graphene Nanocomposites-The Influence of Hard-Segment Content and Nanoparticle Loading on Mechanical Properties.

Polyurethane and polyurea-based adhesives are widely used in various applications, from automotive to electronics and medical applications. The adhesive performance depends strongly on its composition, and developing the formulation-structure-property relationship is crucial to making better products. Here, we investigate the dependence of the linear viscoelastic properties of polyurea nanocomposites, with an IPDI-based polyurea (PUa) matrix and exfoliated graphene nanoplatelet (xGnP) fillers, on the hard-segment weight fraction (HSWF) and the xGnP loading. We characterize the material using scanning electron microscopy (SEM) and dynamic mechanical analysis (DMA). It is found that changing the HSWF leads to a significant variation in the stiffness of the material, from about 10 MPa for 20% HSWF to about 100 MPa for 30% HSWF and about 250 MPa for the 40% HSWF polymer (as measured by the tensile storage modulus at room temperature). The effect of the xGNP loading was significantly more limited and was generally within experimental error, except for the 20% HSWF material, where the xGNP addition led to about an 80% increase in stiffness. To correctly interpret the DMA results, we developed a new physics-based rheological model for the description of the storage and loss moduli. The model is based on the fractional calculus approach and successfully describes the material rheology in a broad range of temperatures (-70 °C-+70 °C) and frequencies (0.1-100 s-1), using only six physically meaningful fitting parameters for each material. The results provide guidance for the development of nanocomposite PUa-based materials.

Combined Description of the Equation of State and Diffusion Coefficient of Liquid Water Using a Two-State Sanchez-Lacombe Approach.

Water is one of the most important compounds on Earth, yet its material properties are still poorly understood. Here, we use a recently developed two-state, two-(time)scale (TS2) dynamic mean-field model combined with the two-state Sanchez-Lacombe (SL) thermodynamic theory in order to describe the equation of state (density as a function of temperature and pressure) and diffusivity of liquid water. In particular, it is shown that in a relatively wide temperature and pressure range (160 K < T < 360 K; 0 < P < 100 MPa), density and self-diffusion obey a special type of dynamic scaling, similar to the "τTV" scaling of Casalini and Roland, but with the negative exponent γ. The model predictions are consistent with experimental data. The new equation of state can be used for various process models and generalized to include multicomponent mixtures.

Surfactant in a Polyol-CO2 Mixture: Insights from a Classical Density Functional Theory Study.

Silicone-polyether (SPE) surfactants, made of a polydimethyl-siloxane (PDMS) backbone and polyether branches, are commonly used as additives in the production of polymeric foams with improved properties. A key step in the production of polymeric foams is the nucleation of gas bubbles in the polymer matrix upon supersaturation of dissolved gas. However, the role of SPE surfactants in the nucleation of gas bubbles is not well understood. In this study, we use classical density functional theory to investigate the effect of an SPE surfactant on the nucleation of CO2 bubbles in a polyol foam formulation. We find that the addition of an SPE surfactant leads to a ∼3-fold decrease in the polyol-CO2 interfacial tension at the surfactant's critical micelle concentration. Additionally, the surfactant is found to reduce the free energy barrier and affect the minimum free energy pathway (MFEP) associated with CO2 bubble nucleation. In the absence of a surfactant, a CO2-rich bubble nucleates from a homogeneous CO2-supersaturated polyol solution by following an MFEP characterized by a single nucleation barrier. Adding a surfactant results in a two-step nucleation process with reduced free energy barriers. The first barrier corresponds to the formation of a spherical aggregate with a liquid-like CO2 core. This spherical aggregate then grows into a CO2-rich bubble (spherical aggregate with a vapor-like CO2 core) of a critical size representing the second barrier. We hypothesize that the stronger affinity of CO2 for PDMS (than polyether) stabilizes the spherical aggregate with the liquid-like CO2 core, leading to a lower free energy barrier for CO2 bubble nucleation. Stabilization of such an aggregate during the early stages of the nucleation may lead to foams with more, smaller bubbles, which can improve their microstrustural features and insulating abilities.

Combined description of pressure-volume-temperature and dielectric relaxation of several polymeric and low-molecular-weight organic glass-formers using SL-TS2 approach.

We apply our recently-developed mean-field "SL-TS2" (two-state Sanchez-Lacombe) model to simultaneously describe dielectric α-relaxation time, τα, and pressure-volume-temperature (PVT) data in four polymers (polystyrene, poly(methylmethacrylate), poly(vinyl acetate) and poly(cyclohexane methyl acrylate)) and four organic molecular glass formers (ortho-terphenyl, glycerol, PCB-62, and PDE). Previously, it has been shown that for all eight materials, the Casalini-Roland thermodynamical scaling, τα = f(Tvγsp) (where T is temperature and vsp is specific volume) is satisfied (R. Casalini and C. M. Roland, Phys. Rev. E, 2004, 69(6), 62501). It has also been previously shown that the same scaling emerges naturally (for sufficiently low pressures) within the "SL-TS2" framework (V. V. Ginzburg, Soft Matter, 2021, 17, 9094-9106). Here, we fit the ambient pressure curves for the relaxation time and the specific volume as functions of temperature for the eight materials and observe a good agreement between theory and experiment. We then use the Casalini-Roland scaling to convert those results into "master curves", thus enabling predictions of relaxation times and specific volumes at elevated pressures. The proposed approach can be used to describe other glass-forming materials, both low-molecular-weight and polymeric.

Combined description of polymer PVT and relaxation data using a dynamic "SL-TS2" mean-field lattice model.

We develop a combined model to describe the pressure-volume-temperature (PVT) thermodynamics and the α- and ß-relaxation time dynamics in glass-forming amorphous materials. The PVT results are described using a two-state modification of the Sanchez-Lacombe equation of state (SL-EoS). The minimization of the Sanchez-Lacombe free energy expression allows one to calculate the density (or specific volume) and the "solid fraction" as a function of temperature and pressure. The solid fraction, ψ(T,P), is then substituted into the TS2 mean-field model (V. V. Ginzburg, Soft Matter, 2020, 16, 810), and the equilibrium relaxation times and the glass transition temperature are calculated. Finally, we use a dynamic model of a Tool-Narayanaswami-Moynihan (TNM)-type to describe the density relaxation as a function of the cooling rate. We applied the new framework to describe experimental data for polystyrene and poly(methylmethacrylate) and found a good qualitative and quantitative agreement.

Nonelectrostatic Adsorption of Polyelectrolytes and Mediated Interactions between Solid Surfaces.

Polymer-mediated interaction between two solid surfaces is directly connected to the properties of the adsorbed polymer layers. Nonelectrostatic interactions with a surface can significantly impact the adsorption of polyelectrolytes to charged surfaces. We use a classical density functional theory to study the effect of various polyelectrolyte solution properties on the adsorption and interaction between two like-charged surfaces. Our results show that nonelectrostatic interactions not only enhance polyelectrolyte adsorption but can also result in qualitatively different salt effects with respect to the adsorbed amount. In particular, we observe decreasing, increasing, and a previously unreported nonmonotonic behavior in the adsorbed amount of polymer with added salt under the conditions studied, although the nonmonotonic regime only occurs for a narrow range in the parameter space. With sufficient nonelectrostatic adsorption, the adsorbed polymer layers produce a long-range repulsive barrier that is strong enough to overcome dispersive interactions that cause surfaces to attract. Concurrently, a short-range bridging attraction is observed when the two polyelectrolyte layers span both the surfaces. Both the repulsive barrier and bridging attraction depend on the charge density of the polymer backbone and the bulk salt concentration but not on the chain length in the semidilute regime studied.

A simple mean-field model of glassy dynamics and glass transition.

We propose a phenomenological model to describe the equilibrium dynamic behavior of amorphous glassy materials. It is assumed that a material can be represented by a lattice of cooperatively re-arranging regions (CRRs), with each CRR having two states, the low-temperature "solid" and the high-temperature "liquid". At low temperatures, the material exhibits two characteristic relaxation times, corresponding to the slow large-scale motion between the "solid" CRRs (α-relaxation) and the faster local motion within individual CRRs (ß-relaxation). At high temperatures, the α- and ß-relaxation times merge, as observed experimentally and suggested by the "Coupling Model" framework. Our new approach is labeled "Two-state, two (time)scale model" or TS2. It is shown that the TS2 treatment can successfully describe the "two-Arrhenius" relaxation time behavior described in several recent experiments. We also apply TS2 to describe the pressure- and molecular-weight dependence of the glass transition temperature in bulk polymers, as well as its dependence on film thickness in thin films.

On the origin of oscillatory interactions between surfaces mediated by polyelectrolyte solution.

We use a numerical implementation of polymer classical density functional theory with an incompressibility condition to study the system consisting of nonadsorbing polyelectrolytes confined by two planar surfaces and quantify the effective interaction between the two planar surfaces as a function of the polyelectrolyte and salt concentrations. Our results indicate that for the uncharged surfaces (and weakly charged surfaces), the effective interaction primarily consists of a short-range attraction due to the depletion followed by repulsion due to the electric double layer overlapping and electrostatic correlations. For salt-free and low salt concentration systems, the magnitude of the repulsion barrier is determined by the overlap between the electric double layers, while at relatively high salt concentrations, the magnitude of the repulsion barrier is determined by the electrostatic correlations. Due to the competition between the electric double layer and the electrostatic correlations, the magnitude of the repulsion barrier varies nonmonotonically. In contrast, a mean-field Poisson-Boltzmann treatment of the electrostatics predicts a monotonically decreasing repulsion barrier with increasing salt concentration. At moderate salt concentrations, our theory predicts oscillatory interaction profiles. A comparison with the mean-field Poisson-Boltzmann treatment of electrostatics suggests that the oscillations are due primarily to electrostatic correlations.

Rheology of Cellulose Ether Excipients Designed for Hot Melt Extrusion.

A new family of cellulosic ether polymeric excipients has been recently engineered for fabrication of amorphous solid dispersions of active pharmaceutical ingredients via hot-melt extrusion (HME). These hydroxypropyl methyl cellulose excipients enable plasticizer-free melt processing at much lower temperatures (135-160 °C) due to their substantially reduced glass transition temperatures ( Tg = 98-110 °C). The novel amorphous cellulose ethers were found to be rheologically solidlike well above their glass transition ( Tg + 70 °C). We demonstrate that in the pharmaceutically relevant HME processing temperature range these polymers behave similarly to yield-stress fluids and flow only when the applied stress exceeds a critical stress value. This critical stress value (0.50 ± 0.05 MPa, 150 °C) is surprisingly high but is easily achieved under typical HME conditions. The origin of their yield-stress fluidlike behavior is hypothesized to arise from hydrogen bonds of the HPMC materials that act as physical cross-links and do not relax within the measured temperature and time window unless the applied stress exceeds the critical stress value. Support for this hypothesis arises from infrared spectroscopic estimates of the free and bound hydrogen bond levels at end-use temperatures.

Oscillatory and Steady Shear Rheology of Model Hydrophobically Modified Ethoxylated Urethane-Thickened Waterborne Paints.

Hydrophobically modified ethoxylated urethane (HEUR) thickeners are widely used as rheology modifiers for waterborne paints. Although the rheology of HEUR solutions in water is fairly well-understood, their impact on the rheology of waterborne latex/pigment suspensions (formulated paints) is more complicated. We study the shear rheology of model HEUR/latex/TiO2 suspensions in water and investigate the dependence of both oscillatory and steady shear behaviors on the strength of the HEUR hydrophobes. We observe that in both oscillatory and steady shear experiments, rheological curves could be shifted onto a single master curve, demonstrating a "time-hydrophobe superposition". We also note that the oscillatory shear behavior exhibits a power-law spectrum of relaxation times, unlike the single-Maxwellian behavior of pure HEUR solutions. On the basis of these results and earlier experimental and theoretical findings, we propose that the rheology of the HEUR-thickened latex/TiO2 suspensions is mainly determined by the transient network of HEUR-bridged latex particles, with a broad distribution of the characteristic lifetimes of the bridge. The model is found to be in good qualitative and semiquantitative agreement with the experiments for both steady shear and oscillatory shear.

Density functional theory for charged fluids.

An improved density functional theory (DFT) for an inhomogeneous charged system (including electrolyte and/or polyelectrolyte) is proposed based on fundamental measure theory, thermodynamic perturbation theory and mean-spherical approximation. Our DFT combines the existing treatment of hard-sphere contributions using fundamental measure theory (FMT) with a new treatment of the electrostatic correlations for the non-bonded ions and chain connectivity that are approximated by employing a first-order Taylor expansion, with the reference fluid density determined using the technique from Gillespie et al. [D. Gillespie et al., J. Phys.: Condens. Matter, 2002, 14, 12129]. We show that the first-order Taylor expansion for the non-bonded electrostatic correlations yields numerically comparable results to the more involved second-order expansion. Furthermore, we find that the existing treatment of the chain connectivity correlation predicts a spurious layer-by-layer phase at moderately large Bjerrum lengths, which is avoided in our new treatment. These simplifications and improvements should significantly facilitate the implementation and reduce the computational cost.

Formulation-Controlled Positive and Negative First Normal Stress Differences in Waterborne Hydrophobically Modified Ethylene Oxide Urethane (HEUR)-Latex Suspensions.

Hydrophobically modified ethylene oxide urethane (HEUR) associative thickeners are widely used to modify the rheology of waterborne paints. Understanding the normal stress behavior of the HEUR-based paints under high shear is critical for many applications such as brush drag and spreading. We observed that the first normal stress difference, N1, at high shear (large Weissenberg number) can be positive or negative depending on the HEUR hydrophobe strength and concentration. We propose that the algebraic sign of the N1 is primarily controlled by two factors: (a) adsorption of HEURs on the latex surface and (b) the ability of HEURs to form transient molecular bridges between latex particles. Such transient bridges are favored for dispersions with small interparticle distances and dense surface coverages; in these systems; HEUR-bridged latex microstructures flow-align in high shear and exhibit positive N1. In the absence of transient bridges (large interparticle distances, low surface coverage), the dispersion rheology is similar to that of weakly interacting spheres, exhibiting negative N1. The results are summarized in a simplified phase diagram connecting formulation, microstructure, and the N1 behavior.

Modeling the interfacial tension in oil-water-nonionic surfactant mixtures using dissipative particle dynamics and self-consistent field theory.

We use two mesoscale simulation methods, dissipative particle dynamics (DPD) and self-consistent field theory (SCFT), to model interfacial tension in ternary oil/water/surfactant mixtures. For model systems where the oil is dodecane and the surfactant is a linear alkyl ethoxylate, the two methods show a good semiquantitative agreement among themselves and with experimental data. We further discuss the advantages and limitations of the two methods and their possible use for surfactant screening in specific industrial applications.

Modeling polymer-induced interactions between two grafted surfaces: comparison between interfacial statistical associating fluid theory and self-consistent field theory.

The interaction between two polymer grafted surfaces is important in many applications, such as nanocomposites, colloid stabilization, and polymer alloys. In our previous work [Jain et al., J. Chem. Phys. 128, 154910 (2008)], we showed that interfacial statistical associating fluid density theory (iSAFT) successfully calculates the structure of grafted polymer chains in the absence/presence of a free polymer. In the current work, we have applied this density functional theory to calculate the force of interaction between two such grafted monolayers in implicit good solvent conditions. In particular, we have considered the case where the segment sizes of the free (sigma(f)) and grafted (sigma(g)) polymers are different. The interactions between the two monolayers in the absence of the free polymer are always repulsive. However, in the presence of the free polymer, the force either can be purely repulsive or can have an attractive minimum depending upon the relative chain lengths of the free (N(f)) and grafted polymers (N(g)). The attractive minimum is observed only when the ratio alpha = N(f)/N(g) is greater than a critical value. We find that these critical values of alpha satisfy the following scaling relation: rho(g) square root(N(g)) beta(3) proportional to alpha(-lambda), where beta = sigma(f)/sigma(g) and lambda is the scaling exponent. For beta = 1 or the same segment sizes of the free and grafted polymers, this scaling relation is in agreement with those from previous theoretical studies using self-consistent field theory (SCFT). Detailed comparisons between iSAFT and SCFT are made for the structures of the monolayers and their forces of interaction. These comparisons lead to interesting implications for the modeling of nanocomposite thermodynamics.

Modeling the thermodynamics of the interaction of nanoparticles with cell membranes.

Interactions between nanoparticles and cell membranes may play a crucial role in determining the cytotoxicity of nanoparticles as well as their potential application as drug delivery vehicles or therapeutic agents. It has been shown that such interactions are often determined not by biochemical but by physicochemical factors (e.g., nanoparticle size, hydrophobicity, and surface charge density). Here, we propose a mesoscale thermodynamic model describing the transitions in membrane morphology observed after exposure to various types of nanoparticles. Our simulations demonstrate under which conditions (determined by particle size and hydrophilic/hydrophobic interactions) the particles can adsorb into the membrane or compromise the membrane integrity to result in the formation of nanosized holes. The model could be refined to include a more accurate description of various phospholipid membranes, and its results could be applied in the design of specific nanoparticles for various biomedical applications.

Effect of hydrodynamic interactions on the evolution of chemically reactive ternary mixtures.

We investigate the structural evolution of an A/B/C ternary mixture in which the A and B components can undergo a reversible chemical reaction to form C. We developed a lattice Boltzmann model for this ternary mixture that allows us to capture both the reaction kinetics and the hydrodynamic interactions within the system. We use this model to study a specific reactive mixture in which C acts as a surfactant, i.e., the formation of C at the A/B interface decreases the interfacial tension between the A and B domains. We found that the dynamics of the system is different for fluids in the diffusive and viscous regimes. In the diffusive regime, the formation of a layer of C at the interface leads to a freezing of the structural evolution in the fluid; the values of the reaction rate constants determine the characteristic domain size in the system. In the viscous regime, where hydrodynamic interactions are important, interfacial reactions cause a slowing down of the domain growth, but do not arrest the evolution of the mixture. The results provide guidelines for controlling the morphology of this complex ternary fluid.