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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Extending Orbital-Optimized Density Functional Theory to L-Edge XPS and Beyond: Spin–Orbit Coupling via Nonorthogonal Quasi-Degenerate Perturbation Theory

Quantum mechanical calculations of core electron binding energies (CEBEs) are relevant to interpreting X-ray photoelectron spectroscopy (XPS). Orbital-optimized density functional theory (OO-DFT) accurately predicts K-edge CEBEs but is challenged by the presence of significant spin–orbit coupling (SOC) at L- and higher edges involving inner-shell orbitals with nonzero angular momentum. Here, to extend OO-DFT to L-edges and higher, our method utilizes scalar-relativistic, spin-restricted open-shell OO-DFT to construct a minimal, quasi-degenerate basis of core-hole states corresponding to a chosen inner-shell (e.g., ionizing all six possible 2p spin orbitals). Nonorthogonal configuration interaction (NOCI) is then used to obtain the matrix elements of the full Hamiltonian including SOC in this quasi-degenerate model space of determinants. Using a screened 1-electron SOC operator parametrized with the Dirac-Coulomb-Breit (DCB) Hamiltonian results in doublet splitting (DS) values for third row elements that are nearly in quantitative agreement with experiment. The resulting NOCI eigenvalues are shifted by the average of the (scalar) OO-DFT CEBEs to yield CEBEs (split by SOC) corrected for dynamic correlation. Comparing calculations on gas phase molecules with experimental results establishes that NO-QDPT with the SCAN functional (NO-QDPT/SCAN), using the DCB screened 1-electron SOC operator is accurate to about 0.2 eV for L-edge CEBEs of molecules containing third row atoms. However, this NO-QDPT approach becomes less accurate for fourth-row elements starting in the middle of the 3d transition metal series, with errors increasing as atomic number increases.

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

Many-Body Benchmark of Electronic Charge and Spin Densities for Li 1–x NiO 2

Accurate benchmarks are particularly important for highly correlated oxides as mean-field approximations often fail to describe the subtle balance of charge transfer and magnetism in these materials with an accuracy comparable to experimental needs. Here we present accurate diffusion Monte Carlo (DMC) results of the electronic charge and spin densities for the tunable highly correlated oxide Li 1–x NiO 2 for x = 0, 1/2, and 1. To enable quantitative comparisons, we introduce a robust density-partitioning scheme, extending Voronoi analysis to assign atomic charges from spatially noisy DMC densities. We then benchmark common approximations used in density functional theory (DFT). Comparison against DMC shows that r 2 SCAN delivers the most balanced performance across charge, spin, and radial density descriptors, nearly reproducing DMC results for LiNiO 2 and apical Ni sites in Li 0.5 NiO 2 . Hybrid functionals (PBE0, SCAN0) perform unexpectedly poorly, and PBE + U + V yields inconsistent trends between charge and spin densities. Therefore, the r 2 SCAN functional minimizes errors relative to DMC while capturing the variable valence of the Ni ion and also retaining the computational efficiency of DFT for large-scale simulations of the tunable structural and electronic phases of Li1−xNiO2. Our study highlights the importance of accurate benchmarking of the fundamental quantities involved in DFT to select appropriate DFT approximations in order to advance the predictive modeling of charge-transfer-driven phenomena in correlated electron systems.

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

Projection-Based Correlated Wave Function in Density Functional Theory Embedding for Periodic Systems

In this report we present a level shift projection operator-based embedding method for systems with periodic boundary conditions—where the “active” subsystem can be described using either density functional theory (DFT) or correlated wave function (WF) methods and the “environment” is described using DFT. Our method allows for k-point sampling, is shown to be exactly equal to the canonical DFT solution of the full system under the limit that we use the full system basis to describe each subsystem, and can treat the active subsystem either with periodic boundary conditions—in what we term “periodic-in-periodic” embedding—or as a molecular cluster—in “cluster-in-periodic” embedding. We explore each of these methods and show that cluster WF-in-periodic DFT embedding can accurately calculate the absorption energy of CO on to a Si(100)-2×1 surface.

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

Application of the ChIMES Force Field to Nonreactive Molecular Systems: Water at Ambient Conditions

