Dissociative and molecular oxygen chemisorption channels on reduced rutile TiO2(110): An STM and TPD study

High-resolution scanning tunneling microscopy (STM) and temperature-programmed desorption (TPD) were used to study the interaction of O₂ with reduced TiO₂(110)–(1 × 1) crystals. STM is the technique of choice to unravel the relation between vacancy and non-vacancy assisted O₂ dissociation channels as a function of temperature. It is revealed that the vacancy-assisted, first O₂ dissociation channel is preferred at low temperature (~ 120 K), whereas the non-vacancy assisted, second O₂ dissociation channel operates at temperatures higher than 150 K–180 K. Based on the STM results on the two dissociative O₂ interaction channels and the TPD data, a new comprehensive model of the O₂ chemisorption on reduced TiO₂(110) is proposed. The model explains the relations between the two dissociative and the molecular O₂ interaction channels. The experimental data are interpreted by considering the available charge in the near-surface region of reduced TiO₂(110) crystals, the kinetics of the two O₂ dissociation channels as well as the kinetics of the diffusion and reaction of Ti interstitials.     

The interaction between hexahydroxytriphenylene and the rutile TiO2(110)-(1 × 1) surface at UHV conditions

The interaction between 2,3,6,7,10,11-hexahydroxytriphenylene (HHTP) and the rutile TiO₂(110)–(1 × 1) surface under ultrahigh vacuum (UHV) conditions was investigated using X-ray photoemission spectroscopy (XPS), ultraviolet photoemission spectroscopy (UPS), near-edge X-ray absorption fine structure (NEXAFS) spectroscopy, and density functional theory (DFT) calculations. The NEXAFS results showed that HHTP molecules formed a submonolayer and a monolayer that aligned along the [001]-direction with, respectively, a more or less flat downward orientation and a more upright orientation to the TiO₂ surface. The HHTP molecules that aligned along the [001]-direction were most likely grafted onto the TiO₂(110) surface by a bidentate bridge between each of the oxygen atoms of one of the catechol units within the HHTP molecule and two adjacent Ti(5f)⁴⁺ ions on the TiO₂(110) surface. The coordination is non-dissociative in the case of the submonolayer, but dissociative in the monolayer, according to the analysis of the C1s XPS, UPS, C1s NEXAFS data and complementary DFT calculations.     

Methylene bromide chemistry and photochemistry on rutile TiO2(110)

The chemistry and photochemistry of methylene bromide (CD₂Br₂) on the rutile TiO₂(110) surface was probed using temperature programmed desorption (TPD). CD₂Br₂ desorbed in three desorption states at 145, 160 and 250 K tentatively assigned to desorption from the multilayer, from an η¹-CD₂Br₂ species and a bridging η²-CD₂Br₂ species, respectively. The latter two TPD states presumably involve binding of CD₂Br₂ molecules to the surface through Br coordination at five-coordinate Ti⁴⁺ surface sites. The 160 and 250 K TPD states saturated at coverages of 1.0 and 0.33 ML, respectively, where 1 ML is equivalent to the surface Ti⁴⁺ site density (5.2 × 10¹⁴ cm^? 2). No thermal decomposition of CD₂Br₂ was observed on either the clean surface or with preadsorbed O₂. UV irradiation of CD₂Br₂ on TiO₂(110) resulted in predominately photodesorption, with trace amounts of photodecomposition evidenced in TPD. The rate of CD₂Br₂ photodesorption from TiO₂(110) occurred with a low cross section (~ 2 × 10^? 21 cm²) similar to that expected from direct optical excitation of CD₂Br₂. This observation suggests that charge carriers generated in TiO₂(110) were no more effective in activating adsorbed CD₂Br₂ molecules than would be expected through direct molecular excitation. These findings suggest that photocatalytic destruction of halocarbons such as CD₂Br₂ on TiO₂ may preferentially occur though indirect processes (such as OH radical attack) as opposed to direct electron transfer processes involving charge carriers generated in TiO₂ by bandgap excitation.     

Electronic energy alignment at the PbSe quantum dots/ZnO(101̅0) interface

