Mechanistic Insights into the Origins of Selectivity in a Cu-Catalyzed C-H Amidation Reaction.

The catalytic transformation of C-H to C-N bonds offers rapid access to fine chemicals and high-performance materials, but achieving high selectivity from undirected aminations of unactivated C(sp3)-H bonds remains an outstanding challenge. We report the origins of the reactivity and selectivity of a Cu-catalyzed C-H amidation of simple alkanes. Using a combination of experimental and computational mechanistic studies and energy decomposition techniques, we uncover a switch in mechanism from inner-sphere to outer-sphere coupling between alkyl radicals and the active Cu(II) catalyst with increasing substitution of the alkyl radical. The combination of computational predictions and detailed experimental validation shows that simultaneous minimization of both Cu-C covalency and alkyl radical size increases the rate of reductive elimination and that both strongly electron-donating and electron-withdrawing substituents on the catalyst accelerate the selectivity-determining C-N bond formation process as a result of a change in mechanism. These findings offer design principles for the development of improved catalyst scaffolds for radical C-H functionalization reactions.

Chemical Bonding and the Role of Node-Induced Electron Confinement.

The chemical bond is the cornerstone of chemistry, providing a conceptual framework to understand and predict the behavior of molecules in complex systems. However, the fundamental origin of chemical bonding remains controversial and has been responsible for fierce debate over the past century. Here, we present a unified theory of bonding, using a separation of electron delocalization effects from orbital relaxation to identify three mechanisms [node-induced confinement (typically associated with Pauli repulsion, though more general), orbital contraction, and polarization] that each modulate kinetic energy during bond formation. Through analysis of a series of archetypal bonds, we show that an exquisite balance of energy-lowering delocalizing and localizing effects are dictated simply by atomic electron configurations, nodal structure, and electronegativities. The utility of this unified bonding theory is demonstrated by its application to explain observed trends in bond strengths throughout the periodic table, including main group and transition metal elements.

Femtosecond Core-Level Spectroscopy Reveals Involvement of Triplet States in the Gas-Phase Photodissociation of Fe(CO)5.

Excitation of iron pentacarbonyl [Fe(CO)5], a prototypical photocatalyst, at 266 nm causes the sequential loss of two CO ligands in the gas phase, creating catalytically active, unsaturated iron carbonyls. Despite numerous studies, major aspects of its ultrafast photochemistry remain unresolved because the early excited-state dynamics have so far eluded spectroscopic observation. This has led to the long-held assumption that ultrafast dissociation of gas-phase Fe(CO)5 proceeds exclusively on the singlet manifold. Herein, we present a combined experimental-theoretical study employing ultrafast extreme ultraviolet transient absorption spectroscopy near the Fe M2,3-edge, which features spectral evolution on 100 fs and 3 ps time scales, alongside high-level electronic structure theory, which enables characterization of the molecular geometries and electronic states involved in the ultrafast photodissociation of Fe(CO)5. We assign the 100 fs evolution to spectroscopic signatures associated with intertwined structural and electronic dynamics on the singlet metal-centered states during the first CO loss and the 3 ps evolution to the competing dissociation of Fe(CO)4 along the lowest singlet and triplet surfaces to form Fe(CO)3. Calculations of transient spectra in both singlet and triplet states as well as spin-orbit coupling constants along key structural pathways provide evidence for intersystem crossing to the triplet ground state of Fe(CO)4. Thus, our work presents the first spectroscopic detection of transient excited states during ultrafast photodissociation of gas-phase Fe(CO)5 and challenges the long-standing assumption that triplet states do not play a role in the ultrafast dynamics.

Geometric Tuning of Coordinatively Unsaturated Copper(I) Sites in Metal-Organic Frameworks for Ambient-Temperature Hydrogen Storage.

Porous solids can accommodate and release molecular hydrogen readily, making them attractive for minimizing the energy requirements for hydrogen storage relative to physical storage systems. However, H2 adsorption enthalpies in such materials are generally weak (-3 to -7 kJ/mol), lowering capacities at ambient temperature. Metal-organic frameworks with well-defined structures and synthetic modularity could allow for tuning adsorbent-H2 interactions for ambient-temperature storage. Recently, Cu2.2Zn2.8Cl1.8(btdd)3 (H2btdd = bis(1H-1,2,3-triazolo-[4,5-b],[4',5'-i])dibenzo[1,4]dioxin; CuI-MFU-4l) was reported to show a large H2 adsorption enthalpy of -32 kJ/mol owing to π-backbonding from CuI to H2, exceeding the optimal binding strength for ambient-temperature storage (-15 to -25 kJ/mol). Toward realizing optimal H2 binding, we sought to modulate the π-backbonding interactions by tuning the pyramidal geometry of the trigonal CuI sites. A series of isostructural frameworks, Cu2.7M2.3X1.3(btdd)3 (M = Mn, Cd; X = Cl, I; CuIM-MFU-4l), was synthesized through postsynthetic modification of the corresponding materials M5X4(btdd)3 (M = Mn, Cd; X = CH3CO2, I). This strategy adjusts the H2 adsorption enthalpy at the CuI sites according to the ionic radius of the central metal ion of the pentanuclear cluster node, leading to -33 kJ/mol for M = ZnII (0.74 Å), -27 kJ/mol for M = MnII (0.83 Å), and -23 kJ/mol for M = CdII (0.95 Å). Thus, CuICd-MFU-4l provides a second, more stable example of optimal H2 binding energy for ambient-temperature storage among reported metal-organic frameworks. Structural, computational, and spectroscopic studies indicate that a larger central metal planarizes trigonal CuI sites, weakening the π-backbonding to H2.

