Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “generalized functions”

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

Longitudinal Surveillance for Chronic Health Conditions in Former United States Department of Energy Site Workers

The aim of the study was to determine (1) the rate at which rescreening former Department of Energy site workers identifies noncommunicable chronic diseases and (2) the development of comorbid conditions. Incidence and prevalence of hypertension, diabetes, reduced kidney function, and peripheral neuropathy at both initial and return screenings were calculated. Risk ratio of chronic disease development at return screening based on the presence of other conditions at initial screening were estimated with generalized linear regression. Prevalence of reduced kidney function was 19% at initial examination and 30% at return examination. The screening program was responsible for identifying 81% of reduced kidney function cases. Similar findings were present for the other chronic conditions examined. As a result, former worker health surveillance programs help identify significant health conditions among DOE workers, subcontractors, and visitors. Longitudinal screening of participants detects additional chronic conditions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Phase Equilibrium of Water with Hexagonal and Cubic Ice Using the SCAN Functional

Machine learning models are rapidly becoming widely used to simulate complex physicochemical phenomena with ab initio accuracy. Here, we use one such model as well as direct density functional theory (DFT) calculations to investigate the phase equilibrium of water, hexagonal ice (Ih), and cubic ice (Ic), with an eye toward studying ice nucleation. The machine learning model is based on deep neural networks and has been trained on DFT data obtained using the SCAN exchange and correlation functional. We use this model to drive enhanced sampling simulations aimed at calculating a number of complex properties that are out of reach of DFT-driven simulations and then employ an appropriate reweighting procedure to compute the corresponding properties for the SCAN functional. This approach allows us to calculate the melting temperature of both ice polymorphs, the driving force for nucleation, the heat of fusion, the densities at the melting temperature, the relative stability of ices Ih and Ic, and other properties. We find a correct qualitative prediction of all properties of interest. In some cases, quantitative agreement with experiment is better than for state-of-the-art semiempirical potentials for water. Furthermore, our results also show that SCAN correctly predicts that ice Ih is more stable than ice Ic.

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

Analytic and numerical vibronic spectra from quasi-classical trajectory ensembles

The truncated Wigner approximation to quantum dynamics in phase space is explored in the context of computing vibronic line shapes for monomer linear optical spectra. We consider multiple model potential forms including a shifted harmonic oscillator with both equal and unequal frequencies on the ground and excited state potentials as well as a shifted Morse potential model. For the equal-frequency shifted harmonic oscillator model, we derive an analytic expression for the exact vibronic line shape that emphasizes the importance of using a quantum mechanical distribution of phase space initial conditions. For the unequal-frequency shifted harmonic oscillator model, we are no longer able to obtain an exact expression for the vibronic line shape in terms of independent deterministic classical trajectories. We show how one can rigorously account for corrections to the truncated Wigner approximation through nonlinear responses of the line shape function to momentum fluctuations along a classical trajectory and demonstrate the qualitative improvement in the resulting spectrum when the leading-order quantum correction is included. Finally, we numerically simulate absorption spectra of a highly anharmonic shifted Morse potential model. We find that, while finite quantization and the dissociation limit are captured with reasonable accuracy, there is a qualitative breakdown of the quasi-classical trajectory ensemble’s ability to describe the vibronic line shape when the relative shift in Morse potentials becomes large. The work presented here provides clarity on the origin of unphysical negative features known to contaminate absorption spectra computed with quasi-classical trajectory ensembles.

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

Martensitic transformation in superlattices of two non-transforming metals

A number of studies have shown that the coherent integration of a second non-transforming phase in martensitic materials results in unprecedented thermo-mechanical properties. We demonstrate the possibility of martensitic transformation by epitaxially combining two non-transforming materials, Mg (hexagonal closed packed) and MgLi (body centered cubic). Density functional theory calculations predict that MgLi/Mg superlattices undergo a transformation crystallographically equivalent to MgSc alloys and can be expected to result in finite temperature shape memory behavior. We predict the zero-temperature energy of the martensite phase to be 9.3 meV/atom lower than austenite for a 50–50 at. % MgLi/Mg nanolaminate, which is expected to result in a martensitic transformation temperature of approximately 150 K. This energy difference can be tuned between 4.9 and 18.1 meV/atom by varying the fraction of the pure Mg phase in the superlattice, potentially enabling the transformation from cryogenic to room and higher temperatures. While Li diffusion might limit the applications of the systems proposed, this work opens the door to a new approach to create martensitic materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

