Formation and evolution of orientation-specific CO2 chains on nonpolar ZnO(10͞10) surfaces.

Clarifying the fundamental adsorption and diffusion process of CO2 on single crystal ZnO surfaces is critical in understanding CO2 activation and transformation over ZnO-based catalysts. By using ultrahigh vacuum-Fourier transform infrared spectroscopy (UHV-FTIRS), we observed the fine structures of CO2 vibrational bands on ZnO(100) surfaces, which are the combinations of different vibrational frequencies, originated from CO2 monomer, dimer, trimer and longer polymer chains along [0001] direction according to the density functional theory calculations. Such novel chain adsorption mode results from the relatively large attractive interaction between CO2 and Zn3c atoms in [0001] direction. Further experiments indicate that the short chains at low coverage evolve into long chains through Ostwald ripening by annealing. At higher CO2 coverage (0.7 ML), the as-grown local (2 × 1) phase of chains first evolve into an unstable local (1 × 1) phase below 150 K, and then into a stable well-defined (2 × 1) phase above 150 K.

CO2 activation and transformation are the key steps in CO2 utilizations such as in environment protections and renewable energy fields. One famous example about CO2 utilization is the methanol synthesis over Cu/ZnO/Al2O3 catalysts by syngas (CO/CO2/H2) industrially1234. In the three-way catalyst, ZnO plays a pivotal role in CO2 activation and stabilization2, however, to elucidate the underlying mechanism of CO2 activation at the molecular level, studies on well-defined single crystal ZnO surfaces under ultrahigh vacuum (UHV) conditions are essential5. The nonpolar mixed-terminated ZnO(100) surface is the energetically most favorable surface and dominates the exposed surfaces of ZnO particles in applications37. Therefore, the investigation of CO2 adsorption and activation behavior on ZnO(100) surfaces is typical and helpful in understanding of the ZnO-catalyzed CO2 chemistry.

The mixed-terminated ZnO(100) surface is composed of rows of ZnO “dimers” separated by trenches along [110] direction. The ZnO “dimer” consists of one threefold coordinated surface Zn cation (Zn3c) and the adjacent surface O anion (O3c), running along the crystallographic [0001] direction89. An early near edge X-ray absorption fine structure (NEXAFS) study suggested a bidentate adsorption configuration of CO2 on ZnO(100) surfaces10. Later, the high resolution electron energy loss spectroscopy (HREELS) results together with density functional theory (DFT) calculations5 supported an unusual tridentate carbonate configuration: the middle C-atom bound to the surface O3c anion and the two end O-atoms of CO2 molecule bound to two surface Zn3c cations along [0001] direction. Besides, experiments also observed two ordered carbonate adlayers: the close packed (1 × 1) phase corresponds to 1 ML CO2 coverage, and the open (2 × 1) phase to 0.5 ML5. (Here, 1 ML is defined as the density of the surface Zn3c cations on the clean surface.) In the (2 × 1) phase, the free surface Zn3c sites were found surprisingly to bind CO more strongly due to charge transfer, which suggests the potential importance of such open phase in the polybasic catalytic reactions such as methanol synthesis from syngas1112.

Up to date, only the tridentate carbonate configuration, and the ordered (2 × 1) and (1 × 1) phases are reported for CO2 on ZnO(100)567891011. We wonder how the initial CO2 tridentate carbonates evolve into the ordered adlayers? What is the mechanism behind the evolution from tridentate carbonates to the ordered adlayers? Monitoring the molecular vibration of adsorbed CO2 is an important way to clarify such evolution process, since the molecular vibration is very sensitive to the chemical surrounding change of the adsorbed molecules. For instance, by using infrared spectroscopy, Heidberg et al. studied the dynamic dipole-dipole-coupling of adjacent CO2 molecules on MgO(100) and NaCl(100) surfaces13141516. Similar results were reported for CO2 adsorption on rutile TiO2(110) surfaces in our previous work17. On the ZnO(100) surface, using the high resolution ultrahigh vacuum-Fourier transform infrared spectroscopy (UHV-FTIRS), Buchholz et al. characterized the ordered (2 × 1) phase at high CO2 coverage6, but they did not observe the evolution from the initial tridentate carbonates to the ordered adlayers.

