Molecular dynamics simulation study of the mechanisms of water diffusion in a hydrated, amorphous polyamide

The physical basis for solvent flow in organic solvent nanofiltration

Organic solvent nanofiltration (OSN) is an emerging membrane technology that could revolutionize chemical separations in numerous vital industries. Despite its significance, there remains a lack of fundamental understanding of solvent transport mechanisms in OSN membranes. Here, we use an extended Flory-Rehner theory, nonequilibrium molecular dynamic simulations, and organic solvent transport experiments to demonstrate that solvent flow in OSN membranes is driven by a pressure gradient. We show that solvent molecules migrate as clusters through interconnected pathways within the membrane pore structure, challenging the widely accepted diffusion-based view of solvent transport in OSN. We further reveal that solvent permeance is dependent on solvent affinity to the OSN membrane, which, in turn, controls the membrane pore structure. Our fundamental insights lay the scientific groundwork for the development of next-generation OSN membranes.

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.

Interplay of the forces governing steroid hormone micropollutant adsorption in vertically-aligned carbon nanotube membrane nanopores

Vertically-aligned carbon nanotube (VaCNT) membranes allow water to conduct rapidly at low pressures and open up the possibility for water purification and desalination, although the ultralow viscous stress in hydrophobic and low-tortuosity nanopores prevents surface interactions with contaminants. In this experimental investigation, steroid hormone micropollutant adsorption by VaCNT membranes is quantified and explained via the interplay of the hydrodynamic drag and friction forces acting on the hormone, and the adhesive and repulsive forces between the hormone and the inner carbon nanotube wall. It is concluded that a drag force above 2.2 × 10-3 pN overcomes the friction force resulting in insignificant adsorption, whereas lowering the drag force from 2.2 × 10-3 to 4.3 × 10-4 pN increases the adsorbed mass of hormones from zero to 0.4 ng cm-2. At a low drag force of 1.6 × 10-3 pN, the adsorbed mass of four hormones is correlated with the hormone-wall adhesive (van der Waals) force. These findings explain micropollutant adsorption in nanopores via the forces acting on the micropollutant along and perpendicular to the flow, which can be exploited for selectivity.

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.

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.

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.

Finely tuned polyamide structure with green plasticizers to construct ultrafast water channels for effective desalination

Highly permeable reverse osmosis (RO) membranes are desirable for alleviating the energy burden and ensuring future water sustainability. Herein, the effectiveness of green plasticizer-assisted interfacial polymerization (GPAIP) for preparing polyamide thin-film composite (TFC) RO membranes with significantly enhanced water permeability was demonstrated. The presence of green citrate plasticizers, namely tributyl citrate (TBC) or acetyl tributyl citrate (ATBC), led to the formation of new hydrogen bonds and inhibited the formation of the initial interchain amide-amide bonding, thus markedly reducing chain rigidity as demonstrated by the decreased elasticity modulus. More flexible polyamide chains resulted in the creation of more ultrafast water channels during filtration. Furthermore, TBC-modified membranes exhibited more elastic polyamide layers and higher water flux than that of ATBC-modified membranes on account of the presence of both hydrogen bond acceptors (OH) and hydrogen bond donors (C=O) in TBC molecules. Specifically, water flux of 0.6 wt% TBC-modified and 0.6 wt% ATBC-modified membranes was 83.6 L m^-2 h^-1 and 49.7 L m^-2 h^-1 respectively, more than 5 times and 3 times that of the pristine membrane. The excellent performance of TFC RO membranes fabricated via GPAIP together with the facile membrane manufacturing process offered the possibility of breaking the predicament in desalination field, which could eventually help ease the current freshwater crisis.

Water Transport and Ion Diffusion Investigation of an Amphotericin B-Based Channel Applied to Forward Osmosis: A Simulation Study

