SEARCH · Engineering Papers

Results for “first principles density functional theory”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,027 records · Page 57

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Advanced thermal/environmental barrier coatings of high-entropy rare earth disilicates tuned by strong anharmonicity of Eu 2 Si 2 O 7

Advancing thermal/environmental barrier coating (TEBC) materials with integrated thermal-mechanical functions is paramount for safeguarding SiC-based ceramic matrix composites (CMCs) in high-efficiency gas turbines. Herein, we employ a synergistic approach, combining density functional theory (DFT) methods and combinatorial chemistry techniques, to design high-performance and low-cost RE 2 Si 2 O 7 (RE = rare earth elements) TEBC materials tailored for enhanced compatibility with SiC-based CMCs. Expanding on phase stability of alloying pure RE 2 Si 2 O 7 , the investigation extends to the mechanical and thermal properties of solid solution systems, including Er 1/2 Y 3/4 Yb 3/4 Si 2 O 7 , Gd 1/4 Er 1/4 Y 3/4 Yb 3/4 Si 2 O 7 , and Eu 1/4 Er 1/4 Y 3/4 Yb 3/4 Si 2 O 7 . The solid solution systems exhibit a major reduction in lattice thermal conductivity relative to their pure counterparts, achieving ultralow values of 0.25 to 0.39 W m −1 K −1 at 1500 K. Furthermore, the coefficients of thermal expansion (CTE) of these solid solutions are precisely tuned within the desired range for SiC (4.4 to 5.5 × 10 −6 K −1 ), while maintaining good mechanical properties. Here, in particular, the addition of Eu 2 Si 2 O 7 demonstrates to be an important variable to the tuning of CTE and lattice thermal conductivity by leveraging its strong anharmonicity, presenting a pioneering avenue for fine-tuning material properties. In summary, this research not only identifies promising TEBC materials with superior thermal properties, but also introduces a valuable computational material design methodology for the rapid discovery of complex materials for harsh environments.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Anisotropy, frustration, and saddle point in the twisted kagome compound ErPdPb

The kagome lattice, with its inherent geometric frustration, provides a rich platform for exploring intriguing magnetic phenomena and topological electronic structures. In reduced-symmetry structures, such as twisted kagome systems involving rare earth elements, additional anisotropy can arise, enabling intriguing properties including spin-ice states, magnetocaloric effects, noncollinear magnetic ordering, and the anomalous Hall effect. Here, we report the synthesis of single crystals of ErPdPb, which features a twisted kagome lattice net of Er atoms within the hexagonal ZrNiAl-type structure, and we investigate its magnetic, electronic, and thermal properties. The material exhibits a highly anisotropic, correlated magnetic state below 2.7 K, as evidenced by magnetic, transport, and heat capacity measurements. Density functional theory (DFT) calculations indicate strong easy-axis anisotropy, consistent with experiment and crystal-field expectations, as well as quasi-one-dimensional bands and a spin-split saddle point at the zone center. Here, the coexistence of competing magnetic interactions along different crystallographic directions suggests an inherent degree of frustration. ErPdPb thus provides a promising platform for exploring the interplay of frustration, anisotropy, and electronic structure in a twisted kagome lattice.

Antiferromagnets↗

DOE OSTI2021

Chemically and electrically tunable spin polarization in ferroelectric Cd-based hybrid organic-inorganic perovskites

Density functional theory computations are used to predict the electronic structure of Cd-based hybrid organic-inorganic high-T C ferroelectric perovskites with TMCM-CdCl 3 being one representative. Here we report Rashba-Dresselhaus spin splitting in the valence band of these nonmagnetic compounds. Interestingly, we find in computations that the splitting is not necessarily sensitive to the polarization of the material but to the organic molecule itself which opens a way to its chemical tunability through the choice of the molecule. Further chemical tunability of spin splitting is shown to be possible through a substitution of Cl in the CdCl 3 chains as the valence band was found to originate from Cl-Cl weekly bonding orbitals. For example, the substitution of Cl with Br in TMCM-CdCl 3 resulted in a ten times increase of spin splitting. Furthermore, the spin polarization in these materials gives origin to persistent spin textures which are coupled to the polarization direction, and, therefore, can be controlled by the electric field. This is promising for spintronics applications.

