Non-fullerene electron acceptors have attracted enormous attention of the research community owing to their advantages of optoelectronic and chemical tunabilities for promoting high-performance polymer solar cells (PSCs). Among them, fused-ring electron acceptors (FREAs) are the most popular ones with the good structural planarity and rigidity, which successfully boost the power conversion efficiencies (PCEs) of PSCs to over 14%. In considering the cost-control of future scale-up applications, it is also worthwhile to explore novel structures that are easy to synthesize and still maintain the advantages of FREAs. In this work, we design and synthesize a new electron acceptor with an unfused backbone, 5, 5'-((2, 5-bis((2-hexyldecyl)oxy)-1, 4-phenylene)bis(thiophene-2-yl))bis(methanylylidene)) bis(3-oxo-2, 3-dihydro-1H-indene-2, 1-diylidene))dimal-ononitrile (ICTP), which contains two thiophenes and one alkoxy benzene as the core and 2-(3-oxo-2, 3-dihydroinden-1-ylidene) malononitrile (IC) as the terminal groups. The synthetic route to ICTP involves only three steps, with high yields. Density functional theory calculations indicate that the non-covalent interactions, O…H and O…S, help reinforce the space conformation between the central core and the terminals. ICTP shows broad and strong absorption in the long-wavelength range between 500 and 760 nm. The highest occupied molecular orbital and lowest unoccupied molecular orbital levels of ICTP were measured to be -5.56 and -3.84 eV by cyclic voltammetry. The suitable absorption and energy levels make ICTP a good acceptor candidate for medium bandgap polymer donors. The best devices based on PBDB-T:ICTP showed a PCE of 4.43%, with an open circuit voltage (V_OC) of 0.97 V, a short circuit current density (J_SC) of 8.29 mA∙cm^-2, and a fill factor (FF) of 0.55, after adding 1% 1, 8-diiodooctane (DIO) as the solvent additive. Atomic force microscopy revealed that DIO could ameliorate the strong aggregation in the blended film and lead to a smoother film surface. The hole and electron mobilities of the optimized device were measured to be 9.64 and 2.03 × 10^-5 cm²∙V^-1∙s^-1, respectively, by the space-charge-limited current method. The relatively low mobilities might be responsible for the moderate PCE. Further studies can be performed to enlarge the conjugation length by including more aromatic rings. This study provides a simple strategy to design non-fullerene acceptors and a valuable reference for the future development of PSCs.

