Challenges in Multiscale Modeling of Polymer Dynamics

Multiscale modelling of transport in polymer-based reverse-osmosis/nanofiltration membranes: present and future

Nanofiltration (NF) and reverse osmosis (RO) processes are physical separation technologies used to remove contaminants from liquid streams by employing dense polymer-based membranes with nanometric voids that confine fluids at the nanoscale. At this level, physical properties such as solvent and solute permeabilities are intricately linked to molecular interactions. Initially, numerous studies focused on developing macroscopic transport models to gain insights into separation properties at the nanometer scale. However, continuum-based models have limitations in nanoconfined situations that can be overcome by force field molecular simulations. Continuum-based models heavily rely on bulk properties, often neglecting critical factors like liquid structuring, pore geometry, and molecular/chemical specifics. Molecular/mesoscale simulations, while encompassing these details, often face limitations in time and spatial scales. Therefore, achieving a comprehensive understanding of transport requires a synergistic integration of both approaches through a multiscale approach that effectively combines and merges both scales. This review aims to provide a comprehensive overview of the state-of-the-art in multiscale modeling of transport through NF/RO membranes, spanning from the nanoscale to continuum media.

Filled Elastomers: Mechanistic and Physics-Driven Modeling and Applications as Smart Materials

Elastomers are made of chain-like molecules to form networks that can sustain large deformation. Rubbers are thermosetting elastomers that are obtained from irreversible curing reactions. Curing reactions create permanent bonds between the molecular chains. On the other hand, thermoplastic elastomers do not need curing reactions. Incorporation of appropriated filler particles, as has been practiced for decades, can significantly enhance mechanical properties of elastomers. However, there are fundamental questions about polymer matrix composites (PMCs) that still elude complete understanding. This is because the macroscopic properties of PMCs depend not only on the overall volume fraction (ϕ) of the filler particles, but also on their spatial distribution (i.e., primary, secondary, and tertiary structure). This work aims at reviewing how the mechanical properties of PMCs are related to the microstructure of filler particles and to the interaction between filler particles and polymer matrices. Overall, soft rubbery matrices dictate the elasticity/hyperelasticity of the PMCs while the reinforcement involves polymer-particle interactions that can significantly influence the mechanical properties of the polymer matrix interface. For ϕ values higher than a threshold, percolation of the filler particles can lead to significant reinforcement. While viscoelastic behavior may be attributed to the soft rubbery component, inelastic behaviors like the Mullins and Payne effects are highly correlated to the microstructures of the polymer matrix and the filler particles, as well as that of the polymer-particle interface. Additionally, the incorporation of specific filler particles within intelligently designed polymer systems has been shown to yield a variety of functional and responsive materials, commonly termed smart materials. We review three types of smart PMCs, i.e., magnetoelastic (M-), shape-memory (SM-), and self-healing (SH-) PMCs, and discuss the constitutive models for these smart materials.

Chemically Specific Systematic Coarse-Grained Polymer Model with Both Consistently Structural and Dynamical Properties

The coarse-grained (CG) model serves as a powerful tool for the simulation of polymer systems; its reliability depends on the accurate representation of both structural and dynamical properties. However, strong correlations between structural and dynamical properties on different scales and also a strong memory effect, enforced by chain connectivity between monomers in polymer systems, render developing a chemically specific systematic CG model a formidable task. In this study, we report a systematic CG approach that combines the iterative Boltzmann inversion (IBI) method and the generalized Langevin equation (GLE) dynamics. Structural properties are ensured by using conservative CG potentials derived from the IBI method. To retrieve the correct dynamical properties in the system, we demonstrate that using a combination of a Rouse-type delta function and a time-dependent short-time kernel in the GLE simulation is practically efficient. The former can be used to adjust the long-time diffusion dynamics, and the latter can be reconstructed from an iterative procedure according to the velocity autocorrelation function (ACF) from all-atomistic (AA) simulations. Taking the polystyrene as an example, we show that not only structural properties of radial distribution function, intramolecular bond, and angle distributions can be reproduced but also dynamical properties of mean-square displacement, velocity ACF, and force ACF resulted from our CG model have quantitative agreement with the reference AA model. In addition, reasonable agreements are observed in other collective properties between our GLE-CG model and the AA simulations as well.

Insight into the Density-Dependence of Pair Potentials for Predictive Coarse-Grained Models

We investigate the temperature- and density-dependence of effective pair potentials for 1-site coarse-grained (CG) models of two industrial solvents, 1,4-dioxane and tetrahydrofuran. We observe that the calculated pair potentials are much more sensitive to density than to temperature. The generalized-Yvon-Born-Green framework reveals that this striking density-dependence reflects corresponding variations in the many-body correlations that determine the environment-mediated indirect contribution to the pair mean force. Moreover, we demonstrate, perhaps surprisingly, that this density-dependence is not important for accurately modeling the intermolecular structure. Accordingly, we adopt a density-independent interaction potential and transfer the density-dependence of the calculated pair potentials into a configuration-independent volume potential. Furthermore, we develop a single global potential that accurately models the intermolecular structure and pressure-volume equation of state across a very wide range of liquid state points. Consequently, this work provides fundamental insight into the density-dependence of effective pair potentials and also provides a significant step toward developing predictive CG models for efficiently modeling industrial solvents.