36 MATERIALS SCIENCE↗

DOE OSTI2021

J eff = $\frac{3}{2}$ metallic phase and unconventional superconductivity in GaTa 4 Se 8

By means of density functional theory plus dynamical mean-field theory (DFT+DMFT) calculations and resonant inelastic x-ray scattering (RIXS) experiments, we investigate the high-pressure phases of the spin-orbit-coupled Jeff = 3/2 insulator GaTa 4 Se 8 . Its metallic phase, derived from the Mott state by applying pressure, is found to carry J eff = 3/2 moments. The characteristic excitation peak in the RIXS spectrum maintains its destructive quantum interference of J eff at the Ta L2-edge up to 10.4 GPa. Our exact diagonalization based DFT+DMFT calculations including spin-orbit coupling also reveal that the J eff = 3/2 character can be clearly identified under high pressure. These results establish the intriguing nature of the correlated metallic magnetic phase, which represents the first confirmed example of J eff =3/2 moments residing in a metal. They also indicate that the pressure-induced superconductivity is likely unconventional and influenced by these J eff = 3/2 moments. Based on a self-energy analysis, we furthermore propose the possibility of doping-induced superconductivity related to a spin-freezing crossover.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

High Thermoelectric Performance and Defect Energetics of Multipocketed Full Heusler Compounds

We report a first-principles density-functional study of electron-phonon interactions in and thermoelectric transport properties of the full Heusler compounds Sr 2 BiAu and Sr 2 SbAu. Our results show that ultrahigh intrinsic bulk thermoelectric performance across a wide range of temperatures is physically possible and point to the presence of multiply degenerate and highly dispersive carrier pockets as the key factor for achieving this. Sr 2 BiAu, which features ten energy-aligned low-effective-mass pockets (six along Γ-X and four at L), is predicted to deliver n-type zT=0.4-4.9 at T=100-700 K. Comparison with the previously investigated compound Ba 2 BiAu shows that the additional L pockets in Sr 2 BiAu significantly increase its low-temperature power factor to a maximum value of 12 mW m -1 K -2 near T=300 K. However, at high temperatures the power factor of Sr 2 BiAu drops below that of Ba 2 BiAu because the L states are heavier and subject to strong scattering by phonon deformation, as opposed to the lighter Γ-X states, which are limited by polar-optical scattering. Sr2SbAu is predicted to deliver a lower n-type zT=3.4 at T=750 K due to appreciable misalignment between the L and Γ-X carrier pockets, generally heavier scattering, and a slightly higher lattice thermal conductivity. Soft acoustic modes, which are responsible for the low lattice thermal conductivity, also increase the vibrational entropy and high-temperature stability of these Heusler compounds, suggesting that their experimental synthesis may be feasible. The dominant intrinsic defects are found to be Au vacancies, which drive the Fermi level towards the conduction band and work in favor of n-doping.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Structural and mechanical properties of monolayer amorphous carbon and boron nitride