The use of an Amphotericin B_Ergosterol (AmBEr) channel as an artificial water channel in forward osmosis filtration (FO) was studied via molecular dynamics (MD) simulation. Three channel models were constructed: a common AmBEr channel and two modified C3deOAmB_Ergosterol (C3deOAmBEr) channels with different diameters (12 Å and 18 Å). During FO filtration simulation, the osmotic pressure of salt-water was a driving force for water permeation. We examined the effect of the modified C3deOAmBEr channel on the water transport performance. By tracing the change of the number of water molecules along with simulation time in the saltwater region, the water permeability of the channel models could be calculated. A higher water permeability was observed for a modified C3deOAmBEr channel, and there was no ion permeation during the entire simulation period. The hydrated ions and water molecules were placed into the channel to explore the ion leakage behavior of the channels. The mean squared displacement (MSD) of ions and water molecules was obtained to study the ion leakage performance. The Amphotericin B-based channels showed excellent selectivity of water molecules against ions. The results obtained on an atomistic scale could assist in determining the properties and the optimal filtration applications for Amphotericin B-based channels.

Nanocomposite desalination membranes made of aromatic polyamide with cellulose nanofibers: synthesis, performance, and water diffusion study

Reverse osmosis membranes of aromatic polyamide (PA) reinforced with a crystalline cellulose nanofiber (CNF) were synthesized and their desalination performance was studied. Comparison with plain PA membranes shows that the addition of CNF reduced the matrix mobility resulting in a molecularly stiffer membrane because of the attractive forces between the surface of the CNFs and the PA matrix. Fourier transform-infrared spectroscopy and X-ray photoelectron spectroscopy results showed complex formation between the carboxy groups of the CNF surface and the m- phenylenediamine monomer in the CNF-PA composite. Molecular dynamics simulations showed that the CNF-PA had higher hydrophilicity which was key for the higher water permeability of the synthesized nanocomposite membrane. The CNF-PA reverse osmosis nanocomposite membranes also showed enhanced antifouling performance and improved chlorine resistance. Therefore, CNF shows great potential as a nanoreinforcing material towards the preparation of nanocomposite aromatic PA membranes with longer operation lifetime due to its antifouling and chlorine resistance properties.

Ion partitioning and permeation in charged low-T* membranes

Understanding ion transport in membrane materials is key to engineering and development of desalination and water purification technologies as well as electro-membrane applications. To date, modeling of ion transport has mainly relied on mean-field approaches, originally intended for weak inter-ionic interactions, i.e., high reduced temperature T*. This condition is violated in many membranes, which could explain disagreement between predicted trends and experiments. The paper highlights observed discrepancies and develops a new approach based on the concept of ion association, more adequate in the low-T^⁎ limit. The new model addresses ion binding and mobility consistently within the same physical picture, applied to different types of single and mixed salts. The resulting relations show a significantly weaker connection between ion partitioning and permeability than the standard ones. Estimates using primitive model (PM) of ions in a homogeneous dielectric suggest that non-PM mechanisms, originating from the molecular structure of the ion-solvating environment, might enhance ion association in membranes. PM analysis also predicts that ion solvation and association must be rigidly related, yet non-PM effects may decouple these phenomena and allow a crossover to non-trivial regimes consistent with experiments and simulations. Despite the crude nature of the presented approach and some questions remaining open, it appears to explain most available experimental data and presents a step towards predictive modeling of ion-selective membrane separations in water-, environment- and energy-related applications.

In silico study of structure and water dynamics in CNT/polyamide nanocomposite reverse osmosis membranes

A comprehensive understanding of the role of CNTs in the reverse osmosis process is disclosed through in silico study.

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.

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.

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.

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.

A probabilistic approach for estimating water permeability in pressure-driven membranes

A probabilistic approach is proposed to estimate water permeability in a cellulose triacetate (CTA) membrane. Water transport across the membrane is simulated in reverse osmosis mode by means of non-equilibrium molecular dynamics (MD) simulations. Different membrane configurations obtained by an annealing MD simulation are considered and simulation results are analyzed by using a hierarchical Bayesian model to obtain the permeability of the different membranes. The estimated membrane permeability is used to predict full-scale water flux by means of a process-level Monte Carlo simulation. Based on the results, the parameters of the model are observed to converge within 5-ns total simulation time. The results also indicate that the use of unique structural configurations in MD simulations is essential to capture realistic membrane properties at the molecular scale. Furthermore, the predicted full-scale water flux based on the estimated permeability is within the same order of magnitude of bench-scale experimental measurement of 1.72×10(-5) m/s.

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.

Water Desalination through Zeolitic Imidazolate Framework Membranes: Significant Role of Functional Groups

