Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Mechanical and Thermal Properties of Low-Density Al 20+x Cr 20-x Mo 20-y Ti 20 V 20+y Alloys

Refractory high-entropy alloys (RHEAs) Al 20+x Cr 20-x Mo 20-y Ti 20 V 20+y ((x, y) = (0, 0), (0, 10), and (10, 15)) were computationally studied to obtain a low density and a better mechanical property. The density functional theory (DFT) method was employed to compute the structural and mechanical properties of the alloys, based on a large unit cell model of randomly distributed elements. Debye–Grüneisen theory was used to study the thermal properties of Al 20+x Cr 20-x Mo 20-y Ti 20 V 20+y . The phase diagram calculation shows that all three RHEAs have a single body-centered cubic (BCC) structure at high temperatures ranging from 1000 K to 2000 K. The RHEA Al 30 Cr 10 Mo 5 Ti 20 V 35 has shown a low density of 5.16 g/cm 3 and a hardness of 5.56 GPa. The studied RHEAs could be potential candidates for high-temperature application materials where high hardness, ductility, and low density are required.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Metal-Mediated Directional Capping of Rod-Packing Metal–Organic Frameworks

We report two new rod-packing metal–organic frameworks (RPMOF) are constructed by regulating the in situ formation of the capping agent. In CPM-s7, carboxylate linkers extend 1D manganese-oxide chains in four additional directions, forming 3D RPMOF. The substitution of Mn 2+ with a stronger Lewis acidic Co 2+ , leads to an acceleration of the hydrolysis-prone sulfonate linker, resulting in presence of sulfate ions to reduce two out of the four carboxylate-extending directions, and thus forming a new 2D rod-packing CPM-s8. Density functional theory calculations and magnetization measurements reveal ferrimagnetic ordering of CPM-s8, signifying the potential of exploring 2D RPMOF for effective low-dimensional magnetic materials.

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

Structural anisotropy in Sb thin films

Sb thin films have attracted wide interest due to their tunable band structure, topological phases, high electron mobility, and thermoelectric properties. We successfully grow epitaxial Sb thin films on a closely lattice-matched GaSb(001) surface by molecular beam epitaxy. We find a novel anisotropic directional dependence on their structural, morphological, and electronic properties. The origin of the anisotropic features is elucidated using first-principles density functional theory (DFT) calculations. The growth regime of crystalline and amorphous Sb thin films was determined by mapping the surface reconstruction phase diagram of the GaSb(001) surface under Sb 2 flux, with confirmation of structural characterizations. Crystalline Sb thin films show a rhombohedral crystal structure along the rhombohedral (211) surface orientation parallel to the cubic (001) surface orientation of the GaSb substrate. At this coherent interface, Sb atoms are aligned with the GaSb lattice along the [1̄10] crystallographic direction but are not aligned well along the [110] crystallographic direction, which results in anisotropic features in reflection of high-energy electron diffraction patterns, misfit dislocation formation, surface morphology, and transport properties. Our DFT calculations show that the preferential orientation of the rhombohedral Sb (211) plane may originate from the GaSb surface, where Sb atoms align with the Ga and Sb atoms on the reconstructed surface. The formation energy calculations confirm the stability of the experimentally observed structures. Our results provide optimal film growth conditions for further studies of novel properties of Bi 1-x Sb x thin films with similar lattice parameters and an identical crystal structure, as well as functional heterostructures of them with III–V semiconductor layers along the (001) surface orientation, supported by a theoretical understanding of the anisotropic film orientation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Predicting phase transitions in PbTi⁢O 3 using zentropy through quasiharmonic phonon calculations