Predicting Natural Rubber Crystallinity by a Novel Machine Learning Algorithm Based on Molecular Dynamics Simulation Data

Natural rubber (NR) with excellent mechanical properties, mainly attributed to its strain-induced crystallization (SIC), has garnered significant scientific and technological interest. With the aid of molecular dynamics (MD) simulations, we can investigate the impacts of crucial structural elements on SIC on the molecular scale. Nonetheless, the computational complexity and time-consuming nature of this high-precision method constrain its widespread application. The integration of machine learning with MD represents a promising avenue for enhancing the speed of simulations while maintaining accuracy. Herein, we developed a crystallinity algorithm tailored to the SIC properties of natural rubber materials. With the data enhancement algorithm, the high evaluation value of the prediction model ensures the accuracy of the computational simulation results. In contrast to the direct utilization of small sample prediction algorithms, we propose a novel concept grounded in feature engineering. The proposed machine learning (ML) methodology consists of (1) An eXtreme Gradient Boosting (XGB) model to predict the crystallinity of NR; (2) a generative adversarial network (GAN) data augmentation algorithm to optimize the utilization of the limited training data, which is utilized to construct the XGB prediction model; (3) an elaboration of the effects induced by phospholipid and protein percentage (ω), hydrogen bond strength (εH), and non-hydrogen bond strength (εNH) of natural rubber materials with crystallinity prediction under dynamic conditions are analyzed by employing weight integration with feature importance analysis. Eventually, we succeeded in concluding that εH has the most significant effect on the strain-induced crystallinity, followed by ω and finally εNH.

Visualization of Self-Assembly and Hydration of a β-Hairpin through Integrated Small and Wide-Angle Neutron Scattering

Fundamental understanding of the structure and assembly of nanoscale building blocks is crucial for the development of novel biomaterials with defined architectures and function. However, accessing self-consistent structural information across multiple length scales is challenging. This limits opportunities to exploit atomic scale interactions to achieve emergent macroscale properties. In this work we present an integrative small- and wide-angle neutron scattering approach coupled with computational modeling to reveal the multiscale structure of hierarchically self-assembled β hairpins in aqueous solution across 4 orders of magnitude in length scale from 0.1 Å to 300 nm. Our results demonstrate the power of this self-consistent cross-length scale approach and allows us to model both the large-scale self-assembly and small-scale hairpin hydration of the model β hairpin CLN025. Using this combination of techniques, we map the hydrophobic/hydrophilic character of this model self-assembled biomolecular surface with atomic resolution. These results have important implications for the multiscale investigation of aqueous peptides and proteins, for the prediction of ligand binding and molecular associations for drug design, and for understanding the self-assembly of peptides and proteins for functional biomaterials.

Interaction of Macromolecular Chain with Phospholipid Membranes in Solutions: A Dissipative Particle Dynamics Simulation Study

The interaction between macromolecular chains and phospholipid membranes in aqueous solution was investigated using dissipative particle dynamics simulations. Two cases were considered, one in which the macromolecular chains were pulled along parallel to the membrane surfaces and another in which they were pulled vertical to the membrane surfaces. Several parameters, including the radius of gyration, shape factor, particle number, and order parameter, were used to investigate the interaction mechanisms during the dynamics processes by adjusting the pulling force strength of the chains. In both cases, the results showed that the macromolecular chains undergo conformational transitions from a coiled to a rod-like structure. Furthermore, the simulations revealed that the membranes can be damaged and repaired during the dynamic processes. The role of the pulling forces and the adsorption interactions between the chains and membranes differed in the parallel and perpendicular pulling cases. These findings contribute to our understanding of the interaction mechanisms between macromolecules and membranes, and they may have potential applications in biology and medicine.

Integrating Machine Learning in the Coarse-Grained Molecular Simulation of Polymers

Machine learning (ML) is having an increasing impact on the physical sciences, engineering, and technology and its integration into molecular simulation frameworks holds great potential to expand their scope of applicability to complex materials and facilitate fundamental knowledge and reliable property predictions, contributing to the development of efficient materials design routes. The application of ML in materials informatics in general, and polymer informatics in particular, has led to interesting results, however great untapped potential lies in the integration of ML techniques into the multiscale molecular simulation methods for the study of macromolecular systems, specifically in the context of Coarse Grained (CG) simulations. In this Perspective, we aim at presenting the pioneering recent research efforts in this direction and discussing how these new ML-based techniques can contribute to critical aspects of the development of multiscale molecular simulation methods for bulk complex chemical systems, especially polymers. Prerequisites for the implementation of such ML-integrated methods and open challenges that need to be met toward the development of general systematic ML-based coarse graining schemes for polymers are discussed.

Approaches and Perspective of Coarse-Grained Modeling and Simulation for Polymer-Nanoparticle Hybrid Systems

