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,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Helium impurities and interactions in lithium

Here we investigate helium interactions with lithium using density functional theory. Like other body-centered cubic (bcc) metals, the lowest-energy site for interstitial helium is a tetrahedral site. However, helium in lithium shows a higher substitutional formation energy than either the tetrahedral or octahedral interstitial formation energies, which is unique. The calculated migration energies of helium in lithium are also very low ( ~ 1 meV), an order of magnitude lower than those in other bcc metals. We find an increase in the binding energy to a vacancy as the number of helium atoms bound to it increases. The extremely low migration energies suggest that helium transport in lithium will be very fast.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Strong Magnetocrystalline Anisotropy Arising from Metal–Ligand Covalency in a Metal–Organic Candidate for 2D Magnetic Order

Layered metal–organic frameworks are promising candidates for new two-dimensional magnets, as their synthetic programmability of these materials can provide a route to diverse structural and electronic properties. However, such framework materials typically lack the heavy elements that engender magnetocrystalline anisotropy in the monolayer ferromagnets reported to date. Alternative sources of magnetic anisotropy are therefore needed in these materials. Here, we report the synthesis of single crystals of the framework material (NMe 4 ) 2 [Fe 2 L 3 ] (H 2 L = 3,6-dichloro-2,5-dihydroxybenzoquinone) and evaluate the angular dependence of its magnetic properties. Oriented-crystal magnetization measurements reveal strong uniaxial anisotropy, where the easy axis is aligned with the crystallographic c axis. While the spin carriers of this structure are isotropic $S = ^5/_2$ Fe III metal centers and $S = ^1/_2$ organic linkers, the anisotropy energy of the framework material is comparable to that of reported 2D ferromagnets. Density functional theory calculations indicate that the observed magnetocrystalline anisotropy arises from ligand-to-metal charge transfer that enhances the magnetic anisotropy of the otherwise isotropic Fe centers, suggesting that metal–ligand covalency can be utilized as a general additive for the development of 2D magnets. These results in this study show the possibility for (NMe 4 ) 2 [Fe 2 L 3 ] to retain magnetic order down to the 2D monolayer limit. In addition, the combination of large magnetic anisotropy and semiconducting character in (NMe 4 ) 2 [Fe 2 L 3 ] highlights its potential as a new 2D magnetic semiconductor.

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

Computational Predictions of the Hydrolysis of 2,4,6-Trinitrotoluene (TNT) and 2,4-Dinitroanisole (DNAN)

Hydrolysis is a common transformation reaction that can affect the environmental fate of many organic compounds. In this study, three proposed mechanisms of alkaline hydrolysis of 2,4,6-trinitrotoluene (TNT) and 2,4-dinitroaniline (DNAN) were investigated with plane-wave density functional theory (DFT) combined with ab initio and classical molecular dynamics (AIMD/MM) free energy simulations, Gaussian basis set DFT calculations, and correlated molecular orbital theory calculations. Most of the computations in this study were carried out using the Arrows web based tools. For each mechanism: Meisenheimer complex formation, nucleophilic aromatic substitution, and proton abstraction reaction energies and activation barriers were calculated for the reaction at each relevant site. For TNT, it was found that the most kinetically favorable first hydrolysis steps involve Meisenheimer complex formation by attachment of OH$^{-}$ at the C1 and C3 arene carbons and proton abstraction from the methyl group. The nucleophilic aromatic substitution reactions at the C2 and C4 arene carbons were found to be thermodynamically favorable. However, the calculated activation barriers were slightly lower than previous studies, but still found to be $\Delta$G$^{\ddagger} \approx 18$ kcal/mol using PBEo AIMD/MM free energy simulations, suggesting the reactions are not kinetically significant. For DNAN, the barriers of nucleophilic aromatic substitution were even greater ($\Delta$G$^{\ddagger} > 29$ kcal/mol PBEo AIMD/MM). Further, the most favorable hydrolysis reaction for DNAN was found to be a two-step process in which the hydroxyl first attacks the C1 carbon to form a Meisenheimer complex at the C1 arene carbon C1-(OCH$_3$)OH$^{-}$, subsequently, the methoxy anion (-OCH$_3$) at the C1 arene carbon dissociates and the proton shuttles from the C1-OH to the dissociated methoxy group, resulting in methanol and an aryloxy anion.

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

