INTRODUCTION
Photochromism is a very useful property for modulating solar radiation.1 Two types of photochromic materials have been widely investigated: organic materials and inorganic transition metal oxides. Although organic photochromic materials exhibit relatively large modulation and fast response, their lack of durability makes them less useful in practical applications. The transition metal oxides such as WO3,2 MoO3,3 Nb2O5,4 V2O55 and TiO26 exhibit photochromism by band gap irradiation. Their photochromic characteristics are so widely variable depending on crystal size, structure, morphology and impurities. Among the photochromic transition metal oxides, WO3 and MoO3 have been studied most extensively. The photochromic mechanisms of these oxides have been described in several different ways, such as color center model,7 double charge injection model8 and small polaron model.9 Despite considerable research, the most important factor determining the overall photochromic (PC) properties of the oxides is still unclear. This is partly due to the lack of uniform preparation method for comparison. Here, we prepared various sizes and crystal structures of nanocrystalline WO3 and MoO3 through simple acid precipitation and only subsequent heat treatment with the same type of precursor. And the crystallite-size and crystal structural dependence of their photochromic properties are comparatively investigated.
EXPERIMENTAL
Materials and Preparation
Ammonium paratungstate pentahydrate (APT, (NH4)10W12O·5H2O, 99%), was purchased from Nippon Chemical Co. Ammonium heptamolybdate tetrahydrate (AHM, (NH4)6Mo7O24·H2O, >98%) and ~66% HNO3 were obtained from Dae-Jung Chemical Co. Both hexagonal nanocrystalline WO3 and MoO3 were prepared by an acid precipitation method using APT or AHM and HNO3 as precursor and acid, respectively.10 For hexagonal WO (h-WO3) 35.0 g of APT was suspended in 250 mL H2O at 90 ℃ and 100 mL of ca. 0.03 M HNO3 solution was added slowly and the reaction mixture was vigorously stirred for two hours. The yellowish precipitate was filtered and washed three times with a copious amount of water. The resulting precipitate was dried at 80 ℃ overnight. This sample has hexagonal crystal structure and denoted as W- 1. In order to increase the crystallinity of hexagonal WO3 without structural transformation, W-1 was annealed at 350 ℃ in air for 3 h and the annealed sample was denoted as W-2. Since increasing crystallinity of W-1 by annealing is limited to the maximum temperature of about 370 ℃, hexagonal WO3 with further improved crystallinity was prepared by this way: the product mixture in the W-1 sample preparation was evaporated to dryness without filtering and washing and the remaining powder was calcined at 350 ℃ in air for 3 h. The resulting powder was washed with water and dried at 80 ℃ overnight. This sample was denoted as W-3. The monoclinic WO3 (γ-WO3) was obtained by thermal treatment of W-1 sample at 400 ℃ in air for 3 h and denoted as W-4. For orthorhombic WO3·H2O (denoted as W-5), the procedure was the same as for W-1, except that the concentration of HNO3 was about 0.1 M. Hexagonal MoO3 (h-MoO3) was prepared by the same way as in W-1 except that the precursor of Mo was AHM and the half amount of nitric acid was used. Orthorhombic MoO3 (α-MoO3) was obtained by annealing h-MoO3 at 500 ℃ in air for 3 h. The resulting samples were labeled Mo-1 and Mo-2 for the former and the latter, respectively.
Characterization
Powder X-ray diffraction (XRD) patterns were collected by Phillips X’Pert MPD diffractometer using CuKα radiation source at a scan rate of 0.02° s-1. Field emission electron microscopy (FE-SEM, Jeol JEM-6700F) was used to get the morphologies of the prepared WO3 and MoO3 powders. Diffuse reflectance spectra of the powder samples were measured by Shimadzu UV-2600 with integrating sphere attachment that has two detectors of PMT and PbSe. The photochromic color changes were performed by illuminating the powder samples with 300 W Xe-arc lamp (ILC Tech). The incident light power was measured about 300 mW/cm2 by power meter (LI Model 7100).
RESULTS AND DISCUSSION
Five different samples were prepared for WO3 and two samples were prepared for MoO3. WO3 is known to have various polymorphic phases such as tetragonal, orthorhombic, monoclinic, triclinic, cubic and hexagonal phases.11 Among them, monoclinic WO3, known as γ-WO3 is the most stable phase at ambient condition. Metastable hexagonal WO3 (h-WO3) has received particular attention for electrochromism12 and photochromism13 as well as γ-WO3. Here we prepared 3 different phases of WO3: hexagonal, monoclinic and orthorhombic hydrated ones for comparison. And three different crystallite sizes for h-WO3 were prepared.
Fig. 1 shows the XRD patterns of all samples of WO3 and MoO3. The WO3 samples labeled W-1 to W-3 show hexagonal patterns in good agreement with JCPDS 33-1387 data. The crystallite size was estimated by Scherrer equation using the highest intensity peak for all samples. W-1 obtained by acid precipitation alone had a crystallite size of about 5 nm. Annealing of W-1 at 350 ℃ in air increased the crystallite size to about 8 nm without structural transformation (W-2 sample). Since the crystal structure changes from hexagonal to monoclinic when annealing temperature increased higher, h -WO3 with crystallite larger than W-2 could not be prepared only by increasing annealing temperature further. Thus, as mentioned in experimental section, the largest crystalline h-WO3 (W-3) was obtained by thermal treatment of W-1 at 350 ℃ with the decomposable byproduct salt NH4NO3 in the precipitation reaction. This thermal treatment increased the crystallite size from 5 nm to 25 nm, whereas only simple annealing increased it from 5 nm to 8 nm at the same temperature. h-WO3 with good crystallinity has been reported by hydrothermal reaction at elevated temperature,12(b),14 however, the thermal treatment method used here for h-WO3 has never been reported previously to our knowledge. The reason for increasing the crystallinity is not clear at this point. Monoclinic WO3 (γ-WO3, W-4) was obtained by thermal transformation of h-WO3 by simple annealing at 400 ℃ in air. The crystallite size of W-4 was estimated to be ca. 29 nm. Orthorhombic hydrated WO3 (W-5) was prepared by the same precipitation reaction as W-1, but with a much larger amount of acid. Interestingly, only the amount of acid gave a dramatic effect on the crystal structure of the product WO3. The crystallinity of W-5 as shown in the lower part of Fig. 1(b) is much better than those of W-1 to W-4 shown in Fig. 1(a). The crystallite size of W-5 was estimated to be about 55 nm from (111) peak.
Figure 1. XRD patterns of the synthesized WO3 and MoO3 samples (Some major peaks are marked).
MoO3, which has the same octahedral motif as WO3 is known to have three main polymorphic structures: orthorhombic, monoclinic and hexagonal. Orthorhombic MoO3 (α-MoO3) is thermodynamically preferred and the other two phases are metastable. Hexagonal MoO3 (h-MoO3) labeled Mo-1 was prepared by the same acid precipitation method as h-WO3. XRD pattern of Mo-1 was shown in the upper part of Fig. 1(b) and coincident with JCPDS 47-0871. The crystallinity of Mo-1 is very high even without annealing, as noticeable from the intensity of XRD peaks. Although the preparation method of MoO3 is identical to that for WO3, there are significant differences in the results. The results of this precipitation reaction are quite interesting considering that good crystalline h-MoO3 was often obtained through hydrothermal reaction,15 because the precipitation reaction even at room temperature gave fairly good crystalline h-MoO3, although the results were not shown here. The estimated crystallite size of Mo-1 was about 49 nm. Orthorhombic MoO3 (Mo-2) was prepared by thermal conversion of Mo-1 at 500 ℃ in air. Partial conversion occurred at annealing temperature below 500 ℃, while complete conversion was achieved above 500 ℃. The XRD pattern of Mo-2 shown in the upper part of Fig. 1(b) shows pure orthorhombic phase coincident with JCPDS 75-0912.
Fig. 2 shows FE-SEM images of the synthesized WO3 and MoO3 samples, showing the particle morphology. Hexagonal WO3 samples with different crystallite sizes do not differ greatly from each other, but appear to be composed of agglomerates composed of fine nanoparticles, and appear to be slightly larger as the particles are fused after heat treatment. Monoclinic WO3 (W-4) prepared by annealing W-1 at 400 ℃ shows slightly larger aggregates without appreciable morphological change, too. The orthorhombic hydrated WO3 (W-5) exhibits a platy morphology and is significantly different and larger than the other WO3 samples, as expected from the peak intensities of the XRD pattern. Hexagonal MoO3 sample (Mo-1) is rod-shaped with a length of 5-10 μm, much larger than all WO3 samples. (See the difference between the scale bars in the figure) The orthorhombic MoO3 (Mo-2) prepared by sintering the Mo-1 sample is more plate-shaped than rod-shaped and has a smaller size than Mo-1.
Figure 2. FE-SEM images of the synthesized WO3 and MoO3 samples (×80,000).
Diffuse reflectance spectra were collected for all samples to obtain bandgap energy information for these samples. Fig. 3 shows the UV-Vis diffuse reflectance spectra of all samples as Kubelka-Munk function (F(R)). As shown in Fig. 3(a), the absorption edges of all the WO3 samples are different. The edge wavelength increased with increasing sample number. Since the crystallite size of WO3 samples increased in order of sample number, the band gap energy of all samples decreased as the crystallite size increased regardless of the crystal structure. The bandgap energy (Eg) of WO3, which is an indirect bandgap semiconductor,16 can be estimated by Tauc plot.17 From the extrapolation of the absorption edge in the plot of (αhν)1/2 or (F(R)hν)1/2 vs. hv (Fig. 4), the bandgap energies of all the samples were estimated and listed in Table 1 along with the unit cell structure and crystallite size. The bandgap energy values of all the WO3 samples systematically decreased from 3.4 to 2.5 eV with increasing crystallite size, regardless of crystal structure. W-5 sample, which is orthorhombic monohydrated WO3, had the smallest Eg as expected from its distinct yellow color. Although Eg value of bulk crystalline WO3 is not well reported for each crystal structure, some reported Eg values were 3.3,18 2.7619 and 2.75 eV20 for hexagonal, monoclinic and orthorhombic hydrated WO3, respectively. The values we obtained are in some agreement with these values. Bandgap energies of the MoO3 samples were estimated as 3.2 and 3.1 eV for Mo-1 and Mo-2, respectively. After the thermal conversion of h-MoO3 to α-MoO3, Eg decreased as the crystallite size increased regardless of crystal structure in the same way as WO3. The reported Eg values for h-MoO3 and α-MoO3 that were prepared in a method similar to ours were 3.05 and 3.18 eV, respectively.21 The result was in the opposite direction. The calculated Eg values for h-MoO3 and α-MoO3 were in line with our result.22
Figure 3. UV-Vis diffuse reflectance spectra of the synthesized WO3 and MoO3 powder samples.
Figure 4. Tauc plots for (a) WO3 and (b) MoO3 samples.
Table 1. Structural and photochromic properties of WO3 and MoO3 samples
When all these WO3 and MoO3 samples were irradiated with UV light, their color change was visually perceptible within a minute, although it varied from sample to sample. Fig. 5 shows the pictures of some samples showing distinct color change before and after UV irradiation. As shown in pictures, W-1 sample shows the most intense visual color change from pale yellow to deep blue under UV irradiation. Fig. 6 shows the diffuse reflectance spectral changes of all powder samples of WO3 under UV irradiation in an ambient condition. For h-WO3 samples with different crystallite sizes (W-1~W-3), The crystallite size dependence of the photochromic properties of these samples is clearly revealed. The smaller the size, the better the color change. W-3 has the smallest Eg among the h-WO3 samples and can absorb more light, but it appears to have the worst photochromic effect. Monoclinic WO3 (W-4), which is larger in size and has smaller Eg than all h-WO3 samples, showed a worse result than W-3 consistently in size dependence regardless of crystal structure. This size dependence does not appear to hold for the W-5 sample, which is the largest in crystallite size and has the smallest Eg among all WO3, as shown in Fig. 5 and 6. W-5 differs from other WO3 samples not only in its crystal structure but also in being hydrated. This color change of W-5 is very interesting because it has been reported that the orthorhombic WO3·H2O did not show photochromism.23 He and Zhao reported the color change of orthorhombic WO3·H2O after irradiation of NIR laser light.24 They attributed this color change to a photoinduced thermochromic effect accompanied by structural change from orthorhombic WO3·H2O to monoclinic WO3. Our result may be related to those of He and Zhao. Further investigation is required. Fig. 7 shows the time dependence of relative reflectance changes at 660 nm due to the photochromic change of all the WO3 samples. The color change of all samples except W-5 were saturated in a few min. W-5 sample showed slightly slower color change than other WO3 samples. This fact implies that the photochromic mechanism of the W-5 sample may be different from that of other samples as mentioned earlier. The initial colorization efficiencies of all samples at 660 nm were estimated with the relative diffuse reflectance and light intensity using eq. (1).23(a)
Figure 5. Photochromic color changes of the hexagonal (W-1), monoclinic WO3 (W-4), orthorhombic WO3·H2O (W-5) and hexagonal MoO3 (Mo1) samples. Upper and lower pictures are before and after UV irradiation for 10 min, respectively.
Figure 6. Diffuse reflectance spectral changes of all WO3 samples (W-1 ~ W-5) after UV irradiation.
Figure 7. The time dependence of relative reflectance changes at 660 nm for all WO3 samples after UV irradiation.
CE = ΔOD/It = log(Ro/R)/It (1)
where CE is the photochromic colorization efficiency, ΔOD is optical density difference, I is the light intensity, t is irradiation time, Ro and R are the reflectance before and after irradiation. The initial photochromic colorization efficiencies were obtained as 48, 7.2, 5.5, 2.0 and 5.7 cm2/kJ for W-1, W-2, W- 3, W-4 and W-5, respectively. The PC response times of all WO3 samples were roughly estimated by fitting the decay of relative reflectance to single exponential decay function and the estimated time constants (τ) and CE were together listed in Table 1 with crystallite size and bandgap energy.
The large difference in photochromic effect between W-1 and other samples cannot be explained by the difference in crystallite size alone. One difference between the W-1 sample and the other samples except W-5 is whether it was thermally treated in air or not. Thermal treatment in air is thought to not only reduce the surface area by increasing the crystallite size, but also play a role in reducing defect sites caused by the lack of oxygen atoms on the crystal surface, so the difference in photochromic effect seems to be related to this matter. Among the proposed mechanisms of photochromism of WO3, briefly mentioned in the introduction, the double charge injection model has been frequently cited. In this model, surfaceadsorbed H2O is oxidized by photogenerated holes in WO3 and electron and proton simultaneously are injected into WO3 to form hydrogen tungsten bronze (HxWO3). The blue color of the bronze is explained by the absorption of visible and near infrared light to cause intervalence charge transfer between W5+ and W6+ in the bronze. This model necessarily requires an oxidizable proton source like adsorbed H2O. In order to check the role of the water adsorbed at ambient condition, W-1 sample was dried at 130 ℃, overnight and the photochromic effect was immediately rechecked. The drying at 130 ℃ overnight didn’t change any crystallinity of the sample as confirmed from XRD patterns (not shown here). Fig. 8 shows the photochromic spectral changes before and after drying. The photochromic effect after drying slightly decreased. But this is not enough to say the adsorbed H2O is critical to the photochromism. Therefore, it is still unclear whether H+ intercalates into the WO3 lattice to form the hydrogen bronze claimed by the double charge injection mechanism. However, our results clearly indicate that heat treatment, which changes crystallinity, reduces the photochromic effect, regardless of size or structure, as shown by the large difference between W-1 and W-2 to W-4.
Figure 8. Photochromic reflectance spectra changes of W-1 before and after drying.
MoO3, like WO3, is well known for its photochromic properties and has been described with the same photochromic mechanisms.3 Fig. 9 shows photochromic spectral changes of two MoO3 samples with different crystal structures shown in XRD patterns. While all WO3 samples show the similar photochromic spectral change, the spectral changes of two MoO3 samples are quite different. Mo-1 (h-MoO3) showed a maximum absorption at 660 nm and rapidly decreased absorption at wavelengths above 700 nm, whereas Mo-2 (α-MoO3) showed a similar spectral change to those of the WO3 samples, showing stronger absorption in the longer wavelength region. Fig. 10 shows the time dependence of relative reflectance change at 660 nm of two MoO3 samples. The photochromic changes of both MoO3 are slower than W-1 to W-4 samples but closer W-5. Since Mo-1, Mo-2 and W-5 have a similar crystallite size that are much larger than W-1 to W-4, the photochromic response time must be definitely related to the crystallite size. Mo-2 showed a lower response time than Mo-1. Since Mo-2 was thermally converted from Mo-1 and had a larger crystallite size than Mo-1 as shown in Table 1, this can be understood as the relationship between response time and crystallite size as mentioned before. The initial colorization efficiencies of Mo-1 and Mo-2 at 660 nm were estimated as 8.3 and 3.9 cm2/kJ, respectively. Despite that bandgap energy of Mo-1 is larger than that of Mo-2 and Mo-2 can absorb more light, Mo-1 showed more than twice better photochromic effect than Mo-2. Xie et al. reported previously the difference of photochromic effect, which is similar to ours, between aqueous α-MoO3 and h-MoO3 nanobelt suspensions prepared by hydrothermal method.22(a) They also explained the photochromic effect of MoO3 as a double charge injection mechanism like WO3. The holes formed by photoexcitation oxidized water and the produced hydrogen ions entered the crystal to form hydrogen molybdenum bronze. The difference in photochromic effect between α-MoO3 and h-MoO3 was explained by the difference in crystalline structure: the hexagonal channel structure of h-MoO3 provides a better pathway for hydrogen ion than zigzag channel of α-MoO3. This explanation may be plausible for their result where the experiment was conducted in an aqueous suspension, However, since all our experiments were performed with powder samples in ambient conditions and showed that adsorbed water is not critical for the photochromic effect, there are considerable questions about the double charge injection mechanism. We therefore believe that the photochromic difference between α-MoO3 that was thermally converted from h-MoO3 and h-MoO3 originates from heat treatment like the WO3 samples. All our results consistently support that heat treatment causing crystallinity change reduces the photochromic effect. As mentioned earlier, the thermal treatment at more than 300 ℃ in air usually decreases the surface oxygen defect sites as well as changing crystal structure. Thus, the photochromic effect strongly related to surface oxygen defect sites to play a role of deep trap sites for photoexcited electrons and holes that are oxygen radicals.
Figure 9. Diffuse reflectance spectral changes of Mo-1 and Mo-2 after UV irradiation.
Figure 10. The time dependence of relative reflectance changes at 550 nm for Mo-1 and Mo-2 samples.
Reversibility in photochromism is very important not only in terms of applications but also in understanding the mechanism of photochromism. In this study we only focused on the coloring process of WO3 and MoO3. Their bleaching process in ambient atmosphere is so slow that comparative studies are not easy. In order to understand the photochromic mechanism of these materials more clearly, the bleaching process under various conditions will be conducted in future studies.
CONCLUSION
In summary, we have investigated comparatively the photochromic (PC) effect of nanocrystalline WO3 and MoO3 with various crystallite size and crystal structures prepared by simple acid precipitation and subsequent thermal treatment. The photochromic effect of h-WO3 and γ-WO3 is highly dependent on crystallite size rather than crystalline structure. The smaller the crystallite size, the better the PC effect.
However, the orthorhombic WO3·H2O and MoO with hexagonal an orthorhombic crystal structures did not follow this trend. One consistent result for all WO3 and MoO3 samples is that the thermal treatment in air, which changes crystallinity, whether it changes the crystal structure or only the crystallite size, reduces the PC effect. Since the thermal treatment reduces the surface oxygen defect sites, we believe that PC effect of WO3 and MoO3 depends critically on the surface oxygen defect sites and these sites serve as deep trap sites for photogenerated electrons and oxygen radical holes. In this case the proton insertion for charge compensation claimed by double charge injection mechanism may not be necessary.
Acknowledgments
This work was supported by Pukyong National University Research Fund (Year 2021-2022).
References
- (a) "Photochromism: Molecules and Systems" Ed. Durr, H.; Bouas-Laurent, H. 2nd ed. 2003, Elsevier.
- (b) Ashrit, P. Chapter 5 in Transition Metal Oxide Thin Film-Based Chromogenics and Devices, 2017, Elsevier.
- Wang, S.; Fan, W.; Liu, Z.; Yu, A.; Jiang, X. J. Mater. Chem. C, 2018, 6, 191.
- (a) Yao, J. N.; Hashimoto, K.; Fujishima, A. Nature, 1992, 355, 624.
- (b) He, T.; Yao, J. J. Photochem. Photobiol. C, 2003, 4, 125. https://doi.org/10.1016/S1389-5567(03)00025-X
- (a) Yao, J. N.; Loo, B. H.; Hashimoto, K.; Fujishima, A. Ber. Bunsen-Ges. 1992, 96, 699.
- (b) Pan, L.; Wang, Y.; Wang, X.; Qu, H.; Zhao, J.; Li, Y.; Gavrilyuk, A. Phys. Chem. Chem. Phys. 2014, 16, 20828. https://doi.org/10.1039/c4cp02834a
- Gavrilyuk, A.; Tritthart, U.; Gey, W. Phys. Chem. Chem. Phys. 2011, 13, 9490.
- Highfield, J. G.; Gratzel, M. J. Phys. Chem. 1988, 92, 464.
- (a) Deb, S. K. Philos. Mag. 1973, 27, 801.
- (b) Deb, S. K.; Chopoorian, D. M. J. Appl. Phys. 1966, 37, 4818. https://doi.org/10.1063/1.1708145
- (c) Deb, S. K. Proc. R. Soc. A, 1968, 304, 211.
- (a) Duncan, R. C.; Faughnan, B. W.; Phillips, W. Appl. Opt. 1970, 9, 2236.
- (b) Crandall, R. S.; Faughnan, B. W; Appl. Phys. Lett. 1975, 26, 120. https://doi.org/10.1063/1.88087
- (a) Schirmer, O. F.; Wittwer, V.; Baur, G.; Brandt, G. J. Electrochem. Soc. 1977, 124, 749.
- (b) Gerard, P.; Deneuville, A.; Hollinger, G. J. App. Phys. 1978, 49, 409. https://doi.org/10.1063/1.324349
- Supothina, S. S.; Seeharaj, P.; Yoriya, S.; Sriyudthsak, M. Ceramics. Int. 2007, 33, 931.
- Wriedt, H. Bull. Alloy Phase Diag. 1989, 10, 368.
- (a) Calber, G.-; Bardeleben, H. J. V.; Fleury, B.; Buissette, V.; Sanchez, C.; Mercier, T. L.; Cassaignon, S.; Portehault, D. Nature Comm. 2019, 327.
- (b) Sun, W.; Yeung, M. T.; Lech, A. T.; Lin, C.-W.; Lee, C.; Li, T.; Dung, X.; Zhou, J.; Kaner, R. Nano Lett. 2015, 15, 4834. https://doi.org/10.1021/acs.nanolett.5b02013
- Li, N.; Zhao, Y.; wang, Yi; Lu, Y.; Song, Y.; Hung, Z.; Li, Y.; Zhano, J. Eur. J. Inorg. Chem. 2015, 2804.
- Gerand, B.; Nowogrocki, G.; Guenot, J.; Figlarz, M.; Chimie, D.; De Picardie, U.; De Leu, S.; Chimie, L.; De Paris, U.; De Cedex, J. P.; Al, E. T. J. Solid State Chem. 1979, 29, 429. https://doi.org/10.1016/0022-4596(79)90199-3
- Chithambararaj, A; Yogamalar, R. N.; Bose, A. C. Cryst. Growth Des. 2016, 16, 1984.
- Kopp, L.; Harman, B. N.; Liu, S. H. Solid State Commun. 1977, 22, 677. https://doi.org/10.1016/0038-1098(77)90248-4
- Tauc, J. Mater. Res. Bull. 1968, 3, 37.
- Huang, K. H.; Zhang, Q.; Yang, F.; He, D. Nano Res. 2010, 3, 281. https://doi.org/10.1007/s12274-010-1031-3
- Kalanur, S.; Seo, H. J. Coll. Interf. Sci. 2018, 509, 440.
- Suzuki, H.; Tomita, O.; Higashi, M.; Abe, R. J. Mater. Chem. A 2017, 5, 10280.
- Wongkrua, P.; Thongtem, T.; Thongtem, S. J. Nanomater. 2013, 2013, 79.
- (a) Zheng, L.; Xu, Y.; Jin, D.; Xie, Y. Chem. Mater. 2009, 21, 5681.
- (b) Yun, J.; Jang, W.; Lee, T.; Lee, Y.; Soon, A. Phys Rev. Lett. 2017, 7, 024025.
- (a) Songara, S.; Gupta, V.; Patra, M. K.; Singh, J.; Saini, L.; Gows, G. S.; Vadera, S. R.; Kumar, N. J. Phys. Chem. Solids 2012, 73, 851.
- (b) Li, N.; Z hao, Y. Wang, Y.; Lu, Y.; Song, Y.; Huang, Z.; Li, Y.; Zhao, J. Eur. J. Inorg. Chem. 2015, 2804.
- He, Y.; Zhao, Y. J. Phys. Chem. C. 2008, 112, 61.