On the quantum dynamical treatment of surface vibrational modes for reactive scattering of H2 from Cu(111) at 925 K

We construct the effective Hartree potential for H2 on Cu(111) as introduced in our earlier work [Dutta et al., J. Chem. Phys. 154, 104103 (2021), and Dutta et al., J. Chem. Phys. 157, 194112 (2022)] starting from the same gas-metal interaction potential obtained for 0 K. Unlike in that work, we now explicitly account for surface expansion at 925 K and investigate different models to describe the surface vibrational modes: (i) a cluster model yielding harmonic normal modes at 0 K and (ii) slab models resulting in phonons at 0 and 925 K according to the quasi-harmonic approximation-all consistently calculated at the density functional theory level with the same exchange-correlation potential. While performing dynamical calculations for the H2(v = 0, j = 0)-Cu(111) system employing Hartree potential constructed with 925 K phonons and surface temperature, (i) the calculated chemisorption probabilities are the highest compared to the other approaches over the energy domain and (ii) the threshold for the reaction probability is the lowest, in close agreement with the experiment. Although the survival probabilities (v' = 0) depict the expected trend (lower in magnitude), the excitation probabilities (v' = 1) display a higher magnitude since the 925 K phonons and surface temperature are more effective for the excitation process compared to the phonons/normal modes obtained from the other approaches investigated to describe the surface.

Control of proton transport and hydrogenation in double-gated graphene

The basal plane of graphene can function as a selective barrier that is permeable to protons1,2 but impermeable to all ions3,4 and gases5,6, stimulating its use in applications such as membranes1,2,7,8, catalysis9,10 and isotope separation11,12. Protons can chemically adsorb on graphene and hydrogenate it13,14, inducing a conductor-insulator transition that has been explored intensively in graphene electronic devices13-17. However, both processes face energy barriers1,12,18 and various strategies have been proposed to accelerate proton transport, for example by introducing vacancies4,7,8, incorporating catalytic metals1,19 or chemically functionalizing the lattice18,20. But these techniques can compromise other properties, such as ion selectivity21,22 or mechanical stability23. Here we show that independent control of the electric field, E, at around 1 V nm-1, and charge-carrier density, n, at around 1 × 1014 cm-2, in double-gated graphene allows the decoupling of proton transport from lattice hydrogenation and can thereby accelerate proton transport such that it approaches the limiting electrolyte current for our devices. Proton transport and hydrogenation can be driven selectively with precision and robustness, enabling proton-based logic and memory graphene devices that have on-off ratios spanning orders of magnitude. Our results show that field effects can accelerate and decouple electrochemical processes in double-gated 2D crystals and demonstrate the possibility of mapping such processes as a function of E and n, which is a new technique for the study of 2D electrode-electrolyte interfaces.

Spiers Memorial Lecture: New directions in molecular scattering

The field of molecular scattering is reviewed as it pertains to gas-gas as well as gas-surface chemical reaction dynamics. We emphasize the importance of collaboration of experiment and theory, from which new directions of research are being pursued on increasingly complex problems. We review both experimental and theoretical advances that provide the modern toolbox available to molecular-scattering studies. We distinguish between two classes of work. The first involves simple systems and uses experiment to validate theory so that from the validated theory, one may learn far more than could ever be measured in the laboratory. The second class involves problems of great complexity that would be difficult or impossible to understand without a partnership of experiment and theory. Key topics covered in this review include crossed-beams reactive scattering and scattering at extremely low energies, where quantum effects dominate. They also include scattering from surfaces, reactive scattering and kinetics at surfaces, and scattering work done at liquid surfaces. The review closes with thoughts on future promising directions of research.

Optimization of Quantum Nuclei Positions with the Adaptive Nuclear-Electronic Orbital Approach

The use of multicomponent methods has become increasingly popular over the last years. Under this framework, nuclei (commonly protons) are treated quantum mechanically on the same footing as the electronic structure problem. Under the use of atomic-centered orbitals, this can lead to some complications as the ideal location of the nuclear basis centers must be optimized. In this contribution, we propose a straightforward approach to determine the position of such centers within the self-consistent cycle of a multicomponent calculation, making use of individual proton charge centroids. We test the method on model systems including the water dimer, a protonated water tetramer, and a porphine system. Comparing to numerical gradient calculations, the adaptive nuclear-electronic orbital (NEO) procedure is able to converge the basis centers to within a few cents of an Ångström and with less than 0.1 kcal/mol differences in absolute energies. This is achieved in one single calculation and with a small added computational effort of up to 80% compared to a regular NEO- self-consistent field run. An example application for the human transketolase proton wire is also provided.

Ultrafast photoinduced C-H bond formation from two small inorganic molecules

The formation of carbon-hydrogen (C-H) bonds via the reaction of small inorganic molecules is of great significance for understanding the fundamental transition from inorganic to organic matter, and thus the origin of life. Yet, the detailed mechanism of the C-H bond formation, particularly the time scale and molecular-level control of the dynamics, remain elusive. Here, we investigate the light-induced bimolecular reaction starting from a van der Waals molecular dimer composed of two small inorganic molecules, H2 and CO. Employing reaction microscopy driven by a tailored two-color light field, we identify the pathways leading to C-H photobonding thereby producing HCO+ ions, and achieve coherent control over the reaction dynamics. Using a femtosecond pump-probe scheme, we capture the ultrafast formation time, i.e., 198 ± 16 femtoseconds. The real-time visualization and coherent control of the dynamics contribute to a deeper understanding of the most fundamental bimolecular reactions responsible for C-H bond formation, thus contributing to elucidate the emergence of organic components in the universe.

