Engineering Papers⌕ Search

Engineering topics

Verma, Pragya

Publications and source records attributed to Verma, Pragya.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Reversible dehydrogenation and rehydrogenation of cyclohexane and methylcyclohexane by single-site platinum catalyst

Abstract Developing highly efficient and reversible hydrogenation-dehydrogenation catalysts shows great promise for hydrogen storage technologies with highly desirable economic and ecological benefits. Herein, we show that reaction sites consisting of single Pt atoms and neighboring oxygen vacancies (V O ) can be prepared on CeO 2 (Pt 1 /CeO 2 ) with unique catalytic properties for the reversible dehydrogenation and rehydrogenation of large molecules such as cyclohexane and methylcyclohexane. Specifically, we find that the dehydrogenation rate of cyclohexane and methylcyclohexane on such sites can reach values above 32,000 mol H2 mol Pt −1 h −1 , which is 309 times higher than that of conventional supported Pt nanoparticles. Combining of DRIFTS, AP-XPS, EXAFS, and DFT calculations, we show that the Pt 1 /CeO 2 catalyst exhibits a super-synergistic effect between the catalytic Pt atom and its support, involving redox coupling between Pt and Ce ions, enabling adsorption, activation and reaction of large molecules with sufficient versatility to drive abstraction/addition of hydrogen without requiring multiple reaction sites.

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

Stabilized open metal sites in bimetallic metal–organic framework catalysts for hydrogen production from alcohols

Liquid organic hydrogen carriers such as alcohols and polyols are a high-capacity means of transporting and reversibly storing hydrogen that demands effective catalysts to drive the (de)hydrogenation reactions under mild conditions. We employed a combined theory/experiment approach to develop MOF-74 catalysts for alcohol dehydrogenation and examine the performance of the open metal sites (OMS), which have properties analogous to the active sites in high-performance single-site catalysts and homogeneous catalysts. Methanol dehydrogenation was used as a model reaction system for assessing the performance of five monometallic M-MOF-74 variants (M = Co, Cu, Mg, Mn, Ni). Co-MOF-74 and Ni-MOF-74 give the highest H 2 productivity. However, Ni-MOF-74 is unstable under reaction conditions and forms metallic nickel particles. To improve catalyst activity and stability, bimetallic (NixMg 1-x )-MOF-74 catalysts were developed that stabilize the Ni OMS and promote the dehydrogenation reaction. An optimal composition exists at (Ni 0.32 Mg 0.68 )-MOF-74 that gives the greatest H2 productivity, up to 203 mL gcat -1 min -1 at 300 °C, and maintains 100% selectivity to CO and H 2 between 225–275 °C. The optimized catalyst is also active for the dehydrogenation of other alcohols. DFT calculations reveal that synergistic interactions between the open metal site and the organic linker lead to lower reaction barriers in the MOF catalysts compared to the open metal site alone. This work expands the suite of hydrogen-related reactions catalyzed by MOF-74 which includes recent work on hydroformulation and our earlier reports of aryl-ether hydrogenolysis. Moreover, it highlights the use of bimetallic frameworks as an effective strategy for stabilizing a high density of catalytically active open metal sites.

08 HYDROGEN↗
DOE OSTI2020

Status and Challenges of Density Functional Theory

Here, we discuss some of the challenges facing density functional theory (DFT) and recent progress in DFT for both ground and excited electronic states. We discuss key aspects of the results we have been able to obtain with the strategy of designing density functionals to have various ingredients and functional forms that are then optimized to accurately predict various types of properties and systems with as much universality as possible. Finally, we make specific recommendations of approximate density functionals that are well suited for particular kinds of applications.

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

M06-SX screened-exchange density functional for chemistry and solid-state physics

Significance Density functionals with Hartree–Fock exchange have been widely used for a wide range of chemical applications, but the nonlocal character of exchange makes long-range exchange computationally expensive for solid-state calculations with periodic boundary conditions, and full exchange is nonphysical for condensed-phase systems. Here, we present a screened-exchange (SX) density functional, M06-SX, that is especially designed to have good accuracy for both solid-state physics and chemical applications with less computational cost than full Hartree–Fock exchange. The M06-SX functional gives accuracy comparable to functionals with full Hartree–Fock exchange for predicting chemical properties of molecules, while also being practical with good accuracy for plane wave calculations on band gaps and lattice constants of solids.

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

M11plus: A Range-Separated Hybrid Meta Functional with Both Local and Rung-3.5 Correlation Terms and High Across-the-Board Accuracy for Chemical Applications

