2D Ionic Liquid-Like State of Charged Rare-Earth Clusters on a Metal Surface.

Rare-earth complexes are vital for separation chemistry and useful in many advanced applications including emission and energy upconversion. Here, 2D rare-earth clusters having net charges are formed on a metal surface, enabling investigations of their structural and electronic properties on a one-cluster-at-a-time basis using scanning tunneling microscopy. While these ionic complexes are highly mobile on the surface at ≈100 K, their mobility is greatly reduced at 5 K and reveals stable and self-limiting clusters. In each cluster, a pair of charged rare-earth complexes formed by electrostatic and dispersive interactions act as a basic unit, and the clusters are chiral. Unlike other non-ionic molecular clusters formed on the surfaces, these rare-earth clusters show mechanical stability. Moreover, their high mobility on the surface suggests that they are in a 2D liquid-like state.

Accurate Prediction of Adiabatic Ionization Potentials of Organic Molecules using Quantum Chemistry Assisted Machine Learning.

In previous work (Dandu et al., J. Phys. Chem. A, 2022, 126, 4528-4536), we were successful in predicting accurate atomization energies of organic molecules using machine learning (ML) models, obtaining an accuracy as low as 0.1 kcal/mol compared to the G4MP2 method. In this work, we extend the use of these ML models to adiabatic ionization potentials on data sets of energies generated using quantum chemical calculations. Atomic specific corrections that were found to improve atomization energies from quantum chemical calculations have also been used in this study to improve ionization potentials. The quantum chemical calculations were performed on 3405 molecules containing eight or fewer non-hydrogen atoms derived from the QM9 data set, using the B3LYP functional with the 6-31G(2df,p) basis set for optimization. Low-fidelity IPs for these structures were obtained using two density functional methods: B3LYP/6-31+G(2df,p) and ωB97XD/6-311+G(3df,2p). Highly accurate G4MP2 calculations were performed on these optimized structures to obtain high-fidelity IPs to use in ML models based on the low-fidelity IPs. Our best performing ML methods gave IPs of organic molecules within a mean absolute deviation of 0.035 eV from the G4MP2 IPs for the whole data set. This work demonstrates that ML predictions assisted by quantum chemical calculations can be used to successfully predict IPs of organic molecules for use in high throughput screening.

Characterization of just one atom using synchrotron X-rays.

Since the discovery of X-rays by Roentgen in 1895, its use has been ubiquitous, from medical and environmental applications to materials sciences1-5. X-ray characterization requires a large number of atoms and reducing the material quantity is a long-standing goal. Here we show that X-rays can be used to characterize the elemental and chemical state of just one atom. Using a specialized tip as a detector, X-ray-excited currents generated from an iron and a terbium atom coordinated to organic ligands are detected. The fingerprints of a single atom, the L2,3 and M4,5 absorption edge signals for iron and terbium, respectively, are clearly observed in the X-ray absorption spectra. The chemical states of these atoms are characterized by means of near-edge X-ray absorption signals, in which X-ray-excited resonance tunnelling (X-ERT) is dominant for the iron atom. The X-ray signal can be sensed only when the tip is located directly above the atom in extreme proximity, which confirms atomically localized detection in the tunnelling regime. Our work connects synchrotron X-rays with a quantum tunnelling process and opens future X-rays experiments for simultaneous characterizations of elemental and chemical properties of materials at the ultimate single-atom limit.

Microkinetic insights into the role of catalyst and water activity on the nucleation, growth, and dissolution during COF-5 synthesis.

The chemical pathway for synthesizing covalent organic frameworks (COFs) involves a complex medley of reaction sequences over a rippling energy landscape that cannot be adequately described using existing theories. Even with the development of state-of-the-art experimental and computational tools, identifying primary mechanisms of nucleation and growth of COFs remains elusive. Other than empirically, little is known about how the catalyst composition and water activity affect the kinetics of the reaction pathway. Here, for the first time, we employ time-resolved in situ Fourier transform infrared spectroscopy (FT-IR) coupled with a six-parameter microkinetic model consisting of ∼10 million reactions and over 20 000 species. The integrated approach elucidates previously unrecognized roles of catalyst pKa on COF yield and water on growth rate and size distribution. COF crystalline yield increases with decreasing pKa of the catalysts, whereas the effect of water is to reduce the growth rate of COF and broaden the size distribution. The microkinetic model reproduces the experimental data and quantitatively predicts the role of synthesis conditions such as temperature, catalyst, and precursor concentration on the nucleation and growth rates. Furthermore, the model also validates the second-order reaction mechanism of COF-5 and predicts the activation barriers for classical and non-classical growth of COF-5 crystals. The microkinetic model developed here is generalizable to different COFs and other multicomponent systems.