In this paper, based on the high resolution UHV-FTIRS and DFT calculations, we reported the fine structures combined by of CO2 vibrational levels on ZnO(100) surfaces with increasing CO2 coverage, which are attributed to the formation of [0001]-oriented short CO2 polymer chains consisting of monomer, dimer, trimer and so on. The as-grown chains at high CO2 coverage further evolve into the unstable local (1 × 1) phase below 150 K, and then relax into the stable (2 × 1) phase above 150 K.

Results and Discussion

Figure 1 presents the polarized infrared reflection absorption spectroscopy (IRRAS) data of adsorbed CO2 on ZnO(100) surfaces at 90 K with IR light incident along [110] direction. At the initial adsorption with 0.1 L (1 L = 1.33 × 10−6 mbar·s) CO2 dosage, one vibrational band at 1622 cm−1 first appears in s-polarized spectra (Fig. 1a); in p-polarized spectra, two bands at 1297 and 978 cm−1 (Fig. 1b) emerge simultaneously. These dramatically lowered vibration frequencies of adsorbed CO2 relative to that of gas-phase CO2 (2349 cm−1) demonstrate evidently that CO2 have chemically adsorbed on ZnO(100) surfaces. Based on the IR judgement principle on dielectric substrates1718, the vibrational bands at 1622, 1297 and 978 cm−1 for CO2 on ZnO(100) surfaces are assigned respectively to the asymmetrical stretching mode (νas(OCO), in-plane), symmetrical stretching mode (νs(OCO), out-of-plane) and the stretching vibration between the carbon atom and the underneath surface O3c atom (ν(CO3c), out-of-plane) for tridentate carbonates. (See Figure S1 for the detailed judgements of the CO2 vibration direction through the polarized IRRAS.) Our assignment is in accordance with previous HREELS and FTIR reports5678910111213141516171819.

Figure 1

IRRA spectra of CO2 adsorbed on ZnO(100) surfaces as a function of CO2 dosage by using (a) s-polarized and (b) p-polarized IR beams, respectively. The IR light incident is along [110] direction. All spectra were acquired at 90 K.

It is interesting that we observed the fine structures of adsorbed CO2 vibrational levels with increasing CO2 coverage. For νas vibration, besides the 1622 cm−1 band (Fig. 1a), a new band appears at 1582 cm−1, and they finally converge into one intense band at 1590 cm−1 for saturated CO2 coverage (2 L). Simultaneously, for νs vibration (Fig. 1b), besides the 1297 cm−1 band, 1313 and 1337 cm−1 bands appear and finally evolve to one sharp 1340 cm−1 band. The ν(CO3c) band evolves from 978 cm−1 to 1008 cm−1 gradually. Therefore, we can divide the fine structures of tridentate carbonate vibration into four different groups as shown in Fig 1: I 1622, 1297, 978 cm−1; II 1582, 1313, 978 cm−1; III 1582, 1337, 1008 cm−1; IV 1590, 1345, 1008 cm−1, which correspond to four configurations of tridentate carbonates. It is clear that the frequencies of the same vibrational modes in different CO2 configurations are well distinguished except a few vibrational modes at 978, 1008 and 1582 cm−1.

To examine whether the fine structures of CO2 vibration is associated with the surface defects, we treated the ZnO(100) surface in atomic oxygen atmosphere of 2 × 10−6 mbar at 750 K and 10 L O2 at 90 K, but no change in the fine structures is observed. Therefore, the prepared ZnO(100) surface can be regarded as the stoichiometric surface with negligible surface defects. Recent STM studies also reported that no apparent oxygen vacancies or miss zinc-oxygen dimers were observed on the vacuum-annealed ZnO(100) surface20. Actually, the fine structures of molecular vibrations have already been observed on some single crystal substrates, such as CO and CO2 on MgO(100) and NaCl(100), which were attributed to the dipole-dipole-coupling in ordered molecular layers13141516 instead of the influence of the surface defects.