Amorphous materials exhibit various characteristics that are not featured by crystals and can sometimes be tuned by their degree of disorder (DOD). Here, we report results on the mechanical properties of monolayer amorphous carbon (MAC) and monolayer amorphous boron nitride (maBN) with different DOD. The pertinent structures are obtained by kinetic-Monte-Carlo (kMC) simulations using machine-learning potentials (MLP) with density-functional-theory (DFT)-level accuracy. An intuitive order parameter, namely the areal fraction F x occupied by crystallites within the continuous random network, is proposed to describe the DOD. We find that F x captures the essence of the DOD: Samples with the same F x but different sizes and arrangements of crystallites, obtained using two distinct kMC procedures, have virtually identical radial distributions functions as well as bond-length and bond-angle distributions. Furthermore, by simulating the fracture process with molecular dynamics, we found that the mechanical responses of MAC and maBN before fracture are mainly determined by F x and are insensitive to the sizes and specific arrangements and to some extent the numbers and area distributions of the crystallites. The behavior of cracks in the two materials is analyzed and found to mainly propagate in meandering paths in the CRN region and to be influenced by crystallites in distinct ways that toughen the material. Furthermore, the present results reveal the relation between structure and mechanical properties in amorphous monolayers and may provide a universal toughening strategy for 2D materials.

2-dimensional systems↗

DOE OSTI2020

Recent developments in the PySCF program package

PySCF is a Python-based general-purpose electronic structure platform that supports first-principles simulations of molecules and solids as well as accelerates the development of new methodology and complex computational workflows. Here, we explain the design and philosophy behind PySCF that enables it to meet these twin objectives. With several case studies, we show how users can easily implement their own methods using PySCF as a development environment. We then summarize the capabilities of PySCF for molecular and solid-state simulations. Finally, we describe the growing ecosystem of projects that use PySCF across the domains of quantum chemistry, materials science, machine learning, and quantum information science.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Computational Fermi Level Engineering and Doping-Type Conversion of Mg:Ga 2 O 3 via Three-Step Synthesis Process

Gallium oxide (Ga 2 O 3 ) is being actively explored for electronics that can operate at high power, temperature, and frequency as well as for deep-ultraviolet optoelectronics and other applications due to its ultra-wide bandgap (UWBG) and low projected fabrication cost of large-size and high-quality crystals. Efficient n-type doping of monoclinic beta-phase of Ga 2 O 3 has been achieved, but p-type doping faces fundamental obstacles due to compensation, deep acceptor levels, and the polaron transport mechanism of free holes. However, aside from the challenges of achieving p-type conductivity, plenty of opportunity exists to engineer the position of the Fermi level for improved design of Ga 2 O 3 -based devices. We use first-principles defect theory and defect equilibrium calculations to simulate a three-step growth-annealing-quench synthesis protocol for hydrogen-assisted Mg doping in ß-Ga 2 O 3 . The simulations take into account the gas phase equilibrium between H 2 , O 2 , and H 2 O, which determines the H chemical potential. We predict Ga 2 O 3 doping-type conversion to a net p-type regime after growth under reducing conditions in the presence of H2 followed by O-rich annealing, which is a similar process to Mg acceptor activation by H removal in GaN. For equilibrium annealing with re-equilibration of compensating O vacancies, there is an optimal temperature that maximizes the Ga 2 O 3 net acceptor density for a given Mg doping level; the acceptor density is further increased in the non-equilibrium annealing scenario without re-equilibration. After quenching to operating temperature, the Ga 2 O 3 Fermi level drops below mid-gap down to about 1.5 eV above the valence band maximum, creating a significant number of uncompensated neutral Mg Ga 0 acceptors. The resulting free hole concentration in Ga 2 O 3 is very low even at elevated operating temperature (~10 8 cm -3 at 400°C) due to the deep energy level of these Mg acceptors, and hole conductivity is further impeded by the polaron hopping mechanism. However, the Fermi-level reduction and suppression of free electron density in this doping-type converted (NA>ND) Ga 2 O 3 material are important for improved designs of Ga 2 O 3 electronic devices. These results illustrate the power of computational predictions not only for new materials but also for their synthesis science.

36 MATERIALS SCIENCE↗

DOE OSTI2020

A First–Principles–Based Sub–Lattice Formalism for Predicting Off–Stoichiometry in Materials for Solar Thermochemical Applications: The Example of Ceria

