Aromatic Polyamide Reverse-Osmosis Membrane: An Atomistic Molecular Dynamics Simulation

Molecular Dynamics Insights into Water Transport Mechanisms in Polyamide Membranes: Influence of Cross-Linking Degree

Polyamide (PA) membranes are widely utilized in desalination and water treatment applications, yet the mechanisms underlying water transport within these amorphous polymer materials remain insufficiently understood. To gain more insight into these problems on a microscopic scale, we employ molecular dynamics (MD) simulations to analyze the relationship between the structural properties and the water permeation behavior of PA membranes. Two distinct atomistic models of PA membranes are developed by controlling their degrees of cross-linking (DC). We then conducted a comparative analysis on their microscopic structural properties and configurations of water inside the membranes and investigated how these differences lead to different water diffusion coefficients. Our results reveal that the membrane with a lower DC exhibits higher polymer mobility and a more orderly microscopic structure, allowing the formation of pores that can hold larger water clusters as well as more transient passages between pores, both contributing to an increased water diffusion coefficient. From these observations, we can conclude that water permeability within PA membranes is governed by both the morphology of semirigid pores and the oscillatory movements of the polymer chains. Overall, these findings contribute to a deeper understanding of the intricate mechanisms governing water permeation in PA membranes and may inform the design of more efficient membranes for reverse osmosis and other water treatment technologies.

Nonbonded Molecular Interaction Controls Aggregation Kinetics of Hydrophobic Molecules in Water

Molecular aggregation frequently occurs during material synthesis, cellular processes, and drug delivery systems, often resulting in decreased performance and efficiency. One major reason for such aggregation in an aqueous solution is hydrophobicity. While the basic understanding of the aggregation process of hydrophobic molecules from a thermodynamic standpoint is known, the present literature lacks a connection between the aggregation kinetics and the molecular basis of hydrophobicity. This study explores how various fluorescent probes (rhodamine dyes) aggregate in an aqueous solution due to their hydrophobicity. The method employs a combination of modeling and characterization to comprehend the aggregation process by examining the nonbonded intermolecular interactions. The aggregation kinetics was analyzed by measuring the average diffusivity of the molecules using fluorescent correlation spectroscopy and NMR diffusion measurements. Through all-atom molecular dynamics (MD) simulations, it has been observed that the level of hydrophobicity is strongly correlated to the total number of hydrogen bonds between water molecules and dyes. In addition, the aggregation frequency of colliding species, which depends on the concentration, is inversely related to hydrogen bonding and the diffusivity of the molecules. This study of small molecules was applied to predict protein aggregation rates, demonstrating strong alignment with the existing literature. The study has also helped to identify and understand the concentration at which a hydrophobic molecule does not aggregate in an aqueous solution. The method developed here could help investigate the aggregation process and its root causes at the molecular level in aqueous systems to develop strategies to control it.

Machine Learning-Aided Inverse Design and Discovery of Novel Polymeric Materials for Membrane Separation

Polymeric membranes have been widely used for liquid and gas separation in various industrial applications over the past few decades because of their exceptional versatility and high tunability. Traditional trial-and-error methods for material synthesis are inadequate to meet the growing demands for high-performance membranes. Machine learning (ML) has demonstrated huge potential to accelerate design and discovery of membrane materials. In this review, we cover strengths and weaknesses of the traditional methods, followed by a discussion on the emergence of ML for developing advanced polymeric membranes. We describe methodologies for data collection, data preparation, the commonly used ML models, and the explainable artificial intelligence (XAI) tools implemented in membrane research. Furthermore, we explain the experimental and computational validation steps to verify the results provided by these ML models. Subsequently, we showcase successful case studies of polymeric membranes and emphasize inverse design methodology within a ML-driven structured framework. Finally, we conclude by highlighting the recent progress, challenges, and future research directions to advance ML research for next generation polymeric membranes. With this review, we aim to provide a comprehensive guideline to researchers, scientists, and engineers assisting in the implementation of ML to membrane research and to accelerate the membrane design and material discovery process.

Mix-Charged Nanofiltration Membrane for Efficient Organic Removal from High-Salinity Wastewater: The Role of Charge Spatial Distribution

The efficient removal of organic contaminants from high-salinity wastewater is crucial for resource recovery and achieving zero discharge. Nanofiltration (NF) membranes are effective in separating organic compounds and monovalent salts, but they typically exhibit an excessive rejection of divalent salts. Modifying the charge characteristics of NF membranes can improve salt permeation; however, the role of charge spatial distribution in governing salt transport behavior is not fully understood. In this study, we developed a mix-charged NF membrane with a horizontal charge distribution by employing interfacial polymerization combined with a polyester template etching and solvent-induced polyamine intercalation strategy. The ratio of positive to negative charge domains in the membrane can be precisely controlled by adjusting the aqueous monomer ratio and polyamine modifier type. X-ray photoelectron spectroscopy (XPS) depth profiling and separation layer thickness analysis confirmed the complete penetration of polyamines into the separation layer, providing direct evidence of the formation of horizontally distributed charge domains. This unique charge distribution results in a high charge density and a near-electroneutral surface, which facilitates the permeation of the divalent salts. The size-dependent "plug-in" modification and covalent cross-linking further reduce pore size, enhancing rejection of small organic molecules. Additionally, the membrane demonstrated exceptional antifouling performance against both negatively and positively charged pollutants, attributed to its unique charge distribution and smooth surface. Molecular dynamics (MD) simulations further revealed that weak electrostatic interactions and a tightly bound hydration layer contribute to the membrane's superior antifouling properties. This work provides valuable insights into the design of NF membranes with tailored microstructures and charge distributions for improved water treatment performance.

Interfacial Polymerization of Aromatic Polyamide Reverse Osmosis Membranes

