Characterizing the Mechanisms of Ca and Mg Carbonate Ion-Pair Formation with Multi-Level Molecular Dynamics/Quantum Mechanics Simulations

The carbonate minerals of Ca and Mg are abundant throughout the lithosphere and have recently garnered significant research interest as possible long-term carbon sinks in the sequestration of atmospheric carbon dioxide. Nonetheless, an understanding of the atomic-level processes comprising their mineralization remains limited. Here, we characterize and contrast the mechanisms of contact ion-pair formation in aqueous Ca and Mg carbonate systems, which represents the most fundamental step leading to the formation of their mineral solids. Utilizing multilevel embedded correlated wavefunction-based ab initio molecular dynamics/quantum mechanics simulations, we characterize not only the dynamics of these processes but also factors arising from the electronic structure of the involved species, revealing further details of the fundamentally different mechanisms for the interconversion between the contact ion-pairs and solvent-shared ion-pairs of Ca versus Mg carbonate.

Introducing the embedded random phase approximation: H2 dissociative adsorption on Cu(111) as an exemplar

The random phase approximation (RPA) as a means of treating electron correlation recently has been shown to outperform standard density functional theory (DFT) approximations in a variety of cases. However, the computational cost of the RPA is substantially more than DFT, especially when aiming to study extended surfaces. Properly accounting for sufficient surface ensemble size, Brillouin zone sampling, and vacuum separation of periodic images in standard periodic-planewave-based DFT code raises the cost to achieve converged results. Here, we show that sub-system embedding schemes enable use of the RPA for modeling heterogeneous reactions at reduced computational cost. We explore two different embedded RPA (emb-RPA) approaches, periodic emb-RPA and cluster emb-RPA. We use the (experimentally and theoretically) well-studied H2 dissociative adsorption on Cu(111) as our exemplar, and first perform full periodic RPA calculations as a benchmark. The full RPA results match well the semi-empirical barrier fit to experimental observables and others derived from high-level computations, e.g., from recent embedded n-electron valence second order perturbation theory [Zhao et al., J. Chem. Theory Comput. 16(11), 7078-7088 (2020)] and quantum Monte Carlo [Doblhoff-Dier et al., J. Chem. Theory Comput. 13(7), 3208-3219 (2017)] simulations. Among the two emb-RPA approaches tested, the cluster emb-RPA accurately reproduces the energy profile (maximum error of 50 meV along the reaction pathway) while reducing the computational cost by approximately two orders of magnitude. We therefore expect that the embedded cluster approach will enable wider RPA implementation in heterogeneous catalysis.

Minimizing the impacts of the ammonia economy on the nitrogen cycle and climate

Ammonia (NH3) is an attractive low-carbon fuel and hydrogen carrier. However, losses and inefficiencies across the value chain could result in reactive nitrogen emissions (NH3, NOx, and N2O), negatively impacting air quality, the environment, human health, and climate. A relatively robust ammonia economy (30 EJ/y) could perturb the global nitrogen cycle by up to 65 Mt/y with a 5% nitrogen loss rate, equivalent to 50% of the current global perturbation caused by fertilizers. Moreover, the emission rate of nitrous oxide (N2O), a potent greenhouse gas and ozone-depleting molecule, determines whether ammonia combustion has a greenhouse footprint comparable to renewable energy sources or higher than coal (100 to 1,400 gCO2e/kWh). The success of the ammonia economy hence hinges on adopting optimal practices and technologies that minimize reactive nitrogen emissions. We discuss how this constraint should be included in the ongoing broad engineering research to reduce environmental concerns and prevent the lock-in of high-leakage practices.

Probing pH-Dependent Dehydration Dynamics of Mg and Ca Cations in Aqueous Solutions with Multi-Level Quantum Mechanics/Molecular Dynamics Simulations

