Molecular View of Polymer/Water Interfaces in Latex Paint

Functional Lignin Building Blocks: Reactive Vinyl Esters with Acrylic Acid

Introducing vinyl groups onto the backbone of technical lignin provides an opportunity to create highly reactive renewable polymers suitable for radical polymerization. In this work, the chemical modification of softwood kraft lignin was pursued with etherification, followed by direct esterification with acrylic acid (AA). In the first step, phenolic hydroxyl and carboxylic acid groups were derivatized into aliphatic hydroxyl groups using ethylene carbonate and an alkaline catalyst. The lignin was subsequently fractionated using a downward precipitation method to recover lignin of defined molar mass and solubility. After recovery, the resulting material was then esterified with AA, resulting in lignin with vinyl functional groups. The first step resulted in approximately 90% of phenolic hydroxyl groups being converted into aliphatic hydroxyls, while the downward fractionation resulted in three samples of lignin with defined molar masses. For the esterification reaction, the weight ratio of reagents, reaction temperature, and reaction time were evaluated as factors that would influence the modification efficacy. ¹³C NMR spectroscopy analysis of lignin samples before and after esterification showed that the optimized reaction conditions could reach approximately 40% substitution of aliphatic hydroxyl groups. Both steps only used lignin and the modifying reagent (no solvent), with the possibility of recovery and reuse of the reagent by dilution and distillation. An additional second esterification step of the resulting lignin sample with acetic acid or propionic acid converted 90% of remaining hydroxyl groups into short-chain carbon aliphatic esters, making a hydrophobic material suitable for further copolymerization with synthetic hydrophobic monomers.

Adsorption Energetics of Amino Acid Analogs on Polymer/Water Interfaces Studied by All-Atom Molecular Dynamics Simulation and a Theory of Solutions

Energetics of adsorption was addressed with all-atom molecular dynamics simulation on the interfaces of poly(2-methoxyethyl acrylate) (PMEA), poly(methyl methacrylate) (PMMA), and poly(butyl acrylate) (PBA) with water. A wide variety of adsorbate solutes were examined, and the free energy of adsorption was computed with the method of energy representation. It was found that the adsorption free energy was favorable (negative) for all the combinations of solute and polymer, and among PMEA, PMMA, and PBA, the strongest adsorption was observed on PMMA for the hydrophobic solutes and on PMEA for the hydrophilic ones. According to the decomposition of the adsorption free energy into the contributions from polymer and water, it was seen that the polymer contribution is larger in magnitude with the solute size. The total free energy of adsorption was correlated well with the solvation free energy in bulk water only for hydrophobic solutes. The roles of the intermolecular interaction components such as electrostatic, van der Waals, and excluded-volume were further studied, and the electrostatic component was influential only in determining the polymer dependences of the adsorption propensities of hydrophilic solutes. The extent of adsorption was shown to be ranked by the van der Waals component in the solute-polymer interaction separately over the hydrophilic and hydrophobic solutes, with the excluded-volume effect from water pointed out to also drive the adsorption.

Molecular dynamics study of wetting of alkanes on water: from high temperature to the supercooled region and the influence of second inflection points of interfacial tensions

To explore the wetting behavior of alkanes on bulk water interfaces, molecular dynamics (MD) simulations were carried out for united-atom PYS alkane models, and for SPC/E and TIP4P/2005 water models over a wide temperature range. The MD results at each temperature were used to find (1) the surface tension of the alkanes (octane, nonane) and water, and (2) the interfacial tensions of the alkane-water systems. These quantities were then used to calculate the spreading coefficient (S) and contact angle (θc) for each alkane on water. At higher temperatures, the contact angle of octane and nonane on water is found to behave in accord with conventional expectations, i.e., it decreases with increasing temperature for both water models as each system approaches the usual high-temperature transition to perfect wetting. At lower temperatures, we found an unusual temperature dependence of S and θc for each PYS alkane on SPC/E water. In contrast to conventional expectations, θc decreases with a decrease in the temperature. For octane-SPC/E water, this unusual behavior of θc occurs due to the presence of second inflection points (SIP) in the vapor-water and the octane-water interfacial tensions, whereas the SIP effect is much less important for the nonane-water system. The unusual temperature dependence of θc observed for nonane on SPC/E water is also found for nonane on TIP4P/2005 water. On the other hand, such unusual wetting behavior has not been observed in the PYS octane-TIP4P/2005 water system, except possibly for the two lowest temperatures studied.

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/TiO₂ 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/TiO₂ 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.

Hydrogen bond directed surface dynamics at tactic poly(methyl methacrylate)/water interface

The complexity of induced ordering for tactic poly(methyl methacrylate) (PMMA) thin films in contact with water is examined through all-atom molecular dynamics with validated potentials. We observe that for the water molecules that are hydrogen bonded to the PMMA surface, the isotactic and atactic PMMA show a 33% longer relaxation time compared to syndiotactic PMMA. Almost 94% of hydrogen bonds are with the carbonyl groups of PMMA, irrespective of temperature and tacticity. The stability in re-orientation and nature of hydrogen bond participation for the carbonyl groups as well as about 20% higher interaction energies of carbonyl group hydrogen bonded with water for atactic form indicates existence of cooperative effects. Quantifying the dynamics of hydrogen bond at the tactic interface is important in understanding the role tacticity plays in controlling adhesion and biocompatibility, a design choice that has been gaining ground in the soft material science community.