Recently, non-fullerene polymer solar cells (NPSCs) have been developed rapidly because of the flexible energy-level variability and excellent optical absorption properties of non-fullerene electron acceptors. Among them, fused-ring electron acceptors (FREAs) with acceptor-donor- acceptor (A-D-A) structures have been extensively exploited in high-performance NPSCs. These FREAs often employ central aromatic fused rings attached to several rigid side-chains and flanked by two electron-deficient terminals. Many efforts have focused on the modification of the central flat conjugated backbone in order to gain broad and strong absorption and dense stacking. However, the preparation of such FREAs is relatively complex, especially for large fused-ring structures. In a previous work, we provided a simple and useful method to extend the effective conjugation length and broaden the absorption spectrum of the acceptor by noncovalent intramolecular interactions. On this basis, in this work, we have designed and synthesized a new A-D-A-type FREA (ITOIC-2Cl) that uses 4, 9-dihydro-s-indaceno[1, 2-b:5, 6-b']dithiophene (IDT) as a central donor unit, bis(alkoxy)-substituted thiophene rings as conformational locking π-bridges between the donor and acceptor units, and cyanoindanones modified with two high-electron-affinity chlorine atoms as end-capping acceptor units. On one hand, we can attain good backbone planarity in the solid state via the noncovalent conformational locking induced by sulfur−oxygen (S···O) and oxygen−hydrogen (CH···O) interactions, which are not strong enough to lock the coplanar conformation in solution, thus simultaneously endowing ITOIC-2Cl with good solubility. On the other hand, we can enhance the intramolecular charge transfer by enhancing the electron deficiency of the terminal groups. The optical and electrochemical properties of ITOIC-2Cl were systematically explored. Moreover, in combination with the donor polymer of [(2, 6-(4, 8-bis(5-(2-ethylhexyl)thiophen-2-yl)-benzo[1, 2-b:4, 5-b']dithiophene))-alt-(5, 5-(1', 3'-di-2-thienyl-5', 7'-bis(2-ethylhexyl)benzo[1', 2'-c:4', 5'-c']dithiophene-4, 8-dione))] (PBDB-T), the photovoltaic performances of the devices and the corresponding blend morphologies were studied. ITOIC-2Cl exhibited a broad absorption spectrum up to 900 nm, which is beneficial for broad harvesting of photons across the visible and NIR region. The PBDB-T:ITOIC-2Cl-based blend films exhibited favorable fibrous nanostructures with appropriate nanoscale phase separation, verified by atomic force microscopy and transmission electron microscopy characterizations. This morphology is beneficial for charge transport. Through the space-charge-limited current measurement, the PBDB-T:ITOIC-2Cl-based device exhibited the high hole/electron mobility of 1.85 × 10⁻⁴/1.19 × 10⁻⁴ cm²∙V⁻¹∙s⁻¹. The PBDB-T:ITOIC-2Cl-based devices obtained a high power conversion efficiency of 9.37%, with an open-circuit voltage (V_oc) of 0.886 V, short-circuit current (J_sc) of 17.09 mA cm⁻², and a fill factor (FF) of 61.8%. These results thus demonstrate the efficacy of the proposed strategy for designing high-performance non-fullerene FREAs.

Ternary blends have been considered as an effective approach to improve power conversion efficiency (PCE) of organic solar cells (OSCs). Among them, the fullerene-containing ternary OSCs have been studied extensively, and their PCEs are as high as over 14%. However, all non-fullerene acceptor ternary OSCs are still limited by their relatively lower PCEs. In this work, we used wide-bandgap benzodithiophene-difluorobenzotriazole copolymer FTAZ as the donor, low-bandgap fused-ring electron acceptor (FREA), fused tris(thieno- thiophene) end-capped by fluorinated 1, 1-dicyanomethylene-3-indanone (FOIC) as acceptor, and two medium-bandgap FREAs, indaceno-dithiophene end- capped by 1, 1-dicyanomethylene-3-indanone (IDT-IC) and indacenodithiophene end-capped by 1, 1-dicyanomethylene-3-benzoindanone (IDT-NC), as the third components to fabricate the ternary blends FTAZ:FOIC:IDT-IC and FTAZ:FOIC:IDT-NC, and investigated the effects of the third components on the performance of ternary OSCs. Both IDT-IC and IDT-NC are based on the same indacenodithiophene core but contain different terminal groups (phenyl and naphthyl). Relative to IDT-IC with phenyl terminal groups, IDT-NC with naphthyl terminal groups has extended π-conjugation, down-shifted lowest unoccupied molecular orbital (LUMO), red-shifted absorption and higher electron mobility. The binary devices based on the FTAZ:FOIC, FTAZ:IDT-IC and FTAZ:IDT-NC blends exhibit PCEs of 9.73%, 7.48% and 7.68%, respectively. Compared with corresponding binary devices, both ternary devices based on FTAZ:FOIC:IDT-IC and FTAZ:FOIC:IDT-NC exhibit better photovoltaic performances. When the IDT-IC weight ratio in acceptors is 50%, the FTAZ:FOIC:IDT-IC ternary devices exhibit the best PCE of 11.2%. The ternary-blend OSCs yield simultaneously improved open-circuit voltage (V_OC), short-circuit current density (J_SC) and fill factor (FF) compared with the binary devices based on FTAZ:FOIC. The higher V_OC is attributed to the higher LUMO energy level of IDT-IC compared with FOIC. The improved J_SC is attributed to the complementary absorption of FOIC and IDT-IC. The introduction of IDT-IC improves blend morphology and charge transport, leading to higher FF. The FTAZ:FOIC:IDT-NC system yields a higher PCE of 10.4% relative to the binary devices based on FTAZ:FOIC as the active layer. However, the PCE of the FTAZ:FOIC:IDT-NC-based ternary devices is lower than that of the FTAZ:FOIC:IDT-IC-based ternary devices. Compared with the binary devices based on FTAZ:FOIC, in FTAZ:FOIC:IDT-NC-based ternary devices, as the ratio of the third component increases, the V_OC increases due to the higher LUMO energy level of IDT-NC, the FF increases due to optimized morphology and improved charge transport, while the J_SC decreases due to the overlapped absorption of FOIC and IDT-NC. The terminal groups in the third components affect the performance of the ternary OSCs. The lower LUMO. energy level of IDT-NC is responsible for the lower V_OC of the FTAZ:FOIC:IDT-NC devices. The red-shifted absorption of IDT-NC leads to the overlapping of the absorption spectra of IDT-NC and FOIC and lower J_SC. On the other hand, replacing the phenyl terminal groups by the naphthyl terminal groups influences the π-π packing and charge transport. The FTAZ:FOIC:IDT-NC blend exhibits higher electron mobility and more balanced charge transport than FTAZ:FOIC:IDT-IC, leading to a higher FF.

Recently, non-fullerene small molecular acceptors (NFSMAs) have received great attention because of their broad and strong absorption spectra and stable active layer morphology when compared with traditional fullerene acceptors. The most widely used strategy to design NFSMAs is through A-D-A type molecules, in which an electron-rich core unit (D) is flanked by two electron-deficient units (A). In order to fine-tune the absorption spectra, energy levels, and photovoltaic properties of NFSMAs, great efforts have been made to modify the conjugated backbone of A-D-A type molecule acceptors. In a previous work, we developed a small molecular electron acceptor, namely MBN-Ph, with an A-D-A structure and an organoboron core unit. MBN-Ph exhibited distinctive absorption spectra with two absorption bands in short- and long-wavelength regions. It is known that side chains or substituents on small molecular electron acceptors can also play an important role in the molecular properties and photovoltaic performance of bulk heterojunction organic solar cells (OSCs). In this work, we report an A-D-A type organoboron compound (MBN-Th) bearing a thienyl substituent on the boron atom, which can be used as an electron acceptor for OSCs. The lowest unoccupied molecular orbital (LUMO) of MBN-Th delocalized on the entire backbone, while the highest occupied molecular orbital (HOMO) localized on the core unit. The unique electronic structure of MBN-Th resulted in two strong absorption peaks at 490 and 726 nm, which indicate a wide absorption spectrum and superior sunlight harvesting capability. Compared with the phenyl substituent, the thienyl group led to an unchanged LUMO energy level, low-lying HOMO energy level by 0.1 eV, and blue-shifted absorption spectrum by 20 nm. OSCs with MBN-Th as an electron acceptor showed a power conversion efficiency of 4.21% and a wide photoresponse from 300 to 850 nm. Our results indicate that the substitution of the boron atom with a thienyl group is an effective strategy to tune the electronic structure of organoboron compounds for applications as electron acceptors in OSCs.

The development of the photocatalytic production of hydrogen from water splitting has attracted immense attention in recent years. CdS is a potential photocatalyst with a visible light response, though it still suffers from a limited activity for hydrogen production due to the fast recombination of photo-induced electron/hole pairs and the low reaction rate of hydrogen evolution on the surface. Studies on the effect of CdS surface structure and properties on hydrogen production are still very limited. In this work, we prepared three CdS nanocrystals with different morphologies: long rod, short rod, and triangular plate. The prepared samples were well characterized by transmission electron microscopy (TEM), X-ray diffraction (XRD), Brunauer-Emmett-Teller (BET) specific surface area analysis, X-ray photoelectron spectroscopy (XPS), photoluminescence (PL) spectroscopy, and UV-Vis diffuse reflectance spectroscopy (UV-Vis DRS). From the results of TEM, XRD and XPS, we find that the three CdS nanocrystals with different morphologies were successfully synthesized. From the PL spectra, we conclude that the area of exposed nonpolar surface and degree of surface defects increase with an increase in aspect ratio. We also performed the photocatalytic hydrogen production reaction using the three CdS crystals. Long rod-like CdS (lr-CdS) exhibits the highest photocatalytic activity, with a hydrogen production rate of 482 μmol·h^-1·g^-1, which is 2.6 times that of short rod-like CdS (sr-CdS) (183 μmol·h^-1·g^-1) and 8.8 times that of triangular plate-like CdS (tp-CdS, 55 μmol h^-1·g^-1). It is found that lr-CdS shows a higher hydrogen production rate than sr-CdS and tp-CdS. We find that the hydrogen production rate is related to the degree of surface defects. Surface defects can trap the photo-induced electrons/holes, thus decreasing their probability of recombination. In addition, these defects can be used to anchor Pd particles to form a heterojunction structure that facilitates the separation of photo-induced charges. Therefore, we also compared three CdS/Pd nanocrystals synthesized with the three abovementioned morphologies with respect to hydrogen production. With 1% (w, mass fraction) Pd, the hydrogen production rate was greatly enhanced compared to all the CdS catalysts. Compared to the unpromoted CdS, the reaction rate is enhanced 43.1, 10.7 and 6.0 times over those of sr-CdS, lr-CdS and tp-CdS, respectively. Notably, the hydrogen production rate with short rod-like CdS/Pd reaches 7884 μmol·h^-1·g^-1, which can be favorably compared with the ever-increasing values reported in the literature. Hopefully, this work provides knowledge on the effect of crystal surface structure and properties on photocatalysis.

Over the past decade, significant progress has been made in theoretical and experimental research in the field of chemical reaction dynamics, moving from triatomic reactions to larger polyatomic reactions. This has challenged the theoretical and computational approaches to polyatomic reaction dynamics in two major areas: the potential energy surface and the dynamics. Highly accurate potential energy surfaces are essential for achieving accurate dynamical information in quantum dynamics calculations. The increased number of degrees of freedom in larger systems poses a significant challenge to the accurate construction of potential energy surfaces. Recently, there has been substantial progress in the development of potential energy surfaces for polyatomic reactive systems. In this article, we review the recent developments made by our group in constructing highly accurately fitted potential energy surfaces for polyatomic reactive systems, based on a neural network approach. A key advantage of the neural network approach is its more faithful representation of the ab initio points. We recently proposed a systematic procedure, based on neural network fitting, for the construction of accurate potential energy surfaces with very small root mean square errors. Based on the neural network approach, we successfully developed potential energy surfaces for polyatomic reactions in the gas phase, including the reactive systems OH₃, HOCO, and CH₅, and the dissociation of gas-phase molecules on metal surfaces, such as H₂O on the Cu(111) surface. These potential energy surfaces were fitted to an unprecedented level of accuracy, representing the most accurate potential energy surfaces calculated for these systems, and were rigorously tested using quantum dynamics calculations. The quantum dynamics calculations based on these potential energy surfaces produce accurate results, which are in good agreement with experiments. We have also proposed a new method for developing permutationally invariant potential energy surfaces, named fundamental-invariant neural networks. Mathematically, fundamental invariants are used to finitely generate the permutation-invariant polynomial ring; thus, fundamental-invariant neural networks can approximate any function to arbitrary accuracy. The use of fundamental invariants minimizes the size of the input permutation-invariant polynomials, which reduces the evaluation time for potential energy calculations, especially for polyatomic systems. Potential energy surfaces for OH₃ and CH₄ were constructed using fundamental-invariant neural networks, with their accuracies confirmed by full-dimensional quantum dynamics and bound-state calculations. These developments in the construction of highly accurate potential energy surfaces are expected to extend the theoretical study of reaction dynamics to larger and more complex systems.

The structure of low-temperature α-Cu₂Se, which is of great importance for understanding the mechanism of the significant increase in thermoelectric performance during the α-β phase transition of Cu₂Se, has still not been fully solved. Because it is restricted by the quality of polycrystal and powder specimens and the accuracy of characterization methods such as the conventional transmission electron microscopy (TEM) and X-ray diffraction (XRD), direct observation with atomic-scale resolution to reveal the structural details has not been realized, although electron diffraction and high-resolution transmission electron microscopy (HRTEM) studies have indicated the complexity of the α-Cu₂Se layered structure. Owing to developments in the focused ion beam (FIB) milling preparation method, high-quality single crystalline specimens with specific crystallographic orientations can be prepared to ensure that atomic-resolution images along a specific orientation can be acquired. Furthermore, the developments in aberration correction technology in TEM and scanning transmission electron microscopy (STEM) allow us to observe the subtle details of structural variation and evolution. Herein, we report, for the first time, the atomic-resolution high-angle annular dark field (HAADF) images acquired along the $ {{\left[ \bar{1}\bar{1}2 \right]}_{\text{c}}}$ axis of α-Cu₂Se using spherical-aberration (C_s)-corrected STEM from FIB-prepared single crystalline specimens. The observations revealed that the complex structure is generated by ordered fluctuations of Se atoms with various forms, including that some of the Se atoms on the two sides of the Cu deficiency layer get closer to each other than the others and the neighboring Cu deficiency layers have different forms of ordered Se fluctuations. These characteristics can only be observed along the $ {{\left[ \bar{1}\bar{1}2 \right]}_{\text{c}}}$ axis, while these details were not visible in a previous study along the $ <\bar{1}10{{>}_{\text{c}}} $ axis or in our results obtained along the ${{\left[ \bar{1}01 \right]}_{\text{c}}} $ axis. By combining the electron diffraction patterns, several models of the unit cell variants were established, including the two-layer and four-layer cells (both have two different shapes) and the two-layer variants with and without central symmetry. These variants can also transform into each other, and an α-Cu₂Se crystal can be formed through the random assembly of these variants. Using the program QSTEM, the corresponding HAADF images of these variants were simulated. The simulation results were similar to the experimental HAADF images and reflected most of the observed details, including the different forms of the ordered fluctuations of Se atoms and the dispersion of Cu atoms, which indicates that our structure models of α-Cu₂Se are reasonable. This work provides new critical information for thoroughly understanding the structure of α-Cu₂Se and the α-β phase transition of Cu₂Se.

In recent years, organic solar cells (OSCs) have attracted increasing attention, and the power conversion efficiency (PCE) of OSCs has markedly improved. To enhance the photovoltaic properties of OSCs, it is important to develop the donor materials in the light-harvesting layer, which mainly include conjugated polymers and small molecules (SMs). Compared with polymeric materials, small-molecule materials have been widely investigated for their superior characteristics, such as well-defined molecular structures that can provide good batch-to-batch reproducibility. In this work, we synthesized three SM donor materials with theacceptor-donor-acceptor (A-D-A) structure by employing the trialkylthienyl-substituted benzodithiophene (TriBDT-T) unit as the D-core unit, and rhodanine (RN), cyano-rhodanine (RCN), and 1, 3-indanone (IDO) as the A end groups, respectively. The optical properties, molecular energy levels, and thermogravimetic characteristics of the three SMs were studied; moreover, the blend morphologies and photovoltaic properties of the devices by employing the non-fullerene (NF) acceptor, IT-4F, were systematically investigated. The results showed that 1) the three SMs exhibit good thermal stabilities as evinced by thermogravimetric analysis (TGA), and all decomposition temperatures exceeded 410 ℃; 2) They all exhibit strong and broad absorption in the visible light range (300–700 nm), and show similar molar extinction coefficients; 3) the HOMO levels are -5.47 eV, -5.54 eV, and -5.44 eV for RN, RCN, and IDO, respectively, implying the clear influence of the different end groups for the energy levels of the A-D-A-type SMs; the slight differences in the optical and electrochemical properties of the corresponding donor material could be attributed to the different electron-withdrawing ability of the A-type end groups. When studying the photovoltaic properties, interestingly, the RN:IT-4F blend was found to form fibrillar-like aggregates with appropriate size, and the corresponding devices exhibited desirable short circuit current (J_sc) and thus the highest PCE value of 9.25%; however, large-size aggregates were formed in the RCN:IT-4F and IDO:IT-4F blend films, resulting in a much lower J_sc and fill factor (FF), and the PCE of the corresponding devices were only around 6%. In summary, by introducing RN, RCN, and IDO as the A units, we synthesized a class of TriBDT-T based A-D-A type SMs. This study shows that terminal A units exert influence on the absorption spectra, molecular energy level, and morphologies after blending with the acceptor material, and hence, the corresponding devices exhibit significant differences in photovoltaic performance. This work also provides useful information for the molecular design of SM donor materials.

