Antimony nitride discovered by theoretical structure prediction, rapid thermal annealing, and in situ X-ray diffraction
Not Available
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.
Not Available
The interplay of cooperative Jahn–Teller (JT) distortions and charge-disproportionation (CD) with a strong electronic correlation in transition metal oxides leads to structural symmetry breaking. Both JT and CD often manifest in the form of significant modifications in electronic and structural properties such as band splitting, metal–insulator transitions (MIT), and enhanced electron lattice interactions. Notably, the charge-disproportionation is a key electronic feature that drives the MIT. Here, we demonstrate and quantify it using first-principles calculations combining density-functional theory, dynamical mean-field theory, and spin–lattice Monte Carlo simulations. Taking Ca 2 FeMnO 6 as a prototypical example of a correlated oxide, our ab initio study shows that MIT in Ca 2 FeMnO 6 arises from the partial localization of oxygen ligand holes at alternate Fe sites that control both charge and magnetic ordering. Interestingly, the band gap was found to be fundamentally controlled by the strength of the charge-transfer energy and not by the Mott–Hubbard interactions. The novel physics and insights presented in this work reveal promising routes for tuning the electronic functionality in transition-metal oxides.
Successful synthesis of some perovskite oxyhydrides and their unique catalytic properties have recently attracted researchers’ attention. However, their surface structure remains unclear. Here we identify stable surface terminations of a prototypical perovskite oxyhydride, BaTiO 2.5 H 0.5 , under catalytically relevant temperatures and pressures by using first-principles thermodynamics based on density functional theory. Five low-index facets, including (100), (010), (210), (011), and (211), and their various terminations for a total of 47 different surfaces have been examined for relative stability at different temperatures (700, 500, 300 K) and gas environments (10 –15 ≤ P O2 ≤ 1 atm, 10–15 ≤ PH2 ≤ 100 atm). The most stable ones are found to be (010)-Ba 2 O 2 , (210)-Ti 2 O 2 , and (211)-Ba 2 O 4 H surface terminations. These polar surfaces are stabilized by charge compensation. This work provides important insights into the stable surfaces of perovskite oxyhydrides for future studies of their catalytic properties.
The natural abundance of potassium in the earth’s crust is 1000 times higher than that of lithium, so energy technologies built on potassium are more sustainable. Potassium-ion batteries have attracted considerable attention because of their relatively low cost and high operating potential, but questions remain about the best anode material for such batteries. Here, in this paper, we report first-principles computations based on density functional theory to investigate the performance of the UiO-66 metal–organic framework as an anode material for potassium-ion batteries; the goal is to provide a fundamental understanding of metal–organic framework (MOF)-based electrodes to guide the design and development of high-performance potassium-ion batteries. Our study includes the stability and electronic properties of potassiated structures and the mechanisms of potassium intercalation and diffusion in the framework. The results indicate that UiO-66 has a maximum specific capacity of 644 mAh/g as the anode of a potassium-ion battery. During potassiation, we observe charge transfer from potassium to carbon or oxygen of UiO-66 near the intercalated K. During K diffusion, the K migrates along the UiO-66 framework with a maximal migration energy barrier of 0.377 eV in the optimal pathway, which is much lower than the barriers for Li and Na diffusion in UiO-66. The diffusion coefficient of K in the anode is several orders of magnitude larger than those of Li and Na. This favors potassium ions over lithium ions or sodium ions when UiO-66 is the anode.
Zinc porphyrin solar cell dyes with donor-π–acceptor architectures combine light absorber (π), electron-donor, and electron-acceptor moieties inside a single molecule with atomic precision. The donor-π–acceptor design promotes the separation of charge carriers following optical excitation. Here, we probe the excited-state electronic structure within such molecules by combining time-resolved X-ray absorption spectroscopy at the N K-edge with first-principles time-dependent density functional theory (TD-DFT) calculations. Furthermore, customized Zn porphyrins with strong-donor triphenylamine groups or weak-donor tri-tert-butylbenzene groups were synthesized. Energetically well-separated N K-edge absorption features simultaneously probe the excited-state electronic structure from the perspectives of the macrocycle and triphenylamine N atoms. New absorption transitions between the macrocycle N atoms and the excited-state HOMO vacancy are observed, and the triphenylamine associated absorption feature blue-shifts, consistent with partial oxidation of the donor groups in the excited state.
Surface transfer doping is proposed to be a potential solution for doping diamond, which is hard to dope for applications in high-power electronics. While MoO 3 is found to be an effective surface electron acceptor for hydrogen-terminated diamond with a negative electron affinity, the effects of commonly existing oxygen vacancies remain elusive. We have performed reactive molecular dynamics simulations to study the deposition of MoO 3–x on a hydrogenated diamond (111) surface and used first-principles calculations based on density functional theory to investigate the electronic structures and charge transfer mechanisms. We find that MoO 3–x is an effective surface electron acceptor and the spatial extent of doped holes in hydrogenated diamond is extended, promoting excellent transport properties. Charge transfer is found to monotonically decrease with the level of oxygen vacancy, providing guidance for engineering of the surface transfer doping process.
Imposing additional confinement in two-dimensional (2D) materials can yield further control over the associated electronic, optical, and topological properties. However, synthesis of ultra-narrow nanoribbons (NRs) remains a challenge, particularly for the transition metal dichalcogenides (TMDs), and synthesizing TMD NRs narrower than 50 nm has remained elusive. In this paper, we report the vapor-phase synthesis of ultra-narrow TaS 2 NRs. The NRs are grown within the hollow cavity of carbon nanotubes, thereby limiting their lateral dimensions and layer number, while simultaneously stabilizing them against the environment. The NRs reach the monolayer (ML) limit and exhibit widths as low as 2.5 nm. Atomic-resolution scanning transmission electron microscopy (STEM) reveals the detailed atomic structure of the ultra-narrow NRs and we observe a hitherto unseen atomic structure supermodulation phenomenon of ordered defect arrays within the NRs. First-principles calculations based on density functional theory (DFT) show the presence of flat bands, as well as edge- and boundary-localized states, and help identify the atomic configuration of the supermodulation. Nanotube-templated synthesis represents a unique, transferable, and broadly deployable route toward ultra-narrow TMD NR growth.
Low-cost transition-metal oxide is regarded as a promising electrocatalyst family for an oxygen evolution reaction (OER). The classic design principle for an oxide electrocatalyst believes that point defect engineering, such as oxygen vacancies (VO .. ) or heteroatom doping, offers the opportunities to manipulate the electronic structure of material toward optimal OER activity. Oppositely, in this work, we discover a counterintuitive phenomenon that both VO .. and an aliovalent dopant (i.e., proton (H + )) in perovskite nickelate (i.e., NdNiO 3 (NNO)) have a considerably detrimental effect on intrinsic OER performance. Detailed characterizations unveil that the introduction of these point defects leads to a decrease in the oxidative state of Ni and weakens Ni-O orbital hybridization, which triggers the local electron-electron correlation and a more insulating state. Evidenced by first-principles calculation using the density functional theory (DFT) method, the OER on nickelate electrocatalysts follows the lattice oxygen mechanism (LOM). The incorporation of point defect increases the energy barrier of transformation from OO*(V O ) to OH*(V O ) intermediates, which is regarded as the rate-determining step (RDS). Furthermore, this work offers a new and significant perspective of the role that lattice defects play in the OER process.
Single-component halide perovskites (HPs) rarely satisfy all the necessary criteria for optoelectronic applications, such as achieving an optimal band gap while maintaining high chemical and structural stability. Alloying halide perovskites has emerged as a promising strategy, not only to enhance stability but also to fine-tune their electronic and optical properties. In this work, we explore multiple degrees of freedom in alloy design, considering different substitution sublattices sites (A, B, or X in ABX3 perovskites), various chemical species (isovalent and hetero-valent elements), and multi-component compositions on a given sublattice. Using first-principles calculations based on density functional theory (DFT), we investigate how compositional variations influence the electronic (band gap) and structural properties (mixing enthalpy) of HP alloys. Our approach employs the polymorphous cell model, allowing full local relaxation which breaks local symmetry while preserving global cubic symmetry—an essential framework for accurately modeling HPs. Our results reveal that X-site mixing (halogen substitution) primarily affects the valence band maximum, allowing target band gap engineering. Additionally, variations in halogen radii introduce internal strain through octahedral distortions, influencing the mixing enthalpy. A-site substitution, while not directly contributing to the band edge states, modifies structural stability via volume effects, indirectly impacting the band gap. B-site alloying plays a dominant role in band gap modulation, leading to either positive or negative band gap bowing. Specifically, isovalent B-site mixing (Sn–Pb) induces strong positive bowing, where the alloy band gap is smaller than the average gap of parent compounds, whereas hetero-valent mixing (Cd–Pb) results in pronounced negative bowing. As an aside, we investigate the competition between the excess energy of disordered alloys vs. that of long-range ordered double perovskites of the same compositions, seeking examples of ordered phases emerging from disordered alloys. Furthermore, our findings provide fundamental insights into the electronic and structural behavior of HP alloys, offering valuable design principles for the development of stable and efficient materials for next-generation photovoltaic and optoelectronic devices.
The II-IV-N2 class of heterovalent ternary nitrides has gained significant interest as alternatives to the III-nitrides for electronic and optoelectronic applications. In this study, we apply first-principles calculations based on density functional theory to systematically investigate the effects of structural distortions due to cation size mismatch on the configurational disorder of the cation sublattice and the valence band structure in this class of materials. We find that larger size mismatch between the group-II and the group-IV cations results in stronger lattice distortions from the ideal hexagonal ratio, which in turn inhibits the propensity of these materials toward octet-rule violating cation disorder. We also demonstrate that the formation energy of a single cation antisite pair, which is fast and simple to calculate, is a strong indicator of a material's propensity toward disorder. Furthermore, the breaking of in-plane symmetry leads to a splitting of the top three valence bands at Γ, which is also directly related to the magnitude of structural distortions. Our work demonstrates that the structural and functional properties of the II-IV-N2 materials can be finely tuned through controllable structural distortions that stem from the choice of cations.
Simulations and diagnostics of high-energy-density plasmas and warm dense matter rely on models of material response properties, both static and dynamic (frequency-dependent). Here, in this work, we systematically investigate variations in dynamic electron–ion collision frequencies ν(ω) in warm dense matter using data from a self-consistent-field average-atom model. We show that including the full quantum density of states, strong collisions, and inelastic collisions lead to significant changes in ν(ω). These changes result in red shifts and broadening of the plasmon peak in the dynamic structure factor, an effect observable in x-ray Thomson scattering spectra, and modify stopping powers around the Bragg peak. These changes improve the agreement of computationally efficient average-atom models with first-principles time-dependent density functional theory in warm dense aluminum, carbon, and deuterium.
Many important chemical processes involve reactivity and dynamics in complex solutions. Gaining a fundamental understanding of these reaction mechanisms is a challenging goal that requires advanced computational and experimental approaches. However, important techniques such as molecular simulation have limitations in terms of scales of time, length, and system complexity. Furthermore, among the currently available solvation models, there are very few designed to describe the interaction between the molecular scale and the mesoscale. To help address this challenge, here, we establish a novel hybrid approach that couples first-principles plane-wave density functional theory with classical density functional theory (cDFT). In this approach, a region of interest described by ab initio molecular dynamics (AIMD) interacts with the surrounding medium described using cDFT to arrive at a self-consistent ground state. cDFT is a robust but efficient mesoscopic approach to accurate thermodynamics of bulk electrolyte solutions over a wide concentration range (up to 2M concentrations). Benchmarking against commonly used continuum models of solvation, such as SMD, as well as experiments, demonstrates that our hybrid AIMD–cDFT method is able to produce reasonable solvation energies for a variety of molecules and ions. With this model, we also examined the solvent effects on a prototype S N 2 reaction of the nucleophilic attack of a chloride ion on methyl chloride in the solution. The resulting reaction pathway profile and the solution phase barrier agree well with experiment, showing that our AIMD/cDFT hybrid approach can provide insight into the specific role of the solvent on the reaction coordinate.
The massive amount of wasted heat energy from industry has pushed the development of thermoelectric (TE) materials that directly convert heat into electricity to a new level of concern. Recently, multicomponent alloys such as GeTe-based and PbSe-based high-entropy (HE) chalcogenides have attracted a great deal of attention due to their potential application as TE materials. The nature of the interatomic bonding, lattice distortion (LD), and the electronic structure in this class of materials is not fully understood. Herein, we report a comprehensive computational investigation of nine GeTe-based HE alloys with eight metallic elements (Ag, Pb, Sb, Bi, Cu, Cd, Mn, and Sn) with large supercells of 1080 atoms each; seven PbSe-based HE solid solutions: Pb0.99−ySb0.012SnySe1−2xTexSx (x = 0.1, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45, with y = 0) with supercells of 1000 atoms each; and five Pb0.99−ySb0.012SnySe1−2xTexSx (y = 0.05, 0.1, 0.15, 0.2, 0.25 with x = 0.25) solid solutions with supercells of 1000 atoms each. All these HE models are theoretically investigated for the first time. The electronic structure, interatomic bonding, charge transfer, and lattice distortion (LD) are investigated by first-principles calculations based on density functional theory. Multicomponent HE alloys can cause a significant LD, which affects their mechanical, thermal, and TE properties. The calculations for the GeTe-based HE chalcogenides showed that they are semiconductors with a narrow bandgap, except for m8, which has a semi-metallic characteristic, and this makes them good candidates for TE applications. For most of these models, the Fermi level shifts upward and locates deeply in the conduction bands, resulting in the enhancement of the electrical conductivity (σ). The bonding properties showed that most bonds in m5 are more dispersed, indicating highest LD and lower lattice thermal conductivity. For PbSe-based HE solid solutions, the LD calculations showed that the models Pb0.99Sb0.012Se0.5Te0.25S0.25 and Pb0.89Sb0.012Sn0.1Se0.5Te0.25S0.25 have the higher LD, and thus a lower lattice thermal conductivity. Such investigations are in high demand since it enables us to design new HE chalcogenides for TE applications. We use the novel concept of total bond order density as a single quantum mechanical metric to characterize the internal cohesion of these HE alloys and correlate with calculated properties, especially the mechanical properties. This work provides a solid database for HE chalcogenides and a road map for many potential applications. Moreover, the computational procedure we developed can be used to design new HE chalcogenides for specific TE applications.
Discovery of structure-property relationships in thin film alloys of complex metal oxides enabled by high-throughput materials synthesis and characterization facilities is demonstrated here with a case-study. In this study, thin films of binary transition metal oxides (Ti–Zn) are prepared by pulsed laser deposition with continuously varying Ti:Zn ratio, creating combinatorial samples for exploration of the properties of this material family. The atomic structure and electronic properties are probed by spatially resolved techniques including x-ray absorption near edge structures (XANES) and x-ray fluorescence (XRF) at the Ti and Zn K-edge, x-ray diffraction, and spectroscopic ellipsometry. The observed properties as a function of Ti:Zn ratio are resolved into mixtures of five distinguishable phases by deploying multivariate curve resolution analysis on the XANES spectral series, under constraints set by results from the other characterization techniques. First-principles computations based on density function theory connect the observed properties of each distinct phase with structural and spectral characteristics of crystalline polymorphs of Ti–Zn oxide. Continuous tuning of the optical absorption edge as a function of Ti:Zn ratio, including the unusual observation of negative optical bowing, exemplifies a functional property of the film correlated to the phase evolution.
Fully self-consistent GW (sc-GW) methods are now available to evaluate quasiparticle and spectral properties of various molecular and bulk systems. However, such techniques based on the full matrix of G and W are computationally demanding. Additionally, the routinely used single-shot GW approximation (G 0 W 0 ) has an undesirable dependency on the choice of initial exchange-correlation functional. In the literature, many so-called self-consistent GW methods are based on diagonal approximation of G and low-ranking approximation of W. It is thus worth checking how good such approximations are in comparison with the full matrix method. In this work, we consider AlAs, AlP, GaP, and ZnS as the prototype systems to perform sc-GW calculations by expressing the full G matrix using a plane-wave basis set. We compared our sc-GW results with the diagonal G and subspace W approximated sc-GW results (sc-GW-diagG and sc-GW-subW methods). In the sc-GW-diagG method, interacting G is expanded in the eigenvectors of noninteracting G such that only diagonal elements are retained, whereas the number of eigenmodes is truncated in sc-GW-subW calculations. A systematic analysis of the results obtained from the above techniques is presented. The differences in the quasiparticle band gap between the approximated and the full matrix sc-GW approaches are mostly less than 1.7%, which validates such widely adopted approximations, and also shows how such low-ranking approximation can be used to include higher-order terms such as the vertex correction.
The recent discovery of superconductivity in Sr-doped NdNiO 2 films grown on SrTiO 3 started a novel field within unconventional superconductivity. Therefore, to understand the similarities and differences between nickelate and cuprate layers on the same SrTiO 3 substrate, here based on the density functional theory we have systematically investigated the structural, electronic, and magnetic properties of NdNiO 2 / SrTiO 3 and CaCuO 2 / SrTiO 3 systems. Our results revealed a strong lattice reconstruction in the case of NdNiO 2 / SrTiO 3 , resulting in a polar film, with the surface and interfacial NiO 2 layers presenting opposite displacements. To avoid the “polar catastrophe,” the NiO 2 surface to the vacuum reconstructs as well. However, for CaCuO 2 / SrTiO 3 , the distortions of those same two CuO 2 layers were in the same direction. In addition, we found this distortion to be approximately independent of the studied range of film thickness for the nickelate films. Furthermore, we also observed a two-dimensional electron gas at the interface between NdNiO 2 and SrTiO 3 , caused by the polar discontinuity, in agreement with recent literature. For NdNiO 2 / SrTiO 3 the two-dimensional electron gas extends over several layers, while for CaCuO 2 / SrTiO 3 this electronic rearrangement is very localized at the interface between CaCuO 2 and SrTiO 3 . The electronic reconstruction found at the interface involves a strong occupation of the Ti 3 d x y state. In both cases, there is a significant electronic charge transfer from the surface Ni or Cu layers to the Ti interface layer. The interfacial Ni and Cu layers are hole and electron doped, respectively. By introducing magnetism and electronic correlation, we observed that the d 3 z 2 - r 2 orbital of Ni becomes itinerant while the same orbital for Cu remains doubly occupied, establishing a clear two- vs one-orbital active framework for the description of these systems. Furthermore, we also observed a strong magnetic reconstruction at the NdNiO 2 surface to vacuum layer where magnetism is basically suppressed.
The thermodynamic properties of hexagonal-close-packed iron (ε–Fe) are essential for investigating the internal structure and dynamic properties of planetary cores. Despite their importance to planetary sciences, experimental investigations of ε–Fe at relevant conditions are still challenging. Therefore, ab initio calculations are crucial to elucidating the thermodynamic properties of this system. Here, we use a free energy calculation scheme based on the phonon gas model compatible with temperature-dependent phonon frequencies. We investigate the effects of electronic thermal excitations, which introduces a temperature dependence on phonon frequencies, and the implication for the thermodynamic properties of ε–Fe at extreme pressure (P) and temperature (T) conditions. We disregard phonon-phonon interactions, i.e., anharmonicity and their effect on phonon frequencies. Nevertheless, the current scheme is also applicable to T -dependent anharmonic frequencies. We conclude that the impact of thermal electronic excitations on vibrational properties is not significant up to ~4000 K at 200 GPa but should not be ignored at higher temperatures or pressures. However, the static free energy F st must always include the effect of thermal excitation fully in a continuum of T. Furthermore, our results for isentropic equations of state show good agreement with data from recent ramp compression experiments up to 1400 GPa conducted at the National Ignition Facility.
A Kohn-Sham density functional approach has recently been developed for the fractional quantum Hall effect, which maps the strongly interacting electrons into a system of weakly interacting composite fermions subject to an exchange correlation potential as well as a density dependent gauge field that mimics the "flux quanta" bound to composite fermions. To get a feel for the role of various terms, we study the behavior of the self-consistent solution as a function of the strength of the exchange correlation potential, which is varied through an ad hoc multiplicative factor. Here, we find that a crystal phase is stabilized when the exchange correlation interaction is sufficiently strong relative to the composite-fermion cyclotron energy. Various properties of this crystal are examined.