Peeking into the Femtosecond Hot-Carrier Dynamics Reveals Unexpected Mechanisms in Plasmonic Photocatalysis.

Plasmonic-driven photocatalysis may lead to reaction selectivity that cannot be otherwise achieved. A fundamental role is played by hot carriers, i.e., electrons and holes generated upon plasmonic decay within the metal nanostructure interacting with molecular species. Understanding the elusive microscopic mechanism behind such selectivity is a key step in the rational design of hot-carrier reactions. To accomplish that, we present state-of-the-art multiscale simulations, going beyond density functional theory, of hot-carrier injections for the rate-determining step of a photocatalytic reaction. We focus on carbon dioxide reduction, for which it was experimentally shown that the presence of a rhodium nanocube under illumination leads to the selective production of methane against carbon monoxide. We show that selectivity is due to a (predominantly) direct hole injection from rhodium to the reaction intermediate CHO. Unexpectedly, such an injection does not promote the selective reaction path by favoring proper bond breaking but rather by promoting bonding of the proper molecular fragment to the surface.

Dichroism of plasmonic chiral nanoalloys by rational design.

Time-dependent density functional theory (TDDFT) simulations are conducted on a series of chiral gold/silver alloy nanowires to explore whether silver doping can produce an enhancement of circular dichroism at the plasmon resonance in these systems, and to identify the quantum-mechanical origin of the observed effects. We find a strong plasmonic dichroism when one or two helixes of gold atoms are substituted by silver in a linear chiral nanotube, whose pure gold counterpart does not display any plasmonic dichroism, and we rationalize this finding in terms of "decoupling" the destructive interference of excitations in the pure gold nanotube via alloying. However, further attempts to increase the plasmonic dichroism by considering multi-shell gold nanowires in which one entire shell is doped with silver did not produce the desired effect, but rather a decrease in circular dichroism. We show that this latter result is due to a more severe destructive interference in the dipole excitation contributions, and suggest that further amplification should be possible in principle by properly tuning simultaneously the nanowire structure and chemical ordering.

Unrevealing the Nitrogen Elusive Chirality of 3-Sulfanyl and 3-Sulfinyl N-Tosyl Isoindolinones by ECD Spectra: An Experimental and Theoretical Investigation.

The nitrogen-hybridization/pyramidalization of two solvated N-tosylisoindolinone derivatives having chiral residues in adjacent (I) or adjacent and distal (II) position has been investigated by a theoretical-computational procedure based on Molecular Dynamics simulations and Quantum-Chemical calculations. After validation of our methodology in providing a reliable repertory of conformations by modeling the electronic circular dichroism (EDC) spectra, the electronic features associated with N-pyramidalization were further characterized through Natural Bond Order (NBO) analysis. Comparing against the N-geometry observed in crystal structures as a reference, our findings reveal that the presence of neighbouring chiral centers induces a more pronounced N-pyramidalization in solution than in the solid state, both in I and II. Furthermore, NBO analysis confirms that the N-lactam mostly retains the sp2 character but exhibits slight configurational distortion (ξI=13°; ξII=21°), which significantly influences the chiroptical activities observed in ECD spectra of I and II. This substantiates the N-lactams as configurationally stable chiral centers.

Oxygen-Promoted on-Surface Synthesis of Polyboroxine Molecules.

We present a protocol for the on-surface synthesis of polyboroxine molecules derived from boroxine molecules precursors. This process is promoted by oxygen species present on the Au(111) surface: oxygen atoms facilitate the detachment of naphthalene units of trinaphthyl-boroxine molecules and bridge two unsaturated boroxine centers to form a boroxine-O-boroxine chemical motif. X-ray spectroscopic characterization shows that, as the synthesis process proceeds, it progressively tunes the electronic properties of the interface, thus providing a promising route to control the electron level alignment.

Combined experimental and computational study of the photoabsorption of the monodoped and nondoped nanoclusters Au24Pt(SR)18, Ag24Pt(SR)18, and Ag25(SR)18.

