Orbital-Optimized Versus Time-Dependent Density Functional Calculations of Intramolecular Charge Transfer Excited States.

The performance of time-independent, orbital-optimized calculations of excited states is assessed with respect to charge transfer excitations in organic molecules in comparison to the linear-response time-dependent density functional theory (TD-DFT) approach. A direct optimization method to converge on saddle points of the electronic energy surface is used to carry out calculations with the local density approximation (LDA) and the generalized gradient approximation (GGA) functionals PBE and BLYP for a set of 27 excitations in 15 molecules. The time-independent approach is fully variational and provides a relaxed excited state electron density from which the extent of charge transfer is quantified. The TD-DFT calculations are generally found to provide larger charge transfer distances compared to the orbital-optimized calculations, even when including orbital relaxation effects with the Z-vector method. While the error on the excitation energy relative to theoretical best estimates is found to increase with the extent of charge transfer up to ca. -2 eV for TD-DFT, no correlation is observed for the orbital-optimized approach. The orbital-optimized calculations with the LDA and the GGA functionals provide a mean absolute error of ∼0.7 eV, outperforming TD-DFT with both local and global hybrid functionals for excitations with a long-range charge transfer character. Orbital-optimized calculations with the global hybrid functional B3LYP and the range-separated hybrid functional CAM-B3LYP on a selection of states with short- and long-range charge transfer indicate that inclusion of exact exchange has a small effect on the charge transfer distance, while it significantly improves the excitation energy, with the best-performing functional CAM-B3LYP providing an absolute error typically around 0.15 eV.

Cubic aromaticity in ligand-stabilized doped Au superatoms.

The magnetic response of valence electrons in doped gold-based M@Au8L8 q superatoms (M = Pd, Pt, Ag, Au, Cd, Hg, Ir, and Rh; L = PPh3; and q = 0, +1, +2) is studied by calculating the gauge including magnetically induced currents (GIMIC) in the framework of the auxiliary density functional theory. The studied systems include 24 different combinations of the dopant, total cluster charge, and cluster structure (cubic-like or oblate). The magnetically induced currents (both diatropic and paratropic) are shown to be sensitive to the atomic structure of clusters, the number of superatomic electrons, and the chemical nature of the dopant metal. Among the cubic-like structures, the strongest aromaticity is observed in Pd- and Pt-doped M@Au8L8 0 clusters. Interestingly, Pd- and Pt-doping increases the aromaticity as compared to a similar all-gold eight-electron system Au9L8 +1. With the recent implementation of the GIMIC in the deMon2k code, we investigated the aromaticity in the cubic and butterfly-like M@Au8 core structures, doped with a single M atom from periods 5 and 6 of groups IX-XII. Surprisingly, the doping with Pd and Pt in the cubic structure increases the aromaticity compared to the pure Au case not only near the central atom but encompassing the whole metallic core, following the aromatic trend Pd > Pt > Au. These doped (Pd, Pt)@Au8 nanoclusters show a closed shell 1S21P6 superatom electronic structure corresponding to the cubic aromaticity rule 6n + 2.

Magnetically induced currents and aromaticity in ligand-stabilized Au and AuPt superatoms.

Understanding magnetically induced currents (MICs) in aromatic or metallic nanostructures is crucial for interpreting local magnetic shielding and NMR data. Direct measurements of the induced currents have been successful only in a few planar molecules but their indirect effects are seen in NMR shifts of probe nuclei. Here, we have implemented a numerically efficient method to calculate gauge-including MICs in the formalism of auxiliary density functional theory. We analyze the currents in two experimentally synthesized gold-based, hydrogen-containing ligand-stabilized nanoclusters [HAu9(PPh3)8]2+ and [PtHAu8(PPh3)8]+. Both clusters have a similar octet configuration of Au(6s)-derived delocalized "superatomic" electrons. Surprisingly, Pt-doping in gold increases the diatropic response of the superatomic electrons to an external magnetic field and enhances the aromaticity of [PtHAu8(PPh3)8]+. This is manifested by a stronger shielding of the hydrogen proton in the metal core of the cluster as compared to [HAu9(PPh3)8]2+, causing a significant upfield shift in agreement with experimental proton NMR data measured for these two clusters. Our method allows the determination of local magnetic shielding properties for any component in large 3D nanostructures, opening the door for detailed interpretation of complex NMR spectra.

Solvent-mediated assembly of atom-precise gold-silver nanoclusters to semiconducting one-dimensional materials.

Bottom-up design of functional device components based on nanometer-sized building blocks relies on accurate control of their self-assembly behavior. Atom-precise metal nanoclusters are well-characterizable building blocks for designing tunable nanomaterials, but it has been challenging to achieve directed assembly to macroscopic functional cluster-based materials with highly anisotropic properties. Here, we discover a solvent-mediated assembly of 34-atom intermetallic gold-silver clusters protected by 20 1-ethynyladamantanes into 1D polymers with Ag-Au-Ag bonds between neighboring clusters as shown directly by the atomic structure from single-crystal X-ray diffraction analysis. Density functional theory calculations predict that the single crystals of cluster polymers have a band gap of about 1.3 eV. Field-effect transistors fabricated with single crystals of cluster polymers feature highly anisotropic p-type semiconductor properties with ≈1800-fold conductivity in the direction of the polymer as compared to cross directions, hole mobility of ≈0.02 cm2 V-1 s-1, and an ON/OFF ratio up to ≈4000. This performance holds promise for further design of functional cluster-based materials with highly anisotropic semiconducting properties.

Solubility-Driven Isolation of a Metastable Nonagold Cluster with Body-Centered Cubic Structure.

The conventional synthetic methodology for atomically precise gold nanoclusters by using reduction in solution offers only the thermodynamically most stable nanoclusters. Herein, a solubility-driven isolation strategy is reported to access a metastable gold cluster. The cluster, with the composition of [Au9 (PPh3 )8 ]+ (1), displays an unusual, nearly perfect body-centered cubic (bcc) structure. As revealed by ESI-MS and UV/Vis measurements, the cluster is metastable in solution and converts to the well-known [Au11 (PPh3 )8 Cl2 ]+ (2) within just 90 min. DFT calculations revealed that although both 1 and 2 are eight-electron superatoms, there is a driving force to convert 1 to 2 as shown by the increased cohesion and larger HOMO-LUMO energy gap of 2. The isolation and crystallization of the metastable gold cluster were achieved in a biphasic reaction system in which reduction of gold precursors and crystallization of 1 took place concurrently. This synthetic protocol represents a successful strategy for investigations of other metastable species in metal nanocluster chemistry.