A number of solar energy conversion strategies depend on exciton dissociation across interfaces between semiconductor quantum dots (QDs) and other electron or hole conducting materials. A critical factor governing exciton dissociation and charge transfer in these systems is the alignment of electronic energy levels across the interface. We probe interfacial electronic energy alignment in a model system, sub-monolayer films of PbSe QDs adsorbed on single crystal ZnO(101?0) surfaces using ultraviolet photoemission spectroscopy. We establish electronic energy alignment as a function of quantum dot size and surface chemistry. We find that replacing insulating oleic-acid capping molecules on the QDs by the short hydrazine or ethanedithiol molecules results in pinning of the valence band maximum (VBM) of QDs to ZnO substrate states, independent of QD size. This is in contrast to similar measurements on TiO₂(110) where the alignment of the PbSe QD VBM to that of the TiO₂ substrate depends on QD size. We interpret these findings as indicative of strong electronic coupling of QDs with the ZnO surface but less with the TiO₂ surface. Based on the measured energy alignment, we predict that electron injection from the 1s_e level in photo-excited PbSe QDs to ZnO can occur with small QDs (diameter ? = 3.4 nm), but energetically unfavorably for larger dots (? = 6.7 nm). In the latter, hot electrons above the 1s_e level are necessary for interfacial electron injection.     

First principles and classical modeling of the oxidized titanium (0001) surface

The formation of a native oxide layer on the Ti(0001) surface is studied by means of FPMD simulations. In agreement with experimental findings, at low temperature we observe quick saturation of the surface reactivity, whereas incorporation of further O₂ molecules in the oxide network takes place after thermal annealing at ～ 800 K. At an O coverage of 2 ML the oxide layer presents an amorphous structure, an approximate TiO stoichiometry, and a broad distribution of Ti oxidation states from + 1 to + 4. We find consistency between the computed Bader atomic charges on Ti atoms and the correspondent charges computed classically by simple electrostatic minimization methods. On this basis we develop an analytic potential to simulate Ti/TiO_x interfaces, including only Coulomb interaction and short-range atomic repulsion terms. Two different parameter sets are proposed and their transferability among TiO₂ allomorphs and thin-layer oxide structures is tested both in static relaxations and in room-temperature MD simulations. These show only negligible changes in the topology of the oxide network after annealing and relaxation of the reference quantum model at the classical level. Moreover, we show that superficial oxide layers can be successfully generated purely classically by truncation of a large-scale Ti/amorphous–TiO₂ system. This may enable large-scale applications of our potential to tribology and biomolecular adsorption phenomena.     

Coarsening phenomena of metal nanoparticles and the influence of the support pre-treatment: Pt/TiO2(110)

One of the technologically most important requirements for the application of oxide-supported metal nanoparticles (NPs) in the fields of molecular electronics, plasmonics, and catalysis is the achievement of thermally stable systems. For this purpose, a thorough understanding of the different pathways underlying thermally-driven coarsening phenomena, and the effect of the nanoparticle synthesis method, support morphology, and degree of support reduction on NP sintering is needed. In this study, the sintering of supported metal NPs has been monitored via scanning tunneling microscopy combined with simulations following the Ostwald ripening and diffusion-coalescence models. Modifications were introduced to the diffusion-coalescence model to incorporate the correct temperature dependence and energetics. Such methods were applied to describe coarsening phenomena of physical-vapor deposited (PVD) and micellar Pt NPs supported on TiO₂(110). The TiO₂(110) substrates were exposed to different pre-treatments, leading to reduced, oxidized and polymer-modified TiO₂ surfaces. Such pre-treatments were found to affect the coarsening behavior of the NPs.No coarsening was observed for the micellar Pt NPs, maintaining their as-prepared size of ~ 3 nm after annealing in UHV at 1060 °C. Regardless of the initial substrate pre-treatment, the average size of the PVD-grown NPs was found to increase after identical thermal cycles, namely, from 0.5 ± 0.2 nm to 1.0 ± 0.3 nm for pristine TiO₂, and from 0.8 ± 0.3 nm to 1.3 ± 0.6 nm for polymer-coated TiO₂ after identical thermal treatments. Although no direct real-time in situ microscopic evidence is available to determine the dominant coarsening mechanism of the PVD NPs unequivocally, our simulations following the diffusion-coalescence coarsening route were in significantly better agreement with the experimental data as compared to those based on the Ostwald-ripening model. The enhanced thermal stability of the micellar NPs as compared to the PVD clusters might be related to their initial larger NP size, narrower size distribution, and larger interparticle distances.     

Adsorption of Xe atoms on the TiO2(1 1 0) surface: A density functional study

