Mesoscale simulation of polymer reaction equilibrium: Combining dissipative particle dynamics with reaction ensemble Monte Carlo. II. Supramolecular diblock copolymers

Determination of interaction parameters in a bottom-up approach employed in reactive dissipative particle dynamics simulations for thermosetting polymers

The limitations in previous dissipative particle dynamics (DPD) studies confined simulations to a narrow resin range. This study refines DPD parameter calculation methodology, extending its application to diverse polymer materials. Using a bottom-up approach with molecular dynamics (MD) simulations, we evaluated solubility parameters and bead number density governing nonbonded interactions via the Flory-Huggins parameter and covalent-bonded interactions. Two solubility parameter methods, Hildebrand and Krevelen-Hoftyzer, were compared for DPD simulations. The Hildebrand method, utilizing MD simulations, demonstrates higher consistency and broader applicability in determining solubility parameters for all DPD particles. The DPD/MD curing reaction process was examined in three epoxy systems: DGEBA/4,4'-DDS, DGEBA/MPDA and DGEBA/DETA. Calculations for the curing profile, gelation point, radial distribution function and branch ratio were performed. Compared to MD data for DGEBA/4,4'-DDS, the maximum deviation in secondary reactions between epoxy and amine groups according to DPD simulations with Krevelen-Hoftyzer was 14.8%, while with the Hildebrand method, it was 1.7%. The accuracy of the DPD curing reaction in reproducing the structural properties verifies its expanded application to general polymeric material simulations. The proposed curing DPD simulations, with a short run time and minimal computational resources, contributes to high-throughput screening for optimal resins and investigates mesoscopic inhomogeneous structures in large resin systems.

How fork-length asymmetry affects solvent connectivity and diffusion in grafted polymeric model membranes

The hydrophilic pore morphology and solvent diffusion within model (amphiphilic) polymer membranes are simulated by dissipative particle dynamics (DPD). The polymers are composed of a backbone of 18 covalently bonded A beads to which at regular intervals side chains are attached. The side chains are composed of linear Ap chains (i.e., -A1-A2…Ap) from which two branches, [AsC] and [ArC], split off (s ≤ r). C beads serve as functionalized hydrophilic pendent sites. The branch lengths (s + 1 and r + 1) are varied. Five repeat unit designs (with general formula A3[Ap[AsC][ArC]]) are considered: A2[A3C][A3C] (symmetric branching), A2[A2C][A4C], A2[AC][A5C], A2[C][A6C] (highly asymmetric branching), and A4[AC][A3C]. The distribution of water (W) and W diffusion through nanophase segregated hydrophilic pores is studied. For similar primary length p, an increase in side chain symmetry favors hydrophilic pore connectivity and long-range water transport. C beads located on the longer [ArC] branches reveal the highest C bead mobility and are more strongly associated with water than the C beads on the shorter [AsC] branches. The connectivity of hydrophilic (W and W + C) phases through mapped replica of selected snapshots obtained from Monte Carlo tracer diffusion simulations is in line with trends found from the W bead diffusivities during DPD simulations. The diffusive pathways for protons (H+) in proton exchange membranes and for hydronium (OH-) in anion exchange membranes are the same as for solvents. Therefore, control of the side chain architecture is an interesting design parameter for optimizing membrane conductivities.

Molecular Weight Distribution of Branched Polymers: Comparison between Monte Carlo Simulation and Flory-Stockmayer Theory

It is challenging to predict the molecular weight distribution (MWD) for a polymer with a branched architecture, though such information will significantly benefit the design and development of branched polymers with desired properties and functions. A Monte Carlo (MC) simulation method based on the Gillespie algorithm is developed to quickly compute the MWD of branched polymers formed through step-growth polymerization, with a branched polyetherimide from two backbone monomers (4,4′-bisphenol A dianhydride and m-phenylenediamine), a chain terminator (phthalic anhydride), and a branching agent (tris[4-(4-aminophenoxy)phenyl] ethane) as an example. This polymerization involves four reactions that can be all reduced to a condensation reaction between an amine group and a carboxylic anhydride group. A comparison between the MC simulation results and the predictions of the Flory-Stockmayer theory on MWD shows that the rates of the reactions are determined by the concentrations of the functional groups on the monomers involved in each reaction. It further shows that the Flory-Stockmayer theory predicts MWD well for systems below the gel point but starts to fail for systems around or above the gel point. However, for all the systems, the MC method can be used to reliably predict MWD no matter if they are below or above the gel point. Even for a macroscopic system, a converging distribution can be quickly obtained through MC simulations on a system of only a few hundred to a few thousand monomers that have the same molar ratios as in the macroscopic system.

