Investigating microstructure evolution in block copolymer membranes.

Block copolymer self-assembly in conjunction with nonsolvent-induced phase separation (SNIPS) has been increasingly leveraged to fabricate integral-asymmetric membranes. The large number of formulation and processing parameters associated with SNIPS, however, has prevented the reliable construction of high performance membranes. In this study, we apply dynamical self-consistent field theory to model the SNIPS process and investigate the effect of various parameters on the membrane morphology: solvent selectivity, nonsolvent selectivity, initial film composition, and glass transition composition. We examine how solvent selectivity and concentration of polymers in the film impact the structure of micelles that connect to form the membrane matrix. In particular, we find that preserving the order in the surface layer and forming a connection between the supporting and surface layer are nontrivial and sensitive to each parameter studied. The effect of each parameter is discussed, and suggestions are made for successfully fabricating viable block copolymer membranes.

Coarsening dynamics of ternary polymer solutions with mobility and viscosity contrasts.

Using phase-field simulations, we investigate the bulk coarsening dynamics of ternary polymer solutions undergoing a glass transition for two models of phase separation: diffusion only and with hydrodynamics. The glass transition is incorporated in both models by imposing mobility and viscosity contrasts between the polymer-rich and polymer-poor phases of the evolving microstructure. For microstructures composed of polymer-poor clusters in a polymer-rich matrix, the mobility and viscosity contrasts significantly hinder coarsening, effectively leading to structural arrest. For microstructures composed of polymer-rich clusters in a polymer-poor matrix, the mobility and viscosity contrasts do not impede domain growth; rather, they change the transient concentration of the polymer-rich phase, altering the shape of the discrete domains. This effect introduces several complexities to the coarsening process, including percolation inversion of the polymer-rich and polymer-poor phases-a phenomenon normally attributed to viscoelastic phase separation.

Molecularly informed field theory for estimating critical micelle concentrations of intrinsically disordered protein surfactants.

The critical micelle concentration (CMC) is a crucial parameter in understanding the self-assembly behavior of surfactants. In this study, we combine simulation and experiment to demonstrate the predictive capability of molecularly informed field theories in estimating the CMC of biologically based protein surfactants. Our simulation approach combines the relative entropy coarse-graining of small-scale atomistic simulations with large-scale field-theoretic simulations, allowing us to efficiently compute the free energy of micelle formation necessary for the CMC calculation while preserving chemistry-specific information about the underlying surfactant building blocks. We apply this methodology to a unique intrinsically disordered protein platform capable of a wide variety of tailored sequences that enable tunable micelle self-assembly. The computational predictions of the CMC closely match experimental measurements, demonstrating the potential of molecularly informed field theories as a valuable tool to investigate self-assembly in bio-based macromolecules systematically.

Emergence of a Spin Microemulsion in Spin-Orbit Coupled Bose-Einstein Condensates.

We report the first numerical prediction of a "spin microemulsion"-a phase with undulating spin domains resembling classical bicontinuous oil-water-surfactant emulsions-in two-dimensional systems of spinor Bose-Einstein condensates with isotropic Rashba spin-orbit coupling. Using field-theoretic numerical simulations, we investigated the melting of a low-temperature stripe phase with supersolid character and find that the stripes lose their superfluidity at elevated temperature and undergo a Kosterlitz-Thouless-like transition into a spin microemulsion. Momentum distribution calculations highlight a thermally broadened occupation of the Rashba circle of low-energy states with macroscopic and isotropic occupation around the ring. We provide a finite-temperature phase diagram that positions the emulsion as an intermediate, structured isotropic phase with residual quantum character before transitioning at higher temperature into a structureless normal fluid.

Assessment of the partial saddle point approximation in field-theoretic polymer simulations.