The dehydration of aqueous calcium and magnesium cations is the most fundamental process controlling their reactivity in chemical and biological phenomena, such as the formation of ionic solids or passing through ion channels. It holds particular relevance in light of recent advancements in the development of carbon capture techniques that rely on mineralization for long-term carbon storage. Specifically, dehydration of Ca2+ and Mg2+ is a key step in proposed carbon capture processes aiming to exploit the relatively high concentration of dissolved carbon dioxide in seawater via the formation of carbonate minerals from solvated Ca2+ and Mg2+ cations for sequestration and storage. Nevertheless, atomic-scale understanding of the dehydration of aqueous Ca2+ and Mg2+ cations remains limited. Here, we utilize rare event sampling via density functional theory molecular dynamics and embedded wavefunction theory calculations to elucidate the dehydration dynamics of aqueous Ca2+ and Mg2+. Emphasis is placed on the investigation of the effect pH has on the stability of the different coordination environments. Our results reveal significant differences in the dehydration dynamics of the two cations and provide insight into how they may be modulated by pH changes.

Solvent Dynamics Are Critical to Understanding Carbon Dioxide Dissolution and Hydration in Water

Simulations of carbon dioxide (CO₂) in water may aid in understanding the impact of its accumulation in aquatic environments and help advance technologies for carbon capture and utilization (via, e.g., mineralization). Quantum mechanical (QM) simulations based on static molecular models with polarizable continuum solvation poorly reproduce the energetics of CO₂ hydration to form carbonic acid in water, independent of the level of QM theory employed. Only with density-functional-theory-based molecular dynamics and rare-event sampling, followed by energy corrections based on embedded correlated wavefunction theory (in conjunction with density functional embedding theory), can a close agreement between theory and experiment be achieved. Such multilevel simulations can serve as benchmarks for simpler, less costly models, giving insight into potential errors of the latter. The strong influence of sampling/averaging over dynamical solvent configurations on the energetics stems from the difference in polarity of both the transition state and product (both polar) versus the reactant (nonpolar). When a solute undergoes a change in polarity during reaction, affecting its interaction with the solvent, careful assessment of the energetic contribution of the solvent response to this change is critical. We show that static models (without structural sampling) that incorporate three explicit water molecules can yield far superior results than models with more explicit water molecules because fewer water molecules yield less configurational artifacts. Static models intelligently incorporating both explicit (molecules directly participating in the reaction) and implicit solvation, along with a proper QM theory, e.g., CCSD(T) for closed-shell systems, can close the accuracy gap between static and dynamic models.

Highly Selective Electrochemical Reduction of CO2 into Methane on Nanotwinned Cu

The electrochemical carbon dioxide reduction reaction (CO₂RR) is a promising route to close the carbon cycle by reducing CO₂ into valuable fuels and chemicals. Electrocatalysts with high selectivity toward a single product are economically desirable yet challenging to achieve. Herein, we demonstrated a highly (111)-oriented Cu foil electrocatalyst with dense twin boundaries (TB) (tw-Cu) that showed a high Faradaic efficiency of 86.1 ± 5.3% toward CH₄ at -1.2 ± 0.02 V vs the reversible hydrogen electrode. Theoretical studies suggested that tw-Cu can significantly lower the reduction barrier for the rate-determining hydrogenation of CO compared to planar Cu(111) under working conditions, which suppressed the competing C-C coupling, leading to the experimentally observed high CH₄ selectivity.

Observation of electron orbital signatures of single atoms within metal-phthalocyanines using atomic force microscopy

