Search NASA⌕ Search

SEARCH · Search NASA

Results for “optimized energy density functionals”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Defect Chemistry and Hydrogen Transport in La/Sr-Based Oxyhydrides

Oxyhydrides in the series La 2 – y Sr y LiH 1 + y O 3 – y are attractive materials for solid-state hydride electrolytes. However, their chemical stability is poor, and thus far only La 2 LiHO 3 (y = 0) has been demonstrated in a working fuel cell. Using a first-principles approach to study defect chemistry and stability, we find that La 2 LiHO 3 can be stabilized over a range of H-rich, O-poor chemical conditions, while Sr 2 LiH 3 O has low stability. Here, defects are critical for ionic transport and for stability, particularly in Sr 2 LiH 3 O, in which anion antisites and H vacancies readily form. In fact, we show that the ground-state structure of Sr 2 LiH 3 O is inherently disordered. La 2 LiHO 3 has low conductivity, which we connect to the high formation energy of defects such as VO and OH that are necessary for hydride conduction. When Sr 2 LiH 3 O is grown under carefully selected conditions, V H + defects can be made prevalent; these defects give rise to much higher ionic conductivity than can be obtained in La 2 LiHO 3 . In general, we show that the choice of synthesis conditions is vital when seeking to optimize the stability and ionic conductivity of these oxyhydrides, thereby tailoring them for use in solid-state hydrogen fuel cells.

08 HYDROGEN↗
DOE OSTI2022

Tuning bandgap and energy stability of Organic-Inorganic halide perovskites through surface engineering

Organohalide perovskite with a variety of surface structures and morphologies have shown promising potential owing to the choice of the type of heterostructure dependent stability. We systematically investigate and discuss the impact of 2-dimensional molybdenum-disulphide (MoS 2 ), molybdenum-diselenide (MoSe 2 ), tungsten-disulphide (WS 2 ), tungsten-diselenide (WSe 2 ), boron-nitiride (BN) and graphene monolayers on bandgap and energy stability of organic–inorganic halide perovskites. We found that MAPbI 3 deposited on BN-ML shows room temperature stability (-25 meV ~ 300 K) with an optimal bandgap of ~ 1.68 eV. The calculated absorption coefficient also lies in the visible-light range with a maximum of 4.9 × 10 4 cm –1 achieved at 2.8 eV photon energy. Furthermore, on the basis of our calculations, we suggest that the encapsulation of an organic–inorganic halide perovskite monolayers by semiconducting monolayers potentially provides greater flexibility for tuning the energy stability and the bandgap.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Molecular insights into CO 2 -to-bicarbonate transformation in functionalized anion exchange ionomers for electrochemical separations

Bipolar membrane (BPM) electrochemical processes are a promising platform for carbon dioxide (CO 2 ) separations, but the molecular level thermodynamic and kinetic understanding of CO 2 -to-bicarbonate (HCO 3 − ) transformation remain poorly understood. This study employs a multiscale computational approach to systematically explore the adsorption and reactive transformation of CO 2 in five anion exchange ionomer systems. Classical molecular dynamics (MD) simulation results demonstrate that polymers with imidazolium groups significantly reduce CO 2 diffusion and enhance (OH − )–CO 2 interactions due to stronger electrostatic and π-interactions. Compared to the commonly used quaternary ammonium ionomers, imidazolium-functionalized ionomers show improved CO 2 proximity and interaction strength. Ab initio MD and density functional theory (DFT) calculations reveal that the benzyl-substituted imidazolium (IM-Ben) substantially reduces the energy barrier for HCO 3 − formation (∼72 meV lower) compared to the alkyl-substituted IM-nBu, while also mitigating imidazolium deprotonation under moderate hydration conditions. Transition state analysis shows IM-Ben forms more extensive hydrogen-bonding networks, which stabilize the transition state structure and contribute to a lower energy barrier for bicarbonate formation. These findings highlight the advantage of the adjacent benzyl moiety in enabling efficient CO 2 -to-bicarbonate transformation via hydrated hydroxide ion counterions, offering mechanistic insights and clear molecular design principles for optimizing anion exchange ionomers at bipolar membrane interfaces for electrochemical CO 2 separation applications.

Bipolar membranes, Reactive transformation of CO2,↗
DOE OSTI2022