Polyamide membranes are widely used in reverse osmosis (RO) water treatment, yet the mechanism of interfacial polymerization during membrane formation is not fully understood. In this work, we perform atomistic molecular dynamics simulations to explore the cross-linking of trimesoyl chloride (TMC) and m-phenylenediamine (MPD) monomers at the aqueous-organic interface. Our studies show that the solution interface provides a function of "concentration and dispersion" of monomers for cross-linking. The process starts with rapid cross-linking, followed by slower kinetics. Initially, amphiphilic MPD monomers diffuse in water and accumulate at the solution interface to interact with TMC monomers from the organic phase. As cross-linking progresses, a precross-linked thin film forms, reducing monomer diffusion and reaction rates. However, the structural flexibility of the amphiphilic film, influenced by interfacial fluctuations and mixed interactions with water and the organic solvent at the solution interface, promotes further cross-linking. The solubility of MPD and TMC monomers in different organic solvents (cyclohexane versus n-hexane) affects the cross-linking rate and surface homogeneity, leading to slight variations in the structure and size distribution of subnanopores. Our study of the interfacial polymerization process in explicit solvents is essential for understanding membrane formation in various solvents, which will be crucial for optimal polyamide membrane design.

Nonpreferential Solvent Transport through an Intrinsic Cyclodextrin Pore in a Polyester Film

We performed equilibrium molecular dynamics simulations to study the transport of water and hexane solvents through cyclodextrin(CD)-based membranes (α-/β-/γ-CD/TMC). Although it is known that water and hexane can permeate through the macrocyclic cavity, surprisingly, when it is present in the CD-based membrane (α-/β-/γ-CD/TMC), these solvents are not permeating through the CD cavity. Interactions between membrane functional group atoms with the water and hexane suggest that these solvents primarily permeate through the polar aggregate pores formed via ester-linkage rather than the CD cavity. Our observation reveals that both solvents can permeate through the membrane; however, the hexane flux was one order of magnitude lower than water flux. Our study suggests that further work is needed to confirm the functional significance of the macrocyclic cavity in solvent permeation and the existence of Janus pathways.

Molecular Design of Functional Polymers for Silica Scale Inhibition

Silica polymerization, which involves the condensation reaction of silicic acid, is a fundamental process with wide-ranging implications in biological systems, material synthesis, and scale formation. The formation of a silica-based scale poses significant technological challenges to energy-efficient operations in various industrial processes, including heat exchangers and water treatment membranes. Despite the common strategy of applying functional polymers for inhibiting silica polymerization, the underlying mechanisms of inhibition remain elusive. In this study, we synthesized a series of nitrogen-containing polymers as silica inhibitors and elucidated the role of their molecular structures in stabilizing silicic acids. Polymers with both charged amine and uncharged amide groups in their backbones exhibit superior inhibition performance, retaining up to 430 ppm of reactive silica intact for 8 h under neutral pH conditions. In contrast, monomers of these amine/amide-containing polymers as well as polymers containing only amine or amide functionalities present insignificant inhibition. Molecular dynamics simulations reveal strong binding between the deprotonated silicic acid and a polymer when the amine groups in the polymer are protonated. Notably, an extended chain conformation of the polymer is crucial to prevent proximity between the interacting monomeric silica species, thereby facilitating effective silica inhibition. Furthermore, the hydrophobic nature of alkyl segments in polymer chains disrupts the hydration shell around the polymer, resulting in enhanced binding with ionized silicic acid precursors compared to monomers. Our findings provide novel mechanistic insights into the stabilization of silicic acids with functional polymers, highlighting the molecular design principles of effective inhibitors for silica polymerization.

Mechanisms and models for water transport in reverse osmosis membranes: history, critical assessment, and recent developments

Water scarcity is one of the greatest societal challenges facing humanity. Reverse osmosis (RO) desalination, a widely used membrane-based technology, has proven to be effective to augment water supply in water-stressed regions of our planet. However, progress in the design and development of RO membranes has been limited. To significantly enhance the performance of RO membranes, it is essential to acquire a deep understanding of the membrane separation and transport mechanisms. In this tutorial review, we cover the pivotal historical developments in RO technology, examine the chemical and physical properties of RO membrane materials, and critically review the models and mechanisms proposed for water transport in RO membranes. Based on recent experimental and computational findings, we conduct a thorough analysis of the key transport models-the solution-diffusion and pore-flow models-to assess their validity for accurately describing water transport in RO membranes. Our analysis involves examining the experimental evidence in favor of the solution-diffusion mechanism. Specifically, we explain whether the water content gradient within the membrane, cited as evidence for the key assumption in the solution-diffusion model, can drive a diffusive transport through RO membranes. Additionally, we review the recent molecular dynamics simulations which support the pore-flow mechanism for describing water transport in RO membranes. We conclude by providing future research directions aimed at addressing key knowledge gaps in water transport phenomena in RO membranes, with the goal of advancing the development of next-generation RO membranes.

Quantitatively Unveiling the Structure-Activity Relationship of Polyamide Nanofiltration Membranes with Complex Structures

Nanofiltration polyamide (NF PA) membranes are widely used in seawater desalination and wastewater treatment due to their excellent permeability. The structure-activity relationship of PA membranes has attracted extensive attention in decades. In this work, NF PA membranes with planar structure, nodular structure, and peak-valley structure were constructed, and the pure water permeance was calculated by nonequilibrium molecular dynamics simulation to quantitatively investigate the structure-activity relationship between the microstructure and water permeance. Results showed that the peak-valley structure had the highest effective utilization rate of the membrane surface, which had the highest number of water molecules that passed through membranes per unit cross-sectional area (7.09). Furthermore, with the increase of the specific surface area ratio, the water permeance of the NF PA with peak-valley increased at a rate about 2.5 times than that of the planar NF PA. Therefore, some molecular scale insights were supplied about the structure-activity relationship of NF PA membranes, which is helpful for the fabrication of high-performance NF PA membranes.

