Search NASA⌕ Search

SEARCH · Search NASA

Results for “DFT”

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 55 records · Page 3

CO adsorption on MnO(100): Experimental benchmarks compared to DFT

In this study, CO adsorption on the MnO(100) surface was studied using temperature programmed desorption (TPD) and density functional theory (DFT). TPD results show that CO is weakly-bound on MnO(100), with an experimental adsorption energy of -35.6 ± 2.1 kJ/mol at terrace sites in the low coverage limit. PBE simulations suggest that CO adsorption causes an implausible (2 × 2) surface reconstruction. PBE+U simulations show no signs of surface reconstruction, and provide an accurate estimate of the adsorption energy (-36.4 kJ/mol) when combined with the DFT-D3 method with Becke-Jonson damping to correct for van der Waals interactions. This simulation also shows that CO adsorbs C-down onto the Mn 2+ terrace site in a tilted geometry, which is also observed experimentally and computationally on the similarly-structured NiO(100) transition metal oxide surface. TPD results for large doses show a plateauing of the coverage at about 0.4 monolayers of CO at 85 K, with a defect coverage equivalent to 0.08 monolayers. Adsorption associated with defect sites is indicated by a high-temperature desorption tail which is not satisfactorily explained by DFT simulations of simple step or oxygen vacancy defects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DFT-based mesoscopic interactions of N adatoms on Cu(100)

The first-principles density functional theory (DFT) method is used to directly calculate the interactions of N atoms of various configurations on Cu(100) surfaces. This allows us to derive a mesoscopic substratemediated interaction model. We find that, for these Cu(100)-supported two-dimensional N clusters, oscillating substrate-mediated interactions obtained from an elastic eigenvector approach including a mesoscopic interacting range (called the mesoscopic model) have reasonable agreement with the DFT results. The DFT periodic boundary conditions are properly mapped into the mesoscopic model. The model is used to analyze the energies of c(2×2) square islands independent of their sizes. Previous experimental studies have shown that, in the low coverage region, N atoms agglomerate to isolated square islands of about 5 nm × 5 nm in size with a c(2×2) structure. With increasing coverage, those islands form regular patterns but do not touch below about 0.4 monolayers. From previous experiments and calculations, the square island pattern shows significant elastic effects which are now confirmed.

74 ATOMIC AND MOLECULAR PHYSICS↗

Interaction of Methane with Calcium Ferrite in the Chemical Looping Partial Oxidation Application: Experimental and DFT Study

In this study, experimental studies using mass spectroscopy, thermogravimetric analysis, in situ X-ray diffraction (XRD), and theoretical density functional theory calculations were carried out to investigate the mechanism of methane partial oxidation over a CaFe 2 O 4 oxygen carrier. The XRD analysis revealed the formation of the calcium ferrite phases during methane reduction occurred in the sequence CaFe 2 O 4 →CaFe 3 O 5 + Ca 2 Fe 2 O 5 →FeO + Fe + CaO→Fe + CaO. Each of these phases was modeled using periodic density functional theory to study the activity of these phases for the reactivity with CH 4 , CO, and H 2 along with the adsorbed surface species to help elucidate the effect that phase changes during reduction have on the chemical looping products CO and H 2 . Mass spectroscopic analysis revealed the formation of H 2 O and CO 2 occurred prior to the formation of CO and H 2 , and these products may be due to the removal of lattice oxygen bound in the discrete primary phase, which contributes to a small portion (~16%) of the total transferable oxygen in the solid. The density functional theory (DFT) results suggest that the gaseous CO and H 2 are not thermodynamically favorable to react with the Ca 2 Fe 2 O 5 where the release of CO and H 2 led to an increase in selectivity for these products. The DFT results suggest that some of the adsorbed atomic hydrogen combined to form adsorbed H 2 and then gaseous H 2 on the CaFe 2 O 4 surface, while the carbon deposition and some adsorbed atomic hydrogen react to form the CaFe 3 O 5 and Ca 2 Fe 2 O 5 phases simultaneously. DFT simulations suggest that both the gaseous H 2 and adsorbed atomic H are less thermodynamically favorable to drive the deep reduction of Ca 2 Fe 2 O 5 to CaO and Fe0. The deep reduction of the calcium ferrite therefore favors the generation of syngas products (CO and H 2 ) through mechanisms of methane dehydrogenation and selective oxidation of carbon to CO. The significant difference in mechanisms for partial oxidation of methane to syngas with the CaFe 2 O 4 carrier can be correlated to the differences in the oxide phases present during the reduction process, concentration of oxygen species, and bonding in the lattice structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Cu(II) Stability and UV-Induced Electron Transfer in a Metal–Organic Hybrid: An EPR, DFT, and Crystallographic Characterization of Copper-Doped Zinc Creatininium Sulfate