The interaction of xenon atoms with the TiO₂(1 1 0) surface of rutile has been studied by density functional theory methods. Five different possible adsorption sites on the relaxed and clean TiO₂(1 1 0) surface and on two different type of oxygen vacancies possible on this oxide substrate have been considered. In the case of the defect-free substrate, and when compared with a previous study concerning the adsorption of Ar atoms also on TiO₂(1 1 0), the xenon atom, as a larger and easier polarizable species, is shown to have a deeper physisorption well, as expected. Likewise, Xe atoms prefer to be bounded to positions nearby the outermost titanium atoms as found previously for Ar. This is in agreement with most studies concerning rare gases adsorbed on transition metal surfaces. In the case of the reduced surfaces, it is found that the interaction is more favourable in the protruding rows. The interaction is dominated by dispersion forces and DFT + dispersion energies are 3.5–5 times larger than the non-corrected DFT values and Xe-surface distances are smaller. Finally, an interesting correlation is obtained for the calculated interaction energies and the Xe–Ti distance.     

Computational study of the adsorption of dimethyl methylphosphonate (DMMP) on the (0 1 0) surface of anatase TiO2 with and without faceting

The adsorption of dimethyl methylphosphonate (DMMP) on the (0 1 0) surface of anatase TiO₂, which is isostructural with the (1 0 0), has been studied using density functional theory and two-dimensionally-periodic slab models. The experimentally-observed faceting of this surface has, for the first time, been included in the modeling. The relaxations of bare surfaces both with and without faceting are similar, leading to an atomic-scale roughening due to inward (outward) displacement of fivefold-coordinated Ti_5c (sixfold-coordinated T_6c) atoms together with outward displacement of threefold-coordinated O_3c atoms. Molecular adsorption occurs by formation of a Ti_5c?OP dative bond with one or more CH?O_2c bonds between CH₃ groups and unsaturated, twofold-coordinated (O_2c) sites. The energies for molecular adsorption, obtained using the B3LYP functional, are virtually identical (about ?21.0 kcal/mol) for the two surfaces and are also close to those found elsewhere for the rutile (1 1 0) and anatase (1 0 1) surfaces. A possible first step in the dissociative adsorption of DMMP has also been modeled and is found to be thermodynamically favored over molecular adsorption to a degree which depends on faceting.     

Enhanced performance of bi-layer Nb2O5 coated TiO2 nanoparticles/nanowires composite photoanode in dye-sensitized solar cells

Dye sensitized solar cells (DSSCs) were fabricated based on coumarin NKX-2700 dye sensitized bi-layer photoanode and quasi-solid state electrolyte sandwiched together with cobalt sulfide coated counter electrode. A novel bi-layer photoanode has been prepared using composite mixtures of 90 wt.% TiO₂ nanoparticles + 10 wt.% TiO₂ nanowires (TNPWs) as active layer and Nb₂O₅ is coated on the active layer, which acts as scattering layer. Hafnium oxide (HfO₂) was applied over the TNPWs/Nb₂O₅ photoanode film, as a blocking layer. TiO₂ nanoparticles (TNPs), TiO₂ nanowires (TNWs) and TNPWs/Nb₂O₅ were characterized by X-ray diffractometer (XRD), scanning electron microscope (SEM) and transmission electron microscope (TEM). The sensitizing organic dye coumarin NKX-2700 displayed maximum absorption wavelength (λ_max) at 525 nm, which could be observed from the UV–vis spectrum. DSSC-1 fabricated with composite bi-layer photoanode revealed enhanced photo-current efficiency (PCE) as compared to other DSSCs and illustrated photovoltaic parameters; short-circuit current J_SC = 18 mA/cm², open circuit voltage (V_OC) = 700 mV, fill factor (FF) = 64% and PCE (η) = 8.06%. The electron transport and charge recombination behaviors of DSSCs were investigated by electrochemical impedance spectra (EIS) and the results illustrated that the DSSC-1 showed the lowest charge transport resistance (R_tr) and the longest electron lifetime (τ_eff). Therefore, in the present investigation, it could be concluded that the novel bi-layer photoanode with blocking layer increased the short circuit current, electron transport and suppressed the recombination of charge carriers at the photoanode/dye/electrolyte interface in DSSC-1.     

Electronic charge transfer between Au nano-particles and TiO2-terminated SrTiO3(0 0 1) substrate