Energetic and Spectroscopic Properties of Astrophysically Relevant MgC 4 H Radicals Using High-Level Ab Initio Calculations

Considering the importance of magnesium-bearing hydrocarbon molecules (MgC n H; n = 2, 4, and 6) in the carbon-rich circumstellar envelopes (e.g., IRC+10216), a total of 28 constitutional isomers of MgC 4 H have been theoretically investigated using density functional theory (DFT) and coupled-cluster methods. The zero-point vibrational energy corrected relative energies at the ROCCSD(T)/cc-pCVTZ level of theory reveal that the linear isomer, 1-magnesapent-2,4-diyn-1-yl (1, 2 Σ + ), is the global minimum geometry on the MgC 4 H potential energy surface. The latter has been detected both in the laboratory and in the evolved carbon star, IRC+10216. The calculated spectroscopic data for 1 match well with the experimental observations (error ~ 0.78%) which validates our theoretical methodology. Plausible isomerization processes happening among different isomers are examined using DFT and coupled-cluster methods. CASPT2 calculations have been performed for a few isomers exhibiting multireference characteristics. The second most stable isomer, 1-ethynyl-1λ 3 -magnesacycloprop-2-ene-2,3-diyl (2, 2 A 1 , μ = 2.54 D), is 146 kJ mol –1 higher in energy than 1 and possibly the next promising candidate to be detected in the laboratory or in the interstellar medium in future.

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

Designing the Nanoconfined Environment for Energy-Efficient Metal Nanoparticle/Ligand-Based Electrocatalysts

Nanoconfined electrocatalysts demonstrate enhanced selectivity and activity owing to new phenomena emerging from chemical species confined at the interface. Creating these nanoconfined pockets at the interface is not trivial. We need highly negative potentials to create the pockets in metal/ligand catalysts like Ag-nanoparticle/ordered ligand interlayer catalysts. Lowering this activation potential can make the electrocatalysis more energy-efficient, and changes in the environment can impact this potential. We used Density Functional Theory (DFT) to evaluate the impact of ligand length and density, interfacial protons, and surface defects like Ag vacancies and Cu and Au dopants. The nanoconfined configuration is stabilized when the ligands are longer and more densely packed, leading to better agglomeration. It is also stabilized when the ligands adsorb protons and when the ligands detach from the surface easily owing to lower charge transfer. Au-doped surfaces displayed the lowest charge transfer and decreased the activation potential.

Jana, Asmita↗
DOE OSTI2025

Nanoconfinement of High Hydrogen-to-Metal Ratio Lanthanum Hydrides in Functionalized Carbon Hosts

Metal hydrides with a high hydrogen content are important for materials-based hydrogen storage and high-temperature superconductivity. Nanoconfinement of metal hydrides in porous hosts is a promising strategy to tune the thermodynamic stability and control the hydrogen-to-metal ratio. However, lanthanum hydride (LaH x ) nanoconfinement in porous materials has been limited due to the challenges associated with isolating and stabilizing nanoparticles of La or La-hydrides. Here we successfully demonstrated the chemical reduction of La­(III) salts to La(0), and subsequent infiltration into pure CMK-3 and nitrogen-doped CMK-3 (NCMK-3) porous carbons. Transmission electron microscopy measurements revealed a uniform distribution of LaH x species within the carbon hosts, while X-ray absorption and photoelectron spectroscopy provided detailed information about the local chemical environment. Sieverts measurements indicate that LaH x @NCMK-3 could desorb up to 0.75 wt % hydrogen, which is higher than non-nitrogen-functionalized CMK-3 (0.43 wt % H). Density Functional Theory and ab initio molecular dynamics calculations indicate that host–guest interaction energies are favorable for porous carbon with nitrogen defects, which is supported by experimental evidence. Moreover, high-pressure synchrotron X-ray diffraction measurements were conducted using a diamond anvil cell up to 60 GPa and reveal that the nitrogen-functionalized nanoporous carbon host favors the formation of higher H:La ratios in the presence of ammonia borane compared to the pure CMK-3 host. This approach could serve as a suitable platform for developing nanoscale superconducting materials at lower pressures and temperatures compared to bulk.

Shivanna, Mohana↗
DOE OSTI2021

