Synthesis of RhH-doped Au-Ag alloy nanoclusters and dopant evolution.

Doping atomically precise metal nanoclusters (NCs) with heterometals is a powerful method for tuning the physicochemical properties of the original NCs at the atomic level. While the heterometals incorporated into metal NCs are limited to group 10-12 metals with closed d-shells, the doping of open d-shell metals remains largely unexplored. Herein, we report the synthesis of Rh-doped Au-Ag alloy NCs by a metal-exchange reaction of [RhHAg24(SPhMe2)18]2- NCs with an Au-thiolate complex. Combined experimental and theoretical structural studies revealed that the synthesized product is a dianionic [RhHAuxAg24-x(SPhMe2)18]2- NC (x = 8-12), consisting of RhH dopant, Au-rich kernel, and Ag-thiolate staple motifs, with the superatomic 8-electron configuration (1S21P6). Under aerobic conditions, the synthesized NCs underwent kernel evolution to generate a 6-electron [RhAuxAg24-x(SPhMe2)18]1- NC (1S21P4), which was initiated by the desorption of hydride from the kernel. Structural analysis of the [RhHAuxAg24-x(SPhMe2)18]2- NC suggests that the kernel evolution is induced by the change in chemical bonds surrounding the hydride in the Au-rich kernel.

Density-Corrected Density Functional Theory for Open Shells: How to Deal with Spin Contamination.

Density functional theory (DFT) is usually used self-consistently to predict chemical properties, but the use of the Hartree-Fock (HF) density improves energetics in certain, well-characterized cases. Density-corrected (DC) DFT provides the theory behind this, but unrestricted Hartree-Fock (UHF) densities yield poor energetics in cases of strong spin contamination. Here we compare with restricted open-shell HF (ROHF) across 13 different functionals and two DC-DFT methods. For significant spin contamination, ROHF densities outperform UHF densities by as much as a factor of 3, depending on the energy functional, and ROHF-DFT improves over self-consistent DFT for most of the tested functionals. We refine the DC(HF)-DFT algorithm to use ROHF densities in cases of severe spin contamination.

One-Step Passivation of Both Sulfur Vacancies and SiO2 Interface Traps of MoS2 Device.

Transition metal dichalcogenides (TMDs) benefit electrical devices with spin-orbit coupling and valley- and topology-related properties. However, TMD-based devices suffer from traps arising from defect sites inside the channel and the gate oxide interface. Deactivating them requires independent treatments, because the origins are dissimilar. This study introduces a single treatment to passivate defects in a multilayer MoS2 FET. By applying back-gate bias, protons from an H-TFSI droplet are injected into the MoS2, penetrating deeply enough to reach the SiO2 gate oxide. The characterizations employing low-temperature transport and deep-level transient spectroscopy (DLTS) studies reveal that the trap density of S vacancies in MoS2 drops to the lowest detection level. The temperature-dependent mobility plot on the SiO2 substrate resembles that of the h-BN substrate, implying that dangling bonds in SiO2 are passivated. The carrier mobility on the SiO2 substrate is enhanced by approximately 2200% after the injection.

Superatom-in-Superatom Nanoclusters: Synthesis, Structure, and Photoluminescence.

We report a new strategy in which a thiolate-protected Ag25 nanocluster can be doped with open d-shell group 8 (Ru, Os) and 9 (Ir) metals by forming metal hydride (RuH2 , OsH2 , IrH) superatoms with a closed d-shell. Structural analyses using various experimental and theoretical methods revealed that the Ag25 nanoclusters were co-doped with the open d-shell metal and hydride species to produce superatom-in-superatom nanoclusters, establishing a novel superatom doping phenomenon for open d-shell metals. The synthesized superatom-in-superatom nanoclusters exhibited dopant-dependent photoluminescence (PL) properties. Comparative PL lifetime studies of the Ag25 nanoclusters doped with 8-10 group metals revealed that both radiative and nonradiative processes were significantly dependent on the dopant. The former is strongly correlated with the electron affinity of the metal dopant, whereas the latter is governed predominantly by the kernel structure changed upon the doping of the metal hydride(s).

Extending density functional theory with near chemical accuracy beyond pure water.