The way to improve Kohn–Sham density functional theory is to improve the exchange–correlation functionals, and functionals have been successively improved by adding new ingredients, especially local spin density gradients, nonlocal Hartree–Fock exchange, and local meta terms based on kinetic energy density. Here, we present a new kind of functional obtained by adding rung-3.5 terms to a functional including local gradients, local meta terms, and range-separated Hartree–Fock exchange. A rung-3.5 term has short-range nonlocality designed to account for nondynamic correlation; we add two kinds of rung-3.5 terms, one kind modeled on position-dependent Hartree–Fock exchange and another modeled on the spin density at a point interacting with the opposite-spin exchange hole at the same point. Here, optimization of the functional yields broad accuracy for both ground states and excited states with especially significant improvement for systems with strong correlation.

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

Revised M11 Exchange-Correlation Functional for Electronic Excitation Energies and Ground-State Properties

The ability of Kohn-Sham density functional theory (KS-DFT) to accurately predict various types of electronic excitation energies with (necessarily approximate) exchange-correlation functionals faces several challenges. Chief among these is that valence excitations are usually inherently multiconfigurational and therefore best treated by functionals with local exchange, whereas Rydberg and charge transfer excitations are often better treated with nonlocal exchange. The question arises of whether one can optimize a functional such that all three kinds of excitations (valence, Rydberg, and charge transfer – including long-range charge transfer) are treated in a balanced and accurate way. The goal of the present work is to try to answer that question and then to optimize a functional with the best possible balanced behavior. Of the variety of functional types available, we select range-separated hybrid meta functionals because (i) range separation allows the percentage of Hartree–Fock (HF) exchange to change with interelectronic separation, and therefore, one can have 100% HF exchange at large interelectronic separations, which gives good performance for long-range charge-transfer excitations, while the range separation allows one to simultaneously have smaller values of HF exchange at small and intermediate inter-electronic separations, which give good performance for valence and Rydberg excitations and (ii) meta functionals allow one to obtain better accuracy with high HF exchange than is possible with functionals whose local part depends only on spin densities and their gradients. Furthermore, this work starts with the range-separated hybrid meta functional, M11, and re-optimizes it (with stronger smoothness restraints) against electronic excitation energies and ground-state properties to obtain a new functional called revM11 that gives good performance for all three types of electronic excitations and at the same time gives very good predictions across-the-board for ground-state properties.

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

Is the Inversion of Phosphorus Trihalides (PF 3 , PCl 3 , PBr 3 , and PI 3 ) a Diradical Process?

This work explores possible reaction paths for the inversion of a series of trigonal pyramidal phosphorus trihalides, PF 3 , PCl 3 , PBr 3 , and PI 3 , and it especially addresses the question of whether and when the bonding of the lowest-energy species along the inversion paths should be described as a hyper open-shell diradical. The various paths for inversion are calculated using a single-reference method within the framework of Kohn-Sham density functional theory and also with multi-reference wave function methods. Our calculated results using both kinds of methods show that, for all the halogens studied (F, Cl, Br, and I), the lowestenergy singlet path for the inversion occurs by the formation of a C 2v transition structure rather than a D 3 h transition structure. This geometrical preference agrees with what has been inferred previously based on closed-shell singlet calculations. But in the present study, we examined not only closed-shell singlet transition states but also open-shell singlet states and triplet states for calculating stationary points and inversion paths, and for some of the phosphorus trihalides, we found that paths involving open-shell configurations are lower in energy than those restricted to closed-shell configurations. We analyzed the changes along the paths in terms of hybridization and orientation of the frontier orbitals and in terms of locally avoided crossings, and the extent of diradical character was quantified by calculating the effective number of unpaired electrons. Even for the singlet inversion path that goes via a D 3 h structure, the barrier for PF 3 , PCl 3 , and PBr 3 is higher for a closed-shell singlet spin state than for the open-shell singlet configuration. Furthermore, the energy of the triplet D 3 h structure is below even the open-shell D 3 h singlet for PCl 3 , PBr 3 , and PI 3 . This necessitates rethinking the role of open-shell states in nominally closedshell processes.

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

Combining Wave Function Methods with Density Functional Theory for Excited States