Semiaromatic Polyamide-Based Membrane in Forward Osmosis: Molecular Insights

Despite the increased interest in forward osmosis (FO) in recent years, the technology's advancement in commercial and industrial applications has been hampered by the absence of suitable FO membranes and ideal draw solutes, which demands the exploration of new membranes and novel draw solutes targeted for some specific applications. In this context, we considered a semiaromatic polyamide (SAPA) for an application where monovalent salt can be permeated but has high selectivity toward divalent salt and excellent water permeability. In this regard, we constructed an atomistic model for the membrane via a heuristic approach using an equilibrated mixture of hydrolyzed trimesoyl chloride and piperazine monomers and performed nonequilibrium molecular dynamics simulations on the SAPA membrane in the FO process to understand the structural properties and performance of the membrane at the atomistic level. We used pure water as the feed and Na₂SO₄ as the draw solution. It is observed that the SAPA membrane shows excellent water permeability and no reverse draw solute flux. To further test the dynamics of salt ions inside the membranes, we performed two distinct equilibrium simulations on systems consisting of either monovalent salt, such as NaCl, or divalent salt, such as Na₂SO₄. The atomistic details of the interactions between the functional groups of the membrane and salt ions provided in this work can inspire further experiments on SAPA membranes in the context of separation of monovalent and divalent salts, which have applications in the treatment of textile industry wastewater.

Hydration behaviors of nonfouling zwitterionic materials

Zwitterionic materials have emerged as highly effective ultralow fouling materials for many applications, however the underlying mechanism of fouling resistance remains unclear. Using ab initio molecular dynamics simulations and surface-sensitive sum frequency generation vibrational spectroscopy, we studied the hydration behaviors of zwitterionic materials, including trimethylamine-N-oxide (TMAO) and carboxybetaines of different charge-separation distances, to understand their fouling-resistant mechanism and provide a design principle for improved performance. Our study reveals that the interplay among hydrogen bonding, net charge, and dipole moment is crucial to the fouling-resistant capabilities of zwitterionic materials. Shortening of the zwitterionic spacing strengthens hydrogen bonding with water against biomolecule attachment due to the increased electrostatic and induction interactions, charge transfer, and improved structural stability. Moreover, the shortened charge separation reduces the dipole moment of zwitterionic materials with an intrinsic near-neutral net charge, decreasing their electrostatic and dipole-dipole interactions with biofoulers, and increasing their resistance to fouling. Compared to carboxybetaine compounds, TMAO has the shortest zwitterionic spacing and exhibits the strongest hydrogen bonding, the smallest net charge, and the minimum dipole moment, making it an excellent nonfouling material.

Correlating the Diffusion of Water to Performance in Model Reverse Osmosis Polyamides with Controlled Crosslink Densities

We systematically reduce the cross-link density of a PA network based on m-phenylene diamine by substituting a fraction of the trifunctional trimesoyl chloride cross-linking agent with a difunctional isophthaloyl analog that promotes chain extension, in order to elucidate robust design cues for improving the polyamide (PA) separation layer in reverse osmosis (RO) membranes for desalination. Thin films of these model PA networks are fully integrated into a composite membrane and evaluated in terms of their water flux and salt rejection. By incorporating 15 mol % of the difunctional chain extender, we reduce the cross-link density of the network by a factor of two, which leads to an 80 % increase in the free or unreacted amine content. The resulting swelling of the PA network in liquid water increases by a factor of two accompanied by a 30 % increase in the salt passage through the membrane. Surprisingly, this leads to a 30 % decrease in the overall permeance of water through the membrane. This conundrum is resolved by quantifying the microscopic diffusion coefficient of water inside the PA network with quasi-elastic neutron scattering. In the highest and lowest cross-link density networks, water shows strong signatures of confined diffusion. At short length scales, the water exhibits a translational diffusion that is consistent with the jump-diffusion mechanism. This translational diffusion coefficient is approximately five times slower in the lowest cross-linked density network, consistent with the reduced water permeance. This is interpreted as water molecules interacting more strongly with the increased free amine content. Over longer length scales the water diffusion is confined, exhibiting mobility that is independent of length scale. The length scales of confinement from the quasi-elastic neutron scattering experiments at which this transition from confined to translational diffusion occurs is on the order of (5 to 6) Å , consistent with complementary X-ray scattering, small angle neutron scattering, and positron annihilation lifetime spectroscopy measurements. The confinement appears to come from heterogeneities in the average inter-atomic distances, suggesting that diffusion occurs by water bouncing between chains and occasionally sticking to the polar functional groups. The results obtained here are compared with similar studies of water diffusion through both rigid porous silicates and ion exchange membranes, revealing robust design cues for engineering high-performance RO membranes.

High-Flux Nanofibrous Composite Reverse Osmosis Membrane Containing Interfacial Water Channels for Desalination

