Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “energy density functional”

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 1,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Unoccupied electronic structure of actinide dioxides

Ligand field density functional theory calculations of the dioxides of thorium, uranium, and plutonium have been combined with high-energy-resolution fluorescence detection (HERFD) in x-ray absorption spectroscopy and inverse photoelectron spectroscopy (IPES) measurements to provide powerful insight into the underlying composition of the unoccupied 5f electronic structure in these 5f localized systems. Fine structure in the 5f 5/2 transitions in HERFD can be directly correlated with the fine structure in the leading edge of the IPES. The shapes, intensities, and systematics in HERFD and IPES are explained in a consistent and rigorous fashion in terms of the j-specific 5f electronic structure.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Putting error bars on density functional theory dataset

This dataset contains submission files and raw output files from high-throughput DFT simulations to analyze the systemic errors in lattice constant, bulk moduli and formation energy predictions for a range of binary and ternary oxides using four exchange correlation functionals (LDA, PBE, PBEsol and vdW-DF-C09). This data was then used as the basis for employing materials informatics methods to predict the expected errors in the lattice constants of the studied compounds. Predicted errors were also used to better the DFT-predicted lattice parameters. Our results emphasize the link between the computed errors and the electron density and hybridization errors of a functional. In essence, these results provide “error bars” for choosing a functional for the creation of high-accuracy, high-throughput datasets as well as avenues for the development of XC functionals with enhanced performance, thereby enabling the accelerated discovery and design of new materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Nitrogen-doped char as a catalyst for wet oxidation of phenol-contaminated water

Catalytic wet oxidation (CWO) of aqueous effluents rich in organic compounds is a very promising technology for the treatment of liquid wastes from biomass conversion processes. CWO reactions occur through the formation of free radical species, produced in the presence of an oxidant, which act on organic contaminates in the effluent. Although the reaction is well known, there exists a lack of affordable catalysts to conduct this process at the lower temperatures and pressures in novel bioenergy processes. Here, this study assessed the catalytic effect of nitrogen-doped chars as such an option. Phenol in aqueous solution was used as a model waste effluent. Treatment was conducted at moderate temperatures (190 to 260 °C), oxygen partial pressure of 1 MPa, and reaction times of 15, 30, and 45 min in stainless steel and glass-lined tube reactors. High pressure liquid chromatography (HPLC) analyses of the products quantified phenol and by-product concentrations used in the calculation of reaction activation energy. The char catalyst was studied by X-ray photoelectron spectroscopy (XPS), transmission electron microscopy (TEM), and scanning electron microscopy (SEM) in order to gain insight into its structure and surface composition. The results indicate that nitrogen-doped char catalysts accelerate the oxidation of phenol by decreasing its reaction activation energy from 82.2 kJ/mol (non-catalyzed) to 40.4 kJ/mol (catalyzed). An analysis from first principles using density functional theory (DFT) was conducted to ascertain which N functional group has the most significant impact on free radical formation in the presence of oxygen. Among all the N functional groups studied, the dipyridinic functional groups showed the most promising characteristics to facilitate the formation of hydroxyl free radicals.

Advanced Oxidation Processes for Organic Pollutant↗
DOE OSTI2025

Single-crystalline orthorhombic GdAlGe as a rare-earth magnetic Dirac nodal-line metal