We demonstrate development of the Chebyshev Interaction Model for Efficient Simulation (ChIMES) for molecular systems through application to water under ambient conditions (298 K, 1 g/cm 3 ). These models, which are comprised of linear combinations of Chebyshev polynomials explicitly describing two- and three-body interactions, are largely fit by force matching to Kohn–Sham Density Functional Theory (DFT). Protocols for selecting user-specified parameters and inclusion of stress tensor data are investigated, and structural and dynamical property prediction for resulting models is benchmarked against DFT. We show that the present ChIMES force fields yield excellent agreement with DFT without the need for additional terms such as those for Coulomb interactions. Overall, we show that tractable parametrization and subsequent accuracy of the present models make ChIMES an ideal candidate for extension of DFT dynamics to larger system sizes and longer time scales.

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

Density-Functional Tight-Binding Parameters for Bulk Zirconium: A Case Study for Repulsive Potentials

Density-functional tight-binding (DFTB) parameters are presented for the simulation of the bulk phases of zirconium. Electronic parameters were obtained using a band structure fitting strategy, while two-center repulsive potentials were created by particle swarm optimization. As objective functions for the repulsive potential fitting, we employed the Birch–Murnaghan equations of state for hexagonal close-packed (HCP), body-centered cubic (BCC) and ω phases of Zr from density-functional theory (DFT). When fractional atomic coordinates are not allowed to change in the generation of the equation-of-state curves, long-range repulsive DFTB potentials are able to almost perfectly reproduce equilibrium structures, relative DFT energies of the bulk phases, and bulk moduli. However, the same potentials lead to artifacts in the DFTB potential energy surfaces when atom positions in the unit cell are allowed to fully relax during the change of unit cell parameters. Finally, conventional short-range repulsive DFTB potentials, while inferior in their ability to reproduce DFT bulk energetics, are able to correctly reproduce the qualitative shape of the DFT potential energy surfaces, including the location of global minima, and can therefore be considered more transferable.

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

Metal-to-Ligand Charge-Transfer Spectrum of a Ru-Bipyridine-Sensitized TiO 2 Cluster from Embedded Multiconfigurational Excited-State Theory

Understanding optical properties of the dye molecule in dye-sensitized solar cells (DSSCs) from first-principles quantum mechanics can contribute to improving the efficiency of such devices. While density functional theory (DFT) and time-dependent DFT have been pivotal in simulating optoelectronic properties of photoanodes used in DSSCs at the atomic scale, questions remain regarding DFT’s adequacy and accuracy to furnish critical information needed to understand the various excited-state processes involved. In this work, we simulate the absorption spectra of a dye-sensitized solar cell analogue, comprised of a Ru-bipyridine (Ru-bpy) dye molecule and a small TiO2 cluster via DFT and via an accurate embedded correlated wavefunction (CW) theory. We generated CW spectra for the adsorbed Ru-bpy dye via a recently introduced capped density functional embedding theory or capped-DFET (to generate the embedding potential that accounts for the interaction of the molecule and the TiO 2 cluster). We then combined capped-DFET with the accurate but expensive multiconfigurational complete active space second-order perturbation theory (CASPT2)—embedded CASPT2. Because the CW theory is conducted on only a portion of the total system in the presence of an embedding potential that describes that portion’s interaction with its environment, we efficiently obtain CW-quality predictions that reflect local properties of the entire system. Specifically, for example, with capped-DFET and embedded CW theory, we can simulate accurately a plethora of metal-to-ligand charge-transfer excited properties at a manageable computational cost. Here, we predict detailed electronic spectra within the visible region, featuring the lowest three singlet and triplet excited states, along with predictions of the singlets’ lifetimes. We illustrated these results using a Jablonski diagram that show the relative energy position of the singlet and longer-lived triplet excited states and analyzed and proposed relaxation paths for the excited state corresponding to the most intense but short-lived absorption (interconversion, intersystem crossing, fluorescence, and phosphorescence) that may lead to longer-lived excited states necessary for efficient charge separation required to generate current in solar cells.

14 SOLAR ENERGY↗
DOE OSTI2022