Surface-Enhanced Raman Scattering Optophysiology Nanofibers for the Detection of Heavy Metals in Single Breast Cancer Cells

Mercury(II) ions (Hg2+) and silver ions (Ag+) are two of the most hazardous pollutants causing serious damage to human health. Here, we constructed surface enhanced Raman scattering (SERS)-active nanofibers covered with 4-mercaptopyridine (4-Mpy)-modified gold nanoparticles (AuNPs) to detect Hg2+ and Ag+. Experimental evidence suggests that the observed spectral changes originate from the combined effect of (i) the coordination between the nitrogen on 4-Mpy and the metal ions and (ii) the 4-Mpy molecular orientation (from flatter to more perpendicular with respect to the metal surface). The relative intensity of a pair of characteristic Raman peaks (at ~428 and ~708 cm-1) was used to quantify the metal ion concentration, greatly increasing the reproducibility of the measurement compared to signal-on or signal-off detection based on a single SERS peak. The detection limit of this method for Hg2+ is lower than that for the Ag+ (5 nM vs 100 nM), which can be explained by the stronger interaction energy between Hg2+ and N compared to Ag+ and N, as demonstrated by density functional theory (DFT) calculations. The Hg2+ and Ag+ ions can be masked by adding ethylenediaminetetraacetate (EDTA) and Cl-, respectively, to the Hg2+ and Ag+ samples. The good sensitivity, high reproducibility and excellent selectivity of these nanosensors were also demonstrated. Furthermore, detection of Hg2+ in living breast cancer cells at the subcellular level is possible thanks to the nanometric size of the herein described SERS nanosensors, allowing high spatial resolution and minimal cell damage.

Zhao, Xingjuan↗
DOE OSTI2025

Strained-Induced Morphological Reconstruction of RuO 2 (110) Thin-Film Electrocatalysts

Strain is a widely used strategy for electrocatalyst engineering. Overstraining, however, can lead to unintentional materials transformation. We investigate the impact of strain on the surface morphology of a rutile RuO 2 (110) film grown on a symmetry-matching rutile TiO 2 (110) substrate. When the film thickness exceeds 9 nm, the RuO 2 surface relaxes by forming step edges that expose the {011} plane. Density functional theory (DFT) calculation shows that the (011) facet is among the lower energy surfaces of rutile RuO 2 , suggesting that this formation incurs minimal energetic penalties. In situ atomic force microscopy (AFM) shows that the film maintains the (110) structure of the terrace during electrochemistry. Inductively coupled plasma–mass spectrometry (ICP-MS) further reveals the insensitivity of the Ru dissolution to strain. Here, our findings show a strain-relieving pathway via surface reconstruction in RuO 2 (110) and provide an example of a strain-relieving mechanism that does not affect dissolution.

Evolution reactions↗
DOE OSTI2023

Electrosynthesis of ethylene glycol from C 1 feedstocks in a flow electrolyzer

Ethylene glycol is a widely utilized commodity chemical, the production of which accounts for over 46 million tons of CO 2 emission annually. Here we report a paired electrocatalytic approach for ethylene glycol production from methanol. Carbon catalysts are effective in reducing formaldehyde into ethylene glycol with a 92% Faradaic efficiency, whereas Pt catalysts at the anode enable formaldehyde production through methanol partial oxidation with a 75% Faradaic efficiency. With a membrane-electrode assembly configuration, we show the feasibility of ethylene glycol electrosynthesis from methanol in a single electrolyzer. The electrolyzer operates a full cell voltage of 3.2 V at a current density of 100 mA cm -2 , with a 60% reduction in energy consumption. Further investigations, using operando flow electrolyzer mass spectroscopy, isotopic labeling, and density functional theory (DFT) calculations, indicate that the desorption of a *CH 2 OH intermediate is the crucial step in determining the selectively towards ethylene glycol over methanol.

09 BIOMASS FUELS↗
DOE OSTI2025

A comprehensive approach for elucidating the interplay between 4f n +1 and 4f n 5d 1 configurations in Ln 2+ complexes

