CC BY
146
УДК 539.8, 532.5, 532.7
Selective filtration of fluids in materials with slit-shaped nanopores
A.A. Tsukanov1, E.V. Shilko1, E. Gutmanas2, and S.G. Psakhie1
1 Institute of Strength Physics and Materials Science SB RAS, Tomsk, 634055, Russia 2 Department of Materials Science and Engineering, Technion-Israel Institute of Technology, Haifa, 32000, Israel
Problems associated with a qualitative increase in the selectivity of fluid filtration remain the major challenge in a variety of areas such as fluid transport through porous materials and media, ion separation, water desalination and purification, and many others. A promising way to solve these problems is to design and develop membranes with slit-shaped nanopores. In the paper, we studied the selectivity and permeability of slit-shaped nanosized pores in the natural mineral (hydroxyapatite) with the use of the nonequilibrium molecular dynamics approach with all-atom models. We showed that the subnanometer-wide slit-shaped pores in hydroxyapatite are able to demonstrate both good salt rejection and relatively high water permeability. In particular, the numerically predicted water permeability of hydroxyapatite with 0.4 nm thick slit-shaped nanopores reaches about 200 L/(m2hbar) that is higher than that of commercial membranes and has the same order of magnitude as the theoretically predicted water permeability through single-layer MoS2 nanoporous membranes. An increase in the nanopore thickness is accompanied by a multiple growth in permeability, which is comparable with advanced 2D-CAP (2D-conjugated aromatic polymer) membranes, but in so doing the filtration selectivity is lost. The results show that nanoporous materials with the connected network of slit-shaped nanopores is a promising filter material for water treatment including seawater desalination and other important technical and environmental applications.
Keywords: slit-shaped nanopore, fluid, filtration, ion rejection, selectivity, membrane, hydroxyapatite, desalination, computer simulation, molecular dynamics
DOI 10.24411/1683-805X-2018-15002
Избирательная фильтрация жидкостей в материалах с щелевидными нанопорами
А.А. Цуканов1, Е.В. Шилько1, Э. Гутманас2, С.Г. Псахье1
1 Институт физики прочности и материаловедения СО РАН, Томск, 634055, Россия 2 Технион - Израильский технологический институт, Хайфа, 32000, Израиль
Проблемы, связанные с недостаточной селективностью при фильтрации жидкостей, в настоящее время являются серьезными вызовами в таких ключевых областях, как транспорт флюидов в пористых средах, разделение ионов, опреснение, очистка воды и многих других. Перспективным способом решения этих проблем является проектирование и разработка мембран с щелевидными нанопорами. Настоящая работа посвящена теоретическому изучению селективности и проницаемости щелевидных наноразмерных пор в природном минерале (гидроксиапатит) с использованием неравновесного молекулярно-динамического подхода и полноатомных моделей. Показано, что субнанометровые щелевидные поры в гидроксиапатите способны демонстрировать одновременно хорошую задержку ионов соли и сравнительно высокую водопроводимость. В частности, численно прогнозируемая водопрово-димость щелевидной нанопоры в гидроксиапатите с раскрытостью 0.4 нм достигает величины порядка 200л/(м2-ч-бар), что выше чем у коммерческих мембран и сравнимо по порядку величины с теоретически предсказанной проводимостью однослойной нанопористой мембраны дисульфида молибдена MoS2. Увеличение раскрытости нанопоры сопровождается многократным увеличением ее проводимости, которая становится сопоставимой с характеристиками передовых 2D-CAP мембран, однако при этом наблюдается значительная потеря в селективности фильтрации. Полученные результаты свидетельствуют, что нанопористые материалы со связанной системой щелевидных нанопор являются перспективным фильтрующим материалом для очистки воды, включая опреснение морской воды и другие важные технические и экологические применения.
Ключевые слова: щелевидная нанопора, жидкость, фильтрация, задержка ионов, селективность, мембрана, гидроксиапатит, опреснение, компьютерное моделирование, молекулярная динамика
© Tsukanov A.A., Shilko E.V., Gutmanas E., Psakhie S.G., 2018
1. Introduction
Fast climate change, rapid growth of the human population, and the permanently increasing man-caused impact on the ecological system of the Earth pose new challenges for mankind. One of the key modern problems is the need for high-quality selective filtration of fluids of various molecular composition on an industrial scale. Topical areas that require high filtration selectivity are, for example, transport of liquids through porous media [1], water desalination and purification [2], hydrogen purification [3, 4], gas separation [5], ion separation [6], DNA sequencing [7], transplantation of islets for noninvasive cure of patients with Type 1 Diabetes [8] and many other. The most important requirements for selective filtration are excellent rejection of impurity ions or organic molecules, minimal resistance to flow of the main fluid, high energy efficiency and the ability to scale the filtering nanosystem up to the scale of practical application. Progress in advanced filtration/encapsulation technologies is determined by the development of new selective/semipermeable membranes with the given geometric and electrostatic characteristics of the pore surfaces. Specific requirements for the characteristics of the pore space and pore surface of membranes are determined by the features of the field of application. In this paper, we consider some of the fundamental questions related to the effect of the pore size and shape on the quality and efficiency of ion-selective filtration of salt water.
The modern widely used membranes for purification and desalination of water are typically characterized by either very high salt rejection and comparatively low water permeability or insufficient salt rejection and very high water permeability [9]. The recent development of the nanotechnology-based filtration membranes with nanosized pores, including nanostructured ceramic membranes [10], nanoporous graphene nanosheets [11], nanotube-based membranes [12, 13], nanoporous boron nitride nanosheets [14], peptide-based membranes [15], and even biologically inspired aquaporin-based membranes [16], promises to provide the breakthrough in this area. Different aspects of the use of nanomaterials and nanostructured materials for water treatment are reviewed elsewhere [17-19].
The origin of the ion rejection in narrow nanopores is the ion size exclusion and/or solvation shell loss mechanisms. The first one takes place in the case of ions that have the atom diameter larger than dimensions of both the water molecule and nanopore entrance. The second mechanism is associated with the energy penalty for the hydrated ion entry to the nanopore, when the ion (partially or fully) loses neighboring water molecules from its solvation shell [20]. To achieve high salt rejection, pores should be of the subnanometer dimension. However, a decrease of the nano-pore size leads to a decrease in its water permeability because the water molecule can also face the energy barrier,
thus rearranging the network of hydrogen bonds with its neighbors when entering the nanopore from the bulk solution. A promising way to solve this problem is to design the nanoporous material with preferentially oriented rectangular (slit-shaped) nanopores, which are characterized by high aspect ratios. These nanopores should have a subnanometer size in one direction (thickness of slit) and a much larger size in perpendicular directions. The slit-like geometry of the nanopore provides a straight path that presents a shorter distance for molecules to travel as compared to the traditional pore geometry.
The recent study on ultrafiltration showed that membranes with slit-shaped pores are able to provide a higher efficiency and selectivity at the given value of permeability than membranes with cylindrical pores [21]. A striking example is the recently synthesized crystalline 2D-conju-gated aromatic polymer (2D-CAP) with nanopores with a rectangular cross section [22]. 2D-CAP is a very prospective base for the production of effective water desalination and ion separation membranes. An almost 100% ion rejection and high water permeability 1172 L/(m2hbar) were numerically predicted for the 2D-CAP membrane [23], which is more efficient than any commercial membrane. The limiting case of rectangular geometry of the nanopore cross section is a slit with the subnanometer width. In addition, the nanoporous material with slit-shaped pores has typically a large specific area of the internal surface, which can be potentially used as adsorbent for different species (e.g. ions, heavy metals, halides/oxyhalides, organics, etc.) depending on the material surface structure and composition.
The results of the above-mentioned studies suggest that slit-shaped nanopores would demonstrate a higher liquid permeability than that of cylindrical (nanotube-like and hole-like) nanopores of the same diameter as the slit width as well as a fairly good ion rejection due to the subnanometer size of the pore entrance (slit width). Therefore, it is important to obtain estimates of the ion rejection and water permeability of a material sample with slit-shaped nanosized pores. The present study is focused on slit-shaped nanopores in the natural biomineral, namely, in hydroxyapatite (HAP) since it is an environmentally-safe, nontoxic and biocompatible material. HAP has a pH-dependent structure and charge distribution on the surface, which further can be important for the pore selectivity control, including selective permeability and competitive adsorption properties.
The study was carried out by means of the nonequi-librium molecular dynamics (NEMD) computer simulation. Note that the current molecular level computer simulations have become a powerful and widely used tool for the study of nanoscale systems and processes, including the molecular structure and property prediction [24], nanoparticle-cell interaction [25], the study of hybrid [26, 27] or heterogeneous nanosystems [28] and nanoconfined matter [29],
ion and molecule adsorption [30, 31], and molecular transport through nanopores [1, 32]. NEMD is the well-known direct numerical way to study salt water filtration through nanopores, including the estimation of fluid fluxes, ion rejection, and ion selectivity [10, 23]. In the present paper, we employed NEMD for numerical simulation of the salt water (~800 mM) filtration through slit-shaped nanopores with the subnanometer width.
2. Computer model for estimating water permeability through slit-shaped nanopores
We developed molecular dynamics models of two systems that contain nanopores with different width:
(i) model 1 with the slit-shaped "narrow" nanopore of the ~0.42 nm width;
(ii) model 2 with the slit-shaped "wide" nanopore of the ~0.78 nm width.
Figure 1 shows equilibrium structures of the systems at 0.1 MPa and T = 300 K. The periodic boundary conditions (PBC) are applied along the y and z directions (z is the vertical axis, and axis y is directed away from the observer, see the coordinate system in Fig. 1). The {y, z} dimensions of the simulation boxes of models 1 and 2 are {3.2621, 3.2545} nm and {3.2621, 3.6161} nm, respectively. The size of the models along axis x is about 20 nm. Due to the use of PBC along the y axis, the nanopore geometry is infinite along axis y.
The length of the nanopore along axis x is 3.7668 nm in both the models. The nanopore walls are {001} surfaces of
hydroxyapatite (HAP). They are formed by protonated phosphate groups H2PO4, which corresponds to pH = 5. The hydroxyapatite model is developed based on the data derived by Lin and Heinz [33, 34]. The INTERFACE FF force field developed by Heinz et al. [35] is utilized to parameterize the models.
Model 1 contains 16916 atoms, including 3372 water molecules (SPC model [36]), 40 sodium and 40 chlorine ions (the average initial salt concentration in feed box [NaCl] is about 800 mM). The hydroxyapatite fragment comprises 2832 atoms according to the brutto-formula 48[Ca9(PO4)6Ca(OH)2Ca(H2PO4)2]. Coordinates of calcium and phosphorus atoms are "frozen" during simulation runs to avoid mineral movements and rotations under loading. The total number of atoms in model 2 is 18689, including 3819 water molecules. The number of ions and HAP atoms is the same as in model 1. Two pistons composed of six-layered fcc copper blocks with the lattice constant a = 0.3616 nm are added to models to control pressure and its gradient along the x axis. To maintain pressure conditions, the external force along the x axis is applied to the left piston Fx = S (p0 + Ap) and to the right piston F2 = -Sp0 (Fig. 1, a) where S is the piston area, p0 = = 0.1 MPa is the atmosphere pressure, and Ap is an applied pressure drop that varies from 50 to 200 MPa. All NEMD simulations are conducted at the temperature T = 300 K. Employment of the SHAKE algorithm [37] for hydrogen atoms allows us to use the time step At = 2 fs. The pairwise interaction cutoff is 1 nm with the 0.8 nm switch-
msM
WjW
'YYYYY
wyyy
0
>>>>>>#6Î
{
î >)>>>>
* KKKKJD*
— Ca of HAP
— P of HAP
— O of HAP
— H of HAP
) — Na+ ion i — Cl- ion
fig m
' m^m&ê
55555c
m
Fig. 1. Representative views of simulation boxes for the studied systems with "narrow" (a) and "wide" (b) slit-shaped nanopores in hydroxyapatite. The molecular configurations correspond to the initial state after equilibration of the models
1600
— Narrow, 50 MPa
— Narrow, 100 MPa / -Narrow, 150 MPa
Narrow, 200 MPa
/ ^
0
2000-
— Wide, 50 MPa
Wide, 100 MPa / - Wide, 150 MPa
0
Na, narrow CI, narrow •o - Na, wide - o - Cl, wide
50 100 150 200 Pressure drop, MPa
250
Fig. 2. Results of the NEMD simulations: the near-linear time dependences of the number of filtered water molecules for "narrow" (a) and "wide" (b) slit-shaped nanopores at different pressure drops from 50 to 200 MPa; average water flux through the nanopores as a function of the pressure drop (c); estimations of Na+ and Cl- ion rejection by nanopores at different pressure regimes (d)
ing distance. Long-range electrostatics is computed using the particle-particle particle-mesh (PPPM) method [38] with a relative accuracy of 10-4. All simulations are performed using the LAMMPS software of the Sandia National Laboratory (USA) [39] at the Lomonosov-2 highperformance supercomputer of the Lomonosov Moscow State University (Russia) [40, 41].
3. Results and discussion
3.1. Water permeability through slit-shaped nanopores
Series of eight NEMD simulations of the salt water filtration through two slit-shaped nanopores of different widths (~0.42 and ~0.78 nm) at the pressure differences 50, 100, 150, and 200 MPa are performed. The total time of each simulation with the narrow nanopore is 100-120 ns, but the filtration process is steady-state only during first 50100 ns because of the closed system and the limited number of water molecules. This steady-state stage can be observed by the near-linear shape of the curves in Figs. 2, a and b. Note that the duration of water flow through the wide nanopore in the steady-state regime (10-20 ns) is 5 times shorter than that through the narrow nanopore because of the about five-fold higher flow rate.
We estimate molecular water fluxes in the steady-state filtration regime for all eight cases (Fig. 2, c). One can see that the water flux for each system is a linear function of the applied hydrostatic pressure. This implies constant water permeability at least within the range of the atmospheric pressure drop up to 200 MPa. Here pure water permeability is conventionally defined as the volume of water that passes through the membrane (or the layer of permeable material) per unit time per unit area and per unit of trans-membrane/translayer pressure. The estimation of water permeability through the narrow slit-shaped nanopore is 4.7 L/(cm2dayMPa) or 196 L/(m2hbar). The obtained value is higher than that of commercial membranes and has the same order of magnitude as the numerically predicted water permeability rate through the single-layer MoS2 nano-porous membrane [42] and the water permeability of na-noporous graphene [43].
The wide slit-shaped nanopore of model 2 demonstrates a water permeability about 5 times higher than that of the nanopore of model 1:27 L/(cm2dayMPa) or 1119 L/(m2 hbar). This value is comparable with that of the 2D-CAP membrane, which is 1172 L/(m2hbar) [23]. However, as it will be seen in the next section, 2D-CAP
t = 125 ns (narrow nanopore, 50 MPa)
Fig. 3. Ion selectivity of slit-shaped nanopores: (a) two chlorine ions (shown by black arrows) passed through the narrow nanopore during 120 ns while sodium ions are completely rejected (notations are the same as in Fig. 1); (b) comparison of the total amounts of passed Na+ (black bars) and Cl- (striped bars) ions
membranes possess a much better salt rejection property than the 0.78 nm wide slit-shaped nanopore in HAP.
For a quantitative comparison, it is important to keep in mind that water permeability depends on the density (concentration) of pores. The direct comparison of pore densities is not exactly correct on account of completely different geometries of hole-like and slit-shaped nanopores. In the examined nanoporous HAP fragments the densities of slit-shaped nanopores are about 0.3 nm/nm2. This means that each 1 nm2 of the working cross section has 0.3 nm of the total length of slits (the value is independent of the pore width). By increasing the nanopore density, it is possible to increase permeability of the membrane (the ion rejection would remain the same). However, this possibility is limited above by features of the synthesis and production processes.
3.2. Salt rejection and ion selectivity
During each simulation run, we analyze the numbers of sodium and chlorine ions that are completely translocated from the feed box to the outlet through the nanopore. With the commonly used formula [23, 44], we estimate the ion rejection for two nanopore widths at four applied pressures (Fig. 2, d).
The narrow (~0.42 nm) slit-shaped nanopore in HAP demonstrates a 100% sodium ion rejection at any considered pressure drop, whereas about 10% of chlorine ions are
permeated (88 to 92% of chlorine ions are rejected depending on the pressure drop Ap). Thus, the narrow nanopore is predicted to possess a high Cl-to-Na ion selectivity (Figs. 3, a and b). We also note that in the considered range of pressure drop Ap the salt ion selectivity of the "narrow" slit-shaped nanopore is practically independent of the Ap value.
The wide nanopore (~0.78 nm) has the maximum ion rejection at the minimum considered pressure drop (Ap = 50 MPa): 86% for sodium ions and 81% for chlorine ions (Fig. 2, d). An increase in the pressure drop up to 150 MPa leads to a small but systematic decrease in the fraction of rejected ions down to 81% for Na+ and 77% for Cl-. A further increase in the pressure drop is accompanied by an accelerated deterioration of the selectivity of the nanopore. In particular, at the pressure drop 200 MPa the selectivity to both ion types becomes approximately the same and is only 72%. This value is insufficient for advanced filtering applications.
Therefore, within the pressure drop range 50-150 MPa the wide nanopore demonstrates low Na-to-Cl selectivity, namely, the difference between the numbers of filtered Cl-and Na+ ions is about 20-25%. The selectivity of the water filtration decreases nonlinearly with increasing pressure drop across the nanopore. At >200 MPa the water permeation through the wide nanopore is nonselective at least to sodium and chlorine ions.
4. Conclusion
In summary, we show that high requirements for a new generation of filter membranes can be met by using materials with slit-shaped nanopores as the basis of such membranes. The key stage in the development of such materials is the determination of the optimum geometric and electrostatic characteristics of nanopore surfaces. We showed that an effective tool for such a design is molecular level computer simulation.
Advantages of slit-shaped nanopores for highly selective filtration were shown in the paper by the example of nonequilibrium molecular dynamics simulations of salt water filtration through subnanometer slit-shaped nanopores (0.42 and 0.78 nm wide) formed by the {001} surfaces of hydroxyapatite, which is an environmentally safe, nontoxic and biocompatible mineral. In particular, it was found that the slit-shaped nanopore of width 0.42 nm demonstrates the 100% Na+ rejection and about 90% Cl- rejection in the considered range of pressure drops (50-200 MPa). At the same time, this nanopore showed the stable Cl-versus Na+ selectivity and comparatively good water permeability. The predicted water permeability of HAP with
0.42.nm wide nanopores and the 0.3 nm-1 density of pores is 196 L/(m2hbar), which is higher than that of commercial membranes and comparable with the predicted water permeability of the MoS2-based membrane and nanoporous graphene. The obtained results allow considering the biomineral with ~0.5 nm wide slit-shaped nanopores as a promising filtering material that could be used in the development of environmentally safe plants for seawater desalination and ion separation. Moreover, the results of the study invite further investigation of permeability and rejection properties of rectangular subnanometer nanopores (especially with a high aspect ratio) in nanosheets of graphene, boron nitride, silicon, MoS2, WS2, and their derivatives.
Acknowledgments
A.A. Tsukanov and E.V. Shilko gratefully acknowledge the financial support of the Russian Science Foundation (Grant No. 17-11-01232). S.G. Psakhie (contributed to the design of the research and analysis of the results) carried out the research within the Fundamental Research Program of the State Academies of Sciences for 2013-2020 (Priority direction III.23). The research was carried out using the equipment of the Research Computing Centre at the Lomo-nosov Moscow State University.
References
1. Yamada T., Matsuzaki R. Effects of slit width on water permeation through graphene membranes by molecular dynamics simulations // Sci. Rep. - 2017. - V. 8. - P. 339-1-339-9.
2. RollingsR.C., Kuan A.T., Golovchenko J.A. Ion selectivity of graphene nanopores // Nat. Commun. - 2016. - V. 7. - P. 11408-1-11408-7.
3. Zhang Y., Shi Q., Liu Y., Wang Y., Meng Z., Xiao C., Deng K., Rao D., Lu R. Hexagonal boron nitride with designed nanopores as a high-
efficiency membrane for separating gaseous hydrogen from methane // J. Phys. Chem. C. - 2015. - V. 119. - P. 19826-19831.
4. Esfandiarpoor S., Fazli M, Ganji M.D. Reactive molecular dynamic simulations on the gas separation performance of porous graphene membrane // Sci. Rep. - 2017. - V. 7. - P. 16561-1-16561-12.
5. Moghadam F., Park H.B. Two-dimensional materials: an emerging platform for gas separation membranes // Current Opinion Chem. Eng. - 2018. - V. 20. - P. 28-38.
6. Garaj S., Hubbard W., Reina A., Kong J., Branton D., Golovchenko J.A. Graphene as a subnanometre trans-electrode membrane // Nature. -2010. - V. 467. - P. 190-193.
7. Venkatesan B.M., Bashir R. Nanopore sensors for nucleic acid analysis // Nat. Nanotech. - 2011. - V 6. - P. 615-624.
8. SongS., Faleo G., YeungR., KantR., Posselt A.M., Desai T.A., Tang Q., Roy S. Silicon nanopore membrane (SNM) for islet encapsulation and immunoisolation under convective transport // Sci. Rep. - 2016. -V. 6. - P. 23679-1-23679-9.
9. Inorganic, Polymeric and Composite Membranes: Structure, Function and Other Correlations. V. 14 / Ed. by S.T. Oyama, S.M. Stagg-Williams. - Elsevier, 2011.
10. Lin J., Murad S. A computer simulation study of the separation of aqueous solutions using thin zeolite membranes // Molec. Phys. -2001. - V. 99. - No. 14. - P. 1175-1181.
11. Konatham D., Yu J., Ho T. A., Striolo A. Simulation insights for graphene-based water desalination membranes // Langmuir. - 2013. -V. 29. - No. 38. - P. 11884-11897.
12. Thomas M., Corry B. A computational assessment of the permeability and salt rejection of carbon nanotube membranes and their application to water desalination // Philos. Trans. Roy. Soc. A. - 2016. -V. 374. - No. 2060. - P. 20150020.
13. Azamat J., Sardroodi J.J. The permeation of potassium and chloride ions through nanotubes: A molecular simulation study // Monatshefte Chem. Chem. - 2014. - V. 145. - No. 6. - P. 881-890.
14. Lvova N.A., Ananina O.Y. Theoretical study of the adsorption properties ofporous boron nitride nanosheets // Comp. Mater. Sci. - 2016. -V. 115. - P. 11-17.
15. Zhao Z., Gupta K.M., He Z., Jiang J. Dipeptide crystals as reverse osmosis membranes for water desalination: Atomistic simulation // J. Phys. Chem. C. - 2018. - V. 122. - No. 11. - P. 6026-6032.
16. Tang C.Y., Zhao Y., Wang R., Hulix-Nielsen C., Fane A.G. Desalination by biomimetic aquaporin membranes: Review of status and prospects // Desalination. - 2013. - V. 308. - P. 34-40.
17. PendergastM.M., HoekE.M. A review of water treatment membrane nanotechnologies // Energy Environ. Sci. - 2011. - V. 4. - No. 6. -P. 1946-1971.
18. Qu X., Alvarez P.J., Li Q. Applications of nanotechnology in water and wastewater treatment // Water Res. - 2013. - V. 47. - No. 12. -P. 3931-3946.
19. Xu P., Zeng G.M., HuangD.L., Feng C.L., Hu S., Zhao M.H., Liu Z.F. Use of iron oxide nanomaterials in wastewater treatment: A review // Sci. Total Env. - 2012. - V 424. - P. 1-10.
20. Corry B. Designing carbon nanotube membranes for efficient water desalination // J. Phys. Chem. B. - 2008. - V. 112. - No. 5. - P. 14271434.
21. Kanani D.M., Fissel W.H., Roy S., Dubnisheva A., Fleischman A., Zydney A.L. Permeability-selectivity analysis for ultrafiltration: Effect of pore geometry // J. Membrane Sci. - 2010. - V. 349. - P. 405-410.
22. Liu W., Luo X., Bao Y., Liu Y.P., Ning G.H., Abdelwahab I., Li L, Nai C.T., Hu Z.G., Zhao D, Liu B, Quek S.Y., Loh K.P A two-dimensional conjugated aromatic polymer via C-C coupling reaction // Nature Chem. - 2017. - V. 9. - No. 6. - P. 563.
23. Yan Y.G., Wang W.S., Li W, Loh K.P., Zhang J. A graphene-like membrane with an ultrahigh water flux for desalination // Nanoscale. -2017. - V. 9. - No. 47. - P. 18951-18958.
24. Wang J., Hou T. Application of molecular dynamics simulations in molecular property prediction. 1. Density and heat of vaporization// J. Chem. Theory Comp. - 2011. - V 7. - No. 7. - P. 21512165.
25. Alshehri M.H. Interactions of boron nitride nanotubes with lipid bi-layer membranes // J. Comput. Theor. Nanosci. - 2018. - V. 15. -No. 1. - P. 311-316.
26. Tsukanov A.A., Psakhie S.G. From the soft matter-hard matter interface to bio-self-organization and hybrid systems // Phys. Mesomech. -2017. - V. 20. - No. 1. - P. 43-54.
27. Tsukanov A.A., Psakhie S.G. Energy and structure of bonds in the interaction of organic anions with layered double hydroxide nano-sheets: A molecular dynamics study // Sci. Rep. - 2016. - V. 6. -P. 19986.
28. Gogoi P., Chandravanshi M., Mandal S.K., Srivastava A., Kanau-jia S.P. Heterogeneous behavior of metalloproteins toward metal ion binding and selectivity: Insights from molecular dynamics studies // J. Biomolec. Struct. Dynamics. - 2016. - V. 34. - No. 7. - P. 14701485.
29. Iakovlev E., Zhilyaev P., Akhatov I. Atomistic study of the solid state inside graphene nanobubbles // Sci. Rep. - 2017. - V. 7. - No. 1. -P. 17906.
30. Schneemilch M., Quirke N. Free energy of adsorption of supported lipid bilayers from molecular dynamics simulation // Chem. Phys. Lett. - 2016. - V. 664. - P. 199-204.
31. Lozhkomoev A.S., Lerner M.I., Tsukanov A.A., Kazantsev S.O., Baki-na O. V., Psakhie S.G. On the possibility of soft matter nanostructure formation based on mesoporous aluminum hydroxide. Biomedical application perspectives // Phys. Mesomech. - 2017. - V. 20. - No. 2. -P. 134-141.
32. Franco L.F., Castier M, Economou I.G. Anisotropic parallel self-diffusion coefficients near the calcite surface: A molecular dynamics study // J. Chem. Phys. - 2016. - V. 145. - No. 8. - P. 084702.
33. Lin T.J. Force Field Parameters and Atomistic Surface Models for Hydroxyapatite and Analysis of Biomolecular Adsorption at Aqueous Interfaces: Doct. Diss. - The University of Akron, 2013.
34. Heinz H. The role of chemistry and pH of solid surfaces for specific adsorption of biomolecules in solution—Accurate computational mod-
els and experiment // J. Phys. Condens. Matt. - 2014. - V. 26. -No. 24. - P. 244105.
35. Heinz H., Lin T.J., Kishore Mishra R., Emami F.S. Thermodynami-cally consistent force fields for the assembly of inorganic, organic, and biological nanostructures: the INTERFACE force field // Lang-muir. - 2013. - V. 29. - No. 6. - P. 1754-1765.
36. Berendsen H.J., Postma J.P., van Gunsteren W.F., Hermans J. Interaction Models for Water in Relation to Protein Hydration // Intermolecular Forces. - Dordrecht: Springer, 1981. - P. 331-342.
37. Ryckaert J.P., Ciccotti G., Berendsen H.J. Numerical integration of the Cartesian equations of motion of a system with constraints: Molecular dynamics of n-alkanes // J. Comp. Phys. - 1977. - V. 23. -No. 3. - P. 327-341.
38. HockneyR.W., EastwoodJ.W. Computer Simulation Using Particles. -Taylor and Francis, Inc., 1988.
39. Plimpton S. Fast parallel algorithms for short-range molecular dynamics // J. Comp. Phys. - 1995. - V. 117. - P. 1-19.
40. Sadovnichy V., Tikhonravov A., Voevodin V, Opanasenko V. "Lomo-nosov": Supercomputing at Moscow State University // Contemporary High Performance Computing: From Petascale toward Exascale. -Chapman and Hall/CRC Comput. Sci. Ser., 2013. - P. 283-307.
41. Adinets A.V., Bryzgalov P.A., Voevodin V.V., Zhumatii S.A.E., Nikiten-ko D.A., Stefanov K.S. Job digest: An approach to dynamic analysis of job characteristics on supercomputers // Num. Meth. Progr. Adv. Comp. - 2012. - V. 13. - No. 4. - P. 160-166.
42. Heiranian M., Farimani A.B., Aluru N.R. Water desalination with a single-layer MoS2 nanopore // Nat. Comm. - 2015. - V. 6. - P. 8616.
43. Cohen-Tanugi D., Grossman J.C. Water desalination across nano-porous graphene // Nano Lett. - 2012. - V. 12. - No. 7. - P. 36023608.
44. Chan W.F., Chen H.Y., Surapathi A., Taylor M.G., Shao X., Marand E, Johnson J.K. Zwitterion functionalized carbon nanotube/polyamide nanocomposite membranes for water desalination // ACS Nano. -2013. - V. 7. - No. 6. - P. 5308-5319.
nocTynma b peaaKUHro 13.08.2018 r.
CeedeHua 06 aemopax
Aleksey A. Tsukanov, Cand. Sci. (Phys.-Math.), Researcher, ISPMS SB RAS, a.a.tsukanov@yandex.ru
Evgeney V. Shilko, Dr. Sci. (Phys.-Math.), Deputy Director, ISPMS SB RAS, shilko@ispms.ru
Elazar Gutmanas, Ph.D., Prof., Technion-Israel Institute of Technology, gutmanas@technion.ac.il
Sergey G. Psakhie, Dr. Sci. (Phys.-Math.), Corresponding Member of RAS, Director, ISPMS SB RAS, sp@ispms.ru