To clarify the interaction between CO2 on ZnO(100), different configurations of two carbonates were designed and checked by DFT calculations. We first calculated the structure of CO2 monomer on the ZnO(100) surface of a (2 × 4) supercell with the long axis along [0001] direction. As shown in Fig. 2a, the tridentate configuration was confirmed, which is in well agreement with previous calculation results5. Then three distinct configurations of two CO2 molecules were calculated. The results reveal that when two CO2 form a chain along [0001] direction, the binding energy per molecule is the lowest, corresponding to the most stable configuration compared to the other distributions of the two CO2. The similar results were also calculated in a recent literature21.

Figure 2

Different distributions for CO2 molecules on ZnO(100) surfaces and their DFT-calculated binding energies.

(a) single CO2 molecule, (b) two CO2 in diagonal, (c) two adjacent CO2 along [110] direction and (d) two adjacent CO2 along [0001] direction. The binding energy for each configuration is given in eV per CO2 molecule, where the negative energy means the adsorption is exothermic. Charge density difference maps for (e) single CO2 molecule and (f) CO2 chains along [0001] direction. Electron accumulation and depletion regions are shown in blue and red, respectively. Panels (a–d) represent the top view and panels (e,f) for the side view.

Based on the above calculation result, we designed CO2 molecular chains with different length on the ZnO(100) surfaces of a (6 × 2) supercell and calculated the corresponding vibrational frequencies. The chain contains, respectively, one, two, three, four, five and infinite carbonates arraying end to end along the long axis [0001] direction. All the results are shown in Table 1. We can see that the calculated 1585 cm−1, 1261 cm−1 and 958 cm−1 for the monomer respectively correspond to the experimental results of Group I: 1622 cm−1, 1297 cm−1 and 978 cm−1. Thus the bands of Group I are assigned to the carbonate monomer vibrations. It is easy to understand that most CO2 are diluted at the initial adsorption on the surface to form carbonate monomers.

Table 1

Calculated vibrational frequencies of CO2 in different configurations and distributions on ZnO(10 0) surfaces.

		Theory			Experimenta
		νas(cm−1)	νs(cm−1)	ν(cm−1)	νas(cm−1)	νs(cm−1)	ν(cm−1)
CO2 chain	Monomer	1585	1261	958	1622	1297	978
	Dimer	1546	1280	984/951	1582	1313	1008/978
	Trimer	1542	1290	984/950	1582	1320	1008/978
	Tetramer	1541	1298	987/981/976/955	1582	1337	1008/978
	Pentamer	1540	1302	986/982/978/976/949
	infinite chain	1563	1310	980	1590	1345	1008
chains distribution	isolated	1567	1303	976
	spacing	1574	1303	968	1590	1345	1008
	neighbouring	1593	1309	966	1618	1345	1001

The corresponding experimental values are also listed for comparison.

As shown in Table 1, the νas dramatically redshifts from 1585 cm−1 of the monomer to 1546 cm−1 of the dimer, and it slightly changes from 1546 to 1540 cm−1 with increasing CO2 from dimer to pentamer. Further lengthening the chain to infinite, the νas blueshifts back to 1563 cm−1. On the contrary, the νs of monomer to pentamer increases monotonously from 1261 to 1302 cm−1, and further to 1310 cm−1 for infinite length. The evolution trend of the calculated results is well consistent with our experimental IR frequencies. Therefore, the fine structures of CO2 vibrations originate from the short CO2 chains composed of monomer, dimer, trimer and so on when dosing CO2 from 0.1 L to 1 L. Accordingly, in Table 1 the strong bands of Group IV measured at saturated CO2 dosage (2 L) are assigned to the infinitely long chain vibrations.

