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

Extending quantum-mechanical benchmark accuracy to biological ligand-pocket interactions

Predicting the binding affinity of ligands to protein pockets is key in the drug design pipeline. The flexibility of ligand-pocket motifs arises from a range of attractive and repulsive electronic interactions during binding. Accurately accounting for all interactions requires robust quantum-mechanical (QM) benchmarks, which are scarce for ligand-pocket systems. Additionally, disagreement between “gold standard” Coupled Cluster (CC) and Quantum Monte Carlo (QMC) methods casts doubt on many benchmarks for larger non-covalent systems. We introduce the “QUantum Interacting Dimer” (QUID) benchmark framework containing 170 non-covalent (non-)equilibrium systems modeling chemically and structurally diverse ligand-pocket motifs. Symmetry-adapted perturbation theory shows that QUID broadly covers non-covalent binding motifs and energetic contributions. Robust binding energies are obtained using complementary CC and QMC methods, achieving agreement of 0.5 kcal/mol. The benchmark data analysis reveals that several dispersion-inclusive density functional approximations provide accurate energy predictions, though their atomic van der Waals forces differ in magnitude and orientation. Contrarily, semiempirical methods and empirical force fields require improvements in capturing non-covalent interactions (NCIs) for out-of-equilibrium geometries. The wide span of NCIs, highly accurate interaction energies, and analysis of molecular properties take QUID beyond the “gold standard” for QM benchmarks of ligand-protein systems.

