SEARCH · Engineering Papers

Results for “Quantum chemistry methods”

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.

SourceResearch connectionFrom yearThrough year

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

At least 883 records · Page 49

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Importance of Electron Correlation on the Geometry and Electronic Structure of [2Fe–2S] Systems: A Benchmark Study of the [Fe 2 S 2 (SCH 3 ) 4 ] 2–,3–,4– , [Fe 2 S 2 (SCys) 4 ] 2– , [Fe 2 S 2 (S- p -tol) 4 ] 2– , and [Fe 2 S 2 (S- o -xyl) 4 ] 2– Complexes

Iron–sulfur clusters are crucial for biological electron transport and catalysis. Obtaining accurate geometries, energetics, manifolds of their excited electronic states, and reduction energies is important to understand their role in these processes. Here, using a [2Fe–2S] model complex with Fe II and Fe III oxidation states, which leads to different charges, i.e., [Fe 2 S 2 (SMe) 4 ] 2–,3–,4– , we benchmarked a variety of computational methodologies ranging from density functional theory (DFT) to post-Hartree–Fock methods, including complete active space self-consistent field (CASSCF), multireference configuration interaction, the second-order N-electron valence state perturbation theory (NEVPT2), and the linearized integrand approximation of adiabatic connection (AC0) approaches. Additionally, we studied three experimentally well-characterized complexes, [Fe 2 S 2 (SCys) 4 ] 2– , [Fe 2 S 2 (S-o-tol) 4 ] 2– , and [Fe 2 S 2 (S-o-xyl) 4 ] 2– , via DFT methods. We conclude that the dynamic electron correlation is important for accurately predicting the geometry of these complexes. Broken symmetry (BS) DFT correctly predicts experimental geometries of low-spin multiplicity, while CASSCF does not. However, BS-DFT significantly overestimates the difference between the low- and high-spin electronic states for a given oxidation state. At the same time, CASSCF underestimates it but provides relative energies closer to the reference NEVPT2 results. Finally, AC0 provides energetics of NEVPT2 quality with the additional advantage of being able to use large CASSCF sizes. NEVPT2 gives the best estimates of the Fe III /Fe III → Fe II /Fe III (4.27 eV) and Fe II /F III → Fe II /F II (7.72 eV) reduction energies. The results provide insight into the electronic structure of these complexes and assist in the understanding of their physical properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Recent Developments in DFTB+, a Software Package for Efficient Atomistic Quantum Mechanical Simulations

DFTB+ is a flexible, open-source software package developed by its community, designed for fast and efficient atomistic quantum mechanical simulations. It employs various methods that approximate density functional theory (DFT), such as density functional-based tight binding (DFTB) and the extended tight binding (xTB) approach allowing simulations of large systems over extended time scales with reasonable accuracy, while being significantly faster than traditional ab initio methods. In recent years, several new extensions of the DFTB method have been developed and implemented in the DFTB+ program package in order to improve the accuracy and generality of the available simulation results. In this paper, we review those enhancements, show several use case examples and discuss the strengths and limitations of its features.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Understanding Reaction Networks through Controlled Approach to Equilibrium Experiments Using Transient Methods

We report a combined experimental/theoretical approach to studying heterogeneous gas/solid catalytic processes using low-pressure pulse response experiments achieving a controlled approach to equilibrium that combined with quantum mechanics (QM)-based computational analysis provides information needed to reconstruct the role of the different surface reaction steps. We demonstrate this approach using model catalysts for ammonia synthesis/decomposition. Polycrystalline iron and cobalt are studied via low-pressure TAP (temporal analysis of products) pulse response, with the results interpreted through reaction free energies calculated using QM on Fe-BCC(110), Fe-BCC(111), and Co-FCC(111) facets. In TAP experiments, simultaneous pulsing of ammonia and deuterium creates a condition where the participation of reactants and products can be distinguished in both forward and reverse reaction steps. This establishes a balance between competitive reactions for D* surface species that is used to observe the influence of steps leading to nitrogen formation as the nitrogen product remains far from equilibrium. Here, the approach to equilibrium is further controlled by introducing delay timing between NH 3 and D 2 which allows time for surface reactions to evolve before being driven in the reverse direction from the gas phase. The resulting isotopic product distributions for NH 2 D, NHD 2 , and HD at different temperatures and delay times and NH 3 /D 2 pulsing order reveal the role of the N 2 formation barrier in controlling the surface concentration of NH x * species, as well as providing information on the surface lifetimes of key reaction intermediates. Conclusions derived for monometallic materials are used to interpret experimental results on a more complex and active CoFe bimetallic catalyst.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Detection of Foreign Materials on Broiler Breast Meat Using a Fusion of Visible Near-Infrared and Short-Wave Infrared Hyperspectral Imaging