Theoretical models that reliably can predict off-stoichiometry in materials via accurate descriptions of underlying thermodynamics are crucial for energy applications. For example, transition-metal and rare-earth oxides that can tolerate a large number of oxygen vacancies, such as CeO 2 and doped CeO 2 , can split water and carbon dioxide via a two-step, oxide-based solar thermochemical (STC) cycle. The search for new STC materials with a performance superior to that of state-of-the-art CeO 2 can benefit from predictions accurately describing the thermodynamics of oxygen vacancies. The sub-lattice formalism, a common tool used to fit experimental data and build temperature-composition phase diagrams, can be useful in this context. Here, sub-lattice models are derived solely from zero-temperature quantum mechanics calculations to estimate fairly accurate temperature- and oxygen-partial-pressure-dependent off-stoichiometries in CeO 2 and Zr-doped CeO 2 . Physical motivations for deriving some of the “excess” sub-lattice model parameters directly from quantum mechanical calculations, instead of fitting to minimize deviations from experimental and/or theoretical data, are identified. As a result, important limitations and approximations of the approach used are specified and extensions to multi-cation oxides are also suggested to help identify novel candidates for water and carbon dioxide splitting and related applications.

08 HYDROGEN↗

DOE OSTI2023

Accuracy of DFT computed oxygen-vacancy formation energies and high-throughput search of solar thermochemical water-splitting compounds

The enthalpy change involved in metal oxide reduction is a key quantity in various processes related to energy conversion and storage, and is of particular interest for computational prediction. Often this prediction involves the simulation of a high temperature reduction process with a 0K methodology like density functional theory (DFT), and it is not infrequent for the high temperature and 0K stable crystal structures to differ. This introduces a conundrum with regards to the choice of crystal structure to utilize in the computation, with approaches in the literature varying and experimental validation remaining scarce. In this work we address both the crystal structure conundrum and the experimental validation, and then apply the insights we gain to guide a high-throughput search for new materials for solar thermochemical water-splitting applications. By computing the DFT+U oxygen vacancy formation energy (ΔE vf ) of a selection of ABO 3 compounds and comparing different crystal structures for each composition, we highlight the issues that arise when the structure utilized in the computation is dynamically unstable at 0K, namely the presence of an artificial lowering of ΔE vf , and the lack of convergence of ΔE vf with cell size. We solve these limitations by identifying and employing a suitable surrogate dynamically stable structure. We then validate the predictive power of our calculations against appositely generated experimental measurements of reduction enthalpy for a series of Hubbard U values, finding an accuracy ranging between 0.2-0.6 eV/O. In light of such conclusions, we revise and expand a previous a high-throughput DFT study on ABO 3 perovskite oxides. As a result, we provide a list of candidate STCH materials, highlight trends with redox-active cation and structural distortion, and identify Mn 4+ , Mn 3+ and Co 3+ as the most promising redox-active cations.

08 HYDROGEN↗

DOE OSTI2023

Symmetry-driven persistent spin texture for the two-dimensional nonsymmorphic CdTe and ZnTe crystal structures

In this paper, two nonsymmorphic two-dimensional structures of CdTe and ZnTe are modeled, and using state-of-the-art density functional theory with the group theory of solids, their symmetry-enforced electronic properties are studied. The in-plane ferroelectricity coupled with strong spin-orbit coupling induces a unidirectional out-of-plane Rashba spin-orbit field that can host a momentum-independent uniform spin configuration known as persistent spin texture (PST) at the Brillouin zone center. PST in these structures is found to be robust against external perturbations such as strain, structural distortion, and independent of layer thickness. These unprecedented intrinsic spin transport properties hold utmost importance in spintronics, as the experimental stringent condition of equal Rashba and Dresselhaus constants [Phys. Rev. Lett. 90, 146801 (2003)] is eliminated. The calculated persistent spin helix wavelength of <~5 nm paves the way for developing next-generation nanosized nonballistic spin field-effect transistors compared with micrometer-sized GaAs/AlGaAs quantum wells. Further, these materials exhibit finite spin Hall conductivity at the band edges and hence can be used in ferromagnet-free spin Hall transistors. Although CdTe and ZnTe systems have been widely studied for photocatalysis and solar cell applications over the past few decades, their potential application in spintronic devices has not been explored. Mono/few layers of CdTe and ZnTe synthesized from (110) facets of bulk zinc-blende crystals [Nat. Commun. 3, 1057 (2012)] satisfy all symmetry operations of the nonsymmorphic space group and hence can be considered ideal materials to verify our theoretical results experimentally.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Facet-dependent structure and dissociation of water at pristine IrO 2 /water interfaces