Lanthanides (Ln) are typically found in the +3 oxidation state. However, in recent decades, their chemistry has been expanded to include the less stable +2 oxidation state across the entire series except promethium (Pm), facilitated by the coordination of ligands such as trimethylsilylcyclopentadienyl, C 5 H 4 SiMe 3 (Cp'). The [LnCp' 3 ] complexes have been the workhorse for the synthesis and theoretical study of the fundamental aspects of divalent lanthanide chemistry, where experimental and computational evidence have suggested the existence of different ground state (GS) configurations, 4f n+1 or 4f n 5d 1 , depending on the specific metal. Standard reduction potentials and 4f n+1 to 4f n 5d 1 promotion energies have been two factors usually considered to rationalize the occurrence of these variable GS configurations, however the driving force behind this phenomenon is still not clear. In this work we present a comprehensive theoretical approach to shed light on this matter using the [LnCp 3 ] - model systems. We begin by calculating 4f n+1 to 4f n 5d 1 promotion energies and successfully correlate them with existing experimental data. Furthermore, we analyze how changes in the GS charge distribution between the Ln ions, LnCp 3 and the reduced [LnCp 3 ] - complexes (Ln = La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, Lu) correlate with experimental trends in redox potentials and the calculated promotion energies. For this purpose, a comprehensive theoretical work that includes relativistic ligand field density functional theory (LFDFT) and relativistic ab initio wavefunction methods was performed. This study will help the rational design of suitable environments to tune the different GS configurations as well as modulating the spectroscopic properties of new Ln 2+ complexes.

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

Pb dimerization greatly accelerates charge losses in MAPbI 3 : Time-domain ab initio analysis

Metal halide perovskites constitute a new type of semiconducting materials with long charge carrier lifetimes and efficient light-harvesting. The performance of perovskite solar cells and related devices is limited by nonradiative charge and energy losses, facilitated by defects. Combining nonadiabatic molecular dynamics and time-domain density functional theory, we demonstrate that charge losses depend strongly on the defect chemical state. By considering an extra Pb atom in CH 3 NH 3 PbI 3 , which is a common defect in lead halide perovskites, we investigate its influence on charge trapping and recombination. In a chemically inert form as a Pb interstitial, the extra Pb atom has only a mild influence on charge recombination. However, if the extra Pb atom binds to a native Pb atom to form a dimer, the charge trapping and recombination are greatly accelerated because the Pb-dimer creates a localized midgap trap state that couples strongly to the perovskite valence band edge. Holes disappear from the valence band two orders of magnitude faster than in the pristine perovskite and recombine with conduction band electrons one order of magnitude faster. The simulations identify the phonon modes involved in the nonradiative charge trapping and recombination and highlight the importance of rapid decoherence within the electronic subsystem for long carrier lifetimes. The detailed atomistic analysis of the charge trapping and recombination mechanisms enriches the understanding of defect properties and provides theoretical guidance for improving perovskite performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2021

Ionization state and dielectric constant in cold rarefied hydrocarbon plasmas of inertial confinement fusion

A combined approach to study cold rarefied matter is introduced that includes a semianalytical method based on the free-energy minimization and ab initio calculations based on the finite-temperature density-functional theory. The approach is used to calculate the ionization state of hydrocarbon (CH) under the shock-release conditions in inertial confinement fusion. Furthermore, the dielectric constant of CH is calculated using the Kubo-Greenwood formulation and contribution from atomic polarizabilities is found to be as important as the free-electron contribution. Using the ionization state and dielectric constant, the electron density profile in the rarefaction wave of the shock-release plasma is obtained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Sulfur-enhanced dynamics of coinage metal(111) surfaces: Step edges versus terraces as locations for metal-sulfur complex formation

The propensity of trace amounts of sulfur adsorbed on coinage metal(111) surfaces to dramatically enhance surface dynamics has been demonstrated by STM observations of accelerated 2D island decay for Cu and Ag. It is generally accepted that this enhancement is due to the formation of adsorbed metal-sulfur complexes, which facilitate surface mass transport of the metal. These complexes were originally proposed to form on terraces following the extraction of metal atoms from step edges and subsequent combination with sulfur on the terraces. However, even when thermodynamically feasible, this mechanism may not be kinetically viable for some complexes due to limited coupling of the complex concentration to the surface diffusion flux of metal atoms. Focusing on the case of Cu, we assess various scenarios where complexes are formed either on terraces or instead directly at step edges, the latter being a new paradigm. In this work, a new pathway is proposed for the formation on terraces. A rich variety of structures incorporating S at step edges exist, which could provide a viable source for complexes, at least from a thermodynamic perspective. However, it is necessary to also assess the activation barrier for complex formation and detachment from step edges. This is facilitated by the nudged-elastic-band analysis of the minimum energy path for this process utilizing machine-learning derived potentials based on density functional theory energetics for the metal-sulfur system.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Helicity-dependent photocurrents in the chiral Weyl semimetal RhSi