Single-crystal X-ray diffraction and electron paramagnetic resonance (EPR) spectroscopic experiments, complemented by quantum chemical DFT calculations, were carried out on the copper-doped metal–organic hybrid and Tutton salt analogue zinc creatininium sulfate to determine its crystal structure, to characterize the electronic structure of the doped Cu(II) binding site, and to propose a pathway for an excited-state, proton-coupled electron transfer (PCET) process in UV-exposed crystals. The crystal structure is isomorphous to that of cadmium creatininium sulfate, which has the transition ion, not in direct coordination with the creatinine, but forming a hexahydrate complex, which is bridged to a creatininium through an intervening sulfate ion. The EPR g (2.446, 2.112, 2.082) and copper hyperfine (A Cu : -327, -59.6, 10.8 MHz) tensor parameters are consistent with doped copper replacing host zinc in the metal–hexahydrate complex. These parameters are similar to those observed for copper hexahydrate in doped Tutton salt systems at low temperature, where the unpaired electron occupies mainly the copper 3d x 2 –y 2 orbital. At room temperature in the Tutton systems, vibration couplings stemming from a dynamic Jahn–Teller effect cause tensor averaging which results in a reduction in their maximum g-tensor and hyperfine tensor values. However, like for the doped isomorphous Cd creatinine crystal, the Cu(II) EPR exhibits little, or no room temperature averaging compared to its low temperature pattern. Samples exposed to 254 nm UV light generate a carbon-centered free radical species, characterized by an isotropic g-tensor (g = 2.0029) and an alpha-proton hyperfine coupling (-24 -14 +4 G). These parameters identify it as a creatinine radical cation formed by the oxidative release of one of its C2 methylene hydrogens. DFT calculations confirm the unpaired electronic structures of both the Cu(II) site and free radical. The growth in radical concentration with an increase in the UV exposure time coincides with a decrease in the copper EPR signal, indicating a coupled light-induced oxidation reduction process. A comparison of the crystal structure with the EPR parameters and DFT results provides evidence for a UV-induced PCET.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Benchmarking DFT Accuracy in Predicting O 1s Binding Energies on Metals

X-ray photoelectron spectroscopy (XPS) is a powerful tool for probing the electronic structure and composition of materials, particularly metals and metal oxides of relevance to solar cells and catalysis. Density functional theory (DFT) is often used to support XPS peak assignments, but its reliability for predicting oxygen species is not well established. Here, we compile a large data set of experimental oxygen binding energies and evaluate corresponding DFT predictions. We find that as the binding energies of metal-bound atomic oxygen species increase, especially above ≈530 eV, there is a general decrease in the accuracy of DFTpredicted values. Thus, high-binding-energy atomic oxygen species, such as those proposed as active for selective Ag-catalyzed epoxidation, are less well represented. The chemical nature of the oxygen species also influences accuracy, with molecularly bound species more reliably captured across the entire range of energies. These findings illustrate the limitations of DFT for interpreting XPS spectra and provide a benchmark for improving computational methods.

Adsorption↗

DFT simulation of the X-ray diffraction pattern of aluminum-ion-intercalated graphite used as the cathode material of the aluminum-ion battery

Many breakthroughs have been achieved in rechargeable aluminum-ion battery technologies in recent years. Most recently, operando X-ray diffraction (XRD) combined with density functional theory (DFT) calculations was reported to study the chloroaluminate anion (AlCl 4 -)-intercalated graphite cathode of the battery. However, there are quite a few discrepancies between the measured and simulated XRD patterns. This work is focused on the simulation of XRD patterns of graphite intercalation compounds (GICs) with DFT calculations. Our results reveal that both the geometry of AlCl 4 - in graphite and the gallery height of GICs are dependent on the intercalant density. At low intercalant density, the gallery height keeps constant, but at high intercalant densities, the gallery height is linearly related to the intercalant density. Our simulated XRD patterns are highly consistent with the measured operando XRD patterns. Not only do the angles of the peaks match very well, but also the relative intensities and the corresponding electrode capacities show reasonable agreement with the experimental results. The DFT simulation of the XRD pattern provides significant information on the stage index and the charge capacity of the GIC electrode.