Hydrogen atom scattering at the Al2O3(0001) surface: a combined experimental and theoretical study

Investigating atom-surface interactions is the key to an in-depth understanding of chemical processes at interfaces, which are of central importance in many fields - from heterogeneous catalysis to corrosion. In this work, we present a joint experimental and theoretical effort to gain insights into the atomistic details of hydrogen atom scattering at the α-Al2O3(0001) surface. Surprisingly, this system has been hardly studied to date, although hydrogen atoms as well as α-Al2O3 are omnipresent in catalysis as reactive species and support oxide, respectively. We address this system by performing hydrogen atom beam scattering experiments and molecular dynamics (MD) simulations based on a high-dimensional machine learning potential trained to density functional theory data. Using this combination of methods we are able to probe the properties of the multidimensional potential energy surface governing the scattering process. Specifically, we compare the angular distribution and the kinetic energy loss of the scattered atoms obtained in experiment with a large number of MD trajectories, which, moreover, allow to identify the underlying impact sites at the surface.

Quantum and classical molecular dynamics for H atom scattering from graphene

This work presents systematic comparisons between classical molecular dynamics (cMD) and quantum dynamics (QD) simulations of 15-dimensional and 75-dimensional models in their description of H atom scattering from graphene. We use an experimentally validated full-dimensional neural network potential energy surface of a hydrogen atom interacting with a large cell of graphene containing 24 carbon atoms. For quantum dynamics simulations, we apply Monte Carlo canonical polyadic decomposition to transform the original potential energy surface (PES) into a sum of products form and use the multi-layer multi-configuration time-dependent Hartree method to simulate the quantum scattering of a hydrogen or deuterium atom with an initial kinetic energy of 1.96 or 0.96 eV and an incident angle of 0°, i.e., perpendicular to the graphene surface. The cMD and QD initial conditions have been carefully chosen in order to be as close as possible. Our results show little differences between cMD and QD simulations when the incident energy of the H atom is equal to 1.96 eV. However, a large difference in sticking probability is observed when the incident energy of the H atom is equal to 0.96 eV, indicating the predominance of quantum effects. To the best of our knowledge, our work provides the first benchmark of quantum against classical simulations for a system of this size with a realistic PES. Additionally, new projectors are implemented in the Heidelberg multi-configuration time-dependent Hartree package for the calculation of the atom scattering energy transfer distribution as a function of outgoing angles.

Describing the scattering of keV protons through graphene

Implementing two-dimensional materials in technological solutions requires fast, economic, and non-destructive tools to ensure efficient characterization. In this context, scattering of keV protons through free-standing graphene was proposed as an analytical tool. Here, we critically evaluate the predicted effects using classical simulations including a description of the lattice's thermal motion and the membrane corrugation via statistical averaging. Our study shows that the zero-point motion of the lattice atoms alone leads to considerable broadening of the signal that is not properly described by thermal averaging of the interaction potential. In combination with the non-negligible probability for introducing defects, it limits the prospect of proton scattering at 5 keV as an analytic tool.

Machine learning electronic structure methods based on the one-electron reduced density matrix

The theorems of density functional theory (DFT) establish bijective maps between the local external potential of a many-body system and its electron density, wavefunction and, therefore, one-particle reduced density matrix. Building on this foundation, we show that machine learning models based on the one-electron reduced density matrix can be used to generate surrogate electronic structure methods. We generate surrogates of local and hybrid DFT, Hartree-Fock and full configuration interaction theories for systems ranging from small molecules such as water to more complex compounds like benzene and propanol. The surrogate models use the one-electron reduced density matrix as the central quantity to be learned. From the predicted density matrices, we show that either standard quantum chemistry or a second machine-learning model can be used to compute molecular observables, energies, and atomic forces. The surrogate models can generate essentially anything that a standard electronic structure method can, ranging from band gaps and Kohn-Sham orbitals to energy-conserving ab-initio molecular dynamics simulations and infrared spectra, which account for anharmonicity and thermal effects, without the need to employ computationally expensive algorithms such as self-consistent field theory. The algorithms are packaged in an efficient and easy to use Python code, QMLearn, accessible on popular platforms.

Ring Polymer Molecular Dynamics Approach to Quantum Dissociative Chemisorption Rates

A ring polymer molecular dynamics (RPMD) method is proposed for the calculation of the dissociative chemisorption rate coefficient on surfaces. The RPMD rate theory is capable of handling quantum effects such as the zero-point energy and tunneling in dissociative chemisorption, while it relies on classical trajectories for the simulation. Applications to H2 dissociative chemisorption are demonstrated. For the highly activated process on Ag(111), strong deviations from Arrhenius behavior are found at low temperatures and attributed to tunneling. On Pt(111), where the dissociation has a barrierless pathway, the RPMD rate coefficient is found to agree with the experimentally derived thermal sticking coefficient within a factor of 2 over a large temperature range. Significant quantum effects are also found.