Understanding the microscopic structure of water at metal oxide interfaces is crucial for advancing electrocatalysis. IrO 2 , specifically, has shown exceptional activity for electrochemical water oxidation, but we currently lack a fundamental understanding of how the surface structure of IrO 2 impacts water reactivity. In this work, we developed a machine learning potential trained to first-principles accuracy for modeling IrO 2 /water interfaces across different facets: (110), (100), (101), and (001). Using extensive machine learning molecular dynamics simulations, we investigated the spontaneous dissociation of water molecules at these interfaces. Our results reveal a distinct dissociation probability trend: (110) > (100) ≈ (101) > (001), which we attribute primarily to the reaction thermodynamics of surface water dissociation. A strong correlation is observed between the surface Ir–O bond distances and the dissociation probabilities, highlighting the role of surface geometry in modulating reactivity. As a consequence, the interfacial solvation structures and hydrogen bonding environments are dynamically tuned by the varying water dissociation capabilities across facets. This work elucidates how water dissociation energetics depend on surface orientation and interfacial structure, offering atomistic insights into manipulating reaction chemistry at electrocatalytic interfaces.

organic↗

DOE OSTI2022

Accurate Prediction of Voltage of Battery Electrode Materials Using Attention-Based Graph Neural Networks

Performing first-principles calculations to discover electrodes’ properties in the large chemical space is a challenging task. While machine learning (ML) has been applied to effectively accelerate those discoveries, most of the applied methods ignore the materials’ spatial information and only use predefined features: based only on chemical compositions. Here, we propose two attention-based graph convolutional neural network techniques to learn the average voltage of electrodes. Our proposed methods, which combine both atomic composition and atomic coordinates in 3D-space, improve the accuracy in voltage prediction significantly when compared to composition-based ML models. The first model directly learns the chemical reaction of electrodes and metal ions to predict their average voltage, whereas the second model combines electrodes’ ML predicted formation energy (E form ) to compute their average voltage. Our E form -based model demonstrates improved accuracy in transferability from our subset of learned Li ions to Na ions. Moreover, we predicted the theoretical voltage of 10 Na x MPO 4 F (M = Ti, Cr, Fe, Cu, Mn, Co, and Ni) fluorophosphate battery frameworks, which are unavailable in the Material Project database. It could be shown that we can expect average voltages higher than 3.1 V from those Na battery frameworks except from the NaTiPO 4 F and TiPO 4 F pair of electrodes, which offer an average voltage of 1.32 V.

25 ENERGY STORAGE↗

DOE OSTI2026

Lattice vacancy migration barriers in Fe-Ni alloys, and an indication as to why Ni atoms diffuse slowly: A first-principles study

