Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “energy density functional”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Self-interaction error overbinds water clusters but cancels in structural energy differences

Here, we gauge the importance of self-interaction errors in density functional approximations (DFAs) for the case of water clusters. To this end, we used the Fermi–Löwdin orbital self-interaction correction method (FLOSIC) to calculate the binding energy of clusters of up to eight water molecules. Three representative DFAs of the local, generalized gradient, and metageneralized gradient families [i.e., local density approximation (LDA), Perdew–Burke–Ernzerhof (PBE), and strongly constrained and appropriately normed (SCAN)] were used. We find that the overbinding of the water clusters in these approximations is not a density-driven error. We show that, while removing self-interaction error does not alter the energetic ordering of the different water isomers with respect to the uncorrected DFAs, the resulting binding energies are corrected toward accurate reference values from higher-level calculations. In particular, self-interaction–corrected SCAN not only retains the correct energetic ordering for water hexamers but also reduces the mean error in the hexamer binding energies to less than 14 meV/ H 2 O from about 42 meV/ H 2 O for SCAN. By decomposing the total binding energy into many-body components, we find that large errors in the two-body interaction in SCAN are significantly reduced by self-interaction corrections. Higher-order many-body errors are small in both SCAN and self-interaction–corrected SCAN. These findings suggest that orbital-by-orbital removal of self-interaction combined with a proper DFA can lead to improved descriptions of water complexes.

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

Too big, too small, or just right? A benchmark assessment of density functional theory for predicting the spatial extent of the electron density of small chemical systems

Multipole moments are the first-order responses of the energy to spatial derivatives of the electric field strength. The quality of density functional theory prediction of molecular multipole moments thus characterizes errors in modeling the electron density itself, as well as the performance in describing molecules interacting with external electric fields. However, only the lowest non-zero moment is translationally invariant, making the higher-order moments origin-dependent. Therefore, instead of using the 3 × 3 quadrupole moment matrix, we utilize the translationally invariant 3 × 3 matrix of second cumulants (or spatial variances) of the electron density as the quantity of interest (denoted by K). The principal components of K are the square of the spatial extent of the electron density along each axis. A benchmark dataset of the principal components of K for 100 small molecules at the coupled cluster singles and doubles with perturbative triples at the complete basis set limit is developed, resulting in 213 independent K components. The performance of 47 popular and recent density functionals is assessed against this Var213 dataset. Several functionals, especially double hybrids, and also SCAN and SCAN0 predict reliable second cumulants, although some modern, empirically parameterized functionals yield more disappointing performance. The H, Li, and Be atoms, in particular, are challenging for nearly all methods, indicating that future functional development could benefit from the inclusion of their density information in training or testing protocols.

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

Absolute dissociative electron attachment cross-section measurement of difluoromethane

Dissociative electron attachment (DEA) to difluoromethane (CH 2 F 2 ) has been studied in the electron energy range of 0 to 16 eV. Two resonant states around 2 and 11.4 eV leading to three different fragment anions (F - , CHF - , and F 2 - ) are observed. Ion-yield curves of the negative ions help us to locate the position of the resonant states. In the ion-yield curve of F - ions, one small hump near 9.8 eV followed by a peak around 15.2 eV is also observed. Absolute DEA cross sections of the F - ion is measured by using the well-known relative flow technique. Dissociation channels associated with each resonant states are identified by density-functional theory (DFT) calculations of the threshold energies. The theory matches quite well within the experimental uncertainty.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

MnRhBi3: A Cleavable Antiferromagnetic Metal