25 ENERGY STORAGE↗

Evolution of lithium ordering with (de)-lithiation in β-LiVOPO 4 : insights through solid-state NMR and first principles DFT calculations

The lithium-ion battery cathode material β-VOPO 4 is capable of intercalating more than one Li ion per transition metal ion due to the accessibility of both the V 5+ /V 4+ and V 4+ /V 3+ redox couples at ~4.5 V and ~2.3 V vs. Li, respectively, giving a theoretical capacity greater than ~300 mA h g -1 . The ability to perform full and reversible two Li-ion intercalation in this material, however, has been a matter of debate and the poor crystallinity of the fully lithiated phase has thus far precluded its complete structural characterisation by conventional diffraction-based methods. In this work, 7 Li and 31 P NMR spectroscopy, in combination with first principles DFT calculations, indicate that chemical lithiation results in a single phase β-Li 2 VOPO 4 exhibiting a complex Li ordering scheme with lithium ions occupying multiple disordered environments. 2D NMR 7 Li correlation experiments were used to deduce the most likely Li ordering for the β-Li 2 VOPO 4 phase from amongst several DFT optimised structures. In contrast, electrochemically lithiated β-Li 2- x VOPO 4 discharged to 1.6 V exhibits, in addition to β-Li 2 VOPO 4 , a β-Li 1.5 VOPO 4 phase. Additionally, the existence of β-Li 1.5 VOPO 4 is not reflected in the flat galvanostatic charge and discharge curves nor is evident from diffraction-based methods due to the very close structural similarity between the β-Li 1.5 VOPO 4 phase and β-Li 2 VOPO 4 phases. We demonstrate that solid state NMR spectroscopy, in combination with DFT results, provides a powerful tool for identifying intermediate states formed during charge/discharge of these complex phosphates as these phases can be distinguished from the end member phases primarily by the nature of the lithium ordering.

25 ENERGY STORAGE↗

DFT-based QM/MM with Particle-Mesh Ewald for Direct, Long-Range Electrostatic Embedding

In this work, we present a DFT-based, QM/MM implementation with long-range electrostatic embedding achieved by direct real-space integration of the particle mesh Ewald (PME) computed electrostatic potential. The key transformation is the interpolation of the electrostatic potential from the PME grid to the DFT quadrature grid, from which integrals are easily evaluated utilizing standard DFT machinery. We provide benchmarks of the numerical accuracy with choice of grid size and real-space corrections, and demonstrate that good convergence is achieved while introducing nominal computational overhead. Furthermore, the approach requires only small modification to existing software packages, as is demonstrated with our implementation in the OpenMM and Psi 4 software. After presenting convergence benchmarks, we evaluate the importance of long-range electrostatic embedding in three solute/solvent systems modeled with QM/MM. Water and BMIM/BF 4 ionic liquid were considered as "simple" and "complex" solvents respectively, with water and p-phenylenediamine (PPD) solute molecules treated at QM level of theory. While electrostatic embedding with standard real-space truncation may introduce negligible error for simple systems such as water solute in water solvent, errors become more significant when QM/MM is applied to complex solvents such as ionic liquids. An extreme example is the electrostatic embedding energy for oxidized PPD in BMIM/BF 4 for which real-space truncation produces severe error even at 2-3 nm cutoff distances. This latter example illustrates that utilization of QM/MM to compute redox potentials within concentrated electrolytes/ionic media requires carefully chosen long-range electrostatic embedding algorithms, with our presented algorithm providing a general and robust approach.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Assessment of DFT functionals for a minimal nitrogenase [Fe(SH) 4 H]- model employing state-of-the-art ab initio methods