Molecular modeling and simulations have emerged as effective and indispensable tools to characterize polymeric systems. They provide fundamental and essential insights to design a product of the required properties and to improve the understanding of a phenomenon at the molecular level for a particular system. The polymer-nanoparticle hybrids are materials with outstanding properties and correspondingly large applications whose study has benefited from this new paradigm. However, despite the significant expansion of modern day computational powers, investigation of the long time and large length scale phenomenon in polymeric and polymer-nanoparticle systems is still a challenging task to complete through all-atom molecular dynamics (AA-MD) simulations. To circumvent this problem, a variety of coarse-grained (CG) models have been proposed, ranging from the generic CG models for qualitative properties predictions to more realistic chemically specific CG models for quantitative properties predictions. These CG models have already delivered some success stories in the study of several spatial and temporal evolutions of many processes. Some of these studies were beyond the feasibility of traditional atomistic resolution models due to either the size or the time constraints. This review captures the different types of popular CG approaches that are utilized in the investigation of the microscopic behavior of polymer-nanoparticle hybrid systems. The rationale of this article is to furnish an overview of the popular CG approaches and their applications, to review several important and most recent developments, and to delineate the perspectives on future directions in the field.

Polymer solution structure and dynamics within pores of hexagonally close-packed nanoparticles

Using coarse-grained molecular dynamics simulations, we examine structure and dynamics of polymer solutions under confinement within the pores of a hexagonally close-packed (HCP) nanoparticle system with nanoparticle diameter fifty times that of the polymer Kuhn segment size. We model a condition where the polymer chain is in a good solvent (i.e., polymer-polymer interaction is purely repulsive and polymer-solvent and solvent-solvent interactions are attractive) and the polymer-nanoparticle and solvent-nanoparticle interactions are purely repulsive. We probe three polymer lengths (N = 10, 114, and 228 Kuhn segments) and three solution concentrations (1, 10, and 25%v) to understand how the polymer chain conformations and chain center-of-mass diffusion change under confinement within the pores of the HCP nanoparticle structure from those seen in bulk. The known trend of bulk polymer R_g² decreasing with increasing concentration no longer holds when confined in the pores of HCP nanoparticle structure; for example, for the 114-mer, the HCP 〈R_g²〉 at 1%v concentration is lower than HCP 〈R_g²〉 at 10%v concentration. The 〈R_g²〉 of the 114-mer and 228-mer exhibit the largest percent decline going from bulk to HCP at the 1%v concentration and the smallest percent decline at the 25%v concentration. We also provide insight into how the confinement ratio (CR) of polymer chain size to pore size within tetrahedral and octahedral pores in the HCP arrangement of nanoparticles affects the chain conformation and diffusion at various concentrations. At the same concentration, the N = 114 has significantly more movement between pores than the N = 228 chains. For the N = 114 polymer, the diffusion between pores (i.e., inter-pore diffusion) accelerates the overall diffusion rate for the confined HCP system while for the N = 228 polymer, the polymer diffusion in the entire HCP is dominated by the diffusion within the tetrahedral or octahedral pores with minor contributions from inter-pore diffusion. These findings augment the fundamental understanding of macromolecular diffusion through large, densely packed nanoparticle assemblies and are relevant to research focused on fabrication of polymer composite materials for chemical separations, storage, optics, and photonics. We perform coarse-grained molecular dynamics simulations to understand structure and dynamics of polymer solutions under confinement within hexagonal close packed nanoparticles with radii much larger than the polymer chain's bulk radius of gyration.

Self-Assembly of Lipid Mixtures in Solutions: Structures, Dynamics Processes and Mechanical Properties

The self-assembly of lipid mixtures in aqueous solution was investigated by dissipative particle dynamics simulation. Two types of lipid molecules were modelled, where three mixed structures, i.e., the membrane, perforated membrane and vesicle, were determined in the self-assembly processes. Phase behaviour was investigated by using the phase diagrams based on the tail chain lengths for the two types of lipids. Several parameters, such as chain number and average radius of gyration, were employed to explore the structural formations of the membrane and perforated membrane in the dynamic processes. Interface tension was used to demonstrate the mechanical properties of the membrane and perforated membrane in the equilibrium state and dynamics processes. Results help us to understand the self-assembly mechanism of the biomolecule mixtures, which has a potential application for designing the lipid molecule-based bio-membranes in solutions.

Coarse-Grained Model Incorporating Short- and Long-Range Effective Potentials for the Fast Simulation of Micelle Formation in Solutions of Ionic Surfactants

