Towards bridging the structure gap in heterogeneous catalysis: the impact of defects in dissociative chemisorption of methane on Ir surfaces.

A quantitative understanding of the role played by defect sites in heterogeneous catalysis is of great importance in designing new and more effective catalysts. In this work, we report a detailed dynamic study of a key step in methane steam reforming under experimentally relevant conditions on a new high-dimensional potential energy surface determined from first principles data with the aid of machine learning, with which the interactions of CH4 with both the flat Ir(111) and stepped Ir(332) surfaces are described. In particular, we argue based on our simulations that the experimentally observed "negatively activated" dissociative chemisorption of methane on Ir surfaces could be due to a combined effect of defects and high substrate temperature, which lowers the reaction barrier relative to that on terraces. Furthermore, a model based on dynamic information of trapping and reaction channels is proposed, which allows a quantitative prediction of the initial sticking probability for different defect densities, thus helping to close the so-called structure gap in heterogeneous catalysis.

Unified and transferable description of dynamics of H2 dissociative adsorption on multiple copper surfaces via machine learning.

Dynamics of gas-surface reactions is of fundamental importance to various interfacial problems. Accurate modeling of gas-surface reaction dynamics requires a globally accurate reactive potential energy surface (PES), typically specialized for one molecule-surface system with no transferability even from one to another surface. As a proof of concept, we report a novel machine learned PES for H2 reactive scattering from multiple low-index copper surfaces. Trained with limited data, this PES enables a uniformly and chemically accurate description of dissociative adsorption of H2/D2 on Cu(111)/Cu(100)/Cu(110) and offers quantitative insights to the remarkable surface temperature effect. More impressively, this PES is also transferable to describe the dynamics of H2 dissociation on Cu(211) without learning any data on that stepped surface, which can be further improved when adding only a small amount of points. Our work opens a new avenue for studying the dynamics of the structure or step density-sensitive gas-surface reactions relevant to heterogeneous catalysis.

Combined experimental and theoretical study on the ultraviolet photodissociation dynamics of 1-bromo-2,6-difluorobenzene in 267 nm-234 nm.

Thanks to their specific molecular symmetry, aromatic molecules and their derivatives represent ideal model systems in understanding photo-induced chemistry of small molecules. Herein, ultraviolet photodissociation dynamics of the 1-bromo-2,6-difluorobenzene molecule has been visualized via imaging the recoiling velocity distributions of photofragments. The measured recoiling angular distributions of the Br(2P3/2) product vary significantly with the increasing photon energy, arguing against the simple bond-fission mechanism within the C2v symmetry. Ab initio calculations reveal that in addition to the C-Br bond cleavage, two additional internal molecular coordinates that break the molecular symmetry are likely involved. The Br out-of-plane bending opens a direct dissociation pathway on the S1-1A″ (S1-1ππ*) state, while the asymmetric C-F stretching significantly changes the orientation of the transition dipole moment. The present study sheds new light on the effect of symmetry breaking in the photodissociation dynamics of symmetric aryl halides, highlighting the multi-dimensional feature of excited state potential energy surfaces.

A modified generalized Langevin oscillator model for activated gas-surface reactions.

Surface motion has proven to influence the gas-surface reactions in various ways. An adequate model to describe the complex lattice effects in a relatively simple way is therefore highly desirable. In this work, we have modified the widely used Generalized Langevin Oscillator (GLO) model to incorporate the molecule-surface coupling that plays an important role in activated dissociation of polyatomic molecules on metal surfaces. To this end, taking the well-studied CHD3+Ni(111) system as an example, we add a coupling potential linearly dependent on the surface oscillating coordinate, which becomes essential in predicting the dissociative sticking coefficients for reactive scattering. We further scale the mass of the surface oscillator on the basis of a mechanic coupling parameter, which has significantly improved the description of the molecule-surface energy transfer for nonreactive scattering. This so-called modified GLO (MGLO) model retains the simplicity and advantages of the original GLO, while yields much more accurate dynamics results that are in remarkably good agreement with the benchmark data calculated using ab initio molecular dynamics. We argue that the MGLO model is applicable to these highly activated gas-surface reactions with strong molecule-surface couplings.

Eley Rideal recombination of hydrogen atoms on Cu(111): Quantitative role of electronic excitation in cross sections and product distributions.