Crystal engineering is a method for discovering new quantum materials and phases, which may be achieved using external pressure or strain. Chemical pressure is unique in that it generates internal pressure perpetually to the lattice. As an example, GdAlSi from the rare-earth (𝑅) 𝑅⁢Al⁢𝑋 (𝑋=Si or Ge) family of Weyl semimetals is considered. Replacing Si with the larger isovalent element Ge creates sufficiently large chemical pressure to induce a structural transition from the tetragonal structure of GdAlSi, compatible with a Weyl semimetallic state, to an orthorhombic phase in GdAlGe, resulting in an inversion-symmetry-protected nodal-line metal. We find that GdAlGe hosts an antiferromagnetic ground state with two successive orderings, at 𝑇 N ⁢1=35K and 𝑇 N ⁢2=30K. In-plane isothermal magnetization shows a magnetic field induced metamagnetic transition at 6.2 T for 2 K. Furthermore, electron-hole compensation gives rise to a large magnetoresistance of ∼100% at 2 K and 14 T. Angle-resolved photoemission spectroscopy measurements and density functional theory calculations reveal a Dirac-like linear band dispersion over an exceptionally large energy range of ∼1.5eV with a high Fermi velocity of ∼10 6 m/s, a rare feature not observed in any magnetic topological materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Machine learning pipeline to predict defect behavior in metallic alloy systems

The interaction between defect and solute atoms is critical to the thermodynamic and kinetic behavior of metallic alloys under exposure to high-energy radiation, causing irradiation damage in materials. Radiation can generate non-equilibrium concentrations of point defects such as vacancies and interstitials. The excess point defects not only accelerate diffusional processes such as precipitation that cause radiation embrittlement, but also change the pathway of phase transformations, including nucleation processes. Understanding these defect behaviors is complicated by the challenge and complexity of addressing each possible local and discrete distribution of environments and chemical interactions around targeted defects-solute or solute-solute complexes. To resolve the challenge, machine learning regression techniques have emerged as powerful tools that can train and construct an energy model to accurately describe the chemical interactions of solutes and defects. In Fiscal Year 2022, the work focused on the workflow development and demonstration using machine learning regression, density functional theory, cluster expansion, and Monte Carlo simulation to predict the effects of ternary solute elements (e.g., aluminum and molybdenum) and point defects on the Cr-rich $\alpha^{\prime}$ precipitation in multicomponent FeCr model alloys. The computational outcomes include the prediction of the ternary phase diagram, vacancy formation energy for different compositions, and the effect of vacancies on the nucleation of Cr-rich clusters. The simulations predict a pronounced change of Cr solubility in bcc Fe by the addition of Al and the rejection of Al atoms from $\alpha^{\prime}$ precipitates. Additionally, the simulations show the formation of Cr-vacancy clusters as the initial nuclei for stable nucleation and growth of $\alpha^{\prime}$ particles. The results demonstrate important outcomes and applications of using machine learning pipeline to study model or commercial alloys with multicomponent solute species and point defects.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Coupling between Emissive Defects on Carbon Nanotubes: Modeling Insights

Covalent functionalization of single-walled carbon nanotubes (SWCNTs) with organic molecules results in red-shifted emissive states associated with sp 3 -defects in the tube lattice, which facilitate their improved optical functionality, including single-photon emission. The energy of the defect-based electronic excitations (excitons) depends on the molecular adducts, the configuration of the defect, and concentration of defects. In this work, we model the interactions between two sp 3 -defects placed at various distances in the (6,5) SWCNT using time-dependent density functional theory. Calculations reveal that these interactions conform to the effective model of J-aggregates for well-spaced defects (>2 nm), leading to a red-shifted and optically allowed (bright) lowest energy exciton. H-aggregate behavior is not observed for any defect orientations, which is beneficial for emission. The splitting between the lowest energy bright and optically forbidden (dark) excitons and the pristine excitonic band are controlled by the single-defect configurations and their axial separation. These findings enable a synthetic design strategy for SWCNTs with tunable near-infrared emission.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Correlation between Spin and Orbital Dynamics during Laser-Induced Femtosecond Demagnetization