Density functional simulations of condensed phase water are typically inaccurate, due to the inaccuracies of approximate functionals. A recent breakthrough showed that the SCAN approximation can yield chemical accuracy for pure water in all its phases, but only when its density is corrected. This is a crucial step toward first-principles biosimulations. However, weak dispersion forces are ubiquitous and play a key role in noncovalent interactions among biomolecules, but are not included in the new approach. Moreover, naïve inclusion of dispersion in HF-SCAN ruins its high accuracy for pure water. Here we show that systematic application of the principles of density-corrected DFT yields a functional (HF-r2SCAN-DC4) which recovers and not only improves over HF-SCAN for pure water, but also captures vital noncovalent interactions in biomolecules, making it suitable for simulations of solutions.

Conversion between Metavalent and Covalent Bond in Metastable Superlattices Composed of 2D and 3D Sublayers.

Reversible conversion over multimillion times in bond types between metavalent and covalent bonds becomes one of the most promising bases for universal memory. As the conversions have been found in metastable states, an extended category of crystal structures from stable states via redistribution of vacancies, research on kinetic behavior of the vacancies is highly in demand. However, it remains lacking due to difficulties with experimental analysis. Herein, the direct observation of the evolution of chemical states of vacancies clarifies the behavior by combining analysis on charge density distribution, electrical conductivity, and crystal structures. Site-switching of vacancies of Sb2Te3 gradually occurs with diverged energy barriers owing to their own activation code: the accumulation of vacancies triggers spontaneous gliding along atomic planes to relieve electrostatic repulsion. Studies on the behavior can be further applied to multiphase superlattices composed of Sb2Te3 (2D) and GeTe (3D) sublayers, which represent superior memory performances, but their operating mechanisms were still under debate due to their complexity. The site-switching is favorable (suppressed) when Te-Te bonds are formed as physisorption (chemisorption) over the interface between Sb2Te3 (2D) and GeTe (3D) sublayers driven by configurational entropic gain (electrostatic enthalpic loss). Depending on the type of interfaces between sublayers, phases of the superlattices are classified into metastable and stable states, where the conversion could only be achieved in the metastable state. From this comprehensive understanding on the operating mechanism via kinetic behaviors of vacancies and the metastability, further studies toward vacancy engineering are expected in versatile materials.

Optical Duality of Molybdenum Disulfide: Metal and Semiconductor.

The two different light-matter interactions between visible and infrared light are not switchable because control mechanisms have not been elucidated so far, which restricts the effective spectral range in light-sensing devices. In this study, modulation of the effective spectral range is demonstrated using the metal-insulator transition of MoS2. Nondegenerate MoS2 exhibits a photoconductive effect in detecting visible light. In contrast, degenerate MoS2 responds only to mid-infrared (not visible) light by displaying a photoinduced heating effect via free carrier absorption. Depending on the doping level, the optical behavior of MoS2 simulates the photoconductivity of either the semiconductor or the metal, further indicating that the optical metal-insulator transition is coherent with its electrical counterpart. The electrical switchability of MoS2 enables the development of an unprecedented and novel design optical sensor that can detect both visible and mid-IR (wavelength of 9.6 µm) ranges with a singular optoelectronic device.

Improving Results by Improving Densities: Density-Corrected Density Functional Theory.

Density functional theory (DFT) calculations have become widespread in both chemistry and materials, because they usually provide useful accuracy at much lower computational cost than wavefunction-based methods. All practical DFT calculations require an approximation to the unknown exchange-correlation energy, which is then used self-consistently in the Kohn-Sham scheme to produce an approximate energy from an approximate density. Density-corrected DFT is simply the study of the relative contributions to the total energy error. In the vast majority of DFT calculations, the error due to the approximate density is negligible. But with certain classes of functionals applied to certain classes of problems, the density error is sufficiently large as to contribute to the energy noticeably, and its removal leads to much better results. These problems include reaction barriers, torsional barriers involving π-conjugation, halogen bonds, radicals and anions, most stretched bonds, etc. In all such cases, use of a more accurate density significantly improves performance, and often the simple expedient of using the Hartree-Fock density is enough. This Perspective explains what DC-DFT is, where it is likely to improve results, and how DC-DFT can produce more accurate functionals. We also outline challenges and prospects for the field.

Density-Corrected DFT Explained: Questions and Answers.