A coarse-grained model comprising short- and long-range effective potentials, parametrized with the iterative Boltzmann inversion (IBI) method, is presented for capturing micelle formation in aqueous solutions of ionic surfactants using as a model system sodium dodecyl sulfate (SDS). In the coarse-grained (CG) model, each SDS molecule is represented as a sequence of four beads while each water molecule is modeled as a single bead. The proposed CG scheme involves ten potential energy functions: four of them describe bonded interactions and control the distribution functions of intramolecular degrees of freedom (bond lengths, valence angles, and dihedrals) along an SDS molecule while the other six account for intermolecular interactions between pairs of SDS and water beads and control the radial distribution functions. The nonbonded effective potentials between coarse-grained SDS molecules extend up to about 12 nm and capture structural and morphological features of the micellar solution both at short and long distances. The long-range component of these potentials, in particular, captures correlations between surfactant molecules belonging to different micelles and is essential to describe ordering associated with micelle formation. A new strategy is introduced for determining the effective potentials through IBI by using information (target distribution functions) extracted from independent atomistic simulations of a micellar reference system (a salt-free SDS solution at total surfactant concentration c_T equal to 103 mM, temperature T equal to 300 K, and pressure P equal to 1 atm) obtained through a multiscale approach described in an earlier study. It employs several optimization steps for bonded and nonbonded interactions and a gradual parametrization of the short- and long-range components of the latter, followed by reparametrization of the bonded ones. The proposed CG model can reproduce remarkably accurately the microstructure and morphology of the reference system within only a few hours of computational time. It is therefore very promising for future studies of structural and morphological behavior of various liquid surfactant formulations.

Systematic bottom-up molecular coarse-graining via force and torque matching using anisotropic particles

We derive a systematic and general method for parametrizing coarse-grained molecular models consisting of anisotropic particles from fine-grained (e.g. all-atom) models for condensed-phase molecular dynamics simulations. The method, which we call anisotropic force-matching coarse-graining (AFM-CG), is based on rigorous statistical mechanical principles, enforcing consistency between the coarse-grained and fine-grained phase-space distributions to derive equations for the coarse-grained forces, torques, masses, and moments of inertia in terms of properties of a condensed-phase fine-grained system. We verify the accuracy and efficiency of the method by coarse-graining liquid-state systems of two different anisotropic organic molecules, benzene and perylene, and show that the parametrized coarse-grained models more accurately describe properties of these systems than previous anisotropic coarse-grained models parametrized using other methods that do not account for finite-temperature and many-body effects on the condensed-phase coarse-grained interactions. The AFM-CG method will be useful for developing accurate and efficient dynamical simulation models of condensed-phase systems of molecules consisting of large, rigid, anisotropic fragments, such as liquid crystals, organic semiconductors, and nucleic acids.

Predicting Young's Modulus of Linear Polyurethane and Polyurethane-Polyurea Elastomers: Bridging Length Scales with Physicochemical Modeling and Machine Learning

Predicting the properties of complex polymeric materials based on monomer chemistry requires modeling physical interactions that bridge molecular, interchain, microstructure, and bulk length scales. For polyurethanes, a polymer class with global commercial and industrial significance, these multiscale challenges are intrinsic due to the thermodynamic incompatibility of the urethane and polyol-rich domains, resulting in heterogeneities from molecular to microstructural length scales. Machine learning can model patterns in data to establish a relationship between the monomer chemistry and bulk material properties, but this is made difficult by small data sets and a diverse set of monomers. Using a data set of 63 industrially relevant and complex elastomers, we demonstrate that accurate machine learning predictions are possible when monomer chemistry is used to estimate interactions at interchain length scales. Here, these features were used to accurately (r² = 0.91) predict the Young's modulus of polyurethane and polyurethane-urea elastomers. Furthermore, by a query of the trained model for compositions that yield a target modulus within the range of accessible values, the capabilities of using this methodology as a design tool are demonstrated. The presented methodology could become increasingly useful in building models for materials with small data sets and may guide the interpretation of the underlying physicochemical forces.

Polymer crystallization under external flow

The general aspects of polymer crystallization under external flow, i.e., flow-induced crystallization (FIC) from fundamental theoretical background to multi-scale characterization and modeling results are presented. FIC is crucial for modern polymer processing, such as blowing, casting, and injection modeling, as two-third of daily-used polymers is crystalline, and nearly all of them need to be processed before final applications. For academics, the FIC is intrinsically far from equilibrium, where the polymer crystallization behavior is different from that in quiescent conditions. The continuous investigation of crystallization contributes to a better understanding on the general non-equilibrium ordering in condensed physics. In the current review, the general theories related to polymer nucleation under flow (FIN) were summarized first as a preliminary knowledge. Various theories and models, i.e., coil-stretch transition and entropy reduction model, are briefly presented together with the modified versions. Subsequently, the multi-step ordering process of FIC is discussed in detail, including chain extension, conformational ordering, density fluctuation, and final perfection of the polymer crystalline. These achievements for a thorough understanding of the fundamental basis of FIC benefit from the development of various hyphenated rheometer, i.e., rheo-optical spectroscopy, rheo-IR, and rheo-x-ray scattering. The selected experimental results are introduced to present efforts on elucidating the multi-step and hierarchical structure transition during FIC. Then, the multi-scale modeling methods are summarized, including micro/meso scale simulation and macroscopic continuum modeling. At last, we briefly describe our personal opinions related to the future directions of this field, aiming to ultimately establish the unified theory of FIC and promote building of the more applicable models in the polymer processing.

