Virtual Synthesis of Thermally Cross-Linked Copolymers from a Novel Implementation of Polymatic

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.

Molecular Dynamics Simulation Study of Organic Solvents Confined in PIM-1 and P84 Polyimide Membranes

Organic solvent nanofiltration (OSN) has recently proved to be a promising separation process thanks to the development of membrane materials with suitable resistance toward organic solvents. Among those materials, P84 polyimide membranes are currently the most used in OSN while PIM-1 membranes have recently attracted attention due to their high permeance in apolar solvents and alcohols. Both P84 and PIM-1 membranes have nanosized free volumes, and their separation performance is finely connected to polymer/solvent interactions. Consequently, modeling OSN membranes at the molecular scale is highly desirable in order to rationalize experimental observations and gain a deeper insight into the molecular mechanisms ruling solvent and solute permeation. A prerequisite for understanding solvent transport through OSN membranes is therefore to characterize the membrane/solvent interactions at the molecular level. For that purpose, we carried out molecular simulations of three different solvents, acetone, methanol, and toluene in contact with P84 and PIM-1 membranes. The solvent uptake by both membranes was found to be correlated to the degree of confinement of the solvent, the polymer swelling ability and polymer/solvent interactions. The translational dynamics of the solvent molecules in the PIM-1 membrane was found to be correlated with the solvent viscosity due to the relatively large pores of this membrane. That was not the case with the P84 membrane, which has a much denser structure than the PIM-1 membrane and for which it was observed that the translational dynamics of the confined solvent molecules was directly correlated to the affinity between the P84 polymer and the solvent.

General-Purpose Coarse-Grained Toughened Thermoset Model for 44DDS/DGEBA/PES

The objective of this work is to predict the morphology and material properties of crosslinking polymers used in aerospace applications. We extend the open-source dybond plugin for HOOMD-Blue to implement a new coarse-grained model of reacting epoxy thermosets and use the 44DDS/DGEBA/PES system as a case study for calibration and validation. We parameterize the coarse-grained model from atomistic solubility data, calibrate reaction dynamics against experiments, and check for size-dependent artifacts. We validate model predictions by comparing glass transition temperatures measurements at arbitrary degree of cure, gel-points, and morphology predictions against experiments. We demonstrate for the first time in molecular simulations the cure-path dependence of toughened thermoset morphologies.

Routine million-particle simulations of epoxy curing with dissipative particle dynamics

Mesoscale simulation techniques have helped to bridge the length scales and time scales needed to predict the microstructures of cured epoxies, but gaps in computational cost and experimental relevance have limited their impact. In this work, we develop an open-source plugin epoxpy for HOOMD-Blue that enables epoxy crosslinking simulations of millions of particles to be routinely performed on a single modern graphics card. We demonstrate the first implementation of custom temperature-time curing profiles with dissipative particle dynamics and show that reaction kinetics depend sensitively on the stochastic bonding rates. We provide guidelines for modeling first-order reaction dynamics in a classic epoxy/hardener/toughener system and show structural sensitivity to the temperature-time profile during cure. We conclude with a discussion of how these efficient large-scale simulations can be used to evaluate ensembles of epoxy processing protocols to quantify the sensitivity of microstructure on processing.

Hypercrosslinked porous polymer materials: design, synthesis, and applications

Hypercrosslinked polymers (HCPs) are a series of permanent microporous polymer materials initially reported by Davankov, and have received an increasing level of research interest. In recent years, HCPs have experienced rapid growth due to their remarkable advantages such as diverse synthetic methods, easy functionalization, high surface area, low cost reagents and mild operating conditions. Judicious selection of monomers, appropriate length crosslinkers and optimized reaction conditions yielded a well-developed polymer framework with an adjusted porous topology. Post fabrication of the as developed network facilitates the incorporation of various chemical functionalities that may lead to interesting properties and enhance the selection toward a specific application. To date, numerous HCPs have been prepared by post-crosslinking polystyrene-based precursors, one-step self-polycondensation or external crosslinking strategies. The advent of these methodologies has prompted researchers to construct well-defined porous polymer networks with customized micromorphology and functionalities. In this review, we describe not only the basic synthetic principles and strategies of HCPs, but also the advancements in the structural and morphological study as well as the frontiers of potential applications in energy and environmental fields such as gas storage, carbon capture, removal of pollutants, molecular separation, catalysis, drug delivery, sensing etc.

NLDFT Pore Size Distribution in Amorphous Microporous Materials

The pore size distribution (PSD) is one of the most important properties when characterizing and designing materials for gas storage and separation applications. Experimentally, one of the current standards for determining microscopic PSD is using indirect molecular adsorption methods such as nonlocal density functional theory (NLDFT) and N₂ isotherms at 77 K. Because determining the PSD from NLDFT is an indirect method, the validation can be a nontrivial task for amorphous microporous materials. This is especially crucial since this method is known to produce artifacts. In this work, the accuracy of NLDFT PSD was compared against the exact geometric PSD for 11 different simulated amorphous microporous materials. The geometric surface area and micropore volumes of these materials were between 5 and 1698 m²/g and 0.039 and 0.55 cm³/g, respectively. N₂ isotherms at 77 K were constructed using Gibbs ensemble Monte Carlo (GEMC) simulations. Our results show that the discrepancies between NLDFT and geometric PSD are significant. NLDFT PSD produced several artificial gaps and peaks that were further confirmed by the coordinates of inserted particles of a specific size. We found that dominant peaks from NLDFT typically reported in the literature do not necessarily represent the truly dominant pore size within the system. The confirmation provides concrete evidence for artifacts that arise from the NLDFT method. Furthermore, a sensitivity analysis was performed to show the high dependency of PSD as a function of the regularization parameter, λ. A higher value of λ produced a broader and smoother PSD that closely resembles geometric PSD. As an alternative, a new criterion for choosing λ, called here the smooth-shift method (SSNLDFT), is proposed that tuned the NLDFT PSD to better match the true geometric PSD. Using the geometric pore size distribution as our reference, the smooth-shift method reduced the root-mean-square deviation by ∼70% when the geometric surface area of the material is greater than 100 m²/g.

