Theoretical Study of Copper Squarate as a Promising Adsorbent for Small Gases Pollutants

Copper squarate is a metal–organic framework with an oxo-carbonic anion organic linker and a doubly charged metal mode. Its structure features large channels that facilitate the adsorption of relatively small molecules. This study focuses on exploring the potential of adsorbing small pollutants, primarily greenhouse gases, with additional investigations conducted on larger pollutants. The objective is to comprehend the efficacy of this new material in single and multiple molecular adsorption processes using theoretical methods based on density functional theory. Furthermore, we find that the molecular adsorption energies range from 3.4 KJ∙mol−1 to 63.32 KJ∙mol−1 depending on the size and number of adsorbed molecules. An exception is noted with an unfavorable adsorption energy value of 47.94 KJ∙mol−1 for 4-nitrophenol. More importantly, we demonstrate that water exerts an inhibitory effect on the adsorption of these pollutants, distinguishing copper squarate as a rare MOF with hydrophilic properties. The Connolly surface was estimated to give a more accurate idea of the volume and surface accessibility of copper squarate. Finally, using Monte Carlo simulations, we present a study of adsorption isotherms for individual molecules and molecules mixed with water. Our results point out that copper squarate is an efficient adsorbent for small molecular pollutants and greenhouse gases.


Introduction
Nowadays, scientists around the world are actively seeking new porous materials to tackle the pollution crisis caused by urbanization and emissions from motor vehicles and industry [1,2].Extensive research has been conducted, resulting in the publication of in-depth studies on air emissions and the various types of pollutants present in our atmosphere, such as carbon dioxide (CO 2 ), the main contributor to the greenhouse effect, which is largely responsible for global warming.Additionally, methane (CH 4 ), tropospheric ozone (O 3 ), dinitrogen monoxide (N 2 O), nitrogen trifluoride (NF 3 ), and sulfur hexafluoride (SF 6 ) are critical trace gases that contribute to climate change.These six compounds are among the main greenhouse gases identified in the Kyoto Protocol [3][4][5][6].
In addition, the combustion of coal, fossil fuels, or biomass also releases air pollutants, such as sulfur dioxide (SO 2 ) and volatile organic compounds (VOCs), which contribute significantly to industrial and vehicle exhaust pollution.At the same time, VOCs produce a variety of pollutants, including polycyclic aromatic hydrocarbons (PAHs) derivatives such as NPAH nitrates.This category includes, for example, 4-nitrophenol (4NP), which is considered one of the most hazardous pollutants in our daily lives [7][8][9][10][11][12][13][14][15].Nitrophenols are classified as priority pollutants by the U.S. Environmental Protection Agency [16].
Due to their high toxicity, atmospheric pollutants have considerable negative effects on human health, wildlife, and global warming.For this reason, emphasis has been placed on experimental and theoretical studies concerning the application of separation, adsorption, gas capture and storage, and selective catalysis as potential solutions [17][18][19][20][21][22].
In this paper, we will focus on a type of MOF [34] that has not yet been explored for adsorption or capture processes: the copper squarate (C-S) metal-organic framework.The crystal structure of this material was first analyzed and described in detail by Dinnebier et al. [35,36].C-S was subsequently investigated by Colmenero et al. [37] using theoretical solid-state methods based on density functional theory (DFT).This study firmly established the crystal structure of copper squarate, previously determined using a restricted rigid-body Rietveld refinement, which discovered that it exhibits the phenomenon of negative linear compressibility due to the presence of the empty channel structural motif.Colmenero et al. demonstrated that other squarate MOFs, such as zinc [36] and uranyl squarates [38], also exhibit negative linear compressibility.In fact, Colmenero [39] has shown that oxo-carbonic acids in the solid state exhibit extremely anomalous mechanical behavior, and recently, Qiu et al. [40] have revealed that lead and barium squarates exhibit zero linear compressibility and thermal expansion effects over a wide range of pressures.The presence of empty spaces inside these structures and their large pore size have led us to believe in the viability of these materials for the capture and storage of hazardous gases.As far as we know, squarate MOFs have only been recently studied for molecular adsorption applications.Namely, the use of cobalt squarate MOF has recently been investigated for CO 2 /N 2 separation by Zhang et al. [41].
The aim of our work is to characterize the behavior of squarate materials in molecular adsorption applications, which have not yet been explored, using solid simulation techniques.Specifically, we will investigate the ability of copper squarate to absorb a series of different classes of gaseous pollutants in the presence and absence of water.We will closely examine its response to these pollutants, determining its preferred adsorption sites, energies, and mechanisms, as well as its adsorption isotherms.Unfortunately, to our knowledge, no experimental or theoretical studies have been published on adsorption in copper squarate.The present work is therefore the first to provide theoretical results that will serve as a starting point for further experiments.The exploration of new materials for the adsorption of contaminants is crucial to counter the wave of pollution caused by urbanization that we are currently experiencing.