A molecular simulation study is reported for water desalination through five zeolitic imidazolate framework (ZIF) membranes, namely ZIF-25, -71, -93, -96, and -97. The five ZIFs possess identical rho-topology but differ in functional groups. The rejection of salt (NaCl) is found to be around 97% in ZIF-25, and 100% in the other four ZIFs. The permeance ranges from 27 to 710 kg/(m(2)·h·bar), about one∼two orders of magnitude higher compared with commercial reverse osmosis membranes. Due to a larger aperture size da, ZIF-25, -71, and -96 exhibit a much higher water flux than ZIF-93 and -97; however, the flux in ZIF-25, -71, and -96 is governed by the polarity of functional group rather than da. With the hydrophobic CH3 group, ZIF-25 has the highest flux despite the smallest da among ZIF-25, -71, and -96. The lifetime of hydrogen bonding in ZIF-25 is shorter than that in ZIF-71 and -96. Furthermore, water molecules undergo a fast flushing motion in ZIF-25, but frequent jumping in ZIF-96 and particularly in ZIF-97. An Arrhenius-type relationship is found between water flux in ZIF-25 and temperature, and the activation energy is predicted to be 6.5 kJ/mol. This simulation study provides a microscopic insight into water desalination in a series of ZIFs, reveals the key factors (aperture size and polarity of functional group) governing water flux, and suggests that ZIF-25 might be an interesting reverse osmosis membrane for high-performance water desalination.

Molecular Dynamics Study of Carbon Nanotubes/Polyamide Reverse Osmosis Membranes: Polymerization, Structure, and Hydration

Carbon nanotubes/polyamide (PA) nanocomposite thin films have become very attractive as reverse osmosis (RO) membranes. In this work, we used molecular dynamics to simulate the influence of single walled carbon nanotubes (SWCNTs) in the polyamide molecular structure as a model case of a carbon nanotubes/polyamide nanocomposite RO membrane. It was found that the addition of SWCNTs decreases the pore size of the composite membrane and increases the Na and Cl ion rejection. Analysis of the radial distribution function of water confined in the pores of the membranes shows that SWCNT+PA nanocomposite membranes also exhibit smaller clusters of water molecules within the membrane, thus suggesting a dense membrane structure (SWCNT+PA composite membranes were 3.9% denser than bare PA). The results provide new insights into the fabrication of novel membranes reinforced with tubular structures for enhanced desalination performance.

Water diffusion in amorphous hydrophilic systems: a stop and go process

The diffusion of H2O in three amorphous polymer-H2O systems is studied as a function of H2O content using molecular dynamics. A picture of H2O molecule motion comprising alternating steps of being bound at an adsorption site ("stop") and moving ("go") emerges. This picture is made quantitative. The bound time, frequency of stop-go steps, and tortuosity all decrease with H2O content. Fourier analysis of particle motion during bound time segments provides a measure of an attempt frequency that is connected quantitatively to the bound time and an activation energy of a hydrogen bond. For increasing H2O content, the polymer-H2O systems swell, leading to an increase in the diffusion coefficient and porosity and a decrease in activation energy.

Hydrated polyamide membrane and its interaction with alginate: a molecular dynamics study

The properties of the hydrated amorphous polyamide (PA) membrane and its binding with alginate are investigated through molecular dynamics simulations. The density of the hydrated membrane, surface morphology, and water diffusion near and inside the membrane are compared to other studies. Particular focus is given to the steered molecular dynamics (SMD) simulation of the binding between the PA membrane and an alginate model. The PA surface composition is determined on the basis of experimental measurements of the oxygen/nitrogen (O/N) ratio. The surface model is built using a configurational-bias Monte Carlo technique. The consistent valence force field (CVFF) is used to describe the atomic interactions in the membrane-foulant system. Simulation results show that the carboxylate groups in both the PA surface and alginate exhibit strong binding with metal ions. This binding mechanism plays a major role in the PA-alginate fouling through the formation of an ionic binding bridge. Specifically, Ca(2+) ions have stronger binding with the carboxylate group than Na(+) ions, while the binding breakdown time is shorter for Ca(2+) than Na(+) because of the comparably higher hydration free energy of Ca(2+) ions with water molecules.