1. Introduction
Nowadays, there is an urgent problem arising from the continuous increase in global water scarcity [
1,
2,
3]. Based on the United Nations report, over 2.3 billion people are facing significant water stress [
4]. An essential solution capable of providing the required reliable life support is desalination—the process of extracting salt and minerals from seawater [
5]. However, the pursuit of sustainable and economically feasible freshwater extraction methods has become paramount. Freshwater scarcity poses profound challenges, particularly in coastal regions where seawater is abundant. Traditional widespread desalination methods like reverse osmosis or distillation, dominant in arid regions such as the United Arab Emirates, Saudi Arabia, and Israel, pose significant drawbacks due to their substantial energy requirements, including both economic and environmental loads [
6]. Investigating innovative approaches is crucial in addressing these challenges. Solar water evaporation emerges as a promising technology capable of purifying water and desalinating seawater [
7]. Unlike conventional methods, solar water evaporation harnesses renewable energy from the sun, ensuring sustainability while minimizing greenhouse gas emissions [
8].
Recent advancements in materials science have spurred the development of photothermal materials capable of efficiently harvesting the entire solar spectrum and converting it into heat [
9]. The technology of interfacial solar evaporation involves placing photothermal materials at the water–air interface to directly convert solar energy into vapor [
10]. To enhance efficiency, interfacial solar steam generation utilizes a porous thermal barrier with low thermal conductivity, effectively isolating solar absorbers from bulk water to minimize heat loss. Furthermore, porous structures or water channels serve as capillaries, facilitating water transport within solar absorbers.
The efficiency of solar evaporation devices significantly depends on the solar absorber (photothermal materials), substrate (the material that provides water transport from bulk water to absorber), and thermal insulation. The request for photothermal materials that possess both excellent broad solar absorption and desirable hydrophilicity has long been a challenge in water purification technology. The organic-based structures exhibit promising properties post-carbonization and are known for their excellent hydrophilicity [
11,
12,
13]. The current literature presents a noticeable gap in broad inorganic-based intrinsic solar absorbers with high hydrophilicity, making it difficult to find materials that seamlessly combine both attributes.
A promising route to address the global water crisis lies in the investigation of hydrates-containing compositions for solar-driven water evaporation [
14,
15]. Hydrate compounds possess crystal lattices that contain water molecules, making them suitable for application as solar water absorbers. Integrating the hydrates into solar water evaporators could offer an effective solution, harnessing the power of sunlight to facilitate water evaporation and purification, thus potentially alleviating water scarcity issues on a global scale. The combination of efficiency, affordability, and durability makes hydrates a compelling solution for sustainable water desalination. The efficiency of hydrate compositions is notable because they can absorb significant amounts of water molecules, which enhances the evaporation process with minimal solar energy input. The economic viability for large-scale deployment of hydrate compositions as absorbers is based on the low production costs, which not only maximize water output but also underscore the cost-effectiveness. Furthermore, the durability to resist harsh environmental conditions enhances their appeal for widespread implications, offering long-term solutions to water shortage challenges.
Recently, we demonstrated that transition-metal selenite hydrates (earth-abundant metals Ni and Co) possess high solar absorbance (>96%) in the solar spectral range, as well as excellent hydrophilicity, which plays a key role in water transport in the solar steam generation [
16]. The solar water evaporation rate of compositions deposited on the glass–wool membrane achieved 1.83 kg∙m
−2∙h
−1 and 2.34 kg∙m
−2∙h
−1 for nickel selenite hydrate and cobalt selenite hydrate, exceeding the theoretical limit of 1.47 kg∙m
−2∙h
−1 for 2D evaporators (see
Supplementary Note S1) [
17]. The hydrophilic properties observed in selenite hydrates arise from the presence of “water of hydration” trapped within its crystal lattice.
Note that for practical applications of photothermal materials, the materials will be subjected to exposure to sun irradiation and experience a relatively high temperature (~30 °C –> 100 °C) [
18,
19,
20]. In particular, in the case of dry conditions and/or continuous sun irradiation, a higher temperature will be obtained, which may induce further compositional/structural alterations to the materials. Accordingly, the clarification of the thermal stability of the constituent components and possible defect-related chemical reactions, as well as the examination of the applicable threshold temperature of the materials, will be crucial for applications. In this research, we attempted to investigate the intricate mechanisms governing water transportation and to uncover the factors responsible for the difference in evaporation rates between the two selenite hydrate photothermal materials. The thermal stability of the synthesized nickel and cobalt selenite hydrates was first analyzed by thermogravimetric analysis (TGA) to examine the possible compositional variation. Based on the observed compositional mass loss, a series of thermal treatments were performed at several different temperatures, and X-ray diffraction (XRD) and Raman spectroscopy were applied to assess their temperature-induced structural variations. The results reveal distinct temperature-dependent alterations of phase and structure for nickel and cobalt selenite hydrates, providing valuable insights into the factors influencing water transportation and evaporation rates.
3. Results and Discussion
To understand the thermal stability of the samples, TGA analysis was first performed in a relatively large range from RT to 800 °C, as the results presented in
Figure 1a,b, showing relatively similar thermograms of the decomposition process of NiSe and CoSe. To identify the number of decomposition steps, the derivative of each thermogram (DTG, %∙°C
−1) is also presented.
As can be seen from the profiles of the TG curves, both samples undergo several steps of mass reduction, as confirmed by the presence of different DTG peaks. In the low-temperature region, the mass of the samples shows a slight but gradual decrease with increasing temperature up to around 180 °C and 140 °C for NiSe and CoSe, respectively. DTG analysis reveals that the first and second steps are consecutive and overlapped: In the DTG curves, both NiSe and CoSe exhibit a peak at around 80 °C, while the positions of the second peak are different, located at around 160 °C and 120 °C, respectively. Increasing the temperature up to 330–360 °C results in a weight loss of 21.87% and 16.73% for NiSe and CoSe, respectively, along with the appearance of a strong DTG peak at 320 °C. At 800 °C, a total weight loss of 61.52% and 57.38% can be detected for NiSe and CoSe, respectively. In the temperature range between 440 and 620 °C, the most significant mass reduction in the TG curves and the appearance of the most intense DTG peak can be observed, which indicates a weight loss of around 30% in this area, probably owing to the formation of NiSeO
3 and CoSeO
3 [
21,
22].
As can also be seen from
Figure 1c,d, CoSe showed better stability than NiSe during both low-temperature treatments. The weight loss of CoSe during heating from 30 to 40 °C was estimated to be 0.337 %, and after 30 min, the final loss achieved 0.743 %, while the weight loss of NiSe amounted to 1.016%. A similar trend occurred at 70 °C, the CoSe powder lost 2.032% of weight, while NiSe lost 2.406%. The possible chemical and physical transformations that cause the observed decomposition steps will be discussed later.
Further analyses of the compositional/structural alteration of the samples induced by thermal treatments were conducted by XRD.
Figure 2 depicts a comparison of the XRD patterns acquired from freshly synthesized NiSe and heat-treated samples. According to the XRD spectrum for the untreated NiSe powder obtained from the aqueous synthesis at RT (
Figure 2a), the most prominent peaks are observed at 2
θ ≈ 29.71° and 15.64°, which mainly correspond to the (101) plane of the selenium Se phase (PDF #06-0362), and the (011) plane of Ni(SeO
3)(H
2O)
2 (PDF#49-0136), respectively. The predominant phase within the sample is Ni(SeO
3)(H
2O)
2, notwithstanding the observation that the peak at 2
θ ≈ 29.71° exhibits the highest intensity. This heightened intensity arises from the overlapping of the peaks originating from both the Ni(SeO
3)(H
2O)
2 phase and elemental Se, resulting in an overall enhanced peak intensity. Furthermore, a minor fraction of NiSe
2 in the Penroseite structure (PDF#41-1495) might also exist, marked by the primary peak at 2
θ ≈ 33.41°. From the data presented in
Figure 2, the phase composition of the heat-treated NiSe powders remains barely altered at temperatures ≤ 100 °C. This observation underscores the thermal stability of the crystalline phases within the NiSe material under these relatively low-temperature conditions. At an elevated temperature of 200 °C, the primary Se phase and the minor NiSe
2 component exhibited consistent phase characteristics, indicating their stability under these conditions. Concurrently, a notable reduction in the intensity of the Ni(SeO
3)(H
2O)
2 peaks was observed, signifying an indicative measure of the reduction in the hydrated phase and reflecting a thermally induced transformation or decomposition of Ni(SeO
3)(H
2O)
2. This is consistent with the presence of the peak at 160 °C in the DTG curve (
Figure 1a). As depicted in
Figure 2b, upon thermal treatment at 400 °C for 30 min, a pronounced alteration in the phase and structural constitution of NiSe becomes evident. The NiSe sample treated at 400 °C is characterized by the coexistence of NiSe
2 in two distinct structural configurations, Penroseite (PDF#41-1495) and Kullerudite (PDF#18-0886), with approximately equimolar proportions in its composition. This transformation highlights the substantial thermally induced metamorphosis in the crystalline phases and structural attributes of NiSe. In addition, no more signals from the elemental Se phase can be observed, as the melting point of Se is 221 °C. Therefore, both the decomposition of dehydrated Ni(SeO
3)(H
2O)
2 and the evaporation/reaction of Se can contribute to the observed mass loss in the range from 200 to 360 °C and the appearance of the DTG peak at 320 °C shown in
Figure 1a. Annealing of the samples at higher temperatures (>400 °C) may result in a further loss of Se and the formation of nonstoichimetric NiSe
x phases, responsible for the mass reduction in the TG curve (
Figure 1a).
Initially, the XRD profiles of the specimens subjected to low-temperature treatments (T ≤ 100 °C) displayed a resemblance, necessitating a more intricate examination of the peaks associated with distinct phases using LabSpec 5 software to perform spectral deconvolution of the specific peaks. Note that the position and width of an XRD peak are associated with the crystalline structure and crystallinity, as well as the presence of defects and nonstoichiometric alterations to distort the lattice of the material. The selection of specific 2θ values for peak analysis was motivated by the imperative to mitigate potential overlaps among multiple phases and the aspiration to discern these peaks distinctly, thus enhancing the precision and quality of the analytical outcomes.
Figure 3a shows the related spectra of the selected peaks representative for NiSeO
3(H
2O)
2 (15.6° and 16.8°, (011) and (10-1) planes, respectively), and
Figure 3b–d show a comprehensive presentation of spectral deconvolution results, specifically providing detailed data regarding the position, relative intensity, and full width at half maximum as a function of temperature, taking advantages of the selected peaks. These data points provide a thorough and precise depiction of the distinctive characteristics of different heating temperatures.
Note that in the graph for the peak position (
Figure 3b), the lines provided represent the 2
θ position of the peak from the database, and the error bar was estimated based on several XRD measurements of each powder. As the temperature increases from RT to 70 °C, the selected 2
θ peaks show a gradual reduction to lower values but then an increase as the temperature further rises up to 200 °C. In the meanwhile, the FWHM values first increase to 70 °C then decrease at higher temperatures. It is worth noting that in this study, two distinct peaks originating from the same components were used for confirmation, and the observed similar trend of variations in the plot demonstrates the validity of the obtained results. Concerning the intensity variation, the strongest peak at 15.6° of NiSeO
3(H
2O)
2 for the sample thermally treated at 70 °C was used as a normalization factor (the maximum intensity) for comparisons. In general, the relative intensity of NiSeO
3(H
2O)
2 also shows a first increase to 70 °C and then a decrease as temperature increases (
Figure 3c). With respect to the strong peak at 15.6° that corresponds to the (011) plane, the peak at 16.9° shows a relatively small intensity variation, which might indicate a slight structural alteration of NiSeO
3(H
2O)
2 upon thermal treatment.
In general, the as-prepared NiSe sample may have a certain amount of adsorbed water on the particle surfaces because of the hydrophilicity and porous nature of the material, besides the presence of hydrated water in the particles. Consequently, upon thermal treatment or sun irradiation, the increase in surface temperature results in the process of desorption and evaporation of the adsorbed water at relatively lower temperatures. That is to say, at lower temperatures, typically below 100 °C, the material undergoes a remarkable transformation in its surface lattice structure attributed to the removal of surface-bound water molecules, leading to dehydration. This can be supported by the observation of a mass reduction in the TGA curve and the presence of the DTG peak at around 80 °C (cf.
Figure 1a). The dehydration process causes a peak shift and an increase in FWHM of NiSeO
3(H
2O)
2.
As the temperature continues to rise, a distinct shift occurs in the dominant process. At higher temperatures, the removal of internally bound hydrated water becomes increasingly prevalent (cf. strong DTG peak at around 160 °C in
Figure 1a). This transition triggers an intriguing internal transportation of hydrated water within the material. The removal of the hydrated water from the lattice in NiSeO
3(H
2O)
2 can increase the structural order in NiSeO
3, and thus, a higher crystallinity can be expected, showing a decrease in the FWHM of the diffraction peak at higher temperatures. In addition, the 2
θ position will shift to higher values due to a dehydration-induced structural contraction.
Similarly, in the case of untreated CoSe, as illustrated by the XRD spectrum of CoSe RT in
Figure 4a, noteworthy peaks manifest at 2
θ ≈ 29.71°, corresponding to the (101) plane of the selenium Se phase (PDF#06-0362) and at 2
θ ≈ 29.61° attributed to the (012) plane of Co(SeO
3)(H
2O)
2 phase (PDF#25-0125). The heightened intensity of the strongest sample peak arises from the superposition of signals stemming from these two phases.
Figure 4a also presents data demonstrating the stability of the phase composition of heat-treated CoSe powders, with negligible variations observed at temperatures <200 °C and a minor reduction in Co(SeO
3)(H
2O)
2 content at 200 °C. Upon elevating the temperature to 400 °C, as shown in
Figure 4b, substantial phase transformations take place as the signals from Se and Co(SeO
3)(H
2O)
2 disappear. The primary phase is now identified as CoSeO
3 (PDF#47-0903), marked by a dominant peak at 2
θ ≈ 31.65°, and CoSe
2 in the Trogtalite structure (PDF#09-0234), characterized by a principal peak at 2
θ ≈ 37.62°. This alteration underscores the significant annealing-induced shift in the crystalline phases within the CoSe material at temperatures higher than 200 °C.
Comparing the XRD results of CoSe and NiSe, the presence of a significant amount of Co(SeO
3)(H
2O)
2 at 200 °C (cf. strong peak at 15.6° in
Figure 4a) indicates less dehydration of the hydrated water in CoSe. The formation of CoSe
2 from Co(SeO
3)(H
2O)
2 appears only in the range 200–400 °C and the structure of oxyselenide is kept even at 400 °C, while the NiSe sample annealed at 400 °C has completely converted to NiSe
2. The difference in the behavior of phase transition between NiSe and CoSe under identical reaction environments can be explained by different energy barriers. For the transition from Co
2+ (Co(SeO
3)(H
2O)
2) to Co
4+ (CoSe
2), the standard reduction potential reaches +3.34 V, and for Se
0, it is a challenge to capture two electrons from Co
2+. The reduction potential in the case of Ni
2+ → Ni
4+ is only +1.59 V, which explains the easier formation of NiSe
2 [
23].
The formation of NiSe
2 in NiSe and the suppression of CoSe
2 in CoSe can be described based on ligand field theory. This theory suggests that ligands can be perceived as negative charges arranged in specific geometric patterns—octahedral or tetrahedral. In an octahedral crystal field, the 3
d orbitals of the transition metal ion are split into two groups, with the energy difference described by the Δoct crystal-field splitting parameter: t
2g orbital symmetry with 0.4Δoct, and e
g orbital symmetry with 0.6Δoct.
Figure 5 shows an illustration of the 3
d orbital configuration of nickel and cobalt ions in selenite hydrate complex: Co
2+ and Ni
2+ have spin states with 3
d orbital configuration t
2g6e
g1 and t
2g6e
g2, respectively. In NiSe
2, the Ni
4+ owns the 3
d orbital configuration as t
2g6e
g0 because Ni
2+ is easier to donate 2 electrons from the e
g orbitals for becoming NiSe
2 rather than for Co
2+ for the formation of CoSe
2, which needs additional energy to break the electron pairs in the t
2g orbitals. Therefore, the way for the oxidation of Ni
2+ to Ni
4+ is more favorable than that of cobalt. In addition, the presence of asymmetrically occupied
d-orbitals (e
g1) in Co(SeO
3)(H
2O)
2 may cause a distortion of the octahedral structure due to the Jahn-Teller effect. This structural alteration is assumed to partly stabilize or destabilize the hydrated water along different directions in CoSe, so the structural asymmetry may result in the appearance of two peaks at 120 and 250 °C in the DTG curve of CoSe (
Figure 1b), with respect to the single DTG peak at 160 °C for NiSe. Accordingly, although a relatively minor fraction of dehydration occurs upon annealing at 200 °C, the mobility of hydrated water in CoSe can be more efficient in certain directions (cf. the lower DTG peak position).
Figure 6 provides more XRD results pertaining to the CoSe sample, specifically, the enlarged spectra, the 2
θ position, relative intensity, and FWHM for Co(SeO
3)(H
2O)
2. At temperatures below 200 °C, unlike NiSe, with increasing temperature, the 2
θ peaks slightly shift to higher values along with a reduction in FWHM, indicating a surface lattice reconstruction. In addition, the peak at 16.7° (−101 plane) shows an insignificant variation in the relative intensity compared with the strong peak at 15.5° (−110 plane). This might be owing to the competition between the desorption of surface water and the movement/dehydration of hydrated water in certain crystallographic directions induced by asymmetric lattice distortions and might also be associated with the capability of internal hydrated water transport within the material at low temperatures. At higher temperatures (200 °C), although Co(SeO
3)(H
2O)
2 exhibits greater stability, retaining its structural integrity, the movement of internal hydrated water along different directions becomes the predominant process, which might cause a higher degree of structural disorder and result in the observed peak broadening.
Figure 7 shows the variations in the crystallite size and the spacing of the crystal layers
d with temperature for Ni(SeO
3)(H
2O)
2 and Co(SeO
3)(H
2O)
2. The determination of crystallite (grain) size
D was conducted by applying the Scherrer equation as outlined in Equation (1) [
24]:
where
K is the numerical factor (Scherrer constant),
λ is the CuK
α wavelength (0.154056 nm), and
β is the full width at half maximum (in radians) of the (hkl) peaks at the diffraction angle 2
θ. The FWHM of the diffraction peak was calculated by spectral fitting using a Gaussian function.
K is a numerical factor referred to as the crystallite-shape factor. Its value is contingent on the crystallite shape and the definitions of the average crystallite size and width. Within the range of 0.87–1.0,
K’s numerical value can vary significantly, although it is commonly accepted as 0.9 [
25].
And the spacing of the crystal layers
d (Å) was calculated via Bragg’s law (Equation (2)) [
26]:
where
n is the diffraction order (1, 2, 3…),
λ is the wavelength (Å), and
θ is the glancing angle (radian).
With increasing temperature <100 °C, Ni(SeO
3)(H
2O)
2 shows a reduction in the crystallite size owing to the removal/dehydration of surface water, while Co(SeO
3)(H
2O)
2 presents a lattice expansion (reduction in
d) owing to the movement of internal hydrated water in certain crystallographic directions. At higher temperatures, the occurrence of dehydration/movement of internal bound hydrated water Ni(SeO
3)(H
2O)
2 may result in a lattice expansion, while Co(SeO
3)(H
2O)
2 presents a lattice reduction due to the competition of the dehydration of the internally bound hydrated water.
Figure 7c,f show a comparison of the intensity ratios of the peaks of the components for NiSe and CoSe in response to the change in annealing temperature. As can be seen, the NiSe sample shows a slight increase in the ratio,
I(Se)/
I(Ni(SeO
3)(H
2O)
2) in the low-temperature region (<100 °C), and then a significant increase at 200 °C, while for Co(SeO
3)(H
2O)
2, the intensity ratio barely changes at low temperatures, and then slightly increases at 200 °C, confirming minor dehydration.
To further examine the annealing-induced structural alterations, micro-Raman spectroscopic analysis was performed on the samples, as their typical Raman spectra are shown in
Figure 8. For untreated NiSe and CoSe samples, the Raman spectra exhibit similar profiles with the presence of two main peaks despite a difference in peak intensities (
Figure 8a,d). At RT, one main peak originated from the symmetric stretching of (SeO
3)
2− in NiSeO
3(H
2O)
2 can be seen at about 823.5 cm
−1 and at 819.8 cm
−1 for CoSe [
27,
28]. In NiSe, the low-intensive peak of the Se-Se stretching mode is observed at 217 cm
−1, which was close to the value observed in other diselenides [
29,
30]. The stretching of the Se–Se bond in metallic Se is associated with the second main peak at around 244.9 cm
−1 for NiSe and at 244.8 cm
−1 for CoSe [
31,
32].
Upon heating, the spectral morphology shows a slight change in the low-temperature range for both samples, while a marked decrease in the intensity of the Raman bands was found for NiSe but not for CoSe upon annealing at 200 °C. Further annealing at 400 °C results in the appearance of Raman signals from NiSe
2 (~203 cm
−1) or CoSe
2 and Co(SeO
3) (~185 and 201 cm
−1) and the disappearance of the signals of Se and Ni(SeO
3)(H
2O)
2 (Co(SeO
3)(H
2O)
2), in good agreement with the XRD result.
Figure 8b,d–f show an enlargement of the variation in the SeO
32− peak at around 820 cm
−1 for the samples. As can be seen, a clear shift of the peak position toward lower wavenumbers after annealing could be observed (
Figure 8g), which is taken to be associated with the surface dehydration/desorption in NiSe and also the internal movement of hydrated water in CoSe. These Raman spectroscopic observations are in line with the XRD results. The perturbations of XRD and Raman spectra not only underscore the thermally induced changes in composition but also provide further insights into the structural alterations associated with temperature variations, offering a nuanced understanding of the dynamic behavior of NiSe and CoSe under thermal conditions.