Field-theoretic simulations are numerical treatments of polymer field theory models that go beyond the mean-field self-consistent field theory level and have successfully captured a range of mesoscopic phenomena. Inherent in molecularly-based field theories is a "sign problem" associated with complex-valued Hamiltonian functionals. One route to field-theoretic simulations utilizes the complex Langevin (CL) method to importance sample complex-valued field configurations to bypass the sign problem. Although CL is exact in principle, it can be difficult to stabilize in strongly fluctuating systems. An alternate approach for blends or block copolymers with two segment species is to make a "partial saddle point approximation" (PSPA) in which the stiff pressure-like field is constrained to its mean-field value, eliminating the sign problem in the remaining field theory, allowing for traditional (real) sampling methods. The consequences of the PSPA are relatively unknown, and direct comparisons between the two methods are limited. Here, we quantitatively compare thermodynamic observables, order-disorder transitions, and periodic domain sizes predicted by the two approaches for a weakly compressible model of AB diblock copolymers. Using Gaussian fluctuation analysis, we validate our simulation observations, finding that the PSPA incorrectly captures trends in fluctuation corrections to certain thermodynamic observables, microdomain spacing, and location of order-disorder transitions. For incompressible models with contact interactions, we find similar discrepancies between the predictions of CL and PSPA, but these can be minimized by regularization procedures such as Morse calibration. These findings mandate caution in applying the PSPA to broader classes of soft-matter models and systems.

Improved Elastic Recovery from ABC Triblock Terpolymers.

The promise of ABC triblock terpolymers for improving the mechanical properties of thermoplastic elastomers is demonstrated by comparison with symmetric ABA/CBC analogs having similar molecular weights and volume fraction of B and A/C domains. The ABC architecture enhances elasticity (up to 98% recovery over 10 cycles) in part through essentially full chain bridging between discrete hard domains leading to the minimization of mechanically unproductive loops. In addition, the unique phase space of ABC triblocks also enables the fraction of hard-block domains to be higher (fhard ≈ 0.4) while maintaining elasticity, which is traditionally only possible with non-linear architectures or highly asymmetric ABA triblock copolymers. These advantages of ABC triblock terpolymers provide a tunable platform to create materials with practical applications while improving our fundamental understanding of chain conformation and structure-property relationships in block copolymers.

A molecularly informed field-theoretic study of the complexation of polycation PDADMA with mixed micelles of sodium dodecyl sulfate and ethoxylated surfactants.

The self-assembly and phase separation of mixtures of polyelectrolytes and surfactants are important to a range of applications, from formulating personal care products to drug encapsulation. In contrast to systems of oppositely charged polyelectrolytes, in polyelectrolyte-surfactant systems the surfactants micellize into structures that are highly responsive to solution conditions. In this work, we examine how the morphology of micelles and degree of polyelectrolyte adsorption dynamically change upon varying the mixing ratio of charged and neutral surfactants. Specifically, we consider a solution of the cationic polyelectrolyte polydiallyldimethylammonium, anionic surfactant sodium dodecyl sulfate, neutral ethoxylated surfactants (C[Formula: see text]EO[Formula: see text]), sodium chloride salt, and water. To capture the chemical specificity of these species, we leverage recent developments in constructing molecularly informed field theories via coarse-graining from all-atom simulations. Our results show how changing the surfactant mixing ratios and the identity of the nonionic surfactant modulates micelle size and surface charge, and as a result dictates the degree of polyelectrolyte adsorption. These results are in semi-quantitative agreement with experimental observations on the same system.

Design of Soft Material Surfaces with Rationally Tuned Water Diffusivity.

Water structure and dynamics can be key modulators of adsorption, separations, and reactions at soft material interfaces, but systematically tuning water environments in an aqueous, accessible, and functionalizable material platform has been elusive. This work leverages variations in excluded volume to control and measure water diffusivity as a function of position within polymeric micelles using Overhauser dynamic nuclear polarization spectroscopy. Specifically, a versatile materials platform consisting of sequence-defined polypeptoids simultaneously offers a route to controlling the functional group position and a unique opportunity to generate a water diffusivity gradient extending away from the polymer micelle core. These results demonstrate an avenue not only to rationally design the chemical and structural properties of polymer surfaces but also to design and tune the local water dynamics that, in turn, can adjust the local activity for solutes.