Eley-Rideal (ER) reaction between a projectile and an adsorbate at a gas-surface interface represents many interesting dynamic features. Electron-hole pairs (EHPs) could play an important role in this process, given the interaction between the accelerated incident atom/molecule and the surface electrons. Here, we examine the EHP effects in a benchmark ER reaction of an impinging H/D atom and pre-covered D/H atoms on Cu(111) by ab initio molecular dynamics simulations. It is found that energy dissipation to EHPs lowers the reaction cross sections by ∼1/3, mainly due to the suppression of the hot atom mechanism, but has a minor effect on rotational distributions, leading to quantitative agreement with experimental data. More interestingly, electronic excitation is responsible for the isotopic effect in translational energy distributions, but probably not for that in angular distributions. Our results highlight the specific role of electronic excitation in detailed dynamic quantities in the ER process.

QTL mapping of flag leaf-related traits in wheat (Triticum aestivum L.).

KEY MESSAGE: QTL controlling flag leaf length, flag leaf width, flag leaf area and flag leaf angle were mapped in wheat. This study aimed to advance our understanding of the genetic mechanisms underlying morphological traits of the flag leaves of wheat (Triticum aestivum L.). A recombinant inbred line (RIL) population derived from ND3331 and the Tibetan semi-wild wheat Zang1817 was used to identify quantitative trait loci (QTLs) controlling flag leaf length (FLL), flag leaf width (FLW), flag leaf area (FLA), and flag leaf angle (FLANG). Using an available simple sequence repeat genetic linkage map, 23 putative QTLs for FLL, FLW, FLA, and FLANG were detected on chromosomes 1B, 2B, 3A, 3D, 4B, 5A, 6B, 7B, and 7D. Individual QTL explained 4.3-68.52% of the phenotypic variance in different environments. Four QTLs for FLL, two for FLW, four for FLA, and five for FLANG were detected in at least two environments. Positive alleles of 17 QTLs for flag leaf-related traits originated from ND3331 and 6 originated from Zang1817. QTLs with pleiotropic effects or multiple linked QTL were also identified on chromosomes 1B, 4B, and 5A; these are potential target regions for fine-mapping and marker-assisted selection in wheat breeding programs.

Effects of surface motion and electron-hole pair excitations in CO2 dissociation and scattering on Ni(100).

The energy transfer between different channels is an important aspect in chemical reactions at surfaces. We investigate here in detail the energy transfer dynamics in a prototypical system, i.e., reactive and nonreactive scattering of CO2 on Ni(100), which is related to heterogeneous catalytic processes with Ni-based catalysts for CO2 reduction. On the basis of our earlier nine-dimensional potential energy surface for CO2/Ni(100), dynamical calculations have been done using the generalized Langevin oscillator (GLO) model combined with local density friction approximation (LDFA), in which the former accounts for the surface motion and the latter accounts for the low-energy electron-hole pair (EHP) excitation. In spite of its simplicity, it is found that the GLO model yields quite satisfactory results, including the significant energy loss and product energy disposal, trapping, and steering dynamics, all of which agree well with the ab initio molecular dynamics ones where many surface atoms are explicitly involved with high computational cost. However, the GLO model fails to describe the reactivity enhancement due to the lattice motion because it intrinsically does not incorporate the variance of barrier height on the surface atom displacement. On the other hand, in LDFA, the energy transferred to EHPs is found to play a minor role and barely alter the dynamics, except for slightly reducing the dissociation probabilities. In addition, vibrational state-selected dissociative sticking probabilities are calculated and previously observed strong mode specificity is confirmed. Our work suggests that further improvement of the GLO model is needed to consider the lattice-induced barrier lowering.

Ab initio molecular dynamics study of the Eley-Rideal reaction of H + Cl-Au(111) → HCl + Au(111): Impact of energy dissipation to surface phonons and electron-hole pairs.

The reaction between an impinging H atom and a Cl atom adsorbed on Au(111), which is a prototype for the Eley-Rideal mechanism, is investigated using ab initio molecular dynamics at different incidence angles. The reaction yielding gaseous HCl with large internal excitation proceeds via both direct and hot-atom mechanisms. Significant energy exchange with both surface phonons and electron-hole pairs has been observed. However, their impact on the reactivity and final state distributions was found to be limited, thanks to the large exothermicity and small barrier of the reaction.

Dissociative chemisorption of methane on Ni(111) using a chemically accurate fifteen dimensional potential energy surface.

A fifteen-dimensional global potential energy surface for the dissociative chemisorption of methane on the rigid Ni(111) surface is developed by a high fidelity fit of ∼200 000 DFT energy points computed using a specific reaction parameter density functional designed to reproduce experimental data. The permutation symmetry and surface periodicity are rigorously enforced using the permutation invariant polynomial-neural network approach. The fitting accuracy of the potential energy surface is thoroughly investigated by examining both static and dynamical attributes of CHD3 dissociation on the frozen surface. This potential energy surface is expected to be chemically accurate as after correction for surface temperature effects it reproduces the measured initial sticking probabilities of CHD3 on Ni(111) for various incidence conditions.