AdS/BCFT from conformal bootstrap: construction of gravity with branes and particles

We initiate a conformal bootstrap program to study AdS 3 /BCFT 2 with heavy excitations. We start by solving the bootstrap equations associated with two-point functions of scalar/non-scalar primaries under the assumption that one-point functions vanish. These correspond to gravity with a brane and a non-spinning/spinning particle where the brane and the particle do not intersect with each other. From the bootstrap equations, we obtain the energy spectrum and the modified black hole threshold. We then carefully analyze the gravity duals and find the results perfectly match the BCFT analysis. In particular, brane self-intersections, which are usually considered to be problematic, are nicely avoided by the black hole formation. Despite the assumption to solve the bootstrap equations, one-point functions of scalar primaries can be non-zero in general. We construct the holographic dual for a non-vanishing one-point function, in which the heavy particle can end on the brane, by holographically computing the Rényi entropy in AdS/BCFT. As a bonus, we find a refined formula for the holographic Rényi entropy, which appears to be crucial to correctly reproduce the boundary entropy term. On the other hand, we explain why one-point functions of non-scalar primaries always vanish from the gravity dual. The non-sensitivity of the solution for the bootstrap equation to the boundary entropy helps us to construct gravity duals with negative tension branes. We also find a holographic dual of boundary primaries.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

From Anderson localization on random regular graphs to many-body localization

Highlights: • Anderson transition from ergodicity to localization on random regular graphs (RRG). • Analytical, pool method, and exact-diagonalization study of correlations on RRG. • Many-body localization (MBL) to ergodicity transition: quantum dots and spin chains. • Dynamical eigenstate correlation functions in RRG and MBL problems. • Anderson localization on RRG as a toy model for MBL. The article reviews the physics of Anderson localization on random regular graphs (RRG) and its connections to many-body localization (MBL) in disordered interacting systems. Properties of eigenstate and energy level correlations in delocalized and localized phases, as well at criticality, are discussed. In the many-body part, models with short-range and power-law interactions are considered, as well as the quantum-dot model representing the limit of the “most long-range” interaction. Central themes – which are common to the RRG and MBL problems – include ergodicity of the delocalized phase, localized character of the critical point, strong finite-size effects, and fractal scaling of eigenstate correlations in the localized phase.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Molecular Modes of Attosecond Charge Migration

Here, first-principles calculations are employed to elucidate the modes of attosecond charge migration (CM) in halogenated hydrocarbon chains. We use constrained density functional theory (DFT) to emulate the creation of a localized hole on the halogen and follow the subsequent dynamics via time-dependent DFT. We find low-frequency CM modes (~ 1 eV) that propagate across the molecule and study their dependence on length, bond order, and halogenation. We observe that the CM speed (~ 4Å/ fs) is largely independent of molecule length, but is lower for triple-bonded versus double-bonded molecules. Additionally, as the halogen mass increases, the hole travels in a more particlelike manner as it moves across the molecule. These heuristics will be useful in identifying molecules and optimal CM detection methods for future experiments, especially for halogenated hydrocarbons which are promising targets for ionization-triggered CM.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Improvements in the orbitalwise scaling down of Perdew–Zunger self-interaction correction in many-electron regions

The Perdew-Zunger (PZ) method offers a way to remove the self-interaction (SI) error from density functional approximations on an orbital by orbital basis. The PZ method provides significant improvements for the properties such as barrier heights or dissociation energies but results in over-correcting the properties well described by SI-uncorrected semi-local functional. One cure to rectify the over-correcting tendency is to scale down the magnitude of SI-correction of each orbital in the many electron region. We have implemented the orbitalwise scaled down SI-correction (OSIC) scheme of Vydrov et al. [J. Chem. Phys. 124, 094108 (2006)] using the Fermi-L owdin SI-correction method. After validating the OSIC implementation with previously reported OSICLSDA results, we examine its performance with the most successful non-empirical SCAN meta-GGA functional. Using different forms of scaling factors to identify one-electron regions, we assess the performance of OSIC-SCAN for a wide range of properties: total energies, ionization potentials and electron affinities for atoms, atomization energies, dissociation and reaction energies, and reaction barrier heights of molecules. Our results show that OSIC-SCAN provides superior results than the previously reported OSIC-LSDA, -PBE, and -TPSS results. Furthermore, we propose selective scaling of OSIC (SOSIC) to remove its major shortcoming that destroys the -1/r asymptotic behavior of the potentials. The SOSIC method gives the highest occupied orbital eigenvalues practically identical to those in PZSIC and unlike OSIC provides bound atomic anions even with larger powers of scaling factors. Furthermore, SOSIC compared to PZSIC or OSIC provides more balanced description of total energies and barrier heights.

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