A nanofibrous composite reverse osmosis (RO) membrane with a polyamide barrier layer containing interfacial water channels was fabricated on an electrospun nanofibrous substrate via an interfacial polymerization process. The RO membrane was employed for desalination of brackish water and exhibited enhanced permeation flux as well as rejection ratio. Nanocellulose was prepared by sequential oxidations of 2,2,6,6-tetramethylpiperidine-1-oxyl (TEMPO) and sodium periodate systems and surface grafting with different alkyl groups including octyl, decanyl, dodecanyl, tetradecanyl, cetyl, and octadecanyl groups. The chemical structure of the modified nanocellulose was verified subsequently by Fourier transform infrared (FTIR), thermal gravimetric analysis (TGA), and solid NMR measurements. Two monomers, trimesoyl chloride (TMC) and m-phenylenediamine (MPD), were employed to prepare a cross-linked polyamide matrix, i.e., the barrier layer of the RO membrane, which integrated with the alkyl groups-grafted nanocellulose to build up interfacial water channels via interfacial polymerization. The top and cross-sectional morphologies of the composite barrier layer were observed by means of scanning electron microscopy (SEM), atomic force microscopy (AFM), and transmission electron microscopy (TEM) to verify the integration structure of the nanofibrous composite containing water channels. The aggregation and distribution of water molecules in the nanofibrous composite RO membrane verified the existence of water channels, demonstrated by molecular dynamics (MD) simulations. The desalination performance of the nanofibrous composite RO membrane was conducted and compared with that of commercially available RO membranes in the processing of brackish water, where 3 times higher permeation flux and 99.1% rejection ratio against NaCl were accomplished. This indicated that the engineering of interfacial water channels in the barrier layer could substantially increase the permeation flux of the nanofibrous composite membrane while retaining the high rejection ratio as well, i.e., to break through the trade-off between permeation flux and rejection ratio. Antifouling properties, chlorine resistance, and long-term desalination performance were also demonstrated to evaluate the potential applications of the nanofibrous composite RO membrane; remarkable durability and robustness were achieved in addition to 3 times higher permeation flux and a higher rejection ratio against commercial RO membranes in brackish water desalination.

When self-assembly meets interfacial polymerization

Interfacial polymerization (IP) and self-assembly are two thermodynamically different processes involving an interface in their systems. When the two systems are incorporated, the interface will exhibit extraordinary characteristics and generate structural and morphological transformation. In this work, an ultrapermeable polyamide (PA) reverse osmosis (RO) membrane with crumpled surface morphology and enlarged free volume was fabricated via IP reaction with the introduction of self-assembled surfactant micellar system. The mechanisms of the formation of crumpled nanostructures were elucidated via multiscale simulations. The electrostatic interactions among m-phenylenediamine (MPD) molecules, surfactant monolayer and micelles, lead to disruption of the monolayer at the interface, which in turn shapes the initial pattern formation of the PA layer. The interfacial instability brought about by these molecular interactions promotes the formation of crumpled PA layer with larger effective surface area, facilitating the enhanced water transport. This work provides valuable insights into the mechanisms of the IP process and is fundamental for exploring high-performance desalination membranes.

Water transport in reverse osmosis membranes is governed by pore flow, not a solution-diffusion mechanism

We performed nonequilibrium molecular dynamics (NEMD) simulations and solvent permeation experiments to unravel the mechanism of water transport in reverse osmosis (RO) membranes. The NEMD simulations reveal that water transport is driven by a pressure gradient within the membranes, not by a water concentration gradient, in marked contrast to the classic solution-diffusion model. We further show that water molecules travel as clusters through a network of pores that are transiently connected. Permeation experiments with water and organic solvents using polyamide and cellulose triacetate RO membranes showed that solvent permeance depends on the membrane pore size, kinetic diameter of solvent molecules, and solvent viscosity. This observation is not consistent with the solution-diffusion model, where permeance depends on the solvent solubility. Motivated by these observations, we demonstrate that the solution-friction model, in which transport is driven by a pressure gradient, can describe water and solvent transport in RO membranes.

Cation-π Interactions Contribute to Hydrophobic Humic Acid Removal for the Control of Hydraulically Irreversible Membrane Fouling

Hydraulically irreversible membrane fouling is a major problem encountered during membrane-based water purification. Membrane foulants present large hydrophobic fractions, with humic acid (HA) being a prevalent example of hydrophobic natural organic matter. Furthermore, HA contains numerous aromatic rings (π electrons), and its hydrophobic interactions are a major cause of irreversible membrane fouling. To address this issue, in this study, we used the cation-π interaction, which is a strong noncovalent, competitive interaction present in water. Because the strength of cation-π interactions depends on the combination of cations and π molecules, utilizing the appropriate cations will effectively remove irreversible fouling caused by hydrophobic HA. We performed macroscale experiments to determine the cleaning potential of the test cations, nanomechanically analyzed the changes in HA cohesion caused by the test cations using a surface force apparatus and an atomic force microscope, and used molecular dynamics simulations to elucidate the HA removal mechanism of test studied cations. We found that the addition of 1-ethyl-3-methylimidazolium, an imidazolium cation with an aromatic moiety, effectively removed the HA layer by weakening its cohesion, and the size, hydrophobicity, and polarity of the HA layer synergistically affected the HA removal mechanism based on the cation-π interactions.

Ionic Liquid-Mediated Interfacial Polymerization for Fabrication of Reverse Osmosis Membranes

This study revealed the effects of incorporating ionic liquid (IL) molecules: 1-ethyl, 1-butyl, and 1-octyl-3-methyl-imidazolium chlorides with different alkyl chain lengths, in interfacial polymerization (IP) on the structure and property (i.e., permeate-flux and salt rejection ratio) relationships of resulting RO membranes. The IL additive was added in the aqueous meta-phenylene diamine (MPD; 0.1% w/v) phase, which was subsequently reacted with trimesoyl chloride (TMC; 0.004% w/v) in the hexane phase to produce polyamide (PA) barrier layer. The structure of resulting free-standing PA thin films was characterized by grazing incidence wide-angle X-rays scattering (GIWAXS), which results were correlated with the performance of thin-film composite RO membranes having PA barrier layers prepared under the same IP conditions. Additionally, the membrane surface properties were characterized by zeta potential and water contact angle measurements. It was found that the membrane prepared by the longer chain IL molecule generally showed lower salt rejection ratio and higher permeation flux, possibly due to the inclusion of IL molecules in the PA scaffold. This hypothesis was supported by the GIWAXS results, where a self-assembled surfactant-like structure formed by IL with the longest aliphatic chain length was detected.

