SEARCH · Engineering Papers

Results for “Linear response”

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 163 records · Page 9

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Accurate prediction of core-level spectra of radicals at density functional theory cost via square gradient minimization and recoupling of mixed configurations

State-specific orbital optimized approaches are more accurate at predicting core-level spectra than traditional linear-response protocols, but their utility had been restricted due to the risk of "variational collapse" down to the ground state. We employ the recently developed square gradient minimization algorithm to reliably avoid variational collapse and study the effectiveness of orbital optimized density functional theory (DFT) at predicting second period element 1s core-level spectra of open-shell systems. Several density functionals (including SCAN, B3LYP, and ωB97X-D3) are found to predict excitation energies from the core to singly occupied levels with high accuracy (≤0.3 eV RMS error) against available experimental data. Higher excited states are, however, more challenging by virtue of being intrinsically multiconfigurational. We thus present a configuration interaction inspired route to self-consistently recouple single determinant mixed configurations obtained from DFT, in order to obtain approximate doublet states. This recoupling scheme is used to predict the C K-edge spectra of the allyl radical, the O K-edge spectra of CO + , and the N K-edge of NO 2 with high accuracy relative to experiment, indicating substantial promise in using this approach for the computation of core-level spectra for doublet species [vs more traditional time dependent DFT, equation of motion coupled cluster singles and doubles (EOM-CCSD), or using unrecoupled mixed configurations]. Finally, we also present general guidelines for computing core-excited states from orbital optimized DFT.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Polar liquids at charged interfaces: A dipolar shell theory

The structure of polar liquids and electrolytic solutions, such as water and aqueous electrolytes, at interfaces underlies numerous phenomena in physics, chemistry, biology, and engineering. In this work, we develop a continuum theory that captures the essential features of dielectric screening by polar liquids at charged interfaces, including decaying spatial oscillations in charge and mass, starting from the molecular properties of the solvent. The theory predicts an anisotropic dielectric tensor of interfacial polar liquids previously studied in molecular dynamics simulations. We explore the effect of the interfacial polar liquid properties on the capacitance of the electrode/electrolyte interface and on hydration forces between two plane-parallel polarized surfaces. In the linear response approximation, we obtain simple formulas for the characteristic decay lengths of molecular and ionic profiles at the interface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

An optimally tuned range-separated hybrid starting point for ab initio GW plus Bethe–Salpeter equation calculations of molecules