The interaction between Au nano-particles and oxide supports is recently discussed in terms of the catalytic activities. This paper reports the electronic charge transfer between Au nano-particles and TiO₂-terminated SrTiO₃(0 0 1) substrate, which is compared with that for stoichiometric(S)-, pseudo-stoichiometric(S¹)- and reduced(R)-TiO₂(1 1 0) supports. We observed the photoelectron spectra of Au 4f, O 2s, Ti 3p, and Sr 4p lines and also measured the work functions for Au/oxides supports using synchrotron-radiation light. As the results, all the O 2s, Ti 3p, and Sr 4p lines for Au/SrTiO₃(0 0 1) show lower binding energy shifts in a quite same manner and abrupt increase in the work function is seen in an initial stage. This clearly evidences an electronic charge transfer from the substrate to Au probably due to a much larger work function of Au than SrTiO₃(0 0 1), which leads to an upward band bending (0.3 eV) just like a Schottky contact. Electronic charge transfers also take place at Au/S- and Au/S¹-TiO₂(1 1 0) and Au/R-TiO₂(1 1 0) interfaces, where electrons are transferred from Au to S- and S¹-TiO₂ and from R-TiO₂ to Au, as predicted by ab initio calculations.     

First-principles molecular dynamics simulations of uranyl ion interaction at the water/rutile TiO2(110) interface

The effects of temperature and solvation on uranyl ion adsorption at the water/rutile TiO₂(110) interface are investigated by Density Functional Theory (DFT) in both static and Born–Oppenheimer molecular dynamics approaches. According to experimental observations, uranyl ion can form two surface complexes in a pH range from 1.5 to 4.5. Based on these observations, the structures of the complexes at 293 K are first calculated in agreement with vacuum static calculations. Then, an increase in temperature (293 to 425 K) induces the reinforcement of uranyl ion adsorption due to the release of water molecules from the solvation shell of uranyl ion. Finally, temperature can modify the nature of the surface species.     

Ab initio study of shear strain effects on ferroelectricity at PbTiO3 thin films

Ferroelectric thin film with the perovskite ABO₃ structure have been widely used in technology applications, e.g., actuators in MEMS/NEMS and nonvolatile random access memories (FeRAM). In order to clarify the effect of the shear strain on the ferroelectricity, the PbTiO₃ thin film as a typical one is chosen. The focus of this study is to put on the PbO-terminated (1? × ?1) and c(2? × ?2) surfaces and the TiO₂-terminated (1? × ?1) surface. Based on ab initio density functional theory calculations with the local density approximation, we have found out that in both the PbO and TiO₂-terminated (1? × ?1) models, the ferroelectricity in the PbO layers was enhanced under the positive shear strain while it was suppressed under the negative one. For the TiO₂ layers, the ferroelectricity was slightly enhanced and sharply suppressed under the positive and negative shear strains, respectively. In the PbO-terminated (2? × ?2) model, the AFE phase was suppressed by the FE phase under the positive shear strain while the opposite trend was found under the negative shear strain. For the PbO layers, the ferroelectricity was enhanced under the positive and negative shear strains. For the TiO₂ layers, the influence of the negative shear strain on the ferroelectricity was larger than that of the positive one. In addition, the ideal strength of the PbTiO₃ thin film with the different terminations was investigated as well.     

TiO2(110)-(1 × 1) supported Cu particles: An FT-RAIRS investigation

The morphology of TiO₂(110)-(1 × 1) supported Cu particles has been investigated by Fourier Transform Reflection Absorption Infrared Spectroscopy (FT-RAIRS), employing adsorbed CO as a probe molecule sensitive to local surface structure. For Cu coverage (deposited at 300 K) less than 2.85 MLE nucleated Cu particles in the range 2 nm–4 nm are formed, as indicated by a final state shift in the core level Cu(2p_3/2) binding energy and by the existence of only transmission bands in the FT-RAIRS spectra for adsorbed CO. ν_S(CO) indicates that these small particles expose sites similar to those of the stepped Cu surfaces Cu(211), Cu(311), and Cu(755). At Cu coverages in the range of 6 MLE and above, corresponding to particle sizes above 4.6 nm, ν_S(CO) indicates the predominance of (110), (100) and (111) adsorption sites. Annealing the Cu layers to 650 K results in the slight growth of the particle sizes, and transformation of the CO adsorption sites corresponding to the close packed facets. The transformation of the local dielectric from that of titania to that dominated by the Cu particle is shown to take place between 3.7 and 4.2 nm, and this change is also to a smaller extent sensitive to the dispersion of the particles.     