Ultrafast Excited-State Dynamics of Cytosine Aza-Derivative and Analogues.

Excited state dynamics of 5-azacytosine (5-AC), 2,4-diamino-1,3,5-triazine (2,4-DT), and 2-amino-1,3,5-triazine (2-AT) were comprehensively investigated by steady state absorption, fluorescence, and femtosecond transient absorption measurements. Time-dependent density functional theory (TDDFT) calculations were performed to help assign the absorption bands and understand the excited state decay mechanisms. The experimental results of excited singlet state dynamics for 5-AC, 2,4-DT, and 2-AT with femtosecond time resolution were reported for the first time. Two distinct decay pathways, with ∼1 ps and tens of picosecond lifetimes, were observed in 5-AC. Only one decay pathway with 17 ps lifetime was observed in 2,4-DT while an emissive state was found in 2-AT. TDDFT calculations suggest that 5-AC has a dark nπ* (S1) state below the first allowed ππ* (S2) state, which leads to the ultrafast decay of the ππ* state. In 2,4-DT, there is no dark nπ* state below the bright ππ* (S1) state and the 17 ps lifetime is assigned to the relaxation from the ππ* (S1) state to ground state. Two dark nπ* states (S1 and S2) were found in 2-AT, which exhibits much more complex excited state dynamics compared with the other two. Photoluminescence in 2-AT has been confirmed to be fluorescence emission from its bright ππ* (S3) state. Our results strongly suggest that electronic structures are very sensitive to the substitution on the triazine ring and that the photophysical properties of nucleic acid analogues depend highly on their molecular structures.

Non-Adiabatic Effects on Excited States of Vinylidene Observed with Slow Photoelectron Velocity-Map Imaging.

High-resolution slow photoelectron velocity-map imaging spectra of cryogenically cooled X̃2B2 H2CC- and D2CC- in the region of the vinylidene triplet excited states are reported. Three electronic bands are observed and, with the assistance of electronic structure calculations and quantum dynamics on ab initio-based near-equilibrium potential energy surfaces, are assigned as detachment to the [Formula: see text] 3B2 (T1), b̃ 3A2 (T2), and à 1A2 (S1) excited states of neutral vinylidene. This work provides the first experimental observation of the à singlet excited state of H2CC. While regular vibrational structure is observed for the ã and à electronic bands, a number of irregular features are resolved in the vicinity of the b̃ band vibrational origin. High-level ab initio calculations suggest that this anomalous structure arises from a conical intersection between the ã and b̃ triplet states near the b̃ state minimum, which strongly perturbs the vibrational levels in the two electronic states through nonadiabatic coupling. Using the adiabatic electron affinity of H2CC previously measured to be 0.490(6) eV by Ervin and co-workers [J. Chem. Phys. 1989, 91, 5974], term energies for the excited neutral states of H2CC are found to be T0(ã 3B2) = 2.064(6), T0(b̃ 3A2) = 2.738(6), and T0(à 1A2) = 2.991(6) eV.

Probing Water Dissociation and Oxygen Replacement on Partially Oxygen-Covered Cu(111) by Reflection Absorption Infrared Spectroscopy.

The presence of chemisorbed oxygen on the Cu(111) surface is known to strongly reduce the activation barrier for water dissociation as compared to bare Cu(111). Here, we present direct experimental evidence for the hydrogen abstraction mechanism responsible for the facile H2O dissociation on an O/Cu(111) surface using reflection absorption infrared spectroscopy (RAIRS) in combination with isotopically labeled reactants. We also observe that chemisorbed hydroxyl species produced by water dissociation on the O/Cu(111) surface undergo an efficient hydrogen atom transfer from trapped water molecules, leading to the rapid replacement of the initial oxygen isotope coverage and the detection of only a single hydroxyl isotopologue on the surface, in apparent contradiction with the hydrogen abstraction mechanism. In the presence of Cu2O oxide islands on the O/Cu(111) surface, water dissociation occurs selectively at the edges of those islands, leading to the self-assembly of isotopically ordered structures.

Hydrogen atom collisions with a semiconductor efficiently promote electrons to the conduction band.