The ab initio GW plus Bethe–Salpeter equation (GW-BSE, where G is the one particle Green's function and W is the screened Coulomb interaction) approach has emerged as a leading method for predicting excitations in both solids and molecules with a predictive power contingent upon several factors. Among these factors are the (1) generalized Kohn–Sham eigensystem used to construct the GW self-energy and to solve the BSE and (2) the efficacy and suitability of the Tamm–Dancoff approximation. Here, we present a detailed benchmark study of low-lying singlet excitations from a generalized Kohn–Sham (gKS) starting point based on an optimally tuned range-separated hybrid (OTRSH) functional. We show that the use of this gKS starting point with one-shot G0W0 and G0W0-BSE leads to the lowest mean absolute errors (MAEs) and mean signed errors (MSEs), with respect to high-accuracy reference values, demonstrated in the literature thus far for the ionization potentials of the GW100 benchmark set and for low-lying neutral excitations of Thiel’s set molecules in the gas phase, without the need for self-consistency. The MSEs and MAEs of one-shot G0W0-BSE@OTRSH excitation energies are comparable to or lower than those obtained with other functional starting points after self-consistency. Additionally, we compare these results with linear-response time-dependent density functional theory (TDDFT) calculations and find GW-BSE to be superior to TDDFT when calculations are based on the same exchange-correlation functional. This work demonstrates tuned range-separated hybrids used in combination with GW and GW-BSE can greatly suppress starting point dependence for molecules, leading to accuracy similar to that for higher-order wavefunction-based theories for molecules without the need for costlier iterations to self-consistency.

McKeon, Caroline A. (ORCID:0000000217373503)↗

DOE OSTI2023

Size-dependent errors in real-time electron density propagation

Real-time (RT) electron density propagation with time-dependent density functional theory (TDDFT) or Hartree–Fock (TDHF) is one of the most popular methods to model the charge transfer in molecules and materials. However, both RT-TDHF and RT-TDDFT within the adiabatic approximation are known to produce inaccurate evolution of the electron density away from the ground state in model systems, leading to large errors in charge transfer and erroneous shifting of peaks in absorption spectra. Given the poor performance of these methods with small model systems and the widespread use of the methods with larger molecular and material systems, here we bridge the gap in our understanding of these methods and examine the size-dependence of errors in RT density propagation. We analyze the performance of RT density propagation for systems of increasing size during the application of a continuous resonant field to induce Rabi-like oscillations, during charge-transfer dynamics, and for peak shifting in simulated absorption spectra. We find that the errors in the electron dynamics are indeed size dependent for these phenomena, with the largest system producing the results most aligned with those expected from linear response theory. The results suggest that although the RT-TDHF and RT-TDDFT methods may produce severe errors for model systems, the errors in charge transfer and resonantly driven electron dynamics may be much less significant for more realistic, large-scale molecules and materials.

Chemistry↗

DOE OSTI2020

Superdiffusive transport of energy in one-dimensional metals

Metals in one spatial dimension are described at the lowest energy scales by the Luttinger liquid theory. It is well understood that this free theory, and even interacting integrable models, can support ballistic transport of conserved quantities including energy. In contrast, realistic one-dimensional metals, even without disorder, contain integrability-breaking interactions that are expected to lead to thermalization and conventional diffusive linear response. We argue that the expansion of energy when such a nonintegrable Luttinger liquid is locally heated above its ground state shows superdiffusive behavior (i.e., spreading of energy that is intermediate between diffusion and ballistic propagation), by combining an analytical anomalous diffusion model with numerical matrix-product–state calculations on a specific perturbed spinless fermion chain. Different metals will have different scaling exponents and shapes in their energy spreading, but the superdiffusive behavior is stable and should be visible in time-resolved experiments.

Bulchandani, Vir B.↗

DOE OSTI2020

Photoemission from hybrid states of Cl@C 60 before and after a stabilizing charge transfer

Photoionization calculations of the endofullerene molecule Cl@C 60 with an open-shell chlorine atom are performed in the linear response density functional theory based on a spherical jellium model of C 60 . Cross sections for atom–fullerene hybrid photoemission studied show the effects of the hybridization symmetry, the giant plasmon and the molecular cavity. Comparisons with the results of Ar@C 60 provide insights into the role of a shell-closing electron and its influence on the dynamics. The results for Cl@C 60 are further compared with those of a more stable configuration that results after a C 60 electron transfers to Cl forming Cl–@${{\text{C}}_{60}}^{+}$. This comparison reveals noticeable differences in the ionization properties of the antibonding hybrid state while the bonding hybrid remains nearly unaltered showing a magnification covering the entire giant plasmon energy range. Furthermore, these predictions should raise interest in probing the molecular configuration using photoelectron spectroscopy.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2024

First-principles calculation of Hubbard U for Terbium metal under high pressure

Abstract Using density functional theory (DFT) and linear response approaches, we compute the on-site Hubbard interaction U of elemental Terbium (Tb) metal in the pressure range ∼ 0–65 GPa. The resulting first-principles U values with experimental crystal structures enable us to examine the magnetic properties of Tb using a DFT+U method. The lowest-energy magnetic states in our calculations for different high-pressure Tb phases—including hcp, α -Sm, and dhcp—are found to be compatible with the corresponding magnetic ordering vectors reported in experiments. The result shows that the inclusion of Hubbard U substantially improves the accuracy and efficiency in modeling correlated rare-earth materials. Our study also provides the necessary U information for other quantum many-body techniques to study Tb under extreme pressure conditions.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Spin qubit properties of the boron-vacancy/carbon defect in the two-dimensional hexagonal boron nitride

Spin qubit defects in two-dimensional materials have a number of advantages over those in three-dimensional hosts including simpler technologies for defect creation and control, as well as qubit accessibility. In this work, we select the V B C B defect in the hexagonal boron nitride (hBN) as a possible optically controllable spin qubit and explain its triplet ground state and neutrality. In this defect a boron vacancy is combined with a carbon dopant substituting the closest boron atom to the vacancy. Our density-functional-theory calculations confirmed that the system has dynamically stable spin triplet and singlet ground states. As revealed from our linear response GW calculations, the spin-sensitive electronic states are localized around the three undercoordinated N atoms and make local peaks in the density of electronic states within the bandgap. Using the triplet and singlet ground state energies, as well as the energies of the optically excited states, obtained from solution to the Bethe–Salpeter equation, we construct the spin-polarization cycle, which is found to be favorable for the spin qubit initialization. The calculated zero-field splitting parameters ensure that the splitting energy between the spin projections in the triplet ground state is comparable to that of the known spin qubits. We thus propose the V B C B defect in hBN as a promising spin qubit.

2D BN↗

DOE OSTI2026

Perturbative second-order optical susceptibility of bulk materials: a symmetry-enforced return to non-orthogonal localized basis sets

The second-order optical susceptibility of semiconductors $\chi^{(2)}_{ijk}(-2\omega;\omega,\omega)$ finds application in metrology, spectroscopy, telecommunications, material characterization, and quantum information. Pioneering calculations of $\chi^{(2)}_{ijk}(-2\omega;\omega,\omega)$ utilized non-orthogonal Gaussian orbitals centered at atoms. That formulation transitioned into plane-wave-based algorithms as time went by. As of late, nevertheless, multiple tools for calculating optical susceptibilities have recast the problem using Wannier (i.e. localized) orbitals, making a comeback onto frameworks based on localized basis sets. Here, in this work, we present an approach for calculating $\chi^{(2)}_{ijk}(-2\omega;\omega,\omega)$ reliant on numerical pseudo-atomic orbitals (PAOs) within perturbation theory in the velocity gauge. Its salient feature is a calculation of ‘Slater–Koster-like’ two-center integrals of the momentum operator in between PAOs identified by symmetry. The approach was successfully tested on paradigmatic cubic silicon carbide (3C-SiC) and gallium arsenide, for which linear responses are contributed as well.

Huamán, Angiolo [Univ. of Arkansas, Fayetteville, ↗

DOE OSTI2022

Interpreting radial correlation Doppler reflectometry using gyrokinetic simulations

A linear response, local model for the DBS amplitude applied to gyrokinetic simulations shows that radial correlation Doppler reflectometry measurements (RCDR, Schirmer et al 2007 Plasma Phys. Control. Fusion 49 1019) are not sensitive to the average turbulence radial correlation length, but to a correlation length that depends on the binormal wavenumber $k_{\perp}$ selected by the Doppler backscattering (DBS) signal. Nonlinear gyrokinetic simulations show that the turbulence naturally exhibits a nonseparable power law spectrum in wavenumber space, leading to a power law dependence of the radial correlation length with binormal wavenumber $l_r \sim C k_{\perp}^{-\alpha} (\alpha \approx 1)$ which agrees with the inverse proportionality relationship between the measured $l_r$ and $k_{\perp} $ observed in experiments (Fernández-Marina et al 2014 Nucl. Fusion 54 072001). This new insight indicates that RCDR characterizes the eddy aspect ratio in the perpendicular plane to the magnetic field. It also motivates future use of a nonseparable turbulent spectrum to quantitatively interpret RCDR and potentially other turbulence diagnostics. The radial correlation length is only measurable when the radial resolution at the cutoff location $W_\text n$ satisfies $W_\text n \ll l_r$, while the measurement becomes dominated by $W_\text n$ for $W_\text n \gg l_r$. This suggests that $l_r$ is likely to be inaccessible for electron-scale DBS measurements (${k_{\perp}}{\rho_{\text s}} \gt 1$). The effect of $W_\text n$ on ion-scale radial correlation lengths could be nonnegligible.

radial correlation Doppler reflectometry↗

DOE OSTI2022

Application of machine learning and artificial intelligence to extend EFIT equilibrium reconstruction

Recent progress in the application of machine learning (ML)/artificial intelligence (AI) algorithms to improve the Equilibrium Fitting (EFIT) code equilibrium reconstruction for fusion data analysis applications is presented. A device-independent portable core equilibrium solver capable of computing or reconstructing equilibrium for different tokamaks has been created to facilitate adaptation of ML/AI algorithms. A large EFIT database comprising of DIII-D magnetic, motional Stark effect, and kinetic reconstruction data has been generated for developments of EFIT model-order-reduction (MOR) surrogate models to reconstruct approximate equilibrium solutions. Furthermore, a neural-network MOR surrogate model has been successfully trained and tested using the magnetically reconstructed datasets with encouraging results. Other progress includes developments of a Gaussian process Bayesian framework that can adapt its many hyperparameters to improve processing of experimental input data and a 3D perturbed equilibrium database from toroidal full magnetohydrodynamic linear response modeling using the Magnetohydrodynamic Resistive Spectrum - Feedback (MARS-F) code for developments of 3D-MOR surrogate models.

Gaussian process↗

DOE OSTI2022

Hubbard-corrected oxide formation enthalpies without adjustable parameters

A density functional theory (DFT) approach to computing transition metal oxide heat of formation without adjustable parameters is presented. Different degrees of d-electron localization in oxides are treated within the DFT+U approach with site-dependent, first-principles Hubbard U-parameters obtained from linear response theory, and delocalized states in the metallic phases are treated without Hubbard corrections. Comparison of relative stabilities of these differently treated phases is enabled by a local d-electron density matrix-dependent model, which was found by genetic programming against experimental reference formation enthalpies. This mathematically simple model does not explicitly depend on the Hubbard-corrected ionic species and is shown to reproduce the heats of formation of the Mott insulators Ca 2 RuO 4 and Y 2 Ru 2 O 7 within ~3% of experimental results, where the experimental training data did not contain Ru oxides. This newly developed method thus absolves from the need for element-specific corrections fitted to experiments in existing Hubbard-corrected approaches to the prediction of reaction energies of transition metal oxides and metals. The absence of fitting parameters opens up here the possibility to predict relative thermodynamic stabilities and reaction energies involving d-states of varying degree of localization at transition metal oxide interfaces and defects, where site-dependent U-parameters will be particularly important and devising a fitting scheme against experimental data with predictive power would be exceedingly difficult.

transition metal oxides↗

DOE OSTI2025

Interplay between electron localization, magnetic order, and Jahn-Teller distortion dictates LiMnO2 phase stability

The development of manganese (Mn)-rich cathodes for Li-ion batteries promises to alleviate potential supply chain bottlenecks in battery manufacturing. Fundamental challenges in Mn-rich cathodes arise from phenomena such as structural changes due to cooperative Jahn-Teller (JT) distortions of in octahedral environments, Mn migration, and phase transformations to spinel-like order, all of which affect the electrochemical performance. These physically complex phenomena motivate an re-examination of the Li-Mn-O rock-salt space, with a focus on the thermodynamics of the prototypical, polymorphs. It is found that the generalized gradient approximation (GGA-PBEsol) and meta-GGA ( ) density functionals with empirically fitted on-site Hubbard corrections yield spurious stable phases for , such as predicting a phase with -like order ( ) to be the ground state instead of the orthorhombic (Pmmn) phase, which is the experimentally known ground state. Accounting for antiferromagnetic order in each structure is shown to have a substantial effect on the total energies and resulting phase stability. By using hybrid-GGA (HSE06) and GGA with self-consistent Hubbard parameters (on-site and inter-site ) calculated from linear response theory, the experimentally observed phase stability trends are recovered. The calculated on-site between Mn- states in the experimentally observed orthorhombic, layered, and spinel phases are significantly smaller than in and disordered layered structures, by within GGA. The smaller values of are shown to be correlated with a collinear ordering of JT distortions, in which all orbitals are oriented in the same direction. This cooperative JT effect can lead to greater electron delocalization from Mn along the states due to increased Mn-O covalency, which contributes to the greater electronic stability compared to the phases with noncollinear JT arrangements. The structures with collinear ordering of JT distortions also generate greater vibrational entropy, which helps stabilize these phases at high temperature. These phases are shown to be strongly insulating with large calculated band gaps , which are computed using HSE06 and .

Kam, Ronald L↗

DOE OSTI2025

Current-based metrology with two-terminal mesoscopic conductors

The traditional approach to quantum parameter estimation focuses on the quantum state, deriving fundamental bounds on precision through the quantum Fisher information. In most experimental settings, however, performing arbitrary quantum measurements is highly unfeasible. In open quantum systems, an alternative approach to metrology involves the measurement of stochastic currents flowing from the system to its environment. However, the present understanding of current-based metrology is mostly limited to Markovian master equations. Considering a parameter estimation problem in a two-terminal mesoscopic conductor, we identify the key elements that determine estimation precision within the Landauer-Büttiker formalism. Crucially, this approach allows us to address arbitrary coupling and temperature regimes. Furthermore, we obtain analytical results for the precision in linear-response and zero-temperature regimes. For the specific parameter estimation task that we consider, we demonstrate that the boxcar transmission function is optimal for current-based metrology in all parameter regimes.

Landauer formula↗

DOE OSTI2020

Spectral-density estimation with the Gaussian integral transform

The spectral-density operator $\hat{ρ}(ω) = δ(ω–\hat{H})$ plays a central role in linear response theory as its expectation value, the dynamical response function, can be used to compute scattering cross sections. In this work, we describe a near optimal quantum algorithm providing an approximation to the spectral density with energy resolution $\Delta$ and error $\epsilon$ using $O(\sqrt{\text{log}_2 (1/ε)[\text{log}_2 (1 / Δ) + \text{log}_2 (1/ε)]/ Δ)}$ operations. This is achieved without using expensive approximations to the time-evolution operator, but instead exploiting qubitization to implement an approximate Gaussian integral transform of the spectral density. Finally, we also describe appropriate error metrics to assess the quality of the spectral function approximations more generally.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Bridging the gap of storage ring light sources and linac-driven free-electron lasers

High-gain free-electron lasers (FELs) are driven by short, high-charge density electron beams as only produced at dedicated single pass or recirculating linear accelerators. We describe new conceptual, technical, and modeling solutions to produce subpicosecond, up to ~100 μJ-energy extreme ultra-violet and soft x-ray FEL pulses at high- and tunable repetition rates, from diffraction-limited storage ring light source. In contrast to previously proposed schemes, we show that lasing can be simultaneous to the standard multibunch radiation emission from short insertion devices, and that it can be obtained with limited impact on the storage ring infrastructure. By virtue of the high-average power but moderate pulse energy, the storage ring-driven high-gain FEL would open the door to unprecedented accuracy in time-resolved spectroscopic analysis of matter in the linear response regime, in addition to inelastic scattering experiments.

43 PARTICLE ACCELERATORS↗

DOE OSTI2020

Generalized gradient approximations with local parameters

Here we develop and demonstrate the performance of a nonseparable form of the generalized gradient approximation (GGA) for exchange and correlation that includes locally varying parameters that match second-order gradient expansion behavior. Specifically, the high- and low-density limits are included to recover locally the linear response through inclusion of their dependence on the electron density. This local parametrization allows the GGA form to provide varying behavior depending on the density regime. On the basis of an extensive series of property calculations involving both molecules and solids, we show that this nonempirical methodology can lead to a balanced GGA description of both finite and extended systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Calculations of electronic excitation by protons and alpha particles in silicon

This paper presents an ab initio method of calculating the valence electron contribution to the stopping and interaction cross sections of swift ions in matter using linear response calculations of the dielectric function in crystals. We estimate the contribution to the stopping and cross sections by repeated electron capture and stripping of the projectile ion by the medium. We also predict the electronic stopping of protons and alpha particles in silicon and compare to measured values. The predictions generally agree with the data at ion speeds greater than 1 atomic unit. We find that electron capture and stripping is at most a small contribution to the stopping for both ions studied.

electron excitation, silicon, nuclear interactions↗