An official website of the United States government
Modeling the flow of polymer solutions requires knowledge at various length and time scales. The macroscopic behavior is described by the overall velocity, pressure, and stress. The polymeric contribution to the stress requires knowledge of the evolution of polymer chains. In this work, we use a microscopic model, the finitely extensible nonlinear elastic (FENE) model, to capture the polymer’s behavior. The benefit of using microscopic models is that they remain faithful to the polymer dynamics without information loss via averaging. Their downside is the computational cost incurred in solving the thousands to millions of differential equations describing the microstructure. Here, we describe a multiscale flow solver that utilizes GPUs for massively parallel, efficient simulations. We compare and contrast the microscopic model with its macroscopic counterpart under various flow conditions. In particular, significant differences are observed under nonlinear flow conditions, where the polymers become highly stretched and oriented.
more »
« less
This content will become publicly available on April 7, 2026
Efficient dynamical field-theoretic simulations for multi-component systems
Understanding the phase behavior and dynamics of multi-component polymeric systems is essential for designing materials used in applications ranging from biopharmaceuticals to consumer products. While computational tools for understanding the equilibrium properties of such systems are relatively mature, simulation platforms for investigating non-equilibrium behavior are comparatively less developed. Dynamic self-consistent field theory (DSCFT) is a method that retains essential microscopic thermodynamics while enabling a continuum-level understanding of multi-component, multi-phase diffusive transport. A challenge with DSCFT is its high computational complexity and cost, along with the difficulty of incorporating thermal fluctuations. External potential dynamics (EPD) offers a more efficient approach to studying inhomogeneous polymers out of equilibrium, providing similar accuracy to DSCFT but with significantly lower computational cost. In this work, we introduce an extension of EPD to enable efficient and stable simulations of multi-species, multi-component polymer systems while embedding thermodynamically consistent noise. We validate this framework through simulations of a triblock copolymer melt and spinodally decomposing binary and ternary polymer blends, demonstrating its capability to capture key features of phase separation and domain growth. Furthermore, we highlight the role of thermal fluctuations in early stage coarsening. This study provides new insights into the interplay between stochastic and deterministic effects in the dynamic evolution of polymeric fluids, with the EPD framework offering a robust and scalable approach for investigating the complex dynamics of multi-component polymeric materials.
more »
« less
- Award ID(s):
- 2104255
- PAR ID:
- 10633251
- Publisher / Repository:
- American Institute of Physics
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 162
- Issue:
- 13
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Two-dimensional Raman and hybrid terahertz-Raman spectroscopic techniques provide invaluable insight into molecular structures and dynamics of condensed-phase systems. However, corroborating experimental results with theory is difficult due to the high computational cost of incorporating quantum-mechanical effects in the simulations. Here, we present the equilibrium–nonequilibrium ring-polymer molecular dynamics (RPMD), a practical computational method that can account for nuclear quantum effects on the two-time response function of nonlinear optical spectroscopy. Unlike a recently developed approach based on the double Kubo transformed (DKT) correlation function, our method is exact in the classical limit, where it reduces to the established equilibrium-nonequilibrium classical molecular dynamics method. Using benchmark model calculations, we demonstrate the advantages of the equilibrium–nonequilibrium RPMD over classical and DKT-based approaches. Importantly, its derivation, which is based on the nonequilibrium RPMD, obviates the need for identifying an appropriate Kubo transformed correlation function and paves the way for applying real-time path-integral techniques to multidimensional spectroscopy.more » « less
-
The self-assembly of block polymers into well-ordered nanostructures underpins their utility across fundamental and applied polymer science, yet only a handful of equilibrium morphologies are known with the simplest AB-type materials. Here, we report the discovery of the A15 sphere phase in single-component diblock copolymer melts comprising poly(dodecyl acrylate)− block −poly(lactide). A systematic exploration of phase space revealed that A15 forms across a substantial range of minority lactide block volume fractions ( f L = 0.25 − 0.33) situated between the σ-sphere phase and hexagonally close-packed cylinders. Self-consistent field theory rationalizes the thermodynamic stability of A15 as a consequence of extreme conformational asymmetry. The experimentally observed A15−disorder phase transition is not captured using mean-field approximations but instead arises due to composition fluctuations as evidenced by fully fluctuating field-theoretic simulations. This combination of experiments and field-theoretic simulations provides rational design rules that can be used to generate unique, polymer-based mesophases through self-assembly.more » « less
-
Abstract Dynamically cross‐linked polymer networks have attracted significant interest in recent years due to their unique and improved properties including increased toughness, malleability, shape memory, and self‐healing. Here, a computational study on the mechanical behavior of dynamically cross‐linked polymer networks is presented. Coarse grained models for different polymer network configurations are established and their mechanical properties using molecular dynamics (MD) simulations are predicted. Consistent with the experimental measurements, the simulation results show that interpenetrating networks (IPN) withstand considerably higher stress compared to the single networks (SN). Additionally, the MD results demonstrate that the origin of this variation in mechanical behavior of IPN and SN configurations goes back to the difference in spatial degrees of freedom of the noncovalent cross‐linkers, represented by nonbonded interactions within the two system types. The results of this work provide new insight for future studies on the design of systems with dual dynamic cross‐linkers where one linkage exchanges rapidly and provides autonomic dynamic character, while the other is a stimulus responsive dynamic covalent linkage that provides stability with dynamic exchange on‐demand.more » « less
-
Abstract Self‐sorting is commonly observed in complex reaction systems, which has been utilized to guide the formation of single major by‐design molecules. However, most studies have been focused on non‐covalent systems, and using self‐sorting to achieve covalently bonded architectures is still relatively less explored. Herein, we first demonstrated the dynamic nature of spiroborate linkage and systematically studied the self‐sorting behavior observed in the transformation between spiroborate‐linked well‐defined polymeric and molecular architectures, which is enabled by spiroborate bond exchange. The scrambling between a macrocycle and a 1D helical covalent polymer led to the formation of a molecular cage, whose structures are all unambiguously elucidated by single‐crystal X‐ray diffraction. The results indicate that the molecular cage is the thermodynamically favored product in this multi‐component reaction system. This work represents the first example of a 1D polymeric architecture transforming into a shape‐persistent molecular cage, driven by dynamic covalent self‐sorting. This study will further guide the design of spiroborate‐based materials and open the possibilities for the development of novel complex yet responsive dynamic covalent molecular or polymeric systems.more » « less