Spin and orbital angular momenta are two intrinsic properties of an electron and are responsible for the physics of a solid. How the spin and orbital evolve with respect to each other on several hundred femtoseconds is largely unknown, but it is at the center of laser-induced ultrafast demagnetization. In this paper, we introduce a concept of the spin–orbital correlation diagram, where spin angular momentum is plotted against orbital angular momentum, much like the position-velocity phase diagram in classical mechanics. We use four sets of highly accurate time-resolved X-ray magnetic circular dichroism data to construct four correlation diagrams for iron and cobalt. To our surprise, a pattern emerges. The trace on the correlation diagram for iron is an arc, and at the end of demagnetization, it has a pronounced cusp. The correlation diagram for cobalt is different and appears more linear but with kinks. We carry out first-principles calculations with two different methods: time-dependent density functional theory (TDDFT) and time-dependent Liouville density functional theory. These two methods agree that the experimental findings for both Fe and Co are not due to experimental errors. It is the spin–orbit coupling that correlates the spin dynamics to the orbital dynamics. Microscopically, Fe and Co have different orbital occupations, which leads to distinctive correlation diagrams. We believe that this correlation diagram presents a useful tool to better understand spin and orbital dynamics on an ultrafast time scale. A brief discussion on the magnetic anisotropy energy is also provided.

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

A Gaussian field approach to the solvation of spherical ions in electrolyte solutions

In this paper, the electrostatic response of an electrolyte solution to a spherical ion is studied with a Gaussian field theory. In order to capture the ionic correlation effect in concentrated solutions, the bulk dielectric response function is described by a two-Yukawa response function. The modified response function of the solution is solved analytically in the spherical geometry, from which the induced charge density and the electrostatic energy are also derived analytically. Comparisons with results for small ions in electrolyte solutions from the hyper-netted chain theory demonstrate the validity of the Gaussian field theory.

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

Structure and Site Evolution of Framework Ni Species in MIL-127 MOFs for Propylene Oligomerization Catalysis

A mixed-valence oxotrimer metal–organic framework (MOF), Ni-MIL-127, with a fully coordinated nickel atom and two iron atoms in the inorganic node, generates a missing linker defect upon thermal treatment in helium (>473 K) to engender an open coordination site on nickel which catalyzes propylene oligomerization devoid of any cocatalysts or initiators. This catalyst is stable for ~20 h on stream at 500 kPa and 473 K, unprecedented for this chemistry. The number of missing linkers on synthesized and activated Ni-MIL-127 MOFs is quantified using temperature-programmed oxidation, 1 H nuclear magnetic resonance spectroscopy, and X-ray absorption spectroscopy to be ~0.7 missing linkers per nickel; thus, a majority of Ni species in the MOF framework catalyze propylene oligomerization. In situ NO titrations under reaction conditions enumerate ~62% of the nickel atoms as catalytically relevant to validate the defect density upon thermal treatment. Propylene oligomerization rates on Ni-MIL-127 measured at steady state have activation energies of 55–67 kJ mol –1 from 448 to 493 K and are first-order in propylene pressures from 5 to 550 kPa. Density functional theory calculations on cluster models of Ni-MIL-127 are employed to validate the plausibility of the missing linker defect and the Cossee–Arlman mechanism for propylene oligomerization through comparisons between apparent activation energies from steady-state kinetics and computation. Here this study illustrates how MOF precatalysts engender defective Ni species which exhibit reactivity and stability characteristics that are distinct and can be engineered to improve catalytic activity for olefin oligomerization.

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

Properties of La 0.7 Ca 0.3 MnO 3 under extreme tensile strain