Relativistic nonorthogonal configuration interaction: application to L 2,3 -edge X-ray spectroscopy

In this article, we develop a relativistic exact-two-component nonorthogonal configuration interaction (X2C-NOCI) for computing L-edge X-ray spectra. This article to our knowledge is the first time NOCI has been used for relativistic wave functions. A set of molecular complexes, including SF 6 , SiCl 4 and [FeCl 6 ] 3– , are used to demonstrate the accuracy and computational scaling of the X2C-NOCI method. Here, our results suggest that X2C-NOCI is able to satisfactorily capture the main features of the L 2,3 -edge X-ray absorption spectra. Excitations from the core require a large amount of orbital relaxation to yield reasonable energies and X2C-NOCI allows us to treat orbital optimization explicitly. However, the cost of computing the nonorthogonal coupling is higher than in conventional CI. Here, we propose an improved integral screening using overlap-scaled density combined with a continuous measure of the generalized Slater–Condon rules that allows us to estimate if an element is zero before attempting a two-electron integral contraction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Role of opacity at the 9 keV back lighter energy used in measuring the equation of state of boron at pressures up to a Gbar

Many experiments have been conducted at the National Ignition Facility to measure the equation of state Hugoniot of plastic, boron, and diamond at extreme pressures up to a Gbar. The “Gbar” design employs a strong spherically converging shock launched through a solid ball of material using a hohlraum radiation drive. The shock front conditions are characterized using X-ray radiography, typically at energies near 9 keV. Here we examine how the opacity of boron at 9 keV changes at high pressures and temperatures. Understanding this is vital to unfolding the density in the shock front as pressures exceed 100 Mbar. We compare opacity calculations from a number of methods including the legacy XSN opacity tables, super transition array (STA), average atom (AVAT), and detailed configuration accounting (DCA) methods. We examine how the changing opacity is correlated with the K-shell occupation calculated using the opacity methods as well as other electronic structure methods such as MECCA, Purgatorio, and an all-electron Kohn-Sham density functional theory calculation that uses optimized all-electron norm-conserving Vanderbilt pseudopotentials (ONCV). We also examine the relative contribution of free-free opacity to the overall opacity, which is dominated by bound-free absorption.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Unusual Role of Point Defects in Perovskite Nickelate Electrocatalysts

Low-cost transition-metal oxide is regarded as a promising electrocatalyst family for an oxygen evolution reaction (OER). The classic design principle for an oxide electrocatalyst believes that point defect engineering, such as oxygen vacancies (VO .. ) or heteroatom doping, offers the opportunities to manipulate the electronic structure of material toward optimal OER activity. Oppositely, in this work, we discover a counterintuitive phenomenon that both VO .. and an aliovalent dopant (i.e., proton (H + )) in perovskite nickelate (i.e., NdNiO 3 (NNO)) have a considerably detrimental effect on intrinsic OER performance. Detailed characterizations unveil that the introduction of these point defects leads to a decrease in the oxidative state of Ni and weakens Ni-O orbital hybridization, which triggers the local electron-electron correlation and a more insulating state. Evidenced by first-principles calculation using the density functional theory (DFT) method, the OER on nickelate electrocatalysts follows the lattice oxygen mechanism (LOM). The incorporation of point defect increases the energy barrier of transformation from OO*(V O ) to OH*(V O ) intermediates, which is regarded as the rate-determining step (RDS). Furthermore, this work offers a new and significant perspective of the role that lattice defects play in the OER process.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

First-principles investigation of the resistive switching energetics in monolayer MoS 2 : insights into metal diffusion and adsorption

A deeper understanding of resistive switching (RS) in 2D materials is essential for advancing neuromorphic computing. The Dissociation-Diffusion-Adsorption (DDA) model offers a useful framework for probing RS mechanisms in non-volatile memory (NVM) and in-memory computing. We have employed first-principles density functional theory (DFT) to explore dissociation, diffusion, and adsorption phenomena within the DDA model, focusing on the interactions between exemplary metal atoms (Au, Ag, Cu) and monolayer MoS 2 . Nudged elastic band (NEB) calculations evaluated diffusion barriers in pristine and sulfur-vacancy MoS 2 . Charged systems were modeled to assess the impact of applied bias on migration pathways. We also examined metal dissociation from bulk electrodes and adsorption at S vacancies. Ag/MoS 2 shows the lowest dissociation barrier (~0.034 eV), while Au and Cu exhibit similar values (~0.32 eV). These insights highlight Ag as a promising candidate for low-energy RS applications and provide guidance for optimizing switching efficiency in 2D memory devices.