First-principles surface reaction rates by ring polymer molecular dynamics and neural network potential: role of anharmonicity and lattice motion

Elementary gas-surface processes are essential steps in heterogeneous catalysis. A predictive understanding of catalytic mechanisms remains challenging due largely to difficulties in accurately characterizing the kinetics of such steps. Experimentally, thermal rates for elementary surface reactions can now be measured using a novel velocity imaging technique, providing a stringent testing ground for ab initio rate theories. Here, we propose to combine ring polymer molecular dynamics (RPMD) rate theory with state-of-the-art first-principles-determined neural network potential to calculate surface reaction rates. Taking NO desorption from Pd(111) as an example, we show that the harmonic approximation and the neglect of lattice motion in the commonly-used transition state theory overestimates and underestimates the entropy change during the desorption process, respectively, leading to opposite errors in rate coefficient predictions and artificial error cancellations. Including anharmonicity and lattice motion, our results reveal a generally neglected surface entropy change due to significant local structural change during desorption and obtain the right answer for the right reasons. Although quantum effects are found to be less important in this system, the proposed approach establishes a more reliable theoretical benchmark for accurately predicting the kinetics of elementary gas-surface processes.

Innovations in nanosynthesis: emerging techniques for precision, scalability, and spatial control in reactions of organic molecules on solid surfaces

The surface science-based approach to synthesising new organic materials on surfaces has gained considerable attention in recent years, owing to its success in facilitating the formation of novel 0D, 1D and 2D architectures. The primary mechanism used to date has been the catalytic transformation of small organic molecules through substrate-enabled reactions. In this Topical Review, we provide an overview of alternate approaches to controlling molecular reactions on surfaces. These approaches include light, electron and ion-initiated reactions, electrospray ionisation deposition-based techniques, collisions of neutral atoms and molecules, and superhydrogenation. We focus on the opportunities afforded by these alternative approaches, in particular where they may offer advantages in terms of selectivity, spatial control or scalability.

CO Inversion on a NaCl(100) Surface: A Multireference Quantum Embedding Study

We develop a multireference quantum embedding model to investigate a recent experimental observation of the isomerization of vibrationally excited CO molecules on a NaCl(100) surface [Science 2020, 367, 175-178]. To explore this mechanism, we built a reduced potential energy surface of CO interacting with NaCl(100) using a second-order multireference perturbation theory, modeling the adsorbate-surface interaction with our previously developed active space embedding theory (ASET). We considered an isolated CO molecule on NaCl(100) and a high-coverage CO monolayer (1/1), and for both we generated potential energy surfaces parametrized by the CO stretching, adsorption, and inversion coordinates. These surfaces are used to determine stationary points and adsorption energies and to perform a vibrational analysis of the states relevant to the inversion mechanism. We found that for near-equilibrium bond lengths, CO adsorbed in the C-down configuration is lower in energy than in the O-down configuration. Stretching of the C-O bond reverses the energetic order of these configurations, supporting the accepted isomerization mechanism. The vibrational constants obtained from these potential energy surfaces show a small (< 10 cm^-1) blue- and red-shift for the C-down and O-down configurations, respectively, in agreement with experimental assignments and previous theoretical studies. Our vibrational analysis of the monolayer case suggests that the O-down configuration is energetically more stable than the C-down one beyond the 16th vibrational excited state of CO, a value slightly smaller than the one from quasi-classical trajectory simulations (22nd) and consistent with the experiment. Our analysis suggests that CO-CO interactions in the monolayer play an important role in stabilizing highly vibrationally excited states in the O-down configuration and reducing the barrier between the C-down and O-down geometries, therefore playing a crucial role in the inversion mechanism.

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.

Velocity-Selected Rotational State Distributions of Nitric Oxide Scattered off Graphene Revealed by Surface-Velocity Map Imaging

We report velocity-dependent internal energy distributions of nitric oxide molecules, NO, scattered off graphene supported on gold to further explore the dynamics of the collision process between NO radicals and graphene. These experiments were performed by directing a molecular beam of NO onto graphene in a surface-velocity map imaging setup, which allowed us to record internal energy distributions of the NO radicals as a function of their velocity. We do not observe bond formation but (1) major contributions from direct inelastic scattering and (2) a smaller trapping-desorption component where some physisorbed NO molecules have residence times on the order of microseconds. This is in agreement with our classical molecular dynamics simulations which also observe a small proportion of two- and multi-bounce collisions events but likewise a small proportion of NO radicals trapped at the surface for the entire length of the molecular dynamics simulations (a few picoseconds). Despite a collision energy of 0.31 eV, which would be sufficient to populate NO(v = 1), we do not detect vibrationally excited nitric oxide.

Examining the Long-Range Effect in Very Long Graphene Nanoribbons: A First-Principles Study

