Based on our aim to investigate the functionalised carbon-cation interfacial water hydrogen-bond network, we prepared GO/CaCl2 aerogel via solution intercalation method36. The synthesis procedure is shown in Fig.S1a and described in detail in the methods section. The as synthesized samples are marked as GOA for graphene oxide aerogel and Ca-GOA for Ca2+ intercalated GO aerogel respectively. Exemplary images of GOA and Ca-GOA are displayed in Fig. S1b-d.
We characterised the structure and chemical properties of the samples before and after intercalation. Scanning electron microscopy (SEM) with energy dispersive spectroscopy (EDS) images are shown in Fig. 1a-c and Fig. S2b-c (see supplementary note 2), illustrating the sample morphology. The synthesized Ca-GOA samples display a porous structure similar to the GOA sample with typical wrinkled large surface area. X-ray photoelectron spectrometry (XPS) reveals that the carbon to oxygen ratio (C/O) of GO remains constant at ~ 2 before and after intercalation (Fig. 1d). C1s XPS spectra curves further suggest that Ca-GOA (Fig. 1f) have similar carbon-carbon and carbon-oxygen bonds composition compared to GO samples (see Fig.S3a-d in supplementary note 3). The XPS survey and atomic percentage of each element for both GO and Ca-GOA are shown in Fig. 1f. Around 5.1% of calcium atoms are detected on the Ca-GOA surface, while in the original GO, no calcium was detected. This confirms the successful intercalation of Ca-ions into the Ca-GOA samples. It is noted that a small number of chlorine and sulfur atoms are also detected on Ca-GOA however are not the main concern in this study.
To investigate the hydrogen-bond network in the Ca-GOA, we performed AWC measurements as described in the methods section and compared the results to literature values from our own studies and others35,37. Figure 2a-c show the recorded and reported adsorption isotherms. The water uptake (m/m0) of Ca-GOA was recorded with varying relative humidity (RH%) at constant room temperature (298K). The relative humidity is then converted into the partial pressure of water vapor to visualise the isotherm. The results show that AWC ability of Ca-GOA samples is significantly improved to up to 2.1 g/g from pure GO samples (GOM) with 0.5 g/g, which indicates the enriched hydrogen network on Ca-GOA surface. It should be noted that the referred GO adsorption curve was measured for GO membranes (GOM) and not in the form of aerogel. However, we ruled out the influence of materials morphology for AWC by showing that GOM and pure GOA have similar AWC ability (supplementary note 3).
Further analysis of the isotherms allows us to get a deeper understanding of interfacial water molecules during the AWC. All three different materials, CaCl235, GOM37 and our measured data of Ca-GOA, fit well under the BET model38,39. However, the fitted BET equations are notably different from each (see supplementary note 4 for fitting details). The degree of the polynomial function of the BET equation, represents the number of adsorbed water molecule layers (n) on the surface of the material35,39,40. Here, n = 7 for CaCl2, n = 4 for GOM and n = 6 for Ca-GOA. Hence, the intercalation of Ca into Ca-GOA results in a stronger hydrogen network on surface compared to pure GO. In other words, the intercalation of cation enhances the interfacial water hydrogen network at the oxidized carbon surface, however the mechanism of such enhancement is still unclear.
Based on the BET model fitting analysis, we established that the Ca-intercalation offers strong adsorption sites for water molecules on the Ca-GOA surface. If that is the case, the adsorption isotherm should follow the Freundlich model, as it represents a heterogeneous surface adsorption process41. In particular, Freundlich model describes a material with special adsorbing sites which are heterogeneously distributed on the adsorbent surface42–44. In this case, the Ca atoms are the strong water-attracting adsorbing sites, well-dispersed on the GO plane, allowing coverage of the whole surface with several layers of water molecules45. As shown in Fig. 2d-f, the water adsorption isotherm of Ca-GOA can be fitted in high agreement with the typical Freundlich isotherm model, while being less in agreement with GO and CaCl2. The detailed information for Freundlich model fitting is shown in supplementary note 5.
The trend of water molecules per oxygen and calcium atom of GOM, GOA, Ca-GOA and CaCl2 versus the water partial pressure are shown in Fig. 2g-h. The uptake of water molecules per oxygen atom on Ca-GOA is up to 8 times higher than that of GOM and GOA at all ambient environments. Surprisingly, the water uptake per calcium ion, is up to 2.7 times higher in Ca-GOA compared to CaCl2 at the highest water partial pressure (Fig. 2i).
Here, we can conclude that the water uptake per calcium and oxygen atom in Ca-GOA is much higher than in the individual materials, CaCl2 and GOA/GOM. Keeping in mind that the Ca-sites were identified as the main contributor to the water uptake in Ca-GOA, it is striking to see that the water uptake per calcium is 2.7 times higher in the Ca-GOA compared to CaCl2. This strongly suggests that this enhancement is linked to an interplay between the functional groups of GO and the intercalated cations to form a synergistically enhanced hydrogen-bond network. To investigate this hypothesis, we further performed density functional theory (DFT) calculations.
We examine the hydrogen bond properties between oxygen functionalities and water molecules with and without the existence of calcium ions using DFT calculation. We select a graphene plane with a bare epoxide group as a simplified model of a GO nanoflake surface as shown in Fig. 3a. Epoxide and hydroxide groups are typical oxygen functionalities on the GO basal plane and have been experimentally confirmed to have a strong hydrogen bonding interaction with water molecules32,46. The DFT calculations were performed at the PBE0-D3BJ/def2-QZVPP level of theory47–49. Figure 3b-c indicates the process of the epoxide hydrogen bonded to one and two water molecules. We obtain hydrogen bond distances and energies typical of moderate hydrogen bonds50–52. Namely, both water molecules hydrogen bonded to the epoxide group with a distance range of 1.96–1.99 Å and enthalpy of − 14.8 kJ/mol and − 13.4 kJ/mol for the first and second water molecule, respectively. However, when the hydrated calcium cation is bound to the epoxide group, the hydrogen bonding network presents significantly different properties. We note that when the epoxy oxygen is coordinated to the Ca cation, the oxygen is bound to the GO surface via one C–O bond, as illustrated in Fig. 3d-e. We note that the bonding situation in Fig. 3d-e also represents the bonding of a hydrated Ca cation to a C–O functional group on the GO surface53. Hereinafter, the oxygen connected to the GO surface with a single C–O bond will be referred to as the GO oxygen, rather than the epoxy oxygen.
With the presence of calcium cation, our calculations suggest that the GO oxygen is coordinate to the hydrated calcium cation. This dramatically enhances the hydrogen-bond network surrounding GO oxygen. Figure 3d presents a scenario when one water molecule is hydrogen bonded to both the GO oxygen and the hydrated calcium ion. Compared with the system in Fig. 3c, this water molecule is now hydrogen bonded via three hydrogen bonds, one with the GO oxygen and two with the hydrating water molecules around the calcium ion. As shown in Fig. 3d, the hydrogen bond formed with the GO oxygen (1.546 Å) is significantly shorter than the hydrogen bonds formed with the water molecules hydrating the calcium cation (1.699 and 1.725 Å). Furthermore, this hydrogen bond is also significantly shorter than that formed between the water molecule and bare epoxide group presented in Fig. 3b. Similarly, in Fig. 3e, we further show the scenario when the GO oxygen is hydrogen bonded to a second water molecule. The calculation shows that the length of the hydrogen bond between the GO oxygen and the second water molecule is also largely shortened, from 1.982 Å (Fig. 3c) to 1.689 Å (Fig. 3e).
From the hydrogen bond enthalpy point of view, our calculations show that the hydrogen bonding between the GO oxygen and water molecules are highly reinforced by the hydrated calcium ion. In the scenario of the bare epoxide group on the GO surface, the hydrogen bond enthalpy between the oxygen and water molecules is around − 14 kJ/mol, as mentioned above. However, in the presence of the hydrated calcium ion, the hydrogen bond enthalpy between the GO and water molecules in the system shown in Fig. 3d increases to as much − 66.7 kJ/mol and in Fig. 3e comes to − 52.9 kJ/mol. Thus, the H-bond enthalpy in the presence of the hydrated calcium ion is about 3–4 times higher than that of a bare system. Both the increased hydrogen bond enthalpy and shorter hydrogen bond lengths are attributed to both the larger hydrogen bond network and stronger hydrogen bonding acceptor strength of the calcium ion decorated GO oxygen.
We further investigate the hydrogen bonding acceptor properties of the GO oxygen based on the atomic polar tensor (APT) charges54,55. The results are shown in Table 1. We obtain the following APT charges of the epoxide oxygen, − 0.40 e, − 0.47 e, and − 0.51 e in the system shown in Fig. 3a-c, respectively. As expected, the negative charge on the oxygen increases with the number of hydrogen bonds it is involved in. The charge on the hydrogen and oxygen atoms of the water molecules remains relatively constant for systems 3b and 3c. Namely, they range between + 0.23 e and + 0.27 e for the hydrogen and between − 0.45 e and − 0.49 e for the oxygen. Coordination of the hydrated calcium cation to the GO oxygen results in a dramatic change in the oxygen charge. In particular, the coordination of the hydrated calcium cation increases the negative charge on the GO oxygen from − 0.40 e to as much as − 0.86 e. This significant increase in negative charge on the GO oxygen makes it a much stronger hydrogen bond acceptor. Accordingly, the H-bond distance with the water molecule is shortened from 1.965 Å (Fig. 3b) to 1.546 Å (Fig. 3d). We note that the later hydrogen bond distance represents an exceptionally short hydrogen bond for a HOH•••O system (i.e., a water molecule coordinated to an oxygen atom)50,51. It also reveals that coordination of calcium ion to the GO oxygen alters the atomic charges on the hydrogen and oxygen of the hydrogen-bonded water molecules. For example, for the systems depicted in Fig. 3b and 3d, the positive charge on the hydrogen increases from + 0.27 e to + 0.37 e and the negative charge on the oxygen increases from − 0.49 e to − 0.63 e (see Table 1). Remarkably, this demonstrates significant medium-range effects of the hydrated calcium cation on the charge of an oxygen centre to which it is bound via a hydrogen bond network (i.e., not covalently bound). The above results indicate that the hydrogen bonding ability of one epoxide group on the GO surface is enhanced by the coordination of the hydrated calcium ion. Such enhancement was shown via the increasing hydrogen bond acceptor strength of GO oxygen and additional hydrogen bonding interactions with the first hydration sphere of calcium cation. It is evident, both theoretically and experimentally, that epoxide groups on a GO surface tend to cluster in islands rather than be uniformly distributed across the surface26,56. This leads to a natural question, whether this single hydrated calcium cation can interact with more than one epoxide group on GO surface.
Table 1
Atomic polar tensor (APT) charge (q) in a.u. involved in the hydrogen bonds and on the Ca atom for the systems in Fig. 3.
Model
|
q (Owater)
|
q (Hwater)
|
q (Oepoxy)
|
q (Ca)
|
3a
|
N/A*
|
N/A
|
–0.40
|
N/A
|
3b
|
–0.49
|
+ 0.27
|
–0.47
|
N/A
|
3c
|
–0.48, − 0.45
|
+ 0.26, + 0.23
|
–0.51
|
N/A
|
3d
|
–0.63
|
+ 0.37
|
–0.84
|
+ 1.34
|
3e
|
–0.61, − 0.60
|
+ 0.40, + 0.35
|
–0.84
|
+ 1.25
|
* N/A corresponds to the absence of atomic polar tensor charge.
Table 2
Comparison between the hydrogen binding enthalpies at 298 K (∆H298,bind, in kJ/mol) for the systems in Fig. 3 obtained in the solid state and in bulk aqueous solution.
Model
|
∆H298, bind (kJ/mol)
|
Solid state
|
Aqueous solution
|
3a
|
N/A*
|
N/A
|
3b
|
–14.8
|
–6.8
|
3c
|
–13.4
|
–5.3
|
3d
|
–66.7
|
–42.7
|
3e
|
–52.9
|
–32.9
|
* N/A corresponds to the absence of hydrogen binding enthalpy.
We further investigate the system with two epoxides on a graphene plane coordinated to a single hydrated calcium ion (see supplementary note 7). Each of the two oxygens on the GO surface is able to hydrogen bond to two water molecules. For the first three water molecules, we obtain binding enthalpies that are similar to those obtained from the functionalized systems in Fig. 3, namely 70.2, 58.5, and 63.6 kJ/mol, respectively. For the fourth water molecule, we obtain a lower binding enthalpy of 29.9 kJ/mol; this reduction is partly attributed to two (rather than three) hydrogen bonds in which this water molecule is involved in (supplementary note 7). Importantly, all these binding enthalpies are significantly larger than those obtained from the systems with the absence of hydrated calcium cation (Fig. 3b-c).
With the results above, we can now further optimize the hydrogen bond enthalpies from the perspective of the experimental conditions. This is because the calcium intercalated GO surface in the experimental settings is expected to be intermediate between the solid state and a bulk aqueous solution. It is well-established that hydrogen bond strengths are influenced by the effect of the solvent. In particular, the H-bond strength with the GO surface decreases with the polarity of the medium in the order of solid-state > non-polar solvents > polar solvents. Thus, the calculated hydrogen bond enthalpies above, which do not include solvent corrections, are expected to represent the upper limits for the experimental setting. Therefore, it is instructive to calculate the hydrogen bond enthalpies in bulk aqueous solution to obtain lower limits for the hydrogen-bond enthalpies. For this purpose, we use the conductor-like polarizable continuum model (CPCM)57, which has been found to provide good performance for aqueous solution58,59. The solvation corrections reduce the hydrogen bond enthalpies for the unfunctionalized GO models by a factor of \(\sim\)2.5, whereas they reduce the bond enthalpies for the Ca-functionalized GO models by a factor of \(\sim\)1.5. As shown in Table 2, the inclusion of the solvation corrections widens the gap between the hydrogen bond enthalpies of the GO models with and without coordination of the calcium ion. Considering that the GO surface in the experimental settings is expected to be intermediate between the solid state and a bulk aqueous solution, this is strong evidence to explain our experimental observation.