DeePKS + ABACUS as a Bridge between Expensive Quantum Mechanical Models and Machine Learning Potentials

Recently, the development of machine learning (ML) potentials has made it possible to perform large-scale and long-time molecular simulations with the accuracy of quantum mechanical (QM) models. However, for different levels of QM methods, such as density functional theory (DFT) at the meta-GGA level and/or with exact exchange, quantum Monte Carlo, etc., generating a sufficient amount of data for training an ML potential has remained computationally challenging due to their high cost. In this work, we demonstrate that this issue can be largely alleviated with Deep Kohn–Sham (DeePKS), an ML-based DFT model. DeePKS employs a computationally efficient neural network-based functional model to construct a correction term added upon a cheap DFT model. Upon training, DeePKS offers closely matched energies and forces compared with high-level QM method, but the number of training data required is orders of magnitude less than that required for training a reliable ML potential. As such, DeePKS can serve as a bridge between expensive QM models and ML potentials: one can generate a decent amount of high-accuracy QM data to train a DeePKS model and then use the DeePKS model to label a much larger amount of configurations to train an ML potential. Further, this scheme for periodic systems is implemented in a DFT package ABACUS, which is open source and ready for use in various applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Modeling Pb(II) Adsorption on Mineral Surfaces: Bridging Density Functional Theory and Experiment with Thermodynamic Insights

Despite decades of work on aqueous lead (Pb) adsorption on a-Fe2O3 (hematite) and a-Al2O3 (alumina), gaps between measurements and modeling obscure molecular-level understanding. Achieving well-matched geometries between theory and experimental for mineral-water interfaces is a hurdle, as surface functional group type and distribution must be accounted for in determining mechanisms. Additionally, computational methods that can describe the substrate are often not appropriate to capture aqueous effects. Progress requires focusing on well-studied and relevant systems, such as key facets (001),(012), and (110) of hematite and alumina, and ubiquitous contaminants such as aqueous Pb. In the past, bulk-parametrized bond-valence principles were used to rationalize Pb(II) adsorption trends. These approaches can break down at surfaces, where flexible bonding environments and adsorption-induced surface relaxations play a critical role. Here, we adapt and apply a density functional theory (DFT) and thermodynamics framework, integrating DFT-calculated energies with experimental data and electrochemical principles, to predict Pb(II) adsorption. Our model results capture trends across the full set of surfaces and predict that inner-sphere Pb(II) sorption on (001) alumina varies from unfavorable to weakly favorable across a range of pH conditions. This aligns with experiment insights that Pb(II) interacts at that surface through outer-sphere interactions. Extending to Fe(II) adsorption, we demonstrate a coverage-dependent site preference, potentially explaining disorder in overlayers grown by the oxidative adsorption of Fe(II) on hematite (001).

lead contamination↗
DOE OSTI2020

Binding of Formic Acid on Anatase TiO2(101)

The adsorption of formic acid (FA) and the formation of formate species on metal oxide surfaces are of great interest in catalysis. Formic acid is used to probe the adsorption properties of surface sites and formates are common intermediates in many catalytic reactions. Here we focus on the interaction of FA with a prototypical anatase TiO2(101) surface by using a combination of scanning tunneling microscopy (STM), infrared reflection absorption spectroscopy (IRAS), electron stimulated desorption (ESD), and density functional theory (DFT) to assess the coverage dependent evolution of different FA-derived surface species at low temperatures (80–240 K). We find that isolated FA adsorbs at 80 K in both monodentate (MD) and bidentate (BD) configurations on top of undercoordinated Ti sites (Ti5c). The MD form is likely deprotonated and readily converts upon annealing to bridging BD on two neighboring Ti5c sites. DFT calculations show that molecularly-bound MD FA is metastable and readily converts into a deprotonated state in the proximity of subsurface oxygen vacancy. The stability of the MD species increases as the availability of the paired Ti5c sites for BD formation diminishes at high coverages. Upon surface saturation, a mixed configuration of alternating bidentate and monodentate species with coverage of 2/3 FA/Ti5c represents the most favorable configuration. This is in contrast with the adsorption of FA on rutile TiO2(110) and many other oxides, where bidentate formate species dominate. Y.W., A.D., G.A.K., R.R., N.G.P., and Z.D. were supported by the U.S. Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences (BES), Division of Chemical Sciences, Geosciences and Biosciences (CSGB) under FWP 47319. B.W. and A.S. acknowledge the support of DOE BES, CSGB Division under Award DESC0007347. The experimental studies were performed in EMSL, a national scientific user facility sponsored by the Department of Energy’s Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory (PNNL). PNNL is a multiprogram national laboratory operated for DOE by Battelle. The authors also acknowledge computational resources from the TIGRESS high-performance computer center at Princeton University.