Photochemistry on TiO2: Mechanisms behind the surface chemistry

Photochemistry from TiO₂ surfaces is described for two cases: The UV-induced photodesorption of O₂ from TiO₂(1 1 0) – 1 × 1; and the hydrophilic effect caused by UV irradiation on TiO₂. In both cases fundamental information about how these processes occur has been found. In the case of the O₂ photodesorption kinetics, it has been found that the rate of the process is proportional to the square root of the UV flux, showing that second-order electron–hole pair recombination is dominant in governing the photodesorption rate. In addition these measurements provide an estimate of the concentration of hole traps in the TiO₂ crystal. In other measurements of the UV-induced hydrophilicity, starting with the atomically-clean TiO₂ surface, it has been shown that the effect occurs suddenly at a critical point during irradiation as a result of photooxidation of a monolayer of hydrocarbon (n-hexane) at equilibrium with ppm concentration of n-hexane in O₂ at 1 atmosphere pressure.     

Periodic density functional theory studies of the VOx/TiO2 (anatase) catalysts: Structure and stability of monomeric species

Periodic density functional theory has been utilized to investigate the structure and stability of monomeric HVO_x species on anatase support. The three most stable surfaces of anatase were investigated, namely the (001), (100) and (101) surfaces. Unlike previous theoretical studies it was found that on the (001) surface vanadia species with five-coordinated vanadium atom are more stable than those with tetrahedrally coordinated vanadium atom. On the other hand, on the (100) and (101) surfaces, the vanadium atom in the vanadia species is still tetrahedrally coordinated. The stability of different VO_x/TiO₂ structures which are not fully dehydrated has been systematically studied and the results show that the vanadia species on the three surfaces follow an order of TiO₂ (001) > TiO₂ (100) > TiO₂ (101). This can be understood from the acidity and basicity of the three anatase surfaces. The results suggest that monomeric VO_x species may be better stabilized if the support exposes more (001) surfaces. Our analyses on electronic structure of the most stable VO_x/TiO₂ structure (D₀₀₁) suggest that its bridging V–O–Ti oxygen atoms may have higher reactivity than the terminal vanadyl oxygen atoms.     

Acute effects of sono-activated photocatalytic titanium dioxide nanoparticles on oral squamous cell carcinoma

Sonodynamic therapy (SDT) is a new treatment modality using ultrasound to activate certain chemical sensitizers for cancer therapy. In this study, effects of high intensity focused ultrasound (HIFU) combined with photocatalytic titanium dioxide (TiO₂) nanoparticles on human oral squamous cell line HSC-2 were investigated. Viability of HSC-2 cells after 0, 0.1, 1, or 3 s of HIFU irradiation with 20, 32, 55 and 73 W cm⁻² intensities in the presence or absence of TiO₂ was measured immediately after the exposures in vitro. Immediate effects of HIFU (3 s, 73 W cm⁻²) combined with TiO₂ on solid tumors were also examined by histological study. Cytotoxic effect of HIFU + TiO₂ in vitro was significantly higher than that of TiO₂ or HIFU alone with the tendency to increase for higher HIFU intensity, duration, and TiO₂ concentration in the suspension. In vivo results showed significant necrosis and tissue damage in HIFU and HIFU + TiO₂ treated samples. However, penetration of TiO₂ nanoparticles into the cell cytoplasm was only observed in HIFU + TiO₂ treated tissues. In this study, our findings provide a rational basis for the development of an effective HIFU based sonodynamic activation method. This approach offers an attractive non-invasive therapy technique for oral cancer in future.     

TiO2 chemical vapor deposition on Si(111) in ultrahigh vacuum: Transition from interfacial phase to crystalline phase in the reaction limited regime

The interaction between the metal organic precursor molecule titanium(IV) isopropoxide (TTIP) and three different surfaces has been studied: Si(111)-(7 × 7), SiO_x/Si(111) and TiO₂. These surfaces represent the different surface compositions encountered during TTIP mediated TiO₂ chemical vapor deposition on Si(111). The surface chemistry of the titanium(IV) isopropoxide precursor and the film growth have been explored by core level photoelectron spectroscopy and x-ray absorption spectroscopy using synchrotron radiation. The resulting film morphology has been imaged with scanning tunneling microscopy. The growth rate depends on both surface temperature and surface composition. The behavior can be rationalized in terms of the surface stability of isopropoxy and isopropyl groups, confirming that growth at 573 K is a reaction limited process.     