Generally, the one-dimensional chain formation requires the symmetry loss of substrate surfaces22, such as trenches2324 or steps2526 on surfaces along specific direction, i.e., the space restriction plays a major role in the chain formation. However, on ZnO(100) surfaces, the CO2 chains are along [0001] direction, rather than along the surface trench direction [110]. To explore the formation mechanism of such CO2 chains, we performed the charge transfer analysis of CO2 monomer and CO2 chains on ZnO(100) surfaces by DFT calculations. The calculated charge density difference maps are shown in Fig. 2e and f, respectively. The bonding formation obviously bents the linear CO2 and induces the charge redistribution. In the single carbonate, as shown in Fig. 2e, the two O atoms of CO2 get more electrons while the C and Zn3c atoms lose more electrons. The charge redistribution induces the extra Coulomb attraction between Zn3c atoms and O atoms of CO2. For the carbonate chain, as shown in Fig. 2f, two O atoms bond to one Zn3c atom, and the induced positive electricity of Zn3c atoms is evidently enhanced. As a result, the extra attractive Coulomb interaction between Zn3c atoms and O atoms of CO2 is strongly enhanced. Such enhanced attractive interaction makes the chain configuration of adsorbed CO2 along [0001] direction more stable. Along [110] direction, on the contrary, the enhanced electrostatic repulsion between CO2 molecules causes the CO2 alignment along [110] direction less stable.

To understand the phase evolution of CO2 adlayers on ZnO(100), the temperature dependence of CO2 chains was studied for fixed CO2 coverages. For the low CO2 coverage of 0.2 ML (corresponding to 0.2 L), the IRRAS results are shown in Fig. 3. Slowly annealing to 230 K, the νas band at 1623 cm−1 gradually converts to 1583 cm−1 in Fig. 3a. At the same time, the three close peaks (1297, 1313, 1325 cm−1) of νs finally convert to one peak at 1337 cm−1, and the 978 cm−1 band to 1008 cm−1 band, as shown in Fig. 3b. (The corresponding p-polarized spectra with IR light incident along [0001] direction can be seen in Figure S2 in the SI.) Such band conversions reveal that the chain conversions from the monomer to long chains happened upon annealing through Ostwald ripening. Our present study provides an effective way to synthetize long CO2 chains along [0001] direction on ZnO(100) surfaces.

Figure 3

IRRA spectra of 0.2 ML CO2 adsorbed on ZnO(100) surfaces with annealing by using (a) s-polarized and (b) p-polarized IR beams, respectively. The IR light incident is along [110] direction. All spectra were acquired at 90 K.

For the high coverage of 0.7 ML (corresponding to 2 L), the IRRAS results with annealing are shown in Fig. 4a and b. As mentioned before we have assigned the 1590 cm−1 to the νas and 1345 cm−1 to the νs of long CO2 chains at 90 K. Slowly annealing to 150 K, the νas band at 1590 cm−1 unexpectedly converts to 1618 cm−1 gradually, as shown in Fig. 4a. Further annealing to 240 K, the band gradually redshifts back to 1590 cm−1. On the other hand, the νs band (1345 cm−1) keeps constant from 90 K to 150 K. Over 150 K, its intensity slightly decreases with a weak redshift. (The corresponding p-polarized spectra with IR light incident along [0001] direction can be seen in Figure S3 in the SI). It is easy to know that for the high CO2 coverage and relative high annealing temperature, the length change of a single long chain will not induce such obvious changes of the CO2 vibration frequencies. But the change of the separated distance along [110] direction between two long chains may induce significant changes of the CO2 vibration frequencies in Fig. 4 due to the interchain interaction.

Figure 4

IRRA spectra of saturated adsorbed CO2 (0.7 ML) on ZnO(100) surfaces with annealing by using (a) s-polarized and (b) p-polarized IR beams, respectively. The IR light incident is along [110] direction. All spectra were acquired at 90 K. (c) Schematic kinetic model to illustrate the mechanism of the phase evolution during annealing. The zigzag lines denote the different lengths of CO2 chains.

Thereafter, we calculated a series of CO2 long chains with different interchain distances, such as the isolated long chain, two neighbouring chains (corresponding to (1 × 1) phase) with the shortest distance of a0, and two spacing chains (corresponding to (2 × 1) phase) with 2a0. Here a0 represents the lattice constant along [110] direction of ZnO(100) surfaces. The calculated results are shown in Table 1. We found that the calculated νas of the spacing chains is 1574 cm−1, which is consistent with the experimental 1590 cm−1 at 90 K and 240 K. The calculated νs of the spacing chains is 1303 cm−1, which is in agreement with the vibration frequencies of 1345 cm−1 at 90 K and 240 K. These obviously indicate that the experimentally observed vibration bands at both 90 K and 240 K belong to the spacing chains. Similarly, the experimentally observed vibration frequencies at 150 K belong to the neighbouring chains.