Results and Discussion
The copper squarate crystal structure was re-optimized, and the results obtained are fairly consistent with those obtained by Colmenero et al. [36].The powder X-ray diffraction (PXRD) patterns of copper squarate were derived from the computed and experimental structures using the program REFLEX included in the Materials Studio program suite [42].The results are displayed in Figure S1 of the Supplementary Materials and, as can be seen, they are consistent.Because of the small size of the unit cell, it was doubled along the axis to facilitate the adsorption of the pollutant, as was conducted in previous works studying the adsorption of similar molecules in materials like zeolite [43].The calculations were then continued with a series of optimizations of the double unit cell containing the pollutant gases.These gases were introduced inside using the Adsorption Locator tool of Materials Studio [42] in order to conduct the adsorption processes inside this material.The results are presented in Table 1, where the lattice parameters of the crystal unit and the double cell of copper squarate are provided, as well as the parameters of the cells distorted due to the capture of one or more pollutants and water molecules.The optimized geometries of different species are compared with the corresponding experimental data in Table S1 of the supporting information.The optimized structure of copper squarate is shown in Figure 1.The optimized copper squarate unit and double cell can be downloaded from the Supplementary Information (Table S2).
Table 1.Structural parameters of copper squarate (C-S) are in the pure form and contain adsorbed pollutant molecules and water separately.The experimental unit cell parameters are from Ref. [36].