Puleva, Mirela [University of Luxembourg, Luxembou↗
DOE OSTI2023

Sensitivity measurements for a 250 MHz quartz shear-horizontal surface acoustic wave biosensor under liquid viscous loading

Surface acoustic wave (SAW) devices have been used in biochemical assays due to their high sensitivity. The device sensitivity is a function of changes in the density and viscosity of the liquid. Here, we studied the effect of fluid viscosity using a 250 MHz quartz shear-horizontal (SH)-SAW biosensor by monitoring different concentrations of binary aqueous/glycerol solutions. In this study, the sensitivity of the biosensor was determined by fitting the data to models derived from perturbation theory. Measurements in water were used as the reference. For a 0% to 50% glycerol solution, an 87°–204° separation in the phase shift was observed. The slope of the plot of the phase shift vs (ηρ)0.5 was used to indicate the sensor’s sensitivity. The sensitivity for our 250 MHz quartz SH-SAW sensors was calculated to be 3.7×10−3m2sKg. The corresponding mass sensitivity was determined to be 9.25 × 105m2Kg. The limit of detection was calculated to be 36 picograms (pg), while the limit of quantification or LOQ was calculated to be 109 pg. Traditionally, liquid phase measurements have been challenging for SAW devices because liquids dampen the vibrating sensors severely. This problem has been largely solved using a transverse (shear) wave instead of the more popular longitudinal or Rayleigh waves. Liquid measurements are now possible using transverse waves, also known as shear waves, because transverse waves are only minimally attenuated by liquids. Shear-horizontal SAW sensors (SH-SAW) show great promise as label-free biosensors because of their ability to handle liquid samples. However, the viscosity of the liquid still induces loading effects and can be measured when the liquid is loaded onto the SH-SAW propagating surface (delay line). When the liquid above the delay line is perturbed by physical or chemical changes, such as binding to a receptor, it alters the propagating acoustic wave. The SH-SAW device can measure these changes in liquid properties as a change in the wave’s phase compared to the original wave. The device’s phase shift was recorded as a function of the changes in the density and viscosity of the binary glycerol solution and used to determine the sensitivity in the linear dynamic range of responses.

Materials Science↗
DOE OSTI2023

The connection between nonzero density and spontaneous symmetry breaking for interacting scalars

We consider U(1)-symmetric scalar quantum field theories at zero temperature. At nonzero charge densities, the ground state of these systems is usually assumed to be a superfluid phase, in which the global symmetry is spontaneously broken along with Lorentz boosts and time translations. We show that, in d > 2 spacetime dimensions, this expectation is always realized at one loop for arbitrary non-derivative interactions, confirming that the physically distinct phenomena of nonzero charge density and spontaneous symmetry breaking occur simultaneously in these systems. We quantify this result by deriving universal scaling relations for the symmetry breaking scale as a function of the charge density, at low and high density. Moreover, we show that the critical value of μ above which a nonzero density develops coincides with the pole mass in the unbroken, Poincaré invariant vacuum of the theory. The same conclusions hold non-perturbatively for an O(N) theory with quartic interactions in d = 3 and 4, at leading order in the 1/N expansion. We derive these results by computing analytically the zero-temperature, finite-μ one-loop effective potential, paying special attention to subtle points related to the iε terms. We check our results against the one-loop low-energy effective action for the superfluid phonons in λΦ 4 theory in d = 4 previously derived by Joyce and ourselves, which we further generalize to arbitrary potential interactions and arbitrary dimensions. As a byproduct, we find analytically the one-loop scaling dimension of the lightest charge-n operator for the λΦ 6 conformal superfluid in d = 3, at leading order in 1/n, reproducing a numerical result of Badel et al. For a λΦ 4 superfluid in d = 4, we also reproduce the Lee-Huang-Yang relation and compute relativistic corrections to it. Finally, we discuss possible extensions of our results beyond perturbation theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Using Light to Enhance Charge Transfer in Battery Materials: Increasing Charging Rates in LiNi 0.5 Mn 1.5 O 4

Electrochemical, spectroscopic, and computational studies on the interplay between light perturbation and coupled electron and Li-ion transfer in LiNi 0.5 Mn 1.5 O 4 (LNMO) cathodes identify that photon perturbation moves the system out of equilibrium, increases entropy, lowers the impedance of the battery, and increases the charging capacity by at least 15% during fast charging by promoting the oxidation of Ni 3+ to Ni 4+ . Here, the correlation with density functional theory calculations points out that oxidation of Ni 3+ faces a 30% higher barrier compared with the oxidation of Ni 2+ and hence exhibits greater responses to the energy transferred from photons. Structural analysis elucidates that photon energy also assists the transition of LNMO from the Jahn−Teller distorted asymmetric lattice with Ni 3+ to the nondistorted Ni 4+ . Synergistically, photons energize the removal of electrons from Ni 3+ and create photoinduced Ni 4+ intermediates, followed by the removal of Li + with almost two times faster diffusion rates.

5 V spinel↗
NASA NTRS1973

The dynamic behavior and compliance of a stream of cavitating bubbles.

Study of the dynamic response of streams of cavitating bubbles to imposed pressure fluctuations to determine the role played by turbopump cavitation in the POGO instability of liquid rockets. Both quasi-static and more general linearized dynamic analyses are made of the perturbations to a cavitating flow through a region of reduced pressure in which the bubbles first grow and then collapse. The results, when coupled with typical bubble number density distribution functions, yield compliances which compare favorably with the existing measurements. Since the fluids involved are frequently cryogenic, a careful examination was made of the thermal effects both on the mean flow and on the perturbations. As a result, the discrepancy between theory and experiment for particular engines could be qualitatively ascribed to reductions in the compliance caused either by these thermal effects or by relatively high reduced frequencies.

Brennen, C.↗
NASA NTRS1994

On the Interaction of CO and NH3 with BH3 and BF3

BH3CO, BF3CO, BH3NH3, and BF3NH3 have been studied using density functional theory (DFT). The computed geometries are in reasonable agreement with experiment. BF3CO is found to be weakly bound. The binding energies of BH3CO, BH3NH3 , and BF3NH3 have also been computed using the coupled cluster singles and doubles level of theory, including a perturbational estimate of the connected triples [CCSD(T)]. These CCSD(T) calculations are performed using correlation consistent polarized-valence triple-zeta basis sets. These are expected to be the most accurate binding energies to date. The DFT binding energies are in good agreement with those obtained at the CCSD(T) level. The bonding is analyzed using the constrained space orbital variation (CSOV) technique. In addition to analyzing the bonding in BH3CO, BH3NH3, and BF3NH3, the CSOV approach is used to study why the bonding in BF3CO is so weak.

Bauschlicher, Charles W., Jr.↗
DOE OSTI2021

Metal-to-Ligand Charge-Transfer Spectrum of a Ru-Bipyridine-Sensitized TiO 2 Cluster from Embedded Multiconfigurational Excited-State Theory

Understanding optical properties of the dye molecule in dye-sensitized solar cells (DSSCs) from first-principles quantum mechanics can contribute to improving the efficiency of such devices. While density functional theory (DFT) and time-dependent DFT have been pivotal in simulating optoelectronic properties of photoanodes used in DSSCs at the atomic scale, questions remain regarding DFT’s adequacy and accuracy to furnish critical information needed to understand the various excited-state processes involved. In this work, we simulate the absorption spectra of a dye-sensitized solar cell analogue, comprised of a Ru-bipyridine (Ru-bpy) dye molecule and a small TiO2 cluster via DFT and via an accurate embedded correlated wavefunction (CW) theory. We generated CW spectra for the adsorbed Ru-bpy dye via a recently introduced capped density functional embedding theory or capped-DFET (to generate the embedding potential that accounts for the interaction of the molecule and the TiO 2 cluster). We then combined capped-DFET with the accurate but expensive multiconfigurational complete active space second-order perturbation theory (CASPT2)—embedded CASPT2. Because the CW theory is conducted on only a portion of the total system in the presence of an embedding potential that describes that portion’s interaction with its environment, we efficiently obtain CW-quality predictions that reflect local properties of the entire system. Specifically, for example, with capped-DFET and embedded CW theory, we can simulate accurately a plethora of metal-to-ligand charge-transfer excited properties at a manageable computational cost. Here, we predict detailed electronic spectra within the visible region, featuring the lowest three singlet and triplet excited states, along with predictions of the singlets’ lifetimes. We illustrated these results using a Jablonski diagram that show the relative energy position of the singlet and longer-lived triplet excited states and analyzed and proposed relaxation paths for the excited state corresponding to the most intense but short-lived absorption (interconversion, intersystem crossing, fluorescence, and phosphorescence) that may lead to longer-lived excited states necessary for efficient charge separation required to generate current in solar cells.

14 SOLAR ENERGY↗
DOE OSTI2023

Enhanced Exciton-to-Trion Conversion by Proton Irradiation of Atomically Thin WS 2

Defect engineering of van der Waals semiconductors has been demonstrated as an effective approach to manipulate the structural and functional characteristics toward dynamic device controls, yet correlations between physical properties with defect evolution remain underexplored. Here, using proton irradiation, we observe an enhanced exciton-to-trion conversion of the atomically thin WS 2 . The altered excitonic states are closely correlated with nanopore induced atomic displacement, W nanoclusters, and zigzag edge terminations, verified by scanning transmission electron microscopy, photoluminescence, and Raman spectroscopy. Density functional theory calculation suggests that nanopores facilitate formation of in-gap states that act as sinks for free electrons to couple with excitons. The ion energy loss simulation predicts a dominating electron ionization effect upon proton irradiation, providing further evidence on band perturbations and nanopore formation without destroying the overall crystallinity. This study provides a route in tuning the excitonic properties of van der Waals semiconductors using an irradiation-based defect engineering approach.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Thermal quasiparticle theory

The widely used thermal Hartree–Fock (HF) theory is generalized to include the effect of electron correlation while maintaining its quasi-independent-particle framework. An electron-correlated internal energy (or grand potential) is postulated in consultation with the second-order finite-temperature many-body perturbation theory (MBPT), which then dictates the corresponding thermal orbital (quasiparticle) energies in such a way that all fundamental thermodynamic relations are obeyed. The associated density matrix is of a one-electron type, whose diagonal elements take the form of the Fermi–Dirac distribution functions, when the grand potential is minimized. The formulas for the entropy and chemical potential are unchanged from those of Fermi–Dirac or thermal HF theory. The theory thus stipulates a finite-temperature extension of the second-order Dyson self-energy of one-particle many-body Green’s function theory and can be viewed as a second-order, diagonal, frequency-independent, thermal inverse Dyson equation. At low temperatures, the theory approaches finite-temperature MBPT of the same order, but it may outperform the latter at intermediate temperatures by including additional electron-correlation effects through orbital energies. Here, a physical meaning of these thermal orbital energies is proposed (encompassing that of thermal HF orbital energies, which has been elusive) as a finite-temperature version of Janak’s theorem.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Modeling separation of lanthanides via heterogeneous ligand binding

Individual lanthanide elements have physical/electronic/magnetic properties that make each useful for specific applications. Several of the lanthanides cations (Ln 3+ ) naturally occur together in the same ores. They are notoriously difficult to separate from each other due to their chemical similarity. Predicting the Ln 3+ differential binding energies (ΔΔE) or free energies (ΔΔG) at different binding sites, which are key figures of merit for separation applications, will help design of materials with lanthanide selectivity. We apply ab initio molecular dynamics (AIMD) simulations and density functional theory (DFT) to calculate ΔΔG for Ln 3+ coordinated to ligands in water and embedded in metal–organic frameworks (MOFs), and ΔΔE for Ln 3+ bonded to functionalized silica surfaces, thus circumventing the need for the computational costly absolute binding (free) energies ΔG and ΔE. Perturbative AIMD simulations of water-inundated simulation cells are applied to examine the selectivity of ligands towards adjacent Ln 3+ in the periodic table. Static DFT calculations with a full Ln 3+ first coordination shell, while less rigorous, show that all ligands examined with net negative charges are more selective towards the heavier lanthanides than a charge-neutral coordination shell made up of water molecules. Amine groups are predicted to be poor ligands for lanthanide-binding. Finally, we also address cooperative ion binding, i.e., using different ligands in concert to enhance lanthanide selectivity.

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

Large-scale galaxy flow from a non-gravitational impulse

A theory is presented describing linear perturbations of an expanding universe containing multiple, independently perturbed, collisionless, gravitationally coupled constituents. Solutions are found in the limit where one initially unperturbed component dominates the total density. The theory is applied to perturbations generated by a nongravitational process in one or more of the light components, as would occur in explosive or radiation-pressure-instability theories of galaxy formation. The apparent dynamical density parameter and correlations between density and velocity amplitude for various populations, are evaluated as a function of cosmic scale factor.

Hogan, Craig J.↗
DOE OSTI2021

Too big, too small, or just right? A benchmark assessment of density functional theory for predicting the spatial extent of the electron density of small chemical systems

Multipole moments are the first-order responses of the energy to spatial derivatives of the electric field strength. The quality of density functional theory prediction of molecular multipole moments thus characterizes errors in modeling the electron density itself, as well as the performance in describing molecules interacting with external electric fields. However, only the lowest non-zero moment is translationally invariant, making the higher-order moments origin-dependent. Therefore, instead of using the 3 × 3 quadrupole moment matrix, we utilize the translationally invariant 3 × 3 matrix of second cumulants (or spatial variances) of the electron density as the quantity of interest (denoted by K). The principal components of K are the square of the spatial extent of the electron density along each axis. A benchmark dataset of the principal components of K for 100 small molecules at the coupled cluster singles and doubles with perturbative triples at the complete basis set limit is developed, resulting in 213 independent K components. The performance of 47 popular and recent density functionals is assessed against this Var213 dataset. Several functionals, especially double hybrids, and also SCAN and SCAN0 predict reliable second cumulants, although some modern, empirically parameterized functionals yield more disappointing performance. The H, Li, and Be atoms, in particular, are challenging for nearly all methods, indicating that future functional development could benefit from the inclusion of their density information in training or testing protocols.

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

Second harmonic generation in topological insulators under quantizing magnetic fields

We theoretically investigate the second harmonic generation (SHG) of topological insulator surface states in a perpendicular magnetic field. Our theory is based on the microscopic expression of the second-order magneto-optical conductivity developed from the density matrix formalism, taking into account hexagonal warping effects on the surface states' band structure. Using numerically exact Landau-level energies and wave functions including hexagonal warping, we calculate the spectrum of SHG conductivities under normal incidence for different values of magnetic field and chemical potential. The imaginary parts of the SHG conductivities show prominent resonant peaks corresponding to one-photon and two-photon inter-Landau-level transitions. Treating the hexagonal warping term perturbatively, these transitions are clarified analytically within a perturbation theory from which approximate selection rules for the allowable optical transitions for SHG are determined. Furthermore, our results show extremely high SHG susceptibility that is easily tunable by magnetic field and doping level for topological surface states in the far-infrared regime, exceeding that of many conventional nonlinear materials. This work highlights the key role of hexagonal warping effects in generating second-order optical responses and provides insights into the nonlinear magneto-optical properties of the topological insulators.

Landau levels↗
DOE OSTI2022

Fermi–Löwdin orbital self-interaction correction of adsorption energies on transition metal ions

Density functional theory (DFT)-based descriptions of the adsorption of small molecules on transition metal ions are prone to self-interaction errors. Here, we show that such errors lead to a large over-estimation of adsorption energies of small molecules on Cu + , Zn + , Zn 2+ , and Mn + in local spin density approximation (LSDA) and Perdew, Burke, Ernzerhof (PBE) generalized gradient approximation calculations compared to reference values computed using the coupled-cluster with single, doubles, and perturbative triple excitations method. These errors are significantly reduced by removing self-interaction using the Perdew–Zunger self-interaction correction (PZ-SIC) in the Fermi–Löwdin Orbital (FLO) SIC framework. In the case of FLO-PBE, typical errors are reduced to less than 0.1 eV. Furthermore, analysis of the results using DFT energies evaluated on self-interaction-corrected densities [DFT(@FLO)] indicates that the density-driven contributions to the FLO-DFT adsorption energy corrections are roughly the same size in DFT = LSDA and PBE, but the total corrections due to removing self-interaction are larger in LSDA.

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

Excursion set mass functions for hierarchical Gaussian fluctuations

It is pointed out that most schemes for determining the mass function of virialized objects from the statistics of the initial density perturbation field suffer from the cloud-in-cloud problem of miscounting the number of low-mass clumps, many of which would have been subsumed into larger objects. The paper proposes a solution based on the theory of the excursion sets of F(r, R sub f), the four-dimensional initial density perturbation field smoothed with a continuous hierarchy of filters of radii R sub f.

Bond, J. R.↗
DOE OSTI2026

Influence of Markovianity and self-consistency on time-resolved spectral functions of driven quantum systems

We present a systematic comparison of the real-time Dyson expansion (RTDE) with established nonequilibrium Green's function (GF) approaches for simulating driven, interacting quantum systems. Focusing on density matrix dynamics, time-off-diagonal GFs, and time-resolved photoemission spectra, we benchmark RTDE against fully self-consistent Kadanoff-Baym equation (KBE) calculations, the generalized Kadanoff-Baym ansatz, and exact diagonalization for small systems using second-order many-body perturbation theory. Using a driven two-band Hubbard model, we show that mean-field single-particle density matrix trajectories provide a reliable baseline for RTDE across a broad range of interaction strengths and excited-carrier populations. Further, RTDE accurately captures correlation effects in the GFs, including long-lived oscillations and revivals that are strongly suppressed by the overdamping inherent to self-consistent KBE schemes. As a consequence, RTDE resolves rich nonequilibrium spectral structure in time-resolved photoemission, such as interaction- and population-dependent quasiparticle splittings and band gap renormalization, which are largely washed out in self-consistent approaches yet are present in exact solutions. Furthermore, our results demonstrate that RTDE bridges the gap between mean-field propagation and full two-time KBE simulations, retaining favorable linear scaling while capturing essential dynamical correlations relevant for ultrafast spectroscopy.

Electronic structure↗
DOE OSTI2021

Phonon-Assisted Ballistic Current from First-Principles Calculations

The bulk photovoltaic effect (BPVE) refers to current generation due to illumination by light in a homogeneous bulk material lacking inversion symmetry. In addition to the intensively studied shift current, the ballistic current, which originates from asymmetric carrier generation due to scattering processes, also constitutes an important contribution to the overall kinetic model of the BPVE. In this Letter, we use a perturbative approach to derive a formula for the ballistic current resulting from the intrinsic electron-phonon scattering in a form amenable to first-principles calculation. We then implement the theory and calculate the ballistic current of the prototypical BPVE material BaTiO 3 using quantum-mechanical density functional theory. The magnitude of the ballistic current is comparable to that of the shift current, and the total spectrum (shift plus ballistic) agrees well with the experimentally measured photocurrents. Furthermore, we show that the ballistic current is sensitive to structural change, which could benefit future photovoltaic materials design.

14 SOLAR ENERGY↗
DOE OSTI2021

Unimolecular and bimolecular formic acid decomposition routes on dispersed Cu nanoparticles

The elementary steps and site requirements in formic acid (HCOOH) dehydrogenation on Cu surfaces remain of keen interest because formate species act as intermediates or spectators in methanol synthesis and water–gas shift reactions. Steady-state and transient kinetic data, isotopic effects, infrared spectra during catalytic and stoichiometric reactions, and theoretical treatments based on density functional theory (DFT) provide evidence for bimolecular reactions, in which saturated bidentate formate (*HCOO*) adlayers, present at 0.25 ML (0.25 *HCOO* per surface Cu atom), react with undissociated species (HCOOH$^\square$) bound at interstices within formate adlayers ($^\square$) to form H-bonded bimolecular HCOOH$^\square$-*HCOO* adducts. The co-existence of vicinal HCOOH$^\square$ and *HCOO* moieties is evident from antisymmetric infrared bands for *HCOO* that become stronger as a result of their H-bonding that perturbs the induced dipole moment of *HCOO* upon vibration, consistent with DFT-derived vibrational frequencies and intensities for such perturbed species. The *HCOO* moiety in this complex undergoes C-H activation via a transition state that is preferentially stabilized through H-bonding with the vicinal HCOOH$^\square$ relative to its *HCOO* precursor. DFT-derived HCOOH dehydrogenation activation barriers and those determined from the evolution of CO 2 from pre-adsorbed *HCOO* species are about 10 kJ mol -1 smaller in the presence of gaseous HCOOH reactants (because of HCOOH$^\square$-*HCOO* interactions) than those for the unimolecular decomposition of bound *HCOO* species. Such bimolecular routes are consistent with measured effects of HCOOH, H 2 , and CO pressures and of H/D isotopic substitution on dehydrogenation turnover rates and represent the predominant channel for the formation of CO 2 and H 2 during catalytic HCOOH dehydrogenation on Cu nanoparticles. A saturated *HCOO* adlayer that retains binding interstices and the presence of HCOOH(g) enable a sequence of elementary steps unavailable for *HCOO* species, thus circumventing unassisted unimolecular routes that exhibit higher activation barriers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