Assessing the accuracy of first-principles computational approaches is instrumental to predict electronic excitations in metal nanoclusters with quantitative confidence. Here we describe a validation study on the optical response of a set of monolayer-protected clusters (MPC). The photoabsorption spectra of Ag25(DMBT)18-, Ag24Pt(DMBT)182- and Au24Pt(SC4H9)18, where DMBT is 2,4-dimethylbenzenethiolate and SC4H9 is n-butylthiolate, have been obtained at low temperature and compared with accurate TDDFT calculations. An excellent match between theory and experiment, with typical deviations of less than 0.1 eV, was obtained, thereby validating the accuracy and reliability of the proposed computational framework. Moreover, an analysis of the TDDFT simulations allowed us to ascribe all relevant spectral features to specific transitions between occupied/virtual orbital pairs. The doping effect of Pt on the optical response of these ultrasmall MPC systems was identified and discussed.

A Resolution of Identity Technique to Speed up TDDFT with Hybrid Functionals: Implementation and Application to the Magic Cluster Series Au8n+4(SC6H5)4n+8 (n = 3-6).

The Resolution of Identity (RI) technique has been employed to speed up the use of hybrid exchange-correlation (xc) functionals at the TDDFT level using the Hybrid Diagonal Approximation. The RI has been implemented within the polTDDFT algorithm (a complex damped polarization method) in the AMS/ADF suite of programs. A speedup factor of 30 has been obtained with respect to a previous numerical implementation, albeit with the same level of accuracy. Comparison of TDDFT simulations with the experimental photoabsorption spectra of the cluster series Au8n+4(SR)4n+8(n = 3-6; R = C6H5) showed the excellent accuracy and efficiency of the method. Results were compared with those obtained via the more simplified and computationally cheaper TDDFT+TB and sTDDFT methods. The present method represents an accurate as well as computationally affordable approach to predict photoabsorption spectra of complex species, realizing an optimal compromise between accuracy and computational efficiency, and is suitable for applications to large metal clusters with sizes up to several hundreds of atoms.

Electronic Structure of Pentagonal Carbon Nanocones: An ab Initio Study.

In this work, we investigate the electronic structure of a particular class of carbon nanocones having a pentagonal tip and C5v symmetry. The ground-state nature of the wave function for these structures can be predicted by the recently proposed generalized Hückel rule that extends the original Hückel rule for annulenes to this class of carbon nanocones. In particular, the structures here considered can be classified as closed-shell or anionic/cationic closed-shells, depending on the geometric characteristics of the cone. The goal of this work is to assess the relationship between the electronic configuration of these carbon nanocones and their ability to gain or lose an electron as well as their adsorption capability. For this, the geometry of these structures in the neutral or ionic forms, as well as systems containing either one lithium or fluorine atom, was optimized at the DFT/B3LYP level. It was found that the electron affinity, ionization potential, and the Li or F adsorption energy present an intimate connection to the ground-state wave function character predicted by the generalized Hückel rule. In fact, a peculiar oscillatory energy behavior was discovered, in which the electron affinity, ionization energy, and adsorption energies oscillate with an increase in the nanocone size. The reasoning behind this is that if the anion is closed-shell, then the neutral nanocone will turn out to be a good electron acceptor, increasing the electron affinity and lithium adsorption energy. On the other hand, in the case of a closed-shell cation, this means that the neutral nanocone will easily lose an electron, leading to a smaller ionization potential and higher fluorine adsorption energy.

Effects of Oxygen Adsorption on the Optical Properties of Ag Nanoparticles.

Plasmonic metal nanoparticles are efficient light harvesters with a myriad of sensing- and energy-related applications. For such applications, the optical properties of nanoparticles of metals such as Cu, Ag, and Au can be tuned by controlling the composition, particle size, and shape, but less is known about the effects of oxidation on the plasmon resonances. In this work, we elucidate the effects of O adsorption on the optical properties of Ag particles by evaluating the thermodynamic properties of O-decorated Ag particles with calculations based on the density functional theory and subsequently computing the photoabsorption spectra with a computationally efficient time-dependent density functional theory approach. We identify stable Ag nanoparticle structures with oxidized edges and a quenching of the plasmonic character of the metal particles upon oxidation and trace back this effect to the sp orbitals (or bands) of Ag particles being involved both in the plasmonic excitation and in the hybridization to form bonds with the adsorbed O atoms. Our work has important implications for the understanding and application of plasmonic metal nanoparticles and plasmon-mediated processes under oxidizing environments.

Time Evolution of Plasmonic Features in Pentagonal Ag Clusters.