Integration of Machine Learning and Coarse-Grained Molecular Simulations for Polymer Materials: Physical Understandings and Molecular Design

In recent years, the synthesis of monomer sequence-defined polymers has expanded into broad-spectrum applications in biomedical, chemical, and materials science fields. Pursuing the characterization and inverse design of these polymer systems requires our fundamental understanding not only at the individual monomer level, but also considering the chain scales, such as polymer configuration, self-assembly, and phase separation. However, our accessibility to this field is still rudimentary due to the limitations of traditional design approaches, the complexity of chemical space along with the burdened cost and time issues that prevent us from unveiling the underlying monomer sequence-structure-property relationships. Fortunately, thanks to the recent advancements in molecular dynamics simulations and machine learning (ML) algorithms, the bottlenecks in the tasks of establishing the structure-function correlation of the polymer chains can be overcome. In this review, we will discuss the applications of the integration between ML techniques and coarse-grained molecular dynamics (CGMD) simulations to solve the current issues in polymer science at the chain level. In particular, we focus on the case studies in three important topics—polymeric configuration characterization, feed-forward property prediction, and inverse design—in which CGMD simulations are leveraged to generate training datasets to develop ML-based surrogate models for specific polymer systems and designs. By doing so, this computational hybridization allows us to well establish the monomer sequence-functional behavior relationship of the polymers as well as guide us toward the best polymer chain candidates for the inverse design in undiscovered chemical space with reasonable computational cost and time. Even though there are still limitations and challenges ahead in this field, we finally conclude that this CGMD/ML integration is very promising, not only in the attempt of bridging the monomeric and macroscopic characterizations of polymer materials, but also enabling further tailored designs for sequence-specific polymers with superior properties in many practical applications.

Influence of the Filler Particles' Surface Morphology on the Polyurethane Matrix's Structure Formation in the Composite

This article presents the surface morphology effect of silicon carbide (SiC) particles on the polyurethane binder’s structure formation in a dispersed-filled composite. The difference in the morphology and surface relief of filler particles was ensured by the implementation of plasma chemical modification. As a result of this modification, the filler consisted of core-shell particles characterized by a SiC core and a carbon shell (SiC@C), as well as a carbon shell decorated with silicon nanoparticles (SiC@C/SiNP) or nanos (SiC@C/SiNW). The study of the relaxation properties of polyurethane composites has shown that the strongest limiting effect on the molecular mobility of boundary layer’s chain segments is exerted by a highly developed surface with a complex relief of SiC@C/SiNP and SiC@C/SiNW particles. An empirical method was proposed to find the polymer fractions spent on the formation of the boundary, transition and bulk layers of the polymer matrix in the composite. It was shown that the morphology of the filler particles’ surface does not affect the dependence of the boundary layer thickness on the filler’s volume fraction. However, with an increase in the degree of surface development, the boundary layer thickness decreases.

Machine-Learning-Guided Discovery of 19F MRI Agents Enabled by Automated Copolymer Synthesis

Modern polymer science suffers from the curse of multidimensionality. The large chemical space imposed by including combinations of monomers into a statistical copolymer overwhelms polymer synthesis and characterization technology and limits the ability to systematically study structure-property relationships. To tackle this challenge in the context of 19F magnetic resonance imaging (MRI) agents, we pursued a computer-guided materials discovery approach that combines synergistic innovations in automated flow synthesis and machine learning (ML) method development. A software-controlled, continuous polymer synthesis platform was developed to enable iterative experimental-computational cycles that resulted in the synthesis of 397 unique copolymer compositions within a six-variable compositional space. The nonintuitive design criteria identified by ML, which were accomplished by exploring <0.9% of the overall compositional space, lead to the identification of >10 copolymer compositions that outperformed state-of-the-art materials.

Wetting characteristics of polymer adhesives with different chain bending stiffness

Polymer adhesives are widely used in daily applications and in industry owing to their flexibility and overall non-toxicity, particularly in interfacial adhesion. The spreading of polymer adhesives on adherend is one of the essential considerations for the interfacial adhesion of polymer adhesives, which is strongly related to their wetting behaviors. While relationships between polymer microstructure and adhesion have been investigated in previous studies, it remains challenging to unveil the effect of polymer microstructure on wettability. To address this issue, here we utilize coarse-grained molecular dynamics (CGMD) simulations to systematically elucidate how the wettability of a polymer adhesive droplet on a surface depends on bending stiffness. The wetting dynamics and the contact angle are studied to show the evolution of morphology of droplets during the wetting process. The results indicate the wettability is weakened by the increase of bending stiffness of polymer chain. Detailed thermodynamic property analysis is further conducted, revealing that the adhesion between the polymer droplet and substrate deteriorates due to the decline of wettability. Interestingly, we observe such deterioration becomes more significant by both increasing the temperature and decreasing the bending stiffness. Our study sheds light on the dependence of chain bending stiffness and temperature on the wetting behavior of polymer adhesive droplets, and offers insights, which, upon experimental validation can then be used for the design of adhesives or hydrogels.

