- 1Key Laboratory of Urban Pollutant Conversion, Institute of Urban Environment, Chinese Academy of Sciences, Xiamen, China
- 2Department of Civil Engineering, and ReNEWACT Laboratory, Landmark University, Omu-Aran, Kwara State, Nigeria
- 3Key Laboratory of Marine Environment and Ecology, Ministry of Education, Ocean University of China, Qingdao, China
- 4Department of Civil and Environmental Engineering, Federal University of Technology, Akure, Nigeria
- 5Department of Agricultural and Biosystems Engineering, Landmark University, Omu-Aran, Kwara State, Nigeria
- 6Department of Geomatics, University of Benin, Benin City, Nigeria
Graphyne and its family members (GFMs) are allotropes of carbon (a class of 2D materials) having unique properties in form of structures, pores and atom hybridizations. Owing to their unique properties, GFMs have been widely utilized in various practical and theoretical applications. In the past decade, GFMs have received considerable attention in the area of water purification and desalination, especially in theoretical and computational aspects. More recently, GFMs have shown greater prospects in achieving optimal separation performance than the experimentally derived commercial polyamide membranes. In this review, recent theoretical and computational advances made in the GFMs research as it relates to water purification and desalination are summarized. Brief details on the properties of GFMs and the commonly used computational methods were described. More specifically, we systematically reviewed the various computational approaches employed with emphasis on the predicted permeability and selectivity of the GFM membranes. Finally, the current challenges limiting their large-scale practical applications coupled with the possible research directions for overcoming the challenges are proposed.
1 Introduction
The importance of water to lives and livelihood cannot be overemphasized. It is not only essential for life but it also constitute fundamental human rights of humanity. Despite its abundant availability on the earth making up to 71% of the universe, only a minute percentage of about 2.5% is available as freshwater for human consumption (Tan et al., 2020). Unfortunately, little or no access to this freshwater is gained for immediate consumption as most are trapped in glaciers and snowfields, and this ultimately results in insufficient potable water also known as water scarcity (Ahmed et al., 2020; Tan et al., 2020). This water stress condition will continue because of the increasing world population, variations in climatic condition, urbanization in conjunction with improved living standards as projected by UN and others (Gude, 2017; WWAP, 2018; United Nations, 2020). Therefore, it is pertinent that new and state-of-the-art technologies are employed to recover water fit for human use from the abundant seawater and wastewater generated by both human and industrial activities. This action will help in conserving and mitigating the future effect of water scarcity. Among many proposed solutions is desalination, a state-of-the-art technology that stems from membrane separation technology. Traditional treatment processes such as screening, filtration, chlorination/fluoridation, coagulation/flocculation, oxidation, and ion exchange (Gregory and Dhond, 1972; Birjandi et al., 2013; Abass et al., 2017; Dos Santos et al., 2017; Cescon and Jiang, 2020; Vecino and Reig, 2022) have been reported as inefficient in treating pollutants such as organic carbon and heavy metals at high levels. Removal of smaller contaminants of ionic size (<5 nm) is equally challenging due to low removal efficiencies, which is mostly achieved at high operating cost (Abu Hasan et al., 2020). On the contrary, nanoporous membrane materials processes such as, electrodialysis, reverse osmosis (RO) and nanofiltration (NF) filters at molecular level. Thus, are able to remove small molecules in form of pollutants from wastewater with removal efficiency up to 95%, at low capital/operating cost with ease of operation (Davenport et al., 2018; Yang et al., 2019; Gurreri et al., 2020; Lasisi and Zhang, 2022).
Largely, RO technique has been employed worldwide over the past decade as a preferred membrane desalination technology over the thermal desalination technologies owing to its higher efficiency and lower energy consumption (Mayor, 2019; Abdullah et al., 2021). Notwithstanding, commonly used polymeric membranes material for RO processes are replete with drawbacks of low water flux and fouling (Heiranian et al., 2015; Mahmoud et al., 2015; Dervin et al., 2016; Kim et al., 2016). However, some techniques have been utilized for the modification of RO membrane to improve its performance. For instance, the preparation of novel membrane materials (Werber et al., 2016; Li et al., 2020a). Novel nanomaterial-based membranes have proven to be uniquely suitable for desalination (Li et al., 2020a), as they allow water molecules to pass through their nanostructures at high rate while limiting the passage of dissolved salts and other solutes (Daer et al., 2015). Some nanomaterials considered as good candidates for desalination membrane include carbon nanotubes (CNTs), zeolites, graphene, and 2D transition-metal carbides, nitrides, and carbonitrides (Mxenes) (Li et al., 2007; Dumée et al., 2010; Gethard et al., 2011; Cohen-Tanugi and Grossman, 2012; Cohen-Tanugi et al., 2016; Karahan et al., 2020). Recent reports show that this new class of membranes possess exceptional flow rate compared to the conventional commercial membranes (Karahan et al., 2020; Ajibade et al., 2021; Wang et al., 2021; Damptey et al., 2022).
Graphene as an allotrope of carbon has been used extensively as a material in membrane technology for water desalination since patented in 2013 by Lockheed Martin as “Perforene” (Boretti et al., 2018). It has been found to have an atomic thickness, which guarantees high water permeability even better than most commercial NF membranes. In the construction and fabrication of desalination membranes, graphene has been used as pristine graphene, graphene oxide (GO) and in some cases as reduced graphene oxide (rGO) (Perreault et al., 2015). Interestingly, a new family of carbon allotropes called graphyne has recently been discovered as an essential material for effective water desalination. It is a one-atom-thick planar sheet with different forms of extended conjunction between acetylene and benzene groups (Kou et al., 2013; Azamat et al., 2020). In their early predictions, Baughman et al. (1987) described graphyne as stable crystalline carbon allotropes with high degrees of sp hybridization and due to the ethynyl units and aromatic moiety rings (which represents the sp- and sp2- hybridized carbon), they were christened after graphite and ethyne (Jia et al., 2017; Mehrdad and Moosavi, 2019). Furthermore, different forms of graphyne were suggested on the strength of inserting triple bonds carbon between benzene rings in graphene including three highly symmetric forms: α, β, and γ-graphyne as displayed in Figure 1. The graphyne is built on the firm connection of benzene rings by acetylene bonds ─C≡C─ (Figure 1). The pore size can be adjusted by changing the number of acetylene bonds between adjacent benzene rings. This acetylene bond is defined as n and when n = 0, the structure becomes graphene but when n = 1, 2, 3, …, n, other γ-graphyne structures are referred to as the graphyne-n which can be graphyne-1 (or graphyne) graphyne-2 (or graphdiyne) and so on. Figure 2 present the atomistic structures of the main different forms of graphyne for better understanding. In addition, graphyne and its family members (GFMs) have good chemical and mechanical properties and they are chemically inert and stable at room temperature (Enyashin and Ivanovskii, 2011; Cranford et al., 2012). They are also flexible enough to withstand deformation caused by high water pressure (Lin and Buehler, 2013) and have been successfully produced in large quantities (Enyashin and Ivanovskii, 2011; Zhu et al., 2013). Furthermore, due to the recent development in radialene and annulene chemistry, some GFMs have been successfully synthesized unto large area multilayer substrate film in form of low dimensional nanostructures, indicating it can be adapted for various applications (Tahara et al., 2013).
FIGURE 1. The structure of (A) Graphene, (B) α-graphyne, (C) β-graphyne, (D) γ-graphyne, (E) Various types of γ-graphyne-n (n = 1, 2, 3, 4, and 5).
Considering the aforementioned properties of GFMs coupled with their uniform conduction pore sizes and pristine pore edge, some authors have recently reported the application of GFMs in water desalination and water purification (Kou et al., 2014; Mehrdad and Moosavi, 2019). However, the number of available reviews focusing on graphyne and its family members for water purification/desalination are still very scanty. It is therefore necessary to look into this special area of interest and review recent literatures reports in order to be able to identify knowledge gaps and research barriers for further improvement of both experimental and theoretical approaches currently being engaged. In this review, the recent advances in the application of GFMs-based membranes for water desalination are summarized and discussed focusing on the time frame from January 2013 till October 2022. Firstly, we discussed the vast and unique properties of GFMs to justify their suitability for desalination technologies. Thereafter, we outlined the common approaches and processes utilized in the investigation of graphyne-based membrane performance used for water desalination as established in literatures. Their outcomes in terms of water permeability and salt rejection were also considered. Lastly, we concluded with sections on current challenges, outlooks and opportunities associated with the continuous and sustainable development of graphyne and its family-based membranes for practical and large-scale applications.
2 Properties of GFMs
The properties of GFMs were briefly described in the introduction section. Basically, GFMs, just as graphene have several unique properties. These properties can be broadly classified into electronics, magnetic, mechanical and atomic, thermal, optical, and structural. The following sections briefly discuss these properties as they relate to various studies.
2.1 Electronics
One of the fundamental properties of materials is the electronic property. To achieve efficient performance in application, it is important to have proper understanding of the material electronic structure in question. Therefore, proper understanding of GFMs electronic structure is key to its performance. Theoretically, GFMs possess better and excellent electronic properties with a natural band gap around 0.44–2.23 eV compare to graphene with band gaps of zero (Malko et al., 2012; Sheng et al., 2019; Sani et al., 2020). The special presence of triple carbon-carbon bond in GFMs allows for reversed chirality properties and momentum shift of their Dirac cones making it possible to achieve tunability of their energy gap (Kim, 2012). These Dirac cones were uniquely found in the electronic structure of graphene, but are now demonstrated with first principles calculations to be present in GFMs (Malko et al., 2012). GFMs show both semiconductor and metallic behavior along zigzag and armchair directions, respectively. One key parameter that is important in the electronic properties of semiconductor materials is electron mobility and GFMs have shown intrinsic electron mobility (Gao et al., 2019; Li et al., 2020b). When different layers of GFMs are stacked using varying methods, they can influence their electronic properties. León and Pacheco (2015) studied the electronic properties of bilayers GFM with different stacking configurations of either its metallic or semiconductor form, it was observed that the band gap in the semiconductor form can be modulated. Zheng et al. (2012) examined bandgap and mobility (two related electric properties parameters) and different stacking behaviors of GFMs. They discovered that benzene ring stacked in the Bernal mode resulted in the most stable stacking mode for bilayer and trilayer GFMs. Their bandgaps are smaller than the intrinsic bandgap of the single-layer GFMs and when the stacking modes of the bilayer and trilayer were revealed, they displayed metallic characteristics. Other detailed studies featuring a natural band gap and superior electric properties of GFMs with different layers and stacking models has been reported (Matsuoka et al., 2017; Gao et al., 2018; Huang et al., 2018).
2.2 Mechanical
The prediction of GFMs structures have attracted many researchers in finding out more about their mechanical properties (Peng et al., 2014). Mostly responsible for their mechanical behavior is the presence of the acetylenic groups. These acetylenic linkages are largely influenced by the Young’s modulus and fracture stress percentages. Cranford and Buehler (2011) investigated the mechanical strength using molecular dynamics (MD) simulations with ReaxFF potential. Appreciable tensile strength and fracture strain were obtained for both arm-chair and zigzag loads. In another study (Yang and Xu, 2012), discovered the disparate impacts of the addition of acetylene groups of graphyne family on the mechanical performance through acting tensile loads on the architectures. For MD simulations in the armchair direction as opposed to zigzag direction, the fracture strain in the armchair loading case remains nearly unchanged whereas the tensile strength gradually reduces with longer acetylenic chains. On the other hand, the tensile strength remains almost the same under the zigzag loading condition, whereas the fracture strain increases minimally due to the presence of fewer bond linkages and low atom density in the molecular plane of GFMs. Ajori et al. (2013) reported the influences of vacancy defects on graphyne mechanical parameters such as strain, stress, Young’s modulus and Poisson’s ratio by employing MD simulations on the basis of Tersoff–Brenner potential function. Their results demonstrated that graphyne has reduced strength and stiffness compared to graphene. Wu et al. (2017b) applied MD simulations on the basis of AIREBO (Adaptive Intermolecular Reactive Empirical Bond Order) potential function to determine fracture stress, fracture strain, and stiffness with due consideration given to both armchair and zigzag load directions. It was revealed that the number of acetylenic chains between the hexagonal rings had great impact on the mechanical properties of GFMs. However, some studies have shown that the mechanical properties of GFMs can be temperature-dependent (Shao et al., 2012; Zhang et al., 2014), although this aspect of research is still ongoing. In-plane stiffness and Poisson’s ratio are also two important elastic parameters in discussing the mechanical properties of GFMs. When compared with graphene, GFMs usually possess smaller in-plane stiffness (Kang et al., 2011; Peng et al., 2012) while the Poisson’s ratios for the GFMs could range from 0.39 to 0.87 and are all larger than that of graphene (Wen et al., 2019). In all the studies considered, the mechanical parameters have proven that GFMs are mechanically stable, which qualifies them as suitable candidates for membrane materials and applications.
2.3 Thermal/thermoelectric
The thermal or heat transport property of materials applied in nanotechnologies and other related application is of utmost importance. Varying theories have been used to explain the thermal conductivity of GFMs such as non-equilibrium MD (Galamba et al., 2007; Font et al., 2021), Green−Kubo formula with equilibrium MD (Manjunatha et al., 2021), Boltzmann transport equation (Tan et al., 2015; Jiang et al., 2017) and so on. Zhang et al. (2012) studied the thermal conductivities (TC) of four GFMs and found that the presence of the acetylenic linkages in the GFMs caused excess reduction in their thermal conductivities due to the associated low atom density in the structures and weak single bonds. 6, 6, 12-graphyne out of the four graphyne-n sheets exhibited anisotropy in the TC, different from other graphyne-n sheets considered. Ouyang et al. (2012) investigated the thermal transport property of γ-graphyne nanoribbons and showed that the thermal conductance of γ-graphyne nanoribbons is approximately 40% lower than that of graphene and the conductance was insensitive to the acetylenic linkages. Hu et al. (2015) in their study on thermal conductivity of graphyne nanotubes revealed that the perfect γ-graphyne nanotube exhibits an unprecedentedly low thermal conductivity that is much lower than those reported for ordinary, defected, and chemically functionalized carbon nanotubes. More recently, Zhang. et al. (2017a) in their study on thermal conductivity of δ-graphyne showed that as the temperature increased, the thermal conductivity of δ-graphyne monotonically decreased because of the presence of the acetylenic linkages. The outcome from all these studies is identical as they all revealed the relative difficulty of heat transport in GFMs, which is largely associated with their low thermal conductivity. Another interesting related concept is the thermoelectric performance of GFMs. As established by past findings, GFMs possess smaller thermal conductivity than graphene. However, the former predictively showed good thermoelectric performance with high figure of merit because of the band gap in GFMs, which can significantly increase its Seebeck coefficient. Some of these findings are reported in literatures (Sun et al., 2015; Tan et al., 2015; Jiang et al., 2017). Generally, an ideal model for thermoelectric material should have a dimensionless figure of merit, ZT >3. A ZT value of 3 and 4.8 was obtained for p-type holes and n-type electrons in GFMs at 293K via combining MD and first-principles simulations together with Boltzmann theory, which indicates that GFMs are potential materials for high-powered thermoelectric devices (Sun et al., 2015).
2.4 Structural and stability
Essential to producing high-performance carbon-based nanoscale materials is the clear understanding of their structural and stability properties. Design optimization and variation of their structure such as edge functionalization and ribbon width can help enhance their structural and stability properties. In addition to this, their mechanical, chemical and electronic properties are also improved in the process (Li et al., 2014). Bai et al. (2011) investigated the stability and structural properties of graphdiyne nanoribbons using density functional theory (DFT). Based on the calculated cohesive energies, one-dimensional graphdiyne nanoribbons displayed more stability than other two-dimensional graphdiyne slab in view of energy. Also, graphdiyne nanoribbons with zigzag edges are more stable than those with the armchair structures. Furthermore, the stability and structural properties of a set of fluorographynes (a typical GFM) were investigated by Enyashin and Ivanovskii (2013). By employing DFT band structure calculations, their results revealed the effect of fluorine on graphyne sheets. The stability of the fluorographynes increases as the ratio of fluorine-carbon becomes higher. Meanwhile, increasing the sp1 atom in the graphyne sheet created some reduction effect in the stability of both the pristine graphyne sheet and the fluorinated derivatives. Solis et al. (2019) examined the structural stability of graphyne and graphdiyne nanoscrolls through MD simulations. They observed that stable nanoscrolls could be created for all the studied structures resulting from the higher structural porosity of the graphyne/graphdiyne compare to graphene, and as a consequence, the π−π stacking interactions decrease. In another study, a relationship was established between the elastic property of the GFM and its structure. As the percentage of acetylenic linkages increases, the in-plane stiffness and layer modulus decrease, while the Poisson’s ratio increases. Also, dependence of structure was attributed to bond density change (Qu et al., 2017).
2.5 Magnetic
The magnetic properties of carbon materials have been considerably studied because of their applications as light non-metallic magnets and promising potential in spintronics (Chen et al., 2017; Gao et al., 2019). Generally, the zigzag direction loading graphyne nanoribbons usually bear a magnetic semiconductor ground state with ferromagnetic order both at the edge and the opposite spin directions, while the armchair direction loading graphyne nanoribbons are non-magnetic semiconductors with its band gap being a function of its width (Li et al., 2020b). The electronic structure of GFMs is well regulated by the adsorption of 3d transition metal atoms and in the process, an excellent magnetic property is impacted via spin polarized semiconductors.
Zhang et al. (2017b) investigate practically the paramagnetic characteristic of a pristine GFM material and GFM material doped with nitrogen. The results revealed that the paramagnetic characteristic of the GFM material doped with Nitrogen was increased, and show a clearer saturated magnetic moment value than the pure GFM material. However, when an asymmetric pyridinic nitrogen was introduced via substitution reaction, a huge local magnetic moment (0.98 μB) was obtained. This suggests that nitrogen of pyridine type is more favorable to improve magnetism of GFMs than just doping N sites. In a subsequent study, the authors obtained a hybrid composite system based on GFM materials and ferrous ions (Fe) doped through a simple and affordable synthetic route showing forth a favorable ferromagnetic characteristic (Zhang et al., 2017b). Despite some achievements made in the explorative studies of the magnetic properties of GFMs, more prospect can be achieved for excellent exhibition of magnetic properties applicable in spintronic devices.
2.6 Water purification and desalination
One of the interesting properties of GFMs besides all the aforementioned qualities is their ability to serve as outstanding separation materials in water purification and desalination. The ability to adjust their pore structures for easy permeation of water molecules makes them great prospect in water desalination and purification. In addition, they have high tunable surface energy and are super hydrophobic (Huang et al., 2018). Compare to graphynes, graphdiynes are semipermeable with several continuous channels that can allow the passage of some reasonable quantity of water at certain applied pressures while effectively rejecting ions (Zhu et al., 2013; Bartolomei et al., 2014). These nanopores channels can be improved via negative charging or functionalization to enhance their water permeability and salt rejection (Azamat et al., 2020). One of the key performances of the special triple carbon-carbon bond present in GFMs is the allowance for the momentum shift of their Dirac cones which makes it possible to achieve tunability in their energy gap; a condition that favors electron mobility hence, their surface functionalization (Kim, 2012; Kang et al., 2019). Some works have highlighted the attractiveness of GFMs as molecular filters for water purification and desalination. These works are elaborately presented in Section 4.
3 GFMs-based membranes: Computational methods
Reports on analysis and performance of graphyne-based membranes have largely been achieved via theoretical and computational approaches. These approaches are usually helpful when the objective is prediction or there is need for understanding an unclear phenomenon, or to explain an astonishing experimental observation (Cohen-Tanugi and Grossman, 2015). Computational approaches are also helpful as guidelines in indicating the best region or area of a material phase in a given application and in making decision as regarding the selective choice of an experimental material especially in the case of graphene sheet. Many prediction studies on graphyne and its family members have been reported. Here we focus on the most common computational modeling tools used in water desalination using membrane technology. The most common computational modeling tools used include molecular dynamics (MD) simulations, quantum mechanics (QM) calculations and combination of QM/MD simulations. DFT computations is the most used quantum mechanics method. In some other applications and subjects of concern, analytical models or finite element approaches can be employed as computational tools.
3.1 Molecular dynamics (MD) simulation
Molecular dynamics (MD) simulation (also known as molecular mechanics, MM) is a significant computational approach that enables the study of particle interactions between atoms by solving Newton’s equations of motion within a given force field. They are useful in examining the implication of the intermolecular interactions (bond stretching and angle bending) and intramolecular interactions (van der Waals and electrostatic interactions) on the behavior of ions and molecules in a system (Hosseini et al., 2019). In addition, the behavior of relatively large molecular systems even up to 109 atoms can be investigated within a physical reasonable timescale of 1 ns and 1 μs. MD simulations works well with several inputs categorized under three factors: i) Initial configuration of the system, ii) set of forcefields and iii) sets of constraints. The positions, atomic elements, bonding states and partial charges of all atoms constitutes the initial configuration of the system while the constraints include the geometric boundary conditions, a thermodynamic ensemble, and a thermostat or barostat (Cohen-Tanugi and Grossman, 2015). MD simulations in recent decade has begun to gain profound attention in water desalination application for nanomaterials such as high-aligned and high-density CNTs, zeolites, boron nitride and molybdenum disulfide (MoS2) (Li et al., 2008; Zhang et al., 2018; Oviroh et al., 2021). In a recent explanation of Cohen-Tanugi and Grossman (2015) on MD simulation study of nanoporous graphene used as an RO membrane, two water reservoirs was divided with a nanoporous graphene layer and a driving force acting as pressure or electric field, was applied to the water and ions. The trajectories provided profound information on the water flux and the solute, which were vital to the membrane water permeability and salt rejection estimation. Insights were also given on the physical dynamics of the system such as the fluid properties, membrane behavior and the physical mechanism for salt permeation.
3.2 Quantum mechanics (QM) calculation/simulation
QM simulations are best and appropriate techniques for structures wherein bonding, chemical reactions or electronic properties are of matter of concerns. The explicit nature of atoms (nuclei and electrons) in QM aids the properties of materials computation to be done with accurate quantum effects. As stated earlier, DFT is the most common quantum mechanics method, and it is an efficient tool in evaluating the interactions between molecules with active consideration of quantum effects. For predicting the ground-state electron density of atomic systems, DFT uses functionals to approximate electron correlations and exchange energies (Cohen-Tanugi and Grossman, 2015). DFT has been widely employed in graphene related studies both within and outside the scope of water desalination (Compagnini et al., 2011; Patra et al., 2016). However, its application in water desalination through graphyne-based membrane is still ongoing.
One major demerit of the QM methods is their expensive computational cost, which mostly results in the consideration of only small systems containing hundreds of atoms. However, compared to MD simulations study in water desalination and gas separation applications, they are effective in overcoming the energy barrier faced by individual molecules permeating through graphyne pores during system operation (Qiu et al., 2019). Other advance QMs, which provide higher accuracy in simulation with extended increased computational time are Quantum Monte Carlo and GW calculations (Hedin and Lundqvist, 1970; Grossman, 2002; Wang, 2011; Golze et al., 2019).
3.3 Hybrid quantum mechanics/molecular dynamics simulation
In an attempt to overcome the individual limitations of MD and QM simulations, the concept of combining these two useful computational tools have emerged. They were first proposed and introduced by Warshel and Levitt as a multiscale computational tool, which allow a reliable QM calculation on a system with a realistic modeling of the complex environment (Warshel and Levitt, 1976) and named it combined quantum mechanical and molecular mechanical (QM/MM) method or hybrid quantum mechanical and molecular dynamics QM/MD methods. The MD simulations which from tens of picoseconds to hundreds of nanoseconds are needed to gain converged statistical sampling for free energy calculations, alongside the change in the electronic structures of the system (in terms of bond forming or breaking processes) requires quantum mechanics tool such as DFT. This therefore confines the application of MD to a small number of atoms for a short time (Shen and Yang, 2018). Similarly, excessive computational cost is incurred with the QM method (Qiu et al., 2019) thus justifying this innovative concept of hybrid QM/MD.
Figure 3 illustrates the hybrid QM/MD strategy as originally introduced by Warshel and Levitt (1976) in which a system, divided into sections and explained at various levels of theory, is considered the local character of most chemical reactions in condensed phases. Therefore, it is easy to make a distinction between a “reaction region” having atoms that are directly and actively involved in the reaction and a “spectator region” where the atoms do not directly participate in the reaction (Groenhof, 2013). Furthermore, three types of interactions occur in hybrid QM/MD potential energy system as depicted in Figure 4. The magenta shaded portion represents the interactions between atoms in the QM region, the light green shaded portion represents the interactions between atoms in the MM region and the orange shaded portion represents the interactions between QM and MD atoms. Individually, the interactions between atoms within the QM and MD regions are unambiguous to explain while those of the two subsystems are relatively complex to describe (Groenhof, 2013). Notwithstanding, a careful look at Figure 3, relatively explains the complexity, such that the QM/MM energy of the total system is assumed to be equal to the energy of the isolated QM subsystem (evaluated at the QM level), in addition to the energy of the complete system (evaluated at the MM level), minus the energy of the isolated QM subsystem (evaluated at the MM level).
In recent times, QM/MD methods, which was firstly conceived in the environment of biomolecular simulation studies (Friesner and Guallar, 2005; Hu and Yang, 2009; Cui, 2016; Lu et al., 2016) has been extended into other areas of solid-state chemistry, solution chemistry, and material science (Gao, 1996; Harschend et al., 2004; Keal et al., 2011; Hofer et al., 2012; Golze et al., 2015; Hofer and Tirler, 2015). Specifically, QM/MD methods has the capacity to model time-dependent behaviors in situation with limited empirical potential functions. However, Newton’s equations of motion is used to determine the movement of each atom while DFT and ab initio (types of QM methods) are used to calculate the interatomic forces (Qiu et al., 2019).
4 GFM-based membranes: Application in water desalination and purification
One of the main areas of GFMs application, which is currently receiving much attention, is in water purification and desalination. For over two decades, advances have been made in the utilization of graphene for several applications of science and engineering, and most especially in water filtration and desalination because of its attractive and advantageous properties compared to some other nanomaterials (Aghigh et al., 2015). However, some steps have been employed in recent years to improve their performances in membrane desalination application. One of such bold steps is the utilization of graphene derivatives, and this is anchored on the minimum possible membrane thickness that can be achieved with GFM, which in turn tend to favor high permeance with appreciable salt rejection (Qiu et al., 2019). This section thus, systematically summarizes the recent research carried out from 2013 till present using GFMs for water purification and desalination from theoretical and computational viewpoints (Table 1). Some essential information on the experimental preparation and simulation process adopted for graphyne-based separation membranes were highlighted. However, emphasis was placed on two most important parameters used to characterize the performance of membrane in desalination process, which are permeability and/or selectivity.
TABLE 1. Summary of reported studies on GFMs membrane material applied in water purification and desalination.
4.1 Pristine graphyne-n membrane water desalination
Graphyne is one-atom-thick carbon allotrope that resemble graphene and it has been engaged in its bare/pristine form (i.e., unfunctionalized) in water desalination. Guo’s group in 2013 (Xue et al., 2013) reported the desalination performance of five various graphynes via MD simulations and these includes a-graphyne, ß-graphyne, and three analogues of γ-graphyne namely graphdiyne, graphyne-3, and graphyne-4. Prior to the experimental analysis, the geometric structures of the graphynes were optimized with DFTs and the constructed rectangular graphyne sheet was fixed at origin during MD simulations. The MD simulations were carried out and visualized by methods established in literatures (Humphrey et al., 1996; Phillips et al., 2005). Other parameters such as water modelling, carbon atoms particles and simulations testing with sp2 carbon in the AMBER force field were established from previous reports (Jorgensen et al., 1983; Pearlman et al., 1995; MacKerell et al., 1998). The water permeability was firstly measured across the five graphynes with ranging pressures from 0–250 MPa, and included in the simulation system is Na+ and Cl− with salt concentration of 1.2 M. The water fluxes measure as a function of the external hydrostatic pressure exhibited a linear relationship with the pressure as depicted in Figure 5A. As revealed in the work, ß-graphyne can transport water about 31% faster than α-graphyne because of its larger effective internal pore area obtained by assuming carbon atom van der Waals radius of 1.7 Å. Meanwhile, graphdiyne was impermeable to water or ions even at pressure of 250 MPa, which connotes that it is a poor candidate for desalination. However, graphyne-3 and graphyne-4 exhibited higher water permeability compared to α-graphyne and ß-graphyne in the study. With respect to salt rejections, α-graphyne, ß-graphyne, and graphyne-3 were reported to reject 100% Na+ and Cl− salt ions. Likewise, at increased pressure (up to 250 Mpa) and salt concentration (from 1.2–3.6 M), their ion rejection ability was still intact thus, indicating their independence of salt concentrations and operating pressures. This observation is in contrast to the behavior of CNTs at a concentration higher than 0.01 M (Fornasiero et al., 2008). Moreover, their rejection of other salt ions present in sea water was exceptional, as none of the ions could permeate through them. Graphyne-4 on the other hand can be permeated by water rapidly than other graphynes, but at a compromised salt rejection. A performance comparison of α-graphyne, ß-graphyne, and graphyne-3 with conventional RO membranes (seawater RO, brackish water RO, high-flux RO and NF) revealed that the permeability of these graphynes at 100% salt rejection can provide two orders of magnitude higher than that of commercial RO membranes as shown in Figure 5B. Therefore, the monolayer graphyne-3 membrane having the best performance was validated to be a promising candidate for application as membranes for water filtration and desalination (Xue et al., 2013).
FIGURE 5. (A) Water permeability across α-graphyne, β-graphyne, and graphyne-3 as a function of external hydrostatic pressure. (B) Performance comparison between α-graphyne, β-graphyne, and graphyne-3 and some conventional RO desalination membranes (C) water molecule motion trajectories without hydrostatic pressure for graphyne-2 (blue points) and graphyne-3 (green points) (D) water molecule motion trajectories with hydrostatic pressure for graphyne-2 (blue points) and graphyne-3 (green points). The black arrows show the direction of water molecule through membrane (E) water permeability for various graphyne-n membranes (N = 3–6) during the contaminants filtration simulation (F) contaminant rejection rates through various graphyne membranes under ΔP = 50 MPa (Figures (A,B) were adapted from Xue et al., 2013; (C,D) from Kou et al., 2013, and (E,F) from Lin and Buehler, 2013).
Further work on the spontaneous and continuous permeation of water molecules through a single-layer graphyne-3 membrane by molecular dynamics simulations was reported by Kou et al. (2013). Graphyne-3 with a dimension of 6.0 by 6.25 nm corresponding to its length and width, was immersed in a water bath. The monolayer graphyne was formed by replacement of carbon atoms and carbon-carbon interactions in the acetylenic linkages, which was described using AIREBO potential. Graphyne-2 and varying CNT diameters with length of 1.35 nm each were also considered for comparison. By using Gromacs 4.0.7 (Hess et al., 2008), the MD simulation was performed and the simulation system was filled with 4018 SPC/E model water molecules as described in Berendsen et al. (1987). Water molecule motion was tracked along the z-axis as a function of time in a simulation system hosting either graphyne-2 and -3 membranes with or without hydrostatic pressure. It was observed that the water molecule without the application of hydrostatic pressure navigate randomly and permeate the graphyne-3 membrane system after sometime, which was not the case in graphyne-2 system (Figure 5C). Meanwhile, after applying hydrostatic pressure, water molecule navigates readily along z-axis and permeate the membrane easily in graphyne-3 system. However, for graphyne-2 system, the water molecules could not permeate the membrane within the simulation period (Figure 5D). Furthermore, the net fluxes of water through graphyne-2, graphyne-3 and the CNTs membranes [(5, 5) CNT, (9, 0) CNT, and (10, 0) CNT] were closely related to their permeabilities. The result was in the order: (10, 0) CNT > (9, 0) CNT > graphyne-3 > (5, 5) CNT >graphyne-2 with the net water flux of graphyne-2 being zero (which is in tandem with the observation made in Figure 5). There was easy movement of water molecules via the monolayer graphyne-3 membrane than the (5, 5) CNT membrane. The net water flux in graphyne-3 membrane is 27.5 ns−1, but 13.5 ns−1 in (5, 5) CNT, even though they have similar nanopore diameters. Although, the net water fluxes of (10, 0) and (9, 0) CNTs were higher than that of graphyne-3 membrane, their work showed that single layer graphyne-3 membrane displayed a higher water permeability than CNT membranes having the same nanopore diameter (Kou et al., 2013).
The study conducted by Lin and Buehler (2013) took a different approach as 2D nanoweb-like graphyne membrane was used for water purification and also to remove contaminants from wastewater and seawater via MD simulations. They examined the relationship between the mechanical ability, filtration mechanism and rejection performance using monolayer graphyne membranes with various acetylenic linkage lengths (N = 3–6). Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) software package was used for the biaxial mechanical tensile tests while the atomistic simulations were done using ReaxFF potential developed by Chenoweth et al. (2008) to model carbon–carbon and carbon–hydrogen interactions. Gromacs 4.0 software package was used for water desalination process simulations (Hess et al., 2008). The representative contaminants examined include .6 M each of copper (II) sulfate (CuSO4), benzene (C6H6), and carbon tetrachloride (CCl4) while NaCl (completely ionized into Na+ and Cl−) was used as the saline water. The biaxial tensile tests of the graphyne membranes exhibited superior mechanical strength with ultimate stress of 16.7–32.3 GPa, which proves their high tolerance to deformations from the membrane installation process. More importantly, the water permeability of the graphyne membranes ranges from 2.9–4.5 × 10−9 mPa−1 s−1, which increases as the acetylenic linkage number increases. However, this flow rate started declining after reaching 5 and 6 acetylenic linkage as shown in Figure 5E. Also, the flow across the graphtriyne membrane exhibited an optimal purification performance which ranges from 3.0–4.0 × 10−9 mPa−1 s−1 and could still reject contaminants excellently under the applied hydrostatic pressures. Furthermore, the graphyne membranes showed higher rejection for CuSO4 and NaCl (which are hydrophilic in nature) compared to C6H6 and CCl4 (which are hydrophobic in nature), which was linked to their different interaction strengths with water molecules. The rejection rate for all contaminants decreases as acetylenic linkage number increases, thereby following the trend CuSO4 > NaCl > CCl4 > C6H6 as shown in Figure 5F.
The work of Zhu et al. (2013) corroborates the work reported by Xue et al., 2013, via the use of extensive MD simulations to examined the performance of γ-graphyne for water desalination at high rate. Parameters used for the graphyne-n sheet with acetylenic linkage number of 1–6 in their study were adapted from Narita and Nagai (1998). SPC/E model described in Berendsen et al. (1987) was used for water in the simulation system. Interaction parameter for graphyne used for Na+ and Cl−was adapted from the Amber99 force field while the non-bonding interactions are modeled by the Lennard-Jones (LJ) and Coulomb potentials. Overall system simulation was achieved using Gromacs 4.0 software. As observed by the authors, no water permeates through graphyne-1 or graphyne-2 membrane because of their small nanopores but at graphyne n with n
Sequel to the findings on the assessment of water permeation through the two graphyne membranes (graphyne-2 and graphyne-3) and CNTs membranes (Kou et al., 2013), Kou et al. (2014) extend their study to examine the potential of water desalination (both water permeability and NaCl salt rejection) via three graphyne-n membranes (n = 3, 4, and 5). The MD simulations was done using Gromacs 4.0.7 as previously established (Hess et al., 2008). As observed during the simulated period, water containing Na+ ions could permeate the three membranes at different rate and with different fresh water production capacity. Graphyne-3 among others yielded only fresh water without Na+. In addition, the salt rejection capacity tested under hydrostatic pressure of 0–350 MPa showed that graphyne-3 membrane had a perfect 100% salt rejection while graphyne-4 and graphyne-5 membranes slightly varied. This outcome contradicts the results reported by Zhu et al. (2013) who claimed that graphyne-4 membrane could also yield 100% NaCl rejection. Though, further investigation by Kou et al. (2014) showed that graphyne-4 could also yield 100% salt rejection. However, Qiu et al. (2019) in their review attributed this inconsistency to different force field parameters used during MD simulations system set-up.
As a continuum to the studies of Lin and Buehler (2013), Kou et al. (2013) and Kou et al. (2014), water and salt permeability of monolayer graphynes-n (n
According to Xue et al. (2013), graphdiyne was impermeable to water or ions even at pressure of 250 MPa, thus, making it an unsuitable candidate for desalination. However, Baghbani et al. (2020) further investigated the suitability of pristine graphdiyne membrane (without functional and chemical groups on its pore surface) for seawater desalination via MD simulation study. DFT was used to optimize the functionalized graphynes using the B3LYP theory level with the aid of GAMESS software (Schmidt et al., 1993) to obtain partial charges of all atoms. NAMD 2.18 software package was used for all the MD analyses (Phillips et al., 2005), while TIP3P model (Jorgensen et al., 1983) was used to model water interactions. Others parameters for ions and graphdiyne sheets were adapted from the CHARMM36 force field. Outcome of the experiment suggest that the graphdiyne membrane was completely impermeable at applied pressure less than 150 MPa, which is in agreement with the work by Xue et al. (2013) who reported no water passage up to 250 MPa. Interestingly at 150 MPa and even up to 650 MPa, Baghbani et al. (2020) showed that water flux could be observed, and continued to increase with increasing pressure gradient. In addition, the water permeability of the pristine graphdiyne nanosheet membrane reached 565.37 L/m2 h bar at 400 MPa, with corresponding salt rejection of 99.41%. At pressure less than 400 MPa, 100% salt rejection was realized and this performance is comparable to conventional RO system. Interestingly, this outcome exceed the performance of both graphyne-4 (90% salt rejection) and graphyne-5 (60% salt rejection) reported by Kou et al. (2014) These results thus reflects the potentials of graphdiyne nanosheet for application in seawater desalination.
Advances was made on the application of graphyne membrane for investigating brine separation performance by Zhang and Gai (2015) in a forward osmosis (FO) system. MD simulations was used to obtain information on graphyne-n (where n = 3, 4, 5, and 6) membranes. Pure water and 5% NaCl were used as the feed and draw solutions, respectively. COMPASS force field and atom-based methods was applied in the calculation of the electrostatic and van der Waals effects. As observed from the simulated FO process, graphyne-3 membrane can yield up to an average water flux of 39.15 Lcm−2 h−1 with a rejection maintained at 100%. Similarly, graphyne-4, graphyne-5, and graphyne-6 membranes all have appreciable water fluxes slightly greater than or less than that of graphyne-3 but have lower salt rejections (Figure 6A). Considering the effect of charge property of each graphyne membrane on water transport, It was observed that the water fluxes of the charged membranes were higher than those of their uncharged counterparts (Figure 6B), thus making charged graphyne-n (n = 3, 4, 5, and 6) membranes advantageous for water transport in FO systems. Overall, their study shows that graphyne-3 possess excellent separation performance for brine separation while graphyne-n (n = 4, 5, and 6) are presumed to be good candidates for brine separation of lower salt rejection demand.
FIGURE 6. (A) Water flux and rejection of Na+ and Cl− of graphyne-n (at n = 3, 4, 5, and 6). (B) Water fluxes of charged and uncharged graphyne-n membranes (at n = 3, 4, 5, and 6) (adapted from Zhang and Gai, 2015).
Studies on graphynes were not limited to only monolayer graphyne membranes as highlighted in the aforementioned works. The performance of bilayer graphyne membranes in water transport and desalination has been investigated. Akhavan et al. (2018) through a non-equilibrium molecular dynamics simulations investigated water desalination performance using double-layer graphyne membrane. Simulation parameter selection and set-up follows recorded procedures (Narita and Nagai, 1998; Cohen-Tanugi and Grossman, 2015). SPC/E model was used for water molecules and adapted OPLS parameters were applied for Na+ and Cl−ions (Berendsen et al., 1987; Lin and Buehler, 2013). All simulations were carried out using Gromacs package version 4.6.1 (Hess et al., 2008). The simulation results show that for both single and double-layer membranes, the water flow rate through the graphyne-4 membrane doubled that of graphyne-3 membrane, even though the addition of second layer membrane lessen the permeation by 50%, which is a factor of the layer spacing. At layer spacing of 0.35 nm, the flow rate reduced minimally while at spacing of 0.6 nm, the flow rate draws close to zero. For the salt rejection, graphyne-3 membranes displayed 100% rejection at all pressure ranges and layer spacing values while graphyne-4 membranes rejection is lower as expected with pressure and layer spacing grossly influencing its rejection performance. This agrees with other simulated desalination results (Lin and Buehler, 2013; Xue et al., 2013; Kou et al., 2014).
A novel concept of efficient water desalination using graphyne membranes was introduced by Nematipour et al. (2021). The authors investigated the feasibility of using anisotropically nanoporous graphyne membranes (ANGMs) to desalinate water via MD simulation process. These membranes were formed from the synthesis of nanosheets from the meta-bromination of the same molecule. In the process, an ANGM with a pore size below 1 nm can be controlled by altering the number of triple bonds between two phenyl rings as shown in Figures 7A–D. Four ANGM-n (n = 1, 2, 3, and 4) membranes family and three γ-graphyne family membranes were investigated. The MD simulation procedures employed were already established in reported literatures (Lin and Buehler, 2013; Xue et al., 2013; Zhu et al., 2013; Kou et al., 2014; Mehrdad and Moosavi, 2019). LAMMPS package (Plimpton, 1995) was used for all MD simulations analysis with Visual Molecular Dynamics (VMD) program (Humphrey et al., 1996). The result obtained by investigating the effect of external pressure as shown in Figure 7E, reveal that a linear relationship exists between the water flux and pressure, and the slope represent the membrane permeability. The order of the water flux performance are graphyne-5 > ANGM-4 > graphyne-4 > ANGM-3 > ANGM-2 > graphyne-3 > ANGM-1. Furthermore, the effect of external pressure on ion-rejection performance was also examined (Figure 7F). A slight decrease was observed in the ion rejection of all the membranes with the exception of graphyne-5. At pressure ≤150 MPa, graphyne-5 membrane exhibits a sharp decline but later increase with increasing pressure above 150 MPa. This action may be due to the rapid movement of water molecules through the membrane at high pressure, which consequently limits the movement of ions across the membrane. Previous studies have reported graphyne-3 membrane as the best membrane for desalination (Lin and Buehler, 2013; Xue et al., 2013; Kou et al., 2014), but recent studies revealed that ANGM-2 and ANGM-3 outperformed graphyne-3 membrane with higher permeabilities and 100% ions rejections. This outstanding performance hence, makes ANGM-2 and ANGM-3 membranes preferred options than graphyne-3 membrane in water desalination process.
FIGURE 7. Schematic representation of (A) ANGM-1, (B) ANGM-2, (C) ANGM-3, and (D) ANGM-4 membranes. (E) Water flux of ANGM and graphyne membranes for various external pressures (F) Ion-rejection rate of ANGM and graphyne membranes as a function of external pressure (Adapted from Nematipour et al., 2021).
4.2 Functionalized graphyne membrane for water desalination
Most preceding studies have largely used MD simulations to investigate the desalination potential of bare graphyne membranes. However, in real water desalination process or contaminated wastewater separation, these bare graphynes usually get functionalized by protons or hydroxyl radicals present in water (Psofogiannakis and Froudakis, 2012). Therefore, functionalization study is essential for its seawater desalination potential application. To address this, Raju et al. (2018) investigated the desalination performance of bare and hydrogenated (H) α-graphyne and γ-{2,3,4}-graphyne membranes using molecular dynamics simulations and upscale continuum analysis. Their results show that water flux through the membranes reduce as a function of applied external pressures with graphyne-4 and H graphyne-4 having the highest permeability. The hydrogenated γ-2-graphyne membrane has the lowest pore size and does not allow the passage of water for pressures lower than 100 MPa. Meanwhile, for the salt rejection, all the membrane tested have rejection values above 75% for pressure up to 2 GPa, with the exception of graphyne-4. Although, γ-graphyne-{2,3}, H γ-graphyne-3, α-graphyne and H α-graphyne membranes yielded high water fluxes, they were able to reject more than 90% of the ions for pressures up to 1 GPa. Moreover, upscale continuum analysis was used to complement the atomistic scale investigations. The analysis showed that the significant increase in permeability in the MD simulations did not match the real-life RO system because of some transport limitations but predict that graphyne membrane higher flux can permit about six times permeate recovery. It was further observed that salt rejection decreases with increasing pore size and applied pressure. The H membranes showed improved salt rejection performance compared to the bare ones owing to their lower pore area. Although, that was achieved with compromised water flux. Their study affirms that pore functionalization is an important factor to be considered while evaluating membranes for desalination.
In real-life desalination systems, carbon atoms are not fixed. However, most of the MD simulation studies applied in graphyne membrane desalination have fixed carbon atoms in the membrane to prevent out-of-plane displacement. Mehrdad and Moosavi (2019), therefore solves this problem by investigating the water permeability via three monolayer graphyne-n (n = 3, 4, and 5) membranes using different functional groups such as hydrogen (–H), fluorine (–F), carboxyl (COO−) and amine (NH3+) introduced into their pores. This they claim will aid the bending out of plane of the functionalized groups which in turn will align with water flow and increase the effective diameter of the pores as also supported by Raju et al. (2018). The MD simulation procedures used have been documented (Lin and Buehler, 2013; Cohen-Tanugi and Grossman, 2014), however, it is important to note that Lennard-Jones (LJ) parameter was used to model non-bonding interactions. Outcome from the study shows that water permeability through graphyne-4 was higher than graphyne-5 membrane (Figure 8A) because of the higher distribution density of nanopores on graphyne-4. This agrees with the work conducted by Zhu et al. (2013). In addition, the effect of the functional groups applied to the pores of all the graphynes-n (n = 3, 4, and 5) were examined. Results from the investigation showed that water permeabilities exhibited a dual nature. Charge distribution on the “–F” and “–H” functionalized graphynes membranes increase their water permeabilities while in NH3+ and COO- functionalized graphynes membranes, the water permeabilities decreased. Meanwhile, for each membrane, graphyne-3 and graphyne-4, ions rejection performance improved (Figure 8B). Optimally performing membranes (bare graphyne-4 and “–F” functionalized graphyne-4), were selected to investigate the influence of applied pressure. At increased pressure values, a corresponding decrease in water permeabilities was observed, although, “–F” functionalized graphyne-4 performs better than bare graphyne-4 membrane at all applied pressures (Figure 8C). Another attempt was made to examine the influence of salt concentration on the system performance using bare and F functionalized graphyne-4 membranes, respectively at 100 MPa. As shown in Figure 8D, there was decline in both the water permeability and the ion rejection of the membranes as salt concentration increases. However, in all, “–F” functionalized graphyne-4 membrane outperformed the bare graphyne-4 and other graphyne membranes considered thus, making it the best candidate for desalination.
FIGURE 8. Water permeability and ion rejection for (A) graphyne-n (n = 3,4,5), Zhu et al. (2013) (B) graphyne-n (n = 3,4,5) with functional groups tested at 100 MPa and 35 g/L (C) bare graphyne-4 and “–F” functionalized graphyne-4 at applied pressures of 100, 150, and 200 Mpa at salt concentration of 35 g/L. (D) Bare graphyne-4 and “–F” functionalized graphyne-4 tested at salt concentration of 0, 20, 35, 50, 65 g/L at applied pressure of 100 Mpa (Adapted from Mehrdad and Moosavi, 2019).
To gain better understanding of graphyne membranes pore functionalization and their influence on water desalting, Azamat et al. (2020) investigated the potential of functionalized γ-graphyne membrane pores and their allotropes for water desalination process. They employed MD simulation technique to examine one-atom-thick γ-graphyne-1 nanosheets with different functional groups on their pore edges. Three functional groups namely hydroxyl (–OH), fluorine (–F), and carboxylic acid (–COOH) were used to create the functionalized γ-graphyne-1 membranes as shown in Figure 9A–C. DFT was used to optimize the functionalized graphynes using the B3LYP theory level with the aid of GAMESS software (Schmidt et al., 1993) in order to obtain partial charges of all atoms. NAMD 2.18 and VMD 1.9.3 software were both used for all the MD runs and analyses (Humphrey et al., 1996; Phillips et al., 2005). The result obtained by investigating the effect of applied pressure showed that the water flux moving through the membrane pores increased with increasing pressure and the water flux of–OH functionalized γ-graphyne-1 membrane was higher than other functionalized γ-graphyne-1 membranes studied (Figure 9D). Meanwhile, at the absence of applied pressure, the water flux was zero. Furthermore, at low pressures (<7.5 MPa), Na+ and Cl−ions could not penetrate all the functionalized γ-graphyne-1 membranes, thereby having salt rejection of 100%. However, few ions could pass at higher applied pressure with corresponding reduction in salt rejection especially in–OH functionalized γ-graphyne-1 membrane. On the other hand, –COOH functionalized γ-graphyne-1 membrane yielded appreciable salt rejection at high applied pressure but has low permeability (Figure 9E). All the same, these functionalized γ-graphyne-1 membranes can better serve as effective tools for water desalination with high water flux and high salt rejection.
FIGURE 9. Functionalization γ-graphyne-1 membrane systems on the pore edge by (A) hydroxyl (-OH), (B) fluorine and hydroxyl (FOH) and (C) carboxylic acid (-COOH). (D) The water flux across the functionalized pores of the three γ-graphyne-1 membranes. (E) The salt rejection of functionalized pores at various applied pressures (Adapted from Azamat et al., 2020).
In addition to previously reported studies, a novel GFM was recently adapted by Jahangirzadeh et al. (2022) to explore their potential for water desalination. The group applied functionalized graphenylene nanosheet membrane for the first time for water desalination using MD simulations, although, it has previously been used for gas separation (Wang et al., 2020). A graphenylene nanosheet was prepared and functionalized with Fluorine (–F), hydrogen (H), combined fluorine & hydroxyl (–F & -OH), and combined fluorine & hydrogen (F & H) at their pore edges. The MD simulation procedures used followed those employed by Azamat et al. (2020). Undertaken at different pressures from 5–100 MPa, their results indicate that water molecules passing per 1 ns for each of the four functionalized graphenylene membrane increased with increasing pressure. The water permeability of the membranes was examined with respect to the applied pressure. Interestingly, graphenylene membrane functionalized with -F and combined “–F” & “-OH” exhibited the highest permeabilities achieved at 10 and 5 MPa, respectively. Although when compared, graphenylene membrane functionalized with “–F” & “-OH” has the highest permeability due to the combined hydrophilicity nature of the two membranes. On the other hand, membrane functionalized with “–H” has the lowest permeability because of the hydrophobic nature of “–H” atoms, which affects its hydrogenated pore. Moreover, the percentage salt rejection decreases consistently with increasing pressure. In contrast to the performance of the membrane permeability, graphenylene membrane functionalized with H exhibited the highest salt rejection percentage. Although, the percentage of salt rejection by membranes functionalized with -F and combined “–F” & “-OH” at 10 MPa and 100 MPa can reach up to 99.4% and 84.5%, respectively. Meanwhile, at pressure <10 MPa, all the functionalized graphenylene membranes had perfect salt rejections of 100%, suggesting their potentials for water desalination.
4.3 Ethanol-water separation and purification
The separation and purification of ethanol–water liquid mixture is an important process in bioethanol production and other chemical industries that deal with biofuels derived from biomass (Nalaparaju et al., 2011). Although in limited accounts, application of GFM membranes has been extended beyond water-salt separation and desalination to separation of ethanol-water mixtures. These GFM membranes could effectively purify the mixtures because of the steep energy barriers of penetration that ethanol molecules usually encounter during separation (Yeo et al., 2019).
Liu et al. (2016) via MD simulations investigated the interfacial adsorption behavior of ethanol-water mixture near the surfaces of graphyne-n sheets with n = 3, 4, and 5. All MD simulations were carried out using LAMMPS package (Plimpton, 1995). Water was modelled by the SPC/E model (Berendsen et al., 1987) while ethanol was modeled based on OPLS-AA force field (Fern et al., 2007; Metya et al., 2014). The intermolecular interactions were described by Lennard-Jones (L-J) 12–6 potential and Coulombic interaction. Other parameters used such as graphynes bond lengths and lattice were obtained from Narita and Nagai (1998). By investigating the interfacial structural properties of ethanol-water mixtures near the single-layer graphyne sheet, the density profile of ethanol and water, for ethanol concentration of 10%–90% in mole fraction, and that of pure solvent using graphyne-4 was presented (Figure 10A). It was observed that the graphyne sheet stimulated long-range ordered interfacial distribution of ethanol with the density profiles displaying sharp adsorption peaks with few lower peaks. Meanwhile, the water density profile near the surface shows a quick decline, hence forming a water depletion region. They reported similar interfacial density profiles behavior for ethanol-water mixtures near the surfaces of graphyne-3 and graphyne-5. Radial distribution function also confirmed the preferential contact of ethanol molecules with the graphyne surfaces as there are obvious peaks for ethanol near the surfaces, which were absent for water. Methyl carbon were found to be closer to graphyne carbon in RDF peak positions than that of oxygen atoms in ethanol (Figure 10B). This ethanol preferential adsorption over water connotes the micro-phase demixing or separation for ethanol-water mixtures near the graphyne surface. This demixing behavior was more pronounced in graphyne-3 surface because of the decrease in pore area compare to graphyne-4 and graphyne-5 surfaces as their nanopores predominantly occupied the ethanol molecules from the mixture. This behavior shows the strong hydrophobic interaction between amphiphilic ethanol molecules and graphyne carbon surfaces, hence making polyporous graphyne surfaces a potential two-dimensional separation membrane (Liu et al., 2016).
FIGURE 10. (A) Density profiles for ethanol and water on the graphyne-4 surface for ethanol concentration of 10%–90% (mole fraction) with the pure solvents (given in the black line). The direction of blue arrow represents the increase of ethanol concentration. The reference z = 0 corresponds to the membrane position. (B) RDF for the graphyne-4 carbon atoms with the carbon atoms of CH3– and CH2– and also the oxygen atoms in ethanol; and with oxygen atoms and hydrogen atoms in water, in different ethanol compositions (from 10 mol% to 90 mol%). (C) Fluxes of water and ethanol as a function of the applied pressures. (D) The filtered molecular number as a function of simulation time for pure water and pure ethanol through graphyne-4 membrane under 100 MPa. (E) Density distributions along z-axis for pure water and pure ethanol (left) and for the 50 mol% mixture (right) on the graphyne-4 surface [Figures (A,B) were adapted from Liu et al., 2016 while (C–E) were adapted from Zhang et al., 2019].
Similar study on ethanol/water mixtures separation using two-dimensional (2D) nanoweb graphynes was carried out by the same group (Zhang et al., 2019). They employed MD simulations to computationally investigate the permeation performance of liquid ethanol–water mixtures across polyporous 2D γ-graphyne sheets. The MD simulation utilized followed same procedures as their initial study (Liu et al., 2016). The flux of ethanol–water mixture through graphyne-4 membrane under various external hydrostatic pressures were examined and their results showed that the flux of ethanol in the mixture is higher than that of water (Figure 10C). Meanwhile, the permeation carried out for individual species over time showed that water permeate faster through the membrane than pure ethanol (Figure 10D). This behavior elucidate the competitive and selective permeation of ethanol relative to water. Furthermore, the interfacial density profiles analysis was done to further understand the permeation performances. A favorable adsorption of ethanol relative to water, on graphyne-4 membrane surface in the ethanol–water mixture with significant adsorption peak for ethanol molecules was observed (Figure 10E). This further strengthen the assertion that amphiphilic ethanol molecules are able to predominantly concentrate/adsorb on hydrophobic carbon surfaces. This is in agreement with reported experimental (Severin et al., 2014), and theoretical results (Metya et al., 2014; Zhao and Yang, 2015).
To gain further understanding of the interfacial adsorption and permeation of ethanol–water mixtures on graphynes, the group (Zhang et al., 2019), combined dispersion-corrected density functional theory (DFT-D) and classical MD simulation strategy was utilized. Force fields such as Amber, OPLS and Charm have been used to represent the interfacial interaction of graphynes (Kou et al., 2013; Lin and Buehler, 2013; Liu et al., 2016), however, these methods preferentially use sp2 carbon parameters to represent graphyne atoms. As graphyne has sp–sp2 hybridized carbon atoms, the authors believes that adopting sp2 instead of sp–sp2 could lead to underestimation of the interaction between molecules and graphyne, which in turn will lead to bias description of interfacial behavior near graphyne surfaces. Therefore, the group (Zhang et al., 2019), utilized sp–sp2 hybridized carbon atoms to cater for this bias. As such, DFT-based quantum computation was initially applied to identify the nature of the interaction of ethanol and water on graphyne-n by parameterizing force field potential, thereafter, MD simulation was used to probe the interfacial adsorption and permeation properties of ethanol–water mixtures on the graphyne surfaces. Similar to their previous study (Liu et al., 2016), the acquired density profile of ethanol display a sharp adsorption peak near the graphyne surface, which represent a strong adsorption behavior while that of water rapidly deplete near the graphyne surface. Furthermore, the permeation of the ethanol–water mixture (50 mol%) displayed considerable ethanol molecules permeation through the graphyne membranes (especially graphyne-4 pores), hence demonstrating the enrichment of ethanol on the permeate side. Furthermore, the adsorption densities on the graphyne surface using the DFT-based force field displayed stronger adsorption densities compared to the OPLS force field. This observation thus, confirms the underestimation of the binding affinity between molecules and graphyne, which has brought about the misinterpretation of the interfacial properties.
Bringing all these three studies on ethanol-water mixtures separation together, it can be inferred that, the energy barrier for each of the graphyne sheet differs depending on the orientation of the ethanol molecule relative to the pore. Additionally, the dispersion attraction of ethanol molecules was stronger as they adsorbed more to the grapyhne surfaces than water molecules. Amongst the three graphyne membrane sheet tested i.e., graphyne-3, graphyne-4 and graphyne-5, graphyne-4 membrane tends to exhibit the most desirable performance for separation of water from ethanol.
5 Challenges and outlooks
The recent progress on the works carried out on GFM in lieu of its application in water purification and desalination have been summarized, and most authors acknowledged the inherent potentials of GFM membranes in serving as ideal candidates and alternative to commercial separation membrane currently used for water desalting. However, these potentials remain latent because of some present and near future challenges such as lack of lab-scale GFM membrane fabrication, diversity in theoretical and computational analysis and desirable rejection of heavy metal/organic rejection contaminants.
5.1 Lack of lab-scale GFM membrane fabrication
It is no doubt that GFM-based membranes are still in their emerging stage of development. There are still several hurdles to cross and one of such big hurdles is to be able to successfully synthesize them at lab-scale. At present, the theoretical and computational research works carried out on the application of GFM membranes for water desalination cannot be adjudge successful until a sure pathway of preparing these membranes, first in the laboratory and eventually at scalable and economical manufacturing process. There is need to extend the characterization of these membranes beyond theoretical and computational predictions into concrete experimental validations. One major factor of consideration in this direction is the quality of GFM material to be used. Although, the synthesis and properties of GFMs have been reported in literatures (Jia et al., 2017; Gao et al., 2019; Kang et al., 2019; Gong et al., 2020), some shortcomings in synthesized GFMs such as defects in the form of wrinkles, tears and vacancies needs to be addressed if it would feature in membrane application. In principle, a defected material would adversely affect the separation performance of the membrane, and most especially the selectivity (Favre, 2022). Therefore, it is essential that a non-destructive approach is exploited. Moreover, the lab-scale prepared GFM-based membranes must possess attractive properties such as high thermal and mechanical strength so that they can withstand high loading pressure during separation applications.
5.2 Diversity in theoretical and computational analysis
Largely, all the investigations carried out to examine the application of graphyne as membrane for water desalination and purification utilized MD simulations. Researchers and research teams were at liberty to choose which computational models and parameter to employ for their works depending on their experimental and system set-up. This resulted in inconsistency in the theoretical and computational determination of the GFM membrane permeabilities and selectivity. The downside of this approach was that, it resulted in inconsistency and variation in the overall theoretical outcome presented in those studies. For instance, some MD computational studies established through their findings (Xue et al., 2013; Kou et al., 2014) that graphyne-3 membrane can achieve 100% salt rejection at higher pressure, which is not attainable with graphyne-4 and graphyne-5 membranes, even though they have higher permeabilities than graphyne-3 membrane because of their larger pores. However, the study conducted by Zhu et al. (2013) argued otherwise that graphyne-4 can also achieve 100% salt rejection even at higher pressure even with the advantage of high permeability, which then made it a better choice over graphyne-3 membrane for water desalination application as claimed in the study. This disagreement could have resulted from factors such as water model type used, force field parameters selected, pressure loading strategy adopted and so on. Another issue with the application of the theoretical and computational methods is how accurate and reliable the reported permeabilities and rejections are. DFT are useful in estimating the interactions between molecules with non-negligible quantum effect role. It uses functions to estimate the electron correlation and exchange energies so as to forecast electron density of atomic systems (Cohen-Tanugi and Grossman, 2015). In theory and computation studies of membrane water desalination, DFT studies in relation to energy landscape could be used to predict very high selectivity, although, this is usually achieved without given due consideration to the influence of some parameters such as molecule population, temperature and pressure fluctuations in real-life systems. Meanwhile, MD simulations can be used to achieve better permeability depending on the number of permeation events of desired species (Qiu et al., 2019). Therefore, to address the inconsistency issue, it is vital to combine and engage both DFT and MD methods simultaneously in a computational study of water desalination through GFMs, to be able to understand how each parameter and model affect the predicted performance. In addition, more studies should focus on combining QM and MD (i.e., hybrid quantum mechanics), to ascertain the previously established performance of the GFM membranes.
5.3 Desirable rejection of heavy metal/organic contaminants
One of the successes of conventional and lab-scale prepared RO, NF and membrane distillation membranes is seen in the efficient separation of heavy metal ions from contaminated water resources (Abass et al., 2016; Qasem et al., 2021; Lasisi et al., 2022) and organic solvents separation (Chen et al., 2019; Dai et al., 2019). Interestingly, their removal efficiency can reach up to 98.75% and 99.3% (Ozaki et al., 2002). In addition, pure and functionalized nanoporous graphenes have been employed as RO membrane either by experimental or MD simulations to examine the separation performance of heavy metal ions (Park et al., 2016; Li et al., 2017), and excellent separation were realized. However, this process (especially heavy metal/organic contaminant rejection by MD simulations) is not only limited to graphene sheets, Lin and Buehler (2013) used pure graphyne membrane to remove heavy metal/organic contaminants from wastewater and seawater via MD simulations. The membrane showed higher rejection for CuSO4 and NaCl compared to C6H6 and CCl4, which were directly linked to their different interaction strengths with water molecules. Notwithstanding, the separation performance was not comparable to those of nanoporous graphenes and as-prepared RO, MD, and NF membranes previously reported. Generally, the ion hydration of heavy metal ions is very strong due to their enlarged Coulombic interaction, and this usually results in weak dehydration effect when passing through the graphene/graphyne pores (Ferchmin, 2002) hence, a low heavy metal ions rejection is observed. One way to address this issue is through membrane functionalization. Functionalized graphene sheets have been prepared experimentally (Mishra and Ramaprabhu, 2011) and theoretically (Li et al., 2017) and used to remove high concentrations of metal salt ions. Their findings show that heavy metal ions could gain improved interaction with the pore edge of the functionalized graphene, which in turn influence the ion rejection performance. This same principle of functionalization can be applied for GFM sheets to improve their heavy metal ions rejection performance. Beside the work of Lin and Buehler (2013), which reported the application of bare graphyne membrane, no other study has considered heavy metal ions rejection performance via functionalized graphyne membranes, or even establish further the validity of bare graphyne membranes heavy metal ions rejection performance till date. Therefore, more studies that will investigate the separation performance of functionalized graphynes pores for heavy metal ions, and achieve desirable outcomes that are comparable to conventional and as-prepared RO and NF membranes are highly necessary. Besides, a thorough understanding of the molecule’s behavior inside the graphyne nanopores is also necessary as mass transport via pure graphyne pores is governed by molecular sieving mechanism. This knowledge will help membrane designers to make right selection of graphyne types and chemical functionalization that will meet the needed separation demand.
6 Conclusion
In this review, the different theoretical and computational progresses made in the research of GFMs for water purification and desalination is investigated. In the first instance, the mechanical, structural, electrical conductivity, and thermal properties of graphynes, graphdiynes, and other GFMs were analyzed for their suitability in water desalination. In addition, systematic evaluation of the different computational approaches engaged in water desalination studies were provided. Outcomes of reviewed literature on GFM membranes demonstrated their extraordinary performance in water purification and desalination. For instance, via MD simulations, graphyne-3 and graphyne-4 are reported to achieve permeability up to approximately two orders of magnitude greater than the current best-performed RO membranes experimental values, while simultaneously maintaining almost a complete 100% salt rejection. The interfacial adsorption behavior of liquid ethanol–water mixtures near the surfaces of GFM membrane sheets and their permeation performance were also examined. Although, few reports have published these findings, it was observed that graphyne-4 membrane amidst others exhibited the most desirable performance for separation of water from ethanol, and shows favorable adsorption of ethanol relative to water in the ethanol–water mixture. Moreover, efforts are currently being made in improving the separation performance of these membranes via chemical functionalization of their pore edges. Inarguably, some notable efforts have been made in the efficient application of GFMs in water purification and desalination. However, there are challenges that need to be surmounted in order to make GFM materials more attractive. Approaches that can be exploited include translation of computational studies into lab-scale experimental tool, then subsequently into industrial scale application, optimizing theoretical and computational tools to achieve consistent outcomes, and establishing basic understanding of molecules transport behavior inside graphyne nanopores to achieve improved permeability and selectivity performance. Interestingly, GFMs are endowed with unique pore characteristics that can make them ideal candidates for future water desalination. This work opens up window opportunities to explore new and emerging research strategies/pathways to overcome challenges associated with practical application of GFM membranes, which are still largely in their early stage of development.
Author contributions
The authors hereby declare the following: KL and OA participated in the conception and design, or analysis and interpretation of the data. All authors participated in drafting the article or revising it critically for important intellectual content; and approval of the final version.
Funding
This work was supported by the Landmark University Centre for Research, Innovations and Discoveries. We also acknowledge the support of ReNEWACT Laboratories, Landmark University, headed by OA for furnishing the resource used in the manuscript revision.
Acknowledgments
The authors would like to thank Landmark University and the Institute of Urban Environment, CAS for providing the platform through which this research work was conducted.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Abass, O. K., Ma, T., Kong, S., Wang, Z., and Mpinda, M. T. (2016). A novel MD-ZVI integrated approach for high arsenic groundwater decontamination and effluent immobilization. Process Saf. Environ. Prot. 102, 190–203. doi:10.1016/j.psep.2016.03.007
Abass, O. K., Zhuo, M., and Zhang, K. (2017). Concomitant degradation of complex organics and metals recovery from fracking wastewater: Roles of nano zerovalent iron initiated oxidation and adsorption. Chem. Eng. J. 328, 159–171. doi:10.1016/j.cej.2017.07.030
Abdullah, N., Yusof, N., Ismail, A. F., and Lau, W. J. (2021). Insights into metal-organic frameworks-integrated membranes for desalination process: A review. Desalination 500, 114867. doi:10.1016/j.desal.2020.114867
Abu Hasan, H., Muhammad, M. H., and Ismail, N. I. (2020). A review of biological drinking water treatment technologies for contaminants removal from polluted water resources. J. Water Process Eng. 33, 101035. doi:10.1016/j.jwpe.2019.101035
Aghigh, A., Alizadeh, V., Wong, H. Y., Islam, M. S., Amin, N., and Zaman, M. (2015). Recent advances in utilization of graphene for filtration and desalination of water: A review. Desalination 365, 389–397. doi:10.1016/j.desal.2015.03.024
Ahmed, F. E., Hashaikeh, R., and Hilal, N. (2020). Hybrid technologies: The future of energy efficient desalination – a review. Desalination 495, 114659. doi:10.1016/j.desal.2020.114659
Ajibade, T. F., Tian, H., Hassan Lasisi, K., Xue, Q., Yao, W., and Zhang, K. (2021). Multifunctional PAN UF membrane modified with 3D-MXene/O-MWCNT nanostructures for the removal of complex oil and dyes from industrial wastewater. Sep. Purif. Technol. 275, 119135. doi:10.1016/j.seppur.2021.119135
Ajori, S., Ansari, R., and Mirnezhad, M. (2013). Mechanical properties of defective γ-graphyne using molecular dynamics simulations. Mat. Sci. Eng. A 561, 34–39. doi:10.1016/j.msea.2012.10.094
Akhavan, M., Schofield, J., and Jalili, S. (2018). Water transport and desalination through double-layer graphyne membranes. Phys. Chem. Chem. Phys. 20, 13607–13615. doi:10.1039/c8cp02076k
Azamat, J., Baghbani, N. B., and Erfan-Niya, H. (2020). Atomistic understanding of functionalized γ-graphyne-1 nanosheet membranes for water desalination. J. Memb. Sci. 604, 118079. doi:10.1016/j.memsci.2020.118079
Baghbani, N., Azamat, J., Erfan-Niya, H., Majidi, S., and Khazini, L. (2020). Molecular insights into water desalination performance of pristine graphdiyne nanosheet membrane. J. Mol. Graph. Model. 101, 107729. doi:10.1016/j.jmgm.2020.107729
Bai, H., Zhu, Y., Qiao, W., and Huang, Y. (2011). Structures, stabilities and electronic properties of graphdiyne nanoribbons. RSC Adv. 1, 768–775. doi:10.1039/c1ra00481f
Bartolomei, M., Carmona-Novillo, E., Hernandez, M. I., Campos-Martinez, J., Pirani, F., Giorgi, G., et al. (2014). Penetration barrier of water through graphynes’ pores: First-principles predictions and force field optimization. J. Phys. Chem. Lett. 5, 751–755. doi:10.1021/jz4026563
Baughman, R. H., Eckhardt, H., and Kertesz, M. (1987). Structure-property predictions for new planar forms of carbon: Layered phases containing sp2 and sp atoms. J. Chem. Phys. 87, 6687–6699. doi:10.1063/1.453405
Berendsen, H. J. C., Grigera, J. R., and Straatsma, T. P. (1987). The missing term in effective pair potentials. J. Phys. Chem. 91, 6269–6271. doi:10.1021/j100308a038
Birjandi, N., Younesi, H., Bahramifar, N., Ghafari, S., Zinatizadeh, A. A., and Sethupathi, S. (2013). Optimization of coagulation-flocculation treatment on paper-recycling wastewater: Application of response surface methodology. J. Environ. Sci. Heal. - Part A Toxic/Hazardous Subst. Environ. Eng. 48, 1573–1582. doi:10.1080/10934529.2013.797307
Boretti, A., Al-Zubaidy, S., Vaclavikova, M., Al-Abri, M., Castelletto, S., and Mikhalovsky, S. (2018). Outlook for graphene-based desalination membranes. npj Clean. Water 1, 5–11. doi:10.1038/s41545-018-0004-z
Cescon, A., and Jiang, J. Q. (2020). Filtration process and alternative filter media material in water treatment. WaterSwitzerl. 12, 3377. doi:10.3390/w12123377
Chen, D., Yan, C., Li, X., Liu, L., Wu, D., and Li, X. (2019). A highly stable PBI solvent resistant nanofiltration membrane prepared via versatile and simple crosslinking process. Sep. Purif. Technol. 224, 15–22. doi:10.1016/j.seppur.2019.04.077
Chen, X., Gao, P., Guo, L., Wen, Y., Zhang, Y., and Zhang, S. (2017). Two-dimensional ferromagnetism and spin filtering in Cr and Mn-doped graphdiyne. J. Phys. Chem. Solids 105, 61–65. doi:10.1016/j.jpcs.2017.02.009
Chenoweth, K., Duin, A. C. T. V., William, A., and Goddard, I. (2008). ReaxFF reactive force field for molecular dynamics simulations of hydrocarbon oxidation. J. Phys. Chem. A 112, 1040–1053. doi:10.1021/jp709896w
Cohen-Tanugi, D., and Grossman, J. C. (2015). Nanoporous graphene as a reverse osmosis membrane: Recent insights from theory and simulation. Desalination 366, 59–70. doi:10.1016/j.desal.2014.12.046
Cohen-Tanugi, D., and Grossman, J. C. (2012). Water desalination across nanoporous graphene. Nano Lett. 12, 3602–3608. doi:10.1021/nl3012853
Cohen-Tanugi, D., and Grossman, J. C. (2014). Water permeability of nanoporous graphene at realistic pressures for reverse osmosis desalination. J. Chem. Phys. 141, 074704. doi:10.1063/1.4892638
Cohen-Tanugi, D., Lin, L. C., and Grossman, J. C. (2016). Multilayer nanoporous graphene membranes for water desalination. Nano Lett. 16, 1027–1033. doi:10.1021/acs.nanolett.5b04089
Compagnini, G., Forte, G., Giannazzo, F., Raineri, V., Magna, A. L., and Deretzis, I. (2011). Ion beam induced defects in graphene: Raman spectroscopy and DFT calculations. J. Mol. Struct. 993, 506–509. doi:10.1016/j.molstruc.2010.12.065
Cranford, S. W., Brommer, D. B., and Buehler, M. J. (2012). Extended graphynes: Simple scaling laws for stiffness, strength and fracture. Nanoscale 4, 7797–7809. doi:10.1039/c2nr31644g
Cranford, S. W., and Buehler, M. J. (2011). Mechanical properties of graphyne. Carbon N. Y. 49, 4111–4121. doi:10.1016/j.carbon.2011.05.024
Cui, Q. (2016). Perspective: Quantum mechanical methods in biochemistry and biophysics. J. Chem. Phys. 145, 140901. doi:10.1063/1.4964410
Daer, S., Kharraz, J., Giwa, A., and Hasan, S. W. (2015). Recent applications of nanomaterials in water desalination: A critical review and future opportunities. Desalination 367, 37–48. doi:10.1016/j.desal.2015.03.030
Dai, J., Li, S., Liu, J., He, J., Li, J., Wang, L., et al. (2019). Fabrication and characterization of a defect-free mixed matrix membrane by facile mixing PPSU with ZIF-8 core–shell microspheres for solvent-resistant nanofiltration. J. Memb. Sci. 589, 117261. doi:10.1016/j.memsci.2019.117261
Damptey, L., Jaato, B. N., Ribeiro, C. S., Varagnolo, S., Power, N. P., Selvaraj, V., et al. (2022). Surface functionalized MXenes for wastewater treatment—a comprehensive review. Glob. Challenges 6, 2100120. doi:10.1002/gch2.202100120
Davenport, D. M., Deshmukh, A., Werber, J. R., and Elimelech, M. (2018). High-pressure reverse osmosis for energy-efficient hypersaline brine desalination: Current status, design considerations, and research needs. Environ. Sci. Technol. Lett. 5, 467–475. doi:10.1021/acs.estlett.8b00274
Dervin, S., Dionysiou, D. D., and Pillai, S. C. (2016). 2D nanostructures for water purification: Graphene and beyond. Nanoscale 8, 15115–15131. doi:10.1039/c6nr04508a
Dos Santos, C. C., Santos, É. L., and Gonçalves, F. (2017). Evaluation of contaminants in fluorosilicic acid used for public water fluoridation in the Santos region, Brazil. Water Sci. Technol. Water Supply 17, 921–928. doi:10.2166/ws.2016.191
Dumée, L., Sears, K., Schütz, J., Finn, N., Duke, M., and Gray, S. (2010). Carbon nanotube based composite membranes for water desalination by membrane distillation. Desalin. Water Treat. 17, 72–79. doi:10.5004/dwt.2010.1701
Enyashin, A. N., and Ivanovskii, A. L. (2013). Fluorographynes: Stability, structural and electronic properties. Superlattices Microstruct. 55, 75–82. doi:10.1016/j.spmi.2012.11.022
Enyashin, A. N., and Ivanovskii, A. L. (2011). Graphene allotropes. Phys. Status Solidi Basic Res. 248, 1879–1883. doi:10.1002/pssb.201046583
Favre, E. (2022). The future of membrane separation processes: A prospective analysis. Front. Chem. Eng. 4, 1–5. doi:10.3389/fceng.2022.916054
Fern, J. T., Keffer, D. J., and Steele, W. V. (2007). Vapor - liquid equilibrium of ethanol by molecular dynamics simulation and voronoi tessellation. J. Phys. Chem. B 111, 13278–13286. doi:10.1021/jp075414u
Font, F., Micou, W., and Bresme, F. (2021). Non-equilibrium molecular dynamics and continuum modelling of transient freezing of atomistic solids. Int. J. Heat. Mass Transf. 164, 120601. doi:10.1016/j.ijheatmasstransfer.2020.120601
Fornasiero, F., Hyung, G. P., Holt, J. K., Stadermann, M., Grigoropoulos, C. P., Noy, A., et al. (2008). Ion exclusion by sub-2-nm carbon nanotube pores. Proc. Natl. Acad. Sci. U. S. A. 105, 17250–17255. doi:10.1073/pnas.0710437105
Friesner, R. A., and Guallar, V. (2005). Ab initio quantum chemical and mixed quantum mechanics/molecular mechanics (QM/MM) methods for studying enzymatic catalysis. Annu. Rev. Phys. Chem. 56, 389–427. doi:10.1146/annurev.physchem.55.091602.094410
Galamba, N., De Castro, C. A. N., and Ely, J. F. (2007). Equilibrium and nonequilibrium molecular dynamics simulations of the thermal conductivity of molten alkali halides. J. Chem. Phys. 126, 204511–11. doi:10.1063/1.2734965
Gao, J. (1996). Hybrid quantum and molecular mechanical simulations: An alternative avenue to solvent effects in organic chemistry. Acc. Chem. Res. 29, 298–305. doi:10.1021/ar950140r
Gao, X., Liu, H., Wang, D., and Zhang, J. (2019). Graphdiyne: Synthesis, properties, and applications. Chem. Soc. Rev. 48, 908–936. doi:10.1039/c8cs00773j
Gao, X., Zhu, Y., Yi, D., Zhou, J., Zhang, S., Yin, C., et al. (2018). Ultrathin graphdiyne film on graphene through solution-phase van der waals epitaxy. Sci. Adv. 4, eaat6378–8. doi:10.1126/sciadv.aat6378
Gethard, K., Sae-Khow, O., and Mitra, S. (2011). Water desalination using carbon-nanotube-enhanced membrane distillation. ACS Appl. Mat. Interfaces 3, 110–114. doi:10.1021/am100981s
Golze, D., Dvorak, M., and Rinke, P. (2019). The GW compendium: A practical guide to theoretical photoemission spectroscopy. Front. Chem. 7, 377. doi:10.3389/fchem.2019.00377
Golze, D., Hutter, J., and Iannuzzi, M. (2015). Wetting of water on hexagonal boron nitride@Rh(111): A QM/MM model based on atomic charges derived for nano-structured substrates. Phys. Chem. Chem. Phys. 17, 14307–14316. doi:10.1039/c4cp04638b
Gong, Y., Shen, L., Kang, Z., Liu, K., Du, Q., Ye, D., et al. (2020). Progress in energy-related graphyne-based materials: Advanced synthesis, functional mechanisms and applications. J. Mat. Chem. A 8, 21408–21433. doi:10.1039/d0ta08521a
Gregory, J., and Dhond, R. V. (1972). Wastewater treatment by ion exchange. Water Res. 6, 681–694. doi:10.1016/0043-1354(72)90183-2
Groenhof, G. (2013). Introduction to QM/MM simulations. Methods Mol. Biol. 924, 43–66. doi:10.1007/978-1-62703-017-5_3
Grossman, J. C. (2002). Benchmark quantum Monte Carlo calculations. J. Chem. Phys. 117, 1434–1440. doi:10.1063/1.1487829
Gude, V. G. (2017). Desalination and water reuse to address global water scarcity. Rev. Environ. Sci. Biotechnol. 16, 591–609. doi:10.1007/s11157-017-9449-7
Gurreri, L., Tamburini, A., Cipollina, A., and Micale, G. (2020). Electrodialysis applications in wastewater treatment for environmental protection and resources recovery: A systematic review on progress and perspectives. Membr. (Basel) 10, 146–193. doi:10.3390/membranes10070146
Harschend, B., Baudin, M., and Hermansson, K. (2004). A combined molecular dynamics+quantum mechanics method for investigation of dynamic effects on local surface structures. J. Chem. Phys. 120, 4939–4948. doi:10.1063/1.1635802
Hedin, L., and Lundqvist, S. (1970). Effects of electron-electron and electron-phonon interactions on the one-electron states of solids. Solid State Phys. - Adv. Res. Appl. 23, 1–181. doi:10.1016/S0081-1947(08)60615-3
Heiranian, M., Farimani, A. B., and Aluru, N. R. (2015). Water desalination with a single-layer MoS 2 nanopore. Nat. Commun. 6, 8616–6. doi:10.1038/ncomms9616
Hess, B., Kutzner, C., Van Der Spoel, D., and Lindahl, E. (2008). GROMACS 4: algorithms for highly efficient, load-Balanced, and scalable molecular simulation. J. Chem. Theory Comput. 4, 435–447. doi:10.1021/ct700301q
Hofer, T. S., Hitzenberger, M., and Randolf, B. R. (2012). Combining a dissociative water model with a hybrid QM/MM approach-A simulation strategy for the study of proton transfer reactions in solution. J. Chem. Theory Comput. 8, 3586–3595. doi:10.1021/ct300062k
Hofer, T. S., and Tirler, A. O. (2015). Combining 2d-periodic quantum chemistry with molecular force fields: A novel QM/MM procedure for the treatment of solid-state surfaces and interfaces. J. Chem. Theory Comput. 11, 5873–5887. doi:10.1021/acs.jctc.5b00548
Hosseini, M., Azamat, J., and Erfan-Niya, H. (2019). Water desalination through fluorine-functionalized nanoporous graphene oxide membranes. Mat. Chem. Phys. 223, 277–286. doi:10.1016/j.matchemphys.2018.10.063
Hu, H., and Yang, W. (2009). Development and application of ab initio QM/MM methods for mechanistic simulation of reactions in solution and in enzymes. J. Mol. Struct. THEOCHEM 898, 17–30. doi:10.1016/j.theochem.2008.12.025
Hu, M., Jing, Y., and Zhang, X. (2015). Low thermal conductivity of graphyne nanotubes from molecular dynamics study. Phys. Rev. B - Condens. Matter Mat. Phys. 91, 155408–155409. doi:10.1103/PhysRevB.91.155408
Huang, C., Li, Y., Wang, N., Xue, Y., Zuo, Z., Liu, H., et al. (2018). Progress in research into 2D graphdiyne-based materials. Chem. Rev. 118, 7744–7803. doi:10.1021/acs.chemrev.8b00288
Humphrey, W., Dalke, A., and Schulten, K. (1996). Vmd: Visual molecular dynamics. J. Mol. Graph. 14, 33–38. doi:10.1016/0263-7855(96)00018-5
Jahangirzadeh, M., Azamat, J., and Erfan-Niya, H. (2022). Molecular insight into water desalination through functionalized graphenylene nanosheet membranes. Comput. Mat. Sci. 203, 111126. doi:10.1016/j.commatsci.2021.111126
Jia, Z., Li, Y., Zuo, Z., Liu, H., Huang, C., and Li, Y. (2017). Synthesis and properties of 2D carbon - graphdiyne. Acc. Chem. Res. 50, 2470–2478. doi:10.1021/acs.accounts.7b00205
Jiang, P. H., Liu, H. J., Cheng, L., Fan, D. D., Zhang, J., Wei, J., et al. (2017). Thermoelectric properties of γ-graphyne from first-principles calculations. Carbon N. Y. 113, 108–113. doi:10.1016/j.carbon.2016.11.038
Jorgensen, W. L., Chandrasekhar, J., Madura, J. D., Impey, R. W., and Klein, M. L. (1983). Comparison of simple potential functions for simulating liquid water. J. Chem. Phys. 79, 926–935. doi:10.1063/1.445869
Kang, J., Li, J., Wu, F., Li, S. S., and Xia, J. B. (2011). Elastic, electronic, and optical properties of two-dimensional graphyne sheet. J. Phys. Chem. C 115, 20466–20470. doi:10.1021/jp206751m
Kang, J., Wei, Z., and Li, J. (2019). Graphyne and its family: Recent theoretical advances. ACS Appl. Mat. Interfaces 11, 2692–2706. doi:10.1021/acsami.8b03338
Karahan, H. E., Goh, K., Zhang, C., Yang, E., Yıldırım, C., Chuah, C. Y., et al. (2020). MXene materials for designing advanced separation membranes. Adv. Mat. 32, 1906697. doi:10.1002/adma.201906697
Keal, T. W., Sherwood, P., Dutta, G., Sokol, A. A., and Catlow, C. R. A. (2011). Characterization of hydrogen dissociation over aluminium-doped zinc oxide using an efficient massively parallel framework for QM/MM calculations. Proc. R. Soc. A Math. Phys. Eng. Sci. 467, 1900–1924. doi:10.1098/rspa.2010.0613
Kim, B. G., and Choi, H. J. (2012). Graphyne: Hexagonal network of carbon with versatile Dirac cones. Phys. Rev. B 115435, 115435–5. doi:10.1103/PhysRevB.86.115435
Kim, Y. M., Ebro, H., and Kim, J. H. (2016). Molecular dynamics simulation of seawater reverse osmosis desalination using carbon nanotube membranes. Desalin. Water Treat. 57, 20169–20176. doi:10.1080/19443994.2015.1112982
Kou, J., Zhou, X., Chen, Y., Lu, H., Wu, F., and Fan, J. (2013). Water permeation through single-layer graphyne membrane. J. Chem. Phys. 139, 064705. doi:10.1063/1.4817596
Kou, J., Zhou, X., Lu, H., Wu, F., and Fan, J. (2014). Graphyne as the membrane for water desalination. Nanoscale 6, 1865–1870. doi:10.1039/c3nr04984a
Lasisi, K. H., Ajibade, F. T., and Zhang, K. (2022). 3, 3′-diaminodiphenyl sulfone engagement in polysulfonamide-based acid-resistant nanofiltration membrane fabrication for efficient separation performance and heavy metal ions removal from wastewater. J. Memb. Sci. 661, 120909. doi:10.1016/j.memsci.2022.120909
Lasisi, K. H., and Zhang, K. (2022). Polyamine-based thin-film composite nanofiltration membrane embedded with catalytic chemical additive for enhanced separation performance and acid stability. J. Memb. Sci. 644, 120155. doi:10.1016/j.memsci.2021.120155
Lenear, C., Becton, M., and Wang, X. (2016). Computational analysis of hydrogenated graphyne folding. Chem. Phys. Lett. 646, 110–118. doi:10.1016/j.cplett.2016.01.025
León, A., and Pacheco, M. (2015). Electronic properties of β-graphyne bilayers. Chem. Phys. Lett. 620, 67–72. doi:10.1016/j.cplett.2014.12.038
Li, J., Wang, H., Yuan, X., Zhang, J., and Chew, J. W. (2020a). Metal-organic framework membranes for wastewater treatment and water regeneration. Coord. Chem. Rev. 404, 213116. doi:10.1016/j.ccr.2019.213116
Li, X., Li, B. H., He, Y. B., and Kang, F. Y. (2020b). A review of graphynes: Properties, applications and synthesis. Xinxing Tan. Cailiao/New Carbon Mater 35, 619–629. doi:10.1016/S1872-5805(20)60518-2
Li, L., Dong, J., and Nenoff, T. M. (2007). Transport of water and alkali metal ions through MFI zeolite membranes during reverse osmosis. Sep. Purif. Technol. 53, 42–48. doi:10.1016/j.seppur.2006.06.012
Li, L., Liu, N., McPherson, B., and Lee, R. (2008). Influence of counter ions on the reverse osmosis through MFI zeolite membranes: Implications for produced water desalination. Desalination 228, 217–225. doi:10.1016/j.desal.2007.10.010
Li, Y., Xu, L., Liu, H., and Li, Y. (2014). Graphdiyne and graphyne: From theoretical predictions to practical construction. Chem. Soc. Rev. 43, 2572–2586. doi:10.1039/c3cs60388a
Li, Y., Xu, Z., Liu, S., Zhang, J., and Yang, X. (2017). Molecular simulation of reverse osmosis for heavy metal ions using functionalized nanoporous graphenes. Comput. Mat. Sci. 139, 65–74. doi:10.1016/j.commatsci.2017.07.032
Lin, S., and Buehler, M. J. (2013). Mechanics and molecular filtration performance of graphyne nanoweb membranes for selective water purification. Nanoscale 5, 11801–11807. doi:10.1039/c3nr03241h
Liu, F., Yang, J., Xu, Z., and Yang, X. (2016). Selective surface adsorption and pore trapping for ethanol-water mixtures near single-layer polyporous graphynes. Appl. Surf. Sci. 387, 1080–1087. doi:10.1016/j.apsusc.2016.06.195
Lu, X., Fang, D., Ito, S., Okamoto, Y., Ovchinnikov, V., and Cui, Q. (2016). QM/MM free energy simulations: Recent progress and challenges. Mol. Simul. 42, 1056–1078. doi:10.1080/08927022.2015.1132317
MacKerell, A. D., Bashford, D., Bellott, M., Dunbrack, R. L., Evanseck, J. D., Field, M. J., et al. (1998). All-atom empirical potential for molecular modeling and dynamics studies of proteins. J. Phys. Chem. B 102, 3586–3616. doi:10.1021/jp973084f
Mahmoud, K. A., Mansoor, B., Mansour, A., and Khraisheh, M. (2015). Functional graphene nanosheets: The next generation membranes for water desalination. Desalination 356, 208–225. doi:10.1016/j.desal.2014.10.022
Malko, D., Neiss, C., Viñes, F., and Görling, A. (2012). Competition for graphene: Graphynes with direction-dependent Dirac cones. Phys. Rev. Lett. 108, 086804–4. doi:10.1103/PhysRevLett.108.086804
Manjunatha, L., Takamatsu, H., and Cannon, J. J. (2021). Atomic-level breakdown of Green–Kubo relations provides new insight into the mechanisms of thermal conduction. Sci. Rep. 11, 5597. doi:10.1038/s41598-021-84446-9
Matsuoka, R., Sakamoto, R., Hoshiko, K., Sasaki, S., Masunaga, H., Nagashio, K., et al. (2017). Crystalline graphdiyne nanosheets produced at a gas/liquid or liquid/liquid interface. J. Am. Chem. Soc. 139, 3145–3152. doi:10.1021/jacs.6b12776
Mayor, B. (2019). Growth patterns in mature desalination technologies and analogies with the energy field. Desalination 457, 75–84. doi:10.1016/j.desal.2019.01.029
Mehrdad, M., and Moosavi, A. (2019). An efficient graphyne membrane for water desalination. Polym. Guildf. 175, 310–319. doi:10.1016/j.polymer.2019.05.054
Metya, A. K., Khan, S., and Singh, J. K. (2014). Wetting transition of the ethanol-water droplet on smooth and textured surfaces. J. Phys. Chem. C 118, 4113–4121. doi:10.1021/jp4096437
Mishra, A. K., and Ramaprabhu, S. (2011). Removal of metals from aqueous solution and sea water by functionalized graphite nanoplatelets based electrodes. J. Hazard. Mat. 185, 322–328. doi:10.1016/j.jhazmat.2010.09.037
Nalaparaju, A., Zhao, X. S., and Jiang, J. W. (2011). Biofuel purification by pervaporation and vapor permeation in metal-organic frameworks: A computational study. Energy Environ. Sci. 4, 2107–2116. doi:10.1039/c0ee00630k
Narita, N., and Nagai, S. (1998). Optimized geometries and electronic structures of graphyne and its family. Phys. Rev. B - Condens. Matter Mat. Phys. 58, 11009–11014. doi:10.1103/PhysRevB.58.11009
Nematipour, A., Houshfar, E., and Fakhrabadi, M. M. S. (2021). Atomic-level engineering of anisotropically nanoporous graphyne membranes for efficient water desalination. Appl. Surf. Sci. 559, 149977. doi:10.1016/j.apsusc.2021.149977
Ouyang, T., Chen, Y., Liu, L. M., Xie, Y., Wei, X., and Zhong, J. (2012). Thermal transport in graphyne nanoribbons. Phys. Rev. B - Condens. Matter Mat. Phys. 85, 235436–235437. doi:10.1103/PhysRevB.85.235436
Oviroh, P. O., Jen, T., Ren, J., Mohlala, L. M., and Warmbier, R. (2021). Nanoporous MoS2 membrane for water desalination: A molecular dynamics study. Langmuir 37 (23), 7127–7137. doi:10.1021/acs.langmuir.1c00708
Ozaki, H., Sharma, K., and Saktaywin, W. (2002). Performance of an ultra-low-pressure reverse osmosis membrane (ULPROM) for separating heavy metal: Effects of interference parameters. Desalination 144, 287–294. doi:10.1016/S0011-9164(02)00329-6
Park, J., Bazylewski, P., and Fanchini, G. (2016). Porous graphene-based membranes for water purification from metal ions at low differential pressures. Nanoscale 8, 9563–9571. doi:10.1039/c5nr09278g
Patra, N., Esan, D. A., and Kra, P. (2016). ACS PC-2013-Dynamics of ion binding to graphene Nanostructures.pdf.
Pearlman, D. A., Case, D. A., Caldwell, J. W., Ross, W. S., Cheatham, T. E., DeBolt, S., et al. (1995). AMBER, a package of computer programs for applying molecular mechanics, normal mode analysis, molecular dynamics and free energy calculations to simulate the structural and energetic properties of molecules. Comput. Phys. Commun. 91, 1–41. doi:10.1016/0010-4655(95)00041-D
Pendergast, M. M., and Hoek, E. M. V. (2011). A review of water treatment membrane nanotechnologies. Energy Environ. Sci. 4, 1946–1971. doi:10.1039/c0ee00541j
Peng, Q., Dearden, A. K., Crean, J., Han, L., Liu, S., Wen, X., et al. (2014). New materials graphyne, graphdiyne, graphone, and graphane: Review of properties, synthesis, and application in nanotechnology. Nanotechnol. Sci. Appl. 7, 1–29. doi:10.2147/NSA.S40324
Peng, Q., Ji, W., and De, S. (2012). Mechanical properties of graphyne monolayers: A first-principles study. Phys. Chem. Chem. Phys. 14, 13385–13391. doi:10.1039/c2cp42387a
Perreault, F., Fonseca De Faria, A., and Elimelech, M. (2015). Environmental applications of graphene-based nanomaterials. Chem. Soc. Rev. 44, 5861–5896. doi:10.1039/c5cs00021a
Phillips, J. C., Braun, R., Wang, W., Gumbart, J., Tajkhorshid, E., Villa, E., et al. (2005). Scalable molecular dynamics with NAMD. J. Comput. Chem. 26, 1781–1802. doi:10.1002/jcc.20289
Plimpton, S. (1995). Fast parallel algorithms for short-range molecular dynamics. J. Comput. Phys. Available at: https://www.lammps.org.
Psofogiannakis, G. M., and Froudakis, G. E. (2012). Computational prediction of new hydrocarbon materials: The hydrogenated forms of graphdiyne. J. Phys. Chem. C 116, 19211–19214. doi:10.1021/jp306704b
Qasem, N. A. A., Mohammed, R. H., and Lawal, D. U. (2021). Removal of heavy metal ions from wastewater: A comprehensive and critical review. npj Clean. Water 4, 36. doi:10.1038/s41545-021-00127-0
Qiu, H., Xue, M., Shen, C., Zhang, Z., and Guo, W. (2019). Graphynes for water desalination and gas separation. Adv. Mat. 31, 1803772. doi:10.1002/adma.201803772
Qu, J., Zhang, H., Li, J., Zhao, S., and Chang, T. (2017). Structure-dependent mechanical properties of extended beta-graphyne. Carbon N. Y. 120, 350–357. doi:10.1016/j.carbon.2017.05.051
Raju, M., Govindaraju, P. B., Van Duin, A. C. T., and Ihme, M. (2018). Atomistic and continuum scale modeling of functionalized graphyne membranes for water desalination. Nanoscale 10, 3969–3980. doi:10.1039/c7nr07963j
Sani, S. S., Mousavi, H., Asshabi, M., and Jalilvand, S. (2020). Electronic properties of graphyne and graphdiyne in tight-binding model. ECS J. Solid State Sci. Technol. 9, 031003. doi:10.1149/2162-8777/ab801b
Schmidt, M. W., Baldridge, K. K., Boatz, J. A., Elbert, S. T., Gordon, M. S., Jensen, J. H., et al. (1993). General atomic and molecular electronic structure system. J. Comput. Chem. 14, 1347–1363. doi:10.1002/jcc.540141112
Severin, N., Sokolov, I. M., and Rabe, J. P. (2014). Dynamics of ethanol and water mixtures observed in a self-Adjusting molecularly thin slit pore. Langmuir 30, 3455–3459. doi:10.1021/la404818a
Shao, T., Wen, B., Melnik, R., Yao, S., Kawazoe, Y., and Tian, Y. (2012). Temperature dependent elastic constants and ultimate strength of graphene and graphyne. J. Chem. Phys. 137, 194901. doi:10.1063/1.4766203
Shen, L., and Yang, W. (2018). Molecular dynamics simulations with quantum mechanics/molecular mechanics and adaptive neural networks. J. Chem. Theory Comput. 14, 1442–1455. doi:10.1021/acs.jctc.7b01195
Sheng, X. L., Chen, C., Liu, H., Chen, Z., Yu, Z. M., Zhao, Y. X., et al. (2019). Two-dimensional second-order topological insulator in graphdiyne. Phys. Rev. Lett. 123, 256402. doi:10.1103/PhysRevLett.123.256402
Solis, D. A., Borges, D. D., Woellner, C. F., and Galvao, D. S. (2019). Structural and thermal stability of graphyne and graphdiyne nanoscroll structures. ACS Appl. Mat. Interfaces 11, 2670–2676. doi:10.1021/acsami.8b03481
Sun, L., Jiang, P. H., Liu, H. J., Fan, D. D., Liang, J. H., Wei, J., et al. (2015). Graphdiyne: A two-dimensional thermoelectric material with high figure of merit. Carbon N. Y. 90, 255–259. doi:10.1016/j.carbon.2015.04.037
Tahara, K., Yamamoto, Y., Gross, D. E., Kozuma, H., Arikuma, Y., Ohta, K., et al. (2013). Syntheses and properties of graphyne fragments: Trigonally expanded dehydrobenzo[12]annulenes. Chemistry 19 (34), 11251–11260. doi:10.1002/chem.201300838
Tan, H. F., Ooi, B. S., and Leo, C. P. (2020). Future perspectives of nanocellulose-based membrane for water treatment. J. Water Process Eng. 37, 101502. doi:10.1016/j.jwpe.2020.101502
Tan, X., Shao, H., Hu, T., Liu, G., Jiang, J., and Jiang, H. (2015). High thermoelectric performance in two-dimensional graphyne sheets predicted by first-principles calculations. Phys. Chem. Chem. Phys. 17, 22872–22881. doi:10.1039/c5cp03466c
United Nations (2020). Water, world water development report 2020 “water and climate change.”. Available at: https://www.unwater.org/world-water-development-report-%0A2020-water-and-climate-change/(Accessed October 2, 2022).
Vecino, X., and Reig, M. (2022). Wastewater treatment by adsorption and/or ion-exchange processes for resource recovery. WaterSwitzerl. 14, 911–914. doi:10.3390/w14060911
Wang, L., Li, F., Wang, J., Li, Y., Li, W., Yang, Y., et al. (2020). High-efficiency helium separation through an inorganic graphenylene membrane: A theoretical study. Phys. Chem. Chem. Phys. 22, 9789–9795. doi:10.1039/d0cp00154f
Wang, Y., Li, Q., Tang, G., and Zhang, N. (2021). Recent progress on carbon based desalination membranes and carbon nanomaterial incorporated non-polyamide desalination membranes. J. Environ. Chem. Eng. 9, 105762. doi:10.1016/j.jece.2021.105762
Wang, Y., Wheeler, N. V., Couny, F., Roberts, P. J., and Benabid, F. (2011). Low loss broadband transmission in hypocycloid-core Kagome hollow-core photonic crystal fiber. Ann. Appl. Stat. 5, 669–671. doi:10.1364/OL.36.000669
Warshel, A., and Levitt, M. (1976). Theoretical studies of enzymic reactions: Dielectric, electrostatic and steric stabilization of the carbonium ion in the reaction of lysozyme. J. Mol. Biol. 103, 227–249. doi:10.1016/0022-2836(76)90311-9
Wen, Y., Gao, E., Hu, Z., Xu, T., Lu, H., Xu, Z., et al. (2019). Chemically modified graphene films with tunable negative Poisson’s ratios. Nat. Commun. 10, 2446. doi:10.1038/s41467-019-10361-3
Werber, J. R., Osuji, C. O., and Elimelech, M. (2016). Materials for next-generation desalination and water purification membranes. Nat. Rev. Mat. 1, 16018. doi:10.1038/natrevmats.2016.18
Wu, B., Jin, H., Yin, J., Zhang, W., Tang, X., Zhang, P., et al. (2017a). Water and salt permeability of monolayer graph-n-yne: Molecular dynamics simulations. Carbon N. Y. 123, 688–694. doi:10.1016/j.carbon.2017.08.005
Wu, B., Tang, X., Yin, J., Zhang, W., Jiang, Y., Zhang, P., et al. (2017b). Study on mechanical properties of graphyne nanostructures by molecular dynamics simulation. Mat. Res. Express 4, 025603. doi:10.1088/2053-1591/aa5aba
WWAP (2018). The united Nations world water development report 2018: Nature-based solutions. Paris, France.
Xue, M., Qiu, H., and Guo, W. (2013). Exceptionally fast water desalination at complete salt rejection by pristine graphyne monolayers. Nanotechnology 24, 505720. doi:10.1088/0957-4484/24/50/505720
Yang, Y., and Xu, X. (2012). Mechanical properties of graphyne and its family - a molecular dynamics investigation. Comput. Mat. Sci. 61, 83–88. doi:10.1016/j.commatsci.2012.03.052
Yang, J., Xu, Z., and Yang, X. (2017). Multiscale molecular simulations on interfacial adsorption and permeation of nanoporous graphynes. Phys. Chem. Chem. Phys. 19, 21481–21489.
Yang, Z., Zhou, Y., Feng, Z., Rui, X., Zhang, T., and Zhang, Z. (2019). A review on reverse osmosis and nanofiltration membranes for water purification. Polymers 11 (8), 1252. doi:10.3390/polym11081252
Yeo, J., Jung, G. S., Martín-Martínez, F. J., Beem, J., Qin, Z., and Buehler, M. J. (2019). Multiscale design of graphyne-based materials for high-performance separation membranes. Adv. Mat. 31, 1805665. doi:10.1002/adma.201805665
Zhang, J., Cui, Y., and Wang, S. (2017a). Lattice thermal conductivity of δ-graphyne — a molecular dynamics study. Phys. E Low-Dimensional Syst. Nanostructures 90, 116–122. doi:10.1016/j.physe.2017.03.018
Zhang, M., Wang, X., Sun, H., Wang, N., Lv, Q., Cui, W., et al. (2017b). Enhanced paramagnetism of mesoscopic graphdiyne by doping with nitrogen. Sci. Rep. 7, 11535. doi:10.1038/s41598-017-11698-9
Zhang, L., Li, J., Peng, G., Liang, L., Kong, Z., Shen, J. W., et al. (2018). Charge-tunable water transport through boron nitride nanotubes. J. Mol. Liq. 258, 98–105. doi:10.1016/j.molliq.2017.12.065
Zhang, W., Xu, Z., and Yang, X. (2019). Molecular simulation of penetration separation for ethanol/water mixtures using two-dimensional nanoweb graphynes. Chin. J. Chem. Eng. 27, 286–292. doi:10.1016/j.cjche.2018.02.028
Zhang, X., and Gai, J. G. (2015). Single-layer graphyne membranes for super-excellent brine separation in forward osmosis. RSC Adv. 5, 68109–68116. doi:10.1039/c5ra09512c
Zhang, Y. Y., Pei, Q. X., Mai, Y. W., and Gu, Y. T. (2014). Temperature and strain-rate dependent fracture strength of graphynes. J. Phys. D. Appl. Phys. 47, 425301. doi:10.1088/0022-3727/47/42/425301
Zhang, Y. Y., Pei, Q. X., and Wang, C. M. (2012). A molecular dynamics investigation on thermal conductivity of graphynes. Comput. Mat. Sci. 65, 406–410. doi:10.1016/j.commatsci.2012.07.044
Zhao, M., and Yang, X. (2015). Segregation structures and miscellaneous diffusions for ethanol/water mixtures in graphene-based nanoscale pores. J. Phys. Chem. C 119, 21664–21673. doi:10.1021/acs.jpcc.5b03307
Zheng, Q., Luo, G., Liu, Q., Quhe, R., Zheng, J., Tang, K., et al. (2012). Structural and electronic properties of bilayer and trilayer graphdiyne. Nanoscale 4, 3990–3996. doi:10.1039/c2nr12026g
Keywords: GFMS, 2D materials, permeabiilty, membrane, computational analysis and strategies, selectivity, graphyne
Citation: Lasisi KH, Abass OK, Zhang K, Ajibade TF, Ajibade FO, Ojediran JO, Okonofua ES, Adewumi JR and Ibikunle PD (2023) Recent advances on graphyne and its family members as membrane materials for water purification and desalination. Front. Chem. 11:1125625. doi: 10.3389/fchem.2023.1125625
Received: 16 December 2022; Accepted: 04 January 2023;
Published: 20 January 2023.
Edited by:
Jian Li, Jiangnan University, ChinaCopyright © 2023 Lasisi, Abass, Zhang, Ajibade, Ajibade, Ojediran, Okonofua, Adewumi and Ibikunle. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Olusegun K. Abass, YWJhc3Mub2x1c2VndW5AbG11LmVkdS5uZw==, c2VndW5hYmFzczFAZ21haWwuY29t