Resolving the electronic structure of a single atom within a molecule is of fundamental importance for understanding and predicting chemical and physical properties of functional molecules such as molecular catalysts. However, the observation of the orbital signature of an individual atom is challenging. We report here the direct identification of two adjacent transition-metal atoms, Fe and Co, within phthalocyanine molecules using high-resolution noncontact atomic force microscopy (HR-AFM). HR-AFM imaging reveals that the Co atom is brighter and presents four distinct lobes on the horizontal plane whereas the Fe atom displays a "square" morphology. Pico-force spectroscopy measurements show a larger repulsion force of about 5 pN on the tip exerted by Co in comparison to Fe. Our combined experimental and theoretical results demonstrate that both the distinguishable features in AFM images and the variation in the measured forces arise from Co's higher electron orbital occupation above the molecular plane. The ability to directly observe orbital signatures using HR-AFM should provide a promising approach to characterizing the electronic structure of an individual atom in a molecular species and to understand mechanisms of certain chemical reactions.

DFT exchange: sharing perspectives on the workhorse of quantum chemistry and materials science

In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method developers and practitioners. The format of the paper is that of a roundtable discussion, in which the participants express and exchange views on DFT in the form of 302 individual contributions, formulated as responses to a preset list of 26 questions. Supported by a bibliography of 777 entries, the paper represents a broad snapshot of DFT, anno 2022.

Earth-abundant photocatalyst for H 2 generation from NH 3 with light-emitting diode illumination

Catalysts based on platinum group metals have been a major focus of the chemical industry for decades. We show that plasmonic photocatalysis can transform a thermally unreactive, earth-abundant transition metal into a catalytically active site under illumination. Fe active sites in a Cu-Fe antenna-reactor complex achieve efficiencies very similar to Ru for the photocatalytic decomposition of ammonia under ultrafast pulsed illumination. When illuminated with light-emitting diodes rather than lasers, the photocatalytic efficiencies remain comparable, even when the scale of reaction increases by nearly three orders of magnitude. This result demonstrates the potential for highly efficient, electrically driven production of hydrogen from an ammonia carrier with earth-abundant transition metals.

Electrochemical Hydrogenation of CO on Cu(100): Insights from Accurate Multiconfigurational Wavefunction Methods

Copper (Cu) remains the most efficacious electrocatalyst for electrochemical CO₂ reduction (CO₂R). Its activity and selectivity are highly facet-dependent. We recently examined the commonly proposed rate-limiting CO hydrogenation step on Cu(111) via embedded correlated wavefunction (ECW) theory and demonstrated that only this higher-level theory yields predictions consistent with potential-dependent experimental kinetics. Here, to understand the differing activities of Cu(111) and Cu(100) in catalyzing CO₂R, we explore CO hydrogenation on Cu(100) using ECW theory. We predict that the preferred pathway involves the reduction of adsorbed CO (*CO) to *COH via proton-coupled electron transfer (PCET) at working potentials, although *CHO also may form with a kinetically accessible but higher barrier. In contrast, our earlier work on Cu(111) concluded that *COH and *CHO formation via PCET are equally feasible. This work illustrates one possible origin of the facet dependence of CO₂R mechanisms and products on Cu electrodes and sheds light on how the selectivity of CO₂R electrocatalysts can be controlled by the surface morphology.

Plasmonic Photocatalysis with Chemically and Spatially Specific Antenna-Dual Reactor Complexes

Plasmonic antenna-reactor photocatalysts have been shown to convert light efficiently to chemical energy. Virtually all chemical reactions mediated by such complexes to date, however, have involved relatively simple reactions that require only a single type of reaction site. Here, we investigate a planar Al nanodisk antenna with two chemically distinct and spatially separated active sites in the form of Pd and Fe nanodisks, fabricated in 90° and 180° trimer configurations. The photocatalytic reactions H₂ + D₂ → 2HD and NH₃ + D₂ → NH₂D + HD were both investigated on these nanostructured complexes. While the H₂-D₂ exchange reaction showed an additive behavior for the linear (180°) nanodisk complex, the NH₃ + D₂ reaction shows a clear synergistic effect of the position of the reactor nanodisks relative to the central Al nanodisk antenna. This study shows that light-driven chemical reactions can be performed with both chemical and spatial control of the specific reaction steps, demonstrating precisely designed antennas with multiple reactors for tailored control of chemical reactions of increasing complexity.