Foreign material (FM) found on a poultry product lowers the quality and safety of the product. We developed a fusion method combining two hyperspectral imaging (HSI) modalities in the visible-near infrared (VNIR) range of 400–1000 nm and the short-wave infrared (SWIR) range of 1000–2500 nm for the detection of FMs on the surface of fresh raw broiler breast fillets. Thirty different types of FMs that could be commonly found in poultry processing plants were used as samples and prepared in two different sizes (5 × 5 mm 2 and 2 × 2 mm 2 ). The accuracies of the developed Fusion model for detecting 2 × 2 mm 2 pieces of polymer, wood, and metal were 95%, 95%, and 81%, respectively, while the detection accuracies of the Fusion model for detecting 5 × 5 mm 2 pieces of polymer, wood, and metal were all 100%. The performance of the Fusion model was higher than the VNIR- and SWIR-based detection models by 18% and 5%, respectively, when F1 scores were compared, and by 38% and 5%, when average detection rates were compared. The study results suggested that the fusion of two HSI modalities could detect FMs more effectively than a single HSI modality.

36 MATERIALS SCIENCE↗

DOE OSTI2024

A static quantum embedding scheme based on coupled cluster theory

Here, we develop a static quantum embedding scheme that utilizes different levels of approximations to coupled cluster (CC) theory for an active fragment region and its environment. To reduce the computational cost, we solve the local fragment problem using a high-level CC method and address the environment problem with a lower-level Møller–Plesset (MP) perturbative method. This embedding approach inherits many conceptual developments from the hybrid second-order Møller–Plesset (MP2) and CC works by Nooijen [J. Chem. Phys. 111, 10815 (1999)] and Bochevarov and Sherrill [J. Chem. Phys. 122, 234110 (2005)]. We go beyond those works here by primarily targeting a specific localized fragment of a molecule and also introducing an alternative mechanism to relax the environment within this framework. We will call this approach MP-CC. We demonstrate the effectiveness of MP-CC on several potential energy curves and a set of thermochemical reaction energies, using CC with singles and doubles as the fragment solver, and MP2-like treatments of the environment. The results are substantially improved by the inclusion of orbital relaxation in the environment. Using localized bonds as the active fragment, we also report results for N=N bond breaking in azomethane and for the central C–C bond torsion in butadiene. We find that when the fragment Hilbert space size remains fixed (e.g., when determined by an intrinsic atomic orbital approach), the method achieves comparable accuracy with both a small and a large basis set. Additionally, our results indicate that increasing the fragment Hilbert space size systematically enhances the accuracy of observables, approaching the precision of the full CC solver.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Modern Approaches to Exact Diagonalization and Selected Configuration Interaction with the Adaptive Sampling CI Method

Recent advances in selected configuration interaction methods have made them competitive with the most accurate techniques available and, hence, creating an increasingly powerful tool for solving quantum Hamiltonians. In this work, we build on recent advances from the adaptive sampling configuration interaction (ASCI) algorithm. We show that a useful paradigm for generating efficient selected CI/exact diagonalization algorithms is driven by fast sorting algorithms, much in the same way iterative diagonalization is based on the paradigm of matrix vector multiplication. We present several new algorithms for all parts of performing a selected CI, which includes new ASCI search, dynamic bit masking, fast orbital rotations, fast diagonal matrix elements, and residue arrays. The ASCI search algorithm can be used in several different modes, which includes an integral driven search and a coefficient driven search. The algorithms presented here are fast and scalable, and we find that because they are built on fast sorting algorithms they are more efficient than all other approaches we considered. After introducing these techniques, we present ASCI results applied to a large range of systems and basis sets to demonstrate the types of simulations that can be practically treated at the full-CI level with modern methods and hardware, presenting double- and triple-ζ benchmark data for the G1 data set. The largest of these calculations is Si$_2$H$_6$ which is a simulation of 34 electrons in 152 orbitals. We also present some preliminary results for fast deterministic perturbation theory simulations that use hash functions to maintain high efficiency for treating large basis sets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

New procedure for evaluation of U(3) coupling and recoupling coefficients

