Molecular dynamics simulations in membrane-based water treatment processes: A systematic overview

A Comprehensive Analysis of the Impact of Inorganic Matter on Membrane Organic Fouling: A Mini Review

Membrane fouling is a non-negligible issue affecting the performance of membrane systems. Particularly, organic fouling is the most persistent and severe form of fouling. The complexation between inorganic and organic matter may exacerbate membrane organic fouling. This mini review systematically analyzes the role of inorganic matter in membrane organic fouling. Inorganic substances, such as metal ions and silica, can interact with organic foulants like humic acids, polysaccharides, and proteins through ionic bonding, hydrogen bonding, coordination, and van der Waals interactions. These interactions facilitate the formation of larger aggregates that exacerbate fouling, especially for reverse osmosis membranes. Molecular simulations using molecular dynamics (MD) and density functional theory (DFT) provide valuable mechanistic insights complementing fouling experiments. Polysaccharide fouling is mainly governed by transparent exopolymer particle (TEP) formations induced by inorganic ion bridging. Inorganic coagulants like aluminum and iron salts mitigate fouling for ultrafiltration but not reverse osmosis membranes. This review summarizes the effects of critical inorganic constituents on fouling by major organic foulants, providing an important reference for membrane fouling modeling and fouling control strategies.

Thermodynamics and kinetic analysis of membrane: Challenges and perspectives

The ultimate structure of the membrane is determined using two important effects: (i) thermodynamic effect and (ii) kinetic effect. Controlling the mechanism of kinetic and thermodynamic processes in phase separation is essential for enhancing membrane performance. However, the relationship between system parameters and the ultimate membrane morphology is still largely empirical. This review focuses on the fundamental ideas behind thermally induced phase separation (TIPS) and nonsolvent-induced phase separation (NIPS) methods, including both kinetic and thermodynamic elements. The thermodynamic approach to understanding phase separation and the effect of different interaction parameters on membrane morphology has been discussed in detail. Furthermore, this review explores the capabilities and limitations of different macroscopic transport models used for the last four decades to explore the phase inversion process. The application of molecular simulations and phase field to understand phase separation has also been briefly examined. Finally, it discusses the thermodynamic approach to understanding phase separation and the consequence of different interaction parameters on membrane morphology, as well as possible directions for artificial intelligence to fill the gaps in the literature. This review aims to provide comprehensive knowledge and motivation for future modeling work for membrane fabrication via new techniques such as nonsolvent-TIPS, complex-TIPS, non-solvent assisted TIPS, combined NIPS-TIPS method, and mixed solvent phase separation.

Adsorption of Diclofenac and Its UV Phototransformation Products in an Aqueous Solution on PVDF: A Molecular Modeling Study

The presence of pharmaceuticals in drinking water has generated considerable scientific interest in potential improvements to polymeric membranes for water purification at the nanoscale. In this work, we investigate the adsorption of diclofenac and its ultraviolet (UV) phototransformation products on amorphous and crystalline poly(vinylidene difluoride) (PVDF) membrane surfaces at the nanoscale using molecular modeling. We report binding affinities by determining the free energy landscape via the extended adaptive biasing force method. The high binding affinities of the phototransformation products found are consistent with qualitative experimental results. For diclofenac, we found similar or better affinities than those for the phototransformation products, which seems to be in contrast to the experimental findings. This discrepancy can only be explained if the maximum adsorption density of diclofenac is much lower than that of the products. Overall, negligible differences between the adsorption affinities of the crystalline phases are observed, suggesting that no tuning of the PVDF surfaces is necessary to optimize filtration capabilities.

A systematic review of the molecular simulation of hybrid membranes for performance enhancements and contaminant removals