Wang, Yang↗
DOE OSTI2020

Energetic Stability of Freestanding and Metal-Supported Borophenes: Quantum Monte Carlo and Density-Functional Theory Calculations

Quantum Monte Carlo (QMC) and density-functional theory (DFT) calculations were carried out to study cohesion energetics of two-dimensional (2D) sheets of boron atoms called borophenes. Our QMC calculations confirmed the polymorphism among free-standing borophenes that was reported in previous DFT studies. Although Perdew-Burke-Ernzerhof (PBE) calculations significantly overestimate the cohesive energies of 2D boron sheets, DFT-PBE relative energetics with respect to each other among various free-standing borophenes are found to be in quantitative agreement with the corresponding QMC results. This suggests that one can make reliable predictions for relative stability of different boron sheets through DFT-PBE calculations. Our analysis of PBE formation energies of borophenes on metal surfaces shows that the polymorphic range is extended for borophenes on the Ag(111) and the Au(111) surfaces beyond that of free-standing borophenes, reflecting recent experimental synthesis of beta(12) and chi(3) boron sheets on the Ag(111) surface. We have also found that a hexagonal borophene can be stabilized through charge transfer from a metal surface and is energetically favored on the AI(111) surface over other borophene structures. Finally, it is found that the bilayer formation could be energetically favored over its monolayer form for the borophene-Au system, especially for borophenes with hexagonal hole densities eta lower than 1/9. This leads to our prediction that in addition to its monolayer form, a bilayer eta = 1/12 borophene can be synthesized on the Au(111) surface, opening a new possibility for borophene-based electronic devices.

ab initio↗
DOE OSTI2021

Experimental and Theoretical Study of the Electronic Structures of Lanthanide Indium Perovskites LnInO 3

Ternary lanthanide indium oxides LnInO 3 (Ln = La, Pr, Nd, Sm) were synthesized by high-temperature solid-state reaction and characterized by Xray powder diffraction. Rietveld refinement of the powder patterns showed the LnInO 3 materials to be orthorhombic perovskites belonging to the space group Pnma, based on almost-regular InO 6 octahedra and highly distorted LnO 12 polyhedra. Experimental structural data were compared with results from density functional theory (DFT) calculations employing a hybrid Hamiltonian. Valence region X-ray photoelectron and K-shell X-ray emission and absorption spectra of the LnInO 3 compounds were simulated with the aid of the DFT calculations. Photoionization of lanthanide 4f orbitals gives rise to a complex final-state multiplet structure in the valence region for the 4f n compounds PrInO 3 , NdInO 3 , and SmInO 3 , and the overall photoemission spectral profiles were shown to be a superposition of final-state 4f n-1 terms onto the cross-section weighted partial densities of states from the other orbitals. The occupied 4f states are stabilized in moving across the series Pr-Nd-Sm. Band gaps were measured using diffuse reflectance spectroscopy. These results demonstrated that the band gap of LaInO 3 is 4.32 eV, in agreement with DFT calculations. This is significantly larger than a band gap of 2.2 eV first proposed in 1967 and based on the idea that In 4d states lie above the top of the O 2p valence band. However, both DFT and X-ray spectroscopy show that In 4d is a shallow core level located well below the bottom of the valence band. Band gaps greater than 4 eV were observed for NdInO 3 and SmInO 3 , but a lower gap of 3.6 eV for PrInO 3 was shown to arise from the occupied Pr 4f states lying above the main O 2p valence band.

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