Atomistic models↗
DOE OSTI2020

Revisiting Competing Paths in Electrochemical CO 2 Reduction on Copper via Embedded Correlated Wavefunction Theory

We re-evaluate two key steps in the mechanism of CO 2 reduction on copper at a higher level of theory capable of correcting inherent errors in density functional theory (DFT) approximations, namely, embedded correlated wavefunction (ECW) theory. Here, we consider the CO reduction step on Cu(111), which is critical to understanding reaction selectivity. Here, we optimize embedding potentials at the periodic plane-wave DFT level using density functional embedding theory (DFET). All possible adsorption sites (adsites) for each adsorbate then are screened with ECW theory at the catalytically active site to refine the local electronic structure. Unsurprisingly, DFT and ECW theory predict different adsite preferences, largely because of DFT’s inability to properly situate the CO 2π* level. Differing preferred adsites suggest that different reaction pathways could emerge from DFT versus ECW theory. Starting from these preferred ECW theory adsites, we then obtain reaction pathways at the plane-wave DFT level using the climbing-image nudged elastic band method to determine minimum energy paths. Thereafter, we perform ECW calculations at the catalytically active site to correct the energetics at each interpolated structure (image) along the reaction pathways. Via this approach, we confirm that the first step in CO reduction via hydrogen transfer on Cu(111) is to form hydroxymethylidyne (*COH) instead of formyl (*CHO). Although the prediction to preferentially form *COH is consistent with that of DFT, the two theories predict quite different structural and mechanistic behaviors, suggesting that verification is needed for other parts of the mechanism of CO 2 reduction, which is the subject of ongoing work.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2023

Two excited-state datasets for quantum chemical UV-vis spectra of organic molecules

Abstract We present two open-source datasets that provide time-dependent density-functional tight-binding (TD-DFTB) electronic excitation spectra of organic molecules. These datasets represent predictions of UV-vis absorption spectra performed on optimized geometries of the molecules in their electronic ground state. The GDB-9-Ex dataset contains a subset of 96,766 organic molecules from the original open-source GDB-9 dataset. The ORNL_AISD-Ex dataset consists of 10,502,904 organic molecules that contain between 5 and 71 non-hydrogen atoms. The data reveals the close correlation between the magnitude of the gaps between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO), and the excitation energy of the lowest singlet excited state energies quantitatively. The chemical variability of the large number of molecules was examined with a topological fingerprint estimation based on extended-connectivity fingerprints (ECFPs) followed by uniform manifold approximation and projection (UMAP) for dimension reduction. Both datasets were generated using the DFTB+ software on the “Andes” cluster of the Oak Ridge Leadership Computing Facility (OLCF).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Coexisting structural disorder and robust spin-polarization in half-metallic FeMnVAl

We report half-metallic ferromagnets (HMFs) are one of the most promising materials in the field of spintronics due to their unique band structure, consisting of one spin subband having metallic characteristics along with another subband with semiconductor-like behavior. In this work, we report the synthesis of a novel quaternary Heusler alloy—FeMnVAl—and we have studied the structural, magnetic, transport, and electronic properties complemented with first-principles calculations. Among different possible structurally ordered arrangements, the optimal structure is identified by theoretical energy minimization. The corresponding spin-polarized band-structure calculations indicate the presence of a half-metallic ferromagnetic ground state. A detailed and careful investigation of the x-ray diffraction data, Mössbauer spectra, and nuclear magnetic resonance spectra suggests the presence of site disorder between the Fe and Mn atoms in the stable ordered structure of the system. The magnetic susceptibility measurement clearly establishes a ferromagnetic-like transition below ~213 K. The Fe Mössbauer spectrometry measurements suggest that only the Mn spins could be responsible for the magnetic order, which is consistent with our theoretical calculation. Surprisingly, the density-functional-theory calculations reveal that the spin-polarization value is almost immunized (92.4% → 90.4%) from the Mn-Fe structural disorder, even when nonmagnetic Fe and moment-carrying Mn sites are entangled inseparably. The robustness of spin polarization and half-metallicity in the studied FeMnVAl compound comprising structural disorder is thus quite interesting and could provide a new direction to investigate and understand the exact role of disorders on spin polarization in this class of materials, beyond the available knowledge.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