The Born-Oppenheimer approximation is the keystone of modern computational chemistry and there is wide interest in understanding under what conditions it remains valid. Hydrogen atom scattering from insulator, semi-metal and metal surfaces has helped provide such information. The approximation is adequate for insulators and for metals it fails, but not severely. Here we present hydrogen atom scattering from a semiconductor surface: Ge(111)c(2 × 8). Experiments show bimodal energy-loss distributions revealing two channels. Molecular dynamics trajectories within the Born-Oppenheimer approximation reproduce one channel quantitatively. The second channel transfers much more energy and is absent in simulations. It grows with hydrogen atom incidence energy and exhibits an energy-loss onset equal to the Ge surface bandgap. This leads us to conclude that hydrogen atom collisions at the surface of a semiconductor are capable of promoting electrons from the valence to the conduction band with high efficiency. Our current understanding fails to explain these observations.

First-Principles Insights into Adiabatic and Nonadiabatic Vibrational Energy-Transfer Dynamics during Molecular Scattering from Metal Surfaces: The Importance of Surface Reactivity.

Energy transfer is ubiquitous during molecular collisions and reactions at gas-surface interfaces. Of particular importance is vibrational energy transfer because of its relevance to bond forming and breaking. In this Perspective, we review recent first-principles studies on vibrational energy-transfer dynamics during molecular scattering from metal surfaces at the state-to-state level. Taking several representative systems as examples, we highlight the intrinsic correlation between vibrational energy transfer in nonreactive scattering and surface reactivity and how it operates in both electronically adiabatic and nonadiabatic pathways. Adiabatically, the presence of a dissociation barrier softens a bond in the impinging molecule and increases its couplings with other molecular modes and surface phonons. In the meantime, the stronger interaction between the molecule and the surface also changes the electronic structure at the barrier, resulting in an increase of nonadiabatic effects. We further discuss future prospects toward a more quantitative understanding of this important surface dynamical process.

Ultrastretchable, Self-Healable, and Tissue-Adhesive Hydrogel Dressings Involving Nanoscale Tannic Acid/Ferric Ion Complexes for Combating Bacterial Infection and Promoting Wound Healing.

Preventing bacterial infections and accelerating wound closure are essential in the process of wound healing. Current wound dressings lack enough mechanical properties, self-healing ability, and tissue adhesiveness, and the bacterial killing also relies on the use of antibiotic drugs. Herein, a well-designed hybrid hydrogel dressing is constructed by simple copolymerization of acrylamide (AM), 3-acrylamido phenylboronic acid (AAPBA), chitosan (CS), and the nanoscale tannic acid (TA)/ferric ion (Fe3+) complex (TFe). The resulting hydrogel possesses lots of free catechol, phenylboronic acid, amine, and hydroxyl groups and contains many reversible and dynamic bonds such as multiple hydrogen bonds and boronate ester bonds, thereby showing satisfactory mechanical properties, fast self-healing ability, and desirable tissue-adhesive performance. Benefiting from the high photothermal conversion efficiency of the TFe, the hydrogel exhibits satisfactory antibacterial activity against both Gram-positive and Gram-negative bacteria. Moreover, the embedded TFe also endows the hydrogel with good antioxidant activity, anti-inflammatory property, and cell proliferation to promote tissue regeneration. Remarkably, in vivo animal assays reveal that the hybrid hydrogel effectively eliminates biofilm bacteria in the wound sites and accelerates the healing process of infected wounds. Taken together, the developed versatile hydrogels overcome the shortcomings of traditional wound dressings and are expected to become potential antibacterial dressings for future biomedical applications.

Rapid In Situ Deposition of Iron-Chelated Polydopamine Coating on the Polyacrylamide Hydrogel Dressings for Combined Photothermal and Chemodynamic Therapy of Skin Wound Infection.

Pathogenic bacterial infections of skin wounds have caused a significant threat to clinical treatment and human life safety. Here, we develop a bactericidal hydrogel dressing consisting of a polyacrylamide (PAM) hydrogel framework with in situ surface-deposition of iron-dopped polydopamine (FePDA). The prepared hydrogel dressing (FePDA-PAM) has a compact surface, good tensile strength, and excellent elastic recovery ability. The introduction of Fe3+ ions improve the photothermal therapy (PTT) efficiency of the PDA and endow the hydrogel dressing with chemodynamic therapy (CDT) properties. In vitro experiments show that the antibacterial effect of FePDA-PAM hydrogel on Staphylococcus aureus reach nearly 100% under the combined action of H2O2 and 808 nm near-infrared (NIR) laser, indicating an excellent combined antibacterial property of PTT and CDT. Furthermore, the FePDA-PAM + H2O2 + NIR treatment group in the in vivo antibacterial experiments displays lowest relative wound area and optimal wound healing within 5 days of treatment, thereby indicating the intensive skin wound disinfection. To summarize, the FePDA-PAM hydrogel has simple preparation and good biosafety. It may serve as a potential wound dressing for the combined PTT/CDT dual-mode antibacterial therapy.