Lattice vacancy migration barriers in ferromagnetic Fe 𝑥 ⁢Ni 1−𝑥 alloys (0.4 ≤ 𝑥 ≤ 0.6) are accurately quantified within the framework of ab initio electronic structure calculations using the nudged elastic band (NEB) method. Both the atomically disordered (A1) fcc phase, as well as the atomically ordered, tetragonal L⁢1 0 phase—which is under consideration as a material for a rare-earth-free gap magnet for advanced engineering applications—are investigated. Across an ensemble of NEB calculations performed on supercell configurations spanning a range of compositions and containing disordered, partially ordered, and fully ordered structures, we find that Ni-vacancy interchanges encounter significantly higher energetic barriers than do Fe-vacancy interchanges. We contend that this aspect is a key factor in determining the differences in mobility between Fe and Ni atoms in this ferromagnetic alloy. Moreover, we are able to interpret these findings in terms of the ferromagnetic alloy's underlying spin-polarized electronic structure. Specifically, we report a coupling between the size of local lattice distortions and the magnitude of the local electronic spin polarization around vacancies. This causes Fe atoms to relax into lattice vacancies, while Ni atoms remain rigidly fixed to their original lattice positions. These results give atomic-scale insight into the longstanding experimental observation that Ni exhibits remarkably slow atomic diffusion in Fe-Ni alloys.

density functional theory↗

DOE OSTI2025

Defect-induced phonon-resonant scattering and its influence on thermal transport of irradiated thorium-dioxide

Thermal transport in proton irradiated thorium-dioxide (ThO 2 ) is investigated. Using a combination of experiments and first-principles computational framework, the role of lattice defects on thermal conductivity is analyzed. A resonant-phonon scattering mechanism beyond the traditionally considered Rayleigh scattering is found to significantly influence low-temperature thermal transport in the presence of irradiation-induced point defects. The existence of localized phonon modes associated with irradiation-induced defects is suggested by the inability of the first-principles based thermal conductivity model—which considers only three-phonon interactions and phonon-defects scattering using the Tamura formalism—to predict the experimental results, unless a resonant scattering mechanism is included. The emergence of additional peaks in the Raman spectra in the proximity of phonon-resonant frequency provides further evidence for the existence of localized modes. Coupled with a microstructure evolution model, this analysis enables more accurate analysis for contrasting the contributions of different phonon scattering mechanisms across all irradiation doses and temperatures.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

Mechanically and electrically tunable Rashba-Edelstein effect in ferroelectric semiconductors, CsGe⁢𝑋 3 (𝑋 = I, Br, Cl)

The ability of materials to convert charge current into spin current is fundamental to many spintronics applications. One means of realizing this conversion is via Rashba-Edelstein effect (REE). Using density functional theory simulations, we predict that REE can be induced in the recently discovered family of semiconducting ferroelectrics, CsGe⁢𝑋 3 (𝑋 = I, Br, Cl). The effect is quantified through Rashba-Edelstein coefficients, 𝜒 𝑥⁢𝑦 =−𝜒 𝑦⁢𝑥 , which are nonzero in valence and conduction bands. The largest values, obtained for CsGeI 3 , are 3.45 × 10 10 $\frac{ℏ}{Λ cm}$ and 0.97 × 10 10 $\frac{ℏ}{Λ cm}$ in the conduction and valence bands, respectively. The values are comparable to, and sometimes exceed, those of other inorganic materials, although the maximal values occur away from the band edges. The coefficients' sign couples to the direction of spontaneous polarization, offering opportunities for nonvolatile spin current manipulation via external electric field. Furthermore, these coefficients are highly tunable through strain engineering owing to strain-induced variations in spin textures and energy isosurfaces. Specifically, 𝜒 in the valence band of CsGeI 3 is enhanced to 3.61 × 10 10 $\frac{ℏ}{Λ cm}$ under 5% biaxial strain. Given this potential, we synthesized CsGeI 3 to validate our simulation structure and found excellent agreement between experiment and simulation, thereby allowing extrapolation of our findings to practically significant temperatures. Furthermore, our study identifies promising materials for nonvolatile, multifunctional spintronic applications.

Electronic structure↗

DOE OSTI2025

Stability, electronic quantum states, and magnetic interactions of Er 3+ ions in Ga 2 ⁢O 3