FINETUNA: fine-tuning accelerated molecular simulations

Abstract Progress towards the energy breakthroughs needed to combat climate change can be significantly accelerated through the efficient simulation of atomistic systems. However, simulation techniques based on first principles, such as density functional theory (DFT), are limited in their practical use due to their high computational expense. Machine learning approaches have the potential to approximate DFT in a computationally efficient manner, which could dramatically increase the impact of computational simulations on real-world problems. However, they are limited by their accuracy and the cost of generating labeled data. Here, we present an online active learning framework for accelerating the simulation of atomic systems efficiently and accurately by incorporating prior physical information learned by large-scale pre-trained graph neural network models from the Open Catalyst Project. Accelerating these simulations enables useful data to be generated more cheaply, allowing better models to be trained and more atomistic systems to be screened. We also present a method of comparing local optimization techniques on the basis of both their speed and accuracy. Experiments on 30 benchmark adsorbate-catalyst systems show that our method of transfer learning to incorporate prior information from pre-trained models accelerates simulations by reducing the number of DFT calculations by 91%, while meeting an accuracy threshold of 0.02 eV 93% of the time. Finally, we demonstrate a technique for leveraging the interactive functionality built in to Vienna ab initio Simulation Package (VASP) to efficiently compute single point calculations within our online active learning framework without the significant startup costs. This allows VASP to work in tandem with our framework while requiring 75% fewer self-consistent cycles than conventional single point calculations. The online active learning implementation, and examples using the VASP interactive code, are available in the open source FINETUNA package on Github.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Nanometer Resolution Structure‐Emission Correlation of Individual Quantum Emitters via Enhanced Cathodoluminescence in Twisted Hexagonal Boron Nitride

Understanding the atomic structure of quantum emitters, often originating from point defects or impuritie, is essential for designing and optimizing materials for quantum technologies such as quantum computing, communication, and sensing. Despite the availability of atomic-resolution scanning transmission electron microscopy and nanoscale cathodoluminescence microscopy, experimentally determining the atomic structure of individual emitters is challenging due to the conflicting needs for thick samples to generate strong cathodoluminescence signals and thin samples for structural analysis. To overcome this challenge, significantly enhanced cathodoluminescence at twisted interfaces is leveraged to achieve sub-nanometer localization precision for the first time in mapping individual quantum emitters in carbon-implanted hexagonal boron nitride. This unprecedent spatial sensitivity, together with correlative electron energy loss spectroscopy quantitative scanning transmission electron microscopy imaging, and first principles density functional theory calculations, enables the identification of the atomic structure of the 440 nm blue emitter in hexagonal boron nitride as a substituted vertical carbon dimer. Building on the atomic structure insights, nanoscale spatially precise creation of blue emitters is demonstrated by electron beam irradiation of carbon-coated hexagonal boron nitride. This advancement in correlating atomic structures with optical properties lays the foundation for a deeper understanding and precise engineering of quantum emitters, significantly advancing the development of cutting-edge quantum information technologies.

2D material↗
DOE OSTI2026

A scalable superconducting nanowire memory array with row–column addressing

Scalable superconducting memory is required for the development of low-energy superconducting computers and fault-tolerant quantum computers. Conventional superconducting logic-based memory cells possess a large footprint that limits scaling; nanowire-based superconducting memory cells, although more compact, have high error rates, which hinders integration into large arrays. Here we report a 4 × 4 superconducting nanowire memory array that is designed for scalable row–column operations and has a functional density of 2.6 Mbit cm −2 . Each memory cell is based on a nanowire loop consisting of two temperature-dependent superconducting switches and a variable kinetic inductor. The arrays operate at 1.3 K, where we implement and characterize multiflux quanta state storage and destructive read-out. By optimizing the write- and read-pulse sequences, we minimize bit errors and maximize operating margins. We achieve a minimum bit error rate of 10 −5 . Here, we also use circuit-level simulations to understand the memory cell’s dynamics, performance limits and stability under varying pulse amplitudes.

Electrical and electronic engineering↗
DOE OSTI2024

Atomic structure of different surface terminations of polycrystalline ZnPd