Atomically precise control of rotational dynamics in charged rare-earth complexes on a metal surface.

Complexes containing rare-earth ions attract great attention for their technological applications ranging from spintronic devices to quantum information science. While charged rare-earth coordination complexes are ubiquitous in solution, they are challenging to form on materials surfaces that would allow investigations for potential solid-state applications. Here we report formation and atomically precise manipulation of rare-earth complexes on a gold surface. Although they are composed of multiple units held together by electrostatic interactions, the entire complex rotates as a single unit when electrical energy is supplied from a scanning tunneling microscope tip. Despite the hexagonal symmetry of the gold surface, a counterion at the side of the complex guides precise three-fold rotations and 100% control of their rotational directions is achieved using a negative electric field from the scanning probe tip. This work demonstrates that counterions can be used to control dynamics of rare-earth complexes on materials surfaces for quantum and nanomechanical applications.

Improving the Accuracy of Composite Methods: A G4MP2 Method with G4-like Accuracy and Implications for Machine Learning.

G4MP2 theory has proven to be a reliable and accurate quantum chemical composite method for the calculation of molecular energies using an approximation based on second-order perturbation theory to lower computational costs compared to G4 theory. However, it has been found to have significantly increased errors when applied to larger organic molecules with 10 or more nonhydrogen atoms. We report here on an investigation of the cause of the failure of G4MP2 theory for such larger molecules. One source of error is found to be the "higher-level correction (HLC)", which is meant to correct for deficiencies in correlation contributions to the calculated energies. This is because the HLC assumes that the contribution is independent of the element and the type of bonding involved, both of which become more important with larger molecules. We address this problem by adding an atom-specific correction, dependent on atom type but not bond type, to the higher-level correction. We find that a G4MP2 method that incorporates this modification of the higher-level correction, referred to as G4MP2A, becomes as accurate as G4 theory (for computing enthalpies of formation) for a test set of molecules with less than 10 nonhydrogen atoms as well as a set with 10-14 such atoms, the set of molecules considered here, with a much lower computational cost. The G4MP2A method is also found to significantly improve ionization potentials and electron affinities. Finally, we implemented the G4MP2A energies in a machine learning method to predict molecular energies.

Ground-state actinide chemistry with scalar-relativistic multiconfiguration pair-density functional theory.

We have examined the performance of Multiconfiguration Pair-Density Functional Theory (MC-PDFT) for computing the ground-state properties of actinide species. Specifically, we focused on the properties of UN2 and various actinyl species. The properties obtained with MC-PDFT at the scalar-relativistic level are compared to Kohn-Sham DFT (KS-DFT); complete active space self-consistent field theory, CASSCF; coupled-cluster theory, CCSD(T) and CCSDT; as well as multireference perturbation theory (CASPT2). We examine the degree to which MC-PDFT improves over KS-DFT and CASSCF while aligning with CASPT2, CCSD(T), and CCSDT. All properties that we considered were for the CASPT2 electronic ground states. For structural parameters, MC-PDFT confers very little advantage over KS-DFT, especially the B3LYP density functional. For NpO2 3+, MC-PDFT and local KS-DFT functionals excessively favor the bent structure, whereas CCSDT and CASPT2 predict the bent and linear structures as isoenergetic. For this special case, hybrid KS-DFT functionals like PBE0 and B3LYP provide results closer to CASPT2 and CCSDT than MC-PDFT. On a more positive note, MC-PDFT is very close to CASPT2 and CCSD(T) for the redox potentials, energetics of redox chemical reactions, as well as ligand-binding energies. These are encouraging results since MC-PDFT is more affordable. The best MC-PDFT functional is ft-PBE. Our findings suggest that MC-PDFT can be used to study systems and excited states with larger strong electron correlation effects than were considered here. However, for the systems and properties considered here, KS-DFT functionals do well, justifying their usage as the bulwark of computational actinyl chemistry over the last two to three decades.

Nitrogen Reduction by Multimetallic trans-Uranium Actinide Complexes: A Theoretical Comparison of Np and Pu to U.