Here, we report an ab initio study of phase stability, defect formation, electronic structure, and multiple magnetic, Dzyaloshinskii-Moriya, optical, hyperfine, and crystal field interactions in erbium (Er)-doped wide band gap 𝛼- and 𝛽-gallium oxides (Ga 2 ⁢O 3 ), critically important to make a foundation for both optoelectronic and quantum information applications. The chemical, structural, mechanical, and dynamical stabilities of the pristine phases are confirmed from respective negative formation energies, negative cohesive energies, favorable elastic constants, and positive phonon frequencies. The phonon dispersions indicate that the Ga-O bonds are uniform in the 𝛼-phase, while they vary in the 𝛽-phase due to the anisotropic polyhedral movement. The defect formation energy analysis confirms that both Er-doped 𝛼- and 𝛽−Ga 2 ⁢O 3 prefer Er 3+ (neutral) state. The underestimated band gaps of the pristine phases from standard density functional theory (DFT) calculations as compared to experimental values are corrected by employing the hybrid functional calculations, resulting in the indirect band gaps of 5.21 eV in 𝛼−Ga 2 ⁢O 3 and 4.94 eV in 𝛽−Ga 2 ⁢O 3 . The site preference energy analysis indicates partial occupation of Er in the octahedral site of Ga. The anisotropic nature of hyperfine tensor coefficients of Er is similar in both phases, which may be due to the occupation of Er in the same octahedral Ga site. On the other hand, the calculated magnetic exchange interaction between two Er dopants is negative for 𝛼 and positive for 𝛽, indicating an antiferromagnetic (AFM) ground state in the former and a ferromagnetic (FM) ground state in the latter. Large values of Dzyaloshinskii-Moriya interactions (DMIs) are obtained along the 𝑥 direction in the 𝛼 and along the 𝑦 direction in the 𝛽. The large DMI may support exotic magnetic textures, a promising direction for spintronic applications. The analysis of dielectric constants and refractive indices of both pristine and Er-doped phases shows a good agreement with available experimental values. The calculated optical anisotropy is slightly higher in 𝛽 than those in 𝛼, which is due to the involvement of lower symmetry in 𝛽. The crystal field coefficients (CFCs) calculated from DFT are used to analyze 4⁢𝑓 multiplets and 4⁢𝑓 −4⁢𝑓 transitions. Thus calculated lowest energy level of the first excited state to the lowest energy level of the ground state is about 1.53 µ⁢m, which is in a good agreement with available experiments, and it falls within the quantum telecommunication wavelength range.

3-dimensional systems↗

DOE OSTI2021

Computational discovery of stable and metastable ternary oxynitrides

Materials design from first principles enables exploration of uncharted chemical spaces. Extensive computational searches have been performed for mixed-cation ternary compounds, but mixed-anion systems are gaining increased interest as well. Central to computational discovery is the crystal structure prediction, where the trade-off between reliance on prototype structures and size limitations of unconstrained sampling has to be navigated. We approach this challenge by letting two complementary structure sampling approaches compete. We use the kinetically limited minimization approach for high-throughput unconstrained crystal structure prediction in smaller cells up to 21 atoms. On the other hand, ternary—and, more generally, multinary—systems often assume structures formed by atomic ordering on a lattice derived from a binary parent structure. Thus, we additionally sample atomic configurations on prototype lattices with cells up to 56 atoms. Using this approach, we searched 65 different charge-balanced oxide–nitride stoichiometries, including six known systems as the control sample. The convex hull analysis is performed both for the thermodynamic limit and for the case of synthesis with activated nitrogen sources. We identified 34 phases that are either on the convex hull or within a viable energy window for potentially metastable phases. We further performed structure sampling for “missing” binary nitrides whose energies are needed for the convex hull analysis. Among these, we discovered metastable Ce3N4 as a nitride analog of the tetravalent cerium oxide, which becomes stable under slightly activated nitrogen condition µN > +0.07 eV. Given the outsize role of CeO2 in research and application, Ce3N4 is a potentially important discovery.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