Comprehensive review of the interfacial behavior of water/oil/surfactant systems using dissipative particle dynamics simulation

A comprehensive understanding of interfacial behavior in water/oil/surfactant systems is critical to evaluating the performance of emulsions in various industries, specifically in the oil and gas industry. To gain fundamental knowledge regarding this interfacial behavior, atomistic methods, e.g., molecular dynamics (MD) simulation, can be employed; however, MD simulation cannot handle phenomena that require more than a million atoms. The coarse-grained mesoscale methods were introduced to resolve this issue. One of the most effective mesoscale coarse-grained approaches for simulating colloidal systems is dissipative particle dynamics (DPD), which bridges the gap between macroscopic time and length scales and molecular-scale simulation. This work reviews the fundamentals of DPD simulation and its progress on colloids and interface systems, especially surfactant/water/oil mixtures. The effects of temperature, salt content, a water/oil ratio, a shear rate, and a type of surfactant on the interfacial behavior in water/oil/surfactant systems using DPD simulation are evaluated. In addition, the obtained results are also investigated through the lens of the chemistry of surfactants and emulsions. The outcome of this comprehensive review demonstrates the importance of DPD simulation in various processes with a focus on the colloidal and interfacial behavior of surfactants at water-oil interfaces.

Interactions of Crosslinked Polyacrylic Acid Polyelectrolyte Gels with Nonionic and Ionic Surfactants

The morphology and stability of surfactant-loaded polyelectrolyte gels are of great interest for a variety of personal care, cosmetic, and pharmaceutical products. However, the mechanisms of surfactant interactions with gel-forming polymers are poorly understood and experimentally challenging. The aim of this work is to explore in silico the specifics of surfactant absorption within polyelectrolyte gels drawing on the examples of typical non-ionic octaethylene glycol monooctyl ether (C₈E₈) and anionic sodium dodecyl sulfate (SDS) surfactants and polyacrylic acid modified with hydrophobic sidechains mimicking the practically important Carbopol polymer. Using the systematically parameterized coarse-grained dissipative particle dynamics models, we generate and characterize the equilibrium conformations and swelling of the polymer films in aqueous solutions with the surfactant concentrations varied up to the critical micelle concentration (cmc). We discover the striking difference in interactions of Carbopol-like polymers with nonionic and ionic surfactants under mildly acidic conditions. The sorption of C₈E₈ within the polymer film is found substantial. As the surfactant concentration increases, the polymer film swells and, close to cmc, becomes unstable due to the formation and growth of water pockets filled with surfactant micelles. Sorption of SDS at the same bulk concentrations is found much lower, with only about 1% of surfactant mass fraction achieved at cmc. As the SDS concentration increases further, a lamellae structure is formed within the film, which remains stable. Reduced swelling and higher stability indicate better prospects of using SDS-type surfactants with Carbopol-based gels in formulations for detergents and personal care products.

Dissipative particle dynamics simulations in colloid and Interface science: a review

Dissipative particle dynamics (DPD) is one of the most efficient mesoscale coarse-grained methodologies for modeling soft matter systems. Here, we comprehensively review the progress in theoretical formulations, parametrization strategies, and applications of DPD over the last two decades. DPD bridges the gap between the microscopic atomistic and macroscopic continuum length and time scales. Numerous efforts have been performed to improve the computational efficiency and to develop advanced versions and modifications of the original DPD framework. The progress in the parametrization techniques that can reproduce the engineering properties of experimental systems attracted a lot of interest from the industrial community longing to use DPD to characterize, help design and optimize the practical products. While there are still areas for improvements, DPD has been efficiently applied to numerous colloidal and interfacial phenomena involving phase separations, self-assembly, and transport in polymeric, surfactant, nanoparticle, and biomolecules systems.

Structural properties of cationic surfactant-fatty alcohol bilayers: insights from dissipative particle dynamics