HF-DFT, the practice of evaluating approximate density functionals on Hartree-Fock densities, has long been used in testing density functional approximations. Density-corrected DFT (DC-DFT) is a general theoretical framework for identifying failures of density functional approximations by separating errors in a functional from errors in its self-consistent (SC) density. Most modern DFT calculations yield highly accurate densities, but important characteristic classes of calculation have large density-driven errors, including reaction barrier heights, electron affinities, radicals and anions in solution, dissociation of heterodimers, and even some torsional barriers. Here, the HF density (if not spin-contaminated) usually yields more accurate and consistent energies than those of the SC density. We use the term DC(HF)-DFT to indicate DC-DFT using HF densities only in such cases. A recent comprehensive study (J. Chem. Theory Comput. 2021, 17, 1368-1379) of HF-DFT led to many unfavorable conclusions. A reanalysis using DC-DFT shows that DC(HF)-DFT substantially improves DFT results precisely when SC densities are flawed.

Ultrafast Carrier-Lattice Interactions and Interlayer Modulations of Bi2Se3 by X-ray Free-Electron Laser Diffraction.

As a 3D topological insulator, bismuth selenide (Bi2Se3) has potential applications for electrically and optically controllable magnetic and optoelectronic devices. Understanding the coupling with its topological phase requires studying the interactions of carriers with the lattice on time scales down to the subpicosecond regime. Here, we investigate the ultrafast carrier-induced lattice contractions and interlayer modulations in Bi2Se3 thin films by time-resolved diffraction using an X-ray free-electron laser. The lattice contraction depends on the carrier concentration and is followed by an interlayer expansion accompanied by oscillations. Using density functional theory and the Lifshitz model, the initial contraction can be explained by van der Waals force modulation of the confined free carrier layers. Our theoretical calculations suggest that the band inversion, related to a topological phase transition, is modulated by the expansion of the interlayer distance. These results provide insights into the topological phase control by light-induced structural change on ultrafast time scales.

Superatom-in-Superatom [RhH@Ag24 (SPhMe2 )18 ]2- Nanocluster.

Heterometal doping is a powerful method for tuning the physicochemical properties of metal nanoclusters. While the heterometals doped into such nanoclusters predominantly include transition metals with closed d-shells, the doping of open d-shell metals remains largely unexplored. Herein, we report the first synthesis of a [RhHAg24 (SPhMe2 )18 ]2- nanocluster, in which a Rh atom with open d-shells ([Kr]4d8 5s1 ) is incorporated into the Ag24 framework by forming a RhH superatom with closed d-shells ([Kr]4d10 ). Combined experimental and theoretical investigations showed that the Ag24 framework was co-doped with Rh and hydride and that the RhH dopant was a superatomic construct of a Pd atom. Additional studies demonstrated that the [RhHAg24 (SPhMe2 )18 ]2- nanocluster was isoelectronic to the [PdAg24 (SPhMe2 )18 ]2- nanocluster with the superatomic 8-electron configuration (1S2 1P6 ). This study demonstrated for the first time that a superatom could be incorporated into a cluster superatom to generate a stable superatom-in-superatom nanocluster.

KS-pies: Kohn-Sham inversion toolkit.

A Kohn-Sham (KS) inversion determines a KS potential and orbitals corresponding to a given electron density, a procedure that has applications in developing and evaluating functionals used in density functional theory. Despite the utility of KS inversions, application of these methods among the research community is disproportionately small. We implement the KS inversion methods of Zhao-Morrison-Parr and Wu-Yang in a framework that simplifies analysis and conversion of the resulting potential in real-space. Fully documented Python scripts integrate with PySCF, a popular electronic structure prediction software, and Fortran alternatives are provided for computational hot spots.

Explaining and Fixing DFT Failures for Torsional Barriers.

Most torsional barriers are predicted with high accuracies (about 1 kJ/mol) by standard semilocal functionals, but a small subset was found to have much larger errors. We created a database of almost 300 carbon-carbon torsional barriers, including 12 poorly behaved barriers, that stem from the Y═C-X group, where Y is O or S and X is a halide. Functionals with enhanced exchange mixing (about 50%) worked well for all barriers. We found that poor actors have delocalization errors caused by hyperconjugation. These problematic calculations are density-sensitive (i.e., DFT predictions change noticeably with the density), and using HF densities (HF-DFT) fixes these issues. For example, conventional B3LYP performs as accurately as exchange-enhanced functionals if the HF density is used. For long-chain conjugated molecules, HF-DFT can be much better than exchange-enhanced functionals. We suggest that HF-PBE0 has the best overall performance.

Density Sensitivity of Empirical Functionals.