Pressure-induced structural and dielectric changes in liquid water at room temperature

Understanding the pressure-dependent dielectric properties of water is crucial for a wide range of scientific and practical applications. In this study, we employ a deep neural network trained on density functional theory data to investigate the dielectric properties of liquid water at room temperature across a pressure range of 0.1–1000 MPa. We observe a nonlinear increase in the static dielectric constant ɛ 0 with increasing pressure, a trend that is qualitatively consistent with experimental observations. This increase in ɛ 0 is primarily attributed to the increase in water density under compression, which enhances collective dipole fluctuations within the hydrogen-bonding network as well as the dielectric response. Despite the increase in ɛ 0 , our results reveal a decrease in the Kirkwood correlation factor G K with increasing pressure. Furthermore, this decrease in G K is attributed to pressure-induced structural distortions in the hydrogen-bonding network, which weaken dipolar correlations by disrupting the ideal tetrahedral arrangement of water molecules.

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

Transient Triamidoamine Neptunium(V)–Mono(Imido) Complexes: C–H Activations and Hydrogen Atom Transfer Driven by Effective Nuclear Charge

Metal-mono(imido) linkages have been known for seven decades, and they are found in transition metal, main group, lanthanide, thorium, and uranium complexes. However, transuranium-mono(imido) complexes remain unknown in any scenario. Here, we present evidence for transient neptunium(V)–mono(imido) complexes. Treatment of [Np III (Tren TIPS )] (1, Tren TIPS = {N(CH 2 CH 2 NSiPr i 3 ) 3 } 3– ) with N 3 R (R = SiMe 3 ; 1-adamantyl, Ad) results in N 2 evolution and dark purple solutions consistent with the formation of [Np V (Tren TIPS )(NR)] (3NpNR). However, solutions of 3NpNR rapidly turn orange, where for R = SiMe 3 the isolated 1:1 products are [Np IV (Tren TIPS ){N(H)SiMe 3 }] (4a) and [Np IV (Tren TIPS-2H ){N(H)SiMe 3 }] (4b, Tren TIPS-2H = {N(CH 2 CH 2 NS i Pri 3 ) 2 (NCH 2 CH 2 NSiPr i 2 C[Me]=CH 2 )} 3– ). The latter contains a dehydrogenated-Pr i vinyl functionality accounting for the source of the two amido H atoms. The reaction for R = Ad proceeds similarly, but only [Np IV (Tren TIPS ){N(H)Ad}] (5a) could be unequivocally confirmed, though its isolation suggests generality of the imido-to-amido functional group transformation. Complexes 4a/4b exhibit slow relaxation of their magnetization, adding to the small number of transuranium single ion magnets. Experimental and computational analysis suggests that the amido products are formed by C–H activation and two sequential hydrogen atom transfer reactions involving a three-step proton-coupled electron-transfer sequence of H • radical abstraction, electron transfer, then another H • radical abstraction step. In contrast to transient 3NpNR, the 5f 2 uranium(IV)-imido complex [K(2.2.2-cryptand)][U IV (Tren TIPS )(NSiMe 3 )] (8UNSiMe 3 ) is robust, even in boiling THF, suggesting the transience of 5f 2 3NpNR is not due to the 5f n -count but the increased effective nuclear charge of neptunium vs uranium. This work highlights divergence of uranium- and neptunium-imido stabilities, emphasizing that the latter is an inherently challenging synthetic target.

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

Sensitivity of the electronic and magnetic structures of cuprate superconductors to density functional approximations

Abstract We discuss the crystal, electronic, and magnetic structures of La 2− x Sr x CuO 4 (LSCO) for x = 0.0 and x = 0.25 employing 13 density functional approximations, representing the local, semi-local, and hybrid exchange-correlation approximations within the Perdew–Schmidt hierarchy. The meta-generalized gradient approximation (meta-GGA) class of functionals is found to perform well in capturing the key properties of LSCO, a prototypical high-temperature cuprate superconductor. In contrast, the localspin-density approximation, GGA, and the hybrid density functional fail to capture the metal-insulator transition under doping.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Nonadiabatic couplings from a variational excited state method based on constrained DFT