According to x-ray diffraction (XRD) measurements, PbTi⁢O 3 undergoes a phase transition from a tetragonal ferroelectric (FE) phase to a cubic paraelectric phase at 763 K. However, x-ray absorption fine-structure (XAFS) measurements indicate that PbTi⁢O 3 is locally tetragonal even after the phase transition. The difference in these results is because XAFS measurements can probe local features of a structure, while XRD averages over such local features. For both measurements to be consistent, PbTi⁢O 3 is macroscopically cubic but locally tetragonal after the phase transition. Despite this, most models, such as the Landau-Ginsburg-Devonshire theory and effective Hamiltonians, are still unable to explain this phenomenon. Moreover, these methods involve model parameters fitted to experimental or theoretical data and do not consider other tetragonal configurations, such as domain walls, to predict the phase transition. In our previous study, we used our zentropy approach to predict the phase transition by considering the tetragonal FE ground-state configuration and the tetragonal 90° and 180° domain wall configurations with their total energies at 0 K. Here, in this paper, the Helmholtz energies of the three configurations are obtained from density functional theory calculations through energy-volume curves and phonon calculations. The predicted phase transition temperature using the meta-GGA 𝑟 2⁢ SCAN and revised multiplicities of configurations is 716 K, showing good agreement with the experimental value of 763 K.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Hydrogen and deuterium tunneling in niobium

We use density functional methods to identify the atomic configurations of H and D atoms trapped by O impurities embedded in bulk Nb. The O atoms are located at the octahedral position in the Nb body-centered cubic (BCC) lattice, and H (D) atoms tunnel between two degenerate tetrahedral sites separated by a mirror plane. Using nudged elastic band (NEB) methods, we calculate the double-well potential for O-H and O-D impurities and the wave functions and tunnel splittings for H and D atoms. Our results agree with those obtained from analysis of heat capacity and neutron scattering measurements on Nb with low concentrations of O-H and O-D.

density functional theory↗
DOE OSTI2021

ScO x rich surface terminations on lanthanide scandate nanoparticles

We report that the lanthanide scandate materials are widely used substrates for thin film growth and the potential applications of the LnScO 3 materials continue to grow with the recent ability to synthesize them as cuboidal faceted nanoparticles. A comprehensive understanding of the surface structure and chemistry of these oxides is essential for their informed application, either as single crystal or high surface area nanoparticle substrates. Here the {100} pseudocubic surfaces of LnScO 3 (Ln = La, Nd, Sm, Gd) nanoparticles were examined with aberration-corrected electron microscopy and higher-level density functional theory to reveal ScO x rich terminations across all investigated surfaces. Mixed terminations of single and double layer nature were observed, indicating the presence of multiple domains at the surface and introducing the possibility of synthetically controlling the surface reconstruction in the future.

36 MATERIALS SCIENCE↗
DOE OSTI2026

MnRhBi3: A Cleavable Antiferromagnetic Metal

This dataset contains DFT input and output files supporting the theoretical modeling in the associated publication (Chem. Mater. 2024, 36, 11306-11316). The calculations characterize MnRhBi3, an orthorhombic (Cmmm) van der Waals-layered intermetallic compound that cleaves easily between neighboring Bi layers. The dataset is organized into three calculation types: (i) Bulk: Structural relaxations of the periodic MnRhBi3 crystal in antiferromagnetic (AFM) and ferromagnetic (FM) configurations, using the vdW-DF-optB86b functional. These provide the equilibrium lattice constants, magnetic energy differences (AFM is 0.5 meV/f.u. lower than FM), and magnetic moments (4.4 µB/Mn, 0.17 µB/Rh, 0.18 µB/Bi) reported in Table 1 of the main text. (ii) Slab: Same magnetic configurations computed with an 18 Ang vacuum layer introduced between Bi layers, used to calculate the cleavage energy Ec = 0.56 J/m2 (AFM) and 0.57 J/m2 (FM), establishing MnRhBi3 as a van der Waals-layered material comparable to graphite, MoS2, and CrI3. (iii) ELF: Single-point calculation on the relaxed bulk AFM geometry with LELF=.TRUE., producing the ELFCAR file used to generate electron localization function isosurfaces and contour maps (Fig. 2, main text) showing Bi lone pairs directed into the van der Waals gaps. All folders contain CONTCAR, INCAR, KPOINTS, OUTCAR, and POSCAR. The ELF/ folder additionally contains ELFCAR. Calculations were performed using VASP 6.3.2 with PBE + vdW-DF-optB86b, PAW potentials, and an energy cutoff of 800 eV.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Phase Diagrams of Alloys and Their Hydrides via On-Lattice Graph Neural Networks and Limited Training Data