A simple method to calculate Wigner coupling coefficients and Racah recoupling coefficients for U(3) in two group–subgroup chains is presented. While the canonical U(3) coupling and recoupling coefficients are applicable to any system that respects U(3) symmetry, the U(3) coupling coefficients are more specific to nuclear structure studies. This new procedure precludes the use of binomial coefficients and alternating sums which were used in the 1973 formulation of Draayer and Akiyama, and in so doing provides a faster and more accurate determination of any and all required results. The resolution of the outer multiplicity is based on the null space concept of the U(3) generators proposed by Alex et al., whereas the inner multiplicity in the angular momentum subgroup chain is obtained from the dimension of the null space of the SO(3) raising operator. It is anticipated that a C++ library will ultimately be available for determining generic coupling and recoupling coefficients associated with both the canonical and the physical group–subgroup chains of U(3).

Cross-Coupling Reaction↗

DOE OSTI2022

Strongly Correlated Quantum Spin Liquids versus Heavy Fermion Metals: A Review

This review considers the topological fermion condensation quantum phase transition (FCQPT) that explains the complex behavior of strongly correlated Fermi systems, such as frustrated insulators with quantum spin liquid and heavy fermion metals. The review contrasts theoretical consideration with recent experimental data collected on both heavy fermion metals (HF) and frustrated insulators. Such a method allows to understand experimental data. We also consider experimental data collected on quantum spin liquid in Lu 3 Cu 2 Sb 3 O 14 and quasi-one dimensional (1D) quantum spin liquid in both YbAlO 3 and Cu(C 4 H 4 N 2 )(NO 3 ) 2 with the aim to establish a sound theoretical explanation for the observed scaling laws, Landau Fermi liquid (LFL) and non-Fermi-liquid (NFL) behavior exhibited by these frustrated insulators. The recent experimental data on the heavy-fermion metal α-YbAl 1-x Fe x B 4 , with x=0.014, and on its sister compounds β-YbAlB 4 and YbCo 2 Ge 4 , carried out under the application of magnetic field as a control parameter are analyzed. We show that the thermodynamic and transport properties as well as the empirical scaling laws follow from the fermion condensation theory. We explain how both the similarity and the difference in the thermodynamic and transport properties of α-YbAl 1-x Fe x B 4 and in its sister compounds β-YbAlB 4 and YbCo 2 Ge 4 emerge, as well as establish connection of these (HF) metals with insulators Lu 3 Cu 2 Sb 3 O 14 , Cu(C 4 H 4 N 2 )(NO 3 ) 2 and YbAlO 3 . We demonstrate that the universal LFL and NFL behavior emerge because the HF compounds and the frustrated insulators are located near the topological FCQPT or are driven by the application of magnetic fields.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Small arms suppression project (LLNL final report)

US Special Operations Command (USSOCOM) was seeking a technological leap in small firearms weapon suppressor technology, because anticipated enemy capabilities are requiring the operators to have smaller detection cross sections to ensure the safe execution of missions. Suppressors have been developed almost exclusively through trial-and-error methods since the time of the original design by Hiram Maxim over one hundred years ago. Consequently USSOCOM deemed it prudent to perform a physicsbased study of weapon suppression to understand performance limits and possibly identify breakthrough technologies. Lawrence Livermore National Laboratory’s (LLNL’s) high performance production level computational tool called ALE3D (Arbitrary Lagrangian-Eulerian 3D and 2D) has unique physics models and numerical algorithms for modeling suppressor dynamics. The flexible and extendable code framework supports fully integrated hydrodynamics, heat transfer, solid and fluid dynamics, and chemistry that can be applied to simulating propellant-driven motion of a bullet down a gun barrel, the transfer of heat from the burning propellant to the barrel and suppressor, the chemistry of muzzle flash, and the shock/acoustic/optical signatures in the near-field. LLNL’s originally anticipated role was to augment ALE3D for this task, by developing the software and analysis methodologies specific to the simulation of blast and muzzle flash phenomena. It was believed that insights provided by our ALE3D simulations in tandem with a coordinated experimental component by our other team members from Oak Ridge National Laboratory (ORNL) and the U. S. Army Armament, Research, Development and Engineering Center (ARDEC), would have excellent prospects of yielding useful suppressor design improvements that could be transitioned to industry and utilized by US Special Operations Command. The three year effort has come to fruition with the development of revolutionary suppressor designs that far outperform any previous or current design by anyone outside this multi-lab team.

42 ENGINEERING↗

DOE OSTI2024

