Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density functional perturbation theory”

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

Magnetic Ni-N-Ni Centers in N-substituted NiO

Here, to explore how anion substitution modifies the existing magnetism in strongly correlated oxides, we investigate local electronic states and magnetic ordering in nickel oxide (NiO) induced by substituting oxygen (O) with nitrogen (N). Each N introduces an additional N 2p hole and modifies the magnetic moment of a neighboring nickel (Ni) cation site, as the exchange interaction between this hole and the Ni $e_{g}$ electrons exceeds the Ni-O-Ni superexchange interaction. This leads to the formation of Ni-N-Ni centers consisting of five spins, without perturbing the antiferromagnetic NiO lattice. These centers are studied using density functional theory and confirmed through high-resolution spectroscopy on N-substituted NiO thin films grown by molecular beam epitaxy. This type of magnetic design may enable future advances in quantum technologies based on strongly correlated materials, such as quantum sensors and spin-based qubits.

Godin, Simon [Univ. of British Columbia, Vancouver↗
DOE OSTI2023

Studying excited-state-specific perturbation theory on the Thiel set

Here, we explore the performance of a recently introduced N 5 -scaling excited-state-specific second order perturbation theory (ESMP2) on the singlet excitations of the Thiel benchmarking set. We find that, without regularization, ESMP2 is quite sensitive to π system size, performing well in molecules with small π systems but poorly in those with larger π systems. With regularization, ESMP2 is far less sensitive to π system size and shows a higher overall accuracy on the Thiel set than CC2, equation of motion-coupled cluster with singles and doubles, CC3, and a wide variety of time-dependent density functional approaches. Unsurprisingly, even regularized ESMP2 is less accurate than multi-reference perturbation theory on this test set, which can, in part, be explained by the set’s inclusion of some doubly excited states but none of the strong charge transfer states that often pose challenges for state-averaging. Beyond energetics, we find that the ESMP2 doubles norm offers a relatively low-cost way to test for doubly excited character without the need to define an active space.

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

Extrinsic n-type doping of Cd 3 As 2 thin films

Cd 3 As 2 provides an excellent platform for studying the physics of three-dimensional Dirac semimetals due to its stability as well as its compatibility with thin film growth. Crystals made using both bulk and thin film synthesis are unintentionally doped n-type, and other than introducing Zn to reduce the carrier concentration, no efforts have been reported to alter this intrinsic doping without major changes to the band structure. Here, group VI elements Te and Se are introduced during epitaxy to increase the electron concentration of the films. Starting from an unintentionally doped electron concentration of 1–2 × 10 17 cm -3 , concentrations of up to 3 × 10 18 cm -3 are achieved. Analysis of Shubnikov–de Haas oscillations reveals good agreement in calculated effective mass and Fermi velocity of highly doped films with unintentionally doped single crystals with similar electron concentrations. The density functional theory is also performed to study the effects of group VI substitutions and confirms no strong perturbations in the electronic structure. This work ultimately demonstrates tunability in the carrier concentration using extrinsic dopants without substantial changes in the band structure, allowing for intentional design of Fermi-level position for device applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Computation of Auger Electron Spectra in Organic Molecules with Multiconfiguration Pair-Density Functional Theory

Efficient and accurate computation of molecular Auger electron spectra for larger systems is limited by the rapid increase in the number of doubly ionized final states as the system size grows. Here, in this work, we benchmark the application of multiconfiguration pair-density functional theory with a restricted active space (RAS) reference wave function for computing the carbon K-edge decay spectra of 20 organic molecules. Decay rates are computed within the one-center approximation. We evaluate the performance of different basis sets and on-top functionals and find that multiconfiguration pair-density functional theory achieves accuracy comparable to RAS followed by second-order perturbation theory, but at significantly lower computational cost.

Fouda, Adam E. A. [Argonne National Laboratory (AN↗
DOE OSTI2026

Redox Potentials for Hydrogen Species in Acetonitrile and Tetrahydrofuran

The properties of the fundamental hydrogen species, proton (H+), hydrogen atom (H•), and hydride anion (H-) are critical to a vast range of chemical processes, yet their thermodynamic properties in nonaqueous solvents are not well established. A hybrid supermolecule-continuum approach is used to predict the Gibbs free energies of solvation (?G°solv) and standard redox potentials (E°) for the 2H+/H2 and H•/H- couples in acetonitrile (MeCN) and tetrahydrofuran (THF) following the approach previously used for water. Gas phase geometry optimizations were done using density functional theory (DFT) with the ?B97XD functional, Møller-Plesset second-order perturbation theory (MP2) with augmented correlation-consistent basis sets and composite-correlated molecular orbital theory (G3MP2 and G3MP2B3). Additional single point calculations in gas phase were performed for H- affinity for MeCN and THF at the coupled cluster CCSD(T) level using MP2 geometries. Solvation was included using the self-consistent reaction field SMD model. For H+, the ?G°solv values are predicted to be -252.2 kcal/mol in MeCN and -261.0 kcal/mol in THF. The redox potentials for the 2H+/H2 couple are +0.49 V in MeCN and +0.11 V in THF relative to the aqueous SHE, consistent with the available experimental data. The ?G°solv for H- is predicted to be -80 kcal/mol in MeCN and -68 in THF, and for H•, ?G°solv is predicted to be -1.3 kcal/mol in THF and -2.0 kcal/mol in MeCN. These solvation energies yield respective calculated redox potentials for the H•/H- couple of -0.14 V in MeCN and -0.68 V in THF relative to the aqueous SHE. H2 solvation is slightly exergonic in MeCN (-3.3 kcal/mol) and near thermoneutral in THF (+1.3 kcal/mol). The pKa(H2) is predicted to be 45.4 in MeCN and 47.8 in THF. The proton-coupled electron transfer (PCET) accounting term CG was evaluated giving values on the ferrocene scale within approximately 2 kcal/mol of literature in both solvents. These results provide benchmarks for hydrogen redox chemistry in MeCN and THF and improve consistency across E°, ?G°solv, pKa, and CG. These results provide benchmarks for hydrogen redox chemistry in nonaqueous media challenging previous estimates based on indirect extrapolations and assumptions.

Duda, Damian P.↗
DOE OSTI2020

Ligand effects on electronic structure and bonding in U(iii) coordination complexes: a combined MCD, EPR and computational study

The trivalent oxidation state of uranium has been shown to undergo unique reactivity, from its ability to activate a variety of small molecules to its role in the catalytic reduction of ethene to ethane amongst others. Central to this unique reactivity and ability to rationally design ligands for isotope separation is the underlying uranium electronic structure. While electronic structure studies of U(IV), U(V), and U(VI) have been extensive, by comparison, analogous studies of more reduced oxidation states such as U(III) remains underdeveloped. Herein we report a combined MCD and EPR spectroscopic approach along with density functional theory and multireference wavefunction calculations to elucidate the effects of ligand perturbation in three uranium(III) Tp* complexes. Overall, the experimental and computational insight suggests that the change in ligand environment across this series of U(III) complexes resulted in only minor perturbations in the uranium electronic structure. This combined approach was also used to redefine the electronic ground state of a U(III) complex with a redox non-innocent Bipy - ligand. Overall, these studies demonstrate the efficacy of the combined experimental and theoretical approach towards evaluating electronic structure and bonding in U(III) complexes and provide important insight into the challenges in altering ligand environments to modify bonding and reactivity in uranium coordination chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2022

Energy Landscape, and Structural and Spectroscopic Characterization of Diazirine and Its Cyclic Isomers

Identifying new nitrogenated hydrocarbon molecules in the interstellar medium (ISM) is challenging due to the lack of comprehensive spectroscopic data from experiments. In this computational work, we focus on investigating the structures, relative energies, spectroscopic constants, and energy landscape of the cyclic isomers of diazirine (c-CH2N2) using ab initio quantum chemical methods. Density functional theory (DFT) methods and coupled cluster theory with singles and doubles including perturbative triples [CCSD(T)] and CCSD(T) with the explicitly correlated F12 b correction [CCSD(T)-F12b] were employed for this purpose along with large correlation consistent cc-pVTZ, cc-pVQZ, and cc-pV5Z basis sets. Harmonic vibrational frequencies, infrared vibrational intensities, rotational constants, and dipole moments are reported. An harmonic vibrational fundamentals along with centrifugal distortion constants, and vibration-rotation interaction constants are also reported for all the cyclic isomers. The energies computed with the CCSD(T) and CCSD(T)-F12b methods were extrapolated to the one-particle complete basis set (CBS) limit following a three-point formula. At the CCSD(T)-F12b/CBS level of theory the 3,3H-diazirine (c-CH2N2) is the lowest energy cyclic isomer followed by 1,3H-diazirine, (E)-1,2H-Diazirine, and (Z)-1,2H-diazirine which are 20.1, 47.8, and 51.3kcal mol-1above the 3,3H-Diazirine, respectively. Accurate structures and spectroscopic constants that are reported here could be useful for future identification of these cyclic nitrogenated organic molecules in the interstellar medium or circumstellar disks.

Energy Landscape↗
DOE OSTI2023

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↗
DOE OSTI2022

Zero-Field Splitting Calculations by Multiconfiguration Pair-Density Functional Theory

Zero-field splitting (ZFS) is a fundamental molecular property that is especially relevant for single-molecule magnets (SMMs), electron paramagnetic resonance spectra, and quantum computing. Developing a method that can accurately predict ZFS parameters can be very powerful for designing new SMMs. One of the challenges is to include external correlation in an inherently multiconfigurational open-shell species for the accurate prediction of magnetic properties. Previously available methods depend on expensive multireference perturbation theory calculations to include external correlation. In this paper, we present spin-orbit-inclusive multiconfiguration and multistate pair-density functional theory (MC-PDFT) calculations of ZFSs; these calculations have a cost comparable to complete-active-space self-consistent field (CASSCF) theory, but they include correlation external to the active space. We found that combining a multistate formulation of MC-PDFT, namely, compressed-state multistate pair-density functional theory, with orbitals optimized by weighted-state-averaged CASSCF, yields reasonably accurate ZFS results

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

Modeling Liquid Water by Climbing up Jacob’s Ladder in Density Functional Theory Facilitated by Using Deep Neural Network Potentials

Within the framework of Kohn–Sham density functional theory (DFT), the ability to provide good predictions of water properties by employing a strongly constrained and appropriately normed (SCAN) functional has been extensively demonstrated in recent years. Here, we further advance the modeling of water by building a more accurate model on the fourth rung of Jacob’s ladder with the hybrid functional, SCAN0. In particular, we carry out both classical and Feynman path-integral molecular dynamics calculations of water with the SCAN0 functional and the isobaric–isothermal ensemble. To generate the equilibrated structure of water, a deep neural network potential is trained from the atomic potential energy surface based on ab initio data obtained from SCAN0 DFT calculations. For the electronic properties of water, a separate deep neural network potential is trained by using the Deep Wannier method based on the maximally localized Wannier functions of the equilibrated trajectory at the SCAN0 level. The structural, dynamic, and electric properties of water were analyzed. The hydrogen-bond structures, density, infrared spectra, diffusion coefficients, and dielectric constants of water, in the electronic ground state, are computed by using a large simulation box and long simulation time. For the properties involving electronic excitations, we apply the GW approximation within many-body perturbation theory to calculate the quasiparticle density of states and bandgap of water. Compared to the SCAN functional, mixing exact exchange mitigates the self-interaction error in the meta-generalized-gradient approximation and further softens liquid water toward the experimental direction. For most of the water properties, the SCAN0 functional shows a systematic improvement over the SCAN functional. However, some important discrepancies remain. The H-bond network predicted by the SCAN0 functional is still slightly overstructured compared to the experimental results.

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

Point Defects in Two-Dimensional Ruddlesden–Popper Perovskites Explored with Ab Initio Calculations

Two-dimensional Ruddlesden–Popper (RP) halide perovskites stand out as excellent layered materials with favorable optoelectronic properties for efficient light-emitting, spintronic, and other spin-related applications. However, properties often determined by defects are not well understood in these perovskite systems. This work investigates the ground state electronic structure of commonly formed defects in a typical RP perovskite structure by density functional theory. Our study reveals that these 2D perovskites generally retain their defect tolerance with limited perturbation of the electronic structure in the case of neutral-type point defects. In contrast, donor/acceptor defects induce deep midgap states, potentially causing harm to the material’s electronic performance. To retain positive intrinsic properties, the halide vacancies and interstitial defects should be avoided. The observed strong electron localization results in trap states and consequently leads to reduced device performance. Finally, this understanding can guide experimental efforts that aim for improved 2D halide perovskite-based device performance.

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

Spin-polarization anisotropy controlled by bending in tungsten diselenide nanoribbons and tunable excitonic states

A WSe 2 monolayer shows many interesting properties due to its spin-orbit coupling induced spin splitting in bands around the Fermi level and the spin-valley configuration. The orientation of the spin polarization in the relevant bands is crucial for the nature of exciton states and the optical valley selectivity. In this work, we studied the WSe 2 nanoribbons under different mechanical bending curvatures and electron/hole doping with density functional theory and their optical absorption and excitonic states with many-body perturbation GW and BSE (Bethe-Salpeter equation) methods. We found that the WSe 2 nanoribbons can exhibit an enhanced SOC effect and a spatially varying spin polarization in bands around the Fermi level under bending. The spin-polarization can show an anisotropy (or asymmetry) in those nearly degenerate bands, leading to a controllable magnetism via bending and electron/hole doping to the nanoribbons, suggesting a potential application in compact and controllable magnetic nanodevices and spintronics. The optical absorption spectrum of the nanoribbon presents a large tunability with bending within the near infrared region of about 0.4 to 1.5 eV, showing an enhanced absorption at a large bending condition. Finally, the exciton states generally show mixed or various spin configuration in the electron and hole pairs that are controlled by bending, potentially useful for applications in spin-based quantum information processes.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Density functional study of relaxation of adsorbate vibration modes: Dominance of anharmonic interaction

Formulation and density functional workflow for calculating the lifetime of vibrational modes of molecular adsorbates on solid surfaces due to vibration–phonon coupling are presented. The anharmonic coupling is invoked to give the correct description of the origin of temperature dependence. Using pyrrolidine (C 4 H 9 N) absorbed on the Cu(001) surface as a concrete example, here we show that the anharmonic coupling can be one to two orders more significant than the harmonic interaction for the broadening of vibrational spectra, especially as temperature increases. These results challenge the common assumption that the anharmonic interaction is weak and call for attention of considering its effect in quantum relaxation and related problems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1997

The Vibrational Frequencies of CaO2, ScO2, and TiO2: A Comparison of Theoretical Methods

The vibrational frequencies of several states of CaO2, ScO2, and TiO2 are computed at using density functional theory (DFT), the Hatree-Fock approach, second order Moller-Plesset perturbation theory (MP2), and the complete-active-space self-consistent-field theory. Three different functionals are used in the DFT calculations, including two hybrid functionals. The coupled cluster singles and doubles approach including the effect of unlinked triples, determined using perturbation theory, is applied to selected states. The Becke-Perdew 86 functional appears to be the cost effective method of choice, although even this functional does not perform well for one state of CaO2. The MP2 approach is significantly inferior to the DFT approaches.

Rosi, Marzio↗
DOE OSTI2021

Origin of Insulating Ferromagnetism in Iron Oxychalcogenide Ce 2 O 2 FeSe 2

An insulating ferromagnetic (FM) phase exists in the quasi-one-dimensional iron oxychalcogenide Ce 2 O 2 FeSe 2 , but its origin is unknown. To understand the FM mechanism, here a systematic investigation of this material is provided, analyzing the competition between ferromagnetic and antiferromagnetic tendencies and the interplay of hoppings, Coulomb interactions, Hund’s coupling, and crystal-field splittings. Overall, our intuitive analysis based on second-order perturbation theory shows that large entanglements between doubly occupied and half filled orbitals play a key role in stabilizing the FM order in Ce 2 O 2 FeSe 2 . In addition, via many-body computational techniques applied to a multiorbital Hubbard model, the phase diagram confirms the proposed FM mechanism.

1-dimensional spin chains↗
DOE OSTI2023

Accurate prediction of global-density-dependent range-separation parameters based on machine learning

In this work, we develop an accurate and efficient XGBoost machine learning model for predicting the global-density-dependent range-separation parameter, ωGDD, for long-range corrected functional (LRC)-ωPBE. This ωGDDML model has been built using a wide range of systems (11 466 complexes, ten different elements, and up to 139 heavy atoms) with fingerprints for the local atomic environment and histograms of distances for the long-range atomic correlation for mapping the quantum mechanical range-separation values. The promising performance on the testing set with 7046 complexes shows a mean absolute error of 0.001 117 a0−1 and only five systems (0.07%) with an absolute error larger than 0.01 a0−1, which indicates the good transferability of our ωGDDML model. In addition, the only required input to obtain ωGDDML is the Cartesian coordinates without electronic structure calculations, thereby enabling rapid predictions. LRC-ωPBE(ωGDDML) is used to predict polarizabilities for a series of oligomers, where polarizabilities are sensitive to the asymptotic density decay and are crucial in a variety of applications, including the calculations of dispersion corrections and refractive index, and surpasses the performance of all other popular density functionals except for the non-tuned LRC-ωPBE. Finally, LRC-ωPBE (ωGDDML) combined with (extended) symmetry-adapted perturbation theory is used in calculating noncovalent interactions to further show that the traditional ab initio system-specific tuning procedure can be bypassed. The present study not only provides an accurate and efficient way to determine the range-separation parameter for LRC-ωPBE but also shows the synergistic benefits of fusing the power of physically inspired density functional LRC-ωPBE and the data-driven ωGDDML model.

Chemistry↗
DOE OSTI2022

The exact exchange–correlation potential in time-dependent density functional theory: Choreographing electrons with steps and peaks

The time-dependent exchange–correlation potential has the unusual task of directing fictitious non-interacting electrons to move with exactly the same probability density as true interacting electrons. This has intriguing implications for its structure, especially in the non-perturbative regime, leading to step and peak features that cannot be captured by bootstrapping any ground-state functional approximation. Here, we review what has been learned about these features in the exact exchange–correlation potential of time-dependent density functional theory in the past decade or so and implications for the performance of simulations when electrons are driven far from any ground state.

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

Efficient Iridium Catalysts for Formic Acid Dehydrogenation: Investigating the Electronic Effect on the Elementary β-Hydride Elimination and Hydrogen Formation Steps

We report herein a series of Cp*Ir complexes containing a rigid 8-aminoquinolinesulfonamide moiety as highly efficient catalysts for the dehydrogenation of formic acid (FA). The complex [Cp*Ir(L)Cl] (HL = N-(quinolin-8-yl)benzenesulfonamide) displayed a high turnover frequency (TOF) of 2.97 × 10 4 h –1 and a good stability (>100 h) at 60 °C. Comparative studies of [Cp*Ir(L)Cl] with the rigid ligand and [Cp*Ir(L')Cl] (HL' = N-propylpypridine-2-sulfonamide) without the rigid aminoquinoline moiety demonstrated that the 8-aminoquinoline moiety could dramatically enhance the stability of the catalyst. The electron-donating ability of the N,N'-chelating ligand was tuned by functionalizing the phenyl group of the L ligand with OMe, Cl, and CF 3 to have a systematical perturbation of the electronic structure of [Cp*Ir(L)Cl]. Experimental kinetic studies and density functional theory (DFT) calculations on this series of Cp*Ir complexes revealed that (i) the electron-donating groups enhance the hydrogen formation step while slowing down the β-hydride elimination and (ii) the electron-withdrawing groups display the opposite effect on these reaction steps, which in turn leads to lower optimum pH for catalytic activity compared to the electron-donating groups.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