Nonadiabatic processes widely exist in photochemistry and photophysics. The theoretical treatment of nonadiabatic processes is an important challenge in theoretical chemistry. In nonadiabatic dynamics, the well-known "Born-Oppenheimer approximation" breaks down due to the involvement of strong nuclear-electronic coupled motions. Hence, the development of a theoretical framework is required for the proper treatment of nonadiabatic dynamics. The fewest-switches trajectory surface-hopping method developed by Tully is one of the most widely used methods in the treatment of nonadiabatic processes because of its rather simple numerical implementation. In this approach, the nuclear degrees of freedom are propagated on the potential energy surface of an electronic state using the classical equations of motion, while the electronic degrees of freedom are propagated along the same trajectory according to the time-dependent Schr?dinger equation. Nonadiabatic effects are included by allowing sudden hops between different potential energy surfaces. After averaging over many trajectories, a reasonable description of nonadiabatic dynamics is achieved at low computational cost. Particularly, when we combine the trajectory surface-hopping dynamics with on-the-fly molecular dynamics, it is possible to describe the nonadiabatic dynamics of realistic polyatomic molecular systems at a fully atomic level with all degrees of freedom included. After the introduction of hybrid multiscale methods, the simulation of photochemistry in solutions and in biological systems becomes possible. The simulation results provide important information concerning the nonadiabatic dynamics of realistic polyatomic systems, such as excited-state lifetime, major active molecular motions, reaction channels, and branching ratio. This review article summarizes some progresses in this field. After briefly introducing the basic theory of widely used fewest switches surface-hopping dynamics methods, we mainly focus on several numerical details in the implementation of on-the-fly fewest switches surface-hopping dynamics. The seamless combination of surface-hopping dynamics and electronic-structure calculations is emphasized in this review, rather than an exhaustive discussion of rigorous nonadiabatic dynamics theories. Numerical methods to estimate nonadiabatic coupling terms are discussed, which allow us to perform the trajectory surface-hopping calculations when the nonadiabatic coupling vectors are not available in the electronic structure calculations. Several important issues, such as decoherence corrections, diabatic or adiabatic representations, and initial sampling methods are discussed in detail. We also summarize the theoretical treatment of the nonadiabatic dynamics of some interesting molecular systems, which include the photostability of DNA, photo-isomerization of organic systems, photochemistry of transition metal complexes, and photovoltaics. Finally, we discuss the theoretical challenges of this direct dynamics approach and provide an outlook of this field from our personal perspective.