Assessing cathode property prediction via exchange-correlation functionals with and without long-range dispersion corrections

We benchmark calculated interlayer spacings, average topotactic voltages, thermodynamic stabilities, and band gaps in layered lithium transition-metal oxides (TMOs) and their de-lithiated counterparts, which are used in lithium-ion batteries as positive electrode materials, against available experimental data. Specifically, we examine the accuracy of properties calculated within density functional theory (DFT) using eight different treatments of electron exchange-correlation: the strongly constrained and appropriately normed (SCAN) and Perdew-Burke-Ernzerhof (PBE) density functionals, Hubbard-U-corrected SCAN and PBE (i.e., SCAN+U and PBE+U), and SCAN(+U) and PBE(+U) with added long-range dispersion (D) interactions (i.e., DFT(+U)+D). van der Waals interactions are included respectively via the revised Vydrov-Van Voorhis (rVV10) for SCAN(+U) and the DFT-D3 for PBE(+U). We find that SCAN-based functionals predict larger voltages due to an underestimation of stability of the MO₂ systems, while also predicting smaller interlayer spacings compared to their PBE-based counterparts. Furthermore, adding dispersion corrections to PBE has a greater effect on voltage predictions and interlayer spacings than with SCAN, indicating that DFT-SCAN - despite being a ground-state theory - fortuitously captures some short and medium-range dispersion interactions better than PBE. While SCAN-based and PBE-based functionals yield qualitatively similar band gap predictions, there is no significant quantitative improvement of SCAN-based functionals over the corresponding PBE-based versions. Finally, we expect SCAN-based functionals to yield more accurate property predictions than the respective PBE-based functionals for most TMOs, given SCAN's stronger theoretical underpinning and better predictions of systematic trends in interlayer spacings, intercalation voltages, and band gaps obtained in this work.

Breaking a dative bond with mechanical forces

Bond breaking and forming are essential components of chemical reactions. Recently, the structure and formation of covalent bonds in single molecules have been studied by non-contact atomic force microscopy (AFM). Here, we report the details of a single dative bond breaking process using non-contact AFM. The dative bond between carbon monoxide and ferrous phthalocyanine was ruptured via mechanical forces applied by atomic force microscope tips; the process was quantitatively measured and characterized both experimentally and via quantum-based simulations. Our results show that the bond can be ruptured either by applying an attractive force of ~150 pN or by a repulsive force of ~220 pN with a significant contribution of shear forces, accompanied by changes of the spin state of the system. Our combined experimental and computational studies provide a deeper understanding of the chemical bond breaking process.

Factors Governing Oxygen Vacancy Formation in Oxide Perovskites

The control of oxygen vacancy (V_O) formation is critical to advancing multiple metal-oxide-perovskite-based technologies. We report the construction of a compact linear model for the neutral V_O formation energy in ABO₃ perovskites that reproduces, with reasonable fidelity, Hubbard-U-corrected density functional theory calculations based on the state-of-the-art, strongly constrained and appropriately normed exchange-correlation functional. We obtain a mean absolute error of 0.45 eV for perovskites stable at 298 K, an accuracy that holds across a large, electronically diverse set of ABO₃ perovskites. Our model considers perovskites containing alkaline-earth metals (Ca, Sr, and Ba) and lanthanides (La and Ce) on the A-site and 3d transition metals (Ti, V, Cr, Mn, Fe, Co, and Ni) on the B-site in six different crystal systems (cubic, tetragonal, orthorhombic, hexagonal, rhombohedral, and monoclinic) common to perovskites. Physically intuitive metrics easily extracted from existing experimental thermochemical data or via inexpensive quantum mechanical calculations, including crystal bond dissociation energies and (solid phase) reduction potentials, are key components of the model. Beyond validation of the model against known experimental trends in materials used in solid oxide fuel cells, the model yields new candidate perovskites not contained in our training data set, such as (Bi,Y)(Fe,Co)O₃, which we predict may have favorable thermochemical water-splitting properties. The confluence of sufficient accuracy, efficiency, and interpretability afforded by our model not only facilitates high-throughput computational screening for any application that requires the precise control of V_O concentrations but also provides a clear picture of the dominant physics governing V_O formation in metal-oxide perovskites.