Multimodal confined water dynamics in reverse osmosis polyamide membranes

While polyamide (PA) membranes are widespread in water purification and desalination by reverse osmosis, a molecular-level understanding of the dynamics of both confined water and polymer matrix remains elusive. Despite the dense hierarchical structure of PA membranes formed by interfacial polymerization, previous studies suggest that water diffusion remains largely unchanged with respect to bulk water. Here, we employ neutron spectroscopy to investigate PA membranes under precise hydration conditions, and a series of isotopic contrasts, to elucidate water transport and polymer relaxation, spanning ps-ns timescales, and Å-nm lengthscales. We experimentally resolve, for the first time, the multimodal diffusive nature of water in PA membranes: in addition to (slowed down) translational jump-diffusion, we observe a long-range and a localized mode, whose geometry and timescales we quantify. The PA matrix is also found to exhibit rotational relaxations commensurate with the nanoscale confinement observed in water diffusion. This comprehensive ‘diffusion map’ can anchor molecular and nanoscale simulations, and enable the predictive design of PA membranes with tuneable performance.

Polymeric membranes are extensively used in water desalination, but the effect of membrane nanostructure on water transport is still elusive. The authors, using quasi-elastic neutron scattering and contrast variation techniques, provide detailed insight into the dynamics of the polymer network and confined water across a wide range of length and timescales.

Molecular Simulations to Elucidate Transport Phenomena in Polymeric Membranes

Despite decades of dominance in separation technology, progress in the design and development of high-performance polymer-based membranes has been incremental. Recent advances in materials science and chemical synthesis provide opportunities for molecular-level design of next-generation membrane materials. Such designs necessitate a fundamental understanding of transport and separation mechanisms at the molecular scale. Molecular simulations are important tools that could lead to the development of fundamental structure-property-performance relationships for advancing membrane design. In this Perspective, we assess the application and capability of molecular simulations to understand the mechanisms of ion and water transport across polymeric membranes. Additionally, we discuss the reliability of molecular models in mimicking the structure and chemistry of nanochannels and transport pathways in polymeric membranes. We conclude by providing research directions for resolving key knowledge gaps related to transport phenomena in polymeric membranes and for the construction of structure-property-performance relationships for the design of next-generation membranes.

Interfacial synthesis of large-area ultrathin polyimine nanofilms as molecular separation membrane

Thin film membranes of covalent organic frameworks are promising for high-permeance molecular separation. However, their synthesis needs a high temperature or longer reaction time, unsuitable for large-scale fabrication of thin film composite membranes. The ultrathin film of porous organic polymers as a separation layer of the composite membrane could be a close alternative to COF membranes. Here we report transition metal ion-catalyzed room temperature fabrication of the ultrathin (≈12 nm) polyimine nanofilms via interfacial polymerization of melamine and triformylphloroglucinol onto hydrolyzed polyacrylonitrile support within a short reaction time. Composite membranes exhibit high water permeance (≈78 L m⁻² h⁻¹ bar⁻¹), high rejection (99.6%) of brilliant blue R (825.9 g mol⁻¹), low rejection of NaCl (≈1.8%) and Na₂SO₄ (≈17%), and enable efficient molecular separation. The role of metal ion catalysts for large-area fabrication of the ultrathin polyimine nanofilm membranes used for molecular separation is demonstrated.

•

Ultrathin porous polyimine nanofilms could be a close alternative to COF membranes

•

Large-area polyimine nanofilms are formed via interfacial polymerization

•

Transition metal ions favor the formation of the nanofilms at room temperature

•

Polyimine nanofilm membranes display superior permselectivity and tunable MWCO

Molecular Methods for Assessing the Morphology, Topology, and Performance of Polyamide Membranes

The molecular-scale morphology and topology of polyamide composite membranes determine the performance characteristics of these materials. However, molecular-scale simulations are computationally expensive and morphological and topological characterization of molecular structures are not well developed. Molecular dynamics simulation and analysis methods for the polymerization, hydration, and quantification of polyamide membrane structures were developed and compared to elucidate efficient approaches for producing and analyzing the polyamide structure. Polymerization simulations that omitted the reaction-phase solvent did not change the observed hydration, pore-size distribution, or water permeability, while improving the simulation efficiency. Pre-insertion of water into the aggregate pores (radius ≈ 4 Å) of dry domains enabled shorter hydration simulations and improved simulation scaling, without altering pore structure, properties, or performance. Medial axis and Minkowski functional methods were implemented to identify permeation pathways and quantify the polyamide morphology and topology, respectively. Better agreement between simulations and experimentally observed systems was accomplished by increasing the domain size rather than increasing the number of ensemble realizations of smaller systems. The largest domain hydrated was an order of magnitude larger by volume than the largest domain previously reported. This work identifies methods that can enable more efficient and meaningful fundamental modeling of membrane materials.

In Situ Characterization of Dehydration during Ion Transport in Polymeric Nanochannels

The transport of hydrated ions across nanochannels is central to biological systems and membrane-based applications, yet little is known about their hydrated structure during transport due to the absence of in situ characterization techniques. Herein, we report experimentally resolved ion dehydration during transmembrane transport using modified in situ liquid ToF-SIMS in combination with MD simulations for a mechanistic reasoning. Notably, complete dehydration was not necessary for transport to occur across membranes with sub-nanometer pores. Partial shedding of water molecules from ion solvation shells, observed as a decrease in the average hydration number, allowed the alkali-metal ions studied here (lithium, sodium, and potassium) to permeate membranes with pores smaller than their solvated size. We find that ions generally cannot hold more than two water molecules during this sterically limited transport. In nanopores larger than the size of the solvation shell, we show that ionic mobility governs the ion hydration number distribution. Viscous effects, such as interactions with carboxyl groups inside the membrane, preferentially hinder the transport of the mono- and dihydrates. Our novel technique for studying ion solvation in situ represents a significant technological leap for the nanofluidics field and may enable important advances in ion separation, biosensing, and battery applications.