Bilayers, self-assembled by cationic surfactants and fatty alcohols in water, are the basic units of lamellar gel networks - creamy formulations extensively used in cosmetics and pharmaceutics. Mesoscopic modelling and study of the bilayers formed by single- or double-tail cationic surfactants (CTAC or DHDAC), and fatty alcohols (FAs) in the lamellar fluid and gel phases were employed. Fatty alcohols with alkyl tail equal to or greater than the surfactant alkyl tail, i.e., C₁₆FA or C₁₈FA and C₂₂FA, were considered. A model formulation was explored with the FA concentration greater than that of the surfactant and the structure of the fluid and gel bilayers in tensionless state characterised via the density profiles across the bilayers, orientational order parameters of the surfactant and FA chains, intrinsic analysis of the bilayer interfaces, and bending rigidity. The intrinsic analysis allows identification and quantification of the coexistence of the interdigitated and non-interdigitated phases present within the gel bilayers. The FA chains were found to conform the primary scaffolding of the bilayers while the surfactant chains tessellate bilayer monolayers from their water-hydrophobic interface. Further, the overlap of the FA chains from the apposed monolayers of the fluid bilayers rises with increasing FA length. Finally, the prevalence of the non-interdigitated phase over the interdigitated phase within the gel bilayers becomes enhanced upon the FA length increase with a preference of the surfactant chains to reside in the non-interdigitated phase rather than the interdigitated phase.

Thermoset resin curing simulation using quantum-chemical reaction path calculation and dissipative particle dynamics

Thermoset resin, which is commonly used as a matrix in carbon-fiber-reinforced plastic, requires curing procedures. We propose a curing simulation technique involving a dissipative particle dynamics (DPD) simulation, which can simulate a larger system and longer time period than those of conventional all-atom molecular dynamics (AA-MD) simulations. The proposed curing DPD simulation can represent the thermoset resin exothermic reaction process precisely by considering each reactivity according to the reaction types calculated via quantum-chemical reaction path calculations. The cure reaction process given by the curing DPD simulation agrees well with that given by a conventional curing AA-MD simulation, but with run-time and computational-resource reductions of 1/480 and 1/10 times, respectively. We also conduct reverse mapping, through which the AA-MD system can be reconstructed from the DPD system, to evaluate the structural and thermomechanical properties. The X-ray diffraction pattern and thermomechanical properties of the reconstructed system agree well with those of the systems derived from the curing AA-MD simulation and experimental setup. Therefore, a cured-resin AA-MD system can be obtained from a curing DPD simulation at an extremely low computational cost, and the thermomechanical properties can be evaluated precisely using this system. The proposed curing simulation technique can be applied in high-throughput screening for better materials properties and in large system calculations.

Simulation Study of Process-Controlled Supramolecular Block Copolymer Phase Separation with Reversible Reaction Algorithm

A supramolecular diblock copolymer formed by reversible bonds between the two blocks shows a rich microphase separation behavior and has great application potential in stimuli-responsive materials. We propose a novel method to describe supramolecular reactions in dissipative particle dynamics, which includes a reversible reaction to accurately reproduce the strength, saturation, and dynamic properties of the reversible bonds in the simulations. The thermodynamic properties and dynamic processes of the supramolecular diblock copolymer melts in both equilibrium and non-equilibrium states were studied using this method. The simulation results show that the method can faithfully characterize phase behaviors and dynamic properties of supramolecular diblock copolymer melts, especially in a non-equilibrium state, which provides a novel tool to unveil self-assembly mechanism and describe the properties of supramolecular block copolymers.

Bilayer sheet protrusions and budding from bilayer membranes induced by hydrolysis and condensation reactions

Shape transformations of flat bilayer membranes and vesicles induced by hydrolysis and condensation reactions of amphiphilic molecules are studied using coarse-grained molecular dynamics simulations. The hydrolysis and condensation reactions result in the formation and dissociation of amphiphilic molecules, respectively. Asymmetric reactions between the inner and outer leaflets of a vesicle can transport amphiphilic molecules between the leaflets. It is found that the resulting area difference between the two leaflets induces bilayer sheet protrusion (BP) and budding at low reduced volumes of the vesicles, whereas BP only occurs at high reduced volumes. The probabilities of these two types of transformations depend on the shear viscosity of the surrounding fluids compared to the membrane as well as the reaction rates. A higher surrounding fluid viscosity leads to more BP formation. The inhomogeneous spatial distribution of the hydrophobic reaction products forms the nuclei of BP formation, and faster diffusion of the products enhances BP formation. Our results suggest that adjustment of the viscosity is important to control membrane shape transformations in experiments.

One-pot production of porous assemblies by PISA of star architecture copolymers: a simulation study

Porous vesicles can be produced in one-pot by the PISA of star architecture copolymers.

Computer Simulations of Static and Dynamical Properties of Weak Polyelectrolyte Nanogels in Salty Solutions