Machine learning and polymer self-consistent field theory in two spatial dimensions.

A computational framework that leverages data from self-consistent field theory simulations with deep learning to accelerate the exploration of parameter space for block copolymers is presented. This is a substantial two-dimensional extension of the framework introduced in the work of Xuan et al. [J. Comput. Phys. 443, 110519 (2021)]. Several innovations and improvements are proposed. (1) A Sobolev space-trained, convolutional neural network is employed to handle the exponential dimension increase of the discretized, local average monomer density fields and to strongly enforce both spatial translation and rotation invariance of the predicted, field-theoretic intensive Hamiltonian. (2) A generative adversarial network (GAN) is introduced to efficiently and accurately predict saddle point, local average monomer density fields without resorting to gradient descent methods that employ the training set. This GAN approach yields important savings of both memory and computational cost. (3) The proposed machine learning framework is successfully applied to 2D cell size optimization as a clear illustration of its broad potential to accelerate the exploration of parameter space for discovering polymer nanostructures. Extensions to three-dimensional phase discovery appear to be feasible.

Predicting surfactant phase behavior with a molecularly informed field theory.

HYPOTHESIS: The computational study of surfactants and self-assembly is challenging because 1) models need to reflect chemistry-specific interactions, and 2) self-assembled structures are difficult to equilibrate with conventional molecular dynamics. We propose to overcome these challenges with a multiscale simulation approach where relative entropy minimization transfers chemically-detailed information from all-atom (AA) simulations to coarse-grained (CG) models that can be simulated using field-theoretic methods. Field-theoretic simulations are not limited by intrinsic physical time scales like diffusion and allow for rigorous equilibration via free energy minimization. This approach should enable the study of properties that are difficult to obtain by particle-based simulations. SIMULATION WORK: We apply this workflow to sodium dodecylsulfate. To ensure chemical fidelity we present an AA force field calibrated against interfacial tension experiments. We generate CG models from AA simulation trajectories and show that particle-based and field-theoretic simulations of the CG model reproduce AA simulations and experimental measurements. FINDINGS: The workflow captures the complex balance of interactions in a multicomponent system ultimately described by an atomistic model. The resulting CG models can study complex 3D phases like double or alternating gyroids, and reproduce salt effects on properties like aggregation number and shape transitions.

Multiscale modeling of solute diffusion in triblock copolymer membranes.

We develop a multiscale simulation model for diffusion of solutes through porous triblock copolymer membranes. The approach combines two techniques: self-consistent field theory (SCFT) to predict the structure of the self-assembled, solvated membrane and on-lattice kinetic Monte Carlo (kMC) simulations to model diffusion of solutes. Solvation is simulated in SCFT by constraining the glassy membrane matrix while relaxing the brush-like membrane pore coating against the solvent. The kMC simulations capture the resulting solute spatial distribution and concentration-dependent local diffusivity in the polymer-coated pores; we parameterize the latter using particle-based simulations. We apply our approach to simulate solute diffusion through nonequilibrium morphologies of a model triblock copolymer, and we correlate diffusivity with structural descriptors of the morphologies. We also compare the model's predictions to alternative approaches based on simple lattice random walks and find our multiscale model to be more robust and systematic to parameterize. Our multiscale modeling approach is general and can be readily extended in the future to other chemistries, morphologies, and models for the local solute diffusivity and interactions with the membrane.

Modeling Microstructure Formation in Block Copolymer Membranes Using Dynamical Self-Consistent Field Theory.