Ultrathin Molecular-Layer-by-Layer Polyamide Membranes: Insights from Atomistic Molecular Simulations

In this study, we present an atomistic simulation study of several physicochemical properties of polyamide (PA) membranes formed from interfacial polymerization or from a molecular-layer-by-layer (mLbL) on a silicon wafer. These membranes are composed of meta-phenylenediamine (MPD) and benzene-1,3,5-tricarboxylic acid chloride (TMC) for potential reverse osmosis (RO) applications. The mLbL membrane generation procedure and the force field models were validated, by comparison with available experimental data, for hydrated density, membrane swelling, and pore size distributions of PA membranes formed by interfacial polymerization. Physicochemical properties such as density, free volume, thickness, the degree of cross-linking, atomic compositions, and average molecular orientation (which is relevant for the mLbL membranes) are compared for these different processes. The mLbL membranes are investigated systematically with respect to TMC monomer growth rate per substrate surface area, MPD/TMC ratio, and the number of mLbL deposition cycles. Atomistic simulations show that the mLbL deposition generates membranes with a constant film growth if both the TMC monomer growth rate and MPD/TMC monomer ratio are kept constant. The film growth rate increases with TMC monomer growth rate or MPD/TMC ratio. Furthermore, it was found on one hand that the mLbL membrane density and free volume varies significantly with respect to the TMC monomer growth rate, while on the other hand the degree of cross-linking and the atomic composition varies considerably with the MPD/TMC ratio. Additionally, it was found that both TMC and MPD orient at a tilted angle with respect to the substrate surface, where their angular distribution and average angle orientation depend on both the TMC growth rate and the number of deposition cycles. This study illustrates that molecular simulations can play a crucial role in the understanding of structural properties that can empower the design of the next generation RO membranes created from molecular-layer-by-layer (mLbL) on a silicon wafer.

Polymer nanofilms with enhanced microporosity by interfacial polymerization

Highly permeable and selective membranes are desirable for energy-efficient gas and liquid separations. Microporous organic polymers have attracted significant attention in this respect owing to their high porosity, permeability and molecular selectivity. However, it remains challenging to fabricate selective polymer membranes with controlled microporosity that are stable in solvents. Here we report a new approach to designing crosslinked, rigid polymer nanofilms with enhanced microporosity by manipulating the molecular structure. Ultrathin polyarylate nanofilms with thickness down to 20 nm are formed in situ by interfacial polymerization. Enhanced microporosity and higher interconnectivity of intermolecular network voids, as rationalized by molecular simulations, are achieved by using contorted monomers for the interfacial polymerization. Composite membranes comprising polyarylate nanofilms with enhanced microporosity fabricated in situ on crosslinked polyimide ultrafiltration membranes show outstanding separation performance in organic solvents, with up to two orders of magnitude higher solvent permeance than membranes fabricated with nanofilms made from non-contorted planar monomers.

Microscopic Model of the Metal-Organic Framework/Polymer Interface: A First Step toward Understanding the Compatibility in Mixed Matrix Membranes

An innovative computational methodology integrating density functional theory calculations and force field-based molecular dynamics simulations was developed to provide a first microscopic model of the interactions at the metal-organic framework (MOF) surface/polymer interface. This was applied to the case of the composite formed by the polymer of intrinsic microporosity, PIM-1, and the zeolitic imidazolate framework, ZIF-8, as a model system. We found that the structure of the composite at the interface is the result of both the chemical affinity between PIM-1 and ZIF-8 and the rigidity of the polymer. Specifically, there is a preferential interaction between the -CN groups of PIM-1 and the NH terminal functions of the organic linker at the ZIF-8 surface. Additionally, the resulting conformation of the polymer gives rise to interfacial microvoids at the vicinity of the MOF surface. The porosity, rigidity, and density of the interfacial polymer were analyzed and compared to those for the bulk polymer. It was shown that the polymer still feels the impact of the MOF surface even at long distances above 15-20 Å. Further, both the polydispersity of the polymer and the flexibility of the MOF surface were revealed to only slightly affect the properties of the MOF/interface. This work, which delivers a microscopic picture of the MOF surface/polymer interactions at the interface, would lead, in turn, to the understanding of the compatibility in MOF-based mixed-matrix membranes.

In Silico Determination of Gas Permeabilities by Non-Equilibrium Molecular Dynamics: CO2 and He through PIM-1

We study the permeation dynamics of helium and carbon dioxide through an atomistically detailed model of a polymer of intrinsic microporosity, PIM-1, via non-equilibrium molecular dynamics (NEMD) simulations. This work presents the first explicit molecular modeling of gas permeation through a high free-volume polymer sample, and it demonstrates how permeability and solubility can be obtained coherently from a single simulation. Solubilities in particular can be obtained to a very high degree of confidence and within experimental inaccuracies. Furthermore, the simulations make it possible to obtain very specific information on the diffusion dynamics of penetrant molecules and yield detailed maps of gas occupancy, which are akin to a digital tomographic scan of the polymer network. In addition to determining permeability and solubility directly from NEMD simulations, the results shed light on the permeation mechanism of the penetrant gases, suggesting that the relative openness of the microporous topology promotes the anomalous diffusion of penetrant gases, which entails a deviation from the pore hopping mechanism usually observed in gas diffusion in polymers.