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.
The electronic and optical properties of orthorhombic (CH 3 NH 3 )BX 3 (B = Sn, Pb; X = F, Cl, Br, I) were investigated by first-principles density functional theory.
With the goal of developing a Si-based anode for Mg-ion batteries (MIBs) that is both efficient and compatible with the current semiconductor industry, the current research utilized classical Molecular Dynamics (MD) simulation in investigating the intercalation of a Mg 2+ ion under an external electric field (E-field) in a 2D bilayer silicene anode (BSA). First principles density functional theory calculations were used to validate the implemented EDIP potentials. Our simulation shows that there exists an optimum E-field value in the range of 0.2–0.4 V Å –1 for Mg 2+ intercalation in BSA. To study the effect of the E-field on Mg 2+ ions, an exhaustive spread of investigations was carried out under different boundary conditions, including calculations of mean square displacement (MSD), interaction energy, radial distribution function (RDF), and trajectory of ions. Our results show that the Mg 2+ ions form a stable bond with Si in BSA. The effects of E-field direction and operating temperature were also investigated. In the X–Y plane in the 0°–45° range, 15° from the X-direction was found to be the optimum direction for intercalation. The results of this work also suggest that BSA does not undergo drastic structural changes during the charging cycles with the highest operating temperature being ~300 K
Here, in this work, we resolve how proton dynamics and halide mixing enhance or impede ionic conduction in protonated lithium antiperovskites (pLiAP) at compositions near the eutectic points of the halide salts. As a material class, pLiAPs of the form Li 3-x OH x X, (X = Cl, Br) show vast compositional design freedom; however, the resulting properties are susceptible to synthesis and processing methodologies. Proton incorporation and halide mixing stabilize the perovskite cubic phase at low temperatures (<50 °C) and using halide mixtures near the eutectic points (~250 to 300 °C) offer possibilities of lower temperature and faster synthesis and processing conditions (<1 h). Mixed-halide compositions such as Li 2 OHCl 0.37 Br 0.63 lead to a 30-fold improvement in room temperature ionic conductivity of a single halide structure, 1.5 × 10 -6 vs. 4.9 × 10 -8 S cm -1 (Li 2 OHCl). We combine infrared spectroscopy and nuclear magnetic resonance with first-principles density functional theory calculations to deconvolute halide mixing effects from local proton dynamics on Li-ion transport. In contrast to what has been supposed, our findings suggest that the halide sublattice dynamics, besides the OH rotation, correlate strongly with the fast-ion conduction at high temperatures.
The first two-dimensional (2D) polymorphs of antimony dioxide, namely, γ-Sb 2 O 4 and δ-Sb 2 O 4 , are predicted using the evolutionary algorithm combined with first-principles density functional theory (DFT) calculations. Out-of-plane ferroelectricity is found in γ-Sb 2 O 4 , while in-plane ferroelectricity is found in δ-Sb 2 O 4 . The predicted dipole moments of γ-Sb 2 O 4 and δ-Sb 2 O 4 phases are 36.63 and 14.96 eÅ, respectively, implying that they can be good candidates for making ferroelectric devices. Further, the calculations show that doping with other group V elements or applying strain can lower the switching energy barriers and thus facilitate switching. Results from GW calculations show indirect band gaps of 5.51 and 3.39 eV for γ-Sb 2 O 4 and δ-Sb 2 O 4 in their monolayers, respectively. Raman spectra are calculated to facilitate the experimental investigation of the predicted structures. The existence of both in-plane and out-of-plane 2D ferroelectricity and the large band gaps make this material system particularly interesting for potential applications.
Double transition metal (DTM) nitride MXenes offer enhanced electrical conductivity and tunable catalytic properties compared to conventional carbide-based MXenes. In this study, we employed first-principles density functional theory (DFT) calculations to discover and computationally validate a novel DTM nitride MXene, Nb 2 TiN 2 , derived from its MAX phase precursor and investigated its potential as an anchoring material (AM) for Li–Se batteries. This newly proposed MXene expands the compositional landscape of DTM nitrides and opens new avenues for functional material design. We performed a comprehensive analysis of the thermodynamic and electronic properties of Nb 2 TiAlN 2 , and the MAX phase precursor to Nb 2 TiN 2 to assess its structural stability and exfoliation potential. Exfoliation energy calculations confirmed the feasibility of synthesizing Nb 2 TiN 2 from Nb 2 TiAlN 2 . We then explored the functionalized form, Nb 2 TiN 2 S 2 , evaluating its capability to serve as an effective anchoring material (AM) in Li–Se batteries by analyzing the reaction mechanisms and kinetics of the selenium reduction reaction (SeRR). Our results indicate that Nb 2 TiN 2 S 2 exhibits a strong binding affinity for lithium polyselenides (Li 2 Se n ), effectively suppressing the shuttle effect. Gibbs free energy calculations for the rate-limiting step of the SeRR reveal favorable kinetics and reduced reaction barriers. Overall, this study provides a detailed evaluation of the structural and electronic properties of a newly proposed DTM nitride MXene and its S-functionalized derivative and the catalyzing effect of Nb 2 TiN 2 S 2 in accelerating the reaction kinetics in Li–Se batteries. These findings underscore the potential importance of the further exploration of MXenes to address current challenges in high-performance Li–Se batteries.
Interest exists in utilizing boron (B) wall conditioning of fusion tokamaks containing tungsten (W) plasma facing components, in order to improve plasma confinement. To understand the interactions of B with W surfaces, first-principles density functional theory calculations have been performed to model the adsorption, diffusion, and solution of B near the W(100), W(110), and W(111) surfaces. The results show that B within a distance of 0.6 nm above the surfaces is adsorbed to the surfaces without activation barriers. B atoms are strongly adsorbed on the W(100) surface with an adsorption energy of 7.80 eV, which is 1.22 and 1.35 eV larger than on the W(110) and W(111) surfaces. B diffusion on the W(100), W(110), and W(111) surfaces has an activation energies of 2.08, 1.12, and 1.47 eV, respectively; while, diffusion from the adsorption sites into the bulk requires 2.2–2.3 eV. The B solution energy below a clean W(100) surface is the lowest, followed by the W(111) and W(110) surfaces. B clustering and B-induced surface deformation as a function of B coverage has been investigated. B on the W(100) surface occupy epitaxial sites at coverages of 0–1.25 ML, but form clusters at higher coverages. B clustering on the W(110) and W(111) surfaces is expected throughout the adsorption process. Compared to a clean surface, B atoms on the W(100) reduce the surface effect on the B solution energy below the surface, while the presence of B on the W(110) and W(111) surfaces generally decreases or increases the B solution energy below surfaces, respectively.
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.
Recent band structure calculations have suggested the potential for band tuning in the chiral semiconductor Ag 3 AuTe 2 to zero upon application of negative strain. In this study, we report on the synthesis of polycrystalline Ag 3 AuTe 2 and investigate its transport and optical properties and mechanical compressibility. Transport measurements reveal the semiconducting behavior of Ag 3 AuTe 2 with high resistivity and an activation energy E a of 0.2 eV. The optical bandgap determined by diffuse reflectance measurements is about three times wider than the experimental E a . Despite the difference, both experimental gaps fall within the range of predicted bandgaps by our first-principles density functional theory (DFT) calculations employing the Perdew–Burke–Ernzerhof and modified Becke–Johnson methods. Furthermore, our DFT simulations predict a progressive narrowing of the bandgap under compressive strain, with a full closure expected at a strain of –4% relative to the lattice parameter. To evaluate the feasibility of gap tunability at such substantial strain, the high-pressure behavior of Ag 3 AuTe 2 was investigated by in situ high-pressure x-ray diffraction up to 47 GPa. Mechanical compression beyond 4% resulted in a pressure-induced structural transformation, indicating the possibility of substantial gap modulation under extreme compression conditions.
A predicted type-II staggered band alignment with an approximately 1.4 eV valence band offset at the ZnGeN 2 /GaN heterointerface has inspired novel band-engineered III-N/ZnGeN 2 heterostructure-based device designs for applications in high performance optoelectronics. We report on the determination of the valence band offset between metalorganic chemical vapor deposition grown (ZnGe) 1- x Ga 2 x N 2 , for x = 0 and 0.06, and GaN using x-ray photoemission spectroscopy. The valence band of ZnGeN 2 was found to lie 1.45–1.65 eV above that of GaN. This result agrees well with the value predicted by first-principles density functional theory calculations using the local density approximation for the potential profile and quasiparticle self-consistent GW calculations of the band edge states relative to the potential. For (ZnGe) 0.94 Ga 0.12 N 2 the value was determined to be 1.29 eV, ~10%–20% lower than that of ZnGeN 2 . The experimental determination of the large band offset between ZnGeN 2 and GaN provides promising alternative solutions to address challenges faced with pure III-nitride-based structures and devices.
We present a novel deep learning (DL) approach to produce highly accurate predictions of macroscopic physical properties of solid solution binary alloys and magnetic systems. The major idea is to make use of the correlations between different physical properties in alloy systems to improve the prediction accuracy of neural network (NN) models. We use multitasking NN models to simultaneously predict the total energy, charge density and magnetic moment. These physical properties mutually serve as constraints during the training of the multitasking NN, resulting in more reliable DL models because multiple physics properties are correctly learned by a single model. Two binary alloys, copper–gold (CuAu) and iron–platinum (FePt), were studied. Our results show that once the multitasking NN's are trained, they can estimate the material properties for a specific configuration hundreds of times faster than first-principles density functional theory calculations while retaining comparable accuracy. We used a simple measure based on the root-mean-squared errors to quantify the quality of the NN models, and found that the inclusion of charge density and magnetic moment as physical constraints leads to more stable models that exhibit improved accuracy and reduced uncertainty for the energy predictions.
Abstract Local structures play a crucial role in the structural polyamorphism and novel electronic properties of amorphous materials, but their accurate measurement at high pressure remains a formidable challenge. In this article, we use the local structure of network-forming GeO 2 glass as an example, to present our recent approaches and advances in high-energy x-ray diffraction, high-pressure x-ray absorption fine structure, and ab initio first-principles density functional theory calculations and simulations. Although GeO 2 glass is one of the best studied materials in the field of high pressure research due to its importance in glass theory and geophysical significance, there are still some long-standing puzzles, such as the existence of appreciable distinct fivefold [5] Ge coordination at low pressure and the sixfold-plus [6+] Ge coordination at ultrahigh pressure. Our work sheds light on the origin of pressure-induced polyamorphism of GeO 2 glass, and the [5] Ge polyhedral units may be the dominant species in the densification mechanism of network-forming glasses from tetrahedral to octahedral amorphous structures.
Abstract Siliconization of the tokamak walls is a candidate method to improve plasma confinement in fusion tokamaks containing tungsten plasma facing components (W PFCs). To understand the interactions of silicon (Si) with W, the Si behavior in bulk W, and near three low-index W surfaces ((100), (110) and (111)) has been investigated using first-principles density functional theory. In bulk W, Si interstitial atoms have a low solution ability and high mobility, and Si atoms can be strongly trapped by W vacancies. The interaction between two Si adatoms is responsible for the stability of adatom superstructures on W surfaces, consistent with previous experimental observation (Tsong and Casanova 1981 Phys. Rev. Lett. 47 113). Although the coverage dependence of Si adsorption and diffusion energetics on surfaces is related to surface orientation, the W(110) surface has lower Si adsorption affinity and higher Si diffusivity than either the W(111) or W(100) surfaces. The most stable Si adatom superstructure on W surfaces is: square c(2 × 2) pattern on W(100) covered with 0.5 ML Si; rectangular c(4 × 2) pattern on W(110) with 0.25 ML Si; and rhombus p(1 × 1) pattern on W(111) with 1 ML Si. The coverage dependence of Si mobility on/toward W surfaces is generally related to the stability of the Si superstructures as a function of coverage on each surface. Interestingly, Si adatoms prefer to transport below the surface and into W subsurface by an exchange mechanism with W atoms, indicating the likelihood of epitaxial growth of W silicide layers on W surfaces during the operation of W PFCs.
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.
Abstract The unique electronic and mechanical properties of two-dimensional (2D) materials make them promising next-generation candidates for a variety of applications. Large-scale searches for high-performing 2D materials are limited to calculating descriptors with computationally demanding first-principles density functional theory. In this work, we alleviate this issue by extending and generalizing crystal graph convolutional neural networks to systems with planar periodicity and train an ensemble of models to predict thermodynamic, mechanical and electronic properties. We carry out a screening of nearly 45,000 structures for two separate applications: mechanical strength and photovoltaics. By collecting statistics of the screened candidates, we investigate structural and compositional design principles that impact the properties of the structures surveyed. Our approach recovers some well-accepted design rules: hybrid organic-inorganic perovskites with lead and tin tend to be good candidates for solar cell applications and titanium based MXenes usually have high stiffness coefficients. Interestingly, other members of the group 4 elements also contribute to increasing the mechanical strength of MXenes. For all-inorganic perovskites, we discover some compositions that have not been deeply studied in the field of photovoltaics and thus open up paths for further investigation. We open-source the code-base to spur further development in this space.
We report a systematic and accurate approach for deriving the bulk free energy surface (FES), a function of temperature, polarization, and strain, from the first-principles density functional theory (DFT) of proper ferroelectrics. The core of our approach is the metadynamics algorithm that extracts the polarization dependence of the FES from all-atom molecular dynamics simulations without an a priori ansatz. The rest of the FES is derived from the metadynamics trajectories that span the relevant phase space. We demonstrate our approach in the case of lead titanate. The errors across the phase transition, due to DFT numerics, all-atom molecular dynamics, and free energy evaluation by enhanced sampling, can be systematically controlled and are of the order of 1 meV/atom. The accuracy of the resulting ab initio FES is only limited by the adopted functional approximation of DFT.
Understanding the structure–property relationships in layered transition-metal carbides or nitrides, known as MXenes, is of critical importance for their rational design, synthesis, and application. However, the vast chemical and structural diversity of MXenes, stemming from their wide range of M and X elements, surface terminations, and different atomic coordination environments, makes it challenging to clearly understand these structure–property relationships. Here, in this work, we perform first-principles density functional theory (DFT) calculations and molecular dynamics (MD) simulations to comprehensively investigate the stability and a variety of physical properties of MXenes with different coordination environments. Using Ti- and Mo-based carbide MXenes as model systems, energetic calculations reveal that Ti-based MXenes are most stable in octahedral coordination, whereas Mo-based MXenes preferentially adopt prismatic coordination. This fundamental difference in preferred atomic coordination gives rise to markedly distinct properties between these two systems as a function of the fraction of octahedral and prismatic sites. For instance, the in-plane stiffness of Ti-based MXenes increases as octahedral coordination becomes dominant, but it decreases in the Mo-based MXenes under the same conditions. Additional stability analyses based on mechanical, lattice-dynamical, and temperature-dependent thermodynamic properties demonstrate that many metastable MXenes not only satisfy the strict stability criteria but can also undergo phase transitions among different structures and even become stabilized at elevated temperatures. Although surface terminations, such as F and O atoms, do not alter the energetic ordering or the overall stiffness trends among stable and metastable MXenes, they influence other material properties. For instance, O termination can induce semiconducting behavior in both stable and metastable Ti 2 CO 2 MXenes. This study significantly advances the fundamental understanding of structure–property relationships in MXenes and provides valuable guidance for developing coordination-based design principles to precisely engineer MXenes with improved properties.
Increased demand for high-performance permanent magnets in the electric vehicle and wind-turbine industries has prompted the search for cost-effective alternatives. Discovering magnetic materials with the desired intrinsic and extrinsic permanent magnet properties presents a significant challenge to researchers because of issues with the global supply of rare-earth elements, material stability, and a low maximum magnetic energy product BH max . While first-principles density functional theory (DFT) predicts materials’ magnetic moments, magnetocrystalline anisotropy constants, and exchange interactions, it cannot compute extrinsic properties such as coercivity (H c ). Although it is possible to calculate H c theoretically with micromagnetic simulations, the predicted value is larger than the experiment by almost an order of magnitude due to the Brown paradox. To circumvent these issues, we employ machine-learning (ML) methods on an extensive database obtained from experiments, DFT calculations, and micromagnetic modeling. The use of a large experimental dataset enables realistic H c predictions for materials such as Ce-doped Nd 2 Fe 14 B, comparing favorably against micromagnetically simulated coercivities. Remarkably, our ML model accurately identifies uniaxial magneto-crystalline anisotropy as the primary contributor to H c . With DFT calculations, we predict the Nd-site-dependent magnetic anisotropy behavior in Nd 2 Fe 14 B, confirming that Nd 4g sites mainly contribute to uniaxial magnetocrystalline anisotropy, and also calculate the Curie temperature (T c ). Finally, both calculated results are in good agreement with the experiments. The coupled experimental dataset and ML modeling with DFT input predict H c with far greater accuracy and speed than was previously possible using micromagnetic modeling. Further, we reverse engineer the grain-boundary and intergrain exchange coupling with micromagnetic simulations by employing the ML predictions.
PbZrO 3 is regarded as the first antiferroelectric and currently under intense reexamination. We propose a multiscale approach that combines classical and first principles density functional theory based simulations to explore polar phases in this material. Application of this approach to PbZrO 3 predicts three polar phases, none of which is the common $R3c$. The lowest energy polar phase - $Cc$ - is metastable at zero Kelvin but can be stabilized by application of electric field. Its structural and polarization response to the electric field is in good agreement with experimental data from the literature. The other two phases are $Ima2$ and $I4cm$ of orthorhombic and tetragonal symmetry, respectively, and predicted to stabilize at finite temperatures and under application of larger electric fields. These two phases exhibit very large strain which is of technological significance.