Block copolymers have attracted recent interest as candidate materials for ultrafiltration membranes, due to their ability to form isoporous integral-asymmetric membranes by the combined processes of self-assembly and nonsolvent-induced phase separation (SNIPS). However, the dependence of surface layer and substructure morphologies on the processing variables associated with SNIPS is not well understood nor is the interplay between microphase and macrophase separation in block copolymers undergoing such coagulation. Here, we use dynamical self-consistent field theory to simulate the microstructure evolution of block copolymer films during SNIPS and find that such films form the desired sponge-like asymmetric porous substructure only if the solvent and nonsolvent have opposite block selectivities and that otherwise they form a dense nonporous microphase-separated film. Our results could have important implications for the choices of solvent and nonsolvent in the processing of block copolymer membranes.

Field-Theoretic Simulation Method to Study the Liquid-Liquid Phase Separation of Polymers.

Liquid-liquid phase separation (LLPS) is a process that results in the formation of a polymer-rich liquid phase coexisting with a polymer-depleted liquid phase. LLPS plays a critical role in the cell through the formation of membrane-less organelles, but it also has a number of biotechnical and biomedical applications such as drug confinement and its targeted delivery. In this chapter, we present a computational efficient methodology that uses field-theoretic simulations (FTS) with complex Langevin (CL) sampling to characterize polymer phase behavior and delineate the LLPS phase boundaries. This approach is a powerful complement to analytical and explicit-particle simulations, and it can serve to inform experimental LLPS studies. The strength of the method lies in its ability to properly sample a large ensemble of polymers in a saturated solution while including the effect of composition fluctuations on LLPS. We describe the approaches that can be used to accurately construct phase diagrams of a variety of molecularly designed polymers and illustrate the method by generating an approximation-free phase diagram for a classical symmetric diblock polyampholyte.

Phase Separation, Reaction Equilibrium, and Self-Assembly in Binary Telechelic Homopolymer Blends.

We study a binary blend of telechelic homopolymers that can form reversible AB-type bonds at the chain ends. Reversibly bonding polymers display novel material properties, including thermal tunability and self-healing, that are not found in conventional covalently bonded polymers. Previous studies of reversibly bonding polymer systems have been limited by the computational demand of accounting for an infinite number of possible reaction products in a spatially inhomogeneous, self-assembled structure. We demonstrate that newly developed theoretical models and numerical methods enable the simultaneous computation of phase equilibrium, reaction equilibrium, and self-assembly via self-consistent field theory. Phase diagrams are computed at a variety of physically relevant conditions and are compared with nonreactive analogues as well as previous experimental studies of telechelic polymer blends.

Decoupling Ion Transport and Matrix Dynamics to Make High Performance Solid Polymer Electrolytes.

Transport of ions through solid polymeric electrolytes (SPEs) involves a complicated interplay of ion solvation, ion-ion interactions, ion-polymer interactions, and free volume. Nonetheless, prevailing viewpoints on the subject promote a significantly simplified picture, likening ion transport in a polymer to that in an unstructured fluid at low solute concentrations. Although this idealized liquid transport model has been successful in guiding the design of homogeneous electrolytes, structured electrolytes provide a promising alternate route to achieve high ionic conductivity and selectivity. In this perspective, we begin by describing the physical origins of the idealized liquid transport mechanism and then proceed to examine known cases of decoupling between the matrix dynamics and ionic transport in SPEs. Specifically we discuss conditions for "decoupled" mobility that include a highly polar electrolyte environment, a percolated path of free volume elements (either through structured or unstructured channels), high ion concentrations, and labile ion-electrolyte interactions. Finally, we proceed to reflect on the potential of these mechanisms to promote multivalent ion conductivity and the need for research into the interfacial properties of solid polymer electrolytes as well as their performance at elevated potentials.

Ionic Compatibilization of Polymers.

