Hiroki
Nada
*a,
Takeshi
Sakamoto
b,
Masahiro
Henmi
c,
Takafumi
Ogawa
d,
Masahiro
Kimura
d and
Takashi
Kato
*b
aNational Institute of Advanced Industrial Science and Technology (AIST), 16-1 Onogawa, Tsukuba, Ibaraki 305-8569, Japan. E-mail: hiroki.nada@aist.go.jp
bDepartment of Chemistry and Biotechnology, School of Engineering, The University of Tokyo, Hongo, Bunkyo-ku, Tokyo 113-8656, Japan. E-mail: kato@chiral.t.u-tokyo.ac.jp
cTechnology Center, Toray Industries, Inc., 1-1, Nihonbashi-Muromachi 2-chome, Chuo-ku, Tokyo 103-8666, Japan
dGlobal Environment Research Laboratories, Toray Industries, Inc., Sonoyama, Otsu, Shiga 520-0842, Japan
First published on 12th December 2019
Membranes with sub-nano channels formed by self-organization of ionic liquid-crystalline (LC) compounds have great potential as water treatment membranes. In this study, the transport mechanisms of water molecules and ions in the sub-nano channels of the LC membranes are investigated by molecular dynamics simulations for NaCl and NaNO3 solutions. The simulation results suggest that there are different transport mechanisms for water molecules and ions; the transport of water molecules occurs by Brownian diffusion, whereas that of ions occurs by jump diffusion between particular sites in the sub-nano channel. A free-energy landscape of an ion in the channel is analyzed using a metadynamics method, which indicates distinct local minima at particular sites and supports the jump diffusion mechanism. The effects of the LC compounds' structural flexibility and the electrostatic interaction with the wall of the sub-nano channels on the permeability of water molecules and ions are also investigated by molecular dynamics simulations. Simulation results suggest both the structural flexibility and the electrostatic interaction, which are characteristics of the LC membranes, are important factors in determining the water treatment performance. The structural flexibility affects the permeability of the membrane to water molecules and the electrostatic interaction reduces the permeability to ions.
Water impactRecently, we developed liquid crystalline (LC) membranes with sub-nano channels formed by self-organization of thermotropic ionic LC compounds. The LC membranes have great potential as excellent water treat membranes. In this work, we elucidated the transport mechanisms of water molecules and ions in the LC membranes by molecular dynamics simulations. We also elucidated the effects of the structural flexibility and the electrostatic interaction, which are characteristics of the LC membranes, on the water treatment performance. The knowledge obtained in this work has great potential to realize the development of future LC membranes of which the performances to water treatment are significantly improved. |
Conventional membranes contain films of cross-linked polyamides or cellulose acetates as separation functional layers. However, these films contain a disordered structure and the nanopores or sub-nanopores formed in those membranes for the transportation of water molecules and ions are not uniform.3–12 Hence, materials that contain nanopores or sub-nanopores with uniform sizes, such as carbon nanotubes (CNTs)13–20 and polymerized liquid crystals,21–29 have attracted a great deal of interest as candidates for efficient water treatment membranes. It is expected that the alignment of pores, which is required for high water permeability, will be easier for polymerized liquid crystals than for CNTs.
Previously, we reported the development of liquid-crystalline (LC) membranes with sub-nano channels formed by self-organization of thermotropic ionic LC compounds.23,25 These membranes showed high performance for water treatment, including excellent water permeability, salt rejection, and ion selectivity. We also investigated the mobility of water molecules and ions in the channel of the LC membranes by molecular dynamics (MD) simulations for a NaCl solution.25 The simulation results qualitatively explained the experimental results for the permeability of water molecules and ions. Like our LC membranes, some LC membranes developed by other groups have also exhibited selectivity for ions and molecules.22,28,29
The performance of LC membranes can be improved further, for example, by changing the chemical structure of the LC compounds. To develop better LC membranes, it is important to resolve factors determining the water treatment performance. It is essential to understand the transport mechanisms for water molecules and ions in sub-nano channels of LC membranes. Because it is difficult to analyze the transport mechanisms experimentally, MD simulations must be used. Although previous MD simulations have suggested that the mobility of water molecules and ions in sub-nano channels depends on the channel size,25 more detailed studies are needed to elucidate the factors that determine the water treatment performance.
The structures of the ionic LC compounds that form the sub-nano channels are flexible and are involved in strong electrostatic interactions with water molecules. These features of LC membranes are not seen in hard materials that have been considered for use as water treatment materials, such as zeolite30–32 and CNTs.13–20 Elucidating how each of these features affects the permeability of water molecules and ions in sub-nano channels will contribute greatly to elucidating the factors that determine the water treatment performance.
In this study, the transport mechanisms of water molecules and ions in sub-nano channels of LC membranes were examined by using MD simulations of NaCl and NaNO3 solutions. We examined a NaNO3 solution because in addition to the practical importance of removing Na+ and Cl− ions from natural water in desalination, removing NO3− ions is also important for environmental and human health reasons.33–35 Furthermore, our understanding of the transport mechanisms can be improved by examining different solutions for which the water treatment performance has been investigated experimentally.23 In addition to the transport mechanisms, the effects of the structural flexibility of LC compounds and the electrostatic interactions on the mobility of water molecules and ions were also investigated by performing MD simulations for both solutions.
Fig. 1 (a) Chemical structures of 1 and 2. (b) Unit system for system 2. Periodic boundary conditions were imposed in all x, y, and z directions of the unit system. Only the interatomic bonds (light blue columns) and the N atom (magenta spheres) of the LC monomers and the counter Cl− ions (chrome yellow spheres) are shown. The length of the simulation system in the z direction (3.6 nm) was twice that in unit system (1.8 nm). (c) The top view and side view of unit system for systems 1 and 2. The distance between the center of the system and the position of the N atom in the x–y plane was 0.67 nm for system 1 and 0.45 nm for system 2. (a) is reproduced from Fig. 1 of our previous paper.25 (b) is reproduced from Fig. S13 in the supporting information of our previous paper.25 (c) is reproduced from Fig. 5 of our previous paper.25 |
Following the previous MD simulation study,25 the simulation system including a sub-nano channel formed by the LC compounds was constructed using a simplified LC compound model from which the alkoxy chains of the nonionic moiety were eliminated (hereafter, the model is referred to as the LC monomer). For each LC compound, a rectangular parallelepiped system in which a sub-nano channel was created with 16 LC monomers and 16 counter Cl− ions (unit system, Fig. 1b) was prepared for a grand canonical Monte Carlo (GCMC) simulation to determine the number of water molecules in the channel. Then, using the final configuration generated by the GCMC simulation, the system for the MD simulations (simulation system) was constructed by putting a copy of the unit system filled with water molecules onto the original system in the z-axis direction for each LC compound. In the simulation systems, two water molecules in the channel were replaced with single Na+ and Cl− or NO3− ions. Hereafter, the simulation systems for 1 and 2 are referred to as system 1 and system 2, respectively.
The interactions for each LC compound were estimated using the potential model created in the previous MD simulation study.25 The interactions for water molecules were estimated using the TIP3P model,36 those for Na+ and Cl− ions were estimated using a model proposed by Joung and Cheatham,37 and those for NO3− ions were estimated using a model proposed by Banerjee et al.38
A standard MD simulation for the NVT ensemble at 298 K was performed for each simulation system. To investigate the effects of the LC monomer's structural flexibility and the electrostatic interaction with the LC monomers on the mobility of water molecules and ions in the channel, the MD simulation for each simulation system was performed under the following conditions: (1) the atoms constituting the LC monomers were free to move according to the force field, except for the benzene ring atoms, which were fixed at their equilibrium positions; (2) all atoms in the LC monomers were fixed at their equilibrium positions; and (3) all atoms in the LC monomers and the counter Cl− ions were fixed at their equilibrium positions, and the charges of the LC monomers' atoms and the counter Cl− ions were removed. Hereafter, conditions (1)–(3) are referred to as free, fixed, and chargeless conditions, respectively. Notably, the purpose of investigating the chargeless condition was qualitative understanding of the effect of the electrostatic interaction with the channel wall on the transport properties of water molecules and ions, but not reproducing the transport properties of water molecules and ions in the case that the charges were removed from real LC membranes. MD simulations of a NaCl solution for the free conditions were performed in our previous study.25 In the present study, MD simulations of a NaCl solution for the fixed and chargeless conditions, and MD simulations of a NaNO3 solution for all three conditions were performed. A free-energy landscape for the Na+ ion in the channel was also analyzed with a metadynamics (MTD) method39–41 to understand the transport mechanism of ions in the channel. Details of the simulation method are given in ESI.†
Fig. 2 〈dr2〉 as a function of t for the center-of-mass of water molecules, the Na+ ion, and the center-of-mass of NO3− ions in NaNO3 solution for systems 1 and 2 under the free conditions. |
The lower permeability for the NO3− ion than for water molecules originated mainly from the stronger interaction with the wall of the channel for the NO3− ion than for water molecules. We confirmed it by comparing the potential energy, U, between the LC monomer and the NO3− ion (−48.77 kJ mol−1) and U between the LC monomer and a water molecule (the averaged value over all water molecules in the channel was −0.34 kJ mol−1). During the simulation, the NO3− ion was placed near the ammonium moiety of the LC monomer, which was positively charged. Therefore, the permeability of the NO3− ion might depend on the structure of the ammonium moiety.
Fig. 3 shows the number density profiles, ρ, for the O atom of the water molecule (Ow), the Na+ ion, the N atom in the NO3− ion (Nn), and the N atom in the LC monomer (N1) along the z-axis direction for system 1 under the free conditions. Each ρ was created using the data for all three simulations. For comparison, ρ for the Na+ and Cl− ions in the NaCl solution in the channel, which were analyzed using the previous MD simulation data,25 were also shown. ρ for both solutions indicate that water molecules were distributed over the channel. However, ions were localized at particular sites in the channel. This difference in ρ between water molecules and ions suggests that the transport mechanism is different between water molecules and ions; the transport of water molecules occurs by Brownian diffusion, whereas that of ions occurs by jump diffusion between particular sites.42 An example of the jump diffusion mechanism for Na+ observed during the simulation is shown in Fig. S3 (ESI†).
In a bulk solution, each ion forms a hydration structure with the surrounding water molecules, and thus the Brownian diffusion of ions is accompanied by that of water molecules. In the channel, however, ions were stably located at particular positions, which were determined by the interaction with water molecules and the channel wall (i.e., with the LC monomers and the counter Cl− ions). Therefore, the ions in the channel did not undergo Brownian diffusion, even when the surrounding water molecules did. The positions at which ions were stably located in the channel were confirmed by creating a free-energy landscape of the Na+ ion using the MTD method (Fig. 4).
Fig. 5 〈dr2〉 as a function of t for the NaCl and NaNO3 solutions in system 1 under the free (solid lines) and fixed (open circles) conditions. |
Fig. 6 ρ of the N1 atom, Na+ ion, Ow atom, Cl− ion, and Nn atom along the z-axis direction for system 1 under the fixed conditions. The origin of the z-axis is the center of the system. |
However, the mobilities of ions for the NaCl solution for the free and fixed conditions were similar. For the NaNO3 solution, the mobility of the NO3− ion was slightly lower for the fixed conditions than for the free conditions, whereas the mobility of the Na+ ion was apparently higher for the fixed conditions than for the free conditions. For both solutions, the jump diffusion of ions occurred. The occurrence of the jump diffusion irrespective of the mobility of the surrounding water molecules suggests that the sites for the stable localization of ions in the channel were mainly determined by the interaction of ions with the channel wall. The substantial lowering of the mobility of water molecules just by fixing the structure of the LC monomers suggests that the LC monomer's structural flexibility is an important factor in increasing the permeability of water molecules in the channel. Future MD simulations for a sub-nano channel created with a realistic model of the LC compound without fixing the position of the benzene ring will provide detailed quantitative information on the relationship between the LC compound's structural flexibility and the water molecule permeability. Strictly, since the mobility of the Na+ ion for the NaNO3 solution increased by fixing the structure of the LC monomers, the LC monomer's structural flexibility would also influence somehow the permeability of ions. Performing a much larger number of the simulations with different initial positions of ions in the channel than in this study is needed to evaluate the degree of the influence.
Fig. 7 〈dr2〉 as a function of t for NaCl and NaNO3 solutions in system 1 for the chargeless condition. |
Fig. 8 ρ of the N1 atom, Na+ ion, Ow atom, Cl− ion, and Nn atom along the z-axis direction for system 1 for the chargeless conditions. The origin of the z-axis is the center of the system. |
The low mobility of water molecules and ions in the channel under the chargeless conditions does not mean that the permeability of water molecules and ions is low. Because the interaction between the channel wall and each of the water molecules and ions under the chargeless conditions is weak, the interface between the channel wall and the solution is slippery as compared with the case that the electrostatic interaction acts on the water molecules and ions. Therefore, the transport of water molecules and ions along the channel should occur readily under external pressure. We confirmed it by performing a nonequilibrium MD simulation for the NaCl solution of system 1 under the fixed and chargeless conditions (Fig. S4†). The electrostatic interaction with the channel wall, which is required for the stable localization of ions at particular sites in the channel, is an important factor in reducing the permeability of ions. Analogously, the importance of the electrostatic interaction with the wall of CNTs and boron nitride nanotubes for the dynamics of water molecules in those nanotubes was discussed by Wei and Luo.20 The structure of the channel wall, which was geometrically not flat, might also influence the permeability of water molecules. Detailed studies about the effect of the channel wall structure on the permeability of water molecules and ions should be done in the future.
Fig. 9 g OwOw in the channel for the free, fixed, and chargeless conditions as functions of the Ow–Ow distance, rOwOw. For reference, gOwOw for the bulk NaCl and NaNO3 solutions obtained by separate MD simulations are also shown. Details of the MD simulations for the bulk solutions are given in ESI† (S1-2). |
Salt | Condition | U ww (kJ mol−1) |
---|---|---|
NaCl | Free | −16.78 |
Fixed | −8.09 | |
Chargeless | −23.23 | |
NaNO3 | Free | −17.05 |
Fixed | −9.08 | |
Chargeless | −24.45 |
The formation of a highly ordered water structure in the channel for the chargeless conditions is analogous to the formation of an ice-like water structure in a restricted sub-nano space confined by hydrophobic walls, such as in a CNT.43,44Uww was much higher for the fixed conditions than for the free conditions (Table 1) because the arrangement of water molecules in the channel under the fixed conditions was dominated by the interaction with the channel wall. rOwOw at which the first peak of gOwOw appeared was longer for the fixed conditions than for the free and chargeless conditions for the same reason. Water molecules in the channels of polymer membranes are present as free and bound water.4 Our results suggest that the structural flexibility of the LC monomers is a factor in increasing the ratio of free water in the channels.
The transport mechanisms of water molecules and ions in the channel were analyzed for the NaCl and NaNO3 solutions. The results suggested that for both solutions, the water molecule and ion transport mechanisms were different. The transport of water molecules occurred by Brownian diffusion, whereas that of ions occurred by jump diffusion among sites where the ions were stably localized. Jump diffusion of ions occurred because the ions could only be stably located at particular positions in the channel, which were mainly determined by the interactions with the channel wall. The existence of particular positions at which ions were stably located in the channel was confirmed by analyzing the free-energy landscape for the Na+ ion in the channel using the MTD method.
The effects of the LC monomer's structural flexibility and the electrostatic interaction with the channel wall on the mobility of water molecules and ions were also investigated by performing MD simulations for fixed and chargeless conditions. The channel size, structural flexibility, and electrostatic interactions were important factors in determining the permeability of water molecules and ions in the channel. The simulations suggested that the structural flexibility increases the permeability of water molecules, and the electrostatic interactions decrease the permeability of ions.
The transport mechanisms of water molecules and ions, and the effects of the LC monomer's structural flexibility and the electrostatic interactions on the permeability of water molecules and ions in the channel clarified in this study will contribute to the development of superior water treatment LC membranes.
Footnote |
† Electronic supplementary information (ESI) available: Details of the simulation method, MD simulations of the bulk solutions, an example of the jump diffusion mechanism for Na+ observed during the simulation, and the mobility of water molecules in a nonequilibrium MD simulation. See DOI: 10.1039/c9ew00842j |
This journal is © The Royal Society of Chemistry 2020 |