Transport diffusion of liquid water and methanol through membranes

Simulation study on the structural and dynamic properties of ethanol confined in nanochannels

The structural and dynamic properties of ethanol molecules in nanochannels of various diameters have been studied using molecular dynamics simulations.

Molecular Simulation and Analysis of Sorption Process toward Theoretical Prediction for Liquid Permeation through Membranes

The need to understand and describe permeation through membranes has driven the development of many well-established transport models. The modeling parameters such as solubility, diffusivity, and permeability represent the intrinsic nature of molecular interactions between membrane and permeants. In this study, we report a simulation and analysis methodology for liquid permeation. On the basis of a single simulation of liquid sorption process into a membrane, the solubility and diffusivity are estimated simultaneously; then, the permeability is predicted by the solution-diffusion model. The methodology is applied to water permeation through two representative membranes: a polymer of intrinsic microporosity (PIM-1) and a zeolitic imidazolate framework (ZIF-96). For amorphous PIM-1 membrane, the predicted water permeability agrees perfectly with simulation. For crystalline ZIF-96 membrane, water permeability is fairly well predicted. Furthermore, water dynamics in the membranes is analyzed by simulation trajectories and water structure is characterized by hydrogen bonds. Together with these microscopic insights, this study provides a simple theoretical approach to quantitatively describe water sorption, diffusion, and permeation, and it can be further applied to other liquid permeation (e.g., organic solvent nanofiltration).

Structure and energetics of hydrogen-bonded networks of methanol on close packed transition metal surfaces

Methanol is a versatile chemical feedstock, fuel source, and energy storage material. Many reactions involving methanol are catalyzed by transition metal surfaces, on which hydrogen-bonded methanol overlayers form. As with water, the structure of these overlayers is expected to depend on a delicate balance of hydrogen bonding and adsorbate-substrate bonding. In contrast to water, however, relatively little is known about the structures methanol overlayers form and how these vary from one substrate to another. To address this issue, herein we analyze the hydrogen bonded networks that methanol forms as a function of coverage on three catalytically important surfaces, Au(111), Cu(111), and Pt(111), using a combination of scanning tunneling microscopy and density functional theory. We investigate the effect of intermolecular interactions, surface coverage, and adsorption energies on molecular assembly and compare the results to more widely studied water networks on the same surfaces. Two main factors are shown to direct the structure of methanol on the surfaces studied: the surface coverage and the competition between the methanol-methanol and methanol-surface interactions. Additionally, we report a new chiral form of buckled hexamer formed by surface bound methanol that maximizes the interactions between methanol monomers by sacrificing interactions with the surface. These results serve as a direct comparison of interaction strength, assembly, and chirality of methanol networks on Au(111), Cu(111), and Pt(111) which are catalytically relevant for methanol oxidation, steam reforming, and direct methanol fuel cells.

Porous aromatic frameworks impregnated with fullerenes for enhanced methanol/water separation

Molecular simulation techniques have revealed that the incorporation of fullerenes within porous aromatic frameworks (PAFs) remarkably enhances methanol uptake while inhibiting water uptake. The highest selectivity of methanol over water is found to be 1540 at low pressure (1 kPa) and decreases gradually with increasing pressure. The adsorption of water is very small compared to methanol, a useful material property for membrane and adsorbent-based separations. Grand canonical Monte Carlo (GCMC) simulations are utilized to calculate the pure component and mixture adsorption isotherms. The water and methanol mixture simulations show that water uptake is further inhibited above the pure component results because of the dominant methanol adsorption. Molecular dynamics (MD) simulations confirm that water diffusivity is also inhibited by strong methanol adsorption in the mixture. Overall, this study reveals profound hydrophobicity in C60@PAF materials and recommends C60@PAFs as suitable applicants for adsorbent and membrane-based separations of methanol/water mixtures and other alcohol/water separation applications.

Molecular dynamics simulation of a pressure-driven liquid transport process in a cylindrical nanopore using two self-adjusting plates

Fluid transport through a nanopore in a membrane was investigated by using a novel molecular dynamics approach proposed in this study. The advantages of this method, relative to dual-control-volume grand-canonical molecular dynamics method, are that it eliminates disruptions to the system dynamics that are normally created by inserting or deleting particles from control volumes, and that it functions well for dense systems due to the number of particles being fixed in the system. Using the proposed method, we examined liquid argon transport through a nanopore by performing nonequilibrium molecular dynamics (NEMD) simulations under different back pressures. Validation of the code was performed by comparing simulation results to published experimental data obtained under equilibrium conditions. NEMD results show that constant pressure difference across the membrane was readily achieved.

Self-consistent molecular dynamics formulation for electric-field-mediated electrolyte transport through nanochannels

A self-consistent molecular dynamics (SCMD) formulation is presented for electric-field-mediated transport of water and ions through a nanochannel connected to reservoirs or baths. The SCMD formulation is compared with a uniform field MD approach, where the applied electric field is assumed to be uniform, for 2nm and 3.5nm wide nanochannels immersed in a 0.5M KCl solution. Reservoir ionic concentrations are maintained using the dual-control-volume grand canonical molecular dynamics technique. Simulation results with varying channel height indicate that the SCMD approach calculates the electrostatic potential in the simulation domain more accurately compared to the uniform field approach, with the deviation in results increasing with the channel height. The translocation times and ionic fluxes predicted by uniform field MD can be substantially different from those predicted by the SCMD approach. Our results also indicate that during a 2ns simulation time K+ ions can permeate through a 1nm channel when the applied electric field is computed self-consistently, while the permeation is not observed when the electric field is assumed to be uniform.