Selective Adsorption of Oxygen from Humid Air in a Metal-Organic Framework with Trigonal Pyramidal Copper(I) Sites.

High or enriched-purity O2 is used in numerous industries and is predominantly produced from the cryogenic distillation of air, an extremely capital- and energy-intensive process. There is significant interest in the development of new approaches for O2-selective air separations, including the use of metal-organic frameworks featuring coordinatively unsaturated metal sites that can selectively bind O2 over N2 via electron transfer. However, most of these materials exhibit appreciable and/or reversible O2 uptake only at low temperatures, and their open metal sites are also potential strong binding sites for the water present in air. Here, we study the framework CuI-MFU-4l (CuxZn5-xCl4-x(btdd)3; H2btdd = bis(1H-1,2,3-triazolo[4,5-b],[4',5'-i])dibenzo[1,4]dioxin), which binds O2 reversibly at ambient temperature. We develop an optimized synthesis for the material to access a high density of trigonal pyramidal CuI sites, and we show that this material reversibly captures O2 from air at 25 °C, even in the presence of water. When exposed to air up to 100% relative humidity, CuI-MFU-4l retains a constant O2 capacity over the course of repeated cycling under dynamic breakthrough conditions. While this material simultaneously adsorbs N2, differences in O2 and N2 desorption kinetics allow for the isolation of high-purity O2 (>99%) under relatively mild regeneration conditions. Spectroscopic, magnetic, and computational analyses reveal that O2 binds to the copper(I) sites to form copper(II)-superoxide moieties that exhibit temperature-dependent side-on and end-on binding modes. Overall, these results suggest that CuI-MFU-4l is a promising material for the separation of O2 from ambient air, even without dehumidification.

Spin parameter optimization for spin-polarized extended tight-binding methods.

We present an optimization strategy for atom-specific spin-polarization constants within the spin-polarized GFN2-xTB framework, aiming to enhance the accuracy of molecular simulations. We compare a sequential and global optimization of spin parameters for hydrogen, carbon, nitrogen, oxygen, and fluorine. Sensitivity analysis using Sobol indices guides the identification of the most influential parameters for a given reference dataset, allowing for a nuanced understanding of their impact on diverse molecular properties. In the case of the W4-11 dataset, substantial error reduction was achieved, demonstrating the potential of the optimization. Transferability of the optimized spin-polarization constants over different properties, however, is limited, as we demonstrate by applying the optimized parameters on a set of singlet-triplet gaps in carbenes. Further studies on ionization potentials and electron affinities highlight some inherent limitations of current extended tight-binding methods that can not be resolved by simple parameter optimization. We conclude that the significantly improved accuracy strongly encourages the present re-optimization of the spin-polarization constants, whereas the limited transferability motivates a property-specific optimization strategy.

Quantum chemical modeling of hydrogen binding in metal-organic frameworks: validation, insight, predictions and challenges.