Excited Costrained Density Functional Theory (XCDFT) is a variational excited state method that extends ground state DFT to the computation of low-lying excited states. It borrows much of the machinery of Constrained DFT (CDFT) with a crucial difference: the constraint imposes a population of one electron in the Hilbert space spanned by the virtuals of a reference ground state. Here, we present theory and implementation for evaluating nonadiabatic coupling vectors (NACVs) between the first excited state computed with XCDFT and the ground state. Our NACVs are computed analytically using density functional perturbation theory with a formalism that is general enough that could be applied to CDFT diabatic states. We showcase the new method with pilot NACV calculations for the conical intersection in H 3 , the avoided crossing in selenoacrolein, and the NACV magnitudes in azobenzene. Despite complications from the nonorthogonality of the wavefunctions, XCDFT’s energy surfaces and NACVs reproduce benchmark values and respect known sum rules within a reasonable degree. This shows that XCDFT is a viable method for nonadiabatic dynamics simulations.

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

Constraining radio mode feedback in galaxy clusters with the cluster radio AGNs properties to z ∼ 1

ABSTRACT We study the properties of the Sydney University Molonglo Sky Survey (SUMSS) 843 MHz radio active galactic nuclei (AGNs) population in galaxy clusters from two large catalogues created using the Dark Energy Survey (DES): ∼11 800 optically selected RM-Y3 and ∼1000 X-ray selected MARD-Y3 clusters. We show that cluster radio loud AGNs are highly concentrated around cluster centres to $z$ ∼ 1. We measure the halo occupation number for cluster radio AGNs above a threshold luminosity, finding that the number of radio AGNs per cluster increases with cluster halo mass as N ∝ M1.2 ± 0.1 (N ∝ M0.68 ± 0.34) for the RM-Y3 (MARD-Y3) sample. Together, these results indicate that radio mode feedback is favoured in more massive galaxy clusters. Using optical counterparts for these sources, we demonstrate weak redshift evolution in the host broad-band colours and the radio luminosity at fixed host galaxy stellar mass. We use the redshift evolution in radio luminosity to break the degeneracy between density and luminosity evolution scenarios in the redshift trend of the radio AGNs luminosity function (LF). The LF exhibits a redshift trend of the form (1 + $z$)γ in density and luminosity, respectively, of γD = 3.0 ± 0.4 and γP = 0.21 ± 0.15 in the RM-Y3 sample, and γD = 2.6 ± 0.7 and γP = 0.31 ± 0.15 in MARD-Y3. We discuss the physical drivers of radio mode feedback in cluster AGNs, and we use the cluster radio galaxy LF to estimate the average radio-mode feedback energy as a function of cluster mass and redshift and compare it to the core (<0.1R500) X-ray radiative losses for clusters at $z$ < 1.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

First-principles wave-vector- and frequency-dependent exchange-correlation kernel for jellium at all densities

Here we propose a spatially and temporally nonlocal exchange correlation (XC) kernel for the spin-unpolarized fluid phase of ground-state jellium for use in time-dependent density functional and linear response calculations. The kernel is constructed to satisfy known properties of the exact XC kernel to accurately describe the correlation energies of bulk jellium and to satisfy frequency-moment sum rules at a wide range of bulk jellium densities, including those low densities that display strong correlation and symmetry breaking. These effects are easier to understand in the simple jellium model than in real systems. All exact constraints satisfied by the recent MCP07 kernel are maintained in the revised MCP07 (rMCP07) kernel, while others are added. The revision $f^{rMCP07}_{XC}$ (q, ω) differs from MCP07 only for nonzero frequencies ω. Only at densities much lower than those of real bulk metals is the frequency dependence of the kernel important for the correlation energy of jellium. As the wave vector q tends to zero, the kernel has a -4$πα(ω)/q^2$ divergence whose frequency-dependent ultranonlocality coefficient $α(ω)$ vanishes in jellium, and is predicted by rMCP07 to be extremely small for the real metals Al and Na.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Chemically and electrically tunable spin polarization in ferroelectric Cd-based hybrid organic-inorganic perovskites