In the present work, we apply recently developed real-time descriptors to study the time evolution of plasmonic features of pentagonal Ag clusters. The method is based on the propagation of the time-dependent Schrödinger equation within a singly excited TDDFT ansatz. We use transition contribution maps (TCMs) and induced density to characterize the optical longitudinal and transverse response of such clusters, when interacting with pulses resonant with the low-energy (around 2-3 eV, A1) size-dependent or the high-energy (around 4 eV, E1) size-independent peak. TCMs plots on the analyzed clusters, Ag25+ and Ag43+ show off-diagonal peaks consistent with a plasmonic response when a longitudinal pulse resonant at A1 frequency is applied, and dominant diagonal spots, typical of a molecular transition, when a transverse E1 pulse is employed. Induced densities confirm this behavior, with a dipole-like charge distribution in the first case. The optical features show a time delay with respect to the evolution of the external pulse, consistent with those found in the literature for real-time TDDFT calculations on metal clusters.

A computational approach for modeling electronic circular dichroism of solvated chromophores.

The present study consists in a novel computational protocol to model the UV-circular dichroism spectra of solvated species. It makes use of quantum-chemical calculations on a series of conformations of a flexible chromophore or on a series of chromophore/solvent clusters extracted from molecular dynamic simulations. The protocol is described and applied to the aqueous cationic tripeptide GAG+ and to the aqueous neutral decapeptide (GVGVP)2 . The protocol has proven able to: (i) properly consider the conformational motion of solute in the given environment; (ii) give the actual statistical weight of each conformational state; (iii) provide a reliable quantum mechanical method able to reproduce the spectral features. Temperature effects on conformations and spectral properties are properly taken into account. The role of explicit solvent on the conformational analysis and the spectra calculation is discussed. The comparison of the calculated circular dichroism spectra with experimental ones recorded at different temperatures represents a strict validation test of the method.

Optimization of density fitting auxiliary Slater-type basis functions for time-dependent density functional theory.

A new set of auxiliary basis function suitable to fit the induced electron density is presented. Such set has been optimized in order to furnish accurate absorption spectra using the complex polarizability algorithm of time-dependent density functional theory (TDDFT). An automatic procedure has been set up, able, thanks to the definition of suitable descriptors, to evaluate the resemblance of the auxiliary basis-dependent calculated spectra with respect to a reference. In this way, it has been possible to reduce the size of the basis set maximizing the basis set accuracy. Thanks to the choice to employ a collection of molecules for each element, such basis has proven transferable to molecules outside the collection. The final sets are therefore much more accurate and smaller than the previously optimized ones and have been already included in the database of the last release of the AMS suite of programs. The availability of the present new set will allow to improve drastically the applicability range of the polTDDFT method with higher accuracy and less computational effort.

Accurate Vertical Excitation Energies of BODIPY/Aza-BODIPY Derivatives from Excited-State Mean-Field Calculations.

We report a benchmark study of vertical excitation energies and oscillator strengths for the HOMO → LUMO transitions of 17 boron-dipyrromethene (BODIPY) structures, showing a large variety of ring sizes and substituents. Results obtained at the time-dependent density functional theory (TDDFT) and at the delta-self-consistent-field (ΔSCF) by using 13 different exchange correlation kernels (within LDA, GGA, hybrid, and range-separated approximations) are benchmarked against the experimental excitation energies when available. It is found that the time-independent ΔSCF DFT method, when used in combination with hybrid PBE0 and B3LYP functionals, largely outperforms TDDFT and can be quite competitive, in terms of accuracy, with computationally more costly wave function based methods such as CC2 and CASPT2.

Coupling between Plasmonic and Molecular Excitations: TDDFT Investigation of an Ag-Nanorod/BODIPY-Dye Interaction.

A time-dependent density functional theory (TDDFT) computational approach is employed to study the optical coupling between a plasmonic system (a Ag50 nanorod) and a fluorescent dye (BODIPY). It is found that the BODIPY dye can interact with a plasmonic system in a rather different and selective way according to the mutual orientation of the fragments. Indeed, (i) the plasmon excitation turns out to be sensitive to the presence of the BODIPY transition and (ii) this can lead to amplify or suppress the resonance accordingly to the relative orientation of the corresponding transition dipoles. To understand the coupling mechanism, we analyze the shape of the induced density in real space and the Individual Component Map of the Oscillator Strength (ICM-OS) plots and achieve a simple rationalization and insight on the origin and features of the coupling. The resulting possibility of understanding plasmon/fluorophore interactions by simple qualitative arguments opens the way to a rational design of hybrid (plasmon + dye) systems with the desired optical behavior.