Metal-to-Ligand Charge-Transfer Spectrum of a Ru-Bipyridine-Sensitized TiO2 Cluster from Embedded Multiconfigurational Excited-State Theory

Understanding optical properties of the dye molecule in dye-sensitized solar cells (DSSCs) from first-principles quantum mechanics can contribute to improving the efficiency of such devices. While density functional theory (DFT) and time-dependent DFT have been pivotal in simulating optoelectronic properties of photoanodes used in DSSCs at the atomic scale, questions remain regarding DFT's adequacy and accuracy to furnish critical information needed to understand the various excited-state processes involved. Here, we simulate the absorption spectra of a dye-sensitized solar cell analogue, comprised of a Ru-bipyridine (Ru-bpy) dye molecule and a small TiO₂ cluster via DFT and via an accurate embedded correlated wavefunction (CW) theory. We generated CW spectra for the adsorbed Ru-bpy dye via a recently introduced capped density functional embedding theory or capped-DFET (to generate the embedding potential that accounts for the interaction of the molecule and the TiO₂ cluster). We then combined capped-DFET with the accurate but expensive multiconfigurational complete active space second-order perturbation theory (CASPT2)-embedded CASPT2. Because the CW theory is conducted on only a portion of the total system in the presence of an embedding potential that describes that portion's interaction with its environment, we efficiently obtain CW-quality predictions that reflect local properties of the entire system. Specifically, for example, with capped-DFET and embedded CW theory, we can simulate accurately a plethora of metal-to-ligand charge-transfer excited properties at a manageable computational cost. Here, we predict detailed electronic spectra within the visible region, featuring the lowest three singlet and triplet excited states, along with predictions of the singlets' lifetimes. We illustrated these results using a Jablonski diagram that show the relative energy position of the singlet and longer-lived triplet excited states and analyzed and proposed relaxation paths for the excited state corresponding to the most intense but short-lived absorption (interconversion, intersystem crossing, fluorescence, and phosphorescence) that may lead to longer-lived excited states necessary for efficient charge separation required to generate current in solar cells.

Hot carrier multiplication in plasmonic photocatalysis

Light-induced hot carriers derived from the surface plasmons of metal nanostructures have been shown to be highly promising agents for photocatalysis. While both nonthermal and thermalized hot carriers can potentially contribute to this process, their specific role in any given chemical reaction has generally not been identified. Here, we report the observation that the H2-D2 exchange reaction photocatalyzed by Cu nanoparticles is driven primarily by thermalized hot carriers. The external quantum yield shows an intriguing S-shaped intensity dependence and exceeds 100% for high light intensities, suggesting that hot carrier multiplication plays a role. A simplified model for the quantum yield of thermalized hot carriers reproduces the observed kinetic features of the reaction, validating our hypothesis of a thermalized hot carrier mechanism. A quantum mechanical study reveals that vibrational excitations of the surface Cu-H bond is the likely activation mechanism, further supporting the effectiveness of low-energy thermalized hot carriers in photocatalyzing this reaction.

Revisiting Understanding of Electrochemical CO2 Reduction on Cu(111): Competing Proton-Coupled Electron Transfer Reaction Mechanisms Revealed by Embedded Correlated Wavefunction Theory