Surfactant Binding to Polymer-Water Interfaces in Atomistic Simulations

Attractive interactions between additive molecules and particle surfaces are key parameters in the design of waterborne suspensions and coatings. We use atomistic molecular dynamics (MD) simulations to determine the potential of mean force for a commonly used industrial surfactant sodium dodecyl sulfate (SDS) interacting with acrylate latex particles. We investigate how the potential of mean force and binding free energy depend on the amount of SDS adsorbed, solution ionic strength, and presence of other charged groups on the particle surface. We show that the potential of mean force for SDS is a sum of two independent terms, from the hydrophobic surfactant tail and charged headgroup: dragging the surfactant tail into solution contributes a linear potential of about kT per CH2 group, while the headgroup is repelled by like charges on the surface with a potential of about the zeta potential. Commercial acrylate latex particles also bear multivalent charged "hairs" as a remnant of their synthesis. These charged hairs result in a heterogeneously charged surface, for which SDS binds more or less strongly depending on the local environment.

Atomistic Molecular Dynamics Simulations of Charged Latex Particle Surfaces in Aqueous Solution

Charged particles in aqueous suspension form an electrical double layer at their surfaces, which plays a key role in suspension properties. For example, binder particles in latex paint remain suspended in the can because of repulsive forces between overlapping double layers. Existing models of the double layer assume sharp interfaces bearing fixed uniform charge, and so cannot describe aqueous binder particle surfaces, which are soft and diffuse, and bear mobile charge from ionic surfactants as well as grafted multivalent oligomers. To treat this industrially important system, we use atomistic molecular dynamics simulations to investigate a structurally realistic model of commercial binder particle surfaces, informed by extensive characterization of particle synthesis and surface properties. We determine the interfacial profiles of polymer, water, bound and free ions, from which the charge density and electrostatic potential can be calculated. We extend the traditional definitions of the inner and outer Helmholtz planes to our diffuse interfaces. Beyond the Stern layer, the simulated electrostatic potential is well described by the Poisson-Boltzmann equation. The potential at the outer Helmholtz plane compares well to the experimental zeta potential. We compare particle surfaces bearing two types of charge groups, ionic surfactant and multivalent oligomers, with and without added salt. Although the bare charge density of a surface bearing multivalent oligomers is much higher than that of a surfactant-bearing surface at realistic coverage, greater counterion condensation leads to similar zeta potentials for the two systems.

Multiscale Molecular Dynamics Simulations of Model Hydrophobically Modified Ethylene Oxide Urethane Micelles

The flower-like micelles of various aggregation numbers of a model hydrophobically modified ethylene oxide urethane (HEUR) molecule, C16E45C16, and their corresponding starlike micelles, containing the surfactants C16E22 and C16E23, were studied by atomistic and coarse-grained molecular dynamic (MD) simulations. We used free energies from umbrella sampling to calculate the size distribution of micelle sizes and the average time for escape of a hydrophobic group from the micelle. Using the coarse-grained MARTINI force field, the most probable size of the model HEUR molecule was thereby determined to be about 80 hydrophobes per micelle and the average hydrophobe escape time to be about 0.1 s, both of which are consistent with previous experimental studies. Atomistic simulations reveal that hydrogen bond formation and the mean lifetime of hydration waters of the poly(ethylene oxide) (or PEO) groups are location-dependent in the HEUR micelle, with PEO groups immediately adjacent to the C16 groups forming the fewest hydrogen bonds with water and having hydration waters with longer lifetimes than those of the PEO groups located further away from the C16 groups.

Potentials of mean force and escape times of surfactants from micelles and hydrophobic surfaces using molecular dynamics simulations

We calculate potentials of mean force (PMFs) and mean first passage times for a surfactant to escape a micelle, for both ionic sodium dodecyl sulfate (SDS) and nonionic ethoxylated alcohol (C12E5) micelles using both atomistic and coarse-grained molecular dynamics (MD) simulations. The PMFs are obtained by umbrella sampling and used in a Smoluchowski first-passage-time theory to obtain the times for a surfactant to escape a micelle. The calculated mean first passage time for an SDS molecule to break away from a micelle (with an aggregation number of 60) is around 2 μs, which is consistent with previous experimental measurements of the "fast relaxation time" for exchange of surfactants between the micellar phase and the bulk solvent. The corresponding escape time calculated for a nonionic ethoxylated alcohol C12E5, with the same tail length as SDS, is 60 μs, which is significantly longer than for SDS primarily because the PMF for surfactant desorption is about 3kT smaller than for C12E5. We also show that two coarse-grained (CG) force fields, MARTINI and SDK, give predictions similar to the atomistic CHARMM force field for the nonionic C12E5 surfactant, but for the ionic SDS surfactant, the CG simulations give a PMF similar to that obtained with CHARMM only if long-range electrostatic interactions are included in the CG simulations, rather than using a shifted truncated electrostatic interaction. We also calculate that the mean first passage time for an SDS and a C12E5 to escape from a latex binder surface is of the order of milliseconds, which is more than 100 times longer than the time for escape from the micelle, indicating that in latex waterborne coatings, SDS and C12E5 surfactants likely bind preferentially to the latex polymer interface rather than form micelles, at least at low surfactant concentrations.