This dataset contains DFT input and output files supporting the theoretical modeling in the associated publication (Chem. Mater. 2024, 36, 11306-11316). The calculations characterize MnRhBi3, an orthorhombic (Cmmm) van der Waals-layered intermetallic compound that cleaves easily between neighboring Bi layers. The dataset is organized into three calculation types: (i) Bulk: Structural relaxations of the periodic MnRhBi3 crystal in antiferromagnetic (AFM) and ferromagnetic (FM) configurations, using the vdW-DF-optB86b functional. These provide the equilibrium lattice constants, magnetic energy differences (AFM is 0.5 meV/f.u. lower than FM), and magnetic moments (4.4 µB/Mn, 0.17 µB/Rh, 0.18 µB/Bi) reported in Table 1 of the main text. (ii) Slab: Same magnetic configurations computed with an 18 Ang vacuum layer introduced between Bi layers, used to calculate the cleavage energy Ec = 0.56 J/m2 (AFM) and 0.57 J/m2 (FM), establishing MnRhBi3 as a van der Waals-layered material comparable to graphite, MoS2, and CrI3. (iii) ELF: Single-point calculation on the relaxed bulk AFM geometry with LELF=.TRUE., producing the ELFCAR file used to generate electron localization function isosurfaces and contour maps (Fig. 2, main text) showing Bi lone pairs directed into the van der Waals gaps. All folders contain CONTCAR, INCAR, KPOINTS, OUTCAR, and POSCAR. The ELF/ folder additionally contains ELFCAR. Calculations were performed using VASP 6.3.2 with PBE + vdW-DF-optB86b, PAW potentials, and an energy cutoff of 800 eV.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Predicting interface structure using the minima hopping method

Here, we adapt the minima hopping method (MHM) to the problem of interfacial structure prediction and apply it to study a canonical problem, the tilt grain boundaries in SrTiO 3 . Our method employs a hybrid approach by first exploring the potential energy surface (PES) of different grain boundary samplings with an empirical force field, among which the fifteen candidates with lower energies are then refined using ab initio density functional theory (DFT) calculations. During the exploratory stage, we bias the search using a local order parameter to primarily sample various reconstructions in the vicinity of the interface, while preserving the crystallinity of the bulk regions. We further enhance the search by incorporating initial structures with rigid body displacements to account for translational variations between bulk phases, enabling the MHM to effectively generate both stoichiometric and nonstoichiometric SrTiO 3 Σ⁢3(111)[110] and Σ⁢3(112)[110] grain boundaries. From an algorithmic standpoint, MHM outperforms earlier studies based on genetic algorithms (GA) by identifying more stable interfacial structures of several SrTiO 3 grain boundaries. The performance of the present implementation of the MHM approach is primarily limited by exploring an approximate description of the PES with a rather simple Buckingham potential. This limitation leads to variations in performance when compared to approaches utilizing more advanced surrogate PES models, such as direct DFT-PES sampling or GA with the embedded atom method (EAM). Despite the present limitations, the MHM approach is able to yield interfacial structures with comparable or lower interfacial energies in specific cases, such as Σ⁢3(111)[110] Γ=1, ±0.5 and Σ⁢3(112)[110] Γ= ±1, −2, underscoring the robustness of the MHM approach even with a simple approximation of the DFT PES. The MHM interfacial structure prediction method thus offers an efficient approach to understanding the grain boundaries and heterointerfaces at the atomic scale, providing an important prerequisite for effective materials design.

density functional theory↗
DOE OSTI2026

Structured for success: conjugated polymer binders with tailored composition and architecture for lithium-ion batteries

Conjugated polymer binders are replacing conventional binders in lithium-ion batteries. Herein, we examine how molecular engineering and hierarchical nanostructuring govern binder functionality and electrochemical performance. Lithium-ion batteries (LIBs) are the leading energy storage technology, yet enhancing their energy density and cycle life remains critical. Significant progress has been made in high-capacity anodes and high-voltage cathodes, but their performance is hindered by electrode degradation, where it is related to the behaviors of binders at the surface and interface. Conventional non-conductive binders like poly(vinylidene difluoride) (PVDF), combined with conductive additives, often fail to maintain electrical pathways under repeated volume changes. Alternatively, conjugated polymer binders have emerged as a superior alternative, simultaneously offering intrinsic conductivity, mechanical flexibility, and strong adhesion through π-conjugated backbones and functional groups. Their tunable molecular structure enables efficient electron/ion transport while mitigating electrode cracking. Additionally, the development of hierarchically ordered nanostructures in conjugated polymer binder can further enhance their electrochemical performance. This review examines the design principles of conjugated polymer binders, focusing on molecular engineering and nanostructural control to optimize their performance in high-loading electrodes, such as silicon-based anodes. By addressing key challenges in binder functionality, these advanced materials pave the way for next-generation high-energy-density LIBs.