Number of research on molecular simulation and design has emerged recently but there is currently a lack of review to present these studies in an organized manner to highlight the advances and feasibility. This paper aims to review the development, structural, physical properties and separation performance of hybrid membranes using molecular simulation approach. The hybrid membranes under review include ionic liquid membrane, mixed matrix membrane, and functionalized hybrid membrane for understanding of the transport mechanism of molecules through the different structures. The understanding of molecular interactions, and alteration of pore sizes and transport channels at atomistic level post incorporation of different components in hybrid membranes posing impact to the selective transport of desired molecules are also covered. Incorporation of molecular simulation of hybrid membrane in related fields such as carbon dioxide (CO₂) removal, wastewater treatment, and desalination are also reviewed. Despite the limitations of current molecular simulation methodologies, i.e., not being able to simulate the membrane operation at the actual macroscale in processing plants, it is still able to demonstrate promising results in capturing molecule behaviours of penetrants and membranes at full atomic details with acceptable separation performance accuracy. From the review, it was found that the best performing ionic liquid membrane, mixed matrix membrane and functionalized hybrid membrane can enhance the performance of pristine membrane by 4 folds, 2.9 folds and 3.3 folds, respectively. The future prospects of molecular simulation in hybrid membranes are also presented. This review could provide understanding to the current advancement of molecular simulation approach in hybrid membranes separation. This could also provide a guideline to apply molecular simulation in the related sectors.

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.

Graphene oxide membranes with a confined mass transfer effect for Li+/Mg2+ separation: a molecular dynamics study

Molecular dynamics simulations were used to investigate the influence of the confined mass transfer effect on the separation of Mg2+ and Li+ from graphene oxide membranes, both in terms of layer spacing and degree of oxidation.

Efficient Removal of Heavy Metals from Aqueous Solutions through Functionalized γ-Graphyne-1 Membranes under External Uniform Electric Fields: Insights from Molecular Dynamics Simulations

Carbon-based nanosheet membranes with functionalized pores have great potential as water treatment membranes. In this study, the separation of Hg²⁺ and Cu²⁺ as heavy metal ions from aqueous solutions using a functionalized γ-graphyne-1 nanosheet membrane is investigated by molecular dynamics simulations. The simulation systems consist of a γ-graphyne-1 nanosheet with -COOH or -NH₂ functional groups on the edge of pores placed in an aqueous solution containing CuCl₂ and HgCl₂. An external electric field is applied as a driving force across the membrane for the separation of heavy metal ions using these functionalized pores. The ion-membrane and water molecule-membrane interaction energies, the radial distribution function of cations, the retention time and permeation of ions through the membrane, the density profile of water and ions, and the hydrogen bond in the system are investigated, and these results reveal that the performance of -NH₂-functionalized γ-graphyne-1 is better than that of -COOH-functionalized γ-graphyne-1 in the separation of Cu²⁺, while the Hg²⁺ cations encounter a high energy barrier as they pass through the membrane, especially in the -COOH-functionalized pore, due to their larger ionic radius and the smaller pore size of this membrane.

Advanced Materials for Energy-Water Systems: The Central Role of Water/Solid Interfaces in Adsorption, Reactivity, and Transport

The structure, chemistry, and charge of interfaces between materials and aqueous fluids play a central role in determining properties and performance of numerous water systems. Sensors, membranes, sorbents, and heterogeneous catalysts almost uniformly rely on specific interactions between their surfaces and components dissolved or suspended in the water-and often the water molecules themselves-to detect and mitigate contaminants. Deleterious processes in these systems such as fouling, scaling (inorganic deposits), and corrosion are also governed by interfacial phenomena. Despite the importance of these interfaces, much remains to be learned about their multiscale interactions. Developing a deeper understanding of the molecular- and mesoscale phenomena at water/solid interfaces will be essential to driving innovation to address grand challenges in supplying sufficient fit-for-purpose water in the future. In this Review, we examine the current state of knowledge surrounding adsorption, reactivity, and transport in several key classes of water/solid interfaces, drawing on a synergistic combination of theory, simulation, and experiments, and provide an outlook for prioritizing strategic research directions.

Molecular dynamics simulation for membrane separation and porous materials: A current state of art review