A detailed chemical understanding of H2 interactions with binding sites in the nanoporous crystalline structure of metal-organic frameworks (MOFs) can lay a sound basis for the design of new sorbent materials. Computational quantum chemical calculations can aid in this quest. To set the stage, we review general thermodynamic considerations that control the usable storage capacity of a sorbent. We then discuss cluster modeling of H2 ligation at MOF binding sites using state-of-the-art density functional theory (DFT) calculations, and how the binding can be understood using energy decomposition analysis (EDA). Employing these tools, we illustrate the connections between the character of the MOF binding site and the associated adsorption thermodynamics using four experimentally characterized MOFs, highlighting the role of open metal sites (OMSs) in accessing binding strengths relevant to room temperature storage. The sorbents are MOF-5, with no open metal sites, Ni2(m-dobdc), containing Lewis acidic Ni(II) sites, Cu(I)-MFU-4l, containing π basic Cu(I) sites and V2Cl2.8(btdd), also containing π-basic V(II) sites. We next explore the potential for binding multiple H2 molecules at a single metal site, with thermodynamics useful for storage at ambient temperature; a materials design goal which has not yet been experimentally demonstrated. Computations on Ca2+ or Mg2+ bound to catecholate or Ca2+ bound to porphyrin show the potential for binding up to 4 H2; there is precedent for the inclusion of both catecholate and porphyrin motifs in MOFs. Turning to transition metals, we discuss the prediction that two H2 molecules can bind at V(II)-MFU-4l, a material that has been synthesized with solvent coordinated to the V(II) site. Additional calculations demonstrate binding three equivalents of hydrogen per OMS in Sc(I) or Ti(I)-exchanged MFU-4l. Overall, the results suggest promising prospects for experimentally realizing higher capacity hydrogen storage MOFs, if nontrivial synthetic and desolvation challenges can be overcome. Coupled with the unbounded chemical diversity of MOFs, there is ample scope for additional exploration and discovery.

Occupied-Virtual Orbitals for Chemical Valence with Applications to Charge Transfer in Energy Decomposition Analysis.

In this article, we introduce the occupied-virtual orbitals for chemical valence (OVOCV). The OVOCVs can replace or complement the closely related idea of the natural orbitals for chemical valence (NOCV). The input is a difference density matrix connecting any initial single determinant to any final determinant, at a given molecular geometry, and a given one-particle basis. This arises in problems such as orbital rearrangement or charge transfer (CT) in energy decomposition analysis (EDA). The OVOCVs block-diagonalize the density difference operator into 2 × 2 blocks, which are spanned by one level that is filled in the initial state (the occupied OVOCV) and one that is empty (the virtual OVOCV). By contrast, the NOCVs fully diagonalize the density difference matrix and therefore are orbitals with mixed occupied-virtual character. Use of the OVOCVs makes it much easier to identify the donor and acceptor orbitals. We also introduce two different types of EDA methods with the OVOCVs and, most importantly, a charge decomposition analysis method that fixes the unreasonably large CT amount obtained directly from NOCV analysis. The square of the CT amount associated with each NOCV pair emerges as the appropriate value from the OVOCV analysis. When connecting the same initial and final states, this value is identical to the CT amount obtained from the independent absolutely localized molecular orbital (ALMO) complementary occupied-virtual orbital pair (COVP) analysis. The total, summed over all pairs, is also exactly the same as the independently suggested excitation number, as proved herein. Several examples are presented to compare NOCVs and OVOCVs: stretched H2+, a strong halogen bond between tetramethylthiourea and iodine, coordination of ethene in Zeise's salt, and binding in the Cp3La···C≡NCy complex.

Generalization of One-Center Nonorthogonal Configuration Interaction Singles to Open-Shell Singlet Reference States: Theory and Application to Valence-Core Pump-Probe States in Acetylacetone.

We formulate a one-center nonorthogonal configuration interaction singles (1C-NOCIS) theory for the computation of core excited states of an initial singlet state with two unpaired electrons. This model, which we refer to as 1C-NOCIS two-electron open-shell (2eOS), is appropriate for computing the K-edge near-edge X-ray absorption spectra (NEXAS) of the valence excited states of closed-shell molecules relevant to pump-probe time-resolved (TR) NEXAS experiments. With the inclusion of core-hole relaxation effects and explicit spin adaptation, 1C-NOCIS 2eOS requires mild shifts to match experiment, is free of artifacts due to spin contamination, and can capture the high-energy region of the spectrum beyond the transitions into the singly occupied molecular orbitals (SOMOs). Calculations on water and thymine illustrate the different key features of excited-state NEXAS, namely, the core-to-SOMO transitions as well as shifts and spin-splittings in the transitions analogous to those of the ground state. Simulations of the TR-NEXAS of acetylacetone after excitation to its π → π* singlet excited state at the carbon K-edge, an experiment carried out recently, showcase the ability of 1C-NOCIS 2eOS to efficiently simulate NEXAS based on nonadiabatic molecular dynamics simulations.

Real-Space Pseudopotential Method for the Calculation of Third-Row Elements X-ray Photoelectron Spectroscopic Signatures.

X-ray photoelectron spectroscopy (XPS) is a powerful characterization technique that unveils subtle chemical environment differences via core-electron binding energy (CEBE) analysis. We extend the development of real-space pseudopotential methods to calculating 1s, 2s, and 2p3/2 CEBEs of third-row elements (S, P, and Si) within the framework of Kohn-Sham density-functional theory (KS-DFT). The new approach systematically prevents variational collapse and simplifies core-excited orbital selection within dense energy level distributions. However, careful error cancellation analysis is required to achieve accuracy comparable to all-electron methods and experiments. Combined with real-space KS-DFT implementation, this development enables large-scale simulations with both Dirichlet boundary conditions and periodic boundary conditions.