Fabrication of desalination membranes by interfacial polymerization: history, current efforts, and future directions

Membrane desalination is a promising technology for addressing the global challenge of water scarcity by augmenting fresh water supply. Continuous progress in this technology relies on development of membrane materials. The state-of-the-art membranes used in a wide range of desalination applications are polyamide thin-film composite (TFC) membranes which are formed by interfacial polymerization (IP). Despite the wide use of such membranes in desalination, their real-world application is still hampered by several technical obstacles. These challenges of the TFC membranes largely stem from the inherent limitations of the polyamide chemistry, as well as the IP reaction mechanisms. In the past decade, we have witnessed substantial progress in the understanding of polyamide formation mechanisms and the development of new IP strategies that can potentially lead to the redesign of TFC membranes. In this Tutorial, we first present a brief history of the development of desalination membranes and highlight the major challenges of the existing TFC membranes. We then proceed to discuss the pros and cons of emerging IP-based fabrication strategies aiming at improving the performance of TFC membranes. Next, we present technical obstacles and recent efforts in the characterization of TFC membranes to enable fundamental understanding of relevant mechanisms. We conclude with a discussion of the current gap between industrial needs and academic research in designing high-performance TFC membranes, and provide an outlook on future research directions for advancing IP-based fabrication processes.

Effect of an ionic environment on membrane fouling: a molecular dynamics study

The effect of the ionic environment on membrane fouling was investigated for polyamide (PA) and graphene oxide (GO) membranes using equilibrium molecular dynamics (MD) simulations. For each of these membranes, bovine serum albumin (BSA) was considered as the model foulant. The effect of the foulant on the membranes is investigated at seawater concentration and also in a normal aqueous environment. We investigated the translational and rotational motion of the protein relative to the membrane, interaction energy between the protein and the membrane surface, structural changes in the protein, and ion distribution around the protein and the membrane surface for all the systems. We found that the effects of ions were very different on both the membranes. Specifically, with an increase in ionic strength, the repulsion between the protein and membrane was observed in the case of GO, while for PA, no significant changes were observed for the same. Also, the ion distribution around the protein and the membrane surface were found to be different. In particular, for GO, there were more number of chloride ions around the protein and the membrane than that of sodium ions, which was probably the reason for the repulsion in the case of GO. However, in the case of PA, the membrane surface did not exhibit any affinity towards a specific ion, and the protein in the case of PA was surrounded by more number of sodium ions than chloride ions.

Discontinuous Molecular Dynamics Simulations of Biomolecule Interfacial Behavior: Study of Ovispirin-1 Adsorption on a Graphene Surface

Fundamental understanding of biomolecular interfacial behavior, such as protein adsorption at the microscopic scale, is critical to broad applications in biomaterials, nanomedicine, and nanoparticle-based biosensing techniques. The goal of achieving both computational efficiency and accuracy presents a major challenge for simulation studies at both atomistic and molecular scales. In this work, we developed a unique, accurate, high-throughput simulation method which, by integrating discontinuous molecular dynamics (DMD) simulations with the Go-like protein-surface interaction model, not only solves the dynamics efficiently, but also describes precisely the protein intramolecular and intermolecular interactions at the atomistic scale and the protein-surface interactions at the coarse-grained scale. Using our simulation method and in-house developed software, we performed a systematic study of α-helical ovispirin-1 peptide adsorption on a graphene surface, and our study focused on the effect of surface hydrophobic interactions and π-π stacking on protein adsorption. Our DMD simulations were consistent with full-atom molecular dynamics simulations and showed that a single ovispirin-1 peptide lay down on the flat graphene surface with randomized secondary structure due to strong protein-surface interactions. Peptide aggregates were formed with an internal hydrophobic core driven by strong interactions of hydrophobic residues in the bulk environment. However, upon adsorption, the hydrophobic graphene surface can break the hydrophobic core by denaturing individual peptide structures, leading to disassembling the aggregate structure and further randomizing the ovispirin-1 peptide's secondary structures.

Unveiling the Molecular Mechanisms of Thickness-Dependent Water Dynamics in an Ultrathin Free-Standing Polyamide Membrane

Aromatic polyamide (PA) membranes fabricated from interfacial polymerization are widely used for desalination and water treatment. The fabrication of the high-flux PA membrane requires a fundamental understanding of the molecular mechanisms of water dynamics in the PA, which is still obscure due to the limited experimental methods. Herein, molecular dynamics (MD) simulations were employed to establish an atomic model of ultrathin free-standing PA membranes with various thickness and to explore the thickness-dependent dynamics of water molecules in the PA membrane. Simulation results illustrate that the simulated PA membrane has an average pore radius of 3 Å similar to the free volume size of the experimental PA membrane measured by PALS. The PA could be identified as the swelling layer (SL) and the confined layer (CL) based on their water diffusion rates. The diffusivity of water in the confined layer of PA membrane was much lower than that in the swelling layer and thus determined the water flux of the PA membrane. The water diffusivity in the sub-8 nm PA membrane is greatly enhanced due to a very thin confined layer thickness, illustrating the mechanism of the experimentally fabricated sub-8 nm PA membrane having the dramatically enhanced water permeability. Furthermore, results show that water molecules tend to transport rapidly in the free space inside the PA membrane. Our results provide some insights into the thickness-dependent water dynamics in the PA on a molecular level and may help to design the next generation of high-flux PA membranes.