Computational frameworks have been under specific attention within the last two decades. Molecular Dynamics (MD) simulations, identical to the other computational approaches, try to address the unknown question, lighten the dark areas of unanswered questions, to achieve probable explanations and solutions. Owing to their complex microporous structure on one side and the intricate biochemical nature of various materials used in the structure, separative membrane materials possess peculiar degrees of complications. More notably, as nanocomposite materials are often integrated into separative membranes, thin-film nanocomposites and porous separative nanocomposite materials could possess an additional level of complexity with regard to the nanoscale interactions brought to the structure. This critical review intends to cover the recent methods used to assess membranes and membrane materials. Incorporation of MD in membrane technology-related fields such as desalination, fuel cell-based energy production, blood purification through hemodialysis, etc., were briefly covered. Accordingly, this review could be used to understand the current extent of MD applications for separative membranes. The review could also be used as a guideline to use the proper MD implementation within the related fields.

A case study of poly (aryl ether sulfone) hemodialysis membrane interactions with human blood: Molecular dynamics simulation and experimental analyses

Patients with end-stage renal diseases (ESRD) require specific health cares as the accumulation of toxins due to the lack of kidney functionality would affect their lives. However, the mortality rate is still high due to cardiovascular diseases, socks, etc. A majority of patients with chronic kidney disease (CKD) require hemodialysis services. Blood purifying membranes, as the main component of hemodialysis setups, however, still suffer from lack of optimum biocompatibility, which results in morbidity and mortality of hemodialysis service receiving patients. The goal of the present case study is to have an in-depth understanding of the current blood-hemodialysis membrane interactions occurring during hemodialysis sessions using poly (aryl ether sulfone)-poly (vinyl pyrrolidone) (PAES-PVP) membrane. Attenuated total reflectance-Fourier transmission infrared (ATR-FTIR) spectroscopy, Raman spectroscopy, and solid-state nuclear magnetic resonance (SSNMR) spectroscopy were used to assess the initial chemical structure of the PAES-PVP membrane along with the variations after with the infections with human blood. Furthermore, scanning electron microscopy (SEM) and Transition electron microscopy (TEM) were used to visualize the structural variation of the membrane, blood aggregations, and blood clots on the membrane surface. Besides, Molecular dynamics (MD) simulation was used to assess the interaction of PAES-PVP with major human blood proteins, in terms of interaction energy, which is a novel contribution to the area. The macromolecules (human serum albumin (HSA), human serum transferrin (TRF), and human fibrinogen (HFG)) were chosen from the plasma protein component. These protein structures were chosen based on their different molecular size. Three advanced spectroscopy techniques and two advanced visualization techniques were used for the assessment of the membranes. Spectroscopy studies revealed amine related peak displacement and intensity shifts as indices for attachment of biological species to the polymeric membrane surfaces. Raman peaks around 370, 798, and 1299 cm^-1, which experienced significant shifts that were related to carbon-nitrogen and sulfur-oxygen bonds due to protein adhesion. Visualization techniques illustrated blood protein fouling patterns and extracellular vesicles' presence in the pore structures into membranes. The findings highlight the importance of whole structure biocompatibility improvement, rather than only focusing on surface modifications of hemodialysis membranes. Molecular dynamics simulation assessment showed various interaction behaviors for different proteins suggesting molecular weight and active residues of the protein macromolecules play an important role in interacting with polymeric structure. FB had the highest interaction (4,274,749.07 kcal/mol) and binding (10,370.90 kcal/mol) energy with the PAES-PVP structure. TRF owned the lowest interaction energy with respect to its lower molecular weight and fewer active residue count.

Optical Analysis of the Internal Void Structure in Polymer Membranes for Gas Separation

Global warming by greenhouse gas emissions is one of the main threats of our modern society, and efficient CO2 capture processes are needed to solve this problem. Membrane separation processes have been identified among the most promising technologies for CO2 capture, and these require the development of highly efficient membrane materials which, in turn, requires detailed understanding of their operation mechanism. In the last decades, molecular modeling studies have become an extremely powerful tool to understand and anticipate the gas transport properties of polymeric membranes. This work presents a study on the correlation of the structural features of different membrane materials, analyzed by means of molecular dynamics simulation, and their gas diffusivity/selectivity. We propose a simplified method to determine the void size distribution via an automatic image recognition tool, along with a consolidated Connolly probe sensing of space, without the need of demanding computational procedures. Based on a picture of the void shape and width, automatic image recognition tests the dimensions of the void elements, reducing them to ellipses. Comparison of the minor axis of the obtained ellipses with the diameters of the gases yields a qualitative estimation of non-accessible paths in the geometrical arrangement of polymeric chains. A second tool, the Connolly probe sensing of space, gives more details on the complexity of voids. The combination of the two proposed tools can be used for a qualitative and rapid screening of material models and for an estimation of the trend in their diffusivity selectivity. The main differences in the structural features of three different classes of polymers are investigated in this work (glassy polymers, superglassy perfluoropolymers and high free volume polymers of intrinsic microporosity), and the results show how the proposed computationally less demanding analysis can be linked with their selectivities.

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.