Probing Intermolecular H-Bonding Interactions in Cyanuric Acid Networks: Quenching of the N K-Edge Sigma Resonances.

The electronic characterization of the cyanuric acid both in gas phase and when embedded within an H-bonded scheme forming a monolayer on the Au(111) surface has been performed by means of X-ray Photoelectron Spectroscopy (XPS) and Near Edge X-ray Absorption Fine Structure (NEXAFS) spectroscopy. The experimental spectra at the N, O, and C K-edges have been assigned with the support of DFT calculations, and the combination between theory and experiment has allowed to us investigate the effect of the H-bonding intermolecular interaction on the spectra. In particular, the H-bond formation in the monolayer leads to a quenching of the N 1s NEXAFS resonances associated with transitions to the sigma empty orbitals localized on the N-H portion of the imide group. On the other hand, the π* empty states remain substantially unperturbed. From a computational point of view, it has been shown that the DFT-TP scheme is not able to describe the N 1s NEXAFS spectra of these systems, and the configuration mixing has to be included, through the TDDFT approach in conjunction with the range-separated XC CAM-B3LYP functional, to obtain a correct reproduction of the N 1s core spectra.

Continuum Electronic States: The Tiresia Code.

A multicenter (LCAO) B-spline basis is described in detail, and its capabilities concerning affording convergent solutions for electronic continuum states and wavepacket propagation are presented. It forms the core of the Tiresia code, which implements static-DFT and TDDFT hamiltonians, as well as single channel Dyson-DFT and Dyson-TDDFT descriptions to include correlation in the bound states. Together they afford accurate and computationally efficient descriptions of photoionization properties of complex systems, both in the single photon and strong field environments. A number of examples are provided.

Predictive optical photoabsorption of Ag24Au(DMBT)18 - via efficient TDDFT simulations.

We report a computational study via time-dependent density-functional theory (TDDFT) methods of the photo-absorption spectrum of an atomically precise monolayer-protected cluster (MPC), the Ag24Au(DMBT)18 single negative anion, where DMBT is the 2,4-dimethylbenzenethiolate ligand. The use of efficient simulation algorithms, i.e., the complex polarizability polTDDFT approach and the hybrid-diagonal approximation, allows us to employ a variety of exchange-correlation (xc-) functionals at an affordable computational cost. We are thus able to show, first, how the optical response of this prototypical compound, especially but not exclusively in the absorption threshold (low-energy) region, is sensitive to (1) the choice of the xc-functionals employed in the Kohn-Sham equations and the TDDFT kernel and (2) the choice of the MPC geometry. By comparing simulated spectra with precise experimental photoabsorption data obtained from room temperature down to low temperatures, we then demonstrate how a hybrid xc-functional in both the Kohn-Sham equations and the diagonal TDDFT kernel at the crystallographically determined experimental geometry is able to provide a consistent agreement between simulated and measured spectra across the entire optical region. Single-particle decomposition analysis tools finally allow us to understand the physical reason for the failure of non-hybrid approaches.

Optical Activity of Metal Nanoclusters Deposited on Regular and Doped Oxide Supports from First-Principles Simulations.

We report a computational study and analysis of the optical absorption processes of Ag20 and Au20 clusters deposited on the magnesium oxide (100) facet, both regular and including point defects. Ag20 and Au20 are taken as models of metal nanoparticles and their plasmonic response, MgO as a model of a simple oxide support. We consider oxide defects both on the oxygen anion framework (i.e., a neutral oxygen vacancy) and in the magnesium cation framework (i.e., replacing Mg++ with a transition metal: Cu++ or Co++). We relax the clusters' geometries via Density-Functional Theory (DFT) and calculate the photo-absorption spectra via Time-Dependent DFT (TDDFT) simulations on the relaxed geometries. We find that the substrate/cluster interaction induces a broadening and a red-shift of the excited states of the clusters, phenomena that are enhanced by the presence of an oxygen vacancy and its localized excitations. The presence of a transition-metal dopant does not qualitatively affect the spectral profile. However, when it lies next to an oxygen vacancy for Ag20, it can strongly enhance the component of the cluster excitations perpendicular to the surface, thus favoring charge injection.