Role of Dynamic Polarization Interactions in the Electrical Double Layer at Calcite (104) Interfaces with Aqueous Solutions

Reactions at mineral interfaces with aqueous solutions control many geochemical and biogeochemical processes in the Earth’s critical zone. At the molecular level, insights into important properties such as the structure of the electrical double layer (EDL) at specific mineral interfaces continue to improve due to the increasing fidelity of laboratory instrumentation and computational approaches. However, molecular simulation approaches suffer from limited reach into relevant scales of time, length, and system complexity. To span this gap, a novel hybrid approach that couples first principle plane-wave density functional theory (DFT) with classical DFT (cDFT) is demonstrated and applied to calcite (104) interfaces with various electrolytes. In this approach, a region of interest described using DFT interacts with the surrounding medium described using cDFT to arrive at a self-consistent ground state. Benchmarking against experimental observations and entirely first principle DFT simulations demonstrates that this hybrid model efficiently encompasses the key short-range and collective interactions in the EDL. Simulations of calcite (104)/solution interfaces reveal the key static and dynamic polarization interactions that give rise to structuring of ions and water. Ion hydration interactions have the strongest effect on the depth of the first minimum in the density distribution of counterions at the surface, and the position and width of the first density peak is largely determined by the strength of ion-correlation forces. Finer details of ion distributions are controlled by mutual polarization of the calcite surface and interfacial electrolyte. Finally, this new ability to efficiently and rigorously predict EDL structure at mineral surfaces in contact with complex solutions paves the way to accurately modeling sorption, nucleation, dissolution, and growth in realistic systems.

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

An Adsorption Study of CH 4 on ZSM-5, MOR, and ZSM-12 Zeolites

CH 4 adsorption was studied experimentally and theoretically on ZSM-5, MOR, and ZSM-12 zeolites using calorimetric measurements at 195 K and plane wave DFT calculations. Differential heats measured on four different H-ZSM-5 samples were determined to be 22.5 ± 1 kJ/mol, independent of Brønsted site density or defect concentration. However, DFT calculations performed using various functionals and on the most stable Brønsted site indicated that CH4 should bind to this site by an additional 1–7 kJ/mol, a discrepancy that is due to the inability of standard DFT methods to capture hydrogen-bonding effects accurately with CH 4 . Differential heats for CH 4 in MOR were 30±1 kJ/mol at low coverages, falling to 25 kJ/mol for coverages above one molecule per 8-membered-ring side pocket, while differential heats on ZSM-12 were initially 22.5 kJ/mol, decreasing to 21 kJ/mol with coverage. DFT calculations on the siliceous form of the zeolites were able to predict these values within 5 kJ/mol in most cases. The results indicate that CH 4 is an excellent probe molecule for characterizing the pore structure of zeolites.

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

Comparing Classical and Machine Learning Force Fields for Modeling Deformation of Metal–Organic Frameworks Relevant for Direct Air Capture

Deformation of metal–organic frameworks (MOFs) induced by adsorbate molecules can affect adsorption properties such as capacity and selectivity, but most computational studies of MOFs assume framework rigidity to simplify calculations. Although flexible force fields (FFs) for MOFs have been parametrized for specific materials, the generality of FFs for reliably modeling adsorbate-induced deformation to accuracy nearing that of density functional theory (DFT) has not been established. This work confirms using DFT calculations that adsorbate-induced deformation can affect CO 2 and H 2 O adsorption energies in a considerable fraction of MOFs promising for direct air capture (DAC). We then benchmark the efficacy of several general-purpose FFs in describing adsorbate-induced deformation for DAC against DFT. Our results show that current classical FFs are insufficient for describing MOF deformation, especially in cases of interest for DAC where strong interactions exist between adsorbed molecules and MOF frameworks. Some emerging machine learning force fields (MLFFs) we tested, particularly CHGNet, MACE-MP-0, and Equiformer V2, appear to be more promising than the classical FF for emulating the deformation behavior described by DFT. The best performing FF (CHGNet), however, fails to achieve the accuracy required for practical predictions with a mean absolute adsorption energy error of 0.124 eV.

adsorption↗
DOE OSTI2026