The complex phase diagram of manganites with simultaneously active spin, charge, orbital, and lattice degrees of freedom continues providing surprises. In a recent groundbreaking experiment, membranes of the perovskite manganite La 0.7 Ca 0.3 MnO 3 (LCMO) deposited on a flexible polymer layer were strained up to 8% [S. S. Hong et al., Science 368, 71 (2020)], much more than achieved by regular strain induced by a rigid substrate. By increasing this strain, a metal-insulator transition was reported. Here we reproduce the results of the experiments using Monte Carlo simulations of the two-orbital double-exchange model including Jahn-Teller distortions at hole density x=1/3. The full phase diagram with varying temperature and Jahn-Teller coupling λ is presented. When the bandwidth W of mobile electrons is reduced, thus when the effective Jahn-Teller coupling λ/W is increased, a metal-insulator transition is found in our simulations, between a ferromagnetic metallic state with uniform charge distribution and an insulator with diagonal charge stripes that retains its ferromagnetic character. In between the hole-rich diagonals, staggered orbital order occurs. We also report resistivity and magnetization measurements alongside with spin correlations and charge structure factors. Our overall conclusions are in agreement with the recent experimental and density functional theory results by Hong et al., and we confirm much earlier ground state predictions of striped ferromagnetic order using energy optimization techniques by T. Hotta et al. [Phys. Rev. Lett. 86, 4922 (2001)]. The experimental observation of one of the states predicted by theory suggests that diagonal stripes could be achieved at other hole densities as well, such as x=1/4, if LCMO membranes with that hole doping were subject to similar strains.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Free spectator nucleons in ultracentral relativistic heavy-ion collisions as a probe of neutron skin

Besides the yield ratio of free spectator neutrons produced in ultracentral 96 Zr + 96 Zr to 96 Ru + 96 Ru collisions, we propose that the yield ratio N n /N p of free spectator neutrons to protons in a single collision system at the BNL Relativistic Heavy Ion Collider or the CERN Large Hadron Collider can be a more sensitive probe of the neutron-skin thickness Δrnp and the slope parameter L of the symmetry energy. Here, the idea is demonstrated based on the proton and neutron density distributions of colliding nuclei obtained from Skyrme-Hartree-Fock-Bogolyubov calculations, and on a Glauber model that provides information of spectator matter. The final spectator particles are produced from direct emission, clusterization by a minimum spanning tree algorithm or a Wigner function approach, and deexcitation of heavy clusters by gemini. A larger Δr np associated with a larger L value increases the isospin asymmetry of spectator matter and thus leads to a larger N n /N p , especially in ultracentral collisions where the multiplicity of free nucleons are free from the uncertainties of cluster formation and deexcitation. We have further shown that the double ratio of N n /N p in isobaric collision systems or in collisions by isotopes helps to cancel the detecting efficiency for protons. Effects from nuclear deformation and electromagnetic excitation are studied, and they are found to be subdominant compared to the expected sensitivity to Δr np .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Enhanced stability and ordering of 2D MnBi2Se2Te2