Efficient prediction of sampling-intensive thermodynamic properties is needed to evaluate material performance and permit high-throughput materials modeling for a diverse array of technology applications. To alleviate the prohibitive computational expense of high-throughput configurational sampling with density functional theory (DFT), surrogate modeling strategies like cluster expansion are many orders of magnitude more efficient but can be difficult to construct in systems with high compositional complexity. We therefore employ minimal-complexity graph neural network models that accurately predict and can even extrapolate to out-of-train distribution formation energies of DFT-relaxed structures from an ideal (unrelaxed) crystallographic representation. This enables the large-scale sampling necessary for various thermodynamic property predictions that may otherwise be intractable and can be achieved with small training data sets. Two exemplars, optimizing the thermodynamic stability of low-density high-entropy alloys and modulating the plateau pressure of hydrogen in metal alloys, demonstrate the power of this approach, which can be extended to a variety of materials discovery and modeling problems.

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

Argon equation of state data to 1 TPa: Shock compression experiments and simulations

Argon is the most abundant noble gas on Earth and its noble, atomic fluid nature makes it an excellent candidate for comparison of experiment and theory at extreme conditions. For this study, we performed a combined computational and experimental study on shock compressed cryogenic liquid argon. Using Sandia's Z machine, we shock compressed liquid argon to 600 GPa and reshock states up to 950 GPa. Laser shock experiments at the Omega Laser facility extend the principal Hugoniot to 1000 GPa and provided temperature data along the principal Hugoniot. The plate impact experiments and laser shock experiments used well-characterized impedance matching standards and demonstrate consistent results between the two platforms over a common range. Density functional theory based molecular dynamics simulations provided additional data on the Hugoniot to 600 GPa. The combined experimental data and simulation results provide constraints on the development of new equation of state models at extreme conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Thermodynamics and Magnetism of YCo5 Compound Doped with Fe and Ni: An Ab Initio Study

YCo5 permanent magnet exhibits high uniaxial magnetocrystalline anisotropy energy and has a high Curie temperature. These are good properties for a permanent magnet, but YCo5 has a low energy product, which is notably insufficient for a permanent magnet. In order to improve the energy product in YCo5, we suggest replacing cobalt with iron, which has a much bigger magnetic moment. With a combination of density-functional-theory calculations and thermodynamic CALculation of PHAse Diagrams (CALPHAD) modeling, we show that a new magnet, YFe3(Ni1-xCox)2, is thermodynamically stable and exhibits an improved energy product without significant detrimental effects on the magnetocrystalline anisotropy energy or the Curie temperature.

anisotropy↗
DOE OSTI2026

4-Electron Oxygen Reduction Reaction (ORR) with Iron Phthalocyanine (FePc) Functionalized Nanowire Templated-3D Fuzzy Graphene (NT-3DFG)

Iron phthalocyanine (FePc) is a promising alternative to platinum-based catalysts for sustainable energy devices; however, the plane-symmetry of Fe-N 4 sites, random aggregation, and poor conductivity of FePc present major barriers for their application as oxygen reduction reaction (ORR) electrocatalysts. Here, we report the synergistic effects of FePc electrocatalysts supported by a nanowire-templated three-dimensional fuzzy graphene (FePc@NT-3DFG) substrate. The in situ functionalized oxygen groups (iFOGs) at the edge of NT-3DFG localize Fe active sites in FePc under alkaline ORR conditions. With a uniform FePc distribution through many single layers of graphene, the NT-3DFG substrates improve O 2 adsorption and catalytic activity while stabilizing the electrochemical activity during reactions. The FePc@NT-3DFG catalyst exhibits fast ORR kinetics with an extremely low Tafel slope of 28.3 ± 2.7 mV dec −1 , a higher half-wave potential of 0.911 ± 0.004 V (vs RHE), and notable long-term stability at 0.5 V (vs RHE) of 96.0 ± 0.4% retention after 30 h. Surface chemistry spectra validate electronic configuration modification of Fe at the iFOGs. Density functional theory calculations indicate that the extra layers of graphene improve oxygen adsorption. Moreover, additional exploration of other transition metal phthalocyanines supports the effects of iFOGs through the transition toward 4e − ORR. This work offers an expanded strategy for active site modification through edge-based graphene substrates for 4e − ORR.