Copper (Cu) electrodes, as the most efficacious of CO₂ reduction reaction (CO₂RR) electrocatalysts, serve as prototypes for determining and validating reaction mechanisms associated with electrochemical CO₂ reduction to hydrocarbons. As in situ electrochemical mechanism determination by experiments is still out of reach, such mechanistic analysis typically is conducted using density functional theory (DFT). The semilocal exchange-correlation (XC) approximations most often used to model such catalysis unfortunately engender a basic error: predicting the wrong adsorption site for CO (a key CO₂RR intermediate) on the most ubiquitous facet of Cu, namely, Cu(111). This longstanding inconsistency casts lingering doubt on previous DFT predictions of the attendant CO₂RR kinetics. Here, we apply embedded correlated wavefunction (ECW) theory, which corrects XC functional error, to study the CO₂RR on Cu(111) via both surface hydride (*H) transfer and proton-coupled electron transfer (PCET). We predict that adsorbed CO (*CO) reduces almost equally to two intermediates, namely, hydroxymethylidyne (*COH) and formyl (*CHO) at -0.9 V vs the RHE. In contrast, semilocal DFT approximations predict a strong preference for *COH. With increasing applied potential, the dominance of *COH (formed via potential-independent surface *H transfer) diminishes, switching to the competitive formation of both *CHO and *COH (both formed via potential-dependent PCET). Our results also demonstrate the importance of including explicitly modeled solvent molecules in predicting electron-transfer barriers and reveal the pitfalls of overreliance on simple surface *H transfer models of reduction reactions.

First-Principles Insights into Plasmon-Induced Catalysis

The size- and shape-controlled enhanced optical response of metal nanoparticles (NPs) is referred to as a localized surface plasmon resonance (LSPR). LSPRs result in amplified surface and interparticle electric fields, which then enhance light absorption of the molecules or other materials coupled to the metallic NPs and/or generate hot carriers within the NPs themselves. When mediated by metallic NPs, photocatalysis can take advantage of this unique optical phenomenon. This review highlights the contributions of quantum mechanical modeling in understanding and guiding current attempts to incorporate plasmonic excitations to improve the kinetics of heterogeneously catalyzed reactions. A range of first-principles quantum mechanics techniques has offered insights, from ground-state density functional theory (DFT) to excited-state theories such as multireference correlated wavefunction methods. Here we discuss the advantages and limitations of these methods in the context of accurately capturing plasmonic effects, with accompanying examples.

Ionic Layering and Overcharging in Electrical Double Layers in a Poisson-Boltzmann Model

Electrical double layers (EDLs) play a significant role in a broad range of physical phenomena related to colloidal stability, diffuse-charge dynamics, electrokinetics, and energy storage applications. Recently, it has been suggested that for large ion sizes or multivalent electrolytes, ions can arrange in a layered structure inside the EDLs. However, the widely used Poisson-Boltzmann models for EDLs are unable to capture the details of ion concentration oscillations and the effect of electrolyte valence on such oscillations. Here, by treating a pair of ions as hard spheres below the distance of closest approach and as point charges otherwise, we are able to predict ionic layering without any additional parameters or boundary conditions while still being compatible with the Poisson-Boltzmann framework. Depending on the combination of ion valence, size, and concentration, our model reveals a structured EDL with spatially oscillating ion concentrations. We report the dependence of critical ion concentration, i.e., the ion concentration above which the oscillations are observed, on the counter-ion valence and the ion size. More importantly, our model displays quantitative agreement with the results of computationally intensive models of the EDL. Finally, we analyze the nonequilibrium problem of EDL charging and demonstrate that ionic layering increases the total charge storage capacity and the charging timescale.

DEP-Dots for 3D cell culture: low-cost, high-repeatability, effective 3D cell culture in multiple gel systems