We have designed a [Fe(SH)4H]- model with the fifth proton binding either to Fe or S. We show that the energy difference between these two isomers (ΔE) is hard to estimate with quantum-mechanical (QM) methods. For example, different density functional theory (DFT) methods give ΔE estimates that vary by almost 140 kJ/mol, mainly depending on the amount of exact Hartree–Fock included (0%–54%). The model is so small that it can be treated by many high-level QM methods, including coupled-cluster (CC) and multiconfigurational perturbation theory approaches. With extrapolated CC series (up to fully connected coupled-cluster calculations with singles, doubles, and triples) and semistochastic heat-bath configuration interaction methods, we obtain results that seem to be converged to full configuration interaction results within 5 kJ/mol. Our best result for ΔE is 101 kJ/mol. With this reference, we show that M06 and B3LYP-D3 give the best results among 35 DFT methods tested for this system. Brueckner doubles coupled cluster with perturbaitve triples seems to be the most accurate coupled-cluster approach with approximate triples. CCSD(T) with Kohn–Sham orbitals gives results within 4–11 kJ/mol of the extrapolated CC results, depending on the DFT method. Single-reference CC calculations seem to be reasonably accurate (giving an error of ~5 kJ/mol compared to multireference methods), even if the D1 diagnostic is quite high (0.25) for one of the two isomers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Accuracy of DFT computed oxygen-vacancy formation energies and high-throughput search of solar thermochemical water-splitting compounds

The enthalpy change involved in metal oxide reduction is a key quantity in various processes related to energy conversion and storage, and is of particular interest for computational prediction. Often this prediction involves the simulation of a high temperature reduction process with a 0K methodology like density functional theory (DFT), and it is not infrequent for the high temperature and 0K stable crystal structures to differ. This introduces a conundrum with regards to the choice of crystal structure to utilize in the computation, with approaches in the literature varying and experimental validation remaining scarce. In this work we address both the crystal structure conundrum and the experimental validation, and then apply the insights we gain to guide a high-throughput search for new materials for solar thermochemical water-splitting applications. By computing the DFT+U oxygen vacancy formation energy (ΔE vf ) of a selection of ABO 3 compounds and comparing different crystal structures for each composition, we highlight the issues that arise when the structure utilized in the computation is dynamically unstable at 0K, namely the presence of an artificial lowering of ΔE vf , and the lack of convergence of ΔE vf with cell size. We solve these limitations by identifying and employing a suitable surrogate dynamically stable structure. We then validate the predictive power of our calculations against appositely generated experimental measurements of reduction enthalpy for a series of Hubbard U values, finding an accuracy ranging between 0.2-0.6 eV/O. In light of such conclusions, we revise and expand a previous a high-throughput DFT study on ABO 3 perovskite oxides. As a result, we provide a list of candidate STCH materials, highlight trends with redox-active cation and structural distortion, and identify Mn 4+ , Mn 3+ and Co 3+ as the most promising redox-active cations.

08 HYDROGEN↗

Strong Correlation DMRG and DFT

This project developed new ways to improve computer simulations of materials where electrons interact strongly with each other, a challenge for today’s most widely used method, density functional theory (DFT). We used an exact numerical method, the density matrix renormalization group (DMRG), to create highly accurate reference results for simple model systems, and used these to test DFT, prove when it will converge, and even train machine-learned functionals. We also invented new kinds of localized basis functions (“gausslets” and “multi-sliced gausslets”) and a “sliced-basis” approach that make high-accuracy simulations faster and more practical. These methods were applied to extended hydrogen systems, enabling the direct derivation of accurate low-energy models from first-principles calculations. We also introduced a new formalism, Conditional-Probability DFT, which could bypass traditional approximations. The tools and results from this work, including open-source software releases, will help scientists design and understand complex quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

PFOA catalytic reaction mechanisms on zerovalent iron (Fe 0 ): A DFT investigation

In this study, density functional theory (DFT) calculations with dispersion corrections are employed to investigate catalytic reactions of perfluorooctanoic acid (PFOA) onto zerovalent iron (Fe 0 ). The main goal of this investigation is to explain the reaction mechanisms, including bond dissociation energies, activation barriers, and rate-determining steps (RDSs), of PFOA degradation. Additionally, we aim to reveal the catalytic effects of Fe 0 by comparing the reaction energy profile for the degradation of PFOA both in solution (i.e., isolated) and after adsorption on Fe 0 . Along with investigating previously reported pathways, this study proposes a new PFOA degradation pathway on Fe0. Overall, the results revealed that Fe 0 plays a key role in the catalytic degradation of PFOA, significantly affecting the RDSs and increasing defluorination rates. For instance, the activation barrier for the decarboxylation step significantly decreased from 3.08 eV to 1.11 eV over the Fe0 surface, while the activation barrier for F dissociation decreased from ∼5 eV to 0.54 eV. Additionally, the proposed pathway in this study (i.e., successive defluorination) was found to be the most favorable pathway for PFOA degradation on Fe 0 . Taken together, the DFT results are consistent with reported experimental findings and provide critical insights that will help advance the design of catalysts for PFAS degradation.