There is recent interest in organometallic complexes of the trans-uranium elements. However, preparation and characterization of such complexes are hampered by radioactivity and chemotoxicity issues as well as the air-sensitive and poorly understood behavior of existing compounds. As such, there are no examples of small-molecule activation via redox reactivity of organometallic trans-uranium complexes. This contrasts with the situation for uranium. Indeed, a multimetallic uranium(III) nitride complex was recently synthesized, characterized, and shown to be able to capture and functionalize molecular nitrogen (N2) through a four-electron reduction process, N2 → N24-. The bis-uranium nitride, U-N-U core of this complex is held in a potassium siloxide framework. Importantly, the N24- product could be further functionalized to yield ammonia (NH3) and other desirable species. Using the U-N-U potassium siloxide complex, K3U-N-U, and its cesium analogue, Cs3U-N-U, as starting points, we use scalar-relativistic and spin-orbit coupled density functional theory calculations to shed light on the energetics and mechanism for N2 capture and functionalization. The N2 → N24- reactivity depends on the redox potentials of the U(III) centers and crucially on the stability of the starting complex with respect to decomposition into the mixed oxidation U(IV)/U(III) K2U-N-U or Cs2U-N-U species. For the trans-uranium, Np and Pu analogues of K3U-N-U, the N2 → N24- process is endoergic and would not occur. Interestingly, modification of the Np-O and Pu-O bonds between the actinide cores and the coordinated siloxide framework to Np-NH, Pu-NH, Np-CH2, and Pu-CH2 bonds drastically improves the reaction free energies. The Np-NH species are stable and can reductively capture and reduce N2 to N24-. This is supported by analysis of the spin densities, molecular structure, long-range dispersion effects, as well as spin-orbit coupling effects. These findings chart a path for achieving small-molecule activation with organometallic neptunium analogues of existing uranium complexes.

Performance of density functional theory for describing hetero-metallic active-site motifs for methane-to-methanol conversion in metal-exchanged zeolites.

Methane-to-methanol conversion (MMC) can be facilitated with high methanol selectivities by copper-exchanged zeolites. There are however two open questions regarding the use of these zeolites to facilitate the MMC process. The first concerns the possibility of operating the three cycles in the stepwise MMC process by these zeolites in an isothermal fashion. The second concerns the possibility of improving the methanol yields by systematic substitution of some copper centers in these active sites with other earth-abundant transition metals. Quantum-mechanical computations can be used to compare methane activation by copper oxide species and analogous mixed-metal systems. To carry out such screening, it is important that we use theoretical methods that are accurate and computationally affordable for describing the properties of the hetero-metallic catalytic species. We have examined the performance of 47 exchange-correlation density functionals for predicting the relative spin-state energies and chemical reactivities of six hetero-metallic [M-O-Cu]2+ and [M-O2 -Cu]2+ , (where MCo, Fe, and Ni), species by comparison with coupled cluster theory including iterative single, double excitations as well as perturbative treatment of triple excitations, CCSD(T). We also performed multireference calculations on some of these systems. We considered two types of reactions (hydrogen addition and oxygen addition) that are relevant to MMC. We recommend the use of τ-HCTH and OLYP to determine the spin-state energy splittings in the hetero-metallic motifs. ωB97, ωB97X, ωB97X-D3, and MN15 performed best for predicting the energies of the hydrogen and oxygen addition reactions. In contrast, local, and semilocal functionals do poorly for chemical reactivity. Using [Fe-O-Cu]2+ as a test, we see that the nonlocal functionals perform well for the methane CH activation barrier. In contrast, the semilocal functionals perform rather poorly. © 2018 Wiley Periodicals, Inc.

Thickness-Controlled Quasi-Two-Dimensional Colloidal PbSe Nanoplatelets.

We demonstrate controlled synthesis of discrete two-dimensional (2D) PbSe nanoplatelets (NPLs), with measurable photoluminescence, via oriented attachment directed by quantum dot (QD) surface chemistry. Halide passivation is critical to the growth of these (100) face-dominated NPLs, as corroborated by density functional theory studies. PbCl2 moieties attached to the (111) and (110) of small nanocrystals form interparticle bridges, aligning the QDs and leading to attachment. We find that a 2D bridging network is energetically favored over a 3D network, driving the formation of NPLs. Although PbI2 does not support bridging, its presence destabilizes the large (100) faces of NPLs, providing means for tuning NPL thickness. Spectroscopic analysis confirms the predicted role of thickness-dependent quantum confinement on the NPL band gap.