It is known that cells grown in 3D are more tolerant to drug treatment than those grown in dispersion, but the mechanism for this is still not clear; cells grown in 3D have opportunities to develop inter-cell communication, but are also closely packed which may impede diffusion. In this study we examine methods for dielectrophoresis-based cell aggregation of both suspension and adherent cell lines, and compare the effect of various drugs on cells grown in 3D and 2D. Comparing viability of pharmacological interventions on 3D cell clusters against both suspension cells and adherent cells grown in monolayer, as well as against a unicellular organism with no propensity for intracellular communication, we suggest that 3D aggregates of adherent cells, compared to suspension cells, show a substantially different drug response to cells grown in monolayer, which increases as the IC₅₀ is approached. Further, a mathematical model of the system for each agent demonstrates that changes to drug response are due to inherent changes in the system of adherent cells from the 2D to 3D state. Finally, differences in the electrophysiological membrane properties of the adherent cell type suggest this parameter plays an important role in the differences found in the 3D drug response.

The generality of the GUGA MRCI approach in COLUMBUS for treating complex quantum chemistry

The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calculations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g., investigations of π-conjugated biradicaloid compounds, calculations of multitudes of excited states, development of diabatization procedures, and furnishing the electronic structure information for on-the-fly surface nonadiabatic dynamics. With fully variational uncontracted spin-orbit MRCI, COLUMBUS provides a unique possibility of performing high-level calculations on compounds containing heavy atoms up to lanthanides and actinides. Crucial for carrying out all of these calculations effectively is the availability of an efficient parallel code for the CI step. Configuration spaces of several billion in size now can be treated quite routinely on standard parallel computer clusters. Emerging developments in COLUMBUS, including the all configuration mean energy multiconfiguration self-consistent field method and the graphically contracted function method, promise to allow practically unlimited configuration space dimensions. Spin density based on the GUGA approach, analytic spin-orbit energy gradients, possibilities for local electron correlation MR calculations, development of general interfaces for nonadiabatic dynamics, and MRCI linear vibronic coupling models conclude this overview.

Oxidation State of GaP Photoelectrode Surfaces under Electrochemical Conditions for Photocatalytic CO2 Reduction

Illuminated GaP electrodes selectively reduce CO₂ to CH₃OH in aqueous solution. To understand the photoelectrocatalytic mechanism, knowledge of the GaP surface atomic structure in contact with water under relevant electrochemical conditions is essential. However, there remains a debate about the oxidation state of GaP, i.e., whether oxide species are present at the surface. To address this issue, we use density functional theory to investigate the adsorption of oxide species on GaP(110), a stable and active facet for CO₂ reduction. We predict that GaP(110) indeed could be oxidized at the standard reduction potential for CO₂ to CH₃OH. However, we find that unoxidized GaP(110) is stable under illumination, as it corresponds to a highly reducing condition induced by photoexcited electrons. We conclude that an oxidized GaP electrode is very likely unstable thermodynamically under photoelectrochemical conditions for CO₂ reduction, and therefore, the relevant GaP/water interface for catalysis is indeed the unoxidized one.

Surface-Plasmon-Induced Ammonia Decomposition on Copper: Excited-State Reaction Pathways Revealed by Embedded Correlated Wavefunction Theory

Ammonia is a promising hydrogen storage medium; however, its decomposition via conventional thermal catalysis requires a significant amount of thermal energy input in order to overcome the reaction barriers. Here, we use embedded correlated wavefunction (ECW) theory to quantify reaction pathways and energetics for ammonia decomposition (N-H bond dissociation and N₂ and H₂ associative desorption) on copper (Cu) nanoparticles using a Cu (111) surface model. We predict that surface plasmon excitations will be able to facilitate ammonia decomposition by substantially reducing the effective barriers along excited-state pathways. We estimate the reductions in reaction barriers for breaking the first N-H bond and for recombinative desorption of surface-bound nitrogen and hydrogen atoms to be approximately 1.7, 0.8, and 0.5 eV, respectively. Further, by using the experimental N₂ desorption barrier as a reference, we compare the accuracy of various theoretical methods, including plane-wave Kohn-Sham density functional theory calculations with commonly used exchange-correlation functionals, embedded complete active space second-order perturbation theory, and embedded multiconfiguration pair-density functional theory. This work offers further confirmation that the ECW theoretical framework is the most robust for treating highly correlated local electronic structures of solids.