Empirical fitting of parameters in approximate density functionals is common. Such fits conflate errors in the self-consistent density with errors in the energy functional, but density-corrected DFT (DC-DFT) separates these two. We illustrate with catastrophic failures of a toy functional applied to H2+ at varying bond lengths, where the standard fitting procedure misses the exact functional; Grimme's D3 fit to noncovalent interactions, which can be contaminated by large density errors such as in the WATER27 and B30 data sets; and double-hybrids trained on self-consistent densities, which can perform poorly on systems with density-driven errors. In these cases, more accurate results are found at no additional cost by using Hartree-Fock (HF) densities instead of self-consistent densities. For binding energies of small water clusters, errors are greatly reduced. Range-separated hybrids with 100% HF at large distances suffer much less from this effect.

Modulation of the photoelectrochemical behavior of Au nanocluster-TiO2 electrode by doping.

Despite the successful debut of gold nanoclusters (Au NCs) in solar cell applications, Au NCs, compared to dyes and quantum dots, have several drawbacks, such as lower extinction coefficients. Any modulation of the physical properties of NCs can have a significant influence on the delicate control of absorbance, energy levels, and charge separation, which are essential to ensure high power conversion efficiency. To this end, we systematically alter the optoelectronic structure of Au18(SR)14 by Ag doping and explain its influence on solar cell performance. Our in-depth spectroscopic and electrochemical characterization combined with computational study reveals that the performance-dictating factors respond in different manners to the Ag doping level, and we determine that the best compromise is the incorporation of a single Ag atom into an Au NC. This new insight highlights the unique aspect of NCs-susceptibility to atomic level doping-and helps establish a new design principle for efficient NC-based solar cells.

Measuring Density-Driven Errors Using Kohn-Sham Inversion.

Kohn-Sham (KS) inversion, that is, the finding of the exact KS potential for a given density, is difficult in localized basis sets. We study the precision and reliability of several inversion schemes, finding estimates of density-driven errors at a useful level of accuracy. In typical cases of substantial density-driven errors, Hartree-Fock density functional theory (HF-DFT) is almost as accurate as DFT evaluated on CCSD(T) densities. A simple approximation in practical HF-DFT also makes errors much smaller than the density-driven errors being calculated. Two paradigm examples, stretched NaCl and the HO·Cl- radical, illustrate just how accurate HF-DFT is.

Density Functional Analysis: The Theory of Density-Corrected DFT.

Density-corrected density functional theory (DC-DFT) is enjoying substantial success in improving semilocal DFT calculations in a wide variety of chemical problems. This paper provides the formal theoretical framework and assumptions for the analysis of any functional minimization with an approximate functional. We generalize DC-DFT to allow the comparison of any two functionals, not just comparison with the exact functional. We introduce a linear interpolation between any two approximations and use the results to analyze global hybrid density functionals. We define the basins of density space in which this analysis should apply and give quantitative criteria for when DC-DFT should apply. We also discuss the effects of strong correlation on the density-driven error, utilizing the restricted HF Hubbard dimer as an example.

Conformational Heterogeneity in Large Macrocyclic Thiophenes.

For conjugated macrocycles, conformational disorder plays a key role in determining whether the unique form of excitons that are fully delocalized over the cyclic framework (cyclic excitons) is formed by photoexcitation. We have investigated the ring size dependence of conformations and photophysical properties of macrocyclic thiophenes of varying ring sizes (C-5NTNV) by using single-molecule fluorescence spectroscopy. We measured modulation depth, M, values and fluorescence intensities. As the ring size increases, the correlation plots of the two parameters show bimodal distributions, revealing that larger macrocycles exhibit extremely congested linear structures. The size dependence of structural changes in macrocyclic thiophenes have been clearly confirmed by molecular dynamics simulation. The number of torsional defects from simulated structures, in conjunction with survival times from fluorescence intensity trajectories and photon coincidence measurements, demonstrated the existence of multiple acyclic chromophores in the larger macrocycles from the ground state due to complete deformation of circular structures.

Halogen and Chalcogen Binding Dominated by Density-Driven Errors.

Dispersion corrections of various kinds usually improve DFT energetics of weak noncovalent interactions. However, in some cases involving molecules or halides, especially those with σ-hole interactions, the density-driven errors of uncorrected DFT are larger than the dispersion corrections. In these abnormal situations, HF-DFT (using Hartree-Fock densities instead of self-consistent densities) greatly improves bond energies, while dispersion corrections can even worsen the results. On the other hand, pnictogen bonds and the S22 data set are normal and are not improved by this procedure. Such effects should be accounted for when parametrizing dispersion interactions.