Catalysis↗

DFT Mechanism Studies: Biomimetic 1,4-NADH Chemoselective, Co-factor Regeneration with [Cp*Rh(bpy)H] + , in Tandem with the Biocatalysis Pathways of a Core Model of the (HLADH)-Zn(II) Mediated Enzyme, in the Enantioselective Reduction of Achiral Ketones to Chiral S-Alcohols

In this study, Quantum Chemical (QC) calculations, utilizing Density Functional Theory (DFT), were performed to investigate the mechanistic aspects of the chemoselective catalyzed reaction of [Cp*Rh(bpy)H] + with the biomimetic NAD + analogues, N-benzylnicotinamide triflate, 1, and β-nicotinamide ribose-5'-methyl phosphate, 2, in the conversion to their 1,4-NADH analogues, 1,4-dihydro-N-benzylnicotinamide, 4, and β-1,4-dihydronicotinamide-5'-ribose methyl phosphate, 5. This reaction was in tandem with the 1,4-NADH dependent HLADH-Zn(II)- catalyzed reduction of achiral ketones to chiral S-alcohols. The [Cp*Rh(bpy)H] + complex, and not its equilibrium tautomer, [η 4 -Cp*HRh(bpy)] + , was found to control the hydride transfer during the biomimetic NAD + /1,4-NADH conversion, through the non-covalent interactions of the biomimetic co-factors with [Cp*Rh(bpy)H] + . The thermodynamics and kinetics for the chiral reduction of the Zn(II) bound ketones, 2-pentanone and 4-phenyl-2-butanone, with co-factor, 4, catalyzed by Zn(SCH 3 ) 2 (Imidazole), a core model of the Zn(II)-based catalytic center of HLADH, was also investigated by the evaluation of two possible reaction pathways: (1) formation of a ZnH from the C4-H hydride transfer of co-factor, 4, followed by reaction of the postulated ZnH with the bound 2-pentanone or 4-phenyl-2-butanone substrate, and (2), the direct C4-H transfer to the bound achiral ketone substrates, to provide the dominant chiral alcohols, S-2-pentanol or S-4-phenyl-2-butanol. The latter pathway was found most viable, and DFT calculations also revealed an essential η 2 -coordination of the 5,6 double bond of co-factor, 4, to the HLADH-Zn(II) metal ion center, upon imidazole decomplexation, providing an asymmetric differentiation of S-η 2 -5,6-1,4-NADH-Zn(II) binding. A proposed new paradigm for the Zn(II)'s non-innocent role in the HLADH-Zn(II) biocatalysis reduction mechanism, for enantioselective hydride transfer to a Zn(II) bound ketone, providing S-alcohols.

1,4 NADH co-factors↗

Bayesian automated weighting of aggregated DFT, MD, and experimental data for candidate thermodynamic models of aluminum with uncertainty quantification

Atomic-scale modeling methods such as density functional theory (DFT) and molecular dynamics (MD) can predict the thermodynamic properties of materials at a lower cost than experimental measurements. However, their regular usage in thermodynamic model construction is hampered by the lack of quantitative agreement with experimental measurements and the lack of uncertainty estimates on the data. To make regular usage of this atomistic simulation data, it is important to assess whether the atomistic simulation datasets, by themselves or in combination with experimental measurements, result in the same physics-informed models best supported by experimental measurements alone. Here, models of aluminum thermodynamic properties are discussed using three data sources: atomistic calculations (DFT and MD), experiments, and a combination of atomistic calculations and experiments. The study shows that, after ensuring self-consistency in predicting key invariant points, both experimental measurements and atomistic calculations can significantly contribute to an optimal model.

36 MATERIALS SCIENCE↗

Ethylene oxidation on unpromoted silver catalysts: Reaction pathway and selectivity analysis using DFT calculations