Crystal Bulk Structure and Analysis of Connolly Free Volume and Surface Area Properti
In this section, we employ another feature of Materials Studio 22.0 [42], the "At Volume and Surfaces" to approximate the Connolly surface of copper squarate and pollutants separately using the probe radius.This allows us to thoroughly explore cavity and surface features of the internal complex structure after the adsorption pro [44] and to obtain the occupied and free volumes, as well as surface area (Table 2).Th parameters are then used to calculate the pore volume using the expression below [45

Crystal Bulk Structure and Analysis of Connolly Free Volume and Surface Area Properties
In this section, we employ another feature of Materials Studio 22.0 [42], the "Atoms Volume and Surfaces" to approximate the Connolly surface of copper squarate and the pollutants separately using the probe radius.This allows us to thoroughly explore the cavity and surface features of the internal complex structure after the adsorption process [44] and to obtain the occupied and free volumes, as well as surface area (Table 2).These parameters are then used to calculate the pore volume using the expression below [45]: From Table 2, we notice a good agreement between the experimental density (2.427 g•cm −3 ) and the theoretical value (2.263 g•cm −3 ), showing the accuracy of the DFT-D3 method in describing unit cell structures in crystalline materials.The total pore volume of copper squarate (1 × 1 × 2 cell) is expected to be 0.0958 g•cm −3 , i.e., 23% of the total cell volume reported in Table 1.This value is close to the experimental one (24.6%)approximated from available data (density and cell volume) in the work of Colmenero et al. [36] considering the (1 × 1 × 2) cell, as well as using the theoretical cell free volume.All these calculations use a probe radius of 1.0 Å and a grid interval of 0.6 Å [45].
To delve deeper into the understanding of C-S capacity as a filter and elucidate the outcomes of the Connolly analysis, we employed the DFT-D3 veracity method.This method allows us to ascertain the maximum quantity of adsorbed pollutants within the bulk.Table 3 shows the structural parameters of the co-adsorption process.

Adsorption Analysis
Achieving optimal DFT performance in this study involves a hierarchical approach.First, we locate the most stable configuration in conformational space, as was previously described using the Adsorption Locator module with COMPASS III force field [46], by dropping the molecule at various pressures employing the Metropolis algorithm.The DFT geometry is then optimized for the most favorable configuration.(Figure 2) shows the single final adsorption structure of pollutants and water in C-S and (Table 4) gives the energy data for the adsorption process.Negative adsorption energies indicate favorable capture of gases [31].2. For all four species, we observe a short consistent interaction involving C-S oxygen with C (CO 2 ), N (N 2 O), S (SO 2 ), and strong hydrogen [47] bound with H (H 2 O), comprising a bound distance of 2.811, 2.954, 2.775, and 1.966 Å, respectively.In addition, a weaker interaction is evident with C (C-S), interacting with O (CO 2 ), N (N 2 O), O (SO 2 ), and O (H 2 O) at distances ranging from 2.847-3.080Å.
In Table 4, Generally, metal-organic frameworks (MOFs) are known to be unstable in humid environments.Zeolite has been highlighted in the literature for its hydrophobic nature [48][49][50].However, recent experimental studies have identified exceptions, such as ZIF-8 and ZIF-90, which show stability in the presence of water [51,52].In our study, we theoretically demonstrate that C-S exhibits stability in the presence of water, which is consistent with the hydrophilic behavior observed in some MOFs.In particular, we show that O-H bond lengths in water approach 0.966 Å, which closely matches the experimental value of 0.98 Å [52].
Methane is a non-polar molecule and has been the subject of previous adsorption studies, particularly on zeolites [45].In our study, different quantities of CH 4 ranging from one to four were absorbed and optimized to identify the most favorable insertion sites.Our results are consistent with those of a previous study in which four methane molecules were observed inside the ZIF-4 cell [32].By analyzing the modeled structures, we observed that all CH 4 molecules are in close proximity to the C-S, involving a hydrogen interaction with the oxygen atom of C-S, with a bonding distance of 2.737 Å, as illustrated in (Figure 2).This interaction represents a limiting case of hydrogen bonding, characterized by a weak donor to a strong acceptor due to their substantial contribution to dispersion [47].This contributes to system stabilization, resulting in negative adsorption values ranging from −19.42 KJ•mol −1 to −23.02 KJ•mol −1 for one and four CH 4 uptake molecules, respectively.
From Table 4, we can observe high adsorption energies of −30.51, −50.05, and −26.94 KJ•mol −1 for the adsorption of one, three, and five O 3 molecules, respectively, including a strong uptake within the bulk.This is attributed to the reactive nature and high polarity of ozone due to the presence of oxygen atoms.Capture involves intermolecular Van der Waals forces, mainly O O3 -O C-S with an average bond length of 2.8 Å for all the optimized structures.In this case, the charged fluctuation of the oxygen and the electronic density leads to a sensitive instantaneous dipole moment, creating a type of London dispersion.
The vacuum inside the cavity allows the capture of one molecule with fluorine atoms due to the steric effect of the electronic density of the halogen, causing a repulsion with the oxo-carbonic anion of the copper squarate.We observed that the adsorption energy is −14.12KJ•mol −1 for NF 3 , while it decreases consistently for CF 4 and SF 6 , with values of −3.40 and −3.78 KJ•mol −1 , respectively.We deduce that as the number of fluorine atoms increases, the adsorption energy decreases.Additionally, we notice that when the central atom is nitrogen, the adsorption energy is higher, as seen with N 2 O.The optimized lengths are shown in Figure 2, where it can be seen that the carbon of C-S interacts with the fluorine of CF 4 and NF 3 at 2.7 and 2.8 Å, while a shorter value of 1.98 Å is noted with SF 6 because its octahedral form minimizes the distance between the pollutant and the bulk.
A positive adsorption energy of +47.94 KJ•mol −1 was noted for 4NP (Table 4), indicating unfavorable capture of the aromatic pollutant by the C-S.This is in contradiction with the adsorption locator calculation, which allows the uptake of one molecule.This implies the lack of accuracy of the force field and the foremost importance of using the DFT-D3 method to complement and achieve a reliable calculation.

Co-Adsorption of Pollutants and Water
In this section, we study the effect of a humid environment on the adsorption inside copper squarate bulk.The adsorption energy and interaction energy are calculated using Equations ( 2) and (3), respectively, and the results are presented in Table 4.
According to Table 4, the co-adsorption of the pollutant with water gives negative adsorption energy, revealing favorable adsorption energy, except for 4NP.The latter has a positive adsorption value in both dry and humid environments, which can be caused by the unavailability of sufficient space on the pores due to its large size and radius.
For CH 4 , SO 2 , and CO 2 , the adsorption energies increased in the presence of water to reach −33.11, −43.48, and −40.39 KJ•mol −1 , respectively.This implies the same interaction as above between the pollutant and the bulk.In addition, weak hydrogen bonds are formed between the water and the pollutant, which further increases the adsorption energy.The same trend is noted for molecules with fluorine atoms, where the values of the mixed adsorption energy increased to reach −26.83, −41.31, and −6.186 KJ•mol −1 for NF 3 , CF 4 , and SF 6 , respectively.The important impact on the adsorption energetic value is due to the multiple interactions caused by the reactive sites of these types of pollutants.In contrast, for ozone and dinitrogen monoxide, the adsorption energy decreases slightly but remains favorable to the adsorption process.Most of the interaction energies are positive, indicating repulsion between the species and far adsorption between the small pollutants and the water inside the pore.This is due to the electron density, as depicted in Figure 3.
Finally, the adsorption of 4NP in the presence of water does not counterbalance capture towards a favorable process.In previous work on nitrophenol on activated carbon bulk under humid conditions, the authors observed a possibility of adsorption by water as it moved to place itself between 4NP and the surface bulk, creating hydrogen bonds and electrostatic interactions.This suggests a lack of space inside the C-S pore, which can create repulsion due to the high charge density of water, 4NP, and C-S cell, as can be seen in Figure 3 where water is placed in the cell extremity, not inside.We also noticed this in the previous section on the Connolly free volume, where we noted that 4NP and 2NP have almost the same volume.The interaction energy here is negative due to the hydrogen bonding of the aromatic compound to the copper squarate.Finally, the adsorption of 4NP in the presence of water does not counterbalance capture towards a favorable process.In previous work on nitrophenol on activated carbon bulk under humid conditions, the authors observed a possibility of adsorption by water as it moved to place itself between 4NP and the surface bulk, creating hydrogen bonds and electrostatic interactions.This suggests a lack of space inside the C-S pore, which can create repulsion due to the high charge density of water, 4NP, and C-S cell, as can be seen in Figure 3 where water is placed in the cell extremity, not inside.We also noticed this in the previous section on the Connolly free volume, where we noted that 4NP and 2NP have almost the same volume.The interaction energy here is negative due to the hydrogen bonding of the aromatic compound to the copper squarate.

CS-CF4-H2O CS-SF6-H2O CS-O3-H2O
These fluctuations in no way interfere with copper squarate's favorable ability to capture these small pollutants.On the contrary, they positively highlight the fact that the adsorbent can capture them in both dry and humid states, without the presence of water playing a counterbalancing role.It is important to remember that the small pollutants, as well as water, are energetically stable separately in the C-S.
The main geometric parameters of the co-adsorption (water/pollutant) were analyzed and are reported in Table 5 above.The results indicate that there is no significant change These fluctuations in no way interfere with copper squarate's favorable ability to capture these small pollutants.On the contrary, they positively highlight the fact that the adsorbent can capture them in both dry and humid states, without the presence of water playing a counterbalancing role.It is important to remember that the small pollutants, as well as water, are energetically stable separately in the C-S.
The main geometric parameters of the co-adsorption (water/pollutant) were analyzed and are reported in Table 5 above.The results indicate that there is no significant change in the bond distance and valence angle between them compared with the isolated adsorbed form.This is due to the presence of more free space, allowing more water molecules to be captured, confirming the observation that the interaction energy is positive.The presence of additional adsorbed molecules can lead to an increase in distortion due to space limitations.This phenomenon has already been observed in solid simulations, notably in studies of methane, where the angle of increased from 109 • to 114 • [45].

Sorption Isotherms
The isothermal sorption of pollutants and water in singular and mixed forms was studied at 298 K using the "sorption isotherm" module implemented in Materials Studio 22.0 [42].This study is based on random moves to perform energy changes for the new conformation.This is by measuring the chemical potential, expressed as the gas fugacity with pressure inside the cell bulk.The results are shown in Figure 4, indicating the average adsorption capacity of gas molecules by the copper squarate.

Quantum computational Methods
All calculations for the solid-state compounds considered in this work were performed with CASTEP code [53] using theoretical methods based on DFT using norm-conserving pseudopotentials and plane wave basis sets [54].Due to the large pore size of MOFs, which contain hundreds of atoms or more, quantum mechanical methods such as By analyzing the adsorbed singular graphs in Figure 4, we observe that most of the pollutants present a consistent loading profile with molecules inserted around zero or zero per the considered unit cell in low-or high-pressure variations.For example, N 2 O shows zero loaded molecules, while the adsorption study shows high negative adsorption energy of −63.32 kJ/mol, indicating a favorable capture.This aligns also with the Connolly studies, which show available space inside copper squarate.These particular results contradict the precise results obtained by the DFT-D3 method.These strange results can be attributed to the limitations of the force field used in this C-S model, which includes pressure and is unable to accurately describe the unusual properties of the material, in particular its negative linear compressibility [36], highlighting the necessity of complementing it with ab initio calculations.In contrast, water is well described, showing an uptake loading of 15 molecules per cell, which is in good agreement with CASTEP calculations.This agreement may indicate the established hydrophilic nature of the solid.Similarly, ozone molecules show comparable behavior due to their common nature with water.
In the presence of water, the line graph representing the pollutant does not change.There are two reasons for this.The first reason is the same as when they are loaded individually.The second reason may be due to the hydrophilicity of water, which dominates the pollutants and competes with them, reducing their loading capacity in the cell.This phenomenon is evident with ozone, which loads around three molecules instead of six in a dry environment.This is due to the stronger hydrogen bonds between the water and the copper squarate to the intermolecular interactions between the pollutants and the bulk.

Quantum Computational Methods
All calculations for the solid-state compounds considered in this work were performed with CASTEP code [53] using theoretical methods based on DFT using norm-conserving pseudopotentials and plane wave basis sets [54].Due to the large pore size of MOFs, which contain hundreds of atoms or more, quantum mechanical methods such as post-HF are impractical for such large systems [55].The geometry optimizations were performed employing the Broyden-Fletcher-Goldfarb-Shanno scheme (BFGS) [55,56].We treated the exchange-correlation energy term using the generalized gradient approximation (GGA) with Perdew-Burke-Ernzerhof (PBE) density functionals [57].Since the pollutant molecules are packed together inside the unit cell of copper squarate, their description requires the inclusion of Van der Waals forces.For this purpose, Grimme's empirical dispersion correction method (DFT-D3) for dispersive vdW interactions was used in calculations [58,59].Customized atom pseudopotentials [60] generated with the CASTEP code and a plane wave kinetic energy cut-off of 1000 eV were utilized.Convergence criteria were set at 2 × 10 −5 eV/atom for energy, 0.05 eV/Å for interatomic forces, 0.1 GPa for maximum stress, and 0.002 Å for displacements.Atomic forces and charges were evaluated at the minimum in potential energy surface to predict the harmonic vibrational infrared spectrum using density functional perturbation theory [61,62].
After a full geometry optimization of all CH 4 , CO 2 , O 3 , SF 6 , NF 3 , CF 4 , N 2 O, SO 2 , 4NP, water, and copper squarate, the most stable adsorption sites of the pollutants inside the solid material were determined in order to minimize the first adsorption energy.The adsorption sites and mechanisms within the copper squarate were determined by Monte Carlo simulations [63], implemented in the Adsorption Locator module using COMPASS III [46] force field in Materials Studio 22.0 [42] code.Simulations involve a fraction of 10,000 steps following a number of cycles of 3 with 15,000 steps per cycle.Subsequent to the initial results, Density Functional Theory (DFT) calculations were conducted to optimize performance by adjusting the generic force fields within the Adsorption Locator module.More specifically, it was already proved that force fields such as COMPASS III [46] and Dreiding [64] tend to overestimate the number of molecules adsorbed in simulations [65], as was demonstrated in previous studies involving ZIF-4 or ZIF-6 [32,45].
The adsorption process occurs through chemical or physical binding, facilitated by the attractive forces between the solid cavity of the copper squarate and the pollutants, in both dry and humid conditions.The adsorption energies were calculated according to equations [66]: In these equations, ∆E ads refers to the isolated adsorption energy for each individual species (X) and ∆E mult ads denotes the mixed adsorption energy for species (X, Y, • • • . ..).E T stands for the total energy of copper squarate (C-S) with gas molecules included, E C−S corresponds to the energies of pure C-S, and E X is the energy of the respective pollutants.In addition, n, m, t represent the quantities of loaded molecules.The assessment of pollutants and water is conducted using a 10 Å supercell.A negative adsorption energy implies favorable adsorption.
The interaction energies between the co-adsorbed pollutants on the C-S material were calculated with the help of the energetic database obtained from DFT calculations using equation [45]: where E (n)−pollutants+Bulk represents the total energy of the solid material after adsorption of the n pollutants by copper squarate bulk, and ∑ E n−pollutants−Bulk is the sum of the total energy of the previously adsorbed species calculated separately.
The Grand Canonical Monte Carlo (GCMC) method [67] implemented in the Sorption module of Materials Studio 22.0 [42] was used to determine adsorption isotherms with the Metropolis algorithm and COMPASS III.The gas molecules were randomly loaded at a fixed temperature of 298 K into the fixed-volume copper squarate structure based on the gas fugacity.The simulations involve 100,000 equilibration steps followed by 1,000,000 Monte Carlo steps, assuming periodic boundary conditions.Long-range interactions were treated using the Ewald and group summation method with an Ewald accuracy of 0.0001 kcal/mol and a cut-off distance of 15.5 Å.The results are presented as a line graph (N = f (µ)), where f is the density of the number of loaded molecules N, expressed as a function of chemical potential.

Conclusions
In this work, we focus our interest on a new kind of MOF, namely copper squarate, whose adsorption behavior has not been explored in previous studies.For this reason, the adsorption of a series of small pollutants and some aromatic compounds was investigated using DFT calculations in a dry and humid environment.As a first observation, we note that small pollutants such as CH 4 , CO 2 , O 3 , SF 6 , NF 3 , CF 4 , N 2 O, and SO 2 favor the adsorption by different amounts depending on their chemical composition.
Squarate demonstrates excellent stability in the presence of the considered pollutants, exhibiting a favorable capture of the latter that results in a negative adsorption energy value.This phenomenon can be attributed to the existence of electrostatic interaction forces between pollutants and squarate.These interactions are particularly pronounced at short distances, occurring between the oxygen of squarate and the carbon atoms, nitrogen, and sulfur in the molecules CO 2 , N 2 O, SO 2 , and H 2 O, respectively.The presence of other weak parallel interactions further reinforces the binding of these molecules of the C-S.Additionally, the non-polarity of CH 4 gas makes it a favorable insert within squarate, accentuated by a hydrogen bond.Finally, the dipole nature of O 3 imparts characteristic properties, enabling it to react very selectively as an electrophile.These characteristic properties enable them to react very selectively as an electrophile.This characteristic facilitates its strong adsorption with squarate, driven by the intermolecular forces of Van der Waals and London dispersion.
Pollutants with higher volume and surface density are adsorbed in smaller quantities; for example, molecules with fluorine o nitrogen atoms.Indeed, adsorption in NF 3 , CF 4 , and SF 6 molecules is influenced by the number of fluorine atoms used and the adapted geometry.The smaller this number, the more the adsorption process is promoted.Furthermore, the more it adopts a complex structure, the more the interactions between pollutants and squarate are minimized.In contrast to all previous gases, nitrophenols are not able to be adsorbed in this bulk because of the small pore size.
The adsorption energy of the molecule containing nitrogen as N 2 O in a dry environment is revelaed to be the highest one (−63.02),even if the maximal insertion is one molecule per cell compared to four insert molecules of CH 4 with a maximum value of −23.02.This indicates the strong interaction involved in the first adsorption.Globally, the addition of water does not interfere with the adsorption, keeping the adsorption energy negative.
The co-adoption of CH 4 , SO 2 , and CO 2 with H 2 O implies an increase in adsorption energy.The same trend is observed for molecules containing fluorine atoms, such as NF 3 , CF 4 , and SF 6 .This impact on adsorption energy values is attributed to multiple interactions caused by the reactive sites of these pollutants, as well as the formation of weak bonds between water molecules.
In a humid environment, the same order is found for O 3 and N 2 O molecules as in the dry state, with a very slight decrease in adsorption energies.Meanwhile, 4NP is an exception to the other pollutants studied.On the other hand, it is characterized by an unfavorable capture on C-S in both humid and dry environments, resulting in positive adsorption energy.This is attributed to its large size and the unavailability of sufficient space in the pore.However, C-S exhibits a high capacity to insert water inside its cavity, revealing its hydrophilic nature.This behavior stands in contrast to enormous MOFs well-documented in the literature for their highly hydrophobic characteristics.Such hydrophilicity is particularly advantageous, especially considering that gas molecules are typically present in humid conditions, allowing us to regard C-S as a revolutionary MOF for adsorption, making it a novel target in this field.
The number of molecules fixed in the C-S decreases from the dry to the wet state, indicating that water molecules compete with pollutants and form strong hydrogen bonds with the squarate, compared to the interactions between pollutants and the bulk.
In future work, it is recommended to simulate a molecular work with a reduced part of copper squarate to understand deeper what is happening inside the cell cavity and comprehend with better precision to adsorption nature and the type of electrostatic interaction, including Van der Waals or hydrogen bonds in both dry and humid environments, as well as to determine the factors that govern this adsorption phenomenon.In addition, to complete this, a molecular dynamics (MD) trajectory will be helpful to understand dynamic adsorption.

Figure 4 .
Figure 4. Sorption isotherms of single and mixture adsorbed pollutants with and without water at 298 K.

Figure 4 .
Figure 4. Sorption isotherms of single and mixture adsorbed pollutants with and without water at 298 K.

Table 2 .
Connolly free volume, surface area of supercell sizes, and density of copper squarate.

Table 2 .
Connolly free volume, surface area of supercell sizes, and density of copper squarate.

Table 3 .
The structural parameters of copper squarate containing a co-adsorbed mixture of a pollutant with water.

Table 4 .
Adsorption and co-adsorption energies in copper squarate calculated with CASTEP.

Table 4 .
Cont.The results indicate favorable adsorption of CO 2 , N 2 O, SO 2 , and H 2 O molecules into C-S using force field loading, with up to five, one, one, and eighteen species per double unit cell, respectively.The accuracy of these findings was validated by CASTEP calculations, revealing adsorption energies of −18.93, −63.32, −38.04, and −51.16KJ•mol −1 , respectively, for a single captured molecule.This suggests the existence of electrostatic interactions, as illustrated in the zoomed images in Figure