X-ray absorption spectroscopy↗
DOE OSTI2024

Toward first principles-based simulations of dense hydrogen

Accurate knowledge of the properties of hydrogen at high compression is crucial for astrophysics (e.g., planetary and stellar interiors, brown dwarfs, atmosphere of compact stars) and laboratory experiments, including inertial confinement fusion. There exists experimental data for the equation of state, conductivity, and Thomson scattering spectra. However, the analysis of the measurements at extreme pressures and temperatures typically involves additional model assumptions, which makes it difficult to assess the accuracy of the experimental data rigorously. On the other hand, theory and modeling have produced extensive collections of data. They originate from a very large variety of models and simulations including path integral Monte Carlo (PIMC) simulations, density functional theory (DFT), chemical models, machine-learned models, and combinations thereof. At the same time, each of these methods has fundamental limitations (fermion sign problem in PIMC, approximate exchange–correlation functionals of DFT, inconsistent interaction energy contributions in chemical models, etc.), so for some parameter ranges accurate predictions are difficult. Recently, a number of breakthroughs in first principles PIMC as well as in DFT simulations were achieved which are discussed in this review. Here we use these results to benchmark different simulation methods. We present an update of the hydrogen phase diagram at high pressures, the expected phase transitions, and thermodynamic properties including the equation of state and momentum distribution. Furthermore, we discuss available dynamic results for warm dense hydrogen, including the conductivity, dynamic structure factor, plasmon dispersion, imaginary-time structure, and density response functions. We conclude by outlining strategies to combine different simulations to achieve accurate theoretical predictions that are based on first principles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

An Impedance Humidity Sensor Based on CVD Grown WSe 2 2D Films

2D materials-based planar devices have been sought after for gas and humidity sensor applications owing to their high sensitivity due to their ultrathin nature and increased surface area. Herein, 2D WSe 2 films-based planner devices were fabricated to evaluate their performance for potential application in humidity sensors. The structure and morphology of the WSe 2 sensing films were studied by Raman spectroscopy, photoluminescence, optical microscopy, and scanning electron microscopy (SEM) techniques. The relative humidity-dependent sensing performances of the 2D WSe 2 films were evaluated using current–voltage (I–V), relative impedance, and complex total impedance (electrochemical impedance spectroscopy, EIS) techniques. This sensor showed a very stable and repeatable performance investigated over the period of 18 months. The response and recovery times of the WSe 2 -based impedance sensors were found to be 40 and 30 s, respectively. The WSe 2 sensor exhibits a 2 ± 1% minimum hysteresis at a lower humidity level (50% RH) and around 8 ± 2% maximum hysteresis at relative humidity level (70% RH). This value of hysteresis is in the acceptable range during the initial investigation of any material to explore its potential for application in humidity sensors. The density functional theory (DFT) analysis was performed with the interaction from water molecules of the WSe 2 sample for a better understanding of the sensing mechanism of devices.

2D- WSe2↗
DOE OSTI2022

Assessing the Role of Interfacial and Metal Sites in Pt/TiO2-Catalyzed Acetic Acid Hydrodeoxygenation

Thermochemical conversion of biomass to produce drop-in quality biofuels typically involves hydrodeoxygenation (HDO) steps following catalytic fast pyrolysis (CFP) to remove excess oxygen and create a more stable bio-oil product. HDO involves co-feeding the CFP vapor-phase product and H2 gas over a bi-functional catalyst. Reducible metal oxide-supported noble-metal catalysts (e.g., Pt/TiO2) are promising materials for HDO chemistry, with recent work aiming to elucidate the role of various Pt/TiO2 actives sites (i.e., Pt-metal, TiO2-support, Pt-TiO2-interfacial sites) in the competing desired deoxygenation and undesired decarboxylation/decarbonylation reactions for important classes of CFP vapor model compounds. Carboxylic acids are one important class of bio-derived compounds that has not been studied extensively for HDO on Pt/TiO2, particularly at the atomic level. Past experimental work evaluating Pt/C and Pt/TiO2 catalysts for acetic acid HDO (AA-HDO) demonstrated selectivity toward C-C and C-O bond-dissociation products, respectively. This work utilizes atomic-scale modeling to discern the role of Pt-metal and Pt-TiO2-interfacial sites in promoting key C-C bond-breaking, C-O bond-breaking, and (de)hydrogenation steps in AA-HDO. Using Pt(111) and Pt-TiO2-interface surface models, adsorption and reaction energetics calculated by density functional theory provide fundamental insights into the role of interfacial sites and oxygen vacancies in promoting desired deoxygenation pathways over undesired decarboxylation/decarbonylation pathways.