Assessment of the PETase conformational changes induced by poly(ethylene terephthalate) binding

Recently, a bacterium strain of Ideonella sakaiensis was identified with the uncommon ability to degrade the poly(ethylene terephthalate) (PET). The PETase from I. sakaiensis strain 201-F6 (IsPETase) catalyzes the hydrolysis of PET converting it to mono(2-hydroxyethyl) terephthalic acid (MHET), bis(2-hydroxyethyl)-TPA (BHET), and terephthalic acid (TPA). Despite the potential of this enzyme for mitigation or elimination of environmental contaminants, one of the limitations of the use of IsPETase for PET degradation is the fact that it acts only at moderate temperature due to its low thermal stability. Besides, molecular details of the main interactions of PET in the active site of IsPETase remain unclear. Herein, molecular docking and molecular dynamics (MD) simulations were applied to analyze structural changes of IsPETase induced by PET binding. Results from the essential dynamics revealed that the β1-β2 connecting loop is very flexible. This loop is located far from the active site of IsPETase and we suggest that it can be considered for mutagenesis to increase the thermal stability of IsPETase. The free energy landscape (FEL) demonstrates that the main change in the transition between the unbound to the bound state is associated with the β7-α5 connecting loop, where the catalytic residue Asp206 is located. Overall, the present study provides insights into the molecular binding mechanism of PET into the IsPETase structure and a computational strategy for mapping flexible regions of this enzyme, which can be useful for the engineering of more efficient enzymes for recycling plastic polymers using biological systems.

A survey of multiscale modeling: Foundations, historical milestones, current status, and future prospects

Research problems in the domains of physical, engineering, biological sciences often span multiple time and length scales, owing to the complexity of information transfer underlying mechanisms. Multiscale modeling (MSM) and high-performance computing (HPC) have emerged as indispensable tools for tackling such complex problems. We review the foundations, historical developments, and current paradigms in MSM. A paradigm shift in MSM implementations is being fueled by the rapid advances and emerging paradigms in HPC at the dawn of exascale computing. Moreover, amidst the explosion of data science, engineering, and medicine, machine learning (ML) integrated with MSM is poised to enhance the capabilities of standard MSM approaches significantly, particularly in the face of increasing problem complexity. The potential to blend MSM, HPC, and ML presents opportunities for unbound innovation and promises to represent the future of MSM and explainable ML that will likely define the fields in the 21st century.

Optimizing the heterogeneous network structure to achieve polymer nanocomposites with excellent mechanical properties

Designing and optimizing the polymer network structure at the molecular level to manipulate its mechanical properties are of great scientific significance. Although heterogeneous multi-network structures have been extensively investigated, little effort has been devoted to investigating heterogeneous single-networks with a well-defined interface. Herein, through coarse-grained molecular dynamics simulation, we successfully fabricated a heterogeneous single-network, which was divided into several regions with different crosslink densities. Firstly, we found that there is an optimal crosslink density ratio between high and low crosslink density regions to obtain the best stress-strain behavior. Secondly, the effect of the regularity of the network topology (by changing the distribution of two-phase regions) on mechanical properties was also studied. It was clearly observed that the polymer network showed better elastic response and mechanical properties as the distribution of two-phase regions became uniform. Finally, we investigated the effect of the selective distribution of nanoparticles (NPs) on mechanical properties by introducing NPs into a pre-designed multiphase network. Results showed that the selective distribution of NPs in the high crosslink density region had a more significant effect on the mechanical reinforcement. Generally, our simulated results may provide some guidelines to design polymer network structures to achieve high-performance polymer nanocomposites with excellent mechanical properties.

A generalized tube model of rubber elasticity

In the present paper, a new type of micro-mechanically motivated chain network model for rubber-like materials is proposed. The model captures topological constraints of polymer network chains, in particular, entanglements. The model demonstrates how the local molecular packing constraints modify under deformation and shows the impact of these changes on the macroscopic elasticity of the material. To this end, we combine concepts of a confining tube and a slip-link (reptation) model. In these models, entanglements of polymer chains play an important role. The nature of entanglements is discussed, and relationships governing entanglements are formulated in terms of molecular physics. In the context of nonlinear elasticity, we apply a non-affine concept which captures the liquid-like behavior of polymer networks at smaller scales in a more realistic way. Model predictions show good agreement with experimental results from uniaxial and biaxial tension tests.

Multiscale Discrete Element Modeling

A multiscale approach to discrete element modeling is presented. A distinctive feature of the method is that each macroscopic discrete element has an associated atomic sample representing the material’s atomic structure. The dynamics of the elements on macro and micro levels are described by systems of ordinary differential equations, which are solved in a self-consistent manner. A full cycle of multiscale simulations is applied to polycrystalline silicon. Macroscale elastic properties of silicon were obtained only using data extracted from the quantum mechanical properties. The results of computational experiments correspond well to the reference data.

Machine Learning of Coarse-Grained Models for Organic Molecules and Polymers: Progress, Opportunities, and Challenges