The role of long-range effect on the modulation of the electronic structure of graphene nanoribbons has been little studied due to the limitations of existing theoretical and computational methods. By splitting a molecule top-down and calculating and jointing the Fock matrix of fragments, we developed a computational method suitable for large-size molecules with random doping and arbitrary geometry. Utilizing this method, we achieved the study of the effects of dopants and curvature on graphene nanoribbons (GNRs). It reveals that both dopants and curvature can change the charge distribution of GNRs, while the influence of dopants is more significant and can extend up to 1-3 nm. The electronic excitation properties of GNRs are also largely modified by the doping state or nonuniform curvature. Our findings provide not only a feasible approach for studying the electronic structure of large-size molecules but also the possibility to improve the properties of graphene-based materials by dopants and local curvature.

Molecular Dynamics Simulations of Nitric Oxide Scattering Off Graphene

We performed classical molecular dynamics simulations to model the scattering process of nitric oxide, NO, off graphene supported on gold. This is motivated by our desire to probe the energy transfer in collisions with graphene. Since many of these collision systems comprising of graphene and small molecules have been shown to scatter non-reactively, classical molecular dynamics appear to describe such systems sufficiently. We directed thousands of trajectories of NO molecules onto graphene along the surface normal, while varying impact position, but also speed, orientation, and rotational excitation of the nitric oxide, and compare the results with experimental data. While experiment and theory do not match quantitatively, we observe agreement that the relative amount of kinetic energy lost during the collision increases with increasing initial kinetic energy of the NO. Furthermore, while at higher collision energies, all NO molecules lose some energy, and the vast majority of NO is scattered back, in contrast at low impact energies, the fraction of those nitric oxide molecules that are trapped at the surface increases, and some NO molecules even gain some kinetic energy during the collision process. The collision energy seems to preferentially go into the collective motion of the carbon atoms in the graphene sheet.

Generation of Sub-nanosecond H Atom Pulses for Scattering from Single-Crystal Epitaxial Graphene

Pulsed molecular beams allow high-density gas samples to be cooled to low internal temperatures and to produce narrow speed distributions. They are particularly useful in combination with pulsed-laser-based detection schemes and have also been used as pump pulses in pump–probe experiments with neutral matter. The mechanical response of pulsed valves and chopper wheels limits the duration of these pulses typically to about 10–100 μs. Bunch compression photolysis has been proposed as a means to produce atomic pulses shorter than 1 ns—an experimental capability that would allow new measurements to be made on chemical systems. This technique employs a spatially chirped femtosecond duration photolysis pulse that produced an ensemble of H atom photoproducts that rebunches into a short pulse downstream. To date, this technique could not produce strong enough beams to allow new experiments to be carried out. In this paper, we report production of pulsed H atom beams consistent with a 700 ps pulse duration and with sufficient intensity to carry out differentially resolved inelastic H scattering experiments from a graphene surface. We observe surprisingly narrow angular distributions for H atoms incident normal to the surface. At low incidence energies quasi-elastic scattering dominates, and at high incidence energy we observe a strongly inelastic scattering channel. These results provide the basis for future experiments where the H atoms synchronously collide with a pulsed-laser-excited surface.

Near-ambient pressure velocity map imaging

We present a new velocity map imaging instrument for studying molecular beam surface scattering in a near-ambient pressure (NAP-VMI) environment. The instrument offers the possibility to study chemical reaction dynamics and kinetics where higher pressures are either desired or unavoidable, adding a new tool to help close the “pressure gap” between surface science and applied catalysis. NAP-VMI conditions are created by two sets of ion optics that guide ions through an aperture and map their velocities. The aperture separates the high pressure ionization region and maintains the necessary vacuum in the detector region. The performance of the NAP-VMI is demonstrated with results from N2O photodissociation and N2 scattering from a Pd(110) surface, which are compared under vacuum and at near-ambient pressure (1 × 10−3 mbar). NAP-VMI has the potential to be applied to, and useful for, a broader range of experiments, including photoelectron spectroscopy and scattering with liquid microjets.

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.

Surface electronic corrugation of a one-dimensional topological metal: Bi(114)

The surface of Bi(114) is a striking example where the reduced dimensionality gives rise to structural rearrangement and new states at the surface. Here, we present a study of the surface structure and electronic corrugation of this quasi one-dimensional topological metal based on helium atom scattering (HAS) measurements. In contrast to low-index metal surfaces, upon scattering from the stepped (114) truncation of Bi, a large proportion of the incident beam is scattered into higher order diffraction channels which in combination with the large surface unit cell makes an analysis challenging. The surface electronic corrugation of Bi(114) is determined, using measurements upon scattering normal to the steps, together with quantum mechanical scattering calculations. Therefore, minimisation routines that vary the shape of the corrugation are employed, in order to minimise the deviation between the calculations and experimental scans. Furthermore, we illustrate that quantum mechanical scattering calculations can be used to determine the orientation of the in- and outgoing beam with respect to the stepped surface structure.

Quantum Tunnelling Driven H2 Formation on Graphene

It is commonly believed that it is unfavorable for adsorbed H atoms on carbonaceous surfaces to form H₂ without the help of incident H atoms. Using ring-polymer instanton theory to describe multidimensional tunnelling effects, combined with ab initio electronic structure calculations, we find that these quantum-mechanical simulations reveal a qualitatively different picture. Recombination of adsorbed H atoms, which was believed to be irrelevant at low temperature due to high barriers, is enabled by deep tunnelling, with reaction rates enhanced by tens of orders of magnitude. Furthermore, we identify a new path for H recombination that proceeds via multidimensional tunnelling but would have been predicted to be unfeasible by a simple one-dimensional description of the reaction. The results suggest that hydrogen molecule formation at low temperatures are rather fast processes that should not be ignored in experimental settings and natural environments with graphene, graphite, and other planar carbon segments.