Plasmonic Circular Dichroism in Chiral Gold Nanowire Dimers.

We report a computational study at the time-dependent density functional theory (TDDFT) level of the chiro-optical spectra of chiral gold nanowires coupled in dimers. Our goal is to explore whether it is possible to overcome destructive interference in single nanowires that damp chiral response in these systems and to achieve intense plasmonic circular dichroism (CD) through a coupling between the nanostructures. We predict a huge enhancement of circular dichroism at the plasmon resonance when two chiral nanowires are intimately coupled in an achiral relative arrangement. Such an effect is even more pronounced when two chiral nanowires are coupled in a chiral relative arrangement. Individual component maps of rotator strength, partial contributions according to the magnetic dipole component, and induced densities allow us to fully rationalize these findings, thus opening the way to the field of plasmonic CD and its rational design.

The Missing Link: Au191(SPh-tBu)66 Janus Nanoparticle with Molecular and Bulk-Metal-like Properties.

Understanding the evolution of the structure and properties in metals from molecule-like to bulk-like has been a long sought fundamental question in science, since Faraday's 1857 work. We report the discovery of a Janus nanomolecule, Au191(SPh-tBu)66 having both molecular and metallic characteristics, explored crystallographically and optically and modeled theoretically. Au191 has an anisotropic, singly twinned structure with an Au155 core protected by a ligand shell made of 24 monomeric [-S-Au-S-] and 6 dimeric [-S-Au-S-Au-S-] staples. The Au155 core is composed of an 89-atom inner core and 66 surface atoms, arranged as [Au3@Au23@Au63]@Au66 concentric shells of atoms. The inner core has a monotwinned/stacking-faulted face-centered-cubic (fcc) structure. Structural evolution in metal nanoparticles has been known to progress from multiply twinned, icosahedral, structures in smaller molecular sizes to untwinned bulk-like fcc monocrystalline nanostructures in larger nanoparticles. The monotwinned inner core structure of the ligand capped Au191 nanomolecule provides the critical missing link, and bridges the size-evolution gap between the molecular multiple-twinning regime and the bulk-metal-like particles with untwinned fcc structure. The Janus nature of the nanoparticle is demonstrated by its optical and electronic properties, with metal-like electron-phonon relaxation and molecule-like long-lived excited states. First-principles theoretical explorations of the electronic structure uncovered electronic stabilization through the opening of a shell-closing gap at the top of the occupied manifold of the delocalized electronic superatom spectrum of the inner core. The electronic stabilization together with the inner core geometric stability and the optimally stapled ligand-capping anchor and secure the stability of the entire nanomolecule.

An efficient hybrid scheme for time dependent density functional theory.

A hybrid approach able to perform Time Dependent Density Functional Theory (TDDFT) simulations with the same accuracy as that of hybrid exchange-correlation (xc-) functionals but at a fraction of the computational cost is developed, implemented, and validated. The scheme, which we name Hybrid Diagonal Approximation (HDA), consists in employing in the response function a hybrid xc-functional (containing a fraction of the non-local Hartree-Fock exchange) only for the diagonal elements of the omega matrix, while the adiabatic local density approximation is employed for the off-diagonal terms. HDA is especially (but not exclusively) advantageous when using Slater type orbital basis sets and allows one to employ them in a uniquely efficient way, as we demonstrate here by implementing HDA in a local version of the Amsterdam Density Functional code. The new protocol is tested on NH3, C6H6, and the [Au25(SCH3)18]- cluster as prototypical cases ranging from small molecules to ligand-protected metal clusters, finding excellent agreement with respect to both full kernel TDDFT simulations and experimental data. Additionally, a specific comparison test between full kernel and HDA is considered at the Casida level on seven other molecular species, which further confirm the accuracy of the approach for all investigated systems. For the [Au25(SCH3)18]- cluster, a speedup by a factor of seven is obtained with respect to the full kernel. The HDA, therefore, promises to provide a quantitative description of the optical properties of medium-sized systems (nanoclusters) at an affordable cost, thanks to its computational efficiency, especially in combination with a complex polarization algorithm version of TDDFT.