Jin, Xiuyu↗
DOE OSTI2021

Interpretations of ground-state symmetry breaking and strong correlation in wavefunction and density functional theories

Significance The ground state of a quantum mechanical system is the lowest-energy eigenstate of the Hamiltonian. In isolation, it persists unchanged forever, with symmetries dictated by those of the Hamiltonian. But near-eigenstates of broken symmetry can persist for long times, even on the scale of human measurement. The appearance of broken symmetries of the electron density or spin density in a density functional calculation can reveal strong correlations among the electrons that are present in a symmetry-unbroken wavefunction. Symmetry breaking can arise when a wave-like fluctuation drops to zero frequency. The presented examples are the stretched hydrogen molecule, antiferromagnetism in solids, and the static charge-density wave in a low-density jellium, which is shown quantitatively to be a zero-frequency plasma wave.

42 ENGINEERING↗
DOE OSTI2026

Custom-trained Machine-learning Interatomic Potentials: ZnCl2 Aqueous Solution

This dataset was generated using an iterative active-learning strategy implemented in the ArcaNN software package (https://github.com/arcann-chem/arcann_training) to train machine-learning interatomic potentials for aqueous ZnCl2 solutions. Each active-learning cycle consisted of three stages: training, exploration, and labeling. The initial training set combined configurations generated in this work from enhanced-sampling ab initio molecular dynamics simulations with configurations from a previously reported neural-network-potential study of aqueous ZnCl2. The enhanced-sampling ab initio molecular dynamics simulations involved Zn–Cl separation and the chloride coordination number around Zn²? as collective variables. These configurations served as the seed dataset. Subsequent active-learning cycles expanded the training set by identifying and labeling configurations that were poorly represented by the current models, thereby improving coverage of ion-association states and changes in local coordination and charge-state environments relevant to the solution free-energy landscape. For all selected configurations, single-point calculations of the total energies and atomic forces were performed within density functional theory using the CP2K Quickstep module. Reference calculations employed the revPBE-D3 and r2SCAN exchange-correlation functionals. Motivated by recent work on aqueous Zn²?, the main revPBE calculations omitted D3 dispersion contributions involving Zn²?, while retaining the D3 correction for water and chloride. For comparison, fully dispersion-corrected revPBE-D3 reference calculations were also performed, with D3 applied to all species, including Zn²?. Valence electrons were treated explicitly, while core electrons were represented using norm-conserving Goedecker–Teter–Hutter pseudopotentials. The wave functions were expanded using the mixed Gaussian-and-plane-wave scheme with TZV2P-MOLOPT basis sets for all elements and a 600 Ry auxiliary plane-wave cutoff for the electron density. Self-consistent-field convergence was accelerated using the orbital-transformation and Direct Inversion in the Iterative Subspace algorithms, with a convergence threshold of 10?6. All single-point calculations were performed in periodic orthorhombic cells. The CELL_REF keyword in CP2K was used to define a fixed reference cell with a box length of 25 Å. This treatment ensured a consistent reference for configurations extracted from NpT trajectories with fluctuating cell dimensions. The resulting DFT energies and atomic forces constitute the ground-truth labels used to train the MLIPs. The resulting MLIP was trained for aqueous ZnCl2 solutions spanning concentrations from 0 to 30 molal and a broad pH range, from strongly acidic to strongly basic conditions. Representative examples of configurations included in the MLIP training dataset are provided below. These include 1) Representative configurations from the dataset labeled at the revPBE-D3 level, with D3 dispersion interactions involving Zn2+ excluded (revPBE-wo-D3). 2) Representative configurations from the dataset labeled at the fully dispersion-corrected revPBE-D3 level, with D3 interactions applied to all species, including Zn2+ (revPBE-D3). 3) Representative configurations from the dataset labeled at the r2SCAN level of theory (r2SCAN).

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗
DOE OSTI2024

Understanding Twinning and Deformation in High Entropy Alloys

A combination of high strength and high ductility has been observed in multi-principal element alloys due to twin formation attributed to low stacking fault energy (SFE). In the pursuit of low SFE alloys, a key bottleneck is the lack of understanding of the composition–SFE cor- relations that would guide tailoring SFE via alloy composition. Using density functional theory (DFT), we show that dopant radius, which have been postulated as a key descriptor for SFE in dilute alloys, does not fully explain SFE trends across different host metals. Instead, charge density is a much more central descriptor. It allows us to (1) explain contrasting SFE trends in Ni and Cu host metals due to various dopants in dilute concentrations, (2) explain the large SFE variations observed in the literature even within a given alloy composition due to the nearest neighbor environments in “model” concentrated alloys, and (3) develop a machine learning model that can be used to predict SFEs in multi-elemental alloys. This model opens a possibility to use charge density as a descriptor for predicting SFE in alloys. Furthermore, a descriptor-less machine learning (ML) model based only on charge density images extracted from density functional theory (DFT) is developed to predict stacking fault energies (SFE) in concentrated alloys. The model is based on convolutional neural networks (CNNs) as one of the promising ML techniques for dealing with complex images and data. Identification of correct descriptors is a key bottleneck to develop ML models for predicting materials properties. Often, in most ML models, textbook physical descriptors such as atomic radius, valence charge and electronegativity are used as descriptors which have limitations because these properties change in concentrated alloys when multiple elements are mixed to form a solid solution. We illustrate that, within the scope of DFT, the search for descriptors can be circumvented by electronic charge density, which is the backbone of the Kohn-Sham DFT and describes the system completely. The performance of our model is demonstrated by predicting SFE of concentrated alloys with an RMSE and R2 of 6.18 mJ/m2 and 0.87, respectively, validating the accuracy of the proposed approach.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Exploring structural transitions at grain boundaries in Nb using a generalized embedded atom interatomic potential

The advancement in experimental techniques, like the atom probe tomography and high resolution electron microscopy, is fueling interest in studying structural transformations of grain boundaries in metal and alloys to uncover correlations between mechanical properties and solute or impurity segregation to grain boundaries. Atomistic modeling is an important tool that can pinpoint the intricate dynamics of grain boundary phase transitions, but the lack of accurate interatomic potentials needed to simulate the complex dynamics of grain boundary structural transitions and identify different metastable phases has been the bottleneck. To this end, we use niobium as a model body centered cubic (BCC) metal and develop an interatomic potential to study grain boundary phase transitions. The potential for Nb is based on a generalization of the embedded atomic method potential and has sufficient flexibility to learn complex energy landscapes using a small set of training structures. We systematically test and validate the using data from ab initio density functional theory calculations and experiments. Using this potential, we calculate energies of multiple symmetric-tilt grain boundaries spanning a wide range of misorientation angles. Additionally, we explore different metastable structures of the Σ 27(552) [$1\overline{1}0$] grain boundary and use molecular dynamic simulations to study the coexistence of metastable phases and grain boundary transitions at finite temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Correlating Experimentally Determined Hydrogen Binding Energy with Hydrogen Evolution Activity over Metal Monolayers on Molybdenum Nitride

It is well established that hydrogen binding energy (HBE) is a key descriptor for hydrogen evolution reaction (HER) activity, and such a relationship is a useful tool for searching efficient and cost-effective HER catalysts. However, in almost all cases, the HBE values are obtained from density functional theory (DFT) calculations. Here in this study, temperature programmed desorption (TPD) was used to experimentally determine the HBE values of metal monolayers supported on molybdenum nitride (Mo 2 N), and electrochemical measurements were performed on the same surfaces. Combined DFT and kinetic Monte Carlo (kMC) simulations were used to validate the trend observed with TPD and the electrochemical HER activity. Depositing one monolayer of Pt on Mo 2 N led to similar HBE values seen for bulk Pt, and electrochemical measurements showed that monolayer Pt on Mo 2 N had similar HER activity to that of bulk Pt. Similar studies were also performed for monolayer Pd on Mo 2 N to confirm the correlation. This work also demonstrates that Mo 2 N is a promising support to reduce precious metal loading in HER catalysts.

08 HYDROGEN↗
DOE OSTI2020

Compositional phase stability of correlated electron materials within DFT+DMFT

Predicting the compositional phase stability of strongly correlated electron materials is an outstanding challenge in condensed matter physics, requiring precise computations of total energies. In this work, we employ the density functional theory plus dynamical mean-field theory (DFT+DMFT) formalism to address local correlations due to transition metal $\textit{d}$ electrons on compositional phase stability in the prototype rechargeable battery cathode material Li x CoO 2 , and detailed comparisons are made with the simpler DFT+$\textit{U}$ approach (i.e., the Hartree-Fock solution of the DMFT impurity problem). Local interactions are found to strongly impact the energetics of the band insulator LiCoO 2 , most significantly via the $E_g$ orbitals, which are partially occupied via hybridization with $\textit{O p}$ states. We find CoO 2 and Li 1/2 CoO 2 to be moderately correlated Fermi liquids with quasiparticle weights of 0.6–0.8 for the $T_{2g}$ states, which are most impacted by the interactions. Further, as compared to DFT+$\textit{U}$, DFT+DMFT considerably dampens the increase in total energy as $\textit{U}$ is increased, which indicates that dynamical correlations are important to describe this class of materials despite the relatively modest quasiparticle weights. Unlike DFT+$\textit{U}$, which can incorrectly drive Li x CoO 2 toward spurious phase separating or charge-ordered states, DFT+DMFT correctly captures the system's phase stability and does not exhibit a strong charge-ordering tendency. Most importantly, the error within DFT+$\textit{U}$ varies strongly as the composition changes, challenging the common practice of artificially tuning $\textit{U}$ within DFT+$\textit{U}$ to compensate the errors of Hartree-Fock. DFT+DMFT predicts the average intercalation voltage decreases relative to DFT, opposite to the result of DFT+$\textit{U}$, which would yield favorable agreement with experiment in conjunction with the overprediction of the voltage by the strongly constrained and appropriately normed DFT functional.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Small energy gap revealed in CrBr 3 by scanning tunneling spectroscopy

CrBr 3 is a layered van der Waals material with magnetic ordering down to the 2D limit. For decades, based on optical measurements, it is believed that the energy gap of CrBr 3 is in the range of 1.68–2.1 eV. However, controversial results have indicated that the band gap of CrBr 3 is possibly smaller than that. An unambiguous determination of the energy gap is critical to the correct interpretations of the experimental results of CrBr 3 . Here, we present the scanning tunneling microscopy and spectroscopy (STM/S) results of CrBr 3 thin and thick flakes exfoliated onto highly ordered pyrolytic graphite (HOPG) surfaces and density functional theory (DFT) calculations to reveal the small energy gap (peak-to-peak energy gap to be 0.57 ± 0.04 eV; or the onset signal energy gap to be 0.29 ± 0.05 eV from dI/dV spectra). Atomic resolution topography images show the defect-free crystal structure and the dI/dV spectra exhibit multiple peak features measured at 77 K. The conduction band – valence band peak pairs in the multi-peak dI/dV spectrum agree very well with all reported optical transitions. STM topography images of mono- and bi-layer CrBr 3 flakes exhibit edge degradation due to short air exposure (~15 min) during sample transfer. The unambiguously determined small energy gap settles the controversy and is the key in better understanding CrBr 3 and similar materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Electronically induced defect creation at semiconductor/oxide interface revealed by time-dependent density functional theory

Carrier induced defect creation at the semiconductor-oxide interface has been known as the origin of electronic device degradation for a long time, but how exactly the interface lattice can be damaged by carriers (especially low-energy ones) remains unclear. Here we carry out real-time time-dependent density functional theory simulations on concrete Si/SiO 2 interfaces to study the interaction between excited electrons and interface bonds. We show that the normal interface Si-H bonds are generally resistant to electrons due to the delocalized nature and high energy level of the Si-H antibonding states, and due to the high-energy barrier to break the Si-H bond. However, if an additional hydrogen atom exists by attaching to a nearby oxygen atom (forming a “Si-H···H-O” complex), the Si-H bond will be greatly weakened, including the reduction of energy barrier for bond breaking, and the lowering of the antibonding state energy level which favors electron injection. Together with the multiple vibrational excitation process, the corresponding Si-H bond can be broken much more easily. Thus we propose that the Si-H···H-O complex will be the center for defect creation and device degradation. Finally, we also explain why such a center might be relatively easy to form during the hydrogen annealing process.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Predictive Model of Charge Mobilities in Organic Semiconductor Small Molecules with Force-Matched Potentials

Charge mobility of crystalline organic semiconductors (OSC) is limited by local dynamic disorder. Recently, the charge mobility for several high mobility OSCs, including TIPS-pentacene, were accurately predicted from a density functional theory (DFT) simulation constrained by the crystal structure and the inelastic neutron scattering spectrum, which provide direct measures of the structure and the dynamic disorder in the length scale and energy range of interest. However, the computational expense required for calculating all of the atomic and molecular forces is prohibitive. In this work, we demonstrate the use of density functional tight binding (DFTB), a semiempirical quantum mechanical method that is 2 to 3 orders of magnitude more efficient than DFT. We show that force matching a many-body interaction potential to DFT derived forces yields highly accurate DFTB models capable of reproducing the low-frequency intricacies of experimental inelastic neutron scattering (INS) spectra and accurately predicting charge mobility. We subsequently predicted charge mobilities from our DFTB model of a number of previously unstudied structural analogues to TIPS-pentacene using dynamic disorder from DFTB and transient localization theory. The approach we establish here could provide a truly rapid simulation pathway for accurate materials properties prediction, in our vision applied to new OSCs with tailored properties.

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

Stability and Performance of 3d Transition Metal Carbo‐Sulfides: A Density Functional Theory Exploration for Li‐Ion Battery Anodes

As the demand for high-performance and reliable energy storage devices continues to rise, identifying new anode materials is crucial for advancing Li-ion battery (LIB) technology. Inspired by recent experimental breakthroughs in synthesizing two-dimensional transition metal carbo-chalcogenides (2D-TMCCs), density functional theory calculations are performed to systematically explore their sulfide variants (TM 2 S 2 C) spanning all 3d transition metals in three possible phases. Through comprehensive evaluations of thermodynamic, dynamic, mechanical, and thermal stabilities, seven stable 2D-TMCC candidates are identified, four of which exhibit superior battery performance. Notably, V-based 2D-TMCCs across all three phases deliver moderate open-circuit voltages (OCV), efficient Li diffusion, and substantial capacities, making them promising candidates for industrial applications without requiring specific phase controls. A Cr-based 2D-TMCC (with sulfur atoms above carbon atoms) offers the highest capacity of 515.40 mAh g −1 , the lowest Li diffusion barrier, and an optimal OCV, highlighting its appealing potential as an anode material for LIBs. Furthermore, significant Li–Li spacing and pronounced electron delocalization in these four 2D-TMCCs suggest a reduced risk of dendrite formation. This work expands the 2D-TMCC family and identifies up-and-coming candidates for next-generation LIB anodes.

anode materials↗
DOE OSTI2021

Phase Diagram of a Deep Potential Water Model

Using the Deep Potential methodology, we construct a model that reproduces accurately the potential energy surface of the SCAN approximation of density functional theory for water, from low temperature and pressure to about 2400 K and 50 GPa, excluding the vapor stability region. The computational efficiency of the model makes it possible to predict its phase diagram using molecular dynamics. Satisfactory overall agreement with experimental results is obtained. Here, the fluid phases, molecular and ionic, and all the stable ice polymorphs, ordered and disordered, are predicted correctly, with the exception of ice III and XV that are stable in experiments, but metastable in the model. The evolution of the atomic dynamics upon heating, as ice VII transforms first into ice VII" and then into an ionic fluid, reveals that molecular dissociation and breaking of the ice rules coexist with strong covalent fluctuations, explaining why only partial ionization was inferred in experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