Infrared Activities of Adsorbed Species on Metal Surfaces: The Puzzle of Adsorbed Methyl (CH3).

Reflection-absorption infrared spectroscopy (RAIRS) is widely used to identify molecular adsorbates on metals during surface chemical reactions, but the interpretation of RAIRS data can be difficult with experiment alone. Here, we reveal from first-principles calculations the origin of the contrasting RAIRS spectra of methyl adsorbed on Pt(111) and Ni(111). We find that the dynamic dipole associated with the symmetric C-H stretch vibration of CH3 along surface normal is significant on Pt(111) but negligibly small on Ni(111), explaining the strong IR activity in the former and the absence of any RAIRS peaks in the latter. This difference is correlated to different charge transfer patterns between metals and the adsorbate, which are determined by the different preferred adsorption sites of methyl on the two surfaces. This work highlights the need of electronic structure calculations in interpreting RAIRS spectra of adsorbates on metal surfaces.

Bridging the Gap between Direct Dynamics and Globally Accurate Reactive Potential Energy Surfaces Using Neural Networks.

Direct dynamics simulations become increasingly popular in studying reaction dynamics for complex systems where analytical potential energy surfaces (PESs) are unavailable. Yet, the number and/or the propagation time of trajectories are often limited by high computational costs, and numerous energies and forces generated on-the-fly become wasted after simulations. We demonstrate here an example of reusing only a very small portion of existing direct dynamics data to reconstruct a 90-dimensional globally accurate reactive PES describing the interaction of CO2 with a movable Ni(100) surface based on a machine learning approach. In addition to reproducing previous results with much better statistics, we predict scattering probabilities of CO2 at the state-to-state level, which is extremely demanding for direct dynamics. We propose this unified way to investigate gaseous and gas-surface reactions of medium size, initiating with hundreds of preliminary direct dynamics trajectories, followed by low-cost and high-quality simulations on full-dimensional analytical PESs.

Reactive and Nonreactive Scattering of HCl from Au(111): An Ab Initio Molecular Dynamics Study.

The HCl + Au(111) system has recently become a benchmark for highly activated dissociative chemisorption, which presumably is strongly affected by electron-hole pair excitation. Previous dynamics calculations, which were based on density functional theory at the generalized gradient approximation level (GGA-DFT) for the molecule-surface interaction, have all overestimated measured reaction probabilities by at least an order of magnitude. Here, we perform ab initio molecular dynamics (AIMD) and AIMD with electronic friction (AIMDEF) calculations employing a density functional that includes the attractive van der Waals interaction. Our calculations model the simultaneous and possibly synergistic effects of surface temperature, surface atom motion, electron-hole pair excitation, the molecular beam conditions of the experiments, and the van der Waals interaction on the reactivity. We find that reaction probabilities computed with AIMDEF and the SRP32-vdW functional still overestimate the measured reaction probabilities, by a factor 18 for the highest incidence energy at which measurements were performed (≈2.5 eV). Even granting that the experiment could have underestimated the sticking probability by about a factor three, this still translates into a considerable overestimation of the reactivity by the current theory. Likewise, scaled transition probabilities for vibrational excitation from ν = 1, j = 1 to ν = 2 are overestimated by the AIMDEF theory, by factors 3-8 depending on the initial conditions modeled. Energy losses to the surface and translational energy losses are, however, in good agreement with experimental values.

Vibrational control of selective bond cleavage in dissociative chemisorption of methanol on Cu(111).

Controlling product branching ratios in a chemical reaction represents a desired but difficult achievement in chemistry. In this work, we demonstrate the first example of altering the branching ratios in a multichannel reaction, i.e., methanol dissociative chemisorption on Cu(111), via selectively exciting specific vibrational modes. To this end, we develop a globally accurate full-dimensional potential energy surface for the CH3OH/Cu(111) system and perform extensive vibrational state-selected molecular dynamics simulations. Our results show that O-H/C-H/C-O stretching vibrational excitations substantially enhance the respective bond scission processes, representing extraordinary bond selectivity. At a given total energy, the branching ratio of C-O/C-H dissociation can increase by as large as 100 times by exciting the C-O stretching mode which possesses an unprecedentedly strong vibrational efficacy on reactivity. This vibrational control can be realized by the well-designed experiment using a linearly polarized laser.