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.