In Fig. 4c, we give a schematic evolution picture of the CO2 chains with increasing temperature. At low temperature of 90 K, on the one hand, the CO2 chains with various length are randomly distributed on the surface due to the low kinetic energy; on the other hand, most of the interchain spaces along [110] are equal to 2a0 caused by the interchain repulsion, forming the local (2 × 1) phase, as shown in Fig. 4c-I. Annealing to 150 K, the CO2 diffusion is enhanced to induce the Ostwald ripening between CO2 chains: the CO2 molecules detach from the short chains and attach to the long ones. Finally, the lengthened chains become neighbouring with others, forming the local (1 × 1) phase, as shown in Fig. 4c-II. However, due to the strong repulsive interaction between neighbouring chains, the (1 × 1) phase is an unstable intermediate state. Further annealing the intermediate state to 240 K, all the chains will relax to the more stable spacing structure, forming the well-defined stable (2 × 1) phase, as shown in Fig. 4c-III.

Conclusion

In conclusion, the formation and evolution of CO2 chains on ZnO(100) surfaces was studied by employing UHV-FTIRS and DFT calculations. We observed the fine structures of CO2 vibrational levels on ZnO(100) surfaces, which are attributed to the formation of CO2 monomer, dimer, trimer and longer polymer chains along [0001] direction. At high CO2 coverage, the as-grown local (2 × 1) phase composed of chains with various lengths further evolve into the unstable local (1 × 1) phase with annealing to 150 K, and then relax into the well-defined stable (2 × 1) phase above 150 K. The mechanisms of the chain formation and evolution were discussed by DFT calculations and a schematic kinetic model.

Methods

Experimental details

The experiments were carried out using an ultrahigh vacuum (UHV) system17 (the base pressure better than 6 × 10−11 mbar) equipped with a vacuum Fourier transform infrared spectroscopy (FTIR) spectrometer (Bruker, VERTEX 80 v), a low energy electron diffraction (LEED)/Auger (AES) spectrometer with gain power of microchannel plates BDL 600IR-MCP. The clean mix-terminated ZnO(100) (8 × 8 × 1 mm, MTI) surface was prepared by repeated cycles of Ar+ sputtering and annealing at 800 K under UHV conditions until no impurities were detected by AES and clear (1 × 1) LEED patterns were obtained. Then, the clean ZnO(100) was oxidized in oxygen atmosphere (5 × 10−7 mbar) at 750 K for 20 minutes. The IR measurements were performed using infrared reflection absorption spectroscopy (IRRAS) mode with a fixed incidence angle of 80°. The recorded data, i.e., the absorbance is defined as A = log10(R0/R), where R0 and R are the reflected signals from the bare and the adsorbate covered surfaces, respectively. The optical path was evacuated in order to avoid any unwanted IR adsorption from gas phase species. High purity CO2 (99.99%) and O2 (99.999%) were dosed via backfilling in the experiments.

Computational details

First-principles calculations were performed using the Vienna ab-initio simulation package (VASP)2728 with a cut-off energy of 500 eV for the basis set. Γ-point was used for Brillouin zone sampling. The projector-augmented wave method (PAW)29 with the PBE type exchange–correlation potentials30 was adopted. To model the ZnO(100) surface, the optimized lattice parameters of bulk ZnO, a = 3.285 Å and c/a = 1.6131, were used to build slabs with six ZnO layers. Two surface unit cells, which have dimensions of 6 × 2 and 2 × 4 along[0001] and [110] directions, respectively, were employed to perform the calculations. The atomic positions of top three layers were optimized until the forces are less than 0.03 eV/Å, while the bottom layers were fixed at bulk positions. A vacuum layer with a thickness of 15 Å was used to minimize interactions between adjacent slabs. The vibrational frequencies were derived from Hessian matrix calculated by finite-displacement method.

Additional Information

How to cite this article: Cao, Y. et al. Formation and evolution of orientation-specific CO2 chains on nonpolar ZnO(1010) surfaces. Sci. Rep. 7, 43442; doi: 10.1038/srep43442 (2017).

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.