The binding of atomic hydrogen on graphene from density functional theory and diffusion Monte Carlo calculations

In this work, density functional theory (DFT) and diffusion Monte Carlo (DMC) methods are used to calculate the binding energy of a H atom chemisorbed on the graphene surface. The DMC value of the binding energy is about 16% smaller in magnitude than the Perdew–Burke–Ernzerhof (PBE) result. The inclusion of exact exchange through the use of the Heyd–Scuseria–Ernzerhof functional brings the DFT value of the binding energy closer in line with the DMC result. It is also found that there are significant differences in the charge distributions determined using PBE and DMC approaches.

Second-Order Active-Space Embedding Theory

Quantum embedding schemes are a promising way to extend multireference computations to large molecules with strong correlation effects localized on a small number of atoms. This work introduces a second-order active-space embedding theory [ASET(2)] which improves upon mean-field frozen embedding by treating fragment-environment interactions via an approximate canonical transformation. The canonical transformation employed in ASET(2) is formulated using the driven similarity renormalization group. The ASET(2) scheme is benchmarked on the N═N bond dissociation in pentyldiazene, the S₀ to S₁ excitation in 1-octene, and the interaction energy of the O₂-benzene complex. The ASET(2) explicit treatment of fragment-environment interactions beyond the mean-field level generally improves the accuracy of embedded computations, and it becomes necessary to achieve an accurate description of excitation energies of 1-octene and the singlet-triplet gap of the O₂-benzene complex.

Efficiently Computing Excitations of Complex Systems: Linear-Scaling Time-Dependent Embedded Mean-Field Theory in Implicit Solvent

Quantum embedding schemes have the potential to significantly reduce the computational cost of first-principles calculations while maintaining accuracy, particularly for calculations of electronic excitations in complex systems. In this work, I combine time-dependent embedded mean field theory (TD-EMFT) with linear-scaling density functional theory and implicit solvation models, extending previous work within the ONETEP code. This provides a way to perform multilevel calculations of electronic excitations on very large systems, where long-range environmental effects, both quantum and classical in nature, are important. I demonstrate the power of this method by performing simulations on a variety of systems, including a molecular dimer, a chromophore in solution, and a doped molecular crystal. This work paves the way for high accuracy calculations to be performed on large-scale systems that were previously beyond the reach of quantum embedding schemes.

Focusing and spin polarization of atomic hydrogen beam

We have developed a spin-polarized-hydrogen beam with a hexapole magnet. By combining the beam chopper and pulsed laser ionization detection, the time-of-flight of the hydrogen beam was measured, and the dependence of the beam profile on the velocity was acquired, which was consistent with the beam trajectory simulations. The spin polarization of the beam was analyzed by using the Stern-Gerlach-type magnet in combination with the spatial scan of the detection laser. The spin polarization was about 95% at a focusing condition due to the hexapole magnet. The polarization was, on the other hand, reduced to about 70% for the beam at higher velocities, which is consistent with simulation results.

Multi-configuration electron-nuclear dynamics: An open-shell approach

The multi-configuration electron-nuclear dynamics for open shell systems with a spin-unrestricted formalism is described. The mean fields are evaluated using second-order reduced density matrices for electronic and nuclear degrees of freedom. Applications to light-element diatomics including equilibrium geometries, electronic energies, dipole moments, and absorption spectra are presented. The von Neumann entropies for different spin states of a LiH molecule are compared.

Isotopic separation of helium through graphyne membranes: a ring polymer molecular dynamics study

Microscopic-level understanding of the separation mechanism for two-dimensional (2D) membranes is an active area of research due to potential implications of this class of membranes for various technological processes. Helium (He) purification from the natural resources is of particular interest due to the shortfall in its production. In this work, we applied the ring polymer molecular dynamics (RPMD) method to graphdiyne (Gr2) and graphtriyne (Gr3) 2D membranes having variable pore sizes for the separation of He isotopes, and compare for the first time with rigorous quantum calculations. We found that the transmission rate through Gr3 is many orders of magnitude greater than Gr2. The selectivity of either isotope at low temperatures is a consequence of a delicate balance between the zero-point energy effect and tunneling of ⁴He and ³He. In particular, a remarkable tunneling effect is reported on the Gr2 membrane at 10 K, leading to a much larger permeation of the lighter species as compared to the heavier isotope. RPMD provides an efficient approach for studying the separation of He isotopes, taking into account quantum effects of light nuclei motions at low temperatures, which classical methods fail to capture.

Multibounce and Subsurface Scattering of H Atoms Colliding with a van der Waals Solid