Preparation and characterization of low density polystyrene/TiO2 core–shell particles for electronic paper application

In order to reduce the density mismatch between TiO₂ and the low dielectric medium and improve the dispersion stability of the electrophoretic particles in the low dielectric medium for electrophoretic display application, polystyrene/titanium dioxide (PS/TiO₂) core–shell particles were prepared via in-situ sol–gel method by depositing TiO₂ on the PS particle which was positively charged with 2-(methacryloyloxy)ehyl trimethylammonium chloride (DMC). The morphology and average particle size of PS/TiO₂ core–shell particles were observed by transmission electron microscopy (TEM), scanning electron microscope (SEM) and particle size analyzer. It was found that density of PS/TiO₂ core–shell particles were reduced obviously and the particles can suspend in the low dielectric medium of low density. The PS/TiO₂ core–shell particles can endure ultrasonic treatment because of the interaction between TiO₂ and PS. Zeta potential and electrophoretic mobility of the fabricated core–shell particles in a low dielectric medium with charge control agent was measured to be −44.3 mV and −6.07 × 10⁻⁶ cm²/Vs, respectively, which presents potential in electronic paper application.     

Theoretical study of the (110) surface of Sn1 - xTixO2 solid solutions with different distribution and content of Ti

The composition and thermodynamic stability of the (110) surface of Sn_1 - xTi_xO₂ rutile solid solutions was investigated as a function of Ti-distribution and content up to the formation of a full TiO₂ surface monolayer. The bulk and (110) surface properties of Sn_1 - xTi_xO₂ were compared to that of the pure SnO₂ and TiO₂ crystal. A large supercell of 720 atoms and a localized basis set based on the Gaussian and plane wave scheme allowed the investigation of very low Ti-content and symmetry. For the bulk, optimization of the crystal structure confirmed that up to a Ti-content of 3.3 at.%, the lattice parameters (a, c) of SnO₂ do not change. Increasing further the Ti-content decreased both lattice parameters down to those of TiO₂. The surface energy of these solid solutions did not change for Ti-substitution in the bulk of up to 20 at.%. In contrast, substitution in the surface layer rapidly decreased the surface energy from 0.99 to 0.74 J/m² with increasing Ti-content from 0 to 20 at.%. As a result, systems with Ti atoms distributed in the surface (surface enrichment) had always lower energies and thus were thermodynamically more favorable than those with Ti homogeneously distributed in the bulk. This was attributed to the lower energy necessary to break the TiO bonds than SnO bonds in the surface layer. In fact, distributing the Ti atoms homogeneously or segregated in the (110) surface led to the same surface energy indicating that restructuring of the surface bond lengths has minimal impact on thermodynamic stability of these rutile systems. As a result, a first theoretical prediction of the composition of Sn_1 - xTi_xO₂ solid solutions is proposed.     

Structural and electronic properties of cobalt carbide Co2C and its surface stability: Density functional theory study

Density functional theory calculations have been performed to investigate the structural and electronic properties of bulk Co₂C and the stability of low index Co₂C surfaces. We found that the formation of Co₂C is exothermic with the formation energy of ? 0.81 eV/Co₂C with respect to Co under the presence of syngas (mixture of CO and H₂). While formed Co₂C can be decomposed further to metal Co and graphite carbon with modest energy gain of 0.37 eV/Co₂C. This suggests that Co₂C is only metastable in Fischer–Tropsch synthesis, which agrees well with experimental findings. The density of states (DOSs) reveals that the Co₂C is paramagnetic and strong metallic-like. The difference of charge density analysis indicates that the bond of Co₂C is of the mixtures of metallic, covalent, and ionic properties. A variety of low index Co₂C surfaces with different terminations are studied. We find that the surface energy of low index stoichiometric Co₂C highly relies on the surface area, the number of coordination of surface atoms and the surface dipole, with the decreased stability order of (101) > (011) > (010) > (110) > (100) > (001) = (111). Our results indicate that under Co-poor condition, the formation of non ? stoichiometric surface (011) and (111) without terminated cobalt is energetically more favorable, while under Co-rich condition the formation of non ? stoichiometric (111) surface with cobalt overlayer are preferential.