The small specific entropy of mixing of high molecular weight polymers implies that most blends of dissimilar polymers are immiscible with poor physical properties. Historically, a wide range of compatibilization strategies have been pursued, including the addition of copolymers or emulsifiers or installing complementary reactive groups that can promote the in situ formation of block or graft copolymers during blending operations. Typically, such reactive blending exploits reversible or irreversible covalent or hydrogen bonds to produce the desired copolymer, but there are other options. Here, we argue that ionic bonds and electrostatic correlations represent an underutilized tool for polymer compatibilization and in tailoring materials for applications ranging from sustainable polymer alloys to organic electronics and solid polymer electrolytes. The theoretical basis for ionic compatibilization is surveyed and placed in the context of existing experimental literature and emerging classes of functional polymer materials. We conclude with a perspective on how electrostatic interactions might be exploited in plastic waste upcycling.

Direct free energy evaluation of classical and quantum many-body systems via field-theoretic simulation.

Free energy evaluation in molecular simulations of both classical and quantum systems is computationally intensive and requires sophisticated algorithms. This is because free energy depends on the volume of accessible phase space, a quantity that is inextricably linked to the integration measure in a coordinate representation of a many-body problem. In contrast, the same problem expressed as a field theory (auxiliary field or coherent states) isolates the particle number as a simple parameter in the Hamiltonian or action functional and enables the identification of a chemical potential field operator. We show that this feature leads a "direct" method of free energy evaluation, in which a particle model is converted to a field theory and appropriate field operators are averaged using a field-theoretic simulation conducted with complex Langevin sampling. These averages provide an immediate estimate of the Helmholtz free energy in the canonical ensemble and the entropy in the microcanonical ensemble. The method is illustrated for a classical polymer solution, a block copolymer melt exhibiting liquid crystalline and solid mesophases, and a quantum fluid of interacting bosons.

Design of Polymeric Zwitterionic Solid Electrolytes with Superionic Lithium Transport.

Progress toward durable and energy-dense lithium-ion batteries has been hindered by instabilities at electrolyte-electrode interfaces, leading to poor cycling stability, and by safety concerns associated with energy-dense lithium metal anodes. Solid polymeric electrolytes (SPEs) can help mitigate these issues; however, the SPE conductivity is limited by sluggish polymer segmental dynamics. We overcome this limitation via zwitterionic SPEs that self-assemble into superionically conductive domains, permitting decoupling of ion motion and polymer segmental rearrangement. Although crystalline domains are conventionally detrimental to ion conduction in SPEs, we demonstrate that semicrystalline polymer electrolytes with labile ion-ion interactions and tailored ion sizes exhibit excellent lithium conductivity (1.6 mS/cm) and selectivity (t + ≈ 0.6-0.8). This new design paradigm for SPEs allows for simultaneous optimization of previously orthogonal properties, including conductivity, Li selectivity, mechanics, and processability.

Interfacial reaction-induced roughening in polymer thin films.

Reactive blending of immiscible polymers is an important process for synthesizing polymer blends with superior properties. We use a phase-field model to understand reaction dynamics and morphology evolution by diffusive transport in layered films of incompatible, end-reactive polymers. We thoroughly investigate this phenomenon over a large parameter space of interface shapes, layer thicknesses, reaction rates specified by a Damkohler number (Daf), and Flory-Huggins interaction parameter (χ), under static conditions with no external fields. For films of the same thickness, the dynamics of the system is not significantly influenced by the length of the film or the initial shape of the interface. The interface between the polymers is observed to roughen, leading to the formation of a spontaneous emulsion. The reaction progresses slower and the interface roughens later for thicker films, and systems with higher χ. Increasing Daf increases the reaction rate and hastens the onset of roughening. The quasi-static interfacial tension decreases with the extent of reaction, but does not become vanishingly small or negative at the onset of roughening. Simulations with reversible reactions and systems where only a fraction of the homopolymers have reactive end groups show that a critical diblock (reaction product) concentration exists, below which interfacial roughening and spontaneous emulsification is not observed. We also demonstrate that thermal fluctuations accelerate the onset of interfacial roughening, and help sustain the system in an emulsified state.