Examination of Replicate Syntheses of Metal Organic Frameworks as a Window into Reproducibility in Materials Chemistry

Replicate experiments are a useful tool in understanding the repeatability of scientific measurements. In 2019, a systematic search for replicate syntheses of a collection of 130 metal–organic frameworks (MOFs) found that 89% of these materials had no reported replicate syntheses apart from the original publications identifying the material (Agrawal, M. Proc. Natl. Acad. Sci. U.S.A. 2020, 117, 877−88210.1073/pnas.1918484117). A potential weakness of that search was that only 5–11 years had elapsed since the original publication of each material. Here, this analysis is extended to all publications 11–17 years after the original publication. Although this extended time period identifies more repeat syntheses, 83% of the materials still have no reported replicate syntheses. We also consider how appropriately selected Density Functional Theory (DFT) calculations can provide corroboration for the experimentally reported crystal structures. By using data from previous high-throughput DFT studies, corroborating evidence from DFT was available for 17% of the 130 structures for which no replicate syntheses are available. In total, approximately 1/3 of the 130 MOFs have data associated with replicate synthesis experiments and/or directly corroborating DFT calculations.

Sholl, David S. [Oak Ridge National Laboratory (OR↗
DOE OSTI2016

Toward Accurate Adsorption Energetics on Clay Surfaces

Clay minerals are ubiquitous in nature, and the manner in which they interact with their surroundings has important industrial and environmental implications. Consequently, a molecular-level understanding of the adsorption of molecules on clay surfaces is crucial. In this regard computer simulations play an important role, yet the accuracy of widely used empirical force fields (FF) and density functional theory (DFT) exchange-correlation functionals is often unclear in adsorption systems dominated by weak interactions. Herein we present results from quantum Monte Carlo (QMC) for water and methanol adsorption on the prototypical clay kaolinite. To the best of our knowledge, this is the first time QMC has been used to investigate adsorption at a complex, natural surface such as a clay. As well as being valuable in their own right, the QMC benchmarks obtained provide reference data against which the performance of cheaper DFT methods can be tested. Indeed using various DFT exchange-correlation functionals yields a very broad range of adsorption energies, and it is unclear a priori which evaluation is better. QMC reveals that in the systems considered here it is essential to account for van der Waals (vdW) dispersion forces since this alters both the absolute and relative adsorption energies of water and methanol. We show, via FF simulations, that incorrect relative energies can lead to significant changes in the interfacial densities of water and methanol solutions at the kaolinite interface. Despite the clear improvements offered by the vdW-corrected and the vdW-inclusive functionals, absolute adsorption energies are often overestimated, suggesting that the treatment of vdW forces in DFT is not yet a solved problem.

Chemistry↗
DOE OSTI2017

Coupled Cluster and Density Functional Studies of Atomic Fluorine Chemisorption on Coronene as Model Systems for Graphene Fluorination

Potential energy curves (PECs) of atomic fluorine adsorption on coronene as a model for graphene or nanocarbon surfaces have been computed. The PECs were obtained by scanning the fluorine atom distance to one of the center carbon atoms of the coronene molecule as model system in a “top” position from 4.0 to1.0 Å in intervals of 0.1 Å using a variety of quantum chemical methods. Various density functional theory (DFT) functionals, such as B3LYP, PBE, PBE0, CAM-B3LYP, and LC-ωPBE; approximate DFT methods such as several levels of the density-functional tight-binding (DFTB) method, as well as ab initio wave function theory methods such as MP2, CCSD, CCSD(T), and G2MS extrapolations, were used to evaluate energies for B3LYP/cc-pVDZ PEC geometries. G2MS is an approximation to the highly accurate CCSD(T)/cc-pVTZ level of theory in our work. We found that fluorine is chemically adsorbed on coronene with a binding energy of 22.9 and 23.3 kcal/mol at the B3LYP/cc-pVDZ and G2MS levels of theory, respectively, and 18.3 and 19.3 kcal/mol after counterpoise correction. Here, we found that pure DFT functionals and their DFTB approximations fail to predict the correct dissociation limit due to the DFT-inherent self-interaction error and various limitations in the DFTB approximation itself.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