The intermetallic compound ZnPd has been found to have desirable characteristics as a catalyst for the steam reforming of methanol. The understanding of the surface structure of ZnPd is important to optimize its catalytic behavior. However, due to the lack of bulk single-crystal samples and the complexity of characterizing surface properties in the available polycrystalline samples using common experimental techniques, all previous surface science studies of this compound have been performed on surface alloy samples formed through thin-film deposition. In this study, we present findings on the chemical and atomic structure of the surfaces of bulk polycrystalline ZnPd studied by a variety of complementary experimental techniques, including scanning tunneling microscopy (STM), x-ray photoelectron spectroscopy (XPS), low energy electron microscopy (LEEM), photoemission electron microscopy (PEEM), and microspot low-energy electron diffraction ( μ -LEED). These experimental techniques, combined with density functional theory (DFT)-based thermodynamic calculations of surface free energy and detachment kinetics at the step edges, confirm that surfaces terminated by atomic layers composed of both Zn and Pd atoms are more stable than those terminated by only Zn or Pd layers. DFT calculations also demonstrate that the primary contribution to the tunneling current arises from Pd atoms, in agreement with the STM results. The formation of intermetallics at surfaces may contribute to the superior catalyst properties of ZnPd over Zn or Pd elemental counterparts. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗
DOE OSTI2022

Inverse Calculation of Burden Distribution Matrix Using B-spline Model Based PDF control in Blast Furnace Burden Charging Process

The inverse calculation of burden distribution matrix (BDM) is one of the most important challenges in the blast furnace operation in iron-making processes. In general, blast furnace consumes 65% of the total energy for the whole steel-making. Focusing on this practical challenge, this article proposes a new burden distribution spatial model in calculating burden charging process, and develops a B-spline approximation-based probability density function (PDF) control algorithm to assign the expected thickness distribution of burden layer and, thus, develops a new method for the required inverse calculation of BDM. First, a novel method for the thickness distribution of burden layer is given using B-spline model to produce an expected distribution shape subjected to a desired tracking within a specific spatial constraint. Then, according to the coexistence of continuous and bounded discrete variables in BDM, a novel hybrid optimization control method by combining integer programming and PDF tracking is further established for the effective inverse calculation of BDM. Finally, the proposed PDF-based iterative inverse calculation of BDM using B-spline models are tested using various data from industrial examples. Furthermore, the simulation results show that the proposed method is well suited to solve the BDM inverse calculation problem in practice.

42 ENGINEERING↗
NASA NTRS2021

Machine Learning the COSMO Model for Predicting Thermodynamics of Electrolyte Mixtures

Bottom-up design of electrolyte mixtures for battery systems requires predicting macro thermodynamic properties from molecular constituents. For instance, molten salt electrolyte batteries require conditions far above room temperature to operate. Therefore, discovering mixtures with increasingly lower eutectic melting points is desirable. A model that can approximate chemical activity is a valuable tool to search through the vast compositional design space. Machine learning can predict properties of materials such as vibrational free energies, electronic energy gaps, and thermal conductivities. Moreover, they can learn physical models such as interatomic potentials. The COSMO-SAC model uses theory and empirical parameterization to predict liquid-vapor and liquid-solid properties using first-principles calculations. However, obtaining activity coefficients required for parameterizing the COSMO-SAC model is costly and limited to a select chemical space. In this work, we explored if machine learning methods could improve the COSMO-SAC model and bridge density functional theory calculations to liquid phase thermodynamic properties. Our data-driven approach uses existing databases for sigma-profiles of organic solvents and reconciles their methodological differences via ensemble averaging. First, an optimal machine learning model is constructed for each dataset. Our machine learning algorithms use the sigma-profile as an input feature to predict binary mixtures' activity coefficients using multi-output regression. Each dataset uses different choices of functionals, methods, and basis sets. Therefore, our ensemble model attempts to predict corrected activity coefficients given the combination of all the model outputs. The activity coefficients used for training are generated using the COSMO-SAC model. This approach enables the extraction of meaningful information from the existing datasets to improve the COSMO-SAC model for obtaining thermodynamic properties of electrolyte mixtures. With the liquid phase activities, we can identify electrolyte mixtures that meet desired phase equilibria conditions.

Thermodynamics↗
DOE OSTI2023

The nature and suppression strategies of interfacial reactions in all-solid-state batteries