Cu/ZnO/Al₂O₃ is one of the most widely used catalysts in industrial methanol synthesis. However, the reaction mechanism and the nature of the active sites on the catalyst for this reaction are still under debate. Thus, detailed information is needed to understand the catalytic processes occurring on the surface of this catalyst. H₂ is one of the reaction gases in methanol synthesis. Studies of the activation and dissociation behaviors of H₂ on ZnO surfaces are of great importance in understanding the catalytic mechanism of methanol synthesis. In this work, the activation and dissociation processes of H₂ on a ZnO(10${\rm{\bar 1}}$0) single crystal surface were investigated in-situ using ambient-pressure X-ray photoelectron spectroscopy (APXPS) and scanning tunneling microscopy (STM), two powerful surface characterization techniques. In the APXPS experiments, results indicated the formation of hydroxyl (OH) species on the ZnO single crystal surface at room temperature in 0.3 mbar (1 mbar = 100 Pa) H₂ atmosphere. Meanwhile, STM measurements showed that the ZnO surface was reconstructed from a (1×1) to a (2×1) structure upon introduction of H₂. These observations revealed adsorption behaviors of H₂ the same as those of atomic H on a ZnO(10${\rm{\bar 1}}$0) surface as seen in previous studies, which could be evidence of the dissociative adsorption of H₂ on a ZnO surface. However, H₂O adsorption on ZnO surfaces can also result in the formation of OH species, which can be observed using XPS. The STM results show that the exposure of H₂O also leads to the reconstruction from a (1×1) to a (2×1) structure on the ZnO(10${\rm{\bar 1}}$0) surface upon H₂ introduction. Hence, it is necessary to exclude the influence of H₂O in this work, because there may be trace amounts of H₂O in the H₂ gas. Therefore, we performed a comparative study of H₂ and H₂O on ZnO(10${\rm{\bar 1}}$0) single crystal surface. A downward band bending of 0.3 eV was observed on the ZnO surface in 0.3 mbar H₂ atmosphere using APXPS, while negligible band bending was shown in the case of the H₂O atmosphere. Moreover, thermal stability studies revealed that the OH group formed in the H₂ atmosphere desorbed at a higher temperature than the one resulting from H₂O adsorption, meaning that the two OH groups formed on the ZnO surface were different. Results in this work provide evidence of the dissociative adsorption of H₂ on the ZnO(10${\rm{\bar 1}}$0) surface at room temperature and atmospheric pressure. This is in contrast to previous findings, in which no H₂ dissociation on a ZnO(10${\rm{\bar 1}}$0) surface under ultra-high vacuum conditions was observed, indicating that the activation of H₂ on ZnO surfaces is a pressure dependent process.