Chemical conversions in catalytic partial oxidation processes of light hydrocarbons are responsible for the production of numerous industrial chemicals, plastics, and intermediates. These processes are relatively expensive to perform, and are typically operated at high thermodynamic inefficiency, so the development of novel, highly efficient catalysts would prove to be very cost effective. Herein, our study focused on surface catalytic mechanisms of the ethylene oxide (EO) formation process. Periodic plane-wave Density Functional Theory (DFT) methods were used to analyze related reaction mechanisms on the Ag(111) surface facet with low coverage. Energetic changes of related species and pathways were calculated. Key surface species are identified to suggest the factors for the observed selectivity during EO formation. Our results are consistent with previous kinetic modeling efforts in the literature which did not employ DFT analysis. Lastly, our study demonstrates how fundamental theoretical investigations and multi-scale modeling techniques are currently impacting the advancement of rational catalyst design and microkinetic modeling techniques in the light hydrocarbon processing industry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Surface Phase Stability of Fe 2 O 3 (001) in Hydrogen Reducing Environments: A DFT and XPS Analysis

Here, this study combines density functional theory (DFT) and ab initio thermodynamics calculations with X-ray photoelectron spectroscopy (XPS) investigations to identify the reduction properties of the Fe 2 O 3 (001) surface with implications for corrosion resistance, hydrogen transport, and energy safety. Ab initio thermodynamics modeling predicts fully hydroxylated surface stability across a broad range of pressures (1 × 10 –23 to 1 × 10 5 mbar) and temperatures below 700 K, consistent with previous experimental studies. Above 800 K, exposures to 1 × 10 –4 mbar H 2 , 1 × 10 –4 mbar O 2 , or 1 × 10 –4 mbar H 2 + 1 × 10 –4 mbar O 2 each yield unique XPS signals indicating a loss of −OH coverage, aligning with DFT predictions. Insight into the mechanism of reduction as a function of H 2 exposure is provided, as well as conditions that promote further reduction toward Fe 3 O 4 . Theoretical and experimental investigations indicate the ability to maintain the Fe 2 O 3 protective layer of iron oxides that have been exposed to H 2 environments by including trace amounts of aqueous O2.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Revealing the Formation Energy–Exfoliation Energy–Structure Correlation of MAB Phases Using Machine Learning and DFT

MAB phases are popular as high-temperature materials with high damage tolerance and excellent electrical conductivity that are used to exfoliate two-dimensional (2D) transition metal borides (MBenes); a promising material for developing next- generation nanodevices. In this report, we explore the correlation between the formation energy, exfoliation energy and structural factors of MAB phases using DFT and machine learning. In order to do so, we developed three different machine learning models based on the support vector machine, deep neural network, and random forest regressor to study the stability of the MAB phases by calculating their formation energies. Furthermore, our support vector machine and deep neural network models are capable of predicting the formation energies with mean absolute errors less than 0.1 eV/atom. We demonstrated that the stability of a MAB phase for a given transition metal decreases when the A element changes from Al to Tl. The density functional theory (DFT) re- vealed that M-A and M-B bond strength strongly correlates with the stability of MAB phases. In addition, the exfoliation possibility of 2D MBenes becomes higher when the A element changes from Al to Tl due to weakening of M-A and M-B bonds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electronic structure analysis and DFT benchmarking of Rydberg-type alkali-metal-crown ether, -cryptand, and -adamanzane complexes

Density functional theory (DFT) and electron propagator theory (EPT) calculations were performed to study ground and excited electronic structures of alkali-metal (M) coordinated 9-crown-3, 24-crown-8, [2.1.1]cryptand, o-Me 2 -1.1.1, and 3 6 Adamanzane complexes. Each complex bears an expanded electron in the periphery and occupies diffuse 1p-, 1d-, 1f-type molecular orbitals (or superatomic 1P, 1D, 1F orbitals) in excited electronic states. The calculated superatomic shell model of the M(9-crown-3) 2 is 1S, 1P, 1D, 1F, 2S, 2P, 2D, 1G and it is held by all other complexes up to the studied 1F level. Due to the highly diffuse nature of the electron, the ionization energies of these complexes are significantly lower (1.6–2.0 eV) and hence these complexes belong to the superalkali category. Further, the ab initio EPT ionization energy and the excitation energies of the Li(9-crown-3) 2 were used to evaluate DFT errors associated with a series of exchange correlation functionals that span multiple rungs of Jacob's ladder (i.e., GGA, meta-GGA, global GGA hybrid, meta-GGA hybrid, range-separated hybrid, double-hybrid). Among these, the best performing functional is the range-separated hybrid CAM-B3LYP and the errors are within 6% of high-level ab initio EPT results. The accuracy of CAM-B3LYP is indeed transferable to similar complexes and hence the findings are expected to accelerate the progression of studies of Rydberg-type systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