Adjusting the Energy Profile for CH–O Interactions Leads to Improved Stability of RNA Stem-Loop Structures in MD Simulations

The role of ribonucleic acid (RNA) in biology continues to grow, but insight into important aspects of RNA behavior is lacking, such as dynamic structural ensembles in different environments, how flexibility is coupled to function, and how function might be modulated by small molecule binding. In the case of proteins, much progress in these areas has been made by complementing experiments with atomistic simulations, but RNA simulation methods and force fields are less mature. It remains challenging to generate stable RNA simulations, even for small systems where well-defined, thermostable structures have been established by experiments. Further many different aspects of RNA energetics have been adjusted in force fields, seeking improvements that are transferable across a variety of RNA structural motifs. In this work, the role of weak CH···O interactions is explored, which are ubiquitous in RNA structure but have received less attention in RNA force field development. By comparing data extracted from high-resolution RNA crystal structures to energy profiles from quantum mechanics and force field calculations, it is shown that CH···O interactions are overly repulsive in the widely used Amber RNA force fields. A simple, targeted adjustment of CH···O repulsion that leaves the remainder of the force field unchanged was developed. Then, the standard and modified force fields were tested using molecular dynamics (MD) simulations with explicit water and salt, amassing over 300 μs of data for multiple RNA systems containing important features such as the presence of loops, base stacking interactions as well as canonical and noncanonical base pairing. In this work and others, standard force fields lead to reproducible unfolding of the NMR-based structures. Including a targeted CH···O adjustment in an otherwise identical protocol dramatically improves the outcome, leading to stable simulations for all RNA systems tested.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Deep Learning for In-Situ Layer Quality Monitoring during Laser-Based Directed Energy Deposition (LB-DED) Additive Manufacturing Process

Defects are a leading issue for the rejection of parts manufactured through the Directed Energy Deposition (DED) Additive Manufacturing (AM) process. In an attempt to illuminate and advance in situ quality monitoring and control of workpieces, we present an innovative data-driven method that synchronously collects sensing data and AM process parameters with a low sampling rate during the DED process. The proposed data-driven technique determines the important influences that individual printing parameters and sensing features have on prediction at the inter-layer qualification to perform feature selection. Three Machine Learning (ML) algorithms including Random Forest (RF), Support Vector Machine (SVM), and Convolutional Neural Network (CNN) are used. During post-production, a threshold is applied to detect low-density occurrences such as porosity sizes and quantities from CT scans that render individual layers acceptable or unacceptable. This information is fed to the ML models for training. Training/testing are completed offline on samples deemed “high-quality” and “low-quality”, utilizing only features recorded from the build process. CNN results show that the classification of acceptable/unacceptable layers can reach between 90% accuracy while training/testing on a “high-quality” sample and dip to 65% accuracy when trained/tested on “low-quality”/“high-quality” (respectively), indicating over-fitting but showing CNN as a promising inter-layer classifier.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Role of Intermolecular Interactions in Deep Eutectic Solvents for CO 2 Capture: Vibrational Spectroscopy and Quantum Chemical Studies

Recent research and reviews on CO 2 capture methods, along with advancements in industry, have highlighted high costs and energy-intensive nature as the primary limitations of conventional direct air capture and storage (DACS) methods. In response to these challenges, deep eutectic solvents (DESs) have emerged as promising absorbents due to their scalability, selectivity, and lower environmental impact compared to other absorbents. However, the molecular origins of their enhanced thermal stability and selectivity for DAC applications have not been explored before. Therefore, the current study focuses on a comprehensive investigation into the molecular interactions within an alkaline DES composed of potassium hydroxide (KOH) and ethylene glycol (EG). Combining Fourier transform infrared (FT-IR) and quantum chemical calculations, the study reports structural changes and intermolecular interactions induced in EG upon addition of KOH and its implications on CO 2 capture. Experimental and computational spectroscopic studies confirm the presence of noncovalent interactions (hydrogen bonds) within both EG and the KOH-EG system and point to the aggregation of ions at higher KOH concentrations. Additionally, molecular electrostatic potential (MESP) surface analysis, natural bond orbital (NBO) analysis, quantum theory of atoms-in-molecules (QTAIM) analysis, and reduced density gradient-noncovalent interaction (RDG-NCI) plot analysis elucidate changes in polarizability, charge distribution, hydrogen bond types, noncovalent interactions, and interaction strengths, respectively. Evaluation of explicit and hybrid models assesses their effectiveness in representing intermolecular interactions. This research enhances our understanding of molecular interactions in the KOH-EG system, which are essential for both the absorption and desorption of CO 2 . The study also aids in predicting and selecting DES components, optimizing their ratios with salts, and fine-tuning the properties of similar solvents and salts for enhanced CO 2 capture efficiency.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Optical and microstructural characterization of Er 3+ doped epitaxial cerium oxide on silicon

