SEARCH · Engineering Papers
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.
In this study, we present density functional theory calculations within the generalized gradient approximation of U 2 RhC 2 and U 2 PtC 2 . We find the calculated density of states are significantly less than that measured by specific heat indicating the need for electronic correlations. The mass enhancement found for U 2 PtC 2 is m*/m band ≈ 4.
Vacancy-mediated diffusion barriers in 316 stainless steel have been systematically calculated using density functional theory to provide essential parameters for mesoscale microstructure evolution models. A statistical sampling approach employing 210 nudged elastic band calculations across multiple special quasi-random structures captures the effects of local chemical environments in this concentrated alloy. The computational methodology addresses challenges specific to chemically disordered systems, including proper magnetic treatment throughout multi-step calculations and validation against experimental structural properties. The calculated activation barriers reveal clear species-dependent diffusion behavior with the hierarchy Ni >> Fe ˜ Cr >> Mo. Nickel exhibits the highest barriers (0.74–1.31 eV, mean 1.045 eV), confirming its role as the slowest-diffusing major component. Iron and chromium show similar moderate barriers averaging 0.587 eV and 0.522 eV, respectively. Remarkably, molybdenum demonstrates exceptionally low barriers (0.12–0.28 eV, mean 0.194 eV), suggesting much higher mobility than previously recognized and potentially significant implications for precipitation kinetics and microstructure evolution. The barrier ranges remain consistent across different 316 SS compositions, supporting parameter transferability for modeling applications. The overall mean barrier of 0.64 eV provides a practical approximation for phase field simulations, while species-specific values enable detailed treatments of diffusion-controlled processes. This systematic approach establishes a validated framework for generating diffusion parameters in other concentrated alloys where experimental data are limited, while providing the first systematic set of species-specific barriers for predictive modeling of 316 stainless steel microstructure evolution.
The anion photoelectron (PE) spectra of EuH{sup −} and the PE spectrum of overlapping EuOH{sup −} and EuO{sup −} anions are presented and analyzed with supporting results from density functional theory calculations on the various anions and neutrals. Results point to ionically bound, high-spin species. EuH and EuOH anions and neutrals exhibit analogous electronic structures: Transitions from {sup 8}Σ{sup −} anion ground states arising from the 4f{sup 7}σ{sub 6s}{sup 2} superconfiguration to the close-lying neutral {sup 9}Σ{sup −} and {sup 7}Σ{sup −} states arising from the 4f{sup 7}σ{sub 6s} superconfiguration are observed spaced by an energy interval similar to the free Eu{sup +} [4f{sup 7}6s] {sup 9}S - {sup 7}S splitting. The electron affinities (EAs) of EuH and EuOH are determined to be 0.771 ± 0.009 eV and 0.700 ± 0.011 eV, respectively. Analysis of spectroscopic features attributed to EuO{sup −} photodetachment is complicated by the likely presence of two energetically competitive electronic states of EuO{sup −} populating the ion beam. However, based on the calculated relative energies of the close-lying anion states arising from the 4f{sup 7}σ{sub 6s} and 4f{sup 6}σ{sub 6s}{sup 2} configurations and the relative energies of the one-electron accessible 4f{sup 7} and 4f{sup 6}σ{sub 6s} neutral states based on ligand-field theory [M. Dulick, E. Murad, and R. F. Barrow, J. Chem. Phys. 85, 385 (1986)], the remaining features are consistent with the 4f{sup 6}σ{sub 6s}{sup 2} {sup 7}Σ{sup −} and 4f{sup 7}σ{sub 6s}{sup 7}Σ{sup −} anion states lying very close in energy (the former was calculated to be 0.15 eV lower in energy than the latter), though the true anion ground state and neutral EA could not be established unambiguously. Calculations on the various EuO anion and neutral states suggest 4f-orbital overlap with 2p orbitals in species with 4f{sup 6} occupancy.
The formation of intermetallic secondary phases, such as Pt{sub 3}Pb, has been observed experimentally at PbTiO{sub 3}/Pt and Pb(Zr,Ti)O{sub 3}/Pt, or PZT/Pt, interfaces. Density functional theory calculations are used here to calculate the work of adhesion of these interfacial systems with and without the secondary intermetallic phase. The charge density maps of the interfaces reveal the electronic interactions at the interface and the impact of the secondary phase. In addition, Bader charge analysis provides a quantitative assessment of electron transfer from the perovskites to the Pt. Analysis of the band diagrams indicates an increase of the potential barrier associated with electron transfer due to the formation of the Pt{sub 3}Pb at PZT/Pt interfaces.
Oxygen vacancies are a common defect in oxides and play vital roles in many technological applications such as oxygen separation, catalytic reactors, solid oxide fuel cells, and solar thermochemical water splitting. The oxygen vacancy formation energy is directly related to the oxide reduction enthalpy and the vacancy concentration and, therefore, is a key quantity underlying these applications. The complexity of measuring oxygen vacancy formation energies experimentally suggests the utility of calculations based on density functional theory (DFT). However, the calculated results are strongly dependent on the exchange-correlation functionals and other parameters such as the Hubbard U . In this work, we compare the performance of the meta-GGA strongly constrained and appropriately normed (SCAN) functional and the commonly used semilocal generalized gradient approximation (GGA) functionals with experimental values of structural parameters, band gaps, magnetic structures, and oxygen vacancy formation energies for six representative oxides, i.e., CaMnO 3 , SrMnO 3 , LaMnO 3 , YMnO 3 , LaFeO 3 , and CeO 2 . Further, our results show that SCAN usually has better agreement with the experimental lattice constants and band gaps and larger magnetic moments than the commonly used GGA functionals. Although SCAN overestimates the oxygen vacancy formation energies of transition metal oxides and therefore requires unusually large Hubbard U values to reproduce the experimental reduction enthalpies, it does predict the correct oxygen vacancy formation energy of CeO 2 , which challenges the commonly used GGA and hybrid functionals. Our results underscore the challenges that exist in describing these complex oxides by DFT and may shed light on developing more accurate exchange-correlation functionals for oxygen vacancy formation energy calculations.
The mechanical properties of perovskite oxides depend on two metal oxide lattices that are intercalated. This provides an opportunity for separate tuning of hardness, Poisson's ratio (transverse expansion in response to the compression), and shear strength. The elastic constants of series of perovskite oxides were studied by first principles approach. Both A-site and B-site cations were systematically varied in order to see their effects on the elastic parameters. To study the effects of A-site cations, we studied the elastic properties of perovskite ATiO{sub 3} for A being Be, Mg, Ca, Sr, or Ba, one at a time. Similarly, for B-site cations, we studied the elastic properties of PbBO{sub 3} for B being Ti, Zr, or Hf, one at a time. The density functional first principles calculations with local density approximation (LDA) and generalized gradient approximation (GGA) were employed. It is found that the maximum C{sub 11} elastic constant is achieved when the atomic size of the cations at A-site and B-site are comparable. We also found that C{sub 12} elastic constant is sensitive to B-site cations while C{sub 44} elastic constant is more sensitive to A-site cations. Details and explanations for such dependencies are discussed.
The atomic structure and electronic properties of Cu{sub n} nanoclusters (n = 4, 6, 7, and 10) supported on cubic nonpolar δ-MoC(001) and orthorhombic C- or Mo-terminated polar β-Mo{sub 2} C(001) surfaces have been investigated by means of periodic density functional theory based calculations. The electronic properties have been analyzed by means of the density of states, Bader charges, and electron localization function plots. The Cu nanoparticles supported on β-Mo{sub 2} C(001), either Mo- or C-terminated, tend to present a two-dimensional structure whereas a three-dimensional geometry is preferred when supported on δ-MoC(001), indicating that the Mo:C ratio and the surface polarity play a key role determining the structure of supported clusters. Nevertheless, calculations also reveal important differences between the C- and Mo-terminated β-Mo{sub 2} C(001) supports to the point that supported Cu particles exhibit different charge states, which opens a way to control the reactivity of these potential catalysts.
Rate constants for the reactions of O{sub 2}{sup −} and O{sup −} with N and O atoms have been measured for the first time as a function of temperature from 173 to 500 K for O{sup −} reactions and 173 to 400 K for O{sub 2}{sup −} reactions. Room temperature rate constants for O{sub 2}{sup −} reacting with N and O are 3.1 × 10{sup −10} and 1.7 × 10{sup −10} cm{sup 3} s{sup −1}, respectively, and the corresponding O{sup −} rate constants are 1.7 × 10{sup −10} and 1.5 × 10{sup −10} cm{sup 3} s{sup −1}, in good agreement with previous values. Temperature dependences are about T{sup −1.7} for both O{sub 2}{sup −} reactions and T{sup −0.6} and T{sup −1.3} for the reactions of O{sup −} with N and O, respectively. Branching for the O{sub 2}{sup −} reaction with N is found to predominantly form O{sup −} (>85%) in contrast to previous measurements, which reported NO{sub 2}+e{sup −} as the main channel. Calculations point to the present results being correct. The potential energy surface for this reaction was calculated using density functional theory, coupled cluster with singles, doubles (triples), complete active space self-consistent field, and complete active space second-order perturbation methods and is found to be quite complex, with agreement between the calculated surface and the observed kinetic data only possible through the inclusion of dynamical correlation.
The lattice diffusion of silver with various charge states in cubic silicon carbide has been investigated by means of high-fidelity spin-polarized density functional theory calculations. The migration energy barrier of the Ag interstitial diffusion is 1.09 eV and 1.11 eV for neutral and q = +1e charge state, respectively, close to the activation energy of Ag diffusion measured in the German HTR fuel program. A general trend is that the migration energy barrier reduces with respect to the increase of charge on Ag, which is much less than the increment in the oxidation energy, suggesting that the lattice diffusion of silver prefers constant neutral state without redox in transition state. Furthermore, our results indicate a scenario that once Ag is deposited to interstitials via the kickout mechanism, it will highly likely perform lattice diffusion across the SiC layer leading to fast release of Ag in Tristructural-Isotropic fuel induced by irradiation.
We have determined experimentally the magnitude of the difference in the splay and bend flexoelectric coefficients, |e{sub 1} − e{sub 3}|, of an oxadiazole bent-core liquid crystal by measuring the critical voltage for the formation of flexodomains together with their wave number. The coefficient |e{sub 1} − e{sub 3}| is found to be a factor of 2–3 times higher than in most conventional calamitic nematic liquid crystals, varying from 8 pCm{sup −1} to 20 pCm{sup −1} across the ∼60 K—wide nematic regime. We have also calculated the individual flexoelectric coefficients e{sub 1} and e{sub 3}, with the dipolar and quadrupolar contributions of the bent-core liquid crystal by combining density functional theory calculations with a molecular field approach and atomistic modelling. Interestingly, the magnitude of the bend flexoelectric coefficient is found to be rather small, in contrast to common expectations for bent-core molecules. The calculations are in excellent agreement with the experimental values, offering an insight into how molecular parameters contribute to the flexoelectric coefficients and illustrating a huge potential for the prediction of flexoelectric behaviour in bent-core liquid crystals.
Low-temperature standard NH 3 -SCR over copper-exchanged zeolite catalysts occurs on NH 3 -solvated Cu-ion active sites in a quasi-homogeneous manner. As key kinetically relevant reaction steps, the reaction intermediate Cu II (NH 3 ) 4 ion hydrolyzes to Cu II (OH)(NH 3 ) 3 ion to gain redox activity. The Cu II (OH)(NH 3 ) 3 ion also transfers between neighboring zeolite cages to form highly reactive reaction intermediates. Via operando electron paramagnetic resonance spectroscopy and SCR kinetic measurements and density functional theory calculations, we demonstrate here that such kinetically relevant steps become energetically more difficult with lower support Brønsted acid strength and density. Consequently, Cu/LTA displays lower Cu atomic efficiency than Cu/CHA and Cu/AEI, which can also be rationalized by considering differences in their support topology. By carrying out hydrothermal aging to eliminate support Brønsted acid sites, both Cu II (NH 3 ) 4 ion hydrolysis and Cu II (OH)(NH 3 ) 3 ion migration are hindered, leading to a marked decrease in Cu atomic efficiency for all catalysts.
The halide perovskitoid compound Cs 3 Bi 2 I 9 (CBI) has attracted considerable interest as a semiconductor because of its outstanding stability and reduced toxicity compared with lead-based halide perovskites. Here, we report the growth of nuclear radiation detection grade CBI bulk crystals (Φ 15 × 60 mm 3 ) with a high resistivity of over 10 10 Ω cm using a modified vertical Bridgman method. Because of their layered crystal structures we investigated the anisotropy in the optical and electrical properties using different crystal orientations. The CBI (001) sample exhibits a resistivity of ~10 12 Ω cm compared to ~10 10 Ω cm for the CBI (100) sample. This is due to the anisotropic mobility in the two crystallographic directions. Using 425 nm LED (~200 mW cm –2 ) illumination CBI (001) possesses a superior optical response with a switching ratio of over 40, which is critically higher than that of CBI (100) (<2). Detectors of 2 mm thickness show a capability of detecting 241 Am@5.49 MeV α particles, with good peak discrimination. A full width at half maximum (FWHM) of 32% was obtained under a bias of 560 V. Simultaneously, the electron mobility and mobility lifetime ( μτ ) were calculated to be 6.10 cm 2 V –1 s –1 and 2.03 × 10 –5 cm V –1 , respectively. First-principles density functional theory calculations confirm the crystallographic anisotropy of the carrier effective masses. Additionally, a significant X-ray sensitivity of 111.9 μC Gy –1 cm –2 for a CBI detector was observed, under 80 kVp X-rays at an electrical field of 450 V cm –1 .
Refractory high entropy alloys (RHEAs) have been proven to be a potential candidate in the biomedical field due to their balanced mechanical properties and biocompatible composition. Recent experimental findings show that RHEAs like HfNbTaTiZr and Hf 0.5 Nb 0.5 Ta 0.5 Ti 1.5 Zr have good mechanical properties such as high polarization and wear resistance than others which establish them as potential materials for biomedical application. In this work, we performed first-principles density functional theory calculations on the mechanical and thermal properties of HfNbTaTiZr and Hf 0.5 Nb 0.5 Ta 0.5 Ti 1.5 Zr. The predicted lattice constant, density, Young's modulus, and Vickers hardness are consistent with the available experimental report, which verifies the accuracy of the applied model. The thermal coefficient of linear expansion of both RHEAs has been investigated by utilizing the Debye theory. The present methods could be applied to study other future RHEAs on exploration of their physical properties.
Structural phase transition sequence of CdSe has been investigated at pressures up to 60 GPa under quasi-hydrostatic conditions using synchrotron X-ray diffraction. A phase transition from the wurtzite type (B4) to the NaCl-type (B1) structure has been observed, followed by another phase transition to an orthorhombic structure at 27 GPa, in agreement with previous reports. We show that this high-pressure orthorhombic phase has a Pnma symmetry rather than being a Cmcm-symmetric structure as previously suggested. From our observations, the appearance of the new reflections and reflection splitting with increasing pressure is due to the change of atomic relative positions in crystal lattice and the difference in the compression ratio of lattice parameters for the Pnma structure, and we find no evidence for the third phase transition reported previously. The pressure-induced phase transition of CdSe has been further confirmed by the density-functional theory calculations.
We investigate entirely electronic torsional vibrational modes in linear cumulene chains. The carbon nuclei of a cumulene are positioned along the primary axis so that they can participate only in the transverse and longitudinal motions. However, the interatomic electronic clouds behave as a torsion spring with remarkable torsional stiffness. The collective dynamics of these clouds can be described in terms of electronic vibrational quanta, which we name torsitons. It is shown that the group velocity of the wavepacket of torsitons is much higher than the typical speed of sound, because of the small mass of participating electrons compared to the atomic mass. For the same reason, the maximum energy of the torsitons in cumulenes is as high as a few electronvolts, while the minimum possible energy is evaluated as a few hundred wavenumbers and this minimum is associated with asymmetry of zero point atomic vibrations. Theory predictions are consistent with the time-dependent density functional theory calculations. Molecular systems for experimental evaluation of the predictions are proposed.
We present that dipole moments are a simple, global measure of the accuracy of the electron density of a polar molecule. Dipole moments also affect the interactions of a molecule with other molecules as well as electric fields. To directly assess the accuracy of modern density functionals for calculating dipole moments, we have developed a database of 200 benchmark dipole moments, using coupled cluster theory through triple excitations, extrapolated to the complete basis set limit. This new database is used to assess the performance of 88 popular or recently developed density functionals. The results suggest that double hybrid functionals perform the best, yielding dipole moments within about 3.6-4.5% regularized RMS error versus the reference values - which is not very different from the 4% regularized RMS error produced by coupled cluster singles and doubles. Many hybrid functionals also perform quite well, generating regularized RMS errors in the 5-6% range. Lastly, some functionals, however, exhibit large outliers, and local functionals in general perform less well than hybrids or double hybrids.
Shear strength measurements have been carried out on rhenium diboride, ReB 2 , to a pressure of 74 GPa using a Radial X-ray Diffraction (R-XRD) technique in a diamond anvil cell using platinum as an internal x-ray pressure standard. The R-XRD result has provided a unique insight into the deformation of hexagonal crystal lattice under non-hydrostatic compression and variation of shear strength with increasing pressure. From R-XRD data, we have estimated hydrostatic component of compression to determine an equation of state of rhenium diboride yielding a bulk modulus of K 0 = 366 ± 25 GPa with a pressure derivative $K^{'}_{0}$ = 4.3 ± 0.5 in good agreement with hydrostatic density functional theory calculations. Here, the average lower bound of shear strength (τ) from various diffraction planes was then calculated using the measured interplanar d-spacing (d m ) and hydrostatic component of d-spacing (d p ) to be shown to approach 6.7 ± 0.4 GPa at 70 GPa. Our results show that the anisotropic compression effects observed in ReB 2 under hydrostatic compression are correlated to electronic structure changes under compression as predicted by theoretical calculations.
Steady state rate expressions can be derived to distinguish the Cossee–Arlman and metallacycle mechanisms postulated for propylene oligomerization on nickel-based catalysts based on product selectivities, where product selectivities for the former are a function of olefin pressure because sequential coordination and insertion steps lead to independent mechanistic pathways for different hexene isomers. In contrast, the metallacycle mechanism presents pressure-independent product selectivities due to successive coordination prior to the kinetically relevant steps in each mechanism. In this work, steady state propylene oligomerization rates and selectivities were measured in the absence of an activator on nickel functionalized UiO-66 metal organic framework (MOF), Ni/UiO-66, to validate the Cossee–Arlman mechanism for light olefin oligomerization. In situ NO titrations reveal that ~5% of nickel sites were active during the reaction, and thus, not all nickel sites are relevant for catalysis. Propylene dimerization was first order in propylene pressure from 5 to 500 kPa with an apparent activation energy of ~20 kJ mol -1 from 453 to 493 K. Calculated apparent activation energies with density functional theory (DFT) calculations on cluster models of Ni/UiO-66 are in agreement with experiment to corroborate the Cossee–Arlman mechanism. Selectivities of hexene products and the ratio of hexene product selectivities on Ni/UiO-66 are in accordance with selectivity expressions derived from the Cossee–Arlman mechanism. In conclusion, analysis of product selectivities can be used more extensively to demarcate the Cossee–Arlman and metallacycle mechanisms for olefin oligomerization on metal-based catalysts.