This dataset includes density functional theory and quantum Monte Carlo calculations on MnBi(Se(1-x)Te(x))4 system for its phase stability and defect formation energies. The data supports publication "Increased Defect Resistance and Ordering in MnBi2(Se(1–x)Te(x))4 via Accurate Diffusion Monte Carlo" (https://doi.org/10.1021/acs.jpcc.5c04299).

36 MATERIALS SCIENCE↗
DOE OSTI2021

Catalyst Design in Nitrate Removal

Based on the volcano plot developed by Dr. Goldsmith group (Report linked in submission), we utilized DFT (density functional theory) calculations to search for bimetallic materials in the application of catalysts in aqueous nitrate removal. The calculations are conducted via the high-throughput automated workflow package developed by our group (Github linked in submission) using VASP commercial first-principles calculation software.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Static Subspace Approximation for Random Phase Approximation Correlation Energies: Implementation and Performance

Developing theoretical understanding of complex reactions and processes at interfaces requires using methods that go beyond semilocal density functional theory to accurately describe the interactions between solvent, reactants and substrates. Methods based on many-body perturbation theory, such as the random phase approximation (RPA), have previously been limited due to their computational complexity. However, this is now a surmountable barrier due to the advances in computational power available, in particular through modern GPU-based supercomputers. In this work, we describe the implementation of RPA calculations within BerkeleyGW and show its favorable computational performance on large complex systems relevant for catalysis and electrochemistry applications. Our implementation builds off of the static subspace approximation which, by employing a compressed representation of the frequency dependent polarizability, enables the evaluation of the RPA correlation energy with significant acceleration and systematically controllable accuracy. We find that the computational cost of calculating the RPA correlation energy scales only linearly with system size for systems containing up to 50 thousand bands, and is expected to scale quadratically thereafter. We also show excellent strong scaling results across several supercomputers, demonstrating the performance and portability of this implementation.

algorithmic development↗
DOE OSTI2023

Marginal stability constraint on runaway electron distribution

High-frequency kinetic instabilities of the strongly anisotropic runaway electrons (RE) can enhance the pitch-angle scattering of the runaways significantly. This wave-induced scattering can easily prevail over runaway scattering on high-Z impurities. In a steady state, collisional damping balances the kinetic drive of the unstable waves, keeping the RE distribution function at marginal stability. The marginal stability constraint limits the achievable RE densities and the shape of the RE distribution function. In this study, we consider whistler and compressional Alfvén waves as the primary source of enhanced elastic scattering of the runaways. By balancing the anomalous Doppler resonance drive with the collisional wave damping, we find the RE distribution function in the ultra-relativistic range of the phase space. We also derive an expression for the spectral energy density of the waves. We show that the power needed to compensate for the wave dissipation is negligible compared to the work of the electric field. The latter is in balance with the synchrotron losses of the runaway electrons.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Fullerene Negative Ions: Formation and Catalysis

We first explore negative-ion formation in fullerenes C 44 to C 136 through low-energy electron elastic scattering total cross sections calculations using our Regge-pole methodology. Then, the formed negative ions C 44 – to C 136 – are used to investigate the catalysis of water oxidation to peroxide and water synthesis from H 2 and O 2 . The exploited fundamental mechanism underlying negative-ion catalysis involves hydrogen bond strength-weakening/breaking in the transition state. Density Functional Theory transition state calculations found C 60 – optimal for both water and peroxide synthesis, C 100 – increases the energy barrier the most, and C 136 – the most effective catalyst in both water synthesis and oxidation to H 2 O 2 .

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

First-principles study of the surface properties of uranium carbides

Uranium carbides have attracted renewed interest as advanced nuclear fuels for Generation IV reactors. As an important property required for gas bubble modeling in nuclear fuels, the surface energy of uranium carbides is scarce in literature. In this work, we study the surface properties of uranium carbides by first-principles density functional theory calculations. Surface orientations with maximum Miller index up to 3, 2 and 2 are investigated for UC, U 2 C 3 and α-UC 2 , respectively. By studying the effects of surface termination and chemical potential on surface energy, we identify the factors that determines the surface stability. From the calculated surface energies, the surface properties of uranium carbide single crystals are obtained from Wulff construction, including equilibrium morphology, dominant surface orientation and area weighted surface energy.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Solid-State Hydroxide Ion Conductivity in Silver(I) Oxide, Ag 2 O

Silver (I) oxide, Ag 2 O, precipitated as microcrystals by combining aqueous silver(I) nitrate and KOH solutions, was found to be a solid-state hydroxide ion conductor with ionic conductivity on the order of 10 -3 S/cm. The proton chemical shifts at 4.87 and -7.35 ppm measured by solid-state 1H-NMR experiments are attributed to water molecules and in-lattice OH - coordinated to silver, respectively. The lack of spinning sidebands around the 4.87 ppm peak indicates rapid reorientation on the NMR timescale, suggesting that the water molecules are adsorbed to the surface of the Ag 2 O crystals. Pulsed field gradient measurements gave similar diffusion coefficients (2 x 10 -7 cm 2 /s at 298K) for all three proton environments indicating chemical exchange between sites on the millisecond timescale. The activation energy for OH - diffusion measured by NMR (0.18 eV) was comparable to that obtained by conductivity measurements and density functional theory (DFT) electronic structure calculations. Furthermore, the calculated Pourbaix diagram of Ag 2 O is consistent with the slightly lower sample density observed in He pycnometry and thermogravimetric measurements.

25 ENERGY STORAGE↗