We investigate the chemical equilibria of weak polyelectrolyte nanogels with reaction ensemble Monte Carlo simulations. With this method, the chemical identity of the nanogel monomers can change between neutral or charged following the acid-base equilibrium reaction HA ⇌ A- + H⁺. We investigate the effect of changing the chemical equilibria by modifying the dissociation constant K a . These simulations allow for the extraction of static properties like swelling equilibria and the way in which charge-both monomer and ionic-is distributed inside the nanogel. Our findings reveal that, depending on the value of K a , added salt can either increase or decrease the gel size. Using the calculated mean-charge configurations of the nanogel from the reaction ensemble simulation as a quenched input to coupled lattice-Boltzmann molecular dynamics simulations, we investigate dynamical nanogel properties such as the electrophoretic mobility μ and the diffusion coefficient D.

Wang-Landau Reaction Ensemble Method: Simulation of Weak Polyelectrolytes and General Acid-Base Reactions

We present a novel method for the study of weak polyelectrolytes and general acid-base reactions in molecular dynamics and Monte Carlo simulations. The approach combines the advantages of the reaction ensemble and the Wang-Landau sampling method. Deprotonation and protonation reactions are simulated explicitly with the help of the reaction ensemble method, while the accurate sampling of the corresponding phase space is achieved by the Wang-Landau approach. The combination of both techniques provides a sufficient statistical accuracy such that meaningful estimates for the density of states and the partition sum can be obtained. With regard to these estimates, several thermodynamic observables like the heat capacity or reaction free energies can be calculated. We demonstrate that the computation times for the calculation of titration curves with a high statistical accuracy can be significantly decreased when compared to the original reaction ensemble method. The applicability of our approach is validated by the study of weak polyelectrolytes and their thermodynamic properties.

Self-consistent description of electrokinetic phenomena in particle-based simulations

A new computational method is presented for study suspensions of charged particles undergoing fluctuating hydrodynamic and electrostatic interactions. The proposed model is appropriate for polymers, proteins, and porous particles embedded in a continuum electrolyte. A self-consistent Langevin description of the particles is adopted in which hydrodynamic and electrostatic interactions are included through a Green's function formalism. An Ewald-like split is adopted in order to satisfy arbitrary boundary conditions for the Stokeslet and Poisson Green functions, thereby providing a formalism that is applicable to any geometry and that can be extended to deformable objects. The convection-diffusion equation for the continuum ions is solved simultaneously considering Nernst-Planck diffusion. The method can be applied to systems at equilibrium and far from equilibrium. Its applicability is demonstrated in the context of electrokinetic motion, where it is shown that the ionic clouds associated with individual particles can be severely altered by the flow and concentration, leading to intriguing cooperative effects.

Reaction Ensemble Monte Carlo Simulation of Complex Molecular Systems

Acceptance rules for reaction ensemble Monte Carlo (RxMC) simulations containing classically modeled atomistic degrees of freedom are derived for complex molecular systems where insertions and deletions are achieved gradually by utilizing the continuous fractional component (CFC) method. A self-consistent manner in which to utilize statistical mechanical data contained in ideal gas free energy parameters during RxMC moves is presented. The method is tested by applying it to two previously studied systems containing intramolecular degrees of freedom: the propene metathesis reaction and methyl-tert-butyl-ether (MTBE) synthesis. Quantitative agreement is found between the current results and those of Keil et al. (J. Chem. Phys. 2005, 122, 164705) for the propene metathesis reaction. Differences are observed between the equilibrium concentrations of the present study and those of Lísal et al. (AIChE J. 2000, 46, 866-875) for the MTBE reaction. It is shown that most of this difference can be attributed to an incorrect formulation of the Monte Carlo acceptance rule. Efficiency gains using CFC MC as opposed to single stage molecule insertions are presented.

Directed assembly of supramolecular copolymers in thin films: thermodynamic and kinetic advantages

Using computer simulation of a coarse-grained model for supramolecular polymers, we investigate the potential of quasiblock copolymers (QBCPs) assembled on chemically patterned substrates for creating device-oriented nanostructures. QBCPs are comprised of AB diblock copolymers and supramolecular B segments that can reversibly bond to any available B terminus, on either the copolymers or the B oligomers, creating a polydisperse blend of B homopolymers, and AB and ABA copolymers. We demonstrate the defect-free replication of patterns with perpendicularly crossing, A-preferential lines, where the same QBCP can simultaneously replicate patterns differing by up to 50% in their length scales. We demonstrate how the pattern affects the distribution of molecular architectures and the key role of supramolecular associations for replicating patterns with different length scales.