Rare-earth ion dopants in solid-state hosts are ideal candidates for quantum communication technologies, such as quantum memories, due to the intrinsic spin–photon interface of the rare-earth ion combined with the integration methods available in the solid state. Erbium-doped cerium oxide (Er:CeO 2 ) is a particularly promising host material platform for such a quantum memory, as it combines the telecom-wavelength (~ 1.5 μm) 4f–4f transition of erbium, a predicted long electron spin coherence time when embedded in CeO 2 , and a small lattice mismatch with silicon. In this work, we report on the epitaxial growth of Er:CeO 2 thin films on silicon using molecular beam epitaxy, with controlled erbium concentration between 2 and 130 parts per million (ppm). We carry out a detailed microstructural study to verify the CeO 2 host structure and characterize the spin and optical properties of the embedded Er 3+ ions as a function of doping density. In as-grown Er:CeO 2 in the 2–3 ppm regime, we identify an EPR linewidth of 245(1) MHz, an optical inhomogeneous linewidth of 9.5(2) GHz, an optical excited state lifetime of 3.5(1) ms, and a spectral diffusion-limited homogeneous linewidth as narrow as 4.8(3) MHz. We test the annealing of Er:CeO 2 films up to 900 °C, which yields narrowing of the inhomogeneous linewidth by 20% and extension of the excited state lifetime by 40%.

36 MATERIALS SCIENCE↗

DOE OSTI2023

The multichannel i -propyl + O2 reaction system: A model of secondary alkyl radical oxidation

The i-propyl + O2 reaction mechanism has been investigated by definitive quantum chemical methods to establish this system as a benchmark for the combustion of secondary alkyl radicals. Focal point analyses extrapolating to the ab initio limit were performed based on explicit computations with electron correlation treatments through coupled cluster single, double, triple, and quadruple excitations and basis sets up to cc-pV5Z. The rigorous coupled cluster single, double, and triple excitations/cc-pVTZ level of theory was used to fully optimize all reaction species and transition states, thus, removing some substantial flaws in reference geometries existing in the literature. The vital i-propylperoxy radical (MIN1) and its concerted elimination transition state (TS1) were found 34.8 and 4.4 kcal mol−1 below the reactants, respectively. Two β-hydrogen transfer transition states (TS2, TS2′) lie above the reactants by (1.4, 2.5) kcal mol−1 and display large Born–Oppenheimer diagonal corrections indicative of nearby surface crossings. An α-hydrogen transfer transition state (TS5) is discovered 5.7 kcal mol−1 above the reactants that bifurcates into equivalent α-peroxy radical hanging wells (MIN3) prior to a highly exothermic dissociation into acetone + OH. The reverse TS5 → MIN1 intrinsic reaction path also displays fascinating features, including another bifurcation and a conical intersection of potential energy surfaces. An exhaustive conformational search of two hydroperoxypropyl (QOOH) intermediates (MIN2 and MIN3) of the i-propyl + O2 system located nine rotamers within 0.9 kcal mol−1 of the corresponding lowest-energy minima.

Chemistry↗

DOE OSTI2023

Time-Gated Raman Spectroscopy (Intern Final Report)

Time-gated Raman spectroscopy is a powerful tool for rapid non-destructive analysis of chemical compounds, particularly when combined with other investigative methods. Time-gating the signal obtained from a sample during Raman spectroscopy allows for the removal of background emission and black body radiation, improving clarity and accuracy of scans. The main objectives of this project were to prove the lab’s time-gated Raman spectroscopy system works as proof of concept and to make select improvements to the in-house LabVIEW software used to run the time-gated system. This time-gated system will be added to the CTK’s TAP reactors in the future, but this is out of the scope of this internship project. This research is important, as understanding the mechanisms behind catalysts leads to designing better systems/materials and the perfection of catalysts, and energy-efficient chemical manufacturing is key when designing a greener future. The research for this project was done with samples at both ambient conditions in sample holders and at operando conditions (i.e., high temperature, unsteady-state) inside a small-volume chemical reactor. Numerous spectro-kinetic scans of the catalysts were obtained, but the advanced dynamics of these systems are still being interpreted. There is data that demonstrates that the time-gated system is correctly rejecting emission from samples, and one can observe real-time changes in the signal from the sample, due to coking or changes induced by redox chemistry. Chemical manufacturers and the planet are the main beneficiaries of catalysis research, as better catalysts will reduce carbon emissions.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Probe the Localized Electrochemical Environment Effects and Electrode Reaction Dynamics for Metal Batteries using In Situ 3D Microscopy