Here, we review state-of-the-art electronic structure methods based both on wave function theory (WFT) and density functional theory (DFT). Strengths and limitations of both the wave function and density functional based approaches are discussed, and modern attempts to combine these two methods are presented. The challenges in modeling excited-state chemistry using both single-reference and multireference methods are described. Topics covered include background, combining density functional theory with single-configuration wave function theory, generalized Kohn–Sham (KS) theory, global hybrids, range-separated hybrids, local hybrids, using KS orbitals in many-body theory (including calculations of the self-energy and the GW approximation), Bethe–Salpeter equation, algorithms to accelerate GW calculations, combining DFT with multiconfigurational WFT, orbital-dependent correlation functionals based on multiconfigurational WFT, building multiconfigurational wave functions from KS configurations, adding correlation functionals to multiconfiguration self-consistent-field (MCSCF) energies, combining DFT with configuration-interaction singles by means of time-dependent DFT, using range separation to combine DFT with MCSCF, embedding multiconfigurational WFT in DFT, and multiconfiguration pair-density functional theory.

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

Hyper Open-Shell Excited Spin States of Transition-Metal Compounds: FeF 2 , FeF 2 ···Ethane, and FeF 2 ···Ethylene

Spin-state energetics are important for understanding properties that involve more than one spin state, for example, catalysis occurring on two or more potential energy surfaces corresponding to different electronic spins. Very often, multiple-spin processes involve transition-metal compounds, and therefore, it is important to understand the electronic structure and energetics of such compounds in different spin states. In this work, we benchmark relative spin-state energies of FeF 2 with respect to the quintet ground spin state using both singleconfigurational and multi-configurational methods, and we examine how they are affected by the binding of ethane and ethylene to the iron center. We also benchmark the binding energies of the complexes. The single-configurational methods used in this work are the Hartree–Fock method, 32 exchange–correlation functionals, and the CCSD(T) coupled-cluster method in both restricted and unrestricted formalisms. The multi-configurational methods that have been used are CASSCF, CASPT2, CASPT3, MRCI, MRCI+Q, and MR-ACPF. The spin-state splitting energies depend on the functional chosen and of the 32 exchange-correlation functionals investigated here we find that for the septet and spin-projected triplet states of FeF 2 , the M06 functional is the best when compared to our best estimates from multi-reference calculations. If all the nine excitation energies are considered, where there are three excited spin states (singlet, triplet, and septet) for each of the three systems (FeF 2 , FeF 2 ···ethane, and FeF 2 ···ethylene), the three best performing functionals are HLE16, SOGGA11-X, and M06-2X. We find that the binding of ethane perturbs the relative spin-state energy of FeF 2 by only a small amount, but the stronger binding of ethylene has a larger effect. For the spin-state splitting energies of FeF 2 using single-reference CCSD(T), we found the predicted results depend very strongly on precisely how the calculations are done, in particular, on the spin-restricted or spin-unrestricted character of the SCF reference state, which can differ even by around 50 kcal/mol for the SCF reference state and the subsequent CCSD(T) calculations. Furthermore, on analyzing the wave functions of both the spinrestricted or spin-unrestricted formalisms, we find that the lowest-energy singlet and triplet states of the complexes, just like FeF 2 in isolation, can have more unpaired electrons than they are usually assumed to have, i.e. they can be hyper open-shell electronic configurations, and that this can significantly lower the energy.

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

Localizing Holes as Polarons and Predicting Band Gaps, Defect Levels, and Delithiation Energies of Solid-State Materials with a Local Exchange-Correlation Functional

Here, we assess the performance of seven exchange–correlation functionals (some with and some without a Hubbard U correction) for their ability (i) to predict band gaps of silicon, diamond, and Li-ion battery cathode materials, (ii) to localize hole polarons and predict delithiation energies in Li-ion battery cathode materials, and (iii) to predict transition levels of charge carriers of doped silicon and diamond. Both local and hybrid exchange–correlation functionals were tested. The local functionals tend to underestimate band gaps and delocalize polarons. The hybrid functionals very often give a good description of both properties, but they may not be practical for calculations involving large unit cells, large ensembles, or dynamics, and therefore a local functional with a Hubbard U correction is often used (giving the method called DFT+U), where the value of a parameter U is adjusted according to the system and the property being investigated. Keeping in mind the importance of computational cost and the undesirability of having to adjust an empirical parameter for each system or property of interest, we recently developed a local functional, namely HLE17, to try to accurately predict band gaps and excitation energies, and we validated it using mostly main-group solids and molecules. Here we test the performance of HLE17 for its ability to predict band gaps and localize polarons in other solid-state materials, and we compare its performance to that of popular local functionals (PBE, PBEsol, and TPSS), a range-separated hybrid functional with screened exchange (HSE06), and DFT+U. We find that HLE17 predicts more accurate band gaps than other local functionals and can localize holes as polarons, which other local functionals usually fail to do, and for a number of cases it is comparable in performance quality to Hubbard-corrected functionals without the need for system-specific parametrization and to hybrid functionals without the high cost. Because HLE17 does not predict accurate lattice constants, we use the single-point method of quantum chemistry, where the geometry is optimized with one functional and the band gap is calculated with HLE17, or we perform calculations with the lattice constants obtained by TPSS and both the fractional intracell coordinates and the electronic structure obtained by HLE17 (a new method denoted HLE17\\TPSS). In particular, we performed calculations by HLE17//TPSS, HLE17//HSE06, HLE17//DFT+U, and HLE17\\TPSS, and these methods usually agree well with each other and give values similar to experiment.

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