Machine Learning Analysis of the Thermodynamic Responses of In Situ Dielectric Spectroscopy Data in Amino Acids and Inorganic Electrolytes

Dielectric spectroscopy (DS) can be a robust in situ technique for geochemical applications. In this study, we applied deep-learning techniques to DS measurement data to enable rapid science interrogation and identification of electrolyte solutions containing salts and amino acids over a wide temperature range (20 to -60 °C). For the purpose of searching for signs of life, detecting amino acids is a fundamental high priority for field and planetary instruments as amino acids are one of the building blocks for life as we know it. A convolutional neural network (CNN) with channel-wise one-dimensional filters is proposed to fulfill the task, using the DS data of amino acid and inorganic salt solutions. Experimental results show that the CNN with two convolutional layers and one fully connected layer can effectively differentiate solutions containing amino acids from those containing salts in both the liquid and solid (water ice) states. To complement the experimental measurements and CNN analysis, the diffusive behaviors of ions (K⁺, Cl^-, and OH^-) were further discussed with atomistic molecular dynamics simulations performed in this work as well as the quantum simulation published in the literature. Combining DS with machine-learning techniques and simulations will greatly facilitate more real-time decision-making of mobility systems for future exploratory endeavors in other worlds beyond Earth.

Desalination behavior analysis of interior-modified carbon nanotubes doped membrane by dielectric spectrum and molecular simulation

Carbon nanotube (CNT)-doped polyamide (PA) membranes have attracted much attention in reverse osmosis (RO) membranes due to their significant advantages of water flux and desalination. In this study, we synthesized multi-walled carbon nanotube (MWNT)/PA RO membrane by 12-oxodidodecanoic acid methyl ester group interior-modified MWNTs (MWNT-C₁₄H₂₅O₄). Then, their mechanism of desalination behavior was successfully analyzed by combining dielectric relaxation spectrum (DRS) and molecular dynamics (MD) simulation. DRS analysis mainly focuses on two aspects: (1) the water volume fraction, average pore size and dielectric parameters of MWNT-C₁₄H₂₅O₄/PA and PA membranes were obtained by model analysis of DRS data. These data of MWNT-C₁₄H₂₅O₄/PA membrane are higher than PA membrane, which indicates that the water flux of the MWNT-C₁₄H₂₅O₄/PA membrane was higher than that of the PA membrane. (2) Further analysis shows that the MWNT-C₁₄H₂₅O₄/PA membranes have high average charge density, ion solvation barrier and reflection coefficient, which indicates that the added interior-modified MWNT can improve the salt rejection of PA membranes. In the microscopic aspect, the desalination behavior of the MWNT-C₁₄H₂₅O₄/PA and PA membrane was analyzed from the aspects of free volume distribution, the dynamic diffusion process of water and ions. The results show that the microscopic data of dynamic simulation well support the conclusion of the DRS method. This study provides a convenient methodology to characterize the properties of the membrane from the aspect of membrane structure.

Molecular Understanding of Ion Effect on Polyzwitterion Conformation in an Aqueous Environment

Polyzwitterions (PZs) are promising materials for the antifouling in reverse osmosis and nanofiltration membrane technology for water treatment. Fundamental understanding of the structure and molecular interactions involving zwitterions is crucial to the optimal design of antifouling in membrane separation. Here we employ the umbrella sampling and molecular dynamics simulations to investigate molecular interactions between sulfobetaine/carboxybetaine zwitterions and different metal ions (Na⁺, K⁺, and Ca²⁺) in an aqueous solution. The simulation results show that these ions can form stable or metastable contact ionic/solvent-shared-ionic pairs with zwitterions. Simulations at different grafting densities of PZ brush arrays reveal complex competitive association mechanisms, which are attributed to nonbonded electrostatic and van der Waals interactions among zwitterions, water molecules, and different metal ions in an aqueous environment. While the high-grafting density of the PZ brush array leads to a strong branch association between different zwitterions in water, this association is decreased at intermediate- and low-grafting densities due to strong zwitterion-water interactions. More importantly, adding ions into water at intermediate- and low-grafting densities further breaks down the zwitterion branch association, resulting in a randomly oriented and dispersed branch configuration with significant swelling of the polymers. The degree of swelling depends on the type of ions, which further changes the surface electrostatic potential of PZ coatings.

Atomistic Simulations of Biofouling and Molecular Transfer of a Cross-linked Aromatic Polyamide Membrane for Desalination

Reverse osmosis through a polyamide (PA) membrane is an important technique for water desalination and purification. In this study, molecular dynamics simulations were performed to study the biofouling mechanism (i.e., protein adsorption) and nonequilibrium steady-state water transfer of a cross-linked PA membrane. Our results demonstrated that the PA membrane surface's roughness is a key factor of surface's biofouling, as the lysozyme protein adsorbed on the surface's cavity site displays extremely low surface diffusivity, blocking water passage, and decreasing water flux. The adsorbed protein undergoes secondary structural changes, particularly in the pressure-driven flowing conditions, leading to strong protein-surface interactions. Our simulations were able to present water permeation close to the experimental conditions with a pressure difference as low as 5 MPa, while all the electrolytes, which are tightly surrounded by hydration water, were effectively rejected at the membrane surfaces. The analysis of the self-intermediate scattering function demonstrates that the dynamics of water molecules coordinated with hydrogen bonds is faster inside the pores than during the translation across the pores. The pressure difference applied shows a negligible effect on the water structure and content inside the membrane but facilitates the transportation of hydrogen-bonded water molecules through the membrane's sub-nanopores with a reduced coordination number. The linear relationship between the water flux and the pressure difference demonstrates the applicability of continuum hydrodynamic principles and thus the stability of the membrane structure.

