SEARCH · Engineering Papers

Results for “inversion 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.

SourceResearch connectionFrom yearThrough year

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

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Coherent phonon pairs and rotational symmetry breaking of charge density wave order in the kagome superconductor CsV 3 Sb 5

Here, in this work we perform ultrafast time-resolved reflectivity measurements to study the symmetry breaking in the charge-density wave (CDW) phase of CsV 3 ⁢Sb 5 . By extracting the coherent phonon spectrum in the CDW phase of CsV 3 ⁢Sb 5 , we discover close phonon pairs near 1.3 THz and 3.1 THz, as well as a new mode at 1.84 THz. The 1.3 THz phonon pair and the 1.84 THz mode are observed up to the CDW transition temperature. Combining density-functional theory calculations, we show that these phonon pairs arise from the coexistence of Star of David and inverse Star of David distortions combined with sixfold rotational symmetry breaking. An anisotropy in the magnitude of transient reflectivity change is also revealed at the onset of CDW order. Our results thus indicate broken sixfold rotational symmetry in the charge-density wave state of CsV 3⁢ Sb 5 , along with the absence of nematic fluctuation above 𝑇 CDW . Meanwhile, the measured coherent phonon spectrum in the CDW phase of CsV 3 ⁢Sb 5−𝑥 ⁢Sn 𝑥 with 𝑥 = 0.03–0.04 matches the staggered inverse Star of David with interlayer 𝜋 phase shift. This CDW structure contrasts with undoped CsV3⁢Sb5 and explains the evolution from a phonon pair to a single mode at 1.3 THz by 𝑥 = 0.03–0.04 Sn doping.