We report the results of inelastic differential scattering experiments and full-dimensional molecular dynamics trajectory simulations for 2.76 eV H atoms colliding at a surface of solid xenon. The interaction potential is based on an effective medium theory (EMT) fit to density functional theory (DFT) energies. The translational energy-loss distributions derived from experiment and theory are in excellent agreement. By analyzing trajectories, we find that only a minority of the scattering results from simple single-bounce dynamics. The majority comes from multibounce collisions including subsurface scattering where the H atoms penetrate below the first layer of Xe atoms and subsequently re-emerge to the gas phase. This behavior leads to observable energy-losses as large as 0.5 eV, much larger than a prediction of the binary collision model (0.082 eV), which is often used to estimate the highest possible energy-loss in direct inelastic surface scattering. The sticking probability computed with the EMT-PES (0.15) is dramatically reduced (5 × 10-6) if we employ a full-dimensional potential energy surface (PES) based on Lennard-Jones (LJ) pairwise interactions. Although the LJ-PES accurately describes the interactions near the H-Xe and Xe-Xe energy minima, it drastically overestimates the effective size of the Xe atom seen by the colliding H atom at incidence energies above about 0.1 eV.

Inelastic Scattering of H Atoms from Surfaces

We have developed an instrument that uses photolysis of hydrogen halides to produce nearly monoenergetic hydrogen atom beams and Rydberg atom tagging to obtain accurate angle-resolved time-of-flight distributions of atoms scattered from surfaces. The surfaces are prepared under strict ultrahigh vacuum conditions. Data from these experiments can provide excellent benchmarks for theory, from which it is possible to obtain an atomic scale understanding of the underlying dynamical processes governing H atom adsorption. In this way, the mechanism of adsorption on metals is revealed, showing a penetration-resurfacing mechanism that relies on electronic excitation of the metal by the H atom to succeed. Contrasting this, when H atoms collide at graphene surfaces, the dynamics of bond formation involving at least four carbon atoms govern adsorption. Future perspectives of H atom scattering from surfaces are also outlined.

Computational approaches to dissociative chemisorption on metals: towards chemical accuracy

We review the state-of-the-art in the theory of dissociative chemisorption (DC) of small gas phase molecules on metal surfaces, which is important to modeling heterogeneous catalysis for practical reasons, and for achieving an understanding of the wealth of experimental information that exists for this topic, for fundamental reasons. We first give a quick overview of the experimental state of the field. Turning to the theory, we address the challenge that barrier heights (E_b, which are not observables) for DC on metals cannot yet be calculated with chemical accuracy, although embedded correlated wave function theory and diffusion Monte-Carlo are moving in this direction. For benchmarking, at present chemically accurate E_b can only be derived from dynamics calculations based on a semi-empirically derived density functional (DF), by computing a sticking curve and demonstrating that it is shifted from the curve measured in a supersonic beam experiment by no more than 1 kcal mol^-1. The approach capable of delivering this accuracy is called the specific reaction parameter (SRP) approach to density functional theory (DFT). SRP-DFT relies on DFT and on dynamics calculations, which are most efficiently performed if a potential energy surface (PES) is available. We therefore present a brief review of the DFs that now exist, also considering their performance on databases for E_b for gas phase reactions and DC on metals, and for adsorption to metals. We also consider expressions for SRP-DFs and briefly discuss other electronic structure methods that have addressed the interaction of molecules with metal surfaces. An overview is presented of dynamical models, which make a distinction as to whether or not, and which dissipative channels are modeled, the dissipative channels being surface phonons and electronically non-adiabatic channels such as electron-hole pair excitation. We also discuss the dynamical methods that have been used, such as the quasi-classical trajectory method and quantum dynamical methods like the time-dependent wave packet method and the reaction path Hamiltonian method. Limits on the accuracy of these methods are discussed for DC of diatomic and polyatomic molecules on metal surfaces, paying particular attention to reduced dimensionality approximations that still have to be invoked in wave packet calculations on polyatomic molecules like CH₄. We also address the accuracy of fitting methods, such as recent machine learning methods (like neural network methods) and the corrugation reducing procedure. In discussing the calculation of observables we emphasize the importance of modeling the properties of the supersonic beams in simulating the sticking probability curves measured in the associated experiments. We show that chemically accurate barrier heights have now been extracted for DC in 11 molecule-metal surface systems, some of which form the most accurate core of the only existing database of E_b for DC reactions on metal surfaces (SBH10). The SRP-DFs (or candidate SRP-DFs) that have been derived show transferability in many cases, i.e., they have been shown also to yield chemically accurate E_b for chemically related systems. This can in principle be exploited in simulating rates of catalyzed reactions on nano-particles containing facets and edges, as SRP-DFs may be transferable among systems in which a molecule dissociates on low index and stepped surfaces of the same metal. In many instances SRP-DFs have allowed important conclusions regarding the mechanisms underlying observed experimental trends. An important recent observation is that SRP-DFT based on semi-local exchange DFs has so far only been successful for systems for which the difference of the metal work function and the molecule's electron affinity exceeds 7 eV. A main challenge to SRP-DFT is to extend its applicability to the other systems, which involve a range of important DC reactions of e.g. O₂, H₂O, NH₃, CO₂, and CH₃OH. Recent calculations employing a PES based on a screened hybrid exchange functional suggest that the road to success may be based on using exchange functionals of this category.

Atom-surface van der Waals potentials of topological insulators and semimetals from scattering measurements