Machine learning (ML) has emerged as one of the most powerful tools transforming all areas of science and engineering. The nature of molecular dynamics (MD) simulations, complex and time-consuming calculations, makes them particularly suitable for ML research. This review article focuses on recent advancements in developing efficient and accurate coarse-grained (CG) models using various ML methods, in terms of regulating the coarse-graining process, constructing adequate descriptors/features, generating representative training data sets, and optimization of the loss function. Two classes of the CG models are introduced: bottom-up and top-down CG methods. To illustrate these methods and demonstrate the open methodological questions, we survey several important principles in constructing CG models and how these are incorporated into ML methods and improved with specific learning techniques. Finally, we discuss some key aspects of developing machine-learned CG models with high accuracy and efficiency. Besides, we describe how these aspects are tackled in state-of-the-art methods and which remain to be addressed in the near future. We expect that these machine-learned CG models can address thermodynamic consistent, transferable, and representative issues in classical CG models.

Development of Coarse-Grained Models for Poly(4-vinylphenol) and Poly(2-vinylpyridine): Polymer Chemistries with Hydrogen Bonding

In this paper, we identify the modifications needed in a recently developed generic coarse-grained (CG) model that captured directional interactions in polymers to specifically represent two exemplary hydrogen bonding polymer chemistries-poly(4-vinylphenol) and poly(2-vinylpyridine). We use atomistically observed monomer-level structures (e.g., bond, angle and torsion distribution) and chain structures (e.g., end-to-end distance distribution and persistence length) of poly(4-vinylphenol) and poly(2-vinylpyridine) in an explicitly represented good solvent (tetrahydrofuran) to identify the appropriate modifications in the generic CG model in implicit solvent. For both chemistries, the modified CG model is developed based on atomistic simulations of a single 24-mer chain. This modified CG model is then used to simulate longer (36-mer) and shorter (18-mer and 12-mer) chain lengths and compared against the corresponding atomistic simulation results. We find that with one to two simple modifications (e.g., incorporating intra-chain attraction, torsional constraint) to the generic CG model, we are able to reproduce atomistically observed bond, angle and torsion distributions, persistence length, and end-to-end distance distribution for chain lengths ranging from 12 to 36 monomers. We also show that this modified CG model, meant to reproduce atomistic structure, does not reproduce atomistically observed chain relaxation and hydrogen bond dynamics, as expected. Simulations with the modified CG model have significantly faster chain relaxation than atomistic simulations and slower decorrelation of formed hydrogen bonds than in atomistic simulations, with no apparent dependence on chain length.

A coarse-grain force field based on quantum mechanics (CGq FF) for molecular dynamics simulation of poly(ethylene glycol)-block-poly(ε-caprolactone) (PEG-b-PCL) micelles

In order to provide the means to predict from molecular dynamics (MD) simulations the structures of copolymer-based micelles in solution, we developed coarse grain force field (CGq FF) parameters for poly(ethylene glycol) (PEG) and for poly(ε-caprolactone) (PCL). A key advance here is the use of quantum mechanics to train the parameters describing the non-bonded (NB) interactions between the CG beads. The functional forms are the same as the MARTINI CG FF so standard MD codes can be used. Our CGq FF describes well the experimentally observed properties for the polymer-air and polymer-water interfaces, indicating the accuracy of the NB interactions. The structural properties (density, radius of gyration (R_g), and end-to-end distance (h)) match both experiment and all atom (AA) simulations. We illustrate the application of this CGq FF by following the formation of a spherical micelle from 250 chains of PEG₂₃-b-PCL₉ diblock copolymer, each block with molecular weight of 1000 Daltons (10 500 beads, corresponding to 123 250 atoms), in a water box with 119 139 water beads (426 553 water molecules).

Experimental Validation of Injection Molding Simulations of 3D Microparts and Microstructured Components Using Virtual Design of Experiments and Multi-Scale Modeling

The increasing demand for micro-injection molding process technology and the corresponding micro-molded products have materialized in the need for models and simulation capabilities for the establishment of a digital twin of the manufacturing process. The opportunities enabled by the correct process simulation include the possibility of forecasting the part quality and finding optimal process conditions for a given product. The present work displays further use of micro-injection molding process simulation for the prediction of feature dimensions and its optimization and microfeature replication behavior due to geometrical boundary effects. The current work focused on the micro-injection molding of three-dimensional microparts and of single components featuring microstructures. First, two virtual a studies were performed to predict the outer diameter of a micro-ring within an accuracy of 10 µm and the flash formation on a micro-component with mass a 0.1 mg. In the second part of the study, the influence of microstructure orientation on the filling time of a microcavity design section was investigated for a component featuring micro grooves with a 15 µm nominal height. Multiscale meshing was employed to model the replication of microfeatures in a range of 17–346 µm in a Fresnel lens product, allowing the prediction of the replication behavior of a microfeature at 91% accuracy. The simulations were performed using 3D modeling and generalized Navier–Stokes equations using a single multi-scale simulation approach. The current work shows the current potential and limitations in the use of micro-injection molding process simulations for the optimization of micro 3D-part and microstructured components.

