Open Access Article
Tahereh Mohammadi
Hafshejani†
a,
Weijia
Wang†
a,
Jonas
Heggemann
b,
Alexei
Nefedov
a,
Stefan
Heissler
a,
Yuemin
Wang
a,
Philipp
Rahe
b,
Peter
Thissen
a and
Christof
Wöll
*a
aKarlsruher Institut für Technologie (KIT), Institut für Funktionelle Grenzflächen (IFG), Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany. E-mail: christof.woell@kit.edu
bFachbereich Physik, Universität Osnabrück, Barbarastrasse 7, 49076 Osnabrück, Germany
First published on 30th June 2020
Detailed information on structural, chemical, and physical properties of natural cleaved (10.4) calcite surfaces was obtained by a combined atomic force microscopy (AFM) and infrared (IR) study using CO as a probe molecule under ultrahigh vacuum (UHV) conditions. The structural quality of the surfaces was determined using non-contact AFM (NC-AFM), which also allowed assigning the adsorption site of CO molecules. Vibrational frequencies of adsorbed CO species were determined by polarization-resolved infrared reflection absorption spectroscopy (IRRAS). At low exposures, adsorption of CO on the freshly cleaved (10.4) calcite surface at a temperature of 62 K led to the occurrence of a single C–O vibrational band located at 2175.8 cm−1, blue-shifted with respect to the gas phase value. For larger exposures, a slight, coverage-induced redshift was observed, leading to a frequency of 2173.4 cm−1 for a full monolayer. The width of the vibrational bands is extremely small, providing strong evidence that the cleaved calcite surface is well-defined with only one CO adsorption site. A quantitative analysis of the IRRA spectra recorded at different surface temperatures revealed a CO binding energy of −0.31 eV. NC-AFM data acquired at 5 K for sub-monolayer CO coverage reveal single molecules imaged as depressions at the position of the protruding surface features, in agreement with the IRRAS results. Since there are no previous experimental data of this type, the interpretation of the results was aided by employing density functional theory calculations to determine adsorption geometries, binding energies, and vibrational frequencies of carbon monoxide on the (10.4) calcite surface. It was found that the preferred geometry of CO on this surface is adsorption on top of calcium in a slightly tilted orientation. With increased coverage, the binding energy shows a small decrease, revealing the presence of repulsive adsorbate–adsorbate interactions.
Calcite is an ionic crystal consisting of calcium cations (Ca2+) and carbonate (CO32−) anions, which are held together through ionic bonding.7 The structure is trigonal (rhombohedral holohedral
2/m) with a R
c space group.7c Cleavage is easiest for the thermodynamically most stable calcite(10.4) surface with a surface energy calculated around 0.59 J m−2.8 The calcite(10.4) surface is furthermore nonpolar, namely a type I surface according to Tasker's classification, and thus bears zero charge (q = 0) and zero dipole moment (μ = 0).9 Accordingly, most previous surface studies reported for this material focused on this surface. The calcite(10.4) surface (a structural model is reproduced in Fig. 1) contains a total of two calcium atoms and two carbonate groups within the (1 × 1) unit cell (with a unit cell size of 5.0 × 8.1 Å2), both are present in the top layer. Note that the two carbonate groups are rotated with respect to each other and, as a consequence, the upmost oxygen atom points either to the left or to the right along the [42
] direction.7b,10
The microscopic nature of the calcite(10.4) surface has intensively been studied by various experimental10,11 and theoretical methods.12 Two deviations from the bulk-truncated (1 × 1) structure have been reported with their origins and geometry still being under discussion. First, a loss of the pg surface symmetry (also referred to as “row-pairing”) has been observed within atomic force microscopy (AFM) data,10,11,13 yet, this surface relaxation can also be explained by an asymmetric tip probing an unrelaxed surface. Second, the existence of a (2 × 1) reconstruction has been proposed on the basis of low energy electron diffraction (LEED)11d and AFM data.11b,14 This reconstruction is characterised by a rectangular supercell with dimensions of 10.0 Å and 8.1 Å oriented along the [0
0] and [42
] directions, respectively.
Here, the surface structure of the UHV-cleaved calcite(10.4) surface is investigated by polarization-resolved infrared reflection absorption spectroscopy (IRRAS) using the surface-ligand IR (SLIR) spectroscopy approach,15 with CO as the probe molecule, and complemented by non-contact atomic force microscopy (NC-AFM) experiments and density functional theory (DFT) calculations. We would like to point out that there are no previous IRRAS studies of adsorbates on well-defined surfaces of macroscopic calcite single crystals available. This lack of experimental data is due to technical problems in recording IR-data using reflection geometries for dielectric substrates, which could only recently be overcome.15,16 The IRRAS results of CO on calcite(10.4) surfaces obtained at low temperature show only one sharp absorption band, providing strong evidence that this surface exhibits only one type of Ca2+ binding site. To our knowledge, these are the first IRRAS data for any adsorbates bound to a calcite single crystal substrate.
The calculations reported here were carried out for a slab (thickness of 4 layers) exposing a (10.4) surface with a (1 × 1) unit cell and containing a total number of 40 atoms (Ca: 8, C: 8, O: 24), see also Fig. 1. This model is similar to that used in earlier studies.8b,24 A vacuum region of 20 Å was found to be large enough to avoid interactions between the surface and adsorbed molecules and their replica images. We define coverage as the ratio of the number of CO molecules to the number of surface Ca ions. All degrees of freedom, except the bottom layer, were allowed to relax until the forces on the atoms were below 1 meV Å−1. The optimization of the atomic coordinates (and unit cell size for the bulk materials) was performed through a conjugate gradient technique, which utilizes the total energy and the Hellmann–Feynman forces on the atoms (and stresses on the unit cell). The structures were considered to be fully relaxed when the forces on the atoms were smaller than 1 meV Å−1. The Brillouin zone integration was performed using a 4 × 4 × 1 mesh within the Monkhorst–Pack scheme.25 Vibrational eigenmodes and frequencies were calculated by the force-constant approach. The start geometry for the optimization was a CO molecule (bond length 1.14 Å) positioned on top of Ca2+ in an upright orientation.
The surface energy is defined as the difference in total energy between the bulk and the surface per unit area and, for a stoichiometric surface, can be obtained by the following equation:26
| Esurf = (Eslab − NEbulk)/2A |
The adsorption energy (Eads) can be calculated as a function of the coverage of the CO molecules:28
The CO vibration was also observed in IRRAS data recorded with s-polarized light (Fig. 3b, bottom). However, the intensity is much smaller, by an order of magnitude, compared to that excited by p-polarized light (Fig. 3b, top). The occurrence of a vibrational feature also for s-polarization reveals that there must be a static or dynamic tilt angle, because only vibrational modes with a component of their transition dipole moment (TDM) orientated parallel to the surface can be excited by s-polarized light.16
In order to determine the binding energy of the CO adlayers we recorded IRRAS-data after heating to different temperatures, see Fig. 4a. In Fig. 4b we present the variation of the integrated band intensities as a function of surface temperature. The position (Fig. 4a) and the intensity (Fig. 4b) of the CO–Ca2+ band stays constant for temperatures below 92 K. Upon further annealing to 100 K, CO starts to desorb from the surface, yielding a decrease in band intensity accompanied by a slight blue shift of the vibrational frequency. A simple quantitative analysis of the IR band intensity evolution, using the Redhead equation,29 yields a binding energy of −0.31 eV for CO on the calcite(10.4) surface.
Next, we recorded AFM micrographs to investigate the adsorption position and CO arrangement after dosing CO into the cold scan head on calcite(10.4) for three different total times of 25 s, 50 s, and 150 s, resulting in a coverage of about 0.006, 0.013 and 0.067 monolayers (ML), respectively. 1 ML is defined as a coverage of two CO molecules per (1 × 1) unit cell. Note that the pressure was slightly higher in the third experiment. The corresponding results are shown in Fig. 5. The calcite surface structure is in all images apparent in the form of four bright protrusions per (2 × 1) unit cell, having unit cell dimensions of 10 × 8.1 Å2 as indicated by white rectangles. The contrast in Fig. 5a, b and d represents a staggered (zig-zag) structure along the [42
] direction, which can be related to the different orientation of the two carbonate groups within the surface unit cell. Additionally, the (2 × 1) reconstruction is visible as a modulation along the [010] direction of the dark rows in all micrographs. As no row-pairing reconstruction is apparent, the pg surface symmetry is maintained.
Individual dark depressions, which are not observed on pristine calcite(10.4) surfaces, are visible and identified as single CO molecules (examples are marked by white arrows in Fig. 5). These depressions are well separated from each other, but are often manipulated while scanning (see Fig. 5c). We did not observe the formation of clusters or adsorbate islands. As apparent from the detailed micrograph in Fig. 5d, the molecules are located at the position of calcite surface features imaged bright in NC-AFM. Despite the sharp single band measured with IRRAS, we also observe molecules positioned slightly off-center the bright calcite surface feature (examples are marked by white circles in Fig. 5b). However, the difference between these orientations might be thermally averaged at the IRRAS measurement temperature of 62 K, while the barrier is too high at the NC-AFM imaging temperature of 5 K. Based on the IRRAS results we can furthermore conclude that the bright protrusions in these NC-AFM micrographs are representing the locations of the surface Ca2+ cations as the CO molecules are located at these positions. Two possibilities including a “positive” or “negative potential nanotip” imaging model were suggested for NC-AFM measurements of calcite(10.4) surfaces before.30 The imaging was explained to be dominated by electrostatic interactions, where either the oxygen or calcium sublattices are imaged as “bright” due to attraction with a positive or negative ion located at the tip apex, respectively.30 Thus, we conclude that the data herein were acquired with a probing tip with a negative termination.
Finally, we studied two different CO/calcite configurations by DFT, namely (i) a single CO molecule adsorbed within the (1 × 1) cell, corresponding to a coverage of 0.5 ML, and (ii) two CO molecules adsorbed within the (1 × 1) cell, corresponding to a complete monolayer coverage (1.0 ML). In all cases, CO was adsorbed on-top of a surface Ca site and initially positioned perpendicular to the surface with the carbon pointing downwards. This geometry corresponds to the common coordination of CO on metals and metal oxides.31 Optimization with DFT yielded a tilt of 14.42° for 0.5 ML and tilts of 14.57° and 15.14° for the two CO molecules within the (1 × 1) unit cell at 1 ML. The optimized geometries obtained for the different configurations are depicted in Fig. 6 and 7, according parameters are shown in Table 1. The slight increase in the tilt angle with coverage reveals the influence of lateral interactions between the CO molecules.32
![]() | ||
| Fig. 6 Side views of CO/calcite(10.4) after energy minimization for two different CO coverages of 0.5 ML (a) and 1.0 ML (b and c). CO tilt angles are included. | ||
| 0.5 ML | 1.0 ML | ||
|---|---|---|---|
| CO | CO | (CO)2 | |
| C–O/Å | 1.14 | 1.14 | 1.14 |
| x(C–Ca(surface))/Å | 0.04 | 0.06 | 0.05 |
| y(C–Ca(surface))/Å | 0.02 | 0.02 | 0.02 |
| z(C–Ca(surface))/Å | 2.78 | 2.78 | 2.79 |
| Ca(surface)–O(subsurface) | 2.33 | 2.33 | 2.33 |
| Θ/(deg) | 14.42 | 14.57 | 15.14 |
| Φ/(deg) | 36.49 | 36.18 | 38.34 |
The calculated adsorption energies and stretching frequencies of CO molecules at different coverages are reported in Table 2. For the isolated CO gas phase molecule, a stretch frequency of 2123 cm−1 was obtained, showing the frequency shift (Δν) of 20 cm−1 compared to the experimental value of 2143 cm−1.32 In the following, we will use the relative shifts to the free molecule, Δν, when comparing experimental and theoretical results.
For a single CO molecule at 0.5 ML coverage, the adsorption energy is calculated as −0.216 eV, and the CO stretching frequency is blue shifted by 27 cm−1 relative to the gas phase value. For the full monolayer, the binding energy is lowered to −0.210 eV and the blue shift is reduced to 21 cm−1. The lowering of binding energy is attributed to repulsive adsorbate–adsorbate interactions, as observed for CO adsorption on other oxide surfaces (e.g. TiO233). The absence of attractive adsorbate–adsorbate interactions also explains why no island formation is observed when increasing exposures. A simple rationalization of the changes in frequency between low and high coverage are difficult, since the overall blue shift (relative to the free molecule) is a result of a combination of several mechanisms: (1) the interaction between the surface electric field and the CO dipole moment (Stark effect), and (2) the repulsive potential (Pauli repulsion) due to the fact that the CO molecule vibrates against the rigid surface (also called wall effect).31,33 Yet, intermolecular repulsive interactions are in agreement with the NC-AFM data where no clustering of the CO molecules is observed.
DFT calculations for CO adsorbed on the (10.4) surface of calcite were carried out for two different surface coverages (0.5 ML and 1.0 ML). According to the calculations, CO adopts a slightly tilted geometry for the two different coverages studied here. In accord with the IRRAS experiments, a blue shift of the vibrational frequencies relative to the free molecule is obtained. The difference between experimental and theoretical shifts is less than 2 cm−1 for 0.5 ML, revealing a very good agreement. For 1.0 ML the difference between experimental and theoretical shifts increases to 12 cm−1. A missing vdW-potential in the DFT functional can partially explain the discrepancy. The theoretical value for the adsorption energy of −0.210 eV for 1.0 ML is good agreement with the experimental value of −0.31 eV. The binding energy decreased for higher coverages (1.0 ML relative to 0.5 ML), demonstrating the presence of repulsive interactions between adjacent CO adsorbates. This observation is in accord with the NC-AFM data indicating the absence of CO island formation.
An Atomistic Approach, J. Phys. Chem. B, 1998, 102(16), 2914–2922 CrossRef CAS;
(c)
R. Howie; J. Zussman and W. Deer, An introduction to rock-forming minerals. Longman, 1992 Search PubMed.
4), Langmuir, 2010, 26(11), 8295–8300 CrossRef PubMed;
(c) P. Fenter, P. Geissbühler, E. DiMasi, G. Srajer, L. Sorensen and N. Sturchio, Surface speciation of calcite observed in situ by high-resolution X-ray reflectivity, Geochim. Cosmochim. Acta, 2000, 64(7), 1221–1228 CrossRef CAS;
(d) S. L. Stipp and M. F. Hochella Jr, Structure and bonding environments at the calcite surface as observed with X-ray photoelectron spectroscopy (XPS) and low energy electron diffraction (LEED), Geochim. Cosmochim. Acta, 1991, 55(6), 1723–1736 CrossRef CAS.
4), Langmuir, 2010, 26(11), 8295–8300 CrossRef PubMed.Footnote |
| † These authors contributed equally to this work. |
| This journal is © the Owner Societies 2021 |