Uncontrollable dendrite growth is closely related to non-uniform reaction environments. However, there is a lack of understanding and analysis methods to probe the localized electrochemical environment (LEE). Here the effects of the LEE are investigated, including localized ion concentrations, current density, and electric potential, on metal plating/stripping dynamics and dendrite minimization. A novel in situ 3D microscopy technique is developed to image the morphology dynamics and deposition rate of Zn plating/stripping processes on 3D Zn–Mn anodes. Using the in situ 3D microscope, the electrode morphology changes during the reactions are directly imaged and Zn deposition rate maps at different time points are obtained. It is found that reaction kinetics are highly correlated to LEE and electrode morphology. To further quantify the LEE effects, the digital twin technique is employed that allows the accurate calculation of the electrochemical environments, such as localized ion concentrations, current density, and electric potential, which cannot be directly measured from experiments. We found that the curvature of the 3D electrode surface determines the LEE and significantly influences reaction kinetics. This provides a new strategy to minimize the dendrite formation by designing and optimizing the 3D geometry of the electrode to control the LEE.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Comparing GGA, GGA+ U , and meta-GGA functionals for redox-dependent binding at open metal sites in metal–organic frameworks

Metal–organic frameworks (MOFs) with open metal sites have been widely investigated for the selective adsorption of small molecules via redox mechanisms where charge transfer can take place between the binding site and the adsorbate of interest. Quantum-chemical screening methods based on density functional theory have emerged as a promising route to accelerate the discovery of MOFs with enhanced binding affinities toward various adsorbates. However, the success of this approach is linked to the accuracy of the underlying density functional approximations (DFAs). In this work, we compare commonly used generalized gradient approximation (GGA), GGA+U, and meta-GGA exchange-correlation functionals in modeling redox-dependent binding at open metal sites in MOFs using O 2 and N 2 as representative small molecules. We find that the self-interaction error inherent to the widely used Perdew, Burke, and Ernzerhof (PBE) GGA predicts metal sites that are artificially redox-active, as evidenced by their strong binding affinities, short metal–adsorbate bond distances, and large degree of charge transfer. The incorporation of metal-specific, empirical Hubbard U corrections based on the transition metal oxide literature systematically reduces the redox activity of the open metal sites, often improving agreement with experiment. Additionally, the binding behavior shifts from strong chemisorption to weaker physisorption as a function of U. The M06-L meta-GGA typically predicts binding energies between those of PBE-D3(BJ) and PBE-D3(BJ)+U when using empirically derived U values from the transition metal oxide literature. Despite the strong sensitivity of the binding affinities toward a given DFA, the GGA, GGA+U, and meta-GGA approaches often yield the same qualitative trends and structure–property relationships.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

The effect of rare-earth substitution on the Debye temperature of inorganic phosphors

In the quest to predict photoluminescent efficiency in rare-earth substituted inorganic phosphors, research has shown that materials consisting of ordered, rigid crystal structures tend to possess the highest photoluminescent quantum yields. A compound's Debye temperature (Θ D ), which can be calculated using ab initio calculations, is an ideal proxy for quantitatively comparing structural rigidity among different inorganic compounds, allowing potentially efficient phosphors to be selected from large crystal structure databases. However, the high computational cost of these calculations limits estimating Θ D for unsubstituted host crystal structures only. It is assumed that the low substitution concentration of the rare-earth luminescent center does not significantly influence a material's Debye temperature. This work evaluates the validity of this approximation by examining the effect of luminescent center substitution on a host structure's Θ D . Two well-known phosphors, (Y 1– x Ce x ) 3 Al 5 O 12 ( x = 0 – 0.05) and Ba 1-δ Eu δ MgAl 10 O 17 (δ = 0 – 0.15), were synthesized with varying rare-earth concentrations, while Θ D was computationally estimated and then determined by ultrasonic pulse-echo speed-of-sound and low-temperature heat capacity measurements. Overall, the ensuing results provide key implications for using Θ D as a proxy for structural rigidity in substituted inorganic compounds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