The phenomenology of resonant scattering has been known since the earliest experiments upon scattering of atomic beams from surfaces and is a means of obtaining experimental information about the fundamentals of weak adsorption systems in the van der Waals regime. We provide an overview of the experimental approach based on new experimental data for the He-Sb₂Te₃(111) system, followed by a comparative overview and perspective of recent results for topological semimetal and insulator surfaces. Moreover, we shortly discuss the perspectives of calculating helium-surface interaction potentials from ab initio calculations. Our perspective demonstrates that atom-surface scattering provides direct experimental information about the atom-surface interaction in the weak physisorption regime and can also be used to determine the lifetime and mean free path of the trapped atom. We further discuss the effects of elastic and inelastic scattering on the linewidth and lifetime of the trapped He atom with an outlook on future developments and applications.

Small Nuclear Quantum Effects in Scattering of H and D from Graphene

We study nuclear quantum effects in H/D sticking to graphene, comparing scattering experiments at near-zero coverage with classical, quantized, and transition-state calculations. The experiment shows H/D sticking probabilities that are indistinguishable from one another and markedly smaller than those expected from a consideration of zero-point energy shifts of the chemisorption transition state. Inclusion of dynamical effects and vibrational anharmonicity via ring-polymer molecular dynamics (RPMD) yields results that are in good agreement with the experimental results. RPMD also reveals that nuclear quantum effects, while modest, arise primarily from carbon and not from H/D motion, confirming the importance of a C atom rehybridization mechanism associated with H/D sticking on graphene.

Exploration of Graphene Defect Reactivity toward a Hydrogen Radical Utilizing a Preactivated Circumcoronene Model

A preexisting chemisorbed defect is well-known to increase the reactivity of graphene which is normally chemically inert. Specifically, the presence of chemisorbed hydrogen atoms forming an sp³-hybridized C-H bond is known to increase the reactivity of neighboring carbon atoms toward additional hydrogenation with wide-ranging applications from materials science to astrochemistry. In this work, static DFT and DFT-based direct dynamics simulations are used to characterize the reactivity of a graphene sheet around an existing C-H bond defect. The spin density landscape shows how to guide subsequent H atom additions, always bonding most strongly to the carbon atom with greatest spin density. Molecular dynamics of an impinging H atom under thermal conditions with defect graphene was used to determine the statistics of probable reactions. The most frequent outcome is inelastic scattering (48%) and then Eley-Rideal (ER) abstraction of the chemisorbed H atom as vibrationally hot H₂ (40%), while the least likely, but probably most interesting, result is formation of a novel C-H bond (12%). The C-H bonds always form in the β sublattice. The carbon atom in the para position shows to be most reactive toward the incoming H atom, followed by the ortho carbon, in agreement with the spin density computed in the static calculations. Globally, the graphene energy surface is repulsive, but the defects create local channels into this energy surface through which reactants can move locally through and react with the activated surface without a barrier.

Hydrogen abstraction/addition reactions in soot surface growth

The hydrogen abstraction (HB) and addition reactions (HD) by H radicals are examined on a series of polycyclic aromatic hydrocarbon (PAH) monomers and models of quasi-surfaces using quasi-classical trajectory (QCT) method. QCT results reproduce the rate constants of HB reactions on PAH monomers from density functional theory (DFT) in the range of 1500-2700 K. The PAH size has a minor impact on the rates of HB reactions, especially at temperatures beyond 2100 K. In contrast, HD reactions have a clear size dependence, and a larger PAH yields a higher rate. It was also found that the preferred reaction pathway changes from HB to HD reactions at ∼1900 K. The rates of surface HB and HD reactions exceed those in the gas phase by nearly one factor of magnitude. Further analysis of the detailed trajectory of the QCT method reveals that about 50% of surface reactions can be attributed to the events of surface diffusion, which depends on the local energy transfer in gas-surface interactions. However, this phenomenon is not preferred in PAH monomers, as expected. Our finding here questions the treatment of the surface reactions of soot as the product of the first collision between the gaseous species and particle surface. The surface diffusion-induced reactions should be accounted for in the rates of the surface HB and HD reactions. The rate constants of HB and HD reactions on each reactive site (surface zig-zag, surface free-edge and pocket free-edge sites) were calculated by QCT method, and are recommended for the further development of surface chemistry models in soot formation.

Towards free-standing graphane: atomic hydrogen and deuterium bonding to nano-porous graphene

Graphane is formed by bonding hydrogen (and deuterium) atoms to carbon atoms in the graphene mesh, with modification from the pure planar sp² bonding towards an sp³ configuration. Atomic hydrogen (H) and deuterium (D) bonding with C atoms in fully free-standing nano porous graphene (NPG) is achieved, by exploiting low-energy proton (or deuteron) non-destructive irradiation, with unprecedented minimal introduction of defects, as determined by Raman spectroscopy and by the C 1s core level lineshape analysis. Evidence of the H- (or D-) NPG bond formation is obtained by bringing to light the emergence of a H- (or D-) related sp³-distorted component in the C 1s core level, clear fingerprint of H-C (or D-C) covalent bonding. The H (or D) bonding with the C atoms of free-standing graphene reaches more than 1/4 (or 1/3) at% coverage. This non-destructive H-NPG (or D-NPG) chemisorption is very stable at high temperatures up to about 800 K, as monitored by Raman and x-ray photoelectron spectroscopy, with complete healing and restoring of clean graphene above 920 K. The excellent chemical and temperature stability of H- (and D-) NPG opens the way not only towards the formation of semiconducting graphane on large-scale samples, but also to stable graphene functionalisation enabling futuristic applications in advanced detectors for the β-spectrum analysis.