Comparative study between continuum and atomistic approaches of liquid flow through a finite length cylindrical nanopore

Steady state pressure driven flow of liquid argon through a finite length cylindrical nanopore was investigated numerically by classical Navier-Stokes (NS) hydrodynamic models and nonequilibrium molecular dynamics (MD) simulations. In both approaches, the nanopore was nominally 2.2 nm in diameter and 6 nm long. For the MD simulations, the intermolecular properties of the walls were specified independently from the liquid. Comparisons between the approaches were made in terms of the gross feature of total flow rate through the nanopore, as well as the more refined considerations of the spatial distributions of pressure, density, and velocity. The results showed that for the NS equations to predict the same trends in total flow rate with increasing pressure difference as the MD simulation, submodels for variations in density and viscosity with pressure are needed to be included. The classical NS boundary conditions quantitatively agreed with the flow rate predictions from MD simulations only under the condition of having a neutral-like solid-liquid interaction. Under these conditions, the NS and MD models also agreed well in streamwise distributions of pressure, density, and velocity, but not in the radial direction.

Transport of a liquid water and methanol mixture through carbon nanotubes under a chemical potential gradient

In this work, we report a dual-control-volume grand canonical molecular dynamics simulation study of the transport of a water and methanol mixture under a fixed concentration gradient through nanotubes of various diameters and surface chemistries. Methanol and water are selected as fluid molecules since water represents a strongly polar molecule while methanol is intermediate between nonpolar and strongly polar molecules. Carboxyl acid (-COOH) groups are anchored onto the inner wall of a carbon nanotube to alter the hydrophobic surface into a hydrophilic one. Results show that the transport of the mixture through hydrophilic tubes is faster than through hydrophobic nanotubes although the diffusion of the mixture is slower inside hydrophilic than hydrophobic pores due to a hydrogen network. Thus, the transport of the liquid mixture through the nanotubes is controlled by the pore entrance effect for which hydrogen bonding plays an important role.

Predictive Assessment of Surface Resistances in Zeolite Membranes Using Atomically Detailed Models

The diffusive transport of molecules through nanoporous membranes is determined by both intracrystalline diffusion and mass transport resistances associated with entering and leaving the membrane material. We compare two methods for assessing the relative importance of these resistances based on atomically detailed descriptions of the membrane material. For extremely thin membranes, net transport can be assessed using dual control volume grand canonical molecular dynamics (DCV-GCMD). We show that previous implementations of this technique may have been influenced by nonisothermal effects in interfacial regions and suggest a simple remedy to this situation. We also introduce an approximate method that uses information only from equilibrium MD simulations, which avoids the significant computational expense associated with DCV-GCMD. This approximate method can be used to rapidly assess the importance of interface-related resistances to mass transport over broad ranges of membrane operating conditions. This method will be useful in allowing a rapid determination of whether these interface resistances are significant in practical experimental situations. These two methods are compared by considering the transport of CH(4) and CF(4) through defect-free silicalite membranes.

Molecular simulation study of water interactions with oligo (ethylene glycol)-terminated alkanethiol self-assembled monolayers

Molecular simulations were performed to study a system consisting of protein (e.g., lysozyme) and self-assembled monolayers (SAMs) terminating with different chemical groups in the presence of explicit water molecules and ions. Mixed SAMs of oligo (ethylene glycol) [S(CH2)4(OCH2CH2)4OH, (OEG)] and hydroxyl-terminated SAMs [S(CH2)4OH] with a mole fraction of OEG at chiOEG = 0.2, 0.5, 0.8, and 1.0 were used in this study. In addition, methyl-terminated SAMs [S(CH2)11CH3] were also studied for comparison. The structural and dynamic behavior of hydration water, the flexibility and conformation state of SAMs, and the orientation and conformation of protein were examined. Simulation results were compared with those of experiments. It appears that there is a correlation between OEG surface resistance to protein adsorption and the surface density of OEG chains, which leads to a large number of tightly bound water molecules around OEG chains and the rapid mobility of hydrated SAM chains.

Coarse-grained nonequilibrium approach to the molecular modeling of permeation through microporous membranes

We present a modeling technique that combines a statistical-mechanical coarse-graining scheme with a nonequilibrium molecular simulation algorithm to provide an efficient simulation of steady-state permeation across a microporous material. The coarse-graining scheme is based on the mapping of an atomistic model to a lattice using multidimensional free-energy and transition-state calculations. The nonequilibrium simulation algorithm is a stochastic, lattice version of the recently developed atomistic dual-control-volume grand canonical molecular dynamics. We demonstrate the approach on a model of methane permeating through a bulk portion of siliceous zeolite ZK4 at 300 K under imposed fugacity differences. We predict the coarse-grained (cage-level) density profiles and observe the development of nonlinearities as the magnitude of the fugacity difference is increased. From the net flux of methane we also predict a mean permeability coefficient under the various conditions. The simulation results are obtained over time scales on the order of microseconds and length scales on the order of dozens of nanometers.