Plasmonic Photocatalysis of Nitrous Oxide into N2 and O2 Using Aluminum-Iridium Antenna-Reactor Nanoparticles

Photocatalysis with optically active "plasmonic" nanoparticles is a growing field in heterogeneous catalysis, with the potential for substantially increasing efficiencies and selectivities of chemical reactions. Here, the decomposition of nitrous oxide (N₂O), a potent anthropogenic greenhouse gas, on illuminated aluminum-iridium (Al-Ir) antenna-reactor plasmonic photocatalysts is reported. Under resonant illumination conditions, N₂ and O₂ are the only observable decomposition products, avoiding the problematic generation of NO_x species observed using other approaches. Because no appreciable change to the apparent activation energy was observed under illumination, the primary reaction enhancement mechanism for Al-Ir is likely due to photothermal heating rather than plasmon-induced hot-carrier contributions. This light-based approach can induce autocatalysis for rapid N₂O conversion, a process with highly promising potential for applications in N₂O abatement technologies, satellite propulsion, or emergency life-support systems in space stations and submarines.

Response to Comment on “Quantifying hot carrier and thermal contributions in plasmonic photocatalysis”

Sivan et al. claim that the methods used to distinguish thermal from hot carrier effects in our recent report are inaccurate and that our data can be explained by a purely thermal mechanism with a fixed activation energy. This conclusion is invalid, because they substantially misinterpret the emissivity of the photocatalyst and assume a linear intensity–dependent temperature in their model that is unrealistic.

Plasmon damping depends on the chemical nature of the nanoparticle interface

The chemical nature of surface adsorbates affects the localized surface plasmon resonance of metal nanoparticles. However, classical electromagnetic simulations are blind to this effect, whereas experiments are typically plagued by ensemble averaging that also includes size and shape variations. In this work, we are able to isolate the contribution of surface adsorbates to the plasmon resonance by carefully selecting adsorbate isomers, using single-particle spectroscopy to obtain homogeneous linewidths, and comparing experimental results to high-level quantum mechanical calculations based on embedded correlated wavefunction theory. Our approach allows us to indisputably show that nanoparticle plasmons are influenced by the chemical nature of the adsorbates 1,7-dicarbadodecaborane(12)-1-thiol (M1) and 1,7-dicarbadodecaborane(12)-9-thiol (M9). These surface adsorbates induce inside the metal electric dipoles that act as additional scattering centers for plasmon dephasing. In contrast, charge transfer from the plasmon to adsorbates-the most widely suggested mechanism to date-does not play a role here.

Unraveling Oxygen Evolution on Iron-Doped β-Nickel Oxyhydroxide: The Key Role of Highly Active Molecular-like Sites

The active site for electrocatalytic water oxidation on the highly active iron(Fe)-doped β-nickel oxyhydroxide (β-NiOOH) electrocatalyst is hotly debated. Here we characterize the oxygen evolution reaction (OER) activity of an unexplored facet of this material with first-principles quantum mechanics. We show that molecular-like 4-fold-lattice-oxygen-coordinated metal sites on the (1̅21̅1) surface may very well be the key active sites in the electrocatalysis. The predicted OER overpotential (η_OER) for a Fe-centered pathway is reduced by 0.34 V relative to a Ni-centered one, consistent with experiments. We further predict unprecedented, near-quantitative lower bounds for the η_OER, of 0.48 and 0.14 V for pure and Fe-doped β-NiOOH(1̅21̅1), respectively. Our hybrid density functional theory calculations favor a heretofore unpredicted pathway involving an iron(IV)-oxo species, Fe⁴⁺=O. We posit that an iron(IV)-oxo intermediate that stably forms under a low-coordination environment and the favorable discharge of Ni³⁺ to Ni²⁺ are key to β-NiOOH's OER activity.