Solid-state Li batteries are promising energy storage devices owing to their high safety and high theoretical energy density. However, the serious interfacial reaction between solid state electrolytes and cathodes deteriorates the battery performance, impeding the realization of long-term cyclability. The buried nature of the interphase layer presents a significant challenge in achieving a comprehensive understanding of the underlying interfacial reaction mechanisms. Herein, we systematically explore the interfacial reaction evolutions and interphase compositions and electronic properties between the popular oxide cathodes and sulfide solid electrolytes (SSEs). This includes analysis of the chemical and electrochemical reactions between cathodes/coatings and SSEs, as well as the electrochemical self-decomposition of SSEs by thermodynamic phase equilibrium analysis. Herein, we disclose that the driving force of the electrochemical reaction at the chemical potential of μ Li is much stronger than that of the chemical reaction, which dominates the interfacial reaction. Preventing the formation of an electronically conductive interphase is crucial in inhibiting the continuous interfacial degradation during long-term cycling, which can be achieved through the optimized combination of cathodes and SSEs, as well as the introduction of functional coatings between them. Based on these findings, the percentage of molar fraction ( f ) of electronically conductive species in the formed interphase is proposed as a key factor for indicating the interfacial stability for the first time. Furthermore, we propose a specific high-throughput screening scheme to filter the functional coating materials by comprehensively evaluating their functionality. The tiered screening identifies 48 coating materials with optimal properties. The work highlights the significant roles of rational coupling of the cathodes and SSEs, and optimizing interfacial coating materials for solid-state batteries. It opens new avenues for engineering an interphase with improved interfacial compatibility to realize long-term cyclability.

25 ENERGY STORAGE↗
DOE OSTI2021

Gaseous cyclodextrin- closo -dodecaborate complexes χCD·B 12 X 12 2– (χ = α, β, and γ; X = F, Cl, Br, and I): electronic structures and intramolecular interactions

A fundamental understanding of cyclodextrin- closo -dodecaborate inclusion complexes is of great interest in supramolecular chemistry. Herein, we report a systematic investigation on the electronic structures and intramolecular interactions of perhalogenated closo -dodecaborate dianions B 12 X 12 2– (X = F, Cl, Br and I) binding to α-, β-, and γ-cyclodextrins (CDs) in the gas phase using combined negative ion photoelectron spectroscopy (NIPES) and density functional theory (DFT) calculations. The vertical detachment energy (VDE) of each complex and electronic stabilization of each dianion due to the CD binding (ΔVDE, relative to the corresponding isolated B 12 X 12 2– ) are determined from the experiments along α-, β- and γ-CD in the form of VDE (ΔVDE): 4.00 (2.10), 4.33 (2.43), and 4.30 (2.40) eV in X = F; 4.09 (1.14), 4.64 (1.69), and 4.69 (1.74) eV in X = Cl; 4.11 (0.91), 4.58 (1.38), and 4.70 (1.50) eV in X = Br; and 3.54 (0.74), 3.88 (1.08), and 4.05 (1.25) eV in X = I, respectively. All complexes have significantly higher VDEs than the corresponding isolated dodecaborate dianions with ΔVDE spanning from 0.74 eV at (α, I) to 2.43 eV at (β, F), sensitive to both host CD size and guest substituent X. DFT-optimized complex structures indicate that all B 12 X 12 2– prefer binding to the wide openings of CDs with the insertion depth and binding motif strongly dependent on the CD size and halogen X. Dodecaborate anions with heavy halogens, i.e. , X = Cl, Br, and I, are found outside of α-CD, while B 12 F 12 2– is completely wrapped by γ-CD. Partial embedment of B 12 X 12 2– into CDs is observed for the other complexes via multipronged B–X … H–O/C interlocking patterns. The simulated spectra based on the density of states agree well with those of the experiments and the calculated VDEs well reproduce the experimental trends. Molecular orbital analyses suggest that the spectral features at low binding energies originated from electrons detached from the dodecaborate dianion, while those at higher binding energies are derived from electron detachment from CDs. Energy decomposition analyses reveal that the electrostatic interaction plays a dominating role in contributing to the host–guest interactions for the X = F series partially due to the formation of a O/C–H … X–B hydrogen bonding network, and the dispersion forces gradually become important with the increase of halogen size.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