Enantioseparation processes and mechanisms in functionalized graphene membranes: Facilitated or retarded transport?

Up to date, functionalized graphene-based membranes have exhibited a promising potential in the enantioseparation. However, since precisely controlling the interlayer distance of two-dimensional materials is a great challenge in practical experiments, the transport mechanism of chiral guests in such membranes, together with various critical parameters that play a controlling role in the transport behaviors of the preferentially binding enantiomer in narrow channels, remains to be explored. The molecular dynamics (MD) simulation, especially using the steered MD (SMD) method, might be an alternative way to investigate the enantioseparation processes and mechanisms of layered membranes with different interlayer distances. In this work, D-alanine modified graphene sheets with different interlayer distances were built as membrane models, whereas D- and L-phenylalanine were selected as chiral probes. The effect of the interlayer distance and the applied external force on the enantioseparation performance was examined. Results show that such two parameters exert a significant influence on the enantioseparation performance: (a) Increasing the interlayer distance would result in a conversion from the retarded to the facilitated mechanism at a proper external force (medium); (b) both the large and small driving forces would only lead to the appearance of the retarded transport for the preferential enantiomer, unlike the moderate force; (c) the interaction energy of L-phenylalanine with D-isomer selector decreases with the rising interlayer distances studied in this work, regardless of what the external force is. Our findings can provide guidance on the practical applications in the membrane-based chiral separation.

Molecular Insights into the Composition-Structure-Property Relationships of Polyamide Thin Films for Reverse Osmosis Desalination

A molecular-level understanding of the structure-property relationship of polyamide (PA) active layers in thin-film-composite membranes remains unclear. We developed an approach to build and hydrate the PA layer in molecular dynamics simulations and reproduced realistic membrane properties, which enabled us to examine the composition-structure-permeability relationships at the molecular level. We discovered the variation of pore size distributions in the dry PA structures at different monomer compositions, leading to different water cluster distributions and wetting properties of hydrated PA films. Membrane swelling was linearly dependent on the degree of cross-linking (DC), and higher water flux was obtained in the more swelling-prone PA films because of the transition in water transport mechanisms. Continuum-like and jumping transport both occurred in PA films with smaller DC, where visible and more persistent channels existed. In the denser films, water molecules relied only on the on-and-off channels to jump from one cavity to another; however, jumping transport was more pronounced even in the less dense PA films, and all the PA structures exhibited oscillations, which provided evidence for the solution-diffusion model rather than the pore-flow model. The results not only contribute to fundamental understanding but also provide insights into the molecule-level design for next-generation membranes.

Modeling and Simulation Studies Analyzing the Pressure-Retarded Osmosis (PRO) and PRO-Hybridized Processes

Pressure-retarded osmosis (PRO) is viewed as a highly promising renewable energy process that generates energy without carbon emissions in the age of the climate change regime. While many experimental studies have contributed to the quest for an efficiency that would make the PRO process commercially viable, computational modeling and simulation studies have played crucial roles in investigating the efficiency of PRO, particularly the concept of hybridizing the PRO process with reverse osmosis (RO). It is crucial for researchers to understand the implications of the simulation and modeling works in order to promote the further development of PRO. To that end, the authors collected many relevant papers and reorganized their important methodologies and results. This review, first of all, presents the mathematical derivation of the fundamental modeling theories regarding PRO including water flux and concentration polarization equations. After that, those theories and thermodynamic theories are then applied to depict the limitations of a stand-alone PRO process and the effectiveness of an RO-PRO hybridized process. Lastly, the review diagnoses the challenges facing PRO-basis processes which are insufficiently resolved by conventional engineering approaches and, in response, presents alternative modeling and simulation approaches as well as novel technologies.

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.