BIOMASS FUELS,INORGANIC, ORGANIC, PHYSICAL, AND AN↗
DOE OSTI2020

Rapid generation of optimal generalized Monkhorst-Pack grids

In this report computational modeling of the properties of crystalline materials has become an increasingly important aspect of materials research, consuming hundreds of millions of CPU-hours at scientific computing centers around the world each year, if not more. A routine operation in such calculations is the evaluation of integrals over the Brillouin zone. We have previously demonstrated that performing such integrals using generalized Monkhorst-Pack k-point grids can roughly double the speed of these calculations relative to the widely-used traditional Monkhorst-Pack grids. However the generation of optimal generalized Monkhorst-Pack grids is not implemented in most software packages due to the computational cost and difficulty of identifying the best grids. To address this problem, we present new algorithms that allow rapid generation of optimal generalized Monkhorst-Pack grids on the fly. We demonstrate that the grids generated by these algorithms are on average significantly more efficient than those generated using existing algorithms across a range of grid densities. For grids that correspond to a real-space supercell with at least 50 Å between lattice points, which is sufficient to converge density functional theory calculations within 1 meV/atom for nearly all materials, our algorithm finds optimized grids in an average of 0.19 s on a single processing core. To facilitate the widespread adoption of this approach, we present new open-source tools including a library designed for integration with third-party software packages.

36 MATERIALS SCIENCE↗
DOE OSTI2021

2D Ni 0.25 Mn 0.75 O 2 : A high-performance cathode for multivalent ion batteries

Although Li-ion batteries have driven portable energy storage in recent decades, there is increasing concern about their safety, cost, and abundance of constituents. Multivalent ion batteries (MVIBs) have the potential to remedy these issues, but they are limited by the currently known MVIB cathodes, which fail to deliver unanimously favorable voltage, energy density, and diffusion kinetics. Here we used density functional theory (DFT) to model the performance of Li, Na, Mg, Ca, and Al ions when paired with 2D Ni 0.25 Mn 0.75 O 2 , a novel cathode that uses increased layer separation to improve on the kinetics of its 3D analog. Our calculations yielded maximum voltages of 3.38 V for Na and 2.7 V for Ca, outperforming 2D Na x MnO 2 and Na x NiO 2 . Diffusion barriers for Li, Na, and Ca are below 300 meV, comparable to existing battery technology and the endpoint 2D cathodes; meanwhile, Mg and Al have prohibitively high diffusion barriers, implying their incompatibility with this cathode. Lastly, density of states calculations and Bader charge analysis show that the cathode becomes conducting following ion adsorption, which is necessary for high-rate performance. 2D Ni 0.25 Mn 0.75 O 2 maintains performance seen with other 2D transition metal oxides while increasing cathode conductivity, indicating that it is a promising candidate for experimental investigation with Li, Na, and Ca ions.

25 ENERGY STORAGE↗
DOE OSTI2020

Efficient computational search for lanthanide-binding additive dopants for advanced U-Zr based fuels

A rational design strategy combining theory and experiment is highly desirable for the development of U-Zr based fuels that are resistant to lanthanide-induced fuel-cladding chemical interaction (FCCI). Here, we performed rapid computational screening of FCCI-mitigating (via strong lanthanide-binding) elements across the periodic table by utilizing three simple criteria and density functional theory (DFT) calculated total energies of relevant elements and compounds retrieved from high-throughput DFT databases. In addition to successfully identifying previously demonstrated additives, our search leads to the discovery of Bi as a highly effective element among all 68 elements screened, and its effectiveness was experimentally confirmed in this study.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