Understanding ion-transfer reactions in silver electrodissolution and electrodeposition from first-principles calculations and experiments.

The electrified aqueous/metal interface is critical in controlling the performance of energy conversion and storage devices, but an atomistic understanding of even basic interfacial electrochemical reactions challenges both experiment and computation. We report a combined simulation and experimental study of (reversible) ion-transfer reactions involved in anodic Ag corrosion/deposition, a model system for interfacial electrochemical processes generating or consuming ions. With the explicit modeling of the electrode potential and a hybrid implicit-explicit solvation model, the density functional theory calculations produce free energy curves predicting thermodynamics, kinetics, partial charge profiles, and reaction trajectories. The calculated (equilibrium) free energy barriers (0.2 eV), and their asymmetries, agree with experimental activation energies (0.4 eV) and transfer coefficients, which were extracted from temperature-dependent voltage-step experiments on Au-supported, Ag-nanocluster substrates. The use of Ag nanoclusters eliminates the convolution of the kinetics of Ag+(aq.) generation and transfer with those of nucleation or etch-pit formation. The results indicate that the barrier is controlled by the bias-dependent competition between partial solvation of the incipient ion, metal-metal bonding, and electrostatic stabilization by image charge, with the latter two factors weakened by stronger positive biases. We also report simulations of the bias-dependence of defect generation relevant to nucleating corrosion by removing an atom from a perfect Ag(100) surface, which is predicted to occur via a vacancy-adatom intermediate. Together, these experiments and calculations provide the first validated, accurate, molecular model of the central steps that govern the rates of important dissolution/deposition reactions broadly relevant across the energy sciences.

Highly Accurate Prediction of NMR Chemical Shifts from Low-Level Quantum Mechanics Calculations Using Machine Learning.

Theoretical predictions of NMR chemical shifts from first-principles can greatly facilitate experimental interpretation and structure identification of molecules in gas, solution, and solid-state phases. However, accurate prediction of chemical shifts using the gold-standard coupled cluster with singles, doubles, and perturbative triple excitations [CCSD(T)] method with a complete basis set (CBS) can be prohibitively expensive. By contrast, machine learning (ML) methods offer inexpensive alternatives for chemical shift predictions but are hampered by generalization to molecules outside the original training set. Here, we propose several new ideas in ML of the chemical shift prediction for H, C, N, and O that first introduce a novel feature representation, based on the atomic chemical shielding tensors within a molecular environment using an inexpensive quantum mechanics (QM) method, and train it to predict NMR chemical shieldings of a high-level composite theory that approaches the accuracy of CCSD(T)/CBS. In addition, we train the ML model through a new progressive active learning workflow that reduces the total number of expensive high-level composite calculations required while allowing the model to continuously improve on unseen data. Furthermore, the algorithm provides an error estimation, signaling potential unreliability in predictions if the error is large. Finally, we introduce a novel approach to keep the rotational invariance of the features using tensor environment vectors (TEVs) that yields a ML model with the highest accuracy compared to a similar model using data augmentation. We illustrate the predictive capacity of the resulting inexpensive shift machine learning (iShiftML) models across several benchmarks, including unseen molecules in the NS372 data set, gas-phase experimental chemical shifts for small organic molecules, and much larger and more complex natural products in which we can accurately differentiate between subtle diastereomers based on chemical shift assignments.

Precursor Chemistry of Lead Bromide Perovskite Nanocrystals.

We investigate the early stages of cesium lead bromide perovskite formation through absorption spectroscopy of stopped-flow reactions, high-throughput mapping, and direct synthesis and titration of potential precursor species. Calorimetric and spectroscopic measurements of lead bromide complex titrations combined with theoretical calculations suggest that bromide complexes with higher coordination numbers than previously considered for nonpolar systems can better explain observed behaviors. Synthesis mapping of binary lead halides reveals multiple lead bromide species with absorption peaks higher than 300 nm, including a previously observed species with a peak at 313 nm and two species with peaks at 345 and 370 nm that also appear as reaction intermediates during the formation of lead bromide perovskites. Based on theoretical calculations of excitonic energies that match within 50 meV, we give a preliminary assignment of these species as two-dimensional magic-sized clusters with side lengths of 2, 3, and 4 unit cells. Kinetic measurements of the conversion of benzoyl bromide precursor are connected to stopped-flow measurements of product formation and demonstrate that the formation of complexes and magic-sized clusters (i.e., nucleation) is controlled by precursor decomposition, whereas the growth rate of 2D and 3D perovskites is significantly slower.