Unified Analytic Expressions for the Entanglement Length, Tube Diameter, and Plateau Modulus of Polymer Melts

By combining molecular dynamics simulations and topological analyses with scaling arguments, we obtain analytic expressions that quantitatively predict the entanglement length N_{e}, the plateau modulus G, and the tube diameter a in melts that span the entire range of chain stiffnesses for which systems remain isotropic. Our expressions resolve conflicts between previous scaling predictions for the loosely entangled [Lin-Noolandi, Gℓ_{K}^{3}/k_{B}T∼(ℓ_{K}/p)^{3}], semiflexible [Edwards-de Gennes: Gℓ_{K}^{3}/k_{B}T∼(ℓ_{K}/p)^{2}], and tightly entangled [Morse, Gℓ_{K}^{3}/k_{B}T∼(ℓ_{K}/p)^{1+ϵ}] regimes, where ℓ_{K} and p are, respectively, the Kuhn and packing lengths. We also find that maximal entanglement (minimal N_{e}) coincides with the onset of local nematic order.

Influence of Small-Scale Correlation on the Interface Evolution of Semiflexible Homopolymer Blends

Within the framework of a dynamic self-consistent field theory, we study the effect of the correlations in a small scale on polymer dynamics, adopting the semiflexible homopolymer blends as the model system. This is accomplished by taking the pair correlation function of ideal semiflexible chains as the Onsager coefficient and the Debye function as an approximation to the Onsager coefficient. Relying on the difference of the two pair correlation functions in the small-scale region, we can identify the effect of small-scale correlations. In the equilibrium state, with the chain length growing, the interface width has a continuous transition from the contour length to radius of gyration. The investigation of interfacial evolution and chain orientation reveals that strong small-scale correlations would accelerate the small-scale dynamic process. We also expect that such a small-scale effect should be highlighted in the process where microscopic phase separation happens.

Multi-scale modeling of the polymer-filler interaction

We report mesoscopic simulations of the interaction between a silica nanoparticle and cis-1,4-polybutadiene chains with realistic coarse-(CG) grained models. The CG models are obtained with a bottom-up Bayesian method based on trajectory matching of atomistic configurations of the system. We then investigate the structural properties of the interfacial region as a function of the grafting density and polymer chain length. We take advantage of the realistic CG models to explore the dynamics of the nanoparticle over a period of 10 microseconds. We show that the dynamics of the nanoparticle is affected by the grafting density and the polymer chain length of the grafted chains.

Dual-potential approach for coarse-grained implicit solvent models with accurate, internally consistent energetics and predictive transferability

The dual-potential approach promises coarse-grained (CG) models that accurately reproduce both structural and energetic properties, while simultaneously providing predictive estimates for the temperature-dependence of the effective CG potentials. In this work, we examine the dual-potential approach for implicit solvent CG models that reflect large entropic effects from the eliminated solvent. Specifically, we construct implicit solvent models at various resolutions, R, by retaining a fraction 0.10 ≤ R ≤ 0.95 of the molecules from a simple fluid of Lennard-Jones spheres. We consider the dual-potential approach in both the constant volume and constant pressure ensembles across a relatively wide range of temperatures. We approximate the many-body potential of mean force for the remaining solutes with pair and volume potentials, which we determine via multiscale coarse-graining and self-consistent pressure-matching, respectively. Interestingly, with increasing temperature, the pair potentials appear increasingly attractive, while the volume potentials become increasingly repulsive. The dual-potential approach not only reproduces the atomic energetics but also quite accurately predicts this temperature-dependence. We also derive an exact relationship between the thermodynamic specific heat of an atomic model and the energetic fluctuations that are observable at the CG resolution. With this generalized fluctuation relationship, the approximate CG models quite accurately reproduce the thermodynamic specific heat of the underlying atomic model.

Age-aware constitutive materials model for a 3D printed polymeric foam

Traditional open or closed-cell stochastic elastomeric foams have wide-ranging applications in numerous industries: from thermal insulation, shock absorbing/gap-filling support cushions, packaging, to light-weight structural and positional components. Recent developments in 3D printing technologies by direct ink-write have opened the possibility of replacing stochastic foam parts by more controlled printed micro-structures with superior stress-distribution and longer functional life. For successful deployment as mechanical support or structural components, it is crucial to characterize the response of such printed materials to long-term external loads in terms of stress-strain behavior evolution and in terms of irreversible structural and load-bearing capacity changes over time. To this end, here we report a thermal-age-aware constitutive model for a 3D printed close-packed foam structure under compression. The model is based on the Ogden hyperfoam strain-energy functional within the framework of Tobolsky two-network scheme. It accurately describes experimentally measured stress-strain response, compression set, and load retention for various aging times and temperatures. Through the technique of time-temperature-superposition the model enables the prediction of long-term changes along with the quantification of uncertainty stemming from sample-to-sample variation and measurement noise. All aging parameters appear to possess the same Arrhenius activation barrier, which suggests a single dominant aging mechanism at the molecular/network level.