Following the microscopic pathway to adsorption through chemisorption and physisorption wells

Adsorption involves molecules colliding at the surface of a solid and losing their incidence energy by traversing a dynamical pathway to equilibrium. The interactions responsible for energy loss generally include both chemical bond formation (chemisorption) and nonbonding interactions (physisorption). In this work, we present experiments that revealed a quantitative energy landscape and the microscopic pathways underlying a molecule’s equilibration with a surface in a prototypical system: CO adsorption on Au(111). Although the minimum energy state was physisorbed, initial capture of the gas-phase molecule, dosed with an energetic molecular beam, was into a metastable chemisorption state. Subsequent thermal decay of the chemisorbed state led molecules to the physisorption minimum. We found, through detailed balance, that thermal adsorption into both binding states was important at all temperatures.

Regulation of Electronic Structure of Graphene Nanoribbon by Tuning Long-Range Dopant-Dopant Coupling at Distance of Tens of Nanometers

Long-range dopant-dopant coupling in graphene nanoribbon (GNR) has been under intensive study for a very long time. Using a newly developed dopant central insertion scheme (DCIS), we performed first-principles study on multiple H, O, OH, and FeN₄ dopants in long (up to 1000 nm) GNRs and found that, although potential energy of the dopant decays exponentially as a function of distance to the dopant, GNR's electronic density of states (DOS) exhibits wave-like oscillation modulated by dopants separated at a distance up to 100 nm. Such an oscillation strongly infers the purely quantum mechanical resonance states constrained between double quantum wells. This has been unambiguously confirmed by our DCIS study together with a one-dimensional quantum well model study, leading to a proof-of-principle protocol prescribing on-demand GNR-DOS regulation. All these not only reveal the underlining mechanism and importance of long-range dopant-dopant coupling specifically reported in GNR, but also open a novel highway for rationally optimizing and designing two-dimensional materials.

Embedded Mean-Field Theory for Solution-Phase Transition-Metal Polyolefin Catalysis

Decreasing the wall-clock time of quantum mechanics/molecular mechanics (QM/MM) calculations without sacrificing accuracy is a crucial prerequisite for widespread simulation of solution-phase dynamical processes. In this work, we demonstrate the use of embedded mean-field theory (EMFT) as the QM engine in QM/MM molecular dynamics (MD) simulations to examine polyolefin catalysts in solution. We show that employing EMFT in this mode preserves the accuracy of hybrid-functional DFT in the QM region, while providing up to 20-fold reductions in the cost per SCF cycle, thereby increasing the accessible simulation time-scales. We find that EMFT reproduces DFT-computed binding energies and optimized bond lengths to within chemical accuracy, as well as consistently ranking conformer stability. Furthermore, solution-phase EMFT/MM simulations provide insight into the interaction strength of strongly coordinating and bulky counterions.

A zeroth-order active-space frozen-orbital embedding scheme for multireference calculations

Multireference computations of large-scale chemical systems are typically limited by the computational cost of quantum chemistry methods. In this work, we develop a zeroth-order active space embedding theory [ASET(0)], a simple and automatic approach for embedding any multireference dynamical correlation method based on a frozen-orbital treatment of the environment. ASET(0) is combined with the second-order multireference driven similarity renormalization group and tested on several benchmark problems, including the excitation energy of 1-octene and bond-breaking in ethane and pentyldiazene. Finally, we apply ASET(0) to study the singlet-triplet gap of p-benzyne and 9,10-anthracyne diradicals adsorbed on a NaCl surface. Our results show that despite its simplicity, ASET(0) is a powerful and sufficiently accurate embedding scheme applicable when the coupling between the fragment and the environment is in the weak to medium regime.

An experimentally validated neural-network potential energy surface for H-atom on free-standing graphene in full dimensionality

We present a first principles-quality potential energy surface (PES) describing the inter-atomic forces for hydrogen atoms interacting with free-standing graphene. The PES is a high-dimensional neural network potential that has been parameterized to 75 945 data points computed with density-functional theory employing the PBE-D2 functional. Improving over a previously published PES [Jiang et al., Science, 2019, 364, 379], this neural network exhibits a realistic physisorption well and achieves a 10-fold reduction in the RMS fitting error, which is 0.6 meV per atom. The chemisorption barrier is 172 meV, which is lower than that of the REBO-EMFT PES (260 meV). We used this PES to calculate about 1.5 million classical trajectories with carefully selected initial conditions to allow for direct comparison to results of H- and D-atom scattering experiments performed at incidence translational energy of 1.9 eV and a surface temperature of 300 K. The theoretically predicted scattering angular and energy loss distributions are in good agreement with experiment, despite the fact that the experiments employed graphene grown on Pt(111). Compared to previous calculations, the agreement with experiments is improved. The remaining discrepancies between experiment and theory are likely due to the influence of the Pt substrate only present in the experiment.