Adsorption of Bovine Serum Albumin on Poly(vinylidene fluoride) Surfaces in the Presence of Ions: A Molecular Dynamics Simulation

Adsorption of bovine serum albumin (BSA) on poly(vinylidene fluoride) (PVDF) surfaces in an aqueous environment was investigated in the presence and absence of excess ions using molecular dynamics simulations. The adsorption process involved diffusion of protein to the surface and dehydration of surface-protein interactions, followed by adsorption and denaturation. Although adsorption of BSA on PVDF surface was observed in the absence of excess ions, denaturation of BSA was not observed during the simulation (1 μs). Basic and acidic amino acids of BSA were found to be directly interacting with PVDF surface. Simulation in a 0.1 M NaCl solution showed delayed adsorption of BSA on PVDF surfaces in the presence of excess ions, with BSA not observed in close proximity to PVDF surface within 700 ns. Adsorption of Cl^- on PVDF surface increased its negative charge, which repelled negatively charged BSA, thereby delaying the adsorption process. These results will be helpful for understanding membrane fouling phenomena in polymeric membranes, and fundamental advancements in these areas will lead to a new generation of membrane materials with improved antifouling properties and reduced energy demands.

Seawater desalination using pillared graphene as a novel nano-membrane in reverse osmosis process: nonequilibrium MD simulation study

Herein, the applicability and efficiency of two types of pillared graphene nanostructures, namely, (6,6)@G and (7,7)@G, were investigated as membranes in reverse osmosis seawater desalination using extensive nonequilibrium molecular dynamics simulations.

Understanding the temperature effect on transport dynamics and structures in polyamide reverse osmosis system via molecular dynamics simulations

Molecular simulations could disclose the transport dynamics, membrane structures and temperature effect on reverse osmosis process.

Water desalination using graphene nanopores: influence of the water models used in simulations

Water desalination using graphene nanopores was studied using different water models. The water permeation was found to be influenced by the bulk transport properties and the hydrogen-bond dynamics of the simulated water.

Development of Nanostructured Water Treatment Membranes Based on Thermotropic Liquid Crystals: Molecular Design of Sub-Nanoporous Materials

Supply of safe fresh water is currently one of the most important global issues. Membranes technologies are essential to treat water efficiently with low costs and energy consumption. Here, the development of self-organized nanostructured water treatment membranes based on ionic liquid crystals composed of ammonium, imidazolium, and pyridinium moieties is reported. Membranes with preserved 1D or 3D self-organized sub-nanopores are obtained by photopolymerization of ionic columnar or bicontinuous cubic liquid crystals. These membranes show salt rejection ability, ion selectivity, and excellent water permeability. The relationships between the structures and the transport properties of water molecules and ionic solutes in the sub-nanopores in the membranes are examined by molecular dynamics simulations. The results suggest that the volume of vacant space in the nanochannel greatly affects the water and ion permeability.

The Carbon-Water Interface: Modeling Challenges and Opportunities for the Water-Energy Nexus

Providing clean water and sufficient affordable energy to all without compromising the environment is a key priority in the scientific community. Many recent studies have focused on carbon-based devices in the hope of addressing this grand challenge, justifying and motivating detailed studies of water in contact with carbonaceous materials. Such studies are becoming increasingly important because of the miniaturization of newly proposed devices, with ubiquitous nanopores, large surface-to-volume ratio, and many, perhaps most of the water molecules in contact with a carbon-based surface. In this brief review, we discuss some recent advances obtained via simulations and experiments in the development of carbon-based materials for applications in water desalination. We suggest possible ways forward, with particular emphasis on the synergistic combination of experiments and simulations, with simulations now sometimes offering sufficient accuracy to provide fundamental insights. We also point the interested reader to recent works that complement our short summary on the state of the art of this important and fascinating field.

Molecular dynamic simulations of pressure-driven water transport through polyamide nanofiltration membranes at different membrane densities

The amorphous aromatic polyamide membranes with different membrane densities were modeled to study the porous structure of free-volume pores and the pressure-driven water transport by using molecular dynamics simulations.