Deng, Qinwen [Univ. of Pennsylvania, Philadelphia,↗

DOE OSTI2023

Large-Scale Materials Modeling at Quantum Accuracy: Ab Initio Simulations of Quasicrystals and Interacting Extended Defects in Metallic Alloys

Ab initio electronic-structure has remained dichotomous between achievable accuracy and length-scale. Quantum many-body (QMB) methods realize quantum accuracy but fail to scale. Density functional theory (DFT) scales favorably but remains far from quantum accuracy. We present a framework that breaks this dichotomy by use of three interconnected modules: (i) invDFT: a methodological advance in inverse DFT linking QMB methods to DFT; (ii) MLXC: a machine-learned density functional trained with invDFT data, commensurate with quantum accuracy; (iii) DFT-FE-MLXC: an adaptive higher-order spectral finite-element (FE) based DFT implementation that integrates MLXC with efficient solver strategies and HPC innovations in FE-specific dense linear algebra, mixed-precision algorithms, and asynchronous compute-communication. Furthermore, we demonstrate a paradigm shift in DFT that not only provides an accuracy commensurate with QMB methods in ground-state energies, but also attains an unprecedented performance of 659.7 PFLOPS (43.1% peak FP64 performance) on 619,124 electrons using 8,000 GPU nodes of Frontier supercomputer.

density functional theory↗

DOE OSTI2020

Approximate Inverse Chain Preconditioner: Iteration Count Case Study for Spectral Support Solvers

As the growing availability of computational power slows, there has been an increasing reliance on algorithmic advances. However, faster algorithms alone will not necessarily bridge the gap in allowing computational scientists to study problems at the edge of scientific discovery in the next several decades. Often, it is necessary to simplify or precondition solvers to accelerate the study of large systems of linear equations commonly seen in a number of scientific fields. Preconditioning a problem to increase efficiency is often seen as the best approach; yet, preconditioners which are fast, smart, and efficient do not always exist. Following the progress of [1], we present a new preconditioner for symmetric diagonally dominant (SDD) systems of linear equations. These systems are common in certain PDEs, network science, and supervised learning among others. Based on spectral support graph theory, this new preconditioner builds off of the work of [2], computing and applying a V-cycle chain of approximate inverse matrices. This preconditioner approach is both algebraic in nature as well as hierarchically-constrained depending on the condition number of the system to be solved. Due to its generation of an Approximate Inverse Chain of matrices, we refer to this as the AIC preconditioner. We further accelerate the AIC preconditioner by utilizing precomputations to simplify setup and multiplications in the con-text of an iterative Krylov-subspace solver. While these iterative solvers can greatly reduce solution time, the number of iterations can grow large quickly in the absence of good preconditioners. Initial results for the AIC preconditioner have shown a very large reduction in iteration counts for SDD systems as compared to standard preconditioners such as Incomplete Cholesky (ICC) and Multigrid (MG). We further show significant reduction in iteration counts against the more advanced Combinatorial Multigrid (CMG) preconditioner. We have further developed no-fill sparsification techniques to ensure that the computational cost of applying the AIC preconditioner does not grow prohibitively large as the depth of the V-cycle grows for systems with larger condition numbers. Our numerical results have shown that these sparsifiers maintain the sparsity structure of our system while also displaying significant reductions in iteration counts.1 2

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Approximate Inverse Chain Preconditioner: Iteration Count Case Study for Spectral Support Solvers

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

First-Principles Studies on Sc 2 RuZ (Z = Si, Ge, Sn) Inverse Heusler Alloys: Structural, Electronic, and Transport Properties

The continuous demand for efficient, nontoxic, and thermally stable materials for room-temperature energy conversion motivates the exploration of novel thermoelectric systems beyond the traditional magnetic Heusler alloys. While full and half-Heusler compounds, especially Co-, Ni-, and Mn-based systems, have demonstrated promising thermoelectric properties, their typically high operating temperatures and magnetic complexities limit their applicability in ambient thermal management. In this context, we investigate whether Sc-based inverse Heusler alloys can offer a viable nonmagnetic alternative with competitive thermoelectric performance. In this work, we perform a systematic first-principles study of the inverse Heusler compounds Sc 2 RuZ (Z = Si, Ge, Sn), focusing on their structural, electronic, mechanical, and thermodynamic-thermoelectric properties. Density Functional Theory (DFT) was employed to compute optimized lattice structures and band dispersion, while dynamical stability was assessed via phonon calculations. Thermoelectric transport coefficients, including Seebeck coefficient, electrical conductivity, and thermal conductivity, were estimated using the semiclassical Boltzmann transport theory within the constant relaxation time approximation. Our results show that all Sc 2 RuZ compounds are thermodynamically stable semiconductors with indirect band gaps of 0.12–0.16 eV and exhibit high elastic moduli, especially Sc 2 RuSn, which demonstrates superior stiffness and incompressibility. Importantly, all compounds display promising room-temperature thermoelectric characteristics, including high Seebeck coefficients and power factors. These findings reveal that Sc 2 RuZ alloys represent a rare class of stable, nonmagnetic inverse Heusler semiconductors with intrinsic thermoelectric potential at room temperature, unlike many existing Heusler systems optimized for spintronics or high-temperature operation. This work expands the known design space for Heusler-based thermoelectrics and offers a theoretical basis for experimental realization of efficient, low-temperature, nonmagnetic thermoelectric materials.

alloys↗

DOE OSTI2024

Force-Free Identification of Minimum-Energy Pathways and Transition States for Stochastic Electronic Structure Theories

Here, the accurate mapping of potential energy surfaces (PESs) is crucial to our understanding of the numerous physical and chemical processes mediated by atomic rearrangements, such as conformational changes and chemical reactions, and the thermodynamic and kinetic feasibility of these processes. Stochastic electronic structure theories, e.g., Quantum Monte Carlo (QMC) methods, enable highly accurate total energy calculations that in principle can be used to construct the PES. However, their stochastic nature poses a challenge to the computation and use of forces and Hessians, which are typically required in algorithms for minimum-energy pathway (MEP) and transition state (TS) identification, such as the nudged elastic band (NEB) algorithm and its climbing image formulation. Here, we present strategies that utilize the surrogate Hessian line-search method, previously developed for QMC structural optimization, to efficiently identify MEP and TS structures without requiring force calculations at the level of the stochastic electronic structure theory. By modifying the surrogate Hessian algorithm to operate in path-orthogonal subspaces and at saddle points, we show that it is possible to identify MEPs and TSs by using a force-free QMC approach. We demonstrate these strategies via two examples, the inversion of the ammonia (NH 3 ) molecule and the nucleophilic substitution (S N 2) reaction F – + CH 3 F → FCH 3 + F – . We validate our results using Density Functional Theory (DFT)- and Coupled Cluster (CCSD, CCSD(T))-based NEB calculations. We then introduce a hybrid DFT-QMC approach to compute thermodynamic and kinetic quantities, free energy differences, rate constants, and equilibrium constants that incorporates stochastically optimized structures and their energies, and show that this scheme improves upon DFT accuracy. Our methods generalize straightforwardly to other systems and other high-accuracy theories that similarly face challenges computing energy gradients, paving the way for highly accurate PES mapping, transition state determination, and thermodynamic and kinetic calculations at significantly reduced computational expense.

Iyer, Gopal R.↗

DOE OSTI2026

Theoretical Study of High Harmonic Generation in Monolayer NbSe 2

High harmonic generation (HHG) is a powerful probe of electron dynamics on attosecond to femtosecond time scales and has been successfully used to detect electronic and structural changes in solid-state quantum materials, including transition-metal dichalcogenides (TMDs). Among TMDs, bulk NbSe 2 exhibits charge density wave (CDW) order below 33 K and becomes superconducting below 7.3 K. Monolayer NbSe 2 also has superconducting and CDW behavior and is therefore interesting as a material whose different structural and electronic properties could be probed via HHG. Here, we predict the HHG response of the pristine 2H and CDW phases of monolayer NbSe 2 using real-time time-dependent density functional theory under the application of a simulated laser pulse excitation. We find that due to the lack of inversion symmetry in both monolayer phases, it is possible to excite even harmonics and that the even harmonics appear as the transverse components of the current response under excitations polarized along the zigzag direction of the monolayer, while odd harmonics arise from the longitudinal current response in all excitation directions. This suggests that the even and odd harmonic responses can be controlled via the polarization of the probing field, opening an avenue for potentially useful applications in optoelectronic devices.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Spectroscopy of 33 Mg with knockout reactions

In this work, the structure of 33 Mg was investigated by means of two knockout reactions, one-neutron removal from 34 Mg and one-proton removal from 34 Al. Using comparative analysis of the population of observed excited states in the residual 33 Mg, the nature of these states can be deciphered. In addition, the long-standing controversy about the parity of the 33 Mg ground state is resolved using momentum distribution analysis, showing a clear signature for negative parity. Partial cross section measurements are compared with the results of eikonal reaction theory combined with large-scale shell model calculations of this complex nucleus located in the island of inversion, where configuration mixing plays a major role.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Pressure effect on band inversion in AE Cd 2 As 2 ( AE =Ca, Sr, Ba)

Recent studies have predicted that magnetic EuCd 2 As 2 can host several different topological states depending on its magnetic order, including a single pair of Weyl points. Here we report on the bulk properties and band inversion induced by pressure in the nonmagnetic analogs AECd 2 As 2 (AE = Ca, Sr, Ba) as studied with density functional theory calculations. Under ambient pressure we find that these compounds are narrow band gap semiconductors, in agreement with experiment. In this work, the size of the band gap is dictated by both the increasing ionicity across the AE series which tends to increase the band gap, as well as the larger nearest neighbor Cd-As distance from increasing atomic size which can decrease the band gap because the conduction band edge is an antibonding state derived mostly from Cd 5s orbitals. The combination of these two competing effects results in a nonmonotonic change of the band gap size across the AE series with SrCd 2 As 2 having the smallest band gap among the three compounds. The application of negative pressure reduces this band gap and causes the band inversion between the Cd 5s and As 4p orbitals along the Γ-A direction to induce a pair of Dirac points. The topological nature of the Dirac points is then confirmed by finding the closed Fermi arcs on the ($10\bar{1}0$) surface.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Identifying Band Inversions in Topological Materials Using Diffusion Monte Carlo

Topological insulators are characterized by insulating bulk states and robust metallic surface states. Band inversion is a hallmark of topological insulators. At time-reversal invariant points in the Brillouin zone, spin–orbit coupling (SOC) induces a swapping of orbital character at the bulk band edges. Reliably detecting band inversion in solid-state systems with many-body methods would aid in identifying possible candidates for spintronics and quantum computing applications and improve our understanding of the physics behind topologically nontrivial systems. Density functional theory (DFT) methods are a well-established means of investigating these interesting materials due to their favorable balance of computational cost and accuracy but often struggle to accurately model the electron–electron correlations present in the many materials containing heavier elements. In this work, we develop a novel method to detect band inversion within continuum quantum Monte Carlo (QMC) methods that can accurately treat the electron correlation and spin–orbit coupling that are crucial to the physics of topological insulators. Our approach applies a momentum-space-resolved atomic population analysis throughout the first Brillouin zone utilizing the Löwdin method and the one-body reduced density matrix produced with diffusion Monte Carlo (DMC). We integrate this method into QMCPACK, an open source ab initio QMC package, so that these ground-state methods can be used to complement experimental studies and validate prior DFT work on predicting the band structures of correlated topological insulators. Here, we demonstrate this new technique on the topological insulator bismuth telluride, which displays band inversion between its Bi-p and Te-p states at the Γ-point. We show an increase in charge on the bismuth-p orbital and a decrease in charge on the tellurium-p orbital when comparing band structures with and without SOC. Additionally, we use our method to compare the degree of band inversion present in monolayer Bi 2 Te 3 , which has no interlayer van der Waals interactions, to that seen in the bilayer and bulk. The method presented here will enable future many-body studies of band inversion that can shed light on the delicate interplay between correlation and topology in correlated topological materials.

Band structure↗

DOE OSTI2024

Primitive quantum gates for an S U ( 2 ) discrete subgroup: Binary octahedral

We construct a primitive gate set for the digital quantum simulation of the 48-element binary octahedral ( BO ) group. This non-Abelian discrete group better approximates S U ( 2 ) lattice gauge theory than previous work on the binary tetrahedral group at the cost of one additional qubit—for a total of six—per gauge link. The necessary primitives are the inversion gate, the group multiplication gate, the trace gate, and the BO Fourier transform. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

Enhanced Second and Third Harmonic Generations in Metasurfaces Enabled through Point Group Design and Accidental Bound States in the Continuum

Here, we investigate enhanced second‑ and third‑harmonic generation (SHG and THG) in GaAs metasurfaces by combining electromagnetic point‑group designs with accidental bound states in the continuum (BICs). Using group theory, we predict the existence, suppression, and polarization selection rules of SHG and THG for metasurfaces belonging to the D 2h and C 2v point groups. In D 2h metasurfaces, inversion symmetry suppresses SHG while enabling polarization‑selective THG, which is strongly enhanced near accidental BICs with ultrahigh Q‑factors approaching 107. By breaking inversion symmetry to form C 2v metasurfaces through the introduction of nanogaps, SHG is now allowed and is dramatically enhanced near accidental BICs. Importantly, accidental BICs preserve polarization selection rules while providing a robust parameter range for high‑Q resonances, in contrast to quasi‑BICs induced by symmetry breaking. Our results demonstrate that the combination of point‑group theory and accidental‑BIC engineering offers a predictive and robust route toward highly efficient nonlinear metasurfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Synthesis and Characterization of Macrocyclic Ionic Liquids for CO 2 Separation

In this work, a series of ionic liquids with symmetric, macrocyclic cations composed of N-heterocycles connected by alkyl or ether linkages coupled with bis(trifluoromethylsulfonyl)imide or bis(pentafluoroethanesulfonyl)amide anions were synthesized to explore their application in CO 2 separation. These novel ionic liquids showed selective physical absorption of CO 2 compared to N 2 making them potential candidates for carbon capture. We examined the effect of cation structure on the physicochemical and separations properties of the ionic liquid and showed that these ionic liquids exhibit an inverse relationship between molar volume and CO 2 sorption capacity following the trend predicted by models based on regular solution theory. The fractional free volumes of the ionic liquids were lower than typical alkyl-imidazolium ionic liquids despite the cyclic structure of the cations which could be related to the flexibility of the alkyl and ether linkages. These results provide added insight into our understanding of structure–property relationships in CO 2 solubility in ionic liquids and afford a path toward enhancing the separations properties based on structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Local Excitations of a Charged Nitrogen Vacancy in Diamond with Multireference Density Matrix Embedding Theory

Here, we investigate the negatively charged nitrogen-vacancy center in diamond using periodic density matrix embedding theory (pDMET). To describe the strongly correlated excited states of this system, the complete active space self-consistent field (CASSCF) followed by n-electron valence state second-order perturbation theory (NEVPT2) was used as the impurity solver. Since the NEVPT2-DMET energies show a linear dependence on the inverse of the size of the embedding subspace, we performed an extrapolation of the excitation energies to the nonembedding limit using a linear regression. The extrapolated NEVPT2-DMET first triplet–triplet excitation energy is 2.31 eV and that for the optically inactive singlet–singlet transition is 1.02 eV, both in agreement with the experimentally observed vertical excitation energies of ∼2.18 eV and ∼1.26 eV, respectively. This is the first application of pDMET to a charged periodic system and the first investigation of the NV – defect using NEVPT2 for periodic supercell models.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Design-Space Exploration for Inverse-Design of Wind Turbine Blades Using Data-Driven Methods

The state-of-the-practice aerodynamic design methods for wind turbine blades is typically based on Blade-Element Momentum (BEM) theory using a pre-designed frozen family of airfoils. The airfoils are themselves typically designed using panel methods. The design of the next-generation of large flexible rotors will need to capture the non-linear aerodynamics and three-dimensional flow to reduce the levelized cost of wind energy. Data-driven methods for aerodynamic design using data generated by computational fluid dynamics offer an attractive alternative to BEM-based methods that captures the non-linear aerodynamics of the component airfoils as well as the root and tip sections. In this work, we develop and demonstrate a framework to "smartly" explore the relevant design space in combination with an appropriate automated CFD pipeline to evaluate the aerodynamics of each design. The design-space exploration framework uses appropriate perturbations to the airfoil shape and induction profile from a baseline shape in combination with the inverse-design using BEM. The resulting designs are evaluated using an automated CFD pipeline using the in-house CFD solver framework "Mercury". We perform verification and validation to establish the capability of the Mercury framework to predict the aerodynamic performance of wind turbines. The CFD simulation of the perturbed blade shapes are optimized to restart from the converged baseline simulation to reduce the computational time. Finally, we demonstrate the design-space exploration technique for the design of the outboard section and the full rotor using perturbations to the shape and operating conditions of the NREL 5-MW turbine.

aerodynamic design methods↗

DOE OSTI2020

LA17-V-InverseRadTransportProblem-PD2Kb. Final Project Report

In the inverse radiation transport problem, a radiation signal emanating from an unknown object is used to reconstruct geometric and material parameters of the object , e.g., the mass of special nuclear material (SNM), the isotopic composition of the SNM (e.g., the enrichment of 235U), the thickness of shielding, etc. The inverse radiation transport problem is one of the most fundamental and important in the emergency response mission, and has been addressed in previous NA-22 and NA-42/84 investments by using derivative-free optimization methods. Our project was aimed at (i) consolidating previous work by applying polynomial chaos expansion (PCE) for reducing significantly the number of forward model simulations required by derivative-free optimization methods; (ii) developing a general-purpose framework and stand-alone software for solving the inverse radiation transport problem by applying a new “predictive modeling” methodology (predictive modeling of coupled multi-physics systems, PM_CMPS) which enables the fusing of experimental and computational information; and (iii) developing several benchmarks (for unscattered gamma rays) for demonstrating the quantification of non-Gaussian features (e.g., skewness and kurtosis, which characterize asymmetries and long tails) of response distributions by deterministically and efficiently propagating parameter uncertainties (covariances and higher-order moments, where available) using high-order sensitivities computed by the recently developed comprehensive adjoint sensitivity analysis methodology (C-ASAM). With regard to (i), application of PCE, the most significant accomplishment of this project was an implementation of PCE into INVERSE that significantly reduced computation time for the DiffeRential Evolution Adaptive Metropolis (DREAM) and Differential Evolution (DE) methods. For inverse transport problems with measurements of neutron leakage, neutron multiplication, and neutron-induced gamma-rays, the solution times for the DREAM and DE methods were reduced by factors of 30-300, reducing run times from hours to minutes while maintaining the same accuracy as the standard (non-PCE) DREAM and DE. Speedups were further improved by exploring sparse grid quadrature sets using the Toolkit for Adaptive Stochastic Modeling and Non-Intrusive ApproximatioN (Tasmanian) toolset. An initial exploration of sparse grids indicated that further speedups factors of 5-10 are obtainable. The implementation of PCE in INVERSE is not in a form that users can easily access. With regard to (ii), solving the inverse problem using a predictive modeling paradigm, the most significant accomplishment of this project were the derivation and implementation of adjoint-based second derivatives in neutron and photon transport theory. There are two codes that compute these derivatives. SENSPG computes them for uncollided gamma rays and SENSMG computes them for neutron leakage. SENSMG is open-source code available on github. Other significant accomplishments included a stand-alone code, MULTI-PRED, that implements the predictive modeling paradigm, along with a detailed manual. The code was implemented in the existing INVERSE code for inverse problems, but a head-to-head comparison of inverse methods was not achieved. We also derived and implemented the capability to compute sensitivities for multiplicity measurements into the SENSMG code, and this capability was used in a data assimilation application, but it was not used for solving inverse problems. With regard to (iii), developing benchmarks to demonstrate the importance of high-order sensitivity and uncertainty analysis, the most significant accomplishment of this project were two test problems that showed how important second-order sensitivities are.

07 ISOTOPE AND RADIATION SOURCES↗

DOE OSTI2025

Discrepancies between Theory and Experiment in Determining the Ionization Energy of NF 3

High-accuracy ab initio thermochemical predictions for the ionization energy of NF 3 , the barrier height (to inversion) of NF 3 + , and the dissociative ionization threshold of NF 3 to NF 2 + + F are presented and incorporated into Active Thermochemical Tables. The adiabatic ionization energy of the first ionization band of NF 3 , calculated at 12.647 ± 0.010 eV, is at odds with previous experimental interpretations by nearly 0.36 eV due to unfavorable Franck-Condon factors associated with this transition. The barrier (to inversion) height is calculated to be about 0.6 eV lower in energy than the prior interpretation, which instigates a discussion of the supposed vibrational structure of the first ionization band of NF 3 . Updated assignments of the photoelectron spectrum are proposed, and the loss in vibrational spacing on the high-energy side of the experimental ionization band is discussed. Rudimentary anharmonic Franck-Condon simulations qualitatively reproduce the broad spectral features observed in experiment.

Active Thermochemical Tables↗

DOE OSTI2023

Unconventional surface state pairs in a high-symmetry lattice with anti-ferromagnetic band-folding

Many complex magnetic structures in a high-symmetry lattice can arise from a superposition of well-defined magnetic wave vectors. These “multi-q” structures have garnered much attention because of interesting real-space spin textures such as skyrmions. However, the role multi-q structures play in the topology of electronic bands in momentum space has remained rather elusive. Here we show that the type-I anti-ferromagnetic 1q, 2q and 3q structures in an face-centered cubic sublattice with band inversion, such as NdBi, can induce unconventional surface state pairs inside the band-folding hybridization bulk gap. Our density functional theory calculations match well with the recent experimental observation of unconventional surface states with hole Fermi arc-like features and electron pockets below the Neel temperature. We further show that these multi-q structures have Dirac and Weyl nodes. Our work reveals the special role that band-folding from anti-ferromagnetism and multi-q structures can play in developing new types of surface states.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