Assessment of electronic structure methods for the determination of the ground spin states of Fe( ii ), Fe( iii ) and Fe( iv ) complexes

Our ability to understand and simulate the reactions catalyzed by iron depends strongly on our ability to predict the relative energetics of spin states. In this work, we studied the electronic structures of Fe 2+ ion, gaseous FeO and 14 iron complexes using Kohn–Sham density functional theory with particular focus on determining the ground spin state of these species as well as the magnitudes of relevant spin-state energy splittings. The 14 iron complexes investigated in this work have hexacoordinate geometries of which seven are Fe(II), five are Fe(III) and two are Fe(IV) complexes. These are calculated using 20 exchange–correlation functionals. In particular, we use a local spin density approximation (LSDA) – GVWN5, four generalized gradient approximations (GGAs) – BLYP, PBE, OPBE and OLYP, two non-separable gradient approximations (NGAs) – GAM and N12, two meta-GGAs – M06-L and M11-L, a meta-NGA – MN15-L, five hybrid GGAs – B3LYP, B3LYP*, PBE0, B97-3 and SOGGA11-X, four hybrid meta-GGAs – M06, PW6B95, MPW1B95 and M08-SO and a hybrid meta-NGA – MN15. The density functional results are compared to reference data, which include experimental results as well as the results of diffusion Monte Carlo (DMC) calculations and ligand field theory estimates from the literature. For the Fe 2+ ion, all functionals except M11-L correctly predict the ground spin state to be quintet. However, quantitatively, most of the functionals are not close to the experimentally determined spin-state splitting energies. For FeO all functionals predict quintet to be the ground spin state. For the 14 iron complexes, the hybrid functionals B3LYP, MPW1B95 and MN15 correctly predict the ground spin state of 13 out of 14 complexes and PW6B95 gets all the 14 complexes right. Here, the local functionals, OPBE, OLYP and M06-L, predict the correct ground spin state for 12 out of 14 complexes. Two of the tested functionals are not recommended to be used for this type of study, in particular M08-SO and M11-L, because M08-SO systematically overstabilizes the high spin state, and M11-L systematically overstabilizes the low spin state.

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

Can Kohn–Sham density functional theory predict accurate charge distributions for both single-reference and multi-reference molecules?

Dipole moments are the first moment of electron density and are fundamental quantities that are often available from experiments. An exchange–correlation functional that leads to an accurate representation of the charge distribution of a molecule should accurately predict the dipole moments of the molecule. It is well known that Kohn–Sham density functional theory (DFT) is more accurate for the energetics of single-reference systems than for the energetics of multi-reference ones, but there has been less study of charge distributions. Here, in this work, we benchmark 48 density functionals chosen with various combinations of ingredients, against accurate experimental data for dipole moments of 78 molecules, in particular 55 single-reference molecules and 23 multi-reference ones. We chose both organic and inorganic molecules, and within the category of inorganic molecules there are both main-group and transition-metal-containing molecules, with some of them being multi-reference. As one would expect, the multi-reference molecules are not as well described by single-reference DFT, and the functionals tested in this work do show larger mean unsigned errors (MUEs) for the 23 multi-reference molecules than the single-reference ones. Five of the 78 molecules have relatively large experimental error bars and were therefore not included in calculating the overall MUEs. For the 73 molecules not excluded, we find that three of the hybrid functionals, B97-1, PBE0, and TPSSh (each with less than or equal to 25% Hartree–Fock (HF) exchange), the range-separated hybrid functional, HSE06 (with HF exchange decreasing from 25% to 0 as interelectronic distance increases), and the hybrid functional, PW6B95 (with 28% HF exchange) are the best performing functionals with each yielding an MUE of 0.18 D. Perhaps the most significant finding of this study is that there exists great similarity among the success rate of various functionals in predicting dipole moments. In particular, of 39 functionals designed as general-purpose functionals and that do not have a global value of 100% HF exchange, the average MUE is 0.23 D, with a standard deviation of only 0.04 D. Among gradient approximations, which are especially interesting because of their speed and portability, the best overall performance is by PBE, HCTH/407, OLYP, OreLYP, and GAM, each with MUE of 0.22 D.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