Weyl semimetals are crystals in which electron bands cross at isolated points in momentum space. Associated with each crossing point (or Weyl node) is a topological invariant known as the Berry monopole charge. The circular photogalvanic effect (CPGE), whereby circular polarized light generates a helicity-dependent photocurrent, is a notable example of a macroscopic property that emerges directly from the topology of the Weyl semimetal band structure. Recently, it was predicted that the amplitude of the CPGE associated with optical transitions near a Weyl node is proportional to its monopole charge. In chiral Weyl systems, nodes of opposite charge are nondegenerate, opening a window of wavelengths where the CPGE resulting from uncompensated Berry charge can emerge. Here, we report measurements of CPGE in the chiral Weyl semimetal RhSi, revealing a CPGE response in an energy window that closes at 0.65 eV, in agreement with the predictions of density functional theory.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Effects of Salt Aggregation in Perfluoroether Electrolytes

Electrolytes comprised of polymers mixed with salts have great potential for enabling the use of Li metal anodes in batteries for increased safety. Ionic conductivity is one of the key performance metrics of these polymer electrolytes and achieving high room-temperature conductivity remains a challenge to date. For a bottom-up design of the polymer electrolytes, we must first understand how the structure of polyelectrolytes on a molecular level determines their properties. Here, we use classical molecular dynamics to study the solvation structure and ion diffusion in electrolytes composed of a short-chain perfluoroether with LiFSI or LiTFSI salts. Density functional theory is also used to provide some insights into the structures and energies of the salt interactions with the perfluoroether. We observe the formation of aggregates of salts in the fluorinated systems even at low salt concentrations. The fluorine-fluorine attraction in the solvent is the governing factor for creating the salt aggregates. The aggregates’ size and lifetime change with concentration and anion. These simulations provide an insight into the structure and dynamics of perfluoroether based electrolytes that can be used to improve Li-ion batteries.

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

FETS-VFFA Cell using Contour-Based Magnets

A pure-conductor winding arrangement has been found that approximates the optics of the design given in (Figure 1). The tunes are not the same but are reasonably flat and exist through the entire 3–12 MeV energy range, so could conceivably be fine-tuned using the vertical winding density function in the next section. The lattice consists of two magnets: a reverse bending Bd magnet followed by a normal bending Bf magnet. These magnets are both sectors whose midpoints are displaced with respect to the midpoints of the sides of a 20-sided polygon with side length 1.25 m. The ring is said to have 10 cells (2 polygon sides per cell) that are 2.5 m long. The displacements are 17 cm inwards for Bd and 17 cm outwards for Bf.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Semiconductor Nanostructures by Scientific Design (Final Technical Report)

We developed and used quantum simulation techniques to predict structural, electronic, optical and thermal properties affecting the scientific design of semiconductor nanostructures. We focused on the study of surfaces and interfaces at the nanoscale, in particular embedding and solvation of semiconductor nanoparticles, and provided a microscopic description of nanostructured materials in realistic environments, directly comparable with experimental conditions. We addressed three major problems: (i) the design of realistic structural models of integrated nanostructures, e.g. quantum dots embedded in solid matrices and nanostructures immersed in fluids; (ii) the development of efficient and accurate calculations to describe the excited state properties of nanostructures, beyond standard Density Functional Theory, so as to understand and predict electronic properties involved in renewable energy applications, for example solar conversion devices. In particular we developed and used methods based on a many body perturbation theory. (iii) the development of calculations of heat transport properties of nanostructured materials, for thermoelectric applications, including predictive calculations of thermal and electrical conductivity and of Seebeck coefficients. Specifically, we developed and used methods based on the Boltzmann transport equation and molecular dynamics.

77 NANOSCIENCE AND NANOTECHNOLOGY↗