Density functional theory computations are used to predict the electronic structure of Cd-based hybrid organic-inorganic high-T C ferroelectric perovskites with TMCM-CdCl 3 being one representative. Here we report Rashba-Dresselhaus spin splitting in the valence band of these nonmagnetic compounds. Interestingly, we find in computations that the splitting is not necessarily sensitive to the polarization of the material but to the organic molecule itself which opens a way to its chemical tunability through the choice of the molecule. Further chemical tunability of spin splitting is shown to be possible through a substitution of Cl in the CdCl 3 chains as the valence band was found to originate from Cl-Cl weekly bonding orbitals. For example, the substitution of Cl with Br in TMCM-CdCl 3 resulted in a ten times increase of spin splitting. Furthermore, the spin polarization in these materials gives origin to persistent spin textures which are coupled to the polarization direction, and, therefore, can be controlled by the electric field. This is promising for spintronics applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Systematic Evaluation of Counterpoise Correction in Density Functional Theory

A widespread belief persists that the Boys–Bernardi function counterpoise (CP) procedure “overcorrects” supramolecular interaction energies for the effects of basis-set superposition error. To the extent that this is true for correlated wave function methods, it is usually an artifact of low-quality basis sets. The question has not been considered systematically in the context of density functional theory, however, where basis-set convergence is generally less problematic. We present a systematic assessment of the CP procedure for a representative set of functionals and basis sets, considering both benchmark data sets of small dimers and larger supramolecular complexes. The latter include layered composite polymers with ~150 atoms and ligand–protein models with ~300 atoms. Provided that CP correction is used, we find that intermolecular interaction energies of nearly complete-basis quality can be obtained using only double-ζ basis sets. Furthermore, this is less expensive as compared to triple-ζ basis sets without CP correction. CP-corrected interaction energies are less sensitive to the presence of diffuse basis functions as compared to uncorrected energies, which is important because diffuse functions are expensive and often numerically problematic for large systems. Our results upend the conventional wisdom that CP “overcorrects” for basis-set incompleteness. In small basis sets, CP correction is mandatory in order to demonstrate that the results do not rest on error cancellation.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Asymmetric equilibrium core structures of pyramidal-II < c + a > dislocations in ten hexagonal-close-packed metals

The structures of pyramidal-II < c + a > dislocations, one of the most important defects in structural hexagonal-close-packed (HCP) metals, have not been fully characterized for many of the HCP metals in use today. Here, we employ ab initio informed phase-field dislocation dynamics to determine the minimum energy structure of pyramidal {1¯1¯22} < 11¯23 > dislocations in ten HCP metals, including Be, Co, Mg, Re, Ti, Zn, Cd, Hf, Y, and Zr. As input for the simulations, we calculate, using first-principles density functional theory, the {1¯1¯22} generalized stacking fault energy (GSFE) curves for all ten metals. From these calculations, it is found that magnetism in Co is necessary for achieving a local minimum in the GSFE curve. We observe in simulations that edge and screw character dislocations split into two partials separated by a low-energy intrinsic stacking fault. The splitting distance is shown to scale inversely with the local minimum energy normalized by the product of its shear modulus and Burgers vector. Interestingly, some HCP metals exhibit an asymmetric structure, with either unequal partial Burgers vectors or widths, in contrast to the symmetric configuration expected from linear elastic dislocation theory. We explain these structures by properties of the local maxima in their GSFE curves. Metals with larger degrees of elastic anisotropy result in dislocations with larger splitting distances than would be expected under the commonly used assumption of elastic isotropy. Furthermore, these findings on the sizes and asymmetry in the structures of pyramidal-II < c + a > dislocations are fundamental to understanding how these dislocations glide and interact or react with other defects when these metals are mechanically strained.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Spatial correlation functions of paracrystals with radial symmetry

Here, we develop a phenomenological model to describe the structure of radially symmetric paracrystals whose long-range order are destroyed by propagation of particle fluctuations. General expressions are derived for the spatial correlation functions in one-, two-, and three-dimensional spaces. And the spatial correlation in paracrystals in reciprocal space is further discussed and clarified. The developed method can be used to quantitatively analyze the microstructure of paracrystalline materials in both real and reciprocal spaces via scattering experiments and computer simulations.

36 MATERIALS SCIENCE↗