Influence of the presence of cations on the water and salt dynamics inside layered graphene oxide (GO) membranes

Although over the past few years, graphene oxide (GO) has emerged as a promising membrane material, the applicability of layered GO membranes in water purification/seawater desalination is still a challenging issue because of the undesirable swelling of GO laminates in the aqueous environment. One of the ways to tune the interlayer spacing and to arrest the undesirable swelling of layered GO membranes in the aqueous environment is to intercalate the interlayer spacing of the GO laminates with cations. Although the cation intercalation imparts stabilization to GO laminates in the aqueous environment, their effect on the performance of the membrane is yet to be addressed in detail. In the present study we have investigated the effect of cation intercalation on the performance of layered GO membranes using molecular dynamics simulation. For the same interlayer spacing, the cation intercalated layered GO membranes have a higher water flux as compared to the corresponding pristine layered GO membranes. In the presence of the cations, the water molecules inside the interlayer gallery get more compactly packed. The presence of the cations also increases the stability of the hydrogen bond network among the water molecules inside the membrane. This can be attributed to slow water reorientation dynamics inside the interlayer gallery in the presence of the cations. The synergistic effect of all these changes is that the water permeability through the cation intercalated layered GO membranes is higher as compared to that through the corresponding pristine layered GO membranes. On the other hand, the intercalation of the cations (K+, Mg2+) leads to higher rejection of Na+ ions whereas the rejection of Cl- ions slightly decreases.

Data for molecular recognition between polyamide thin film composite on the polymeric subtract by molecular dynamic

This paper focus to examine the best molecular interaction between Polyamide Thin Film Composite (PA TFC) layers with different properties of the support membrane. The support membrane of Nylon 66 (N66) and Polyvinylidene fluoride (PVDF) was chosen to represent the hydrophilic and hydrophobic model respectively in the Molecular Dynamic (MD) simulation. The Condensed-Phase Optimized Molecular Potential for Atomistic Simulation Studies (COMPASS) force field was used with the total simulation runs were set 1000 picoseconds run production ensembles. The temperature and pressure set for both ensembles were 298 K and 1 atm respectively. The validity of our model densities data was check and calculated where the deviation must be less than 6%. The comparison between hydrophobic and hydrophilic of the support membrane data was examined by the distance and magnitude of intensity of the Radial Distribution Function (RDF's) trends.

Molecular Structure of Aromatic Reverse Osmosis Polyamide Barrier Layers

The molecular structures of polyamide barrier layers in reverse osmosis membranes, made by interfacial polymerization of m-phenylenediamine and trimesoyl chloride under different reaction and post-treatment conditions, were characterized by grazing incidence wide-angle X-ray scattering (GIWAXS). The molecular backbone packing is consistent with two different aromatic molecular packing motifs (parallel and perpendicular) with preferential surface-induced orientation. The results suggest that the perpendicular, T-shaped, packing motif (5 Å spacing) might be associated with optimal membrane permeance, compared with the parallel packings (3.5-4.0 Å spacings).

Anomalous Dynamics of Water in Polyamide Matrix

Water in polymer matrixes is likely to show anomalous dynamics, a problem that has not been well understood yet. Here, we performed atomistic molecular dynamics simulations to study the water dynamics in a polyamide (PA) matrix, the bulk phase of well-known reverse osmosis membranes. For time-dependent ensemble average, water molecules experienced ballistic diffusion at a shorter time scale, followed by a crossover from subdiffusion to Brownian diffusion at a time scale ∼10 ns, and non-Gaussian diffusion, an indication of anomalous dynamics, sticks on even in the Brownian diffusion region. The anomalous dynamics mainly originates from two distinct motions including small-step continuous diffusion and jumping diffusion. The jumping motion has a mean length of 3.08 ± 0.31 Å and characteristic relaxation time of 0.218 ± 0.040 ns, which dominates the water diffusion in a fully hydrated PA matrix. It comprised low- and high-frequency jumps; the former is almost unchanged, and the latter remarkably increases with the increase of the hydration level. Surrounding neighbors of water strongly affect the jumping frequency, which exponentially or linearly decays with the increase in the number of atoms from the PA matrix. Although the PA matrix is flexible, associated with the water dynamics, the translocation of water is mainly through either tracing the position of neighboring water or jumping into the adjacent accommodation space.

How Pore Hydrophilicity Influences Water Permeability?

Membrane separation is playing increasingly important role in providing clean water. Simulations predict that membrane pores with strong hydrophobicity produce ultrahigh water permeability as a result of low friction. However, experiments demonstrate that hydrophilic pores favor higher permeability. Herein we simulate water molecules transporting through interlayers of two-dimensional nanosheets with various hydrophilicities using nonequilibrium molecular dynamics. We reveal that there is a threshold pressure drop (ΔP _T), exceeding which stable water permeability appears. Strongly hydrophobic pores exhibit extremely high ΔP _T, prohibiting the achievement of ultrahigh water permeability under the experimentally accessible pressures. Under pressures < ΔP _T, water flows in hydrophobic pores in a running-stop mode because of alternative wetting and nonwetting, thus leading to significantly reduced permeability. We discover that hydrophilic modification to one surface of the nanosheet can remarkably reduce ΔP _T by > 99%, indicating a promising strategy to experimentally realize ultrafast membranes.

Transport mechanism of water molecules passing through polyamide/COF mixed matrix membranes

Mixed matrix membranes (MMMs) have gained significant attention due to their high water permeability without the cost of salt rejection. The mechanism of permeability promotion for PA/COFs MMMs is investigated in this work from molecular insights.

Molecular dynamics simulation studies of the structure and antifouling performance of a gradient polyamide membrane

The structure and the antifouling performance of the first gradient